#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 DDCompactView &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
	SensitiveCaloDetector (const std::string &iname, const DDCompactView &cpv, 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

	SensitiveDetector (const std::string &iname, const DDCompactView &cpv, const SensitiveDetectorCatalog &, edm::ParameterSet const &p)

	~SensitiveDetector () override

Protected Member Functions
HFShowerG4Hit *	createNewHit ()

G4bool	hitExists ()

void	updateHit (HFShowerG4Hit *)

Protected Member Functions inherited from SensitiveDetector
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 19 of file HFWedgeSD.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 19 of file HFWedgeSD.cc.

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

HFWedgeSD::~HFWedgeSD ( )

override

Definition at line 27 of file HFWedgeSD.cc.

References theHC.

                       {
   delete theHC;
 }

Member Function Documentation

void HFWedgeSD::clear ( void )

override

Definition at line 68 of file HFWedgeSD.cc.

Referenced by EndOfEvent().

68 {}

void HFWedgeSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 123 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 91 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 70 of file HFWedgeSD.cc.

70 {}

void HFWedgeSD::EndOfEvent ( G4HCofThisEvent * eventHC )

override

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

133 {}

G4bool HFWedgeSD::hitExists ( )

protected

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

References clearHits(), ecaldqm::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 72 of file HFWedgeSD.cc.

72 {}

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 128 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 113 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;
 }