#include <CaloSD.h>

Inheritance diagram for CaloSD:

Classes
struct	Detector

Public Member Functions
	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

uint32_t	setDetUnitId (const G4Step *step) override=0

	~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 ()

Protected Member Functions
bool	checkHit (int k=0)

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

virtual void	endEvent ()

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

virtual bool	filterHit (CaloG4Hit *, double)

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 double	getEnergyDeposit (const G4Step *step)

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

virtual void	initEvent (const BeginOfEvent *)

virtual void	initRun ()

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...

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

Protected Attributes
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_

Private Member Functions
void	cleanHitCollection ()

bool	saveHit (CaloG4Hit *, int k=0)

void	storeHit (CaloG4Hit *, int k=0)

Private Attributes
std::unordered_map< unsigned int, unsigned int >	boundaryCrossingParentMap_

int	cleanIndex [2]

double	correctT

bool	corrTOFBeam

bool	doFineCalo_

bool	doFineCaloThisStep_

double	eMinFine_

double	eminHitD

std::vector< Detector >	fineDetectors_

int	hcID [2]

std::map< CaloHitID, CaloG4Hit * >	hitMap [2]

bool	ignoreReject

bool	ignoreTrackID

bool	isParameterized

const SimTrackManager *	m_trackManager

std::unique_ptr< CaloMeanResponse >	meanResponse [2]

int	nCheckedHits [2]

int	primAncestor

int	primIDSaved [2]

std::vector< std::unique_ptr< CaloG4Hit > >	reusehit [2]

std::unique_ptr< CaloSlaveSD >	slave [2]

CaloG4HitCollection *	theHC [2]

float	timeSlice

std::map< int, TrackWithHistory * >	tkMap

int	totalHits [2]

bool	useMap

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

Detailed Description

Definition at line 40 of file CaloSD.h.

Constructor & Destructor Documentation

◆ CaloSD()

CaloSD::CaloSD	(	const std::string &	aSDname,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager,
		float	timeSlice = `1.`,
		bool	ignoreTkID = `false`,
		const std::string &	newcolname = `""`
	)

Definition at line 33 of file CaloSD.cc.

