Bcm1fSD Class Reference

#include <Bcm1fSD.h>

Inheritance diagram for Bcm1fSD:

Public Member Functions
	Bcm1fSD (std::string, const DDCompactView &, SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)
virtual void	EndOfEvent (G4HCofThisEvent *)
void	fillHits (edm::PSimHitContainer &, std::string use)
virtual bool	ProcessHits (G4Step *, G4TouchableHistory *)
virtual uint32_t	setDetUnitId (G4Step *)
virtual	~Bcm1fSD ()
Public Member Functions inherited from SensitiveTkDetector
	SensitiveTkDetector (std::string &iname, const DDCompactView &cpv, SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)
Public Member Functions inherited from SensitiveDetector
virtual void	AssignSD (std::string &vname)
Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point)
Local3DPoint	FinalStepPosition (G4Step *s, coordinates)
virtual std::vector< std::string >	getNames ()
virtual void	Initialize (G4HCofThisEvent *eventHC)
Local3DPoint	InitialStepPosition (G4Step *s, coordinates)
std::string	nameOfSD ()
void	NaNTrap (G4Step *step)
void	Register ()
	SensitiveDetector (std::string &iname, const DDCompactView &cpv, SensitiveDetectorCatalog &, edm::ParameterSet const &p)
virtual	~SensitiveDetector ()
Public Member Functions inherited from Observer< const BeginOfEvent * >
	Observer ()
void	slotForUpdate (const BeginOfEvent *iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfTrack * >
	Observer ()
void	slotForUpdate (const BeginOfTrack *iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfJob * >
	Observer ()
void	slotForUpdate (const BeginOfJob *iT)
virtual	~Observer ()

Private Member Functions
virtual void	clearHits ()
virtual bool	closeHit (G4Step *)
virtual void	createHit (G4Step *)
TrackInformation *	getOrCreateTrackInformation (const G4Track *)
virtual bool	newHit (G4Step *)
virtual void	sendHit ()
void	update (const BeginOfEvent *)
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfTrack *)
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfJob *)
	This routine will be called when the appropriate signal arrives. More...
virtual void	updateHit (G4Step *)

Private Attributes
float	energyCut
float	energyHistoryCut
int	eventno
Local3DPoint	globalEntryPoint
Local3DPoint	globalExitPoint
uint32_t	lastId
unsigned int	lastTrack
G4TrackToParticleID *	myG4TrackToParticleID
std::string	myName
UpdatablePSimHit *	mySimHit
G4VPhysicalVolume *	oldVolume
std::string	pname
TrackingSlaveSD *	slave
G4ProcessTypeEnumerator *	theG4ProcessTypeEnumerator

Additional Inherited Members
Public Types inherited from SensitiveDetector
enum	coordinates { WorldCoordinates, LocalCoordinates }

Detailed Description

Definition at line 30 of file Bcm1fSD.h.

Constructor & Destructor Documentation

