#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 edm::EventSetup &cpv, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, 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 edm::EventSetup &es, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p)

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 edm::EventSetup &es, const SensitiveDetectorCatalog &, edm::ParameterSet const &p, 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

const SimTrackManager *	m_trackManager

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 18 of file HFWedgeSD.h.

Constructor & Destructor Documentation

HFWedgeSD::HFWedgeSD	(	const std::string &	iname,
		const edm::EventSetup &	cpv,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

explicit

Definition at line 19 of file HFWedgeSD.cc.

     : SensitiveCaloDetector(iname, es, clg, p),
       m_trackManager(manager),
       hcID(-1),
       theHC(nullptr),
       currentHit(nullptr) {}

HFWedgeSD::~HFWedgeSD ( )

override

Definition at line 30 of file HFWedgeSD.cc.

References theHC.

30 { delete theHC; }

HFWedgeSD::theHC

HFShowerG4HitsCollection * theHC

Definition: HFWedgeSD.h:47

Member Function Documentation

void HFWedgeSD::clear ( void )

override

Definition at line 66 of file HFWedgeSD.cc.

Referenced by EndOfEvent().

66 {}

void HFWedgeSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 115 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 88 of file HFWedgeSD.cc.

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

Referenced by ProcessHits().

                                        {
   LogDebug("FiberSim") << "HFWedgeSD::CreateNewHit for ID " << currentID << " Track " << trackID
                        << " Edep: " << edep / MeV << " MeV; Time: " << time << " ns; Position (local) " << localPos
                        << " (global ) " << globalPos << " direction " << momDir;
   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 68 of file HFWedgeSD.cc.

68 {}

void HFWedgeSD::EndOfEvent ( G4HCofThisEvent * eventHC )

override

Definition at line 61 of file HFWedgeSD.cc.

References clear(), LogDebug, and theHC.

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

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

overridevirtual

Implements SensitiveCaloDetector.

Definition at line 125 of file HFWedgeSD.cc.

125 {}

G4bool HFWedgeSD::hitExists ( )

protected

Definition at line 72 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 32 of file HFWedgeSD.cc.

References clearHits(), bysipixelclustmulteventfilter_cfi::collectionName, hcID, LogDebug, and theHC.

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

void HFWedgeSD::PrintAll ( )

override

Definition at line 70 of file HFWedgeSD.cc.

70 {}

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

overridevirtual

Implements SensitiveDetector.

Definition at line 42 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 120 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 107 of file HFWedgeSD.cc.

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

Referenced by createNewHit(), and hitExists().

                                              {
   if (edep != 0) {
     aHit->updateEnergy(edep);
     LogDebug("FiberSim") << "HFWedgeSD: Add energy deposit in " << currentID << " edep " << edep / MeV << " MeV";
   }
   previousID = currentID;
 }