References cleanIndex, collName_, correctT, corrTOFBeam, currentHit, currentID, hgcalTestNeighbor_cfi::detector, detName_, doFineCalo_, edepositEM, edepositHAD, eMinFine_, eminHit, eminHitD, energyCut, entranceLocal, entrancePoint, fineDetectors_, forceSave, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcID, hcn_, mps_fire::i, ignoreTrackID, incidentEnergy, isParameterized, dqmdumpme::k, kmaxIon, kmaxNeutron, kmaxProton, personalPlayback::level, meanResponse, Skims_PA_cff::name, nCheckedHits, nHC_, AlCaHLTBitMon_ParallelJobs::p, posGlobal, previousID, primAncestor, primIDSaved, reset(), slave, suppressHeavy, theHC, timeSlice, tmaxHit, totalHits, useMap, and useResMap_.

     : SensitiveCaloDetector(name, clg, newcollName),
       G4VGFlashSensitiveDetector(),
       eminHit(0.),
       m_trackManager(manager),
       ignoreTrackID(ignoreTkID),
       timeSlice(timeSliceUnit),
       eminHitD(0.) {
   //Parameters
   nHC_ = (newcollName.empty()) ? 1 : 2;
   for (int k = 0; k < 2; ++k) {
     currentHit[k] = nullptr;
     theHC[k] = nullptr;
     hcID[k] = -1;
   }
   bool dd4hep = p.getParameter<bool>("g4GeometryDD4hepSource");
   int addlevel = dd4hep ? 1 : 0;
   edm::ParameterSet m_CaloSD = p.getParameter<edm::ParameterSet>("CaloSD");
   energyCut = m_CaloSD.getParameter<double>("EminTrack") * CLHEP::GeV;
   tmaxHit = m_CaloSD.getParameter<double>("TmaxHit") * CLHEP::ns;
   std::vector<double> eminHits = m_CaloSD.getParameter<std::vector<double>>("EminHits");
   std::vector<double> tmaxHits = m_CaloSD.getParameter<std::vector<double>>("TmaxHits");
   hcn_ = m_CaloSD.getParameter<std::vector<std::string>>("HCNames");
   useResMap_ = m_CaloSD.getParameter<std::vector<int>>("UseResponseTables");
   std::vector<double> eminHitX = m_CaloSD.getParameter<std::vector<double>>("EminHitsDepth");
   suppressHeavy = m_CaloSD.getParameter<bool>("SuppressHeavy");
   kmaxIon = m_CaloSD.getParameter<double>("IonThreshold") * CLHEP::MeV;
   kmaxProton = m_CaloSD.getParameter<double>("ProtonThreshold") * CLHEP::MeV;
   kmaxNeutron = m_CaloSD.getParameter<double>("NeutronThreshold") * CLHEP::MeV;
   nCheckedHits[0] = m_CaloSD.getUntrackedParameter<int>("CheckHits", 25);
   nCheckedHits[1] = nCheckedHits[0];
   useMap = m_CaloSD.getUntrackedParameter<bool>("UseMap", true);
   int verbn = m_CaloSD.getUntrackedParameter<int>("Verbosity", 0);
   corrTOFBeam = m_CaloSD.getParameter<bool>("CorrectTOFBeam");
   double beamZ = m_CaloSD.getParameter<double>("BeamPosition") * CLHEP::cm;
   correctT = beamZ / CLHEP::c_light / CLHEP::nanosecond;
   doFineCalo_ = m_CaloSD.getParameter<bool>("DoFineCalo");
   eMinFine_ = m_CaloSD.getParameter<double>("EminFineTrack") * CLHEP::MeV;
   std::vector<std::string> fineNames = m_CaloSD.getParameter<std::vector<std::string>>("FineCaloNames");
   std::vector<int> fineLevels = m_CaloSD.getParameter<std::vector<int>>("FineCaloLevels");
   std::vector<int> useFines = m_CaloSD.getParameter<std::vector<int>>("UseFineCalo");
   for (auto& level : fineLevels)
     level += addlevel;
 
   SetVerboseLevel(verbn);
   for (int k = 0; k < 2; ++k)
     meanResponse[k].reset(nullptr);
   for (unsigned int k = 0; k < hcn_.size(); ++k) {
     if (name == hcn_[k]) {
       if (k < eminHits.size())
         eminHit = eminHits[k] * CLHEP::MeV;
       if (k < eminHitX.size())
         eminHitD = eminHitX[k] * CLHEP::MeV;
       if (k < tmaxHits.size())
         tmaxHit = tmaxHits[k] * CLHEP::ns;
       if (k < useResMap_.size() && useResMap_[k] > 0) {
         meanResponse[0] = std::make_unique<CaloMeanResponse>(p);
         break;
       }
     }
   }
   detName_ = name;
   collName_[0] = name;
   collName_[1] = newcollName;
   slave[0] = std::make_unique<CaloSlaveSD>(name);
   slave[1].reset(nullptr);
 
   for (int k = 0; k < 2; ++k) {
     currentID[k] = CaloHitID(timeSlice, ignoreTrackID);
     previousID[k] = CaloHitID(timeSlice, ignoreTrackID);
   }
   isParameterized = false;
 
   entrancePoint.set(0., 0., 0.);
   entranceLocal.set(0., 0., 0.);
   posGlobal.set(0., 0., 0.);
   incidentEnergy = edepositEM = edepositHAD = 0.f;
 
   primAncestor = 0;
   forceSave = false;
   for (int k = 0; k < 2; ++k) {
     cleanIndex[k] = primIDSaved[k] = totalHits[k] = 0;
   }
 
   edm::LogVerbatim("CaloSim") << "CaloSD: Minimum energy of track for saving it " << energyCut / CLHEP::GeV
                               << " GeV\n        Use of HitID Map " << useMap << "\n        Check last "
                               << nCheckedHits[0] << " before saving the hit\n        Correct TOF globally by "
                               << correctT << " ns (Flag =" << corrTOFBeam << ")\n        Save hits recorded before "
                               << tmaxHit << " ns and if energy is above " << eminHit / CLHEP::MeV
                               << " MeV (for depth 0) or " << eminHitD / CLHEP::MeV
                               << " MeV (for nonzero depths);\n        Time Slice Unit " << timeSlice
                               << "\nIgnore TrackID Flag " << ignoreTrackID << " doFineCalo flag " << doFineCalo_
                               << "\nBeam Position " << beamZ / CLHEP::cm << " cm";
   if (doFineCalo_)
     edm::LogVerbatim("DoFineCalo") << "Using finecalo v2";
 
   // Treat fine calorimeters
   edm::LogVerbatim("CaloSim") << "CaloSD: Have a possibility of " << fineNames.size() << " fine calorimeters of which "
                               << useFines.size() << " are selected";
   for (unsigned int k = 0; k < fineNames.size(); ++k)
     edm::LogVerbatim("CaloSim") << "[" << k << "] " << fineNames[k] << " at " << fineLevels[k];
   std::ostringstream st1;
   for (unsigned int k = 0; k < useFines.size(); ++k)
     st1 << " [" << k << "] " << useFines[k] << ":" << fineNames[useFines[k]];
   edm::LogVerbatim("CaloSim") << "CaloSD used calorimeters" << st1.str();
   const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
   std::vector<G4LogicalVolume*>::const_iterator lvcite;
   for (unsigned int i = 0; i < useFines.size(); i++) {
     G4LogicalVolume* lv = nullptr;
     G4String name = static_cast<G4String>(fineNames[useFines[i]]);
     for (lvcite = lvs->begin(); lvcite != lvs->end(); lvcite++) {
       G4String namx(static_cast<std::string>(dd4hep::dd::noNamespace((*lvcite)->GetName())));
       if (namx == name) {
         lv = (*lvcite);
         break;
       }
     }
     if (lv != nullptr) {
       CaloSD::Detector detector;
       detector.name = name;
       detector.lv = lv;
       detector.level = fineLevels[useFines[i]];
       fineDetectors_.emplace_back(detector);
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD::Loads information for " << fineDetectors_.size() << " fine detectors";
   unsigned int k(0);
   for (const auto& detector : fineDetectors_) {
     edm::LogVerbatim("CaloSim") << "Detector[" << k << "] " << detector.name << " at level " << detector.level
                                 << " pointer to LV: " << detector.lv;
   }
 #endif
 }

◆ ~CaloSD()

CaloSD::~CaloSD ( )

override

Definition at line 174 of file CaloSD.cc.

174 {}

Member Function Documentation

◆ checkHit()

bool CaloSD::checkHit ( int k = 0 )

protected

Definition at line 488 of file CaloSD.cc.

References currentHit, currentID, mps_fire::end, newFWLiteAna::found, hitMap, dqmiolumiharvest::j, dqmdumpme::k, SiStripPI::max, nCheckedHits, nhits, theHC, updateHit(), and useMap.

Referenced by hitExists(), processHit(), and ProcessHits().

                            {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: checkHit for " << k;
 #endif
   //look in the HitContainer whether a hit with the same ID already exists:
   bool found = false;
   if (useMap) {
     std::map<CaloHitID, CaloG4Hit*>::const_iterator it = hitMap[k].find(currentID[k]);
     if (it != hitMap[k].end()) {
       currentHit[k] = it->second;
       found = true;
     }
   } else if (nCheckedHits[k] > 0) {
     int nhits = theHC[k]->entries();
     int minhit = std::max(nhits - nCheckedHits[k], 0);
     int maxhit = nhits - 1;
 
     for (int j = maxhit; j > minhit; --j) {
       if ((*theHC[k])[j]->getID() == currentID[k]) {
         currentHit[k] = (*theHC[k])[j];
         found = true;
         break;
       }
     }
   }
 
   if (found) {
     updateHit(currentHit[k], k);
   }
   return found;
 }

◆ cleanHitCollection()

void CaloSD::cleanHitCollection ( )

private

Definition at line 1039 of file CaloSD.cc.

References cleanIndex, correctT, corrTOFBeam, relativeConstraints::empty, cond::serialization::equal(), filterHit(), CaloG4Hit::getTimeSlice(), hitMap, mps_fire::i, cuy::ii, dqmiolumiharvest::j, dqmdumpme::k, nHC_, hltrates_dqm_sourceclient-live_cfg::offset, AlCaHLTBitMon_ParallelJobs::p, reusehit, jetUpdater_cfi::sort, theHC, hcalRecHitTable_cff::time, totalHits, and useMap.

Referenced by EndOfEvent(), and update().

                                 {
   for (int k = 0; k < nHC_; ++k) {
     std::vector<CaloG4Hit*>* theCollection = theHC[k]->GetVector();
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD: collection before merging, size = " << theHC[k]->entries();
 #endif
     if (reusehit[k].empty())
       reusehit[k].reserve(theHC[k]->entries() - cleanIndex[k]);
 
     // if no map used, merge before hits to have the save situation as a map
     if (!useMap) {
       std::vector<CaloG4Hit*> hitvec;
 
       hitvec.swap(*theCollection);
       sort((hitvec.begin() + cleanIndex[k]), hitvec.end(), CaloG4HitLess());
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("CaloSim") << "CaloSD::cleanHitCollection: sort hits in buffer starting from "
                                   << "element = " << cleanIndex[k];
       for (unsigned int i = 0; i < hitvec.size(); ++i) {
         if (hitvec[i] == nullptr)
           edm::LogVerbatim("CaloSim") << i << " has a null pointer";
         else
           edm::LogVerbatim("CaloSim") << i << " " << *hitvec[i];
       }
 #endif
       CaloG4HitEqual equal;
       for (unsigned int i = cleanIndex[k]; i < hitvec.size(); ++i) {
         int jump = 0;
         for (unsigned int j = i + 1; j < hitvec.size() && equal(hitvec[i], hitvec[j]); ++j) {
           ++jump;
           // merge j to i
           (*hitvec[i]).addEnergyDeposit(*hitvec[j]);
           (*hitvec[j]).setEM(0.);
           (*hitvec[j]).setHadr(0.);
           reusehit[k].emplace_back(hitvec[j]);
           hitvec[j] = nullptr;
         }
         i += jump;
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("CaloSim") << "CaloSD: cleanHitCollection merge the hits in buffer ";
       for (unsigned int i = 0; i < hitvec.size(); ++i) {
         if (hitvec[i] == nullptr)
           edm::LogVerbatim("CaloSim") << i << " has a null pointer";
         else
           edm::LogVerbatim("CaloSim") << i << " " << *hitvec[i];
       }
 #endif
       //move all nullptr to end of list and then remove them
       hitvec.erase(std::stable_partition(
                        hitvec.begin() + cleanIndex[k], hitvec.end(), [](CaloG4Hit* p) { return p != nullptr; }),
                    hitvec.end());
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("CaloSim") << "CaloSD::cleanHitCollection: remove the merged hits in buffer, new size = "
                                   << hitvec.size();
 #endif
       hitvec.swap(*theCollection);
       totalHits[k] = theHC[k]->entries();
     }
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD: collection after merging, size= " << theHC[k]->entries()
                                 << " Size of reusehit= " << reusehit[k].size()
                                 << "\n      starting hit selection from index = " << cleanIndex[k];
 #endif
 
     int addhit = 0;
     for (unsigned int i = cleanIndex[k]; i < theCollection->size(); ++i) {
       CaloG4Hit* aHit((*theCollection)[i]);
 
       // selection
 
       double time = aHit->getTimeSlice();
       if (corrTOFBeam)
         time += correctT;
       if (!filterHit(aHit, time)) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("CaloSim") << "CaloSD: dropped CaloG4Hit "
                                     << " " << *aHit;
 #endif
 
         // create the list of hits to be reused
 
         reusehit[k].emplace_back((*theCollection)[i]);
         (*theCollection)[i] = nullptr;
         ++addhit;
       }
     }
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD: Size of reusehit after selection = " << reusehit[k].size()
                                 << " Number of added hit = " << addhit;
 #endif
     if (useMap) {
       if (addhit > 0) {
         int offset = reusehit[k].size() - addhit;
         for (int ii = addhit - 1; ii >= 0; --ii) {
           CaloHitID theID = reusehit[k][offset + ii]->getID();
           hitMap[k].erase(theID);
         }
       }
     }
 
     //move all nullptr to end of list and then remove them
     theCollection->erase(
         std::stable_partition(
             theCollection->begin() + cleanIndex[k], theCollection->end(), [](CaloG4Hit* p) { return p != nullptr; }),
         theCollection->end());
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD: hit collection after selection, size = " << theHC[k]->entries();
     theHC[k]->PrintAllHits();
 #endif
 
     cleanIndex[k] = theHC[k]->entries();
   }
 }

◆ clear()

void CaloSD::clear ( void )

override

Definition at line 428 of file CaloSD.cc.

428 {}

◆ clearHits()

void CaloSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 847 of file CaloSD.cc.

References cleanIndex, dqmdumpme::k, nHC_, previousID, primIDSaved, CaloHitID::reset(), and slave.

Referenced by update().

                        {
   for (int k = 0; k < nHC_; ++k) {
     previousID[k].reset();
     primIDSaved[k] = -99;
     cleanIndex[k] = 0;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD: Clears hit vector for " << GetName()
                                 << " and initialise slave: " << slave[k].get()->name();
 #endif
     slave[k].get()->Initialize();
   }
 }

◆ createNewHit()

CaloG4Hit * CaloSD::createNewHit	(	const G4Step *	aStep,
		const G4Track *	theTrack,
		int	k
	)

protected

Definition at line 623 of file CaloSD.cc.

References SensitiveDetector::cmsTrackInformation(), currentID, CaloHitID::depth(), doFineCaloThisStep_, edepositEM, edepositHAD, relativeConstraints::empty, energyCut, entranceLocal, entrancePoint, f, forceSave, getNumberOfHits(), incidentEnergy, isItFineCalo(), dqmdumpme::k, TrackWithHistory::momentum(), posGlobal, primIDSaved, reusehit, CaloG4Hit::setEM(), CaloG4Hit::setEntry(), CaloG4Hit::setEntryLocal(), CaloG4Hit::setHadr(), CaloG4Hit::setID(), CaloG4Hit::setIncidentEnergy(), CaloG4Hit::setPosition(), TrackInformation::setStoreTrack(), TrackWithHistory::setToBeSaved(), shortreprID(), mathSSE::sqrt(), storeHit(), CaloHitID::timeSliceID(), tkMap, totalHits, CaloHitID::trackID(), CaloHitID::unitID(), updateHit(), and useMap.

Referenced by processHit(), ZdcSD::ProcessHits(), ProcessHits(), and HGCalSD::processSecondHit().

                                                                                    {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD::CreateNewHit " << getNumberOfHits() << " for " << GetName()
                               << " Unit:" << currentID[k].unitID() << " " << currentID[k].depth()
                               << " Edep= " << edepositEM << " " << edepositHAD
                               << " primaryID= " << currentID[k].trackID()
                               << " timeSlice= " << currentID[k].timeSliceID() << " ID= " << theTrack->GetTrackID()
                               << " " << theTrack->GetDefinition()->GetParticleName()
                               << " E(GeV)= " << theTrack->GetKineticEnergy() / CLHEP::GeV
                               << " parentID= " << theTrack->GetParentID() << "\n Ein= " << incidentEnergy
                               << " entranceGlobal: " << entrancePoint << " entranceLocal: " << entranceLocal
                               << " posGlobal: " << posGlobal;
 #endif
 
   CaloG4Hit* aHit;
   if (!reusehit[k].empty()) {
     aHit = reusehit[k].back().release();
     aHit->setEM(0.f);
     aHit->setHadr(0.f);
     reusehit[k].pop_back();
   } else {
     aHit = new CaloG4Hit;
   }
 
   aHit->setID(currentID[k]);
   aHit->setEntry(entrancePoint.x(), entrancePoint.y(), entrancePoint.z());
   aHit->setEntryLocal(entranceLocal.x(), entranceLocal.y(), entranceLocal.z());
   aHit->setPosition(posGlobal.x(), posGlobal.y(), posGlobal.z());
   aHit->setIncidentEnergy(incidentEnergy);
   updateHit(aHit, k);
 
   storeHit(aHit, k);
   TrackInformation* trkInfo = cmsTrackInformation(theTrack);
 
 #ifdef EDM_ML_DEBUG
   if (doFineCaloThisStep_)
     edm::LogVerbatim("DoFineCalo") << "New hit " << shortreprID(aHit) << " using finecalo;"
                                    << " isItFineCalo(post)=" << isItFineCalo(aStep->GetPostStepPoint()->GetTouchable())
                                    << " isItFineCalo(pre)=" << isItFineCalo(aStep->GetPreStepPoint()->GetTouchable());
 #endif
 
   // 'Traditional', non-fine history bookkeeping
   if (!doFineCaloThisStep_) {
     double etrack = 0;
     if (currentID[k].trackID() == primIDSaved[k]) {  // The track is saved; nothing to be done
     } else if (currentID[k].trackID() == theTrack->GetTrackID()) {
       etrack = theTrack->GetKineticEnergy();
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("CaloSim") << "CaloSD: set save the track " << currentID[k].trackID() << " etrack " << etrack
                                   << " eCut " << energyCut << " force: " << forceSave
                                   << " save: " << (etrack >= energyCut || forceSave);
 #endif
       if (etrack >= energyCut || forceSave) {
         trkInfo->setStoreTrack();
       }
     } else {
       TrackWithHistory* trkh = tkMap[currentID[k].trackID()];
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("CaloSim") << "CaloSD : TrackWithHistory pointer for " << currentID[k].trackID() << " is "
                                   << trkh;
 #endif
       if (trkh != nullptr) {
         etrack = sqrt(trkh->momentum().Mag2());
         if (etrack >= energyCut) {
           trkh->setToBeSaved();
 #ifdef EDM_ML_DEBUG
           edm::LogVerbatim("CaloSim") << "CaloSD: set save the track " << currentID[k].trackID() << " with Hit";
 #endif
         }
       }
     }
     primIDSaved[k] = currentID[k].trackID();
   }
 
   if (useMap)
     ++totalHits[k];
   return aHit;
 }

◆ DrawAll()

void CaloSD::DrawAll ( )

override

Definition at line 430 of file CaloSD.cc.

430 {}

◆ endEvent()

void CaloSD::endEvent ( )

protectedvirtual

Reimplemented in HCalSD, and HGCSD.

Definition at line 870 of file CaloSD.cc.

Referenced by update().

870 {}

◆ EndOfEvent()

void CaloSD::EndOfEvent ( G4HCofThisEvent * eventHC )

override

Definition at line 412 of file CaloSD.cc.

References cleanHitCollection(), dqmdumpme::k, nHC_, and theHC.

                                         {
   // clean the hits for the last tracks
 
   cleanHitCollection();
 
 #ifdef EDM_ML_DEBUG
   for (int k = 0; k < nHC_; ++k) {
     if (theHC[k] == nullptr)
       edm::LogVerbatim("CaloSim") << "CaloSD: EndofEvent entered for container " << k << " with no entries";
     else
       edm::LogVerbatim("CaloSim") << "CaloSD: EndofEvent entered for container " << k << " with " << theHC[0]->entries()
                                   << " entries";
   }
 #endif
 }

◆ EnergyCorrected()

double CaloSD::EnergyCorrected	(	const G4Step &	step,
		const G4Track *
	)

protectedvirtual

Reimplemented in ECalSD.

Definition at line 364 of file CaloSD.cc.

Referenced by ProcessHits().

364 { return aStep.GetTotalEnergyDeposit(); }

◆ fillHits()

void CaloSD::fillHits	(	edm::PCaloHitContainer &	cc,
		const std::string &	hname
	)

overridevirtual

Implements SensitiveCaloDetector.

Definition at line 441 of file CaloSD.cc.

References gpuPixelDoublets::cc, dqmdumpme::k, Skims_PA_cff::name, nHC_, and slave.

                                                                       {
   for (int k = 0; k < nHC_; ++k) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD: Tries to transfer " << slave[k].get()->hits().size() << " hits for "
                                 << slave[k].get()->name() << "   " << hname;
 #endif
     if (slave[k].get()->name() == hname) {
       cc = slave[k].get()->hits();
       slave[k].get()->Clean();
     }
   }
 }

◆ filterHit()

bool CaloSD::filterHit	(	CaloG4Hit *	hit,
		double	time
	)

protectedvirtual

Reimplemented in HCalSD, HGCSD, HGCalSD, HFNoseSD, HGCScintSD, and AHCalSD.

Definition at line 916 of file CaloSD.cc.

References eminHit, eminHitD, hcalRecHitTable_cff::time, and tmaxHit.

Referenced by cleanHitCollection().

                                                   {
   double emin(eminHit);
   if (hit->getDepth() > 0)
     emin = eminHitD;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD::filterHit(..) Depth " << hit->getDepth() << " Emin = " << emin << " ("
                               << eminHit << ", " << eminHitD << ")";
 #endif
   return ((time <= tmaxHit) && (hit->getEnergyDeposit() > emin));
 }

◆ findBoundaryCrossingParent()

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

protected

Definition at line 551 of file CaloSD.cc.

References boundaryCrossingParentMap_, SensitiveDetector::cmsTrackInformation(), TrackInformation::crossedBoundary(), Exception, SimTrackManager::getTrackByID(), l1ctLayer2EG_cff::id, m_trackManager, TrackWithHistory::parentID(), cosmicMuonLinks_cfi::parentTrack, printableDecayChain(), TrackInformation::setStoreTrack(), and HLT_2023v12_cff::track.

Referenced by ProcessHits().

                                                                                          {
   TrackInformation* trkInfo = cmsTrackInformation(track);
   unsigned int id = track->GetTrackID();
   // First see if this track is already in the map
   auto it = boundaryCrossingParentMap_.find(id);
   if (it != boundaryCrossingParentMap_.end()) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("DoFineCalo") << "Track " << id << " parent already cached: " << it->second;
 #endif
     return it->second;
   }
   // Then see if the track itself crosses the boundary
   else if (trkInfo->crossedBoundary()) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("DoFineCalo") << "Track " << id << " crosses boundary itself";
 #endif
     boundaryCrossingParentMap_[id] = id;
     trkInfo->setStoreTrack();
     return id;
   }
   // Else, traverse the history of the track
   std::vector<unsigned int> decayChain{id};
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("DoFineCalo") << "Track " << id << ": Traversing history to find boundary-crossing parent";
 #endif
   unsigned int parentID = track->GetParentID();
   while (true) {
     if (parentID == 0)
       throw cms::Exception("Unknown", "CaloSD")
           << "Hit end of parentage for track " << id << " without finding a boundary-crossing parent";
     // First check if this ancestor is already in the map
     auto it = boundaryCrossingParentMap_.find(parentID);
     if (it != boundaryCrossingParentMap_.end()) {
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("DoFineCalo") << "  Track " << parentID
                                      << " boundary-crossing parent already cached: " << it->second;
 #endif
       // Store this parent also for the rest of the traversed decay chain
       for (auto ancestorID : decayChain)
         boundaryCrossingParentMap_[ancestorID] = it->second;
 #ifdef EDM_ML_DEBUG
       // In debug mode, still build the rest of the decay chain for debugging
       decayChain.push_back(parentID);
       while (parentID != it->second) {
         parentID = m_trackManager->getTrackByID(parentID, true)->parentID();
         decayChain.push_back(parentID);
       }
       edm::LogVerbatim("DoFineCalo") << "  Full decay chain: " << printableDecayChain(decayChain);
 #endif
       return it->second;
     }
     // If not, get this parent from the track manager (expensive)
     TrackWithHistory* parentTrack = m_trackManager->getTrackByID(parentID, true);
     if (parentTrack->crossedBoundary()) {
       if (markAsSaveable)
         parentTrack->setToBeSaved();
       decayChain.push_back(parentID);
       // Record this boundary crossing parent for all traversed ancestors
       for (auto ancestorID : decayChain)
         boundaryCrossingParentMap_[ancestorID] = parentID;
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("DoFineCalo") << "  Found boundary-crossing ancestor " << parentID << " for track " << id
                                      << "; decay chain: " << printableDecayChain(decayChain);
 #endif
       return parentID;
     }
     // Next iteration
     decayChain.push_back(parentID);
     parentID = parentTrack->parentID();
   }
 }

