#include <HFWedgeSD.h>

Inheritance diagram for HFWedgeSD:

Public Member Functions
void	clear () override

void	clearHits () override

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

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

	HFWedgeSD (const std::string &, const SensitiveDetectorCatalog &clg, const SimTrackManager *)

void	Initialize (G4HCofThisEvent *HCE) override

void	PrintAll () override

bool	ProcessHits (G4Step step, G4TouchableHistory tHistory) override

uint32_t	setDetUnitId (const G4Step *) override

	~HFWedgeSD () override

Public Member Functions inherited from SensitiveCaloDetector
virtual void	reset ()

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

Protected Member Functions
HFShowerG4Hit *	createNewHit ()

G4bool	hitExists ()

void	updateHit (HFShowerG4Hit *)

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 Attributes
HFShowerG4Hit *	currentHit

int	currentID

double	edep

G4ThreeVector	globalPos

int	hcID

std::map< int, HFShowerG4Hit * >	hitMap

G4ThreeVector	localPos

G4ThreeVector	momDir

int	previousID

HFShowerG4HitsCollection *	theHC

double	time

int	trackID

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

Detailed Description

Definition at line 20 of file HFWedgeSD.h.

Constructor & Destructor Documentation

HFWedgeSD::HFWedgeSD	(	const std::string &	iname,
		const SensitiveDetectorCatalog &	clg,
		const SimTrackManager *	manager
	)

explicit

Definition at line 21 of file HFWedgeSD.cc.

     : SensitiveCaloDetector(iname, clg), hcID(-1), theHC(nullptr), currentHit(nullptr) {
   edm::LogVerbatim("FiberSim") << "HFWedgeSD : Instantiated for " << iname;
 }

HFWedgeSD::~HFWedgeSD ( )

override

Definition at line 28 of file HFWedgeSD.cc.

References theHC.

28 { delete theHC; }

HFWedgeSD::theHC

HFShowerG4HitsCollection * theHC

Definition: HFWedgeSD.h:45

Member Function Documentation

void HFWedgeSD::clear ( void )

override

Definition at line 65 of file HFWedgeSD.cc.

Referenced by EndOfEvent().

65 {}

void HFWedgeSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 118 of file HFWedgeSD.cc.

References hitMap, and previousID.

Referenced by Initialize().

                           {
   hitMap.erase(hitMap.begin(), hitMap.end());
   previousID = -1;
 }

HFShowerG4Hit * HFWedgeSD::createNewHit ( )

protected

Definition at line 87 of file HFWedgeSD.cc.

References currentID, edep, globalPos, hitMap, localPos, MeV, momDir, previousID, HFShowerG4Hit::setGlobalPos(), HFShowerG4Hit::setHitId(), HFShowerG4Hit::setLocalPos(), HFShowerG4Hit::setPrimMomDir(), HFShowerG4Hit::setTime(), HFShowerG4Hit::setTrackId(), theHC, time, trackID, and updateHit().

Referenced by ProcessHits().

                                        {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("FiberSim") << "HFWedgeSD::CreateNewHit for ID " << currentID << " Track " << trackID
                                << " Edep: " << edep / CLHEP::MeV << " MeV; Time: " << time << " ns; Position (local) "
                                << localPos << " (global ) " << globalPos << " direction " << momDir;
 #endif
   HFShowerG4Hit* aHit = new HFShowerG4Hit;
   aHit->setHitId(currentID);
   aHit->setTrackId(trackID);
   aHit->setTime(time);
   aHit->setLocalPos(localPos);
   aHit->setGlobalPos(globalPos);
   aHit->setPrimMomDir(momDir);
   updateHit(aHit);
 
   theHC->insert(aHit);
   hitMap.insert(std::pair<int, HFShowerG4Hit*>(previousID, aHit));
 
   return aHit;
 }

void HFWedgeSD::DrawAll ( )

override

Definition at line 67 of file HFWedgeSD.cc.

