#include <TkAccumulatingSensitiveDetector.h>

Inheritance diagram for TkAccumulatingSensitiveDetector:

Public Member Functions
void	clearHits () override

void	EndOfEvent (G4HCofThisEvent *) override

void	fillHits (edm::PSimHitContainer &, const std::string &) override

bool	ProcessHits (G4Step , G4TouchableHistory ) override

uint32_t	setDetUnitId (const G4Step *) override

	TkAccumulatingSensitiveDetector (const std::string &, const GeometricDet , const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager )

	~TkAccumulatingSensitiveDetector () override

Public Member Functions inherited from SensitiveTkDetector
	SensitiveTkDetector (const std::string &iname, const SensitiveDetectorCatalog &clg)

Public Member Functions inherited from SensitiveDetector
void	EndOfEvent (G4HCofThisEvent *eventHC) override

const std::vector< std::string > &	getNames () const

void	Initialize (G4HCofThisEvent *eventHC) override

bool	isCaloSD () const

	SensitiveDetector (const std::string &iname, const SensitiveDetectorCatalog &, bool calo)

	~SensitiveDetector () override

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 ()

Protected Member Functions
void	update (const BeginOfEvent *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfJob *) override
	This routine will be called when the appropriate signal arrives. More...

void	update (const BeginOfTrack *) override
	This routine will be called when the appropriate signal arrives. More...

Protected Member Functions inherited from SensitiveDetector
TrackInformation *	cmsTrackInformation (const G4Track *aTrack)

Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point) const

Local3DPoint	FinalStepPosition (const G4Step *step, coordinates) const

Local3DPoint	InitialStepPosition (const G4Step *step, coordinates) const

Local3DPoint	LocalPostStepPosition (const G4Step *step) const

Local3DPoint	LocalPreStepPosition (const G4Step *step) const

void	NaNTrap (const G4Step *step) const

void	setNames (const std::vector< std::string > &)

Private Member Functions
bool	closeHit (const G4Step *)

void	createHit (const G4Step *)

bool	newHit (const G4Step *)

void	sendHit ()

void	updateHit (const G4Step *)

Private Attributes
bool	allowZeroEnergyLoss

float	energyCut

float	energyHistoryCut

int	eventno

Local3DPoint	globalEntryPoint

Local3DPoint	globalExitPoint

uint32_t	lastId

int	lastTrack

UpdatablePSimHit *	mySimHit

bool	neverAccumulate

const G4VPhysicalVolume *	oldVolume

const GeometricDet *	pDD_

std::string	pname

bool	printHits

float	px

float	py

float	pz

double	rTracker

double	rTracker2

std::unique_ptr< TrackingSlaveSD >	slaveHighTof

std::unique_ptr< TrackingSlaveSD >	slaveLowTof

std::unique_ptr< const G4ProcessTypeEnumerator >	theG4ProcTypeEnumerator

const SimTrackManager *	theManager

TrackerG4SimHitNumberingScheme *	theNumberingScheme

std::unique_ptr< FrameRotation >	theRotation

float	theTofLimit

double	zTracker

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

Detailed Description

Definition at line 31 of file TkAccumulatingSensitiveDetector.h.

Constructor & Destructor Documentation

◆ TkAccumulatingSensitiveDetector()

