#include <HGCSD.h>

Inheritance diagram for HGCSD:

Public Member Functions
	HGCSD (const std::string &, const HGCalTBDDDConstants , const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager )

uint32_t	setDetUnitId (const G4Step *step) override

	~HGCSD () override=default

Public Member Functions inherited from CaloSD
	CaloSD (const std::string &aSDname, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &p, const SimTrackManager *, float timeSlice=1., bool ignoreTkID=false, const std::string &newcolname="")

void	clear () override

void	clearHits () override

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

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

void	Initialize (G4HCofThisEvent *HCE) override

bool	isItFineCalo (const G4VTouchable *touch)

void	newCollection (const std::string &name, edm::ParameterSet const &p)

void	PrintAll () override

G4bool	ProcessHits (G4Step step, G4TouchableHistory ) override

bool	ProcessHits (G4GFlashSpot aSpot, G4TouchableHistory ) override

void	reset () override

	~CaloSD () override

Public Member Functions inherited from SensitiveCaloDetector
	SensitiveCaloDetector (const std::string &iname, const SensitiveDetectorCatalog &clg, const std::string &newcollname="")

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, const std::string &newcollname="")

	~SensitiveDetector () override

Public Member Functions inherited from Observer< const BeginOfRun *>
	Observer ()

void	slotForUpdate (const BeginOfRun * iT)

virtual	~Observer ()

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 EndOfTrack *>
	Observer ()