Bcm1fSD::Bcm1fSD	(	std::string	name,
		const DDCompactView &	cpv,
		SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

Definition at line 33 of file Bcm1fSD.cc.

References SensitiveDetector::AssignSD(), energyCut, energyHistoryCut, edm::ParameterSet::getParameter(), SensitiveDetectorCatalog::logicalNames(), myG4TrackToParticleID, SensitiveDetector::Register(), slave, and theG4ProcessTypeEnumerator.

                                                       : 
  SensitiveTkDetector(name, cpv, clg, p), myName(name), mySimHit(0),
  oldVolume(0), lastId(0), lastTrack(0), eventno(0) {
  
  edm::ParameterSet m_TrackerSD = p.getParameter<edm::ParameterSet>("Bcm1fSD");
  energyCut           = m_TrackerSD.getParameter<double>("EnergyThresholdForPersistencyInGeV")*GeV; //default must be 0.5 (?)
  energyHistoryCut    = m_TrackerSD.getParameter<double>("EnergyThresholdForHistoryInGeV")*GeV;//default must be 0.05 (?)

  edm::LogInfo("Bcm1fSD") <<"Criteria for Saving Tracker SimTracks: \n "
                       <<" History: "<<energyHistoryCut<< " MeV ; Persistency: "<< energyCut<<" MeV\n"
                       <<" Constructing a Bcm1fSD with ";

  slave  = new TrackingSlaveSD(name);
  
  // Now attach the right detectors (LogicalVolumes) to me
  std::vector<std::string>  lvNames = clg.logicalNames(name);
  this->Register();
  for (std::vector<std::string>::iterator it = lvNames.begin(); it != lvNames.end(); it++)
  {
     edm::LogInfo("Bcm1fSD")<< name << " attaching LV " << *it;
     this->AssignSD(*it);
  }

  theG4ProcessTypeEnumerator = new G4ProcessTypeEnumerator;
  myG4TrackToParticleID = new G4TrackToParticleID;
}

Bcm1fSD::~Bcm1fSD ( )

virtual

Definition at line 64 of file Bcm1fSD.cc.

References myG4TrackToParticleID, slave, and theG4ProcessTypeEnumerator.

                  { 
  delete slave;
  delete theG4ProcessTypeEnumerator;
  delete myG4TrackToParticleID;
}

Member Function Documentation

void Bcm1fSD::clearHits ( )

privatevirtual

Implements SensitiveDetector.

Definition at line 284 of file Bcm1fSD.cc.

References TrackingSlaveSD::Initialize(), and slave.

Referenced by update().

                        {
    slave->Initialize();
}

bool Bcm1fSD::closeHit ( G4Step *  aStep )

privatevirtual

Definition at line 211 of file Bcm1fSD.cc.

References PSimHit::exitPoint(), SensitiveDetector::InitialStepPosition(), SensitiveDetector::LocalCoordinates, LogDebug, mag(), and mySimHit.

Referenced by newHit().

                                     {

  if (mySimHit == 0) return false; 
  const float tolerance = 0.05 * mm; // 50 micron are allowed between the exit 
  // point of the current hit and the entry point of the new hit
  Local3DPoint theEntryPoint = SensitiveDetector::InitialStepPosition(aStep,LocalCoordinates);  
  LogDebug("Bcm1fSD")<< " closeHit: distance = " << (mySimHit->exitPoint()-theEntryPoint).mag();

  if ((mySimHit->exitPoint()-theEntryPoint).mag()<tolerance) return true;
  return false;
}

void Bcm1fSD::createHit ( G4Step *  aStep )

privatevirtual

Definition at line 223 of file Bcm1fSD.cc.

References SensitiveDetector::ConvertToLocal3DPoint(), PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), SensitiveDetector::FinalStepPosition(), globalEntryPoint, globalExitPoint, SensitiveDetector::InitialStepPosition(), lastId, lastTrack, SensitiveDetector::LocalCoordinates, LogDebug, myG4TrackToParticleID, mySimHit, oldVolume, G4TrackToParticleID::particleID(), PV3DBase< T, PVType, FrameType >::phi(), pname, G4ProcessTypeEnumerator::processId(), setDetUnitId(), theG4ProcessTypeEnumerator, PV3DBase< T, PVType, FrameType >::theta(), PSimHit::trackId(), findQualityFiles::v, and SensitiveDetector::WorldCoordinates.

Referenced by ProcessHits().

                                      {

  if (mySimHit != 0) {
    delete mySimHit;
    mySimHit=0;
  }
    
  G4Track * theTrack  = aStep->GetTrack(); 
  G4VPhysicalVolume * v = aStep->GetPreStepPoint()->GetPhysicalVolume();

  Local3DPoint theEntryPoint = SensitiveDetector::InitialStepPosition(aStep,LocalCoordinates);  
  Local3DPoint theExitPoint  = SensitiveDetector::FinalStepPosition(aStep,LocalCoordinates); 
  
  float thePabs             = aStep->GetPreStepPoint()->GetMomentum().mag()/GeV;
  float theTof              = aStep->GetPreStepPoint()->GetGlobalTime()/nanosecond;
  float theEnergyLoss       = aStep->GetTotalEnergyDeposit()/GeV;
  int theParticleType       = myG4TrackToParticleID->particleID(theTrack);
  uint32_t theDetUnitId     = setDetUnitId(aStep);
  
  globalEntryPoint = SensitiveDetector::InitialStepPosition(aStep,WorldCoordinates);
  globalExitPoint = SensitiveDetector::FinalStepPosition(aStep,WorldCoordinates);
  pname = theTrack->GetDynamicParticle()->GetDefinition()->GetParticleName();
  
  unsigned int theTrackID = theTrack->GetTrackID();
  
  G4ThreeVector gmd  = aStep->GetPreStepPoint()->GetMomentumDirection();
  // convert it to local frame
  G4ThreeVector lmd = ((G4TouchableHistory *)(aStep->GetPreStepPoint()->GetTouchable()))->GetHistory()->GetTopTransform().TransformAxis(gmd);
  Local3DPoint lnmd = ConvertToLocal3DPoint(lmd);
  float theThetaAtEntry = lnmd.theta();
  float thePhiAtEntry = lnmd.phi();
    
  mySimHit = new UpdatablePSimHit(theEntryPoint,theExitPoint,thePabs,theTof,
                  theEnergyLoss,theParticleType,theDetUnitId,
                  theTrackID,theThetaAtEntry,thePhiAtEntry,
                  theG4ProcessTypeEnumerator->processId(theTrack->GetCreatorProcess()));  

  LogDebug("Bcm1fSD") << " Created PSimHit: " << pname << " " 
                   << mySimHit->detUnitId() << " " << mySimHit->trackId()
                   << " " << mySimHit->energyLoss() << " " 
                   << mySimHit->entryPoint() << " " << mySimHit->exitPoint();
  lastId = theDetUnitId;
  lastTrack = theTrackID;
  oldVolume = v;
}