TkAccumulatingSensitiveDetector::TkAccumulatingSensitiveDetector	(	const std::string &	name,
		const GeometricDet *	pDD,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

Definition at line 41 of file TkAccumulatingSensitiveDetector.cc.

     : SensitiveTkDetector(name, clg),
       pDD_(pDD),
       theManager(manager),
       rTracker(1200. * CLHEP::mm),
       zTracker(3000. * CLHEP::mm),
       mySimHit(nullptr),
       lastId(0),
       lastTrack(0),
       oldVolume(nullptr),
       px(0.0f),
       py(0.0f),
       pz(0.0f),
       eventno(0),
       pname("") {
   edm::ParameterSet m_TrackerSD = p.getParameter<edm::ParameterSet>("TrackerSD");
   allowZeroEnergyLoss = m_TrackerSD.getParameter<bool>("ZeroEnergyLoss");
   neverAccumulate = m_TrackerSD.getParameter<bool>("NeverAccumulate");
   printHits = m_TrackerSD.getParameter<bool>("PrintHits");
   theTofLimit = m_TrackerSD.getParameter<double>("ElectronicSigmaInNanoSeconds") * 3 * CLHEP::ns;  // 3 sigma
   energyCut =
       m_TrackerSD.getParameter<double>("EnergyThresholdForPersistencyInGeV") * CLHEP::GeV;  //default must be 0.5
   energyHistoryCut =
       m_TrackerSD.getParameter<double>("EnergyThresholdForHistoryInGeV") * CLHEP::GeV;  //default must be 0.05
   rTracker2 = rTracker * rTracker;
  
   // No Rotation given in input, automagically choose one based upon the name
   std::string rotType;
   theRotation = std::make_unique<TrackerFrameRotation>();
   rotType = "TrackerFrameRotation";
  
 #ifdef FAKEFRAMEROTATION
   theRotation.reset(new FakeFrameRotation());
   rotType = "FakeFrameRotation";
 #endif
  
   edm::LogVerbatim("TrackerSim") << " TkAccumulatingSensitiveDetector: "
                                  << " Criteria for Saving Tracker SimTracks: \n"
                                  << " History: " << energyHistoryCut << " MeV; Persistency: " << energyCut
                                  << " MeV;  TofLimit: " << theTofLimit << " ns"
                                  << "\n FrameRotation type " << rotType << " rTracker(cm)= " << rTracker / CLHEP::cm
                                  << " zTracker(cm)= " << zTracker / CLHEP::cm
                                  << " allowZeroEnergyLoss: " << allowZeroEnergyLoss
                                  << " neverAccumulate: " << neverAccumulate << " printHits: " << printHits;
  
   slaveLowTof = std::make_unique<TrackingSlaveSD>(name + "LowTof");
   slaveHighTof = std::make_unique<TrackingSlaveSD>(name + "HighTof");
  
   std::vector<std::string> temp;
   temp.push_back(slaveLowTof.get()->name());
   temp.push_back(slaveHighTof.get()->name());
   setNames(temp);
  
   theG4ProcTypeEnumerator = std::make_unique<G4ProcessTypeEnumerator>();
   theNumberingScheme = nullptr;
 }

References allowZeroEnergyLoss, energyCut, energyHistoryCut, edm::ParameterSet::getParameter(), GeV, Skims_PA_cff::name, neverAccumulate, AlCaHLTBitMon_ParallelJobs::p, printHits, rTracker, rTracker2, SensitiveDetector::setNames(), slaveHighTof, slaveLowTof, AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, theG4ProcTypeEnumerator, theNumberingScheme, theRotation, theTofLimit, and zTracker.

◆ ~TkAccumulatingSensitiveDetector()

TkAccumulatingSensitiveDetector::~TkAccumulatingSensitiveDetector ( )

override

Definition at line 102 of file TkAccumulatingSensitiveDetector.cc.

102 {}

Member Function Documentation

◆ clearHits()

void TkAccumulatingSensitiveDetector::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 354 of file TkAccumulatingSensitiveDetector.cc.

                                                 {
   slaveLowTof.get()->Initialize();
   slaveHighTof.get()->Initialize();
 }

References slaveHighTof, and slaveLowTof.

Referenced by update().

◆ closeHit()

bool TkAccumulatingSensitiveDetector::closeHit ( const G4Step * aStep )

private

Definition at line 331 of file TkAccumulatingSensitiveDetector.cc.

                                                                   {
   const float tolerance2 = 0.0025f;  // (0.5 mm)^2 are allowed between entry and exit
   Local3DPoint theEntryPoint = theRotation.get()->transformPoint(LocalPreStepPosition(aStep));
   LogDebug("TrackerSimDebug") << " closeHit: distance = " << (mySimHit->exitPoint() - theEntryPoint).mag();
  
   return ((mySimHit->exitPoint() - theEntryPoint).mag2() < tolerance2) ? true : false;
 }

References PSimHit::exitPoint(), SensitiveDetector::LocalPreStepPosition(), LogDebug, mag(), mag2(), mySimHit, and theRotation.

Referenced by newHit().

◆ createHit()

void TkAccumulatingSensitiveDetector::createHit ( const G4Step * aStep )

private

Definition at line 204 of file TkAccumulatingSensitiveDetector.cc.

                                                                    {
   // VI: previous hit should be already deleted
   //     in past here was a check if a hit is inside a sensitive detector,
   //     this is not needed, because call to senstive detector happens
   //     only inside the volume
   const G4Track* theTrack = aStep->GetTrack();
   Local3DPoint theExitPoint = theRotation.get()->transformPoint(LocalPostStepPosition(aStep));
   Local3DPoint theEntryPoint;
   //
   //  Check particle type - for gamma and neutral hadrons energy deposition
   //  should be local (VI)
   //
   if (0.0 == theTrack->GetDefinition()->GetPDGCharge()) {
     theEntryPoint = theExitPoint;
   } else {
     theEntryPoint = theRotation.get()->transformPoint(LocalPreStepPosition(aStep));
   }
  
   //
   //    This allows to send he skipEvent if it is outside!
   //
   const G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
   float thePabs = preStepPoint->GetMomentum().mag() / GeV;
   float theTof = preStepPoint->GetGlobalTime() / nanosecond;
   float theEnergyLoss = aStep->GetTotalEnergyDeposit() / GeV;
   int theParticleType = G4TrackToParticleID::particleID(theTrack);
   uint32_t theDetUnitId = setDetUnitId(aStep);
   int theTrackID = theTrack->GetTrackID();
   if (theDetUnitId == 0) {
     edm::LogWarning("TkAccumulatingSensitiveDetector::createHit") << " theDetUnitId is not valid for " << GetName();
     throw cms::Exception("TkAccumulatingSensitiveDetector::createHit")
         << "cannot get theDetUnitId for G4Track " << theTrackID;
   }
  
   // To whom assign the Hit?
   // First iteration: if the track is to be stored, use the current number;
   // otherwise, get to the mother
   unsigned int theTrackIDInsideTheSimHit = theTrackID;
  
   const TrackInformation* temp = cmsTrackInformation(theTrack);
   if (!temp->storeTrack()) {
     // Go to the mother!
     theTrackIDInsideTheSimHit = theTrack->GetParentID();
     LogDebug("TrackerSimDebug") << " TkAccumulatingSensitiveDetector::createHit(): setting the TrackID from "
                                 << theTrackIDInsideTheSimHit << " to the mother one " << theTrackIDInsideTheSimHit
                                 << " " << theEnergyLoss;
   } else {
     LogDebug("TrackerSimDebug") << " TkAccumulatingSensitiveDetector:createHit(): leaving the current TrackID "
                                 << theTrackIDInsideTheSimHit;
   }
  
   const G4ThreeVector& gmd = preStepPoint->GetMomentumDirection();
   // convert it to local frame
   G4ThreeVector lmd =
       ((G4TouchableHistory*)(preStepPoint->GetTouchable()))->GetHistory()->GetTopTransform().TransformAxis(gmd);
   Local3DPoint lnmd = theRotation.get()->transformPoint(ConvertToLocal3DPoint(lmd));
   float theThetaAtEntry = lnmd.theta();
   float thePhiAtEntry = lnmd.phi();
  
   mySimHit = new UpdatablePSimHit(theEntryPoint,
                                   theExitPoint,
                                   thePabs,
                                   theTof,
                                   theEnergyLoss,
                                   theParticleType,
                                   theDetUnitId,
                                   theTrackIDInsideTheSimHit,
                                   theThetaAtEntry,
                                   thePhiAtEntry,
                                   theG4ProcTypeEnumerator.get()->processId(theTrack->GetCreatorProcess()));
   lastId = theDetUnitId;
   lastTrack = theTrackID;
  
   // only for debugging
   if (printHits) {
     // point on Geant4 unit (mm)
     globalEntryPoint = ConvertToLocal3DPoint(preStepPoint->GetPosition());
     globalExitPoint = ConvertToLocal3DPoint(aStep->GetPostStepPoint()->GetPosition());
     // in CMS unit (GeV)
     px = preStepPoint->GetMomentum().x() / CLHEP::GeV;
     py = preStepPoint->GetMomentum().y() / CLHEP::GeV;
     pz = preStepPoint->GetMomentum().z() / CLHEP::GeV;
     oldVolume = preStepPoint->GetPhysicalVolume();
     pname = theTrack->GetDefinition()->GetParticleName();
     LogDebug("TrackerSimDebug") << " Created PSimHit: " << pname << " " << mySimHit->detUnitId() << " "
                                 << mySimHit->trackId() << " " << theTrackID
                                 << " p= " << aStep->GetPreStepPoint()->GetMomentum().mag() << " "
                                 << mySimHit->energyLoss() << " " << mySimHit->entryPoint() << " "
                                 << mySimHit->exitPoint();
   }
 }

References SensitiveDetector::cmsTrackInformation(), SensitiveDetector::ConvertToLocal3DPoint(), PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), Exception, PSimHit::exitPoint(), GeV, globalEntryPoint, globalExitPoint, lastId, lastTrack, SensitiveDetector::LocalPostStepPosition(), SensitiveDetector::LocalPreStepPosition(), LogDebug, mySimHit, oldVolume, G4TrackToParticleID::particleID(), PV3DBase< T, PVType, FrameType >::phi(), pname, printHits, px, py, pz, setDetUnitId(), groupFilesInBlocks::temp, theG4ProcTypeEnumerator, theRotation, PV3DBase< T, PVType, FrameType >::theta(), and PSimHit::trackId().

