HFWedgeSD Class Reference

#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::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::~HFWedgeSD ( )

override

Definition at line 26 of file HFWedgeSD.cc.

References theHC.

{ delete theHC; }

Member Function Documentation

clear()

void HFWedgeSD::clear ( void  )

override

Definition at line 64 of file HFWedgeSD.cc.

Referenced by EndOfEvent().

{}

clearHits()

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

createNewHit()

HFShowerG4Hit * HFWedgeSD::createNewHit ( )

protected

Definition at line 86 of file HFWedgeSD.cc.

References currentID, edep, globalPos, hitMap, localPos, 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;
}

DrawAll()

void HFWedgeSD::DrawAll ( )

override

Definition at line 66 of file HFWedgeSD.cc.

{}

EndOfEvent()

void HFWedgeSD::EndOfEvent ( G4HCofThisEvent *  eventHC )

override

Definition at line 59 of file HFWedgeSD.cc.

References clear(), and theHC.

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

fillHits()

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

overridevirtual

Implements SensitiveCaloDetector.

Definition at line 128 of file HFWedgeSD.cc.

{}

hitExists()

G4bool HFWedgeSD::hitExists ( )

protected

Definition at line 70 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;
}

Initialize()

void HFWedgeSD::Initialize ( G4HCofThisEvent *  HCE )

override

Definition at line 28 of file HFWedgeSD.cc.

References clearHits(), bysipixelclustmulteventfilter_cfi::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();
}

PrintAll()

void HFWedgeSD::PrintAll ( )

override

Definition at line 68 of file HFWedgeSD.cc.

{}

ProcessHits()

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

overridevirtual

Implements SensitiveDetector.

Definition at line 40 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;
}

setDetUnitId()

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));
}

updateHit()

void HFWedgeSD::updateHit ( HFShowerG4Hit *  aHit )

protected

Definition at line 107 of file HFWedgeSD.cc.

References currentID, edep, 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;
}

Member Data Documentation

currentHit

HFShowerG4Hit* HFWedgeSD::currentHit

private

Definition at line 49 of file HFWedgeSD.h.

Referenced by hitExists(), and ProcessHits().

currentID

int HFWedgeSD::currentID

private

Definition at line 46 of file HFWedgeSD.h.

Referenced by createNewHit(), hitExists(), ProcessHits(), and updateHit().

edep

double HFWedgeSD::edep

private

Definition at line 47 of file HFWedgeSD.h.

Referenced by createNewHit(), ProcessHits(), and updateHit().

globalPos

G4ThreeVector HFWedgeSD::globalPos

private

Definition at line 48 of file HFWedgeSD.h.

Referenced by createNewHit(), and ProcessHits().

hcID

int HFWedgeSD::hcID

private

Definition at line 42 of file HFWedgeSD.h.

Referenced by Initialize().

hitMap

std::map<int, HFShowerG4Hit*> HFWedgeSD::hitMap

private

Definition at line 44 of file HFWedgeSD.h.

Referenced by clearHits(), createNewHit(), and hitExists().

localPos

G4ThreeVector HFWedgeSD::localPos

private

Definition at line 48 of file HFWedgeSD.h.

Referenced by createNewHit(), and ProcessHits().

momDir

G4ThreeVector HFWedgeSD::momDir

private

Definition at line 48 of file HFWedgeSD.h.

Referenced by createNewHit(), and ProcessHits().

previousID

int HFWedgeSD::previousID

private

Definition at line 46 of file HFWedgeSD.h.

Referenced by clearHits(), createNewHit(), hitExists(), and updateHit().

theHC

HFShowerG4HitsCollection* HFWedgeSD::theHC

private

Definition at line 43 of file HFWedgeSD.h.

Referenced by createNewHit(), EndOfEvent(), Initialize(), and ~HFWedgeSD().

time

double HFWedgeSD::time

private

Definition at line 47 of file HFWedgeSD.h.

Referenced by createNewHit(), and ProcessHits().

trackID

int HFWedgeSD::trackID

private

Definition at line 46 of file HFWedgeSD.h.

Referenced by createNewHit(), and ProcessHits().