◆ getAttenuation()

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

protected

Definition at line 726 of file CaloSD.cc.

References funct::abs(), HltBtagPostValidation_cff::c, ALCARECOTkAlJpsiMuMu_cff::charge, fastSimProducer_cff::density, and mps_merge::weight.

Referenced by TotemT2ScintSD::getEnergyDeposit(), AHCalSD::getEnergyDeposit(), HGCalTB16SD01::getEnergyDeposit(), HcalTB06BeamSD::getEnergyDeposit(), HGCScintSD::getEnergyDeposit(), EcalTBH4BeamSD::getEnergyDeposit(), DreamSD::getEnergyDeposit(), ECalSD::getEnergyDeposit(), HcalTB02SD::getEnergyDeposit(), and HCalSD::getEnergyDeposit().

                                                                                                  {
   double weight = 1.;
   double charge = aStep->GetPreStepPoint()->GetCharge();
   double length = aStep->GetStepLength();
 
   if (charge != 0. && length > 0.) {
     double density = aStep->GetPreStepPoint()->GetMaterial()->GetDensity();
     double dedx = aStep->GetTotalEnergyDeposit() / length;
     double rkb = birk1 / density;
     double c = birk2 * rkb * rkb;
     if (std::abs(charge) >= 2.)
       rkb /= birk3;  // based on alpha particle data
     weight = 1. / (1. + rkb * dedx + c * dedx * dedx);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD::getAttenuation in " << aStep->GetPreStepPoint()->GetMaterial()->GetName()
                                 << " Charge " << charge << " dE/dx " << dedx << " Birk Const " << rkb << ", " << c
                                 << " Weight = " << weight << " dE " << aStep->GetTotalEnergyDeposit();
 #endif
   }
   return weight;
 }