Referenced by ProcessHits().

◆ EndOfEvent()

void TkAccumulatingSensitiveDetector::EndOfEvent ( G4HCofThisEvent * )

override

Definition at line 339 of file TkAccumulatingSensitiveDetector.cc.

                                                                  {
   LogDebug("TrackerSimDebug") << " Saving the last hit in a ROU " << GetName();
   if (mySimHit != nullptr)
     sendHit();
 }

References LogDebug, mySimHit, and sendHit().

◆ fillHits()

void TkAccumulatingSensitiveDetector::fillHits	(	edm::PSimHitContainer &	cc,
		const std::string &	hname
	)

overridevirtual

Implements SensitiveTkDetector.

Definition at line 359 of file TkAccumulatingSensitiveDetector.cc.

                                                                                               {
   if (slaveLowTof.get()->name() == hname) {
     cc = slaveLowTof.get()->hits();
   } else if (slaveHighTof.get()->name() == hname) {
     cc = slaveHighTof.get()->hits();
   }
 }

References slaveHighTof, and slaveLowTof.

◆ newHit()

bool TkAccumulatingSensitiveDetector::newHit ( const G4Step * aStep )

private

Definition at line 314 of file TkAccumulatingSensitiveDetector.cc.

                                                                 {
   const G4Track* theTrack = aStep->GetTrack();
  
   // for neutral particles do not merge hits (V.I.)
   if (0.0 == theTrack->GetDefinition()->GetPDGCharge())
     return true;
  
   uint32_t theDetUnitId = setDetUnitId(aStep);
   int theTrackID = theTrack->GetTrackID();
  
   LogDebug("TrackerSimDebug") << "newHit: OLD(detID,trID) = (" << lastId << "," << lastTrack << "), NEW = ("
                               << theDetUnitId << "," << theTrackID << ") Step length(mm)= " << aStep->GetStepLength()
                               << " Edep= " << aStep->GetTotalEnergyDeposit()
                               << " p= " << aStep->GetPreStepPoint()->GetMomentum().mag();
   return ((theTrackID == lastTrack) && (lastId == theDetUnitId) && closeHit(aStep)) ? false : true;
 }