void Bcm1fSD::EndOfEvent ( G4HCofThisEvent *  )

virtual

Reimplemented from SensitiveDetector.

Definition at line 153 of file Bcm1fSD.cc.

References LogDebug, myName, mySimHit, and sendHit().

                                          {
  
  LogDebug("Bcm1fSD")<< " Saving the last hit in a ROU " << myName;

  if (mySimHit == 0) return;
  sendHit();
}

void Bcm1fSD::fillHits ( edm::PSimHitContainer &  c, std::string  use  )

virtual

Implements SensitiveTkDetector.

Definition at line 161 of file Bcm1fSD.cc.

References TrackingSlaveSD::hits(), n, TrackingSlaveSD::name(), and slave.

                                                         {
  if (slave->name() == n)  c=slave->hits();
}

TrackInformation * Bcm1fSD::getOrCreateTrackInformation ( const G4Track *  gTrack )

private

Definition at line 288 of file Bcm1fSD.cc.

References info(), and groupFilesInBlocks::temp.

Referenced by ProcessHits().

                                                                             {
  G4VUserTrackInformation* temp = gTrack->GetUserInformation();
  if (temp == 0){
    edm::LogError("Bcm1fSD") <<" ERROR: no G4VUserTrackInformation available";
    abort();
  }else{
    TrackInformation* info = dynamic_cast<TrackInformation*>(temp);
    if (info == 0){
      edm::LogError("Bcm1fSD") <<" ERROR: TkSimTrackSelection: the UserInformation does not appear to be a TrackInformation";
      abort();
    }
    return info;
  }
}

bool Bcm1fSD::newHit ( G4Step *  aStep )

privatevirtual

Definition at line 197 of file Bcm1fSD.cc.

References closeHit(), lastId, lastTrack, LogDebug, mySimHit, and setDetUnitId().

Referenced by ProcessHits().

                                   {

  G4Track * theTrack = aStep->GetTrack(); 
  uint32_t theDetUnitId = setDetUnitId(aStep);
  unsigned int theTrackID = theTrack->GetTrackID();

  LogDebug("Bcm1fSD") << " OLD (d,t) = (" << lastId << "," << lastTrack 
                   << "), new = (" << theDetUnitId << "," << theTrackID << ") return "
                   << ((theTrackID == lastTrack) && (lastId == theDetUnitId));
  if ((mySimHit != 0) && (theTrackID == lastTrack) && (lastId == theDetUnitId) && closeHit(aStep))
    return false;
  return true;
}

bool Bcm1fSD::ProcessHits ( G4Step *  aStep, G4TouchableHistory *    )

virtual

Implements SensitiveDetector.

Definition at line 70 of file Bcm1fSD.cc.