◆ getDepth()

uint16_t CaloSD::getDepth ( const G4Step * )

protectedvirtual

Reimplemented in ECalSD.

Definition at line 912 of file CaloSD.cc.

Referenced by ZdcSD::ProcessHits(), and ProcessHits().

912 { return 0; }

◆ getEnergyDeposit()

double CaloSD::getEnergyDeposit ( const G4Step * step )

protectedvirtual

Reimplemented in HCalSD, ECalSD, HcalTB02SD, CastorSD, HGCSD, DreamSD, HGCalSD, HFNoseSD, EcalTBH4BeamSD, HGCScintSD, HcalTB06BeamSD, HGCalTB16SD01, AHCalSD, and TotemT2ScintSD.

Definition at line 362 of file CaloSD.cc.

Referenced by ProcessHits().

362 { return aStep->GetTotalEnergyDeposit(); }

◆ getFromLibrary()

bool CaloSD::getFromLibrary ( const G4Step * step )

protectedvirtual

Reimplemented in HCalSD, and CastorSD.

Definition at line 366 of file CaloSD.cc.

Referenced by ProcessHits().

366 { return false; }

◆ getNumberOfHits()

int CaloSD::getNumberOfHits ( int k = 0 )

protected

Definition at line 520 of file CaloSD.cc.

References dqmdumpme::k, and theHC.

Referenced by createNewHit(), and HCalSD::getFromLibrary().

520 { return theHC[k]->entries(); }

CaloSD::theHC

CaloG4HitCollection * theHC[2]

Definition: CaloSD.h:174

dqmdumpme.k

k

Definition: dqmdumpme.py:60

◆ getResponseWt()

double CaloSD::getResponseWt	(	const G4Track *	aTrack,
		int	k = `0`
	)

protected

Definition at line 927 of file CaloSD.cc.

References SensitiveDetector::cmsTrackInformation(), TrackInformation::genParticleP(), TrackInformation::genParticlePID(), dqmdumpme::k, and meanResponse.

Referenced by ECalSD::EnergyCorrected(), HGCScintSD::getEnergyDeposit(), HFNoseSD::getEnergyDeposit(), HGCalSD::getEnergyDeposit(), HGCSD::getEnergyDeposit(), ECalSD::getEnergyDeposit(), and HCalSD::getEnergyDeposit().

                                                          {
   double wt = 1.0;
   if (meanResponse[k].get()) {
     TrackInformation* trkInfo = cmsTrackInformation(aTrack);
     wt = meanResponse[k].get()->getWeight(trkInfo->genParticlePID(), trkInfo->genParticleP());
   }
   return wt;
 }

◆ getTrackID()

int CaloSD::getTrackID ( const G4Track * aTrack )

protectedvirtual

Reimplemented in ECalSD.

Definition at line 872 of file CaloSD.cc.

References SensitiveDetector::cmsTrackInformation(), and TrackInformation::getIDonCaloSurface().

Referenced by ECalSD::getTrackID(), ZdcSD::ProcessHits(), and ProcessHits().

                                             {
   int primaryID = 0;
   TrackInformation* trkInfo = cmsTrackInformation(aTrack);
   if (trkInfo) {
     primaryID = trkInfo->getIDonCaloSurface();
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "Track ID: " << trkInfo->getIDonCaloSurface() << ":" << aTrack->GetTrackID() << ":"
                                 << primaryID;
 #endif
   } else {
     primaryID = aTrack->GetTrackID();
 #ifdef EDM_ML_DEBUG
     edm::LogWarning("CaloSim") << "CaloSD: Problem with primaryID **** set by force to TkID **** " << primaryID;
 #endif
   }
   return primaryID;
 }

◆ hitExists()

bool CaloSD::hitExists	(	const G4Step *	aStep,
		int	k
	)

protected

Definition at line 462 of file CaloSD.cc.

References checkHit(), currentHit, currentID, dqmdumpme::k, posGlobal, previousID, resetForNewPrimary(), and updateHit().

Referenced by ZdcSD::ProcessHits(), ProcessHits(), and HGCalSD::processSecondHit().

                                                  {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: hitExists for " << k;
 #endif
   // Update if in the same detector, time-slice and for same track
   if (currentID[k] == previousID[k]) {
     updateHit(currentHit[k], k);
     return true;
   }
 
   // Note T. Klijnsma:
   // This is a rather strange place to set these class variables.
   // The code would be much more readable if all logic for determining
   // whether to update a hit or create a new hit is done in one place,
   // and only then perform the actual updating or creating of the hit.
 
   // Reset entry point for new primary (only for the primary hit)
   if (k == 0) {
     posGlobal = aStep->GetPreStepPoint()->GetPosition();
     if (currentID[k].trackID() != previousID[k].trackID()) {
       resetForNewPrimary(aStep);
     }
   }
   return checkHit(k);
 }

◆ ignoreRejection()

void CaloSD::ignoreRejection ( )

inlineprotected

Definition at line 112 of file CaloSD.h.

References ignoreReject.

Referenced by HGCalSD::setDetUnitId().

112 { ignoreReject = true; }

CaloSD::ignoreReject

bool ignoreReject

Definition: CaloSD.h:178

◆ initEvent()

void CaloSD::initEvent ( const BeginOfEvent * )

protectedvirtual

Reimplemented in HCalSD, and HGCSD.

Definition at line 868 of file CaloSD.cc.

Referenced by update().

868 {}

◆ Initialize()

void CaloSD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 388 of file CaloSD.cc.

References collName_, detName_, hcID, dqmdumpme::k, nHC_, theHC, and totalHits.

                                             {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD : Initialize called for " << detName_;
 #endif
 
   for (int k = 0; k < nHC_; ++k) {
     totalHits[k] = 0;
 
     //This initialization is performed at the beginning of an event
     //------------------------------------------------------------
     theHC[k] = new CaloG4HitCollection(detName_, collName_[k]);
 
     if (hcID[k] < 0) {
       hcID[k] = G4SDManager::GetSDMpointer()->GetCollectionID(collName_[k]);
     }
     //theHC ownership is transfered here to HCE
     HCE->AddHitsCollection(hcID[k], theHC[k]);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD : Add hits collection for " << hcID[k] << " for " << k << ":" << detName_
                                 << ":" << collName_[k];
 #endif
   }
 }

◆ initRun()

void CaloSD::initRun ( )

protectedvirtual

Reimplemented in HCalSD, HGCSD, HGCalSD, HFNoseSD, DreamSD, ZdcSD, and HGCScintSD.

Definition at line 866 of file CaloSD.cc.

Referenced by update().

866 {}

◆ isItFineCalo()

bool CaloSD::isItFineCalo ( const G4VTouchable * touch )

Definition at line 368 of file CaloSD.cc.

References hgcalTestNeighbor_cfi::detector, fineDetectors_, cuy::ii, personalPlayback::level, convertSQLiteXML::ok, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by createNewHit(), and ProcessHits().

                                                    {
   bool ok(false);
   int level = ((touch->GetHistoryDepth()) + 1);
   for (const auto& detector : fineDetectors_) {
     if (level > 0 && level >= detector.level) {
       int ii = level - detector.level;
       G4LogicalVolume* lv = touch->GetVolume(ii)->GetLogicalVolume();
       ok = (lv == detector.lv);
 #ifdef EDM_ML_DEBUG
       std::string name1 = (lv == nullptr) ? "Unknown" : lv->GetName();
       edm::LogVerbatim("CaloSim") << "CaloSD: volume " << name1 << ":" << detector.name << " at Level "
                                   << detector.level << " Flag " << ok;
 #endif
       if (ok)
         break;
     }
   }
   return ok;
 }

◆ newCollection()

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

Definition at line 176 of file CaloSD.cc.

References collName_, detName_, hcn_, dqmdumpme::k, meanResponse, Skims_PA_cff::name, nHC_, AlCaHLTBitMon_ParallelJobs::p, slave, and useResMap_.