References closeHit(), funct::false, lastId, lastTrack, LogDebug, and setDetUnitId().

Referenced by ProcessHits().

◆ ProcessHits()

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

overridevirtual

Implements SensitiveDetector.

Definition at line 104 of file TkAccumulatingSensitiveDetector.cc.

                                                                                     {
   LogDebug("TrackerSimDebug") << " Entering a new Step " << aStep->GetTotalEnergyDeposit() << " "
                               << aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName();
  
   if (aStep->GetTotalEnergyDeposit() > 0. || allowZeroEnergyLoss) {
     if (!mySimHit) {
       createHit(aStep);
     } else if (neverAccumulate || newHit(aStep)) {
       sendHit();
       createHit(aStep);
     } else {
       updateHit(aStep);
     }
     return true;
   }
   return false;
 }

References allowZeroEnergyLoss, createHit(), LogDebug, mySimHit, neverAccumulate, newHit(), sendHit(), and updateHit().

Referenced by LaserAlignmentSimulation::update().

◆ sendHit()

void TkAccumulatingSensitiveDetector::sendHit ( )

private

Definition at line 171 of file TkAccumulatingSensitiveDetector.cc.

                                               {
   if (mySimHit == nullptr)
     return;
   if (printHits) {
     TkSimHitPrinter thePrinter("TkHitPositionOSCAR.dat");
     thePrinter.startNewSimHit(GetName(),
                               oldVolume->GetLogicalVolume()->GetName(),
                               mySimHit->detUnitId(),
                               mySimHit->trackId(),
                               lastTrack,
                               eventno);
     thePrinter.printLocal(mySimHit->entryPoint(), mySimHit->exitPoint());
     thePrinter.printGlobal(globalEntryPoint, globalExitPoint);
     thePrinter.printHitData(pname, mySimHit->pabs(), mySimHit->energyLoss(), mySimHit->timeOfFlight());
     thePrinter.printGlobalMomentum(px, py, pz);
     LogDebug("TrackerSimDebug") << " Storing PSimHit: " << mySimHit->detUnitId() << " " << mySimHit->trackId() << " "
                                 << mySimHit->energyLoss() << " " << mySimHit->entryPoint() << " "
                                 << mySimHit->exitPoint();
   }
  
   if (mySimHit->timeOfFlight() < theTofLimit) {
     slaveLowTof.get()->processHits(*mySimHit);  // implicit conversion (slicing) to PSimHit!!!
   } else {
     slaveHighTof.get()->processHits(*mySimHit);  // implicit conversion (slicing) to PSimHit!!!
   }
   //
   // clean up
   delete mySimHit;
   mySimHit = nullptr;
   lastTrack = 0;
   lastId = 0;
 }

References PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), eventno, PSimHit::exitPoint(), globalEntryPoint, globalExitPoint, lastId, lastTrack, LogDebug, mySimHit, oldVolume, PSimHit::pabs(), pname, TkSimHitPrinter::printGlobal(), TkSimHitPrinter::printGlobalMomentum(), TkSimHitPrinter::printHitData(), printHits, TkSimHitPrinter::printLocal(), px, py, pz, slaveHighTof, slaveLowTof, TkSimHitPrinter::startNewSimHit(), theTofLimit, PSimHit::timeOfFlight(), and PSimHit::trackId().

Referenced by EndOfEvent(), and ProcessHits().

◆ setDetUnitId()

uint32_t TkAccumulatingSensitiveDetector::setDetUnitId ( const G4Step * step )

overridevirtual

Implements SensitiveDetector.

Definition at line 122 of file TkAccumulatingSensitiveDetector.cc.

                                                                          {
   return theNumberingScheme->g4ToNumberingScheme(step->GetPreStepPoint()->GetTouchable());
 }

References TrackerG4SimHitNumberingScheme::g4ToNumberingScheme(), and theNumberingScheme.

Referenced by createHit(), and newHit().

◆ update() [1/3]

void TkAccumulatingSensitiveDetector::update ( const BeginOfEvent * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 345 of file TkAccumulatingSensitiveDetector.cc.

                                                                   {
   clearHits();
   eventno = (*i)()->GetEventID();
   delete mySimHit;
   mySimHit = nullptr;
 }

References clearHits(), eventno, and mySimHit.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

◆ update() [2/3]

void TkAccumulatingSensitiveDetector::update ( const BeginOfJob * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 352 of file TkAccumulatingSensitiveDetector.cc.

352 { theNumberingScheme = &(numberingScheme(*pDD_)); }

References numberingScheme(), pDD_, and theNumberingScheme.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

◆ update() [3/3]

void TkAccumulatingSensitiveDetector::update ( const BeginOfTrack * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 126 of file TkAccumulatingSensitiveDetector.cc.

                                                                     {
   const G4Track* gTrack = (*bot)();
  
 #ifdef DUMPPROCESSES
   if (gTrack->GetCreatorProcess()) {
     edm::LogVerbatim("TrackerSim") << " -> PROCESS CREATOR : " << gTrack->GetCreatorProcess()->GetProcessName();
   } else {
     edm::LogVerbatim("TrackerSim") << " -> No Creator process";
   }
 #endif
  
   //
   //Position
   //
   const G4ThreeVector& pos = gTrack->GetPosition();
   LogDebug("TrackerSimDebug") << " update(..) of " << gTrack->GetDefinition()->GetParticleName()
                               << " trackID= " << gTrack->GetTrackID() << " E(MeV)= " << gTrack->GetKineticEnergy()
                               << " Ecut= " << energyCut << " R(mm)= " << pos.perp() << " Z(mm)= " << pos.z();
  
   //
   // Check if in Tracker Volume
   //
   if (pos.x() * pos.x() + pos.y() * pos.y() < rTracker2 && std::abs(pos.z()) < zTracker) {
     //
     // inside the Tracker
     //
     TrackInformation* info = nullptr;
     if (gTrack->GetKineticEnergy() > energyCut) {
       info = cmsTrackInformation(gTrack);
       info->storeTrack(true);
     }
     //
     // Save History?
     //
     if (gTrack->GetKineticEnergy() > energyHistoryCut) {
       if (!info) {
         info = cmsTrackInformation(gTrack);
       }
       info->putInHistory();
       LogDebug("TrackerSimDebug") << " Track inside the tracker selected for HISTORY"
                                   << " Track ID= " << gTrack->GetTrackID();
     }
   }
 }

References funct::abs(), SensitiveDetector::cmsTrackInformation(), energyCut, energyHistoryCut, info(), LogDebug, rTracker2, and zTracker.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

◆ updateHit()

void TkAccumulatingSensitiveDetector::updateHit ( const G4Step * aStep )

private

Definition at line 296 of file TkAccumulatingSensitiveDetector.cc.

                                                                    {
   // VI: in past here was a check if a hit is inside a sensitive detector,
   //     this is not needed, because call to senstive detector happens
   //     only inside the volume
   Local3DPoint theExitPoint = theRotation.get()->transformPoint(LocalPostStepPosition(aStep));
   float theEnergyLoss = aStep->GetTotalEnergyDeposit() / GeV;
   mySimHit->setExitPoint(theExitPoint);
   mySimHit->addEnergyLoss(theEnergyLoss);
   if (printHits) {
     globalExitPoint = ConvertToLocal3DPoint(aStep->GetPostStepPoint()->GetPosition());
     LogDebug("TrackerSimDebug") << " updateHit: for " << aStep->GetTrack()->GetDefinition()->GetParticleName()
                                 << " trackID= " << aStep->GetTrack()->GetTrackID() << " deltaEloss= " << theEnergyLoss
                                 << "\n Updated PSimHit: " << mySimHit->detUnitId() << " " << mySimHit->trackId() << " "
                                 << mySimHit->energyLoss() << " " << mySimHit->entryPoint() << " "
                                 << mySimHit->exitPoint();
   }
 }