void	slotForUpdate (const EndOfTrack * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const EndOfEvent *>
	Observer ()

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

Public Member Functions inherited from Observer< const BeginOfJob *>
	Observer ()

void	slotForUpdate (const BeginOfJob * iT)

virtual	~Observer ()

Protected Member Functions
void	endEvent () override

bool	filterHit (CaloG4Hit *, double) override

double	getEnergyDeposit (const G4Step *) override

void	initEvent (const BeginOfEvent *) override

void	initRun () override

void	update (const BeginOfRun *) override

void	update (const BeginOfTrack *trk) override

void	update (const EndOfTrack *trk) override

void	update (const BeginOfEvent *) override

void	update (const ::EndOfEvent *) override

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

Protected Member Functions inherited from CaloSD
bool	checkHit (int k=0)

CaloG4Hit *	createNewHit (const G4Step , const G4Track , int k)

virtual double	EnergyCorrected (const G4Step &step, const G4Track *)

unsigned int	findBoundaryCrossingParent (const G4Track *track, bool markParentAsSaveable=true)

double	getAttenuation (const G4Step *aStep, double birk1, double birk2, double birk3) const

virtual uint16_t	getDepth (const G4Step *)

virtual bool	getFromLibrary (const G4Step *step)

int	getNumberOfHits (int k=0)

double	getResponseWt (const G4Track *, int k=0)

virtual int	getTrackID (const G4Track *)

bool	hitExists (const G4Step *, int k)

void	ignoreRejection ()

void	printDetectorLevels (const G4VTouchable *) const

void	processHit (const G4Step *step)

virtual void	processSecondHit (const G4Step , const G4Track )

void	resetForNewPrimary (const G4Step *)

void	setNumberCheckedHits (int val, int k=0)

void	setParameterized (bool val)

G4ThreeVector	setToGlobal (const G4ThreeVector &, const G4VTouchable *) const

G4ThreeVector	setToLocal (const G4ThreeVector &, const G4VTouchable *) const

virtual int	setTrackID (const G4Step *)

void	setUseMap (bool val)

std::string	shortreprID (const CaloHitID &ID)

std::string	shortreprID (const CaloG4Hit *hit)

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

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

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

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

void	update (const ::EndOfEvent *) override

void	updateHit (CaloG4Hit *, int k)

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

Protected Member Functions inherited from Observer< const EndOfEvent *>
virtual void	update (const EndOfEvent *)=0
	This routine will be called when the appropriate signal arrives. More...

Private Member Functions
bool	isItinFidVolume (const G4ThreeVector &)

uint32_t	setDetUnitId (ForwardSubdetector &, int, int, int, int, G4ThreeVector &)

Private Attributes
std::vector< double >	angles_

bool	dd4hep_

double	eminHit_

HGCalGeometryMode::GeometryMode	geom_mode_

const HGCalTBDDDConstants *	hgcons_

int	levelT_

std::unique_ptr< HGCMouseBite >	mouseBite_

double	mouseBiteCut_

ForwardSubdetector	myFwdSubdet_

std::string	nameX_

std::unique_ptr< HGCNumberingScheme >	numberingScheme_

bool	rejectMB_

double	slopeMin_

bool	storeAllG4Hits_

std::vector< double >	t_Angle_

std::vector< double >	t_dEStep1_

std::vector< double >	t_dEStep2_

uint32_t	t_EventID_

std::vector< int >	t_Layer_

std::vector< int >	t_Parcode_

std::vector< double >	t_TrackE_

bool	waferRot_

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

Static Protected Member Functions inherited from CaloSD
static std::string	printableDecayChain (const std::vector< unsigned int > &decayChain)

Protected Attributes inherited from CaloSD
std::string	collName_ [2]

CaloG4Hit *	currentHit [2]

CaloHitID	currentID [2]

std::string	detName_

float	edepositEM

float	edepositHAD

double	eminHit

double	energyCut

G4ThreeVector	entranceLocal

G4ThreeVector	entrancePoint

bool	forceSave

std::vector< std::string >	hcn_

float	incidentEnergy

double	kmaxIon

double	kmaxNeutron

double	kmaxProton

int	nHC_

G4ThreeVector	posGlobal

CaloHitID	previousID [2]

bool	suppressHeavy

double	tmaxHit

std::vector< int >	useResMap_

Detailed Description

Definition at line 26 of file HGCSD.h.

Constructor & Destructor Documentation

◆ HGCSD()

HGCSD::HGCSD	(	const std::string &	name,
		const HGCalTBDDDConstants *	hgc,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

Definition at line 31 of file HGCSD.cc.

References angles_, dd4hep_, CaloSD::eminHit, eminHit_, ForwardEmpty, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HGCEE, HGCHEB, HGCHEF, mouseBite_, mouseBiteCut_, myFwdSubdet_, mergeVDriftHistosByStation::name, nameX_, numberingScheme_, AlCaHLTBitMon_ParallelJobs::p, rejectMB_, CaloSD::setNumberCheckedHits(), CaloSD::setUseMap(), storeAllG4Hits_, funct::tan(), and waferRot_.

     : CaloSD(name,
              clg,
              p,
              manager,
              (float)(p.getParameter<edm::ParameterSet>("HGCSD").getParameter<double>("TimeSliceUnit")),
              p.getParameter<edm::ParameterSet>("HGCSD").getParameter<bool>("IgnoreTrackID")),
       hgcons_(hgc),
       slopeMin_(0),
       levelT_(99) {
 #ifdef plotDebug
   tree_ = nullptr;
 #endif
   numberingScheme_.reset(nullptr);
   mouseBite_.reset(nullptr);
 
   edm::ParameterSet m_HGC = p.getParameter<edm::ParameterSet>("HGCSD");
   eminHit_ = m_HGC.getParameter<double>("EminHit") * CLHEP::MeV;
   storeAllG4Hits_ = m_HGC.getParameter<bool>("StoreAllG4Hits");
   rejectMB_ = m_HGC.getParameter<bool>("RejectMouseBite");
   waferRot_ = m_HGC.getParameter<bool>("RotatedWafer");
   angles_ = m_HGC.getUntrackedParameter<std::vector<double>>("WaferAngles");
   double waferSize = m_HGC.getUntrackedParameter<double>("WaferSize") * CLHEP::mm;
   double mouseBite = m_HGC.getUntrackedParameter<double>("MouseBite") * CLHEP::mm;
   mouseBiteCut_ = waferSize * tan(30.0 * CLHEP::deg) - mouseBite;
   dd4hep_ = p.getParameter<bool>("g4GeometryDD4hepSource");
 
   if (storeAllG4Hits_) {
     setUseMap(false);
     setNumberCheckedHits(0);
   }
   //this is defined in the hgcsens.xml
   G4String myName = name;
   myFwdSubdet_ = ForwardSubdetector::ForwardEmpty;
   nameX_ = "HGCal";
   if (myName.find("HitsEE") != std::string::npos) {
     myFwdSubdet_ = ForwardSubdetector::HGCEE;
     nameX_ = "HGCalEESensitive";
   } else if (myName.find("HitsHEfront") != std::string::npos) {
     myFwdSubdet_ = ForwardSubdetector::HGCHEF;
     nameX_ = "HGCalHESiliconSensitive";
   } else if (myName.find("HitsHEback") != std::string::npos) {
     myFwdSubdet_ = ForwardSubdetector::HGCHEB;
     nameX_ = "HGCalHEScintillatorSensitive";
   }
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "**************************************************"
                              << "\n"
                              << "*                                                *"
                              << "\n"
                              << "* Constructing a HGCSD  with name " << name << "\n"
                              << "*                                                *"
                              << "\n"
                              << "**************************************************";
 #endif
   edm::LogVerbatim("HGCSim") << "HGCSD:: Threshold for storing hits: " << eminHit << " for " << nameX_ << " subdet "
                              << myFwdSubdet_;
   edm::LogVerbatim("HGCSim") << "Flag for storing individual Geant4 Hits " << storeAllG4Hits_;
   edm::LogVerbatim("HGCSim") << "Reject MosueBite Flag: " << rejectMB_ << " Size of wafer " << waferSize
                              << " Mouse Bite " << mouseBite << ":" << mouseBiteCut_ << " along " << angles_.size()
                              << " axes";
 }

◆ ~HGCSD()

HGCSD::~HGCSD ( )

overridedefault

Member Function Documentation

◆ endEvent()

void HGCSD::endEvent ( )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 255 of file HGCSD.cc.

                      {
 #ifdef plotDebug
   if (tree_)
     tree_->Fill();
 #endif
 }

◆ filterHit()

bool HGCSD::filterHit	(	CaloG4Hit *	aHit,
		double	time
	)

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 262 of file HGCSD.cc.

References eminHit_, CaloG4Hit::getEnergyDeposit(), hcalRecHitTable_cff::time, and CaloSD::tmaxHit.

                                                   {
   return ((time <= tmaxHit) && (aHit->getEnergyDeposit() > eminHit_));
 }

◆ getEnergyDeposit()

double HGCSD::getEnergyDeposit ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 99 of file HGCSD.cc.

References funct::abs(), angle(), CaloSD::getResponseWt(), nano_mu_digi_cff::layer, levelT_, alignCSCRings::r, slopeMin_, t_Angle_, t_dEStep1_, t_dEStep2_, t_Layer_, t_Parcode_, t_TrackE_, and z.

                                                   {
   double r = aStep->GetPreStepPoint()->GetPosition().perp();
   double z = std::abs(aStep->GetPreStepPoint()->GetPosition().z());
 
 #ifdef EDM_ML_DEBUG
   G4int parCode = aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
   G4LogicalVolume* lv = aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume();
   edm::LogVerbatim("HGCSim") << "HGCSD: Hit from standard path from " << lv->GetName() << " for Track "
                              << aStep->GetTrack()->GetTrackID() << " ("
                              << aStep->GetTrack()->GetDefinition()->GetParticleName() << ":" << parCode << ") R = " << r
                              << " Z = " << z << " slope = " << r / z << ":" << slopeMin_;
 #endif
 
   // Apply fiductial volume
   if (r < z * slopeMin_) {
     return 0.0;
   }
 
   double wt1 = getResponseWt(aStep->GetTrack());
   double wt2 = aStep->GetTrack()->GetWeight();
   double destep = wt1 * aStep->GetTotalEnergyDeposit();
   if (wt2 > 0)
     destep *= wt2;
 
 #ifdef plotDebug
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   G4double tmptrackE = aStep->GetTrack()->GetKineticEnergy();
   G4int parCodex = aStep->GetTrack()->GetDefinition()->GetPDGEncoding();
   G4double angle = (aStep->GetTrack()->GetMomentumDirection().theta()) / CLHEP::deg;
   G4int layer = ((touch->GetHistoryDepth() == levelT_) ? touch->GetReplicaNumber(0) : touch->GetReplicaNumber(2));
   G4int ilayer = (layer - 1) / 3;
   if (aStep->GetTotalEnergyDeposit() > 0) {
     t_Layer_.emplace_back(ilayer);
     t_Parcode_.emplace_back(parCodex);
     t_dEStep1_.emplace_back(aStep->GetTotalEnergyDeposit());
     t_dEStep2_.emplace_back(destep);
     t_TrackE_.emplace_back(tmptrackE);
     t_Angle_.emplace_back(angle);
   }
 #endif
 
   return destep;
 }

◆ initEvent()

void HGCSD::initEvent ( const BeginOfEvent * g4Event )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 242 of file HGCSD.cc.

References t_Angle_, t_dEStep1_, t_dEStep2_, t_EventID_, t_Layer_, t_Parcode_, and t_TrackE_.

                                                  {
   const G4Event* evt = (*g4Event)();
   t_EventID_ = evt->GetEventID();
 #ifdef plotDebug
   t_Layer_.clear();
   t_Parcode_.clear();
   t_dEStep1_.clear();
   t_dEStep2_.clear();
   t_TrackE_.clear();
   t_Angle_.clear();
 #endif
 }

◆ initRun()

void HGCSD::initRun ( )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 226 of file HGCSD.cc.

References t_Angle_, t_dEStep1_, t_dEStep2_, t_EventID_, t_Layer_, t_Parcode_, t_TrackE_, and compare::tfile.

                     {
 #ifdef plotDebug
   edm::Service<TFileService> tfile;
   if (tfile.isAvailable()) {
     tree_ = tfile->make<TTree>("TreeHGCSD", "TreeHGCSD");
     tree_->Branch("EventID", &t_EventID_);
     tree_->Branch("Layer", &t_Layer_);
     tree_->Branch("ParticleCode", &t_Parcode_);
     tree_->Branch("dEStepOriginal", &t_dEStep1_);
     tree_->Branch("dEStepWeighted", &t_dEStep2_);
     tree_->Branch("TrackEnergy", &t_TrackE_);
     tree_->Branch("ThetaAngle", &t_Angle_);
   }
 #endif
 }

◆ isItinFidVolume()

bool HGCSD::isItinFidVolume ( const G4ThreeVector & )

inlineprivate

Definition at line 48 of file HGCSD.h.

48 { return true; }

◆ setDetUnitId() [1/2]

uint32_t HGCSD::setDetUnitId ( const G4Step * step )

overridevirtual

Implements CaloSD.

Definition at line 143 of file HGCSD.cc.

References TauDecayModes::dec, EcalPhiSymFlatTableProducers_cfi::id, isotrackApplyRegressor::k, nano_mu_digi_cff::layer, levelT_, DTRecHitClients_cfi::local, callgraph::module, mouseBite_, myFwdSubdet_, rejectMB_, HGCalTestNumbering::unpackHexagonIndex(), and z.

                                                 {
   const G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
   const G4VTouchable* touch = preStepPoint->GetTouchable();
 
   //determine the exact position in global coordinates in the mass geometry
   G4ThreeVector hitPoint = preStepPoint->GetPosition();
   float globalZ = touch->GetTranslation(0).z();
   int iz(globalZ > 0 ? 1 : -1);
 
   //convert to local coordinates (=local to the current volume):
   G4ThreeVector localpos = touch->GetHistory()->GetTopTransform().TransformPoint(hitPoint);
 
   //get the det unit id with
   ForwardSubdetector subdet = myFwdSubdet_;
 
   int layer(-1), moduleLev(-1), module(-1), cell(-1);
   if (touch->GetHistoryDepth() == levelT_) {
     layer = touch->GetReplicaNumber(0);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "Depths: " << touch->GetHistoryDepth() << " name " << touch->GetVolume(0)->GetName()
                                << " layer:module:cell " << layer << ":" << moduleLev << ":" << module << ":" << cell;
 #endif
   } else {
     layer = touch->GetReplicaNumber(2);
     module = touch->GetReplicaNumber(1);
     cell = touch->GetReplicaNumber(0);
     moduleLev = 1;
   }
 #ifdef EDM_ML_DEBUG
   const G4Material* mat = aStep->GetPreStepPoint()->GetMaterial();
   edm::LogVerbatim("HGCSim") << "Depths: " << touch->GetHistoryDepth() << " name " << touch->GetVolume(0)->GetName()
                              << ":" << touch->GetReplicaNumber(0) << "   " << touch->GetVolume(1)->GetName() << ":"
                              << touch->GetReplicaNumber(1) << "   " << touch->GetVolume(2)->GetName() << ":"
                              << touch->GetReplicaNumber(2) << "    layer:module:cell " << layer << ":" << moduleLev
                              << ":" << module << ":" << cell << " Material " << mat->GetName() << ":"
                              << mat->GetRadlen();
   for (int k = 0; k < touch->GetHistoryDepth(); ++k)
     edm::LogVerbatim("HGCSim") << "Level [" << k << "] " << touch->GetVolume(k)->GetName() << ":"
                                << touch->GetReplicaNumber(k);
 #endif
   // The following statement should be examined later before elimination
   // VI: this is likely a check if media is vacuum - not needed
   if (aStep->GetPreStepPoint()->GetMaterial()->GetRadlen() > 100000.)
     return 0;
 
   uint32_t id = setDetUnitId(subdet, layer, module, cell, iz, localpos);
   if (rejectMB_ && id != 0) {
     int det, z, lay, wafer, type, ic;
     HGCalTestNumbering::unpackHexagonIndex(id, det, z, lay, wafer, type, ic);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCSim") << "ID " << std::hex << id << std::dec << " Input " << subdet << ":" << layer << ":"
                                << module << ":" << cell << ":" << iz << localpos.x() << ":" << localpos.y()
                                << " Decode " << det << ":" << z << ":" << lay << ":" << wafer << ":" << type << ":"
                                << ic;
 #endif
     G4ThreeVector local =
         ((moduleLev >= 0) ? (touch->GetHistory()->GetTransform(moduleLev).TransformPoint(hitPoint)) : G4ThreeVector());
     if (mouseBite_->exclude(local, z, layer, wafer, 0))
       id = 0;
   }
   return id;
 }

◆ setDetUnitId() [2/2]

uint32_t HGCSD::setDetUnitId	(	ForwardSubdetector &	subdet,
		int	layer,
		int	module,
		int	cell,
		int	iz,
		G4ThreeVector &	pos
	)

private

Definition at line 266 of file HGCSD.cc.

References EcalPhiSymFlatTableProducers_cfi::id, nano_mu_digi_cff::layer, and numberingScheme_.

                                                                                                                     {
   uint32_t id = numberingScheme_ ? numberingScheme_->getUnitID(subdet, layer, module, cell, iz, pos) : 0;
   return id;
 }

◆ update() [1/6]

void CaloSD::update

overrideprotected

Definition at line 781 of file CaloSD.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                      {
   endEvent();
   for (int k = 0; k < nHC_; ++k) {
     slave[k].get()->ReserveMemory(theHC[k]->entries());
 
     int count(0);
     double eEM(0.0);
     double eHAD(0.0);
     double eEM2(0.0);
     double eHAD2(0.0);
 #ifdef EDM_ML_DEBUG
     int wrong(0);
     double tt(0.0);
     double zloc(0.0);
     double zglob(0.0);
     double ee(0.0);
 #endif
     int hc_entries = theHC[k]->entries();
     for (int i = 0; i < hc_entries; ++i) {
 #ifdef EDM_ML_DEBUG
       if (!saveHit((*theHC[k])[i], k)) {
         ++wrong;
       }
 #else
       saveHit((*theHC[k])[i], k);
 #endif
 
       ++count;
       double x = (*theHC[k])[i]->getEM();
       eEM += x;
       eEM2 += x * x;
       x = (*theHC[k])[i]->getHadr();
       eHAD += x;
       eHAD2 += x * x;
 #ifdef EDM_ML_DEBUG
       tt += (*theHC[k])[i]->getTimeSlice();
       ee += (*theHC[k])[i]->getIncidentEnergy();
       zglob += std::abs((*theHC[k])[i]->getEntry().z());
       zloc += std::abs((*theHC[k])[i]->getEntryLocal().z());
 #endif
     }
 
     double norm = (count > 0) ? 1.0 / count : 0.0;
     eEM *= norm;
     eEM2 *= norm;
     eHAD *= norm;
     eHAD2 *= norm;
     eEM2 = std::sqrt(eEM2 - eEM * eEM);
     eHAD2 = std::sqrt(eHAD2 - eHAD * eHAD);
 #ifdef EDM_ML_DEBUG
     tt *= norm;
     ee *= norm;
     zglob *= norm;
     zloc *= norm;
     edm::LogVerbatim("CaloSim") << "CaloSD: " << GetName() << " store " << count << " hits; " << wrong
                                 << " track IDs not given properly and " << totalHits - count
                                 << " hits not passing cuts\n EmeanEM= " << eEM << " ErmsEM= " << eEM2
                                 << "\n EmeanHAD= " << eHAD << " ErmsHAD= " << eHAD2 << " TimeMean= " << tt
                                 << " E0mean= " << ee << " Zglob= " << zglob << " Zloc= " << zloc << " ";
 #endif
     std::vector<std::unique_ptr<CaloG4Hit>>().swap(reusehit[k]);
     if (useMap)
       hitMap[k].erase(hitMap[k].begin(), hitMap[k].end());
   }
   tkMap.erase(tkMap.begin(), tkMap.end());
   boundaryCrossingParentMap_.clear();
 }

◆ update() [2/6]

void CaloSD::update

overrideprotected

Definition at line 1014 of file CaloSD.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                            {
   int primary = -1;
   TrackInformation* trkInfo = cmsTrackInformation((*trk)());
   if (trkInfo->isPrimary())
     primary = (*trk)()->GetTrackID();
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "New track: isPrimary " << trkInfo->isPrimary() << " primary ID = " << primary
                               << " primary ancestor ID " << primAncestor;
 #endif
 
   // update the information if a different primary track ID
 
   if (primary > 0 && primary != primAncestor) {
     primAncestor = primary;
 
     // clean the hits information
 
     bool clean(false);
     for (int k = 0; k < nHC_; ++k)
       if (theHC[k]->entries() > 0)
         clean = true;
     if (clean)
       cleanHitCollection();
   }
 }

◆ update() [3/6]

void CaloSD::update

overrideprotected

Definition at line 750 of file CaloSD.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

750 { initRun(); }

HGCSD::initRun

void initRun() override

Definition: HGCSD.cc:226

◆ update() [4/6]

void CaloSD::update

overrideprotected

Definition at line 752 of file CaloSD.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                                {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: Dispatched BeginOfEvent for " << GetName() << " !";
 #endif
   clearHits();
   initEvent(g4Event);
 }

◆ update() [5/6]

void CaloSD::update

overrideprotected

Definition at line 760 of file CaloSD.cc.

Referenced by MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), and MatrixUtil.Steps::overwrite().

                                          {
   int id = (*trk)()->GetTrackID();
   TrackInformation* trkI = cmsTrackInformation((*trk)());
   int lastTrackID = -1;
   if (trkI)
     lastTrackID = trkI->getIDonCaloSurface();
   if (id == lastTrackID) {
     auto trksForThisEvent = m_trackManager->trackContainer();
     if (!trksForThisEvent->empty()) {
       TrackWithHistory* trkH = trksForThisEvent->back();
       if (trkH->trackID() == id) {
         tkMap[id] = trkH;
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("CaloSim") << "CaloSD: get track " << id << " from Container of size "
                                     << trksForThisEvent->size() << " with ID " << trkH->trackID();
 #endif
       }
     }
   }
 }

◆ update() [6/6]

void HGCSD::update ( const BeginOfJob * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob *>.

Definition at line 206 of file HGCSD.cc.

References angles_, dd4hep_, Exception, geom_mode_, HGCalTBDDDConstants::geomMode(), hgcons_, levelT_, HGCalTBDDDConstants::levelTop(), HGCalTBDDDConstants::minSlope(), mouseBite_, mouseBiteCut_, nameX_, numberingScheme_, rejectMB_, slopeMin_, and waferRot_.

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

                                         {
   if (hgcons_ != nullptr) {
     geom_mode_ = hgcons_->geomMode();
     slopeMin_ = hgcons_->minSlope();
     levelT_ = hgcons_->levelTop();
     if (dd4hep_)
       ++levelT_;
     numberingScheme_ = std::make_unique<HGCNumberingScheme>(*hgcons_, nameX_);
     if (rejectMB_)
       mouseBite_ = std::make_unique<HGCMouseBite>(*hgcons_, angles_, mouseBiteCut_, waferRot_);
   } else {
     edm::LogError("HGCSim") << "HGCSD : Cannot find HGCalTBDDDConstants for " << nameX_;
     throw cms::Exception("Unknown", "HGCSD") << "Cannot find HGCalTBDDDConstants for " << nameX_ << "\n";
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCSim") << "HGCSD::Initialized with mode " << geom_mode_ << " Slope cut " << slopeMin_
                              << " top Level " << levelT_ << " dd4hep " << dd4hep_;
 #endif
 }

Member Data Documentation

◆ angles_

std::vector<double> HGCSD::angles_

private

Definition at line 62 of file HGCSD.h.

Referenced by HGCSD(), and update().

◆ dd4hep_

bool HGCSD::dd4hep_

private

Definition at line 61 of file HGCSD.h.

Referenced by HGCSD(), and update().

◆ eminHit_

double HGCSD::eminHit_

private

Definition at line 55 of file HGCSD.h.

Referenced by filterHit(), and HGCSD().

◆ geom_mode_

HGCalGeometryMode::GeometryMode HGCSD::geom_mode_

private

Definition at line 52 of file HGCSD.h.

Referenced by update().

◆ hgcons_

const HGCalTBDDDConstants* HGCSD::hgcons_

private

Definition at line 50 of file HGCSD.h.

Referenced by update().

◆ levelT_

int HGCSD::levelT_

private

Definition at line 58 of file HGCSD.h.

Referenced by getEnergyDeposit(), setDetUnitId(), and update().

◆ mouseBite_

std::unique_ptr<HGCMouseBite> HGCSD::mouseBite_

private

Definition at line 54 of file HGCSD.h.

Referenced by HGCSD(), setDetUnitId(), and update().

◆ mouseBiteCut_

double HGCSD::mouseBiteCut_

private

Definition at line 60 of file HGCSD.h.

Referenced by HGCSD(), and update().

◆ myFwdSubdet_

ForwardSubdetector HGCSD::myFwdSubdet_

private

Definition at line 56 of file HGCSD.h.

Referenced by HGCSD(), and setDetUnitId().

◆ nameX_

std::string HGCSD::nameX_

private

Definition at line 51 of file HGCSD.h.

Referenced by HGCSD(), and update().

◆ numberingScheme_

std::unique_ptr<HGCNumberingScheme> HGCSD::numberingScheme_

private

Definition at line 53 of file HGCSD.h.

Referenced by HGCSD(), setDetUnitId(), and update().

◆ rejectMB_

bool HGCSD::rejectMB_

private

Definition at line 59 of file HGCSD.h.

Referenced by HGCSD(), setDetUnitId(), and update().

◆ slopeMin_

double HGCSD::slopeMin_

private

Definition at line 57 of file HGCSD.h.

Referenced by getEnergyDeposit(), and update().

◆ storeAllG4Hits_

bool HGCSD::storeAllG4Hits_

private

Definition at line 59 of file HGCSD.h.

Referenced by HGCSD().

◆ t_Angle_

std::vector<double> HGCSD::t_Angle_

private

Definition at line 70 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ t_dEStep1_

std::vector<double> HGCSD::t_dEStep1_

private

Definition at line 69 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ t_dEStep2_

std::vector<double> HGCSD::t_dEStep2_

private

Definition at line 69 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ t_EventID_

uint32_t HGCSD::t_EventID_

private

Definition at line 67 of file HGCSD.h.

Referenced by initEvent(), and initRun().

◆ t_Layer_

std::vector<int> HGCSD::t_Layer_

private

Definition at line 68 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ t_Parcode_

std::vector<int> HGCSD::t_Parcode_

private

Definition at line 68 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ t_TrackE_

std::vector<double> HGCSD::t_TrackE_

private

Definition at line 69 of file HGCSD.h.

Referenced by getEnergyDeposit(), initEvent(), and initRun().

◆ waferRot_

bool HGCSD::waferRot_

private

Definition at line 59 of file HGCSD.h.

Referenced by HGCSD(), and update().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HGCSD()

◆ ~HGCSD()

Member Function Documentation

◆ endEvent()

◆ filterHit()

◆ getEnergyDeposit()

◆ initEvent()

◆ initRun()

◆ isItinFidVolume()

◆ setDetUnitId() [1/2]

◆ setDetUnitId() [2/2]

◆ update() [1/6]

◆ update() [2/6]

◆ update() [3/6]

◆ update() [4/6]

◆ update() [5/6]

◆ update() [6/6]

Member Data Documentation

◆ angles_

◆ dd4hep_

◆ eminHit_

◆ geom_mode_

◆ hgcons_

◆ levelT_

◆ mouseBite_

◆ mouseBiteCut_

◆ myFwdSubdet_

◆ nameX_

◆ numberingScheme_

◆ rejectMB_

◆ slopeMin_

◆ storeAllG4Hits_

◆ t_Angle_

◆ t_dEStep1_

◆ t_dEStep2_

◆ t_EventID_

◆ t_Layer_

◆ t_Parcode_

◆ t_TrackE_

◆ waferRot_