Referenced by HGCalSD::update().

                                                                           {
   for (unsigned int k = 0; k < hcn_.size(); ++k) {
     if (name == hcn_[k]) {
       if (k < useResMap_.size() && useResMap_[k] > 0) {
         meanResponse[1] = std::make_unique<CaloMeanResponse>(p);
         break;
       }
     }
   }
   slave[1] = std::make_unique<CaloSlaveSD>(name);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD:: Initialise a second collection for " << name;
   for (int k = 0; k < nHC_; ++k)
     edm::LogVerbatim("CaloSim") << "Detector[" << k << "]" << detName_ << " collection " << collName_[k];
 
 #endif
 }

◆ printableDecayChain()

std::string CaloSD::printableDecayChain ( const std::vector< unsigned int > & decayChain )

staticprotected

Definition at line 526 of file CaloSD.cc.

References mps_fire::i, and contentValuesCheck::ss.

Referenced by findBoundaryCrossingParent().

                                                                                {
   std::stringstream ss;
   for (long unsigned int i = 0; i < decayChain.size(); i++) {
     if (i > 0)
       ss << " <- ";
     ss << decayChain[i];
   }
   return ss.str();
 }

◆ PrintAll()

void CaloSD::PrintAll ( )

override

Definition at line 432 of file CaloSD.cc.

References dqmdumpme::k, nHC_, and theHC.

                       {
   for (int k = 0; k < nHC_; ++k) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD: Collection " << theHC[k]->GetName();
 #endif
     theHC[k]->PrintAllHits();
   }
 }

◆ printDetectorLevels()

void CaloSD::printDetectorLevels ( const G4VTouchable * touch ) const

protected

Definition at line 1157 of file CaloSD.cc.

References mps_fire::i, cuy::ii, personalPlayback::level, Skims_PA_cff::name, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by HGCScintSD::setDetUnitId(), and HGCalSD::setDetUnitId().

                                                                 {
   //Print name and copy numbers
   int level = ((touch->GetHistoryDepth()) + 1);
   std::ostringstream st1;
   st1 << level << " Levels:";
   if (level > 0) {
     for (int ii = 0; ii < level; ii++) {
       int i = level - ii - 1;
       G4VPhysicalVolume* pv = touch->GetVolume(i);
       std::string name = (pv != nullptr) ? pv->GetName() : "Unknown";
       st1 << " " << name << ":" << touch->GetReplicaNumber(i);
     }
   }
   edm::LogVerbatim("CaloSim") << st1.str();
 }

◆ processHit()

void CaloSD::processHit ( const G4Step * step )

inlineprotected

Definition at line 117 of file CaloSD.h.

References checkHit(), createNewHit(), currentHit, currentID, previousID, and updateHit().

Referenced by HCalSD::getFromHFLibrary(), CastorSD::getFromLibrary(), HCalSD::getFromParam(), HCalSD::getHitFibreBundle(), HCalSD::getHitPMT(), and HCalSD::hitForFibre().

                                              {
     // check if it is in the same unit and timeslice as the previous one
     if (currentID[0] == previousID[0]) {
       updateHit(currentHit[0], 0);
     } else if (!checkHit()) {
       currentHit[0] = createNewHit(step, step->GetTrack(), 0);
     }
   }

◆ ProcessHits() [1/2]

G4bool CaloSD::ProcessHits	(	G4Step *	step,
		G4TouchableHistory *
	)

overridevirtual

Implements SensitiveDetector.

Reimplemented in ZdcSD.

Definition at line 194 of file CaloSD.cc.

References createNewHit(), currentHit, currentID, TauDecayModes::dec, hcalRecHitTable_cff::depth, doFineCalo_, doFineCaloThisStep_, edepositEM, edepositHAD, hcalRecHitTable_cff::energy, findBoundaryCrossingParent(), getDepth(), getEnergyDeposit(), getFromLibrary(), getTrackID(), hitExists(), ignoreReject, G4TrackToParticleID::isGammaElectronPositron(), isItFineCalo(), isParameterized, CaloHitID::markAsFinecaloTrackID(), SensitiveDetector::NaNTrap(), previousID, processSecondHit(), setDetUnitId(), CaloHitID::setID(), shortreprID(), hcalRecHitTable_cff::time, and CaloHitID::unitID().

                                                              {
   NaNTrap(aStep);
   ignoreReject = false;
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD::" << GetName() << " ID= " << aStep->GetTrack()->GetTrackID()
                               << " prID= " << aStep->GetTrack()->GetParentID()
                               << " Eprestep= " << aStep->GetPreStepPoint()->GetKineticEnergy()
                               << " step= " << aStep->GetStepLength() << " Edep= " << aStep->GetTotalEnergyDeposit();
 #endif
 
   // Class variable to determine whether finecalo rules should apply for this step
   doFineCaloThisStep_ = (doFineCalo_ && isItFineCalo(aStep->GetPreStepPoint()->GetTouchable()));
 
   // apply shower library or parameterisation
   // independent on energy deposition at a step
   if (isParameterized) {
     if (getFromLibrary(aStep)) {
       // for parameterized showers the primary track should be killed
       // secondary tracks should be killed if they are in the same volume
       (aStep->GetTrack())->SetTrackStatus(fStopAndKill);
       if (0 < aStep->GetNumberOfSecondariesInCurrentStep()) {
         auto tv = aStep->GetSecondaryInCurrentStep();
         auto vol = aStep->GetPreStepPoint()->GetPhysicalVolume();
         for (auto& tk : *tv) {
           if (tk->GetVolume() == vol) {
             const_cast<G4Track*>(tk)->SetTrackStatus(fStopAndKill);
           }
         }
       }
       return true;
     }
   }
 
   // ignore steps without energy deposit
   edepositEM = edepositHAD = 0.f;
   if (aStep->GetTotalEnergyDeposit() <= 0.0) {
     return false;
   }
 
   // check unitID
   unsigned int unitID = setDetUnitId(aStep);
   auto const theTrack = aStep->GetTrack();
   uint16_t depth = getDepth(aStep);
 
   double time = theTrack->GetGlobalTime() / nanosecond;
   int primaryID = getTrackID(theTrack);
   if (unitID > 0) {
     currentID[0].setID(unitID, time, primaryID, depth);
   } else {
     if (!ignoreReject) {
       const G4TouchableHistory* touch = static_cast<const G4TouchableHistory*>(theTrack->GetTouchable());
       edm::LogVerbatim("CaloSim") << "CaloSD::ProcessHits: unitID= " << unitID
                                   << " currUnit=   " << currentID[0].unitID() << " Detector: " << GetName()
                                   << " trackID= " << theTrack->GetTrackID() << " "
                                   << theTrack->GetDefinition()->GetParticleName()
                                   << "\n Edep= " << aStep->GetTotalEnergyDeposit()
                                   << " PV: " << touch->GetVolume(0)->GetName()
                                   << " PVid= " << touch->GetReplicaNumber(0) << " MVid= " << touch->GetReplicaNumber(1);
     }
     return false;
   }
   double energy = getEnergyDeposit(aStep);
   if (energy <= 0.0) {
     return false;
   }
 
   if (doFineCaloThisStep_) {
     currentID[0].setID(unitID, time, findBoundaryCrossingParent(theTrack), depth);
     currentID[0].markAsFinecaloTrackID();
   }
 
   if (G4TrackToParticleID::isGammaElectronPositron(theTrack)) {
     edepositEM = energy;
   } else {
     edepositHAD = energy;
   }
 #ifdef EDM_ML_DEBUG
   const G4TouchableHistory* touch = static_cast<const G4TouchableHistory*>(theTrack->GetTouchable());
   edm::LogVerbatim("CaloSim") << "CaloSD::" << GetName() << " PV:" << touch->GetVolume(0)->GetName()
                               << " PVid=" << touch->GetReplicaNumber(0) << " MVid=" << touch->GetReplicaNumber(1)
                               << " Unit:" << std::hex << unitID << std::dec << " Edep=" << edepositEM << " "
                               << edepositHAD << " ID=" << theTrack->GetTrackID() << " pID=" << theTrack->GetParentID()
                               << " E=" << theTrack->GetKineticEnergy() << " S=" << aStep->GetStepLength() << "\n "
                               << theTrack->GetDefinition()->GetParticleName() << " primaryID= " << primaryID
                               << " currentID= (" << currentID[0] << ") previousID= (" << previousID[0] << ")";
 #endif
   if (!hitExists(aStep, 0)) {
     currentHit[0] = createNewHit(aStep, theTrack, 0);
   } else {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("DoFineCalo") << "Not creating new hit, only updating " << shortreprID(currentHit[0]);
 #endif
   }
 
   processSecondHit(aStep, theTrack);
 
   return true;
 }

◆ ProcessHits() [2/2]

bool CaloSD::ProcessHits	(	G4GFlashSpot *	aSpot,
		G4TouchableHistory *
	)

override

Definition at line 294 of file CaloSD.cc.