67 {}

void HFWedgeSD::EndOfEvent ( G4HCofThisEvent * eventHC )

override

Definition at line 60 of file HFWedgeSD.cc.

References clear(), and theHC.

                                                {
   edm::LogVerbatim("FiberSim") << "HFWedgeSD: Sees" << theHC->entries() << " hits";
   clear();
 }

void HFWedgeSD::fillHits	(	edm::PCaloHitContainer &	,
		const std::string &
	)

overridevirtual

Implements SensitiveCaloDetector.

Definition at line 128 of file HFWedgeSD.cc.

128 {}

G4bool HFWedgeSD::hitExists ( )

protected

Definition at line 71 of file HFWedgeSD.cc.

References currentHit, currentID, hitMap, previousID, and updateHit().

Referenced by ProcessHits().

                             {
   // Update if in the same detector, time-slice and for same track
   if (currentID == previousID) {
     updateHit(currentHit);
     return true;
   }
 
   std::map<int, HFShowerG4Hit*>::const_iterator it = hitMap.find(currentID);
   if (it != hitMap.end()) {
     updateHit(currentHit);
     return true;
   }
 
   return false;
 }

void HFWedgeSD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 30 of file HFWedgeSD.cc.

References clearHits(), ecaldqm::collectionName, hcID, and theHC.

                                                {
   edm::LogVerbatim("FiberSim") << "HFWedgeSD : Initialize called for " << GetName() << " in collection " << HCE;
   theHC = new HFShowerG4HitsCollection(GetName(), collectionName[0]);
   if (hcID < 0)
     hcID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
   HCE->AddHitsCollection(hcID, theHC);
   edm::LogVerbatim("FiberSim") << "HFWedgeSD : Add hit collectrion for " << collectionName[0] << ":" << hcID << ":" << theHC;
 
   clearHits();
 }

void HFWedgeSD::PrintAll ( )

override

Definition at line 69 of file HFWedgeSD.cc.

69 {}

G4bool HFWedgeSD::ProcessHits	(	G4Step *	step,
		G4TouchableHistory *	tHistory
	)

overridevirtual

Implements SensitiveDetector.

Definition at line 41 of file HFWedgeSD.cc.

References createNewHit(), currentHit, currentID, edep, globalPos, hitExists(), localPos, momDir, setDetUnitId(), time, and trackID.

                                                                 {
   G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
   const G4VTouchable* touch = preStepPoint->GetTouchable();
   currentID = setDetUnitId(aStep);
   trackID = aStep->GetTrack()->GetTrackID();
   edep = aStep->GetTotalEnergyDeposit();
   time = (preStepPoint->GetGlobalTime()) / ns;
 
   globalPos = preStepPoint->GetPosition();
   localPos = touch->GetHistory()->GetTopTransform().TransformPoint(globalPos);
   const G4DynamicParticle* particle = aStep->GetTrack()->GetDynamicParticle();
   momDir = particle->GetMomentumDirection();
 
   if (hitExists() == false && edep > 0.)
     currentHit = createNewHit();
 
   return true;
 }

uint32_t HFWedgeSD::setDetUnitId ( const G4Step * aStep )

overridevirtual

Implements SensitiveDetector.

Definition at line 123 of file HFWedgeSD.cc.

Referenced by ProcessHits().

                                                     {
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   return (touch->GetReplicaNumber(0));
 }

void HFWedgeSD::updateHit ( HFShowerG4Hit * aHit )

protected

Definition at line 108 of file HFWedgeSD.cc.

References currentID, edep, MeV, previousID, and HFShowerG4Hit::updateEnergy().

Referenced by createNewHit(), and hitExists().

                                              {
   if (edep != 0) {
     aHit->updateEnergy(edep);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("FiberSim") << "HFWedgeSD: Add energy deposit in " << currentID << " edep " << edep / CLHEP::MeV << " MeV";
 #endif
   }
   previousID = currentID;
 }