References createHit(), energyCut, energyHistoryCut, getOrCreateTrackInformation(), info(), lastTrack, LogDebug, newHit(), TrackInformation::putInHistory(), sendHit(), TrackInformation::storeTrack(), and updateHit().

                                                               {

  LogDebug("Bcm1fSD") << " Entering a new Step " 
            << aStep->GetTotalEnergyDeposit() << " " 
            << aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName();

  G4Track * gTrack  = aStep->GetTrack(); 
  if ((unsigned int)(gTrack->GetTrackID()) != lastTrack) {
    
    if (gTrack->GetKineticEnergy() > energyCut){
      TrackInformation* info = getOrCreateTrackInformation(gTrack);
      info->storeTrack(true);
    }
    if (gTrack->GetKineticEnergy() > energyHistoryCut){
      TrackInformation* info = getOrCreateTrackInformation(gTrack);
      info->putInHistory();
    }
  }

  if (aStep->GetTotalEnergyDeposit()>0.) {
    if (newHit(aStep) == true) {
      sendHit();
      createHit(aStep);
    } else {
      updateHit(aStep);
    }
    return true;
  }
  return false;
}

void Bcm1fSD::sendHit ( )

privatevirtual

Definition at line 165 of file Bcm1fSD.cc.

References PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), lastId, lastTrack, LogDebug, mySimHit, pname, TrackingSlaveSD::processHits(), slave, and PSimHit::trackId().

Referenced by EndOfEvent(), and ProcessHits().

                      {  
  if (mySimHit == 0) return;
  LogDebug("Bcm1fSD") << " Storing PSimHit: " << pname << " " << mySimHit->detUnitId() 
                   << " " << mySimHit->trackId() << " " << mySimHit->energyLoss() 
                   << " " << mySimHit->entryPoint() << " " << mySimHit->exitPoint();
    
  slave->processHits(*mySimHit); 

  // clean up
  delete mySimHit;
  mySimHit = 0;
  lastTrack = 0;
  lastId = 0;
}

uint32_t Bcm1fSD::setDetUnitId ( G4Step *  aStep )

virtual

Implements SensitiveDetector.

Definition at line 101 of file Bcm1fSD.cc.

References gather_cfg::cout, testEve_cfg::level, and LogDebug.

Referenced by createHit(), and newHit().

                                              {
 
  unsigned int detId = 0;

  LogDebug("Bcm1fSD")<< " DetID = "<<detId; 
  
  //Find number of levels
  const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
  int level = 0;
  if (touch) level = ((touch->GetHistoryDepth())+1);
  
  //Get name and copy numbers
  if ( level > 1 ) {
     
     G4String sensorName   = touch->GetVolume(0)->GetName();
     G4String diamondName  = touch->GetVolume(1)->GetName();
     G4String detectorName = touch->GetVolume(2)->GetName();
     G4String volumeName   = touch->GetVolume(3)->GetName();
     
     if ( sensorName != "BCM1FSensor" )
        std::cout << " Bcm1fSD::setDetUnitId -w- Sensor name not BCM1FSensor " << std::endl;
     if ( detectorName != "BCM1F" )
        std::cout << " Bcm1fSD::setDetUnitId -w- Detector name not BCM1F " << std::endl;
     
     int sensorNo    = touch->GetReplicaNumber(0);
     int diamondNo   = touch->GetReplicaNumber(1);
//     int detectorNo  = touch->GetReplicaNumber(2);
     int volumeNo    = touch->GetReplicaNumber(3);
     
     // Detector ID definition
     // detId = XYYZ
     // X  = volume,  1: +Z, 2: -Z
     // YY = diamond, 01-12, 12: phi = 90 deg, numbering clockwise when looking from the IP
     // Z  = sensor,  1 or 2, clockwise when looking from the IP
     
//     detId = 10*volumeNo + diamondNo;
     detId = 1000*volumeNo + 10*diamondNo + sensorNo;
     
//     for ( int ii = 0; ii < level; ii++ ) {
//       int i      = level - ii - 1;
//       G4String name   = touch->GetVolume(i)->GetName();
//       int copyno = touch->GetReplicaNumber(i);
//       const G4ThreeVector trans = touch->GetVolume(i)->GetTranslation();
//       std::cout << " detId = " << detId << " name = " << name << " " << copyno
//             << ",     translation x, y, z = " << trans.x() << ", " <<trans.y() << ", " <<trans.z()
//             << ", eta = " << trans.eta() << std::endl;
//     }
  }

  return detId;
}

void Bcm1fSD::update ( const BeginOfEvent *  )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 271 of file Bcm1fSD.cc.

References clearHits(), eventno, and mySimHit.

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                           {

  clearHits();
  eventno = (*i)()->GetEventID();
  mySimHit = 0;
}