References checkHit(), createNewHit(), currentHit, currentID, TauDecayModes::dec, hcalRecHitTable_cff::depth, edepositEM, edepositHAD, EnergyCorrected(), entranceLocal, entrancePoint, getDepth(), getTrackID(), incidentEnergy, G4TrackToParticleID::isGammaElectronPositron(), posGlobal, previousID, setDetUnitId(), CaloHitID::setID(), hcalRecHitTable_cff::time, HLT_2023v12_cff::track, CaloHitID::unitID(), and updateHit().

                                                                  {
   edepositEM = edepositHAD = 0.f;
   const G4Track* track = aSpot->GetOriginatorTrack()->GetPrimaryTrack();
 
   double edep = aSpot->GetEnergySpot()->GetEnergy();
   if (edep <= 0.0) {
     return false;
   }
 
   G4Step fFakeStep;
   G4StepPoint* fFakePreStepPoint = fFakeStep.GetPreStepPoint();
   G4StepPoint* fFakePostStepPoint = fFakeStep.GetPostStepPoint();
   fFakePreStepPoint->SetPosition(aSpot->GetPosition());
   fFakePostStepPoint->SetPosition(aSpot->GetPosition());
 
   G4TouchableHandle fTouchableHandle = aSpot->GetTouchableHandle();
   fFakePreStepPoint->SetTouchableHandle(fTouchableHandle);
   fFakeStep.SetTotalEnergyDeposit(edep);
   edep = EnergyCorrected(fFakeStep, track);
 
   // zero edep means hit outside the calorimeter
   if (edep <= 0.0) {
     return false;
   }
 
   if (G4TrackToParticleID::isGammaElectronPositron(track)) {
     edepositEM = edep;
   } else {
     edepositHAD = edep;
   }
 
   unsigned int unitID = setDetUnitId(&fFakeStep);
 
   if (unitID > 0) {
     // time of initial track
     double time = track->GetGlobalTime() / nanosecond;
     int primaryID = getTrackID(track);
     uint16_t depth = getDepth(&fFakeStep);
     currentID[0].setID(unitID, time, primaryID, depth);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD:: GetSpotInfo for Unit 0x" << std::hex << currentID[0].unitID() << std::dec
                                 << " Edeposit = " << edepositEM << " " << edepositHAD;
 #endif
     // Update if in the same detector, time-slice and for same track
     if (currentID[0] == previousID[0]) {
       updateHit(currentHit[0], 0);
     } else {
       posGlobal = aSpot->GetEnergySpot()->GetPosition();
       // Reset entry point for new primary
       if (currentID[0].trackID() != previousID[0].trackID()) {
         entrancePoint = aSpot->GetPosition();
         entranceLocal = aSpot->GetTouchableHandle()->GetHistory()->GetTopTransform().TransformPoint(entrancePoint);
         incidentEnergy = track->GetKineticEnergy();
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("CaloSim") << "CaloSD: Incident energy " << incidentEnergy / CLHEP::GeV << " GeV and"
                                     << " entrance point " << entrancePoint << " (Global) " << entranceLocal
                                     << " (Local)";
 #endif
       }
       if (!checkHit()) {
         currentHit[0] = createNewHit(&fFakeStep, track, 0);
       }
     }
     return true;
   }
   return false;
 }

◆ processSecondHit()

void CaloSD::processSecondHit	(	const G4Step *	,
		const G4Track *
	)

protectedvirtual

Reimplemented in HGCalSD.

Definition at line 914 of file CaloSD.cc.

Referenced by ProcessHits().

914 {}

◆ reset()

void CaloSD::reset ( void )

overridevirtual

Reimplemented from SensitiveCaloDetector.

Definition at line 860 of file CaloSD.cc.

References fpCaloG4HitAllocator.

Referenced by CaloSD().

                    {
   if (fpCaloG4HitAllocator) {
     fpCaloG4HitAllocator->ResetStorage();
   }
 }

◆ resetForNewPrimary()

void CaloSD::resetForNewPrimary ( const G4Step * aStep )

protected

Definition at line 713 of file CaloSD.cc.

References entranceLocal, entrancePoint, incidentEnergy, and setToLocal().

Referenced by HCalSD::getFromHFLibrary(), CastorSD::getFromLibrary(), HCalSD::getHitFibreBundle(), HCalSD::getHitPMT(), hitExists(), ZdcSD::setTrackID(), and setTrackID().

                                                    {
   auto const preStepPoint = aStep->GetPreStepPoint();
   entrancePoint = preStepPoint->GetPosition();
   entranceLocal = setToLocal(entrancePoint, preStepPoint->GetTouchable());
   incidentEnergy = preStepPoint->GetKineticEnergy();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD::resetForNewPrimary for " << GetName()
                               << " ID= " << aStep->GetTrack()->GetTrackID() << " Ein= " << incidentEnergy / CLHEP::GeV
                               << " GeV and"
                               << " entrance point global: " << entrancePoint << " local: " << entranceLocal;
 #endif
 }

◆ saveHit()

bool CaloSD::saveHit	(	CaloG4Hit *	aHit,
		int	k = `0`
	)

private

Definition at line 948 of file CaloSD.cc.

References correctT, corrTOFBeam, TauDecayModes::dec, doFineCalo_, Exception, CaloG4Hit::getDepth(), CaloG4Hit::getEM(), CaloG4Hit::getHadr(), CaloG4Hit::getTimeSlice(), CaloG4Hit::getTrackID(), CaloG4Hit::getUnitID(), SimTrackManager::giveMotherNeeded(), CaloG4Hit::isFinecaloTrackID(), dqmdumpme::k, m_trackManager, convertSQLiteXML::ok, shortreprID(), slave, hcalRecHitTable_cff::time, and SimTrackManager::trackExists().

Referenced by update().

                                            {
   int tkID;
   bool ok = true;
 
   double time = aHit->getTimeSlice();
   if (corrTOFBeam)
     time += correctT;
 
   // More strict bookkeeping for finecalo
   if (doFineCalo_ && aHit->isFinecaloTrackID()) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("DoFineCalo") << "Saving hit " << shortreprID(aHit);
 #endif
     if (!m_trackManager)
       throw cms::Exception("Unknown", "CaloSD") << "m_trackManager not set, needed for finecalo!";
     if (!m_trackManager->trackExists(aHit->getTrackID()))
       throw cms::Exception("Unknown", "CaloSD")
           << "Error on hit " << shortreprID(aHit) << ": Parent track not in track manager";
     slave[k].get()->processHits(aHit->getUnitID(),
                                 aHit->getEM() / CLHEP::GeV,
                                 aHit->getHadr() / CLHEP::GeV,
                                 time,
                                 aHit->getTrackID(),
                                 aHit->getDepth());
   }
   // Regular, not-fine way:
   else {
     if (m_trackManager) {
       tkID = m_trackManager->giveMotherNeeded(aHit->getTrackID());
       if (tkID == 0) {
         if (m_trackManager->trackExists(aHit->getTrackID()))
           tkID = (aHit->getTrackID());
         else {
           ok = false;
         }
       }
     } else {
       tkID = aHit->getTrackID();
       ok = false;
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("DoFineCalo") << "Saving hit " << shortreprID(aHit) << " with trackID=" << tkID;
 #endif
 
     slave[k].get()->processHits(
         aHit->getUnitID(), aHit->getEM() / CLHEP::GeV, aHit->getHadr() / CLHEP::GeV, time, tkID, aHit->getDepth());
   }
 
 #ifdef EDM_ML_DEBUG
   if (!ok)
     edm::LogWarning("CaloSim") << "CaloSD:Cannot find track ID for " << aHit->getTrackID();
   edm::LogVerbatim("CaloSim") << "CalosD: Track ID " << aHit->getTrackID() << " changed to " << tkID
                               << " by SimTrackManager Status " << ok;
 #endif
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: Store Hit at " << std::hex << aHit->getUnitID() << std::dec << " "
                               << aHit->getDepth() << " due to " << tkID << " in time " << time << " of energy "
                               << aHit->getEM() / CLHEP::GeV << " GeV (EM) and " << aHit->getHadr() / CLHEP::GeV
                               << " GeV (Hadr)";
 #endif
   return ok;
 }

◆ setDetUnitId()

uint32_t CaloSD::setDetUnitId ( const G4Step * step )

overridepure virtual

Implements SensitiveDetector.

Implemented in HCalSD, ECalSD, HcalTB02SD, CastorSD, HGCSD, DreamSD, HGCalSD, HFNoseSD, HGCScintSD, EcalTBH4BeamSD, HcalTB06BeamSD, ZdcSD, HGCalTB16SD01, AHCalSD, and TotemT2ScintSD.

Referenced by ProcessHits().

◆ setNumberCheckedHits()

void CaloSD::setNumberCheckedHits	(	int	val,
		int	k = `0`
	)

inlineprotected

Definition at line 126 of file CaloSD.h.

References dqmdumpme::k, nCheckedHits, and heppy_batch::val.

Referenced by HFNoseSD::HFNoseSD(), HGCalSD::HGCalSD(), HGCScintSD::HGCScintSD(), and HGCSD::HGCSD().

126 { nCheckedHits[k] = val; }

CaloSD::nCheckedHits

int nCheckedHits[2]

Definition: CaloSD.h:187

heppy_batch.val

val

Definition: heppy_batch.py:351

dqmdumpme.k

k

Definition: dqmdumpme.py:60

◆ setParameterized()

void CaloSD::setParameterized ( bool val )

inlineprotected

Definition at line 114 of file CaloSD.h.

References isParameterized, and heppy_batch::val.

Referenced by CastorSD::CastorSD(), HCalSD::HCalSD(), and ZdcSD::ZdcSD().

114 { isParameterized = val; }

CaloSD::isParameterized

bool isParameterized

Definition: CaloSD.h:177

heppy_batch.val

val

Definition: heppy_batch.py:351

◆ setToGlobal()

G4ThreeVector CaloSD::setToGlobal	(	const G4ThreeVector &	local,
		const G4VTouchable *	touch
	)		const

protected

Definition at line 458 of file CaloSD.cc.

References DTRecHitClients_cfi::local.

                                                                                              {
   return touch->GetHistory()->GetTopTransform().Inverse().TransformPoint(local);
 }

◆ setToLocal()

G4ThreeVector CaloSD::setToLocal	(	const G4ThreeVector &	global,
		const G4VTouchable *	touch
	)		const

protected

Definition at line 454 of file CaloSD.cc.

Referenced by ZdcSD::calculateCherenkovDeposit(), DreamSD::curve_LY(), HcalTB02SD::curve_LY(), ECalSD::EnergyCorrected(), ECalSD::getDepth(), ECalSD::getRadiationLength(), resetForNewPrimary(), and HcalTB06BeamSD::setDetUnitId().

                                                                                              {
   return touch->GetHistory()->GetTopTransform().TransformPoint(global);
 }

◆ setTrackID()

int CaloSD::setTrackID ( const G4Step * aStep )

protectedvirtual

Reimplemented in ZdcSD.

Definition at line 890 of file CaloSD.cc.

References SensitiveDetector::cmsTrackInformation(), TrackInformation::getIDonCaloSurface(), previousID, and resetForNewPrimary().

Referenced by HCalSD::getFromHFLibrary(), CastorSD::getFromLibrary(), HCalSD::getFromParam(), and HCalSD::hitForFibre().

                                           {
   auto const theTrack = aStep->GetTrack();
   TrackInformation* trkInfo = cmsTrackInformation(theTrack);
   int primaryID = trkInfo->getIDonCaloSurface();
   if (primaryID <= 0) {
     primaryID = theTrack->GetTrackID();
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "Track ID: " << trkInfo->getIDonCaloSurface() << ":" << theTrack->GetTrackID() << ":"
                               << primaryID;
 #endif
 
   if (primaryID != previousID[0].trackID()) {
     resetForNewPrimary(aStep);
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD::setTrackID for " << GetName()
                               << " trackID= " << aStep->GetTrack()->GetTrackID() << " primaryID= " << primaryID;
 #endif
   return primaryID;
 }

◆ setUseMap()

void CaloSD::setUseMap ( bool val )

inlineprotected

Definition at line 115 of file CaloSD.h.

References useMap, and heppy_batch::val.

Referenced by HFNoseSD::HFNoseSD(), HGCalSD::HGCalSD(), HGCScintSD::HGCScintSD(), and HGCSD::HGCSD().

115 { useMap = val; }

CaloSD::useMap

bool useMap

Definition: CaloSD.h:179

heppy_batch.val

val

Definition: heppy_batch.py:351

◆ shortreprID() [1/2]

std::string CaloSD::shortreprID ( const CaloHitID & ID )

protected

Definition at line 537 of file CaloSD.cc.

References contentValuesCheck::ss.

Referenced by createNewHit(), ProcessHits(), saveHit(), and shortreprID().

                                                  {
   std::stringstream ss;
   ss << GetName() << "/" << ID.unitID() << "/trk" << ID.trackID() << "/d" << ID.depth() << "/time" << ID.timeSliceID();
   if (ID.isFinecaloTrackID())
     ss << "/FC";
   return ss.str();
 }

◆ shortreprID() [2/2]

std::string CaloSD::shortreprID ( const CaloG4Hit * hit )

protected

Definition at line 546 of file CaloSD.cc.

References shortreprID().

546 { return shortreprID(hit->getID()); }

hit

Definition: SiStripHitEffFromCalibTree.cc:87

CaloSD::shortreprID

std::string shortreprID(const CaloHitID &ID)

Definition: CaloSD.cc:537

◆ storeHit()

void CaloSD::storeHit	(	CaloG4Hit *	hit,
		int	k = `0`
	)

private

Definition at line 936 of file CaloSD.cc.

References hitMap, dqmdumpme::k, previousID, theHC, CaloHitID::trackID(), and useMap.

Referenced by createNewHit().

                                            {
   if (hit == nullptr || previousID[k].trackID() < 0) {
     edm::LogWarning("CaloSim") << "CaloSD: hit to be stored is nullptr !!"
                                << " previousID.trackID()= " << previousID[k].trackID();
     return;
   }
 
   theHC[k]->insert(hit);
   if (useMap)
     hitMap[k].insert(std::pair<CaloHitID, CaloG4Hit*>(previousID[k], hit));
 }

◆ update() [1/5]

void CaloSD::update ( const BeginOfRun * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 748 of file CaloSD.cc.

References initRun().

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

748 { initRun(); }

CaloSD::initRun

virtual void initRun()

Definition: CaloSD.cc:866

◆ update() [2/5]

void CaloSD::update ( const BeginOfEvent * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 750 of file CaloSD.cc.

References clearHits(), and initEvent().

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

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

◆ update() [3/5]

void CaloSD::update ( const BeginOfTrack * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack *>.

Definition at line 1012 of file CaloSD.cc.

References clean(), cleanHitCollection(), SensitiveDetector::cmsTrackInformation(), TrackInformation::isPrimary(), dqmdumpme::k, nHC_, primAncestor, and theHC.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), 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() [4/5]

void CaloSD::update ( const EndOfTrack * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack *>.

Definition at line 758 of file CaloSD.cc.

References SensitiveDetector::cmsTrackInformation(), TrackInformation::getIDonCaloSurface(), l1ctLayer2EG_cff::id, m_trackManager, tkMap, SimTrackManager::trackContainer(), and TrackWithHistory::trackID().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), 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() [5/5]

void CaloSD::update ( const ::EndOfEvent * )

overrideprotected

Definition at line 779 of file CaloSD.cc.

References funct::abs(), boundaryCrossingParentMap_, submitPVResolutionJobs::count, mps_fire::end, endEvent(), edm::roottree::getEntry(), hitMap, mps_fire::i, dqmdumpme::k, nHC_, reusehit, saveHit(), slave, mathSSE::sqrt(), edm::swap(), theHC, tkMap, totalHits, groupFilesInBlocks::tt, useMap, x, and z.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), 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();
 }

◆ updateHit()

void CaloSD::updateHit	(	CaloG4Hit *	aHit,
		int	k
	)

protected

Definition at line 702 of file CaloSD.cc.

References CaloG4Hit::addEnergyDeposit(), currentID, edepositEM, edepositHAD, dqmdumpme::k, and previousID.

Referenced by checkHit(), createNewHit(), hitExists(), processHit(), and ProcessHits().

                                              {
   aHit->addEnergyDeposit(edepositEM, edepositHAD);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD:" << GetName() << " Add energy deposit in " << currentID[k]
                               << " Edep_em(MeV)= " << edepositEM << " Edep_had(MeV)= " << edepositHAD;
 #endif
 
   // buffer for next steps:
   previousID[k] = currentID[k];
 }

Member Data Documentation

◆ boundaryCrossingParentMap_

std::unordered_map<unsigned int, unsigned int> CaloSD::boundaryCrossingParentMap_

private

Definition at line 197 of file CaloSD.h.

Referenced by findBoundaryCrossingParent(), and update().

◆ cleanIndex

int CaloSD::cleanIndex[2]

private

Definition at line 184 of file CaloSD.h.

Referenced by CaloSD(), cleanHitCollection(), and clearHits().

◆ collName_

std::string CaloSD::collName_[2]

protected

Definition at line 159 of file CaloSD.h.

Referenced by CaloSD(), Initialize(), newCollection(), and HGCalSD::update().

◆ correctT

double CaloSD::correctT

private

Definition at line 191 of file CaloSD.h.

Referenced by CaloSD(), cleanHitCollection(), and saveHit().

◆ corrTOFBeam

bool CaloSD::corrTOFBeam

private

Definition at line 180 of file CaloSD.h.

Referenced by CaloSD(), cleanHitCollection(), and saveHit().

◆ currentHit

CaloG4Hit* CaloSD::currentHit[2]

protected

Definition at line 152 of file CaloSD.h.

Referenced by CaloSD(), checkHit(), hitExists(), processHit(), ZdcSD::ProcessHits(), ProcessHits(), and HGCalSD::processSecondHit().

◆ currentID

CaloHitID CaloSD::currentID[2]

protected

Definition at line 146 of file CaloSD.h.

Referenced by CaloSD(), checkHit(), createNewHit(), ECalSD::curve_LY(), HCalSD::getFromHFLibrary(), CastorSD::getFromLibrary(), HCalSD::getFromParam(), HCalSD::getHitFibreBundle(), HCalSD::getHitPMT(), hitExists(), HCalSD::hitForFibre(), processHit(), ZdcSD::ProcessHits(), ProcessHits(), HGCalSD::processSecondHit(), and updateHit().

◆ detName_

std::string CaloSD::detName_

protected

Definition at line 159 of file CaloSD.h.

Referenced by CaloSD(), Initialize(), and newCollection().

◆ doFineCalo_

bool CaloSD::doFineCalo_

private

Definition at line 192 of file CaloSD.h.

Referenced by CaloSD(), ProcessHits(), and saveHit().

◆ doFineCaloThisStep_

bool CaloSD::doFineCaloThisStep_

private

Definition at line 200 of file CaloSD.h.

Referenced by createNewHit(), and ProcessHits().

◆ edepositEM

float CaloSD::edepositEM

protected

Definition at line 144 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), HCalSD::getFromHFLibrary(), CastorSD::getFromLibrary(), HCalSD::getFromParam(), HCalSD::getHitFibreBundle(), HCalSD::getHitPMT(), HCalSD::hitForFibre(), ZdcSD::ProcessHits(), ProcessHits(), HGCalSD::processSecondHit(), and updateHit().

◆ edepositHAD

float CaloSD::edepositHAD

protected

Definition at line 144 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), HCalSD::getFromHFLibrary(), CastorSD::getFromLibrary(), HCalSD::getFromParam(), HCalSD::getHitFibreBundle(), HCalSD::getHitPMT(), HCalSD::hitForFibre(), ZdcSD::ProcessHits(), ProcessHits(), HGCalSD::processSecondHit(), and updateHit().