void Bcm1fSD::update ( const BeginOfTrack *  )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 278 of file Bcm1fSD.cc.

References pname.

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                            {

  const G4Track* gTrack = (*bot)();
  pname = gTrack->GetDynamicParticle()->GetDefinition()->GetParticleName();
}

void Bcm1fSD::update ( const BeginOfJob *  )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 269 of file Bcm1fSD.cc.

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{ }

void Bcm1fSD::updateHit ( G4Step *  aStep )

privatevirtual

Definition at line 180 of file Bcm1fSD.cc.

References UpdatablePSimHit::addEnergyLoss(), PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), SensitiveDetector::FinalStepPosition(), globalExitPoint, SensitiveDetector::LocalCoordinates, LogDebug, mySimHit, pname, UpdatablePSimHit::setExitPoint(), PSimHit::trackId(), and SensitiveDetector::WorldCoordinates.

Referenced by ProcessHits().

                                      {

  Local3DPoint theExitPoint = SensitiveDetector::FinalStepPosition(aStep,LocalCoordinates); 
  float theEnergyLoss = aStep->GetTotalEnergyDeposit()/GeV;
  mySimHit->setExitPoint(theExitPoint);
  LogDebug("Bcm1fSD")<< " Before : " << mySimHit->energyLoss();
  mySimHit->addEnergyLoss(theEnergyLoss);
  globalExitPoint = SensitiveDetector::FinalStepPosition(aStep,WorldCoordinates);

  LogDebug("Bcm1fSD") << " Updating: new exitpoint " << pname << " " 
                   << theExitPoint << " new energy loss " << theEnergyLoss 
                   << "\n Updated PSimHit: " << mySimHit->detUnitId() 
                   << " " << mySimHit->trackId()
                   << " " << mySimHit->energyLoss() << " " 
                   << mySimHit->entryPoint() << " " << mySimHit->exitPoint();
}

Member Data Documentation

float Bcm1fSD::energyCut

private

Definition at line 68 of file Bcm1fSD.h.

Referenced by Bcm1fSD(), and ProcessHits().

float Bcm1fSD::energyHistoryCut

private

Definition at line 69 of file Bcm1fSD.h.

Referenced by Bcm1fSD(), and ProcessHits().

int Bcm1fSD::eventno

private

Definition at line 76 of file Bcm1fSD.h.

Referenced by update().

Local3DPoint Bcm1fSD::globalEntryPoint

private

Definition at line 71 of file Bcm1fSD.h.

Referenced by createHit().

Local3DPoint Bcm1fSD::globalExitPoint

private

Definition at line 72 of file Bcm1fSD.h.

Referenced by createHit(), and updateHit().

uint32_t Bcm1fSD::lastId

private

Definition at line 74 of file Bcm1fSD.h.

Referenced by createHit(), newHit(), and sendHit().

unsigned int Bcm1fSD::lastTrack

private

Definition at line 75 of file Bcm1fSD.h.

Referenced by createHit(), newHit(), ProcessHits(), and sendHit().

G4TrackToParticleID* Bcm1fSD::myG4TrackToParticleID

private

Definition at line 65 of file Bcm1fSD.h.

Referenced by Bcm1fSD(), createHit(), and ~Bcm1fSD().

std::string Bcm1fSD::myName

private

Definition at line 66 of file Bcm1fSD.h.

Referenced by EndOfEvent().

UpdatablePSimHit* Bcm1fSD::mySimHit

private

Definition at line 67 of file Bcm1fSD.h.

Referenced by closeHit(), createHit(), EndOfEvent(), newHit(), sendHit(), update(), and updateHit().

G4VPhysicalVolume* Bcm1fSD::oldVolume

private

Definition at line 73 of file Bcm1fSD.h.

Referenced by createHit().

std::string Bcm1fSD::pname

private

Definition at line 77 of file Bcm1fSD.h.

Referenced by createHit(), sendHit(), update(), and updateHit().

TrackingSlaveSD* Bcm1fSD::slave

private

Definition at line 63 of file Bcm1fSD.h.

Referenced by Bcm1fSD(), clearHits(), fillHits(), sendHit(), and ~Bcm1fSD().

G4ProcessTypeEnumerator* Bcm1fSD::theG4ProcessTypeEnumerator

private

Definition at line 64 of file Bcm1fSD.h.

Referenced by Bcm1fSD(), createHit(), and ~Bcm1fSD().