◆ eMinFine_

double CaloSD::eMinFine_

private

Definition at line 193 of file CaloSD.h.

Referenced by CaloSD().

◆ eminHit

double CaloSD::eminHit

protected

Definition at line 148 of file CaloSD.h.

Referenced by CaloSD(), filterHit(), and HGCSD::HGCSD().

◆ eminHitD

double CaloSD::eminHitD

private

Definition at line 190 of file CaloSD.h.

Referenced by CaloSD(), and filterHit().

◆ energyCut

double CaloSD::energyCut

protected

Definition at line 148 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), HCalSD::getHitFibreBundle(), and HCalSD::getHitPMT().

◆ entranceLocal

G4ThreeVector CaloSD::entranceLocal

protected

Definition at line 141 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), ProcessHits(), and resetForNewPrimary().

◆ entrancePoint

G4ThreeVector CaloSD::entrancePoint

protected

Definition at line 140 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), ProcessHits(), and resetForNewPrimary().

◆ fineDetectors_

std::vector<Detector> CaloSD::fineDetectors_

private

Definition at line 199 of file CaloSD.h.

Referenced by CaloSD(), and isItFineCalo().

◆ forceSave

bool CaloSD::forceSave

protected

Definition at line 157 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), and ECalSD::getTrackID().

◆ hcID

int CaloSD::hcID[2]

private

Definition at line 182 of file CaloSD.h.

Referenced by CaloSD(), and Initialize().

◆ hcn_

std::vector<std::string> CaloSD::hcn_

protected

Definition at line 149 of file CaloSD.h.

Referenced by CaloSD(), and newCollection().

◆ hitMap

std::map<CaloHitID, CaloG4Hit*> CaloSD::hitMap[2]

private

Definition at line 195 of file CaloSD.h.

Referenced by checkHit(), cleanHitCollection(), storeHit(), and update().

◆ ignoreReject

bool CaloSD::ignoreReject

private

Definition at line 178 of file CaloSD.h.

Referenced by ignoreRejection(), and ProcessHits().

◆ ignoreTrackID

bool CaloSD::ignoreTrackID

private

Definition at line 176 of file CaloSD.h.

Referenced by CaloSD().

◆ incidentEnergy

float CaloSD::incidentEnergy

protected

Definition at line 143 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), ProcessHits(), and resetForNewPrimary().

◆ isParameterized

bool CaloSD::isParameterized

private

Definition at line 177 of file CaloSD.h.

Referenced by CaloSD(), ProcessHits(), and setParameterized().

◆ kmaxIon

double CaloSD::kmaxIon

protected

Definition at line 155 of file CaloSD.h.

Referenced by CaloSD(), ECalSD::ECalSD(), ECalSD::getEnergyDeposit(), HCalSD::getEnergyDeposit(), and HCalSD::HCalSD().

◆ kmaxNeutron

double CaloSD::kmaxNeutron

protected

Definition at line 155 of file CaloSD.h.

Referenced by CaloSD(), ECalSD::ECalSD(), ECalSD::getEnergyDeposit(), HCalSD::getEnergyDeposit(), and HCalSD::HCalSD().

◆ kmaxProton

double CaloSD::kmaxProton

protected

Definition at line 155 of file CaloSD.h.

Referenced by CaloSD(), ECalSD::ECalSD(), ECalSD::getEnergyDeposit(), HCalSD::getEnergyDeposit(), and HCalSD::HCalSD().

◆ m_trackManager

const SimTrackManager* CaloSD::m_trackManager

private

Definition at line 169 of file CaloSD.h.

Referenced by findBoundaryCrossingParent(), saveHit(), and update().

◆ meanResponse

std::unique_ptr<CaloMeanResponse> CaloSD::meanResponse[2]

private

Definition at line 172 of file CaloSD.h.

Referenced by CaloSD(), getResponseWt(), and newCollection().

◆ nCheckedHits

int CaloSD::nCheckedHits[2]

private

Definition at line 187 of file CaloSD.h.

Referenced by CaloSD(), checkHit(), and setNumberCheckedHits().

◆ nHC_

int CaloSD::nHC_

protected

Definition at line 158 of file CaloSD.h.

Referenced by CaloSD(), cleanHitCollection(), clearHits(), EndOfEvent(), fillHits(), HGCalSD::HGCalSD(), Initialize(), newCollection(), PrintAll(), HGCalSD::update(), and update().

◆ posGlobal

G4ThreeVector CaloSD::posGlobal

protected

Definition at line 142 of file CaloSD.h.

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

◆ previousID

CaloHitID CaloSD::previousID[2]

protected

Definition at line 146 of file CaloSD.h.

Referenced by CaloSD(), clearHits(), hitExists(), processHit(), ProcessHits(), ZdcSD::setTrackID(), setTrackID(), storeHit(), and updateHit().

◆ primAncestor

int CaloSD::primAncestor

private

Definition at line 183 of file CaloSD.h.

Referenced by CaloSD(), and update().

◆ primIDSaved

int CaloSD::primIDSaved[2]

private

Definition at line 186 of file CaloSD.h.

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

◆ reusehit

std::vector<std::unique_ptr<CaloG4Hit> > CaloSD::reusehit[2]

private

Definition at line 198 of file CaloSD.h.

Referenced by cleanHitCollection(), createNewHit(), and update().

◆ slave

std::unique_ptr<CaloSlaveSD> CaloSD::slave[2]

private

Definition at line 171 of file CaloSD.h.

Referenced by CaloSD(), clearHits(), fillHits(), newCollection(), saveHit(), and update().

◆ suppressHeavy

bool CaloSD::suppressHeavy

protected

Definition at line 154 of file CaloSD.h.

Referenced by CaloSD(), ECalSD::ECalSD(), ECalSD::getEnergyDeposit(), HCalSD::getEnergyDeposit(), and HCalSD::HCalSD().

◆ theHC

CaloG4HitCollection* CaloSD::theHC[2]

private

Definition at line 174 of file CaloSD.h.

Referenced by CaloSD(), checkHit(), cleanHitCollection(), EndOfEvent(), getNumberOfHits(), Initialize(), PrintAll(), storeHit(), and update().

◆ timeSlice

float CaloSD::timeSlice

private

Definition at line 189 of file CaloSD.h.

Referenced by CaloSD().

◆ tkMap

std::map<int, TrackWithHistory*> CaloSD::tkMap

private

Definition at line 196 of file CaloSD.h.

Referenced by createNewHit(), and update().

◆ tmaxHit

double CaloSD::tmaxHit

protected

Definition at line 148 of file CaloSD.h.

Referenced by CaloSD(), AHCalSD::filterHit(), HGCScintSD::filterHit(), HFNoseSD::filterHit(), HGCalSD::filterHit(), HGCSD::filterHit(), HCalSD::filterHit(), and filterHit().

◆ totalHits

int CaloSD::totalHits[2]

private

Definition at line 185 of file CaloSD.h.

Referenced by CaloSD(), cleanHitCollection(), createNewHit(), Initialize(), and update().

◆ useMap

bool CaloSD::useMap

private

Definition at line 179 of file CaloSD.h.

Referenced by CaloSD(), checkHit(), cleanHitCollection(), createNewHit(), setUseMap(), storeHit(), and update().

◆ useResMap_

std::vector<int> CaloSD::useResMap_

protected

Definition at line 150 of file CaloSD.h.

Referenced by CaloSD(), and newCollection().

Classes

Public Member Functions

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ CaloSD()

◆ ~CaloSD()

Member Function Documentation

◆ checkHit()

◆ cleanHitCollection()

◆ clear()

◆ clearHits()

◆ createNewHit()

◆ DrawAll()

◆ endEvent()

◆ EndOfEvent()

◆ EnergyCorrected()

◆ fillHits()

◆ filterHit()

◆ findBoundaryCrossingParent()

◆ getAttenuation()

◆ getDepth()

◆ getEnergyDeposit()

◆ getFromLibrary()

◆ getNumberOfHits()

◆ getResponseWt()

◆ getTrackID()

◆ hitExists()

◆ ignoreRejection()

◆ initEvent()

◆ Initialize()

◆ initRun()

◆ isItFineCalo()

◆ newCollection()

◆ printableDecayChain()

◆ PrintAll()

◆ printDetectorLevels()

◆ processHit()

◆ ProcessHits() [1/2]

◆ ProcessHits() [2/2]

◆ processSecondHit()

◆ reset()

◆ resetForNewPrimary()

◆ saveHit()

◆ setDetUnitId()

◆ setNumberCheckedHits()

◆ setParameterized()

◆ setToGlobal()

◆ setToLocal()

◆ setTrackID()

◆ setUseMap()

◆ shortreprID() [1/2]

◆ shortreprID() [2/2]

◆ storeHit()

◆ update() [1/5]

◆ update() [2/5]

◆ update() [3/5]

◆ update() [4/5]

◆ update() [5/5]

◆ updateHit()

Member Data Documentation

◆ boundaryCrossingParentMap_

◆ cleanIndex

◆ collName_

◆ correctT

◆ corrTOFBeam

◆ currentHit

◆ currentID

◆ detName_

◆ doFineCalo_

◆ doFineCaloThisStep_

◆ edepositEM

◆ edepositHAD

◆ eMinFine_

◆ eminHit