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

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

bool	ProcessHits (G4GFlashSpot aSpot, G4TouchableHistory ) override

G4bool	ProcessHits (G4Step step, 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)

Public Member Functions inherited from SensitiveDetector
void	EndOfEvent (G4HCofThisEvent *eventHC) override

const std::vector< std::string > &	getNames () const

void	Initialize (G4HCofThisEvent *eventHC) override

bool	isCaloSD () const

	SensitiveDetector (const std::string &iname, const SensitiveDetectorCatalog &, bool calo)

	~SensitiveDetector () override

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

CaloG4Hit *	createNewHit (const G4Step , const G4Track )

virtual void	endEvent ()

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

virtual bool	filterHit (CaloG4Hit *, double)

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

double	getResponseWt (const G4Track *)

virtual int	getTrackID (const G4Track *)

void	hitBookkeepingFineCalo (const G4Step step, const G4Track currentTrack, CaloG4Hit *hit)

bool	hitExists (const G4Step *)

void	ignoreRejection ()

virtual void	initEvent (const BeginOfEvent *)

virtual void	initRun ()

void	printDetectorLevels (const G4VTouchable *) const

void	processHit (const G4Step *step)

void	resetForNewPrimary (const G4Step *)

void	setNumberCheckedHits (int val)

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)

void	update (const ::EndOfEvent *) override

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

void	update (const BeginOfRun *) 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	updateHit (CaloG4Hit *)

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
CaloG4Hit *	currentHit

CaloHitID	currentID

float	edepositEM

float	edepositHAD

double	eminHit

double	energyCut

G4ThreeVector	entranceLocal

G4ThreeVector	entrancePoint

bool	forceSave

float	incidentEnergy

double	kmaxIon

double	kmaxNeutron

double	kmaxProton

G4ThreeVector	posGlobal

CaloHitID	previousID

bool	suppressHeavy

double	tmaxHit

Private Member Functions
void	cleanHitCollection ()

bool	saveHit (CaloG4Hit *)

void	storeHit (CaloG4Hit *)

Private Attributes
int	cleanIndex

double	correctT

bool	corrTOFBeam

bool	doFineCalo_

double	eMinFine_

double	eminHitD

std::vector< Detector >	fineDetectors_

int	hcID

std::map< CaloHitID, CaloG4Hit * >	hitMap

bool	ignoreReject

bool	ignoreTrackID

bool	isParameterized

const SimTrackManager *	m_trackManager

std::unique_ptr< CaloMeanResponse >	meanResponse

int	nCheckedHits

int	primAncestor

int	primIDSaved

std::vector< std::unique_ptr< CaloG4Hit > >	reusehit

std::unique_ptr< CaloSlaveSD >	slave

CaloG4HitCollection *	theHC

float	timeSlice

std::map< int, TrackWithHistory * >	tkMap

int	totalHits

bool	useMap

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

Detailed Description

Definition at line 39 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`
	)

Definition at line 34 of file CaloSD.cc.

     : SensitiveCaloDetector(name, clg),
       G4VGFlashSensitiveDetector(),
       eminHit(0.),
       currentHit(nullptr),
       m_trackManager(manager),
       theHC(nullptr),
       ignoreTrackID(ignoreTkID),
       hcID(-1),
       timeSlice(timeSliceUnit),
       eminHitD(0.) {
   //Parameters
   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");
   std::vector<std::string> hcn = m_CaloSD.getParameter<std::vector<std::string>>("HCNames");
   std::vector<int> 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 = m_CaloSD.getUntrackedParameter<int>("CheckHits", 25);
   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);
   meanResponse.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 = std::make_unique<CaloMeanResponse>(p);
         break;
       }
     }
   }
   slave = std::make_unique<CaloSlaveSD>(name);
  
   currentID = CaloHitID(timeSlice, ignoreTrackID);
   previousID = 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 = cleanIndex = totalHits = primIDSaved = 0;
   forceSave = false;
  
   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
                               << " 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";
  
   // 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
 }

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

◆ ~CaloSD()

CaloSD::~CaloSD ( )

override

Definition at line 157 of file CaloSD.cc.

157 {}

Member Function Documentation

◆ checkHit()

bool CaloSD::checkHit ( )

protected

Definition at line 414 of file CaloSD.cc.

                       {
   //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.find(currentID);
     if (it != hitMap.end()) {
       currentHit = it->second;
       found = true;
     }
   } else if (nCheckedHits > 0) {
     int nhits = theHC->entries();
     int minhit = std::max(nhits - nCheckedHits, 0);
     int maxhit = nhits - 1;
  
     for (int j = maxhit; j > minhit; --j) {
       if ((*theHC)[j]->getID() == currentID) {
         currentHit = (*theHC)[j];
         found = true;
         break;
       }
     }
   }
  
   if (found) {
     updateHit(currentHit);
   }
   return found;
 }

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

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

◆ cleanHitCollection()

void CaloSD::cleanHitCollection ( )

private

Definition at line 986 of file CaloSD.cc.

                                 {
   std::vector<CaloG4Hit*>* theCollection = theHC->GetVector();
  
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: collection before merging, size = " << theHC->entries();
 #endif
   if (reusehit.empty())
     reusehit.reserve(theHC->entries() - cleanIndex);
  
   // 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), hitvec.end(), CaloG4HitLess());
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD::cleanHitCollection: sort hits in buffer starting from "
                                 << "element = " << cleanIndex;
     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; 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.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, 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 = theHC->entries();
   }
  
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: collection after merging, size= " << theHC->entries()
                               << " Size of reusehit= " << reusehit.size()
                               << "\n      starting hit selection from index = " << cleanIndex;
 #endif
  
   int addhit = 0;
   for (unsigned int i = cleanIndex; 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.emplace_back((*theCollection)[i]);
       (*theCollection)[i] = nullptr;
       ++addhit;
     }
   }
  
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: Size of reusehit after selection = " << reusehit.size()
                               << " Number of added hit = " << addhit;
 #endif
   if (useMap) {
     if (addhit > 0) {
       int offset = reusehit.size() - addhit;
       for (int ii = addhit - 1; ii >= 0; --ii) {
         CaloHitID theID = reusehit[offset + ii]->getID();
         hitMap.erase(theID);
       }
     }
   }
  
   //move all nullptr to end of list and then remove them
   theCollection->erase(
       std::stable_partition(
           theCollection->begin() + cleanIndex, theCollection->end(), [](CaloG4Hit* p) { return p != nullptr; }),
       theCollection->end());
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: hit collection after selection, size = " << theHC->entries();
   theHC->PrintAllHits();
 #endif
  
   cleanIndex = theHC->entries();
 }

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

Referenced by EndOfEvent(), and update().

◆ clear()

void CaloSD::clear ( void )

override

Definition at line 363 of file CaloSD.cc.

363 {}

◆ clearHits()

void CaloSD::clearHits ( )

overridevirtual

Implements SensitiveDetector.

Definition at line 798 of file CaloSD.cc.

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

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

Referenced by update().

◆ createNewHit()

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

protected

Definition at line 566 of file CaloSD.cc.

                                                                             {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD::CreateNewHit " << getNumberOfHits() << " for " << GetName()
                               << " Unit:" << currentID.unitID() << " " << currentID.depth() << " Edep= " << edepositEM
                               << " " << edepositHAD << " primaryID= " << currentID.trackID()
                               << " timeSlice= " << currentID.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.empty()) {
     aHit = reusehit.back().release();
     aHit->setEM(0.f);
     aHit->setHadr(0.f);
     reusehit.pop_back();
   } else {
     aHit = new CaloG4Hit;
   }
  
   aHit->setID(currentID);
   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);
  
   storeHit(aHit);
   TrackInformation* trkInfo = cmsTrackInformation(theTrack);
  
   bool currentlyInsideFineVolume = !doFineCalo_ ? false : isItFineCalo(aStep->GetPostStepPoint()->GetTouchable());
  
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("DoFineCalo") << "Creating new hit " << aHit->getUnitID() << " using "
                                  << (currentlyInsideFineVolume ? "FINECALO" : "normal CaloSD")
                                  << "; currentID.trackID=" << currentID.trackID()
                                  << " currentID.fineTrackID=" << currentID.fineTrackID()
                                  << " isItFineCalo(aStep->GetPostStepPoint()->GetTouchable())="
                                  << isItFineCalo(aStep->GetPostStepPoint()->GetTouchable())
                                  << " isItFineCalo(aStep->GetPreStepPoint()->GetTouchable())="
                                  << isItFineCalo(aStep->GetPreStepPoint()->GetTouchable())
                                  << " theTrack=" << theTrack->GetTrackID() << " ("
                                  << " parentTrackId=" << theTrack->GetParentID()
                                  << " getIDonCaloSurface=" << trkInfo->getIDonCaloSurface() << ")"
                                  << " primIDSaved=" << primIDSaved;
 #endif
  
   // If fine calo is activated for the current volume, perform track/hit
   // saving logic for fineCalo
   if (currentlyInsideFineVolume) {
     hitBookkeepingFineCalo(aStep, theTrack, aHit);
   }
   // 'Traditional', non-fine history bookkeeping
   else {
     double etrack = 0;
     if (currentID.trackID() == primIDSaved) {  // The track is saved; nothing to be done
     } else if (currentID.trackID() == theTrack->GetTrackID()) {
       etrack = theTrack->GetKineticEnergy();
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("CaloSim") << "CaloSD: set save the track " << currentID.trackID() << " etrack " << etrack
                                   << " eCut " << energyCut << " force: " << forceSave
                                   << " save: " << (etrack >= energyCut || forceSave);
 #endif
       if (etrack >= energyCut || forceSave) {
         trkInfo->storeTrack(true);
         trkInfo->putInHistory();
       }
     } else {
       TrackWithHistory* trkh = tkMap[currentID.trackID()];
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("CaloSim") << "CaloSD : TrackWithHistory pointer for " << currentID.trackID() << " is " << trkh;
 #endif
       if (trkh != nullptr) {
         etrack = sqrt(trkh->momentum().Mag2());
         if (etrack >= energyCut) {
           trkh->save();
 #ifdef EDM_ML_DEBUG
           edm::LogVerbatim("CaloSim") << "CaloSD: set save the track " << currentID.trackID() << " with Hit";
 #endif
         }
       }
     }
     primIDSaved = currentID.trackID();
   }
  
   if (useMap)
     ++totalHits;
   return aHit;
 }

References SensitiveDetector::cmsTrackInformation(), currentID, CaloHitID::depth(), doFineCalo_, edepositEM, edepositHAD, energyCut, entranceLocal, entrancePoint, f, funct::false, CaloHitID::fineTrackID(), forceSave, TrackInformation::getIDonCaloSurface(), getNumberOfHits(), CaloG4Hit::getUnitID(), GeV, hitBookkeepingFineCalo(), incidentEnergy, isItFineCalo(), TrackWithHistory::momentum(), posGlobal, primIDSaved, TrackInformation::putInHistory(), reusehit, TrackWithHistory::save(), CaloG4Hit::setEM(), CaloG4Hit::setEntry(), CaloG4Hit::setEntryLocal(), CaloG4Hit::setHadr(), CaloG4Hit::setID(), CaloG4Hit::setIncidentEnergy(), CaloG4Hit::setPosition(), mathSSE::sqrt(), storeHit(), TrackInformation::storeTrack(), CaloHitID::timeSliceID(), tkMap, totalHits, CaloHitID::trackID(), CaloHitID::unitID(), updateHit(), and useMap.

Referenced by processHit(), and ProcessHits().

◆ DrawAll()

void CaloSD::DrawAll ( )

override

Definition at line 365 of file CaloSD.cc.

365 {}

◆ endEvent()

void CaloSD::endEvent ( )

protectedvirtual

Reimplemented in HGCSD.

Definition at line 819 of file CaloSD.cc.

819 {}

Referenced by update().

◆ EndOfEvent()

void CaloSD::EndOfEvent ( G4HCofThisEvent * eventHC )

override

Definition at line 350 of file CaloSD.cc.

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

References cleanHitCollection(), and theHC.

◆ EnergyCorrected()

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

protectedvirtual

Reimplemented in ECalSD.

Definition at line 308 of file CaloSD.cc.

308 { return aStep.GetTotalEnergyDeposit(); }

Referenced by ProcessHits().

◆ fillHits()

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

overridevirtual

Implements SensitiveCaloDetector.

Definition at line 374 of file CaloSD.cc.

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

References slave.

◆ filterHit()

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

protectedvirtual

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

Definition at line 863 of file CaloSD.cc.

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

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

Referenced by cleanHitCollection().

◆ getAttenuation()

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

protected

Definition at line 683 of file CaloSD.cc.

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

References funct::abs(), 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().

◆ getDepth()

uint16_t CaloSD::getDepth ( const G4Step * )

protectedvirtual

Reimplemented in ECalSD.

Definition at line 861 of file CaloSD.cc.

861 { return 0; }

Referenced by ProcessHits().

◆ getEnergyDeposit()

double CaloSD::getEnergyDeposit ( const G4Step * step )

protectedvirtual

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

Definition at line 306 of file CaloSD.cc.

306 { return aStep->GetTotalEnergyDeposit(); }

Referenced by ProcessHits().

◆ getFromLibrary()

bool CaloSD::getFromLibrary ( const G4Step * step )

protectedvirtual

Reimplemented in HCalSD, CastorSD, and ZdcSD.

Definition at line 310 of file CaloSD.cc.

310 { return false; }

Referenced by ProcessHits().

◆ getNumberOfHits()

int CaloSD::getNumberOfHits ( )

protected

Definition at line 443 of file CaloSD.cc.

443 { return theHC->entries(); }

References theHC.

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

◆ getResponseWt()

double CaloSD::getResponseWt ( const G4Track * aTrack )

protected

Definition at line 874 of file CaloSD.cc.

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

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

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

◆ getTrackID()

int CaloSD::getTrackID ( const G4Track * aTrack )

protectedvirtual

Reimplemented in ECalSD.

Definition at line 821 of file CaloSD.cc.

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

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

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

◆ hitBookkeepingFineCalo()

void CaloSD::hitBookkeepingFineCalo	(	const G4Step *	step,
		const G4Track *	currentTrack,
		CaloG4Hit *	hit
	)

protected

Definition at line 459 of file CaloSD.cc.

                                                                                                    {
   /*
   Performs bookkeeping: Determines what trackIDs are to be recorded for the hit (typically some
   some parent trackID), and also sets the right flags on either the TrackInformation object or 
   the TrackWithHistory object to make sure the right track is saved to the SimTrack collection.
  
   `currentTrack` is the track that is currently being processed by Geant
   */
   TrackInformation* trkInfo = cmsTrackInformation(currentTrack);
   // Copy the  class's currentID so we can freely modify it without influencing
   // hits created later in possibly non-fine detectors by the same track
   CaloHitID hitID = currentID;
   // First check if the current currentTrack passes criteria
   if (trkInfo->crossedBoundary()) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("DoFineCalo") << "currentTrack " << currentTrack->GetTrackID()
                                    << " itself has crossedBoundary=" << trkInfo->crossedBoundary()
                                    << " ; recording it for hit " << hit->getUnitID();
 #endif
     hitID.setFineTrackID(currentTrack->GetTrackID());
     hit->setID(hitID);  // Actually overwrite the ID for the hit
     trkInfo->storeTrack(true);
     return;
   }
   // currentTrack itself does not pass thresholds / does not cross boundary; go through its history to find a track that does
   TrackWithHistory* recordTrackWithHistory;
   // Keep track of decay chain of this track for debugging purposes
   std::vector<unsigned int> decayChain;
   decayChain.push_back(currentTrack->GetTrackID());
   // Find the first parent of this track that passes the required criteria
   // Start from first parent
   unsigned int recordTrackID = currentTrack->GetParentID();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("DoFineCalo") << "Trying to find the first parent of hit " << hit->getUnitID()
                                  << " that passes saving criterion (crosses boundary or specific criterion)"
                                  << "; starting with first parent track " << recordTrackID;
 #endif
   // Check whether this first parent actually exists
   if (recordTrackID <= 0) {
     // Track ID 0 is not a track;
     // This means the current currentTrack has no parent, but apparently it also didn't fit saving criteria
     throw cms::Exception("Unknown", "CaloSD") << "ERROR: Track " << currentTrack->GetTrackID()
                                               << " has no parent, does not fit saving criteria, but left hit "
                                               << hit->getUnitID() << "; recording it but it's weird!";
   }
   // Start progressing through the track's history
   while (true) {
     // Record the decay chain for debugging purposes
     decayChain.push_back(recordTrackID);
     recordTrackWithHistory = m_trackManager->getTrackByID(recordTrackID);
     if (recordTrackID < (unsigned int)hitID.trackID()) {
       // A parent of the currentTrack has a lower trackID than the current
       // hitID.trackID(). This means the current hitID.trackID() does not point
       // to the earliest ancestor of the currentTrack.
       // The current hitID.trackID() might not be a saved track yet, but the
       // ancestor is *always* a saved track.
       // Fix this by overwriting the hitID's track ID.
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("DoFineCalo") << "History-tracking progressed to track " << recordTrackID
                                      << ", which is an earlier ancestor than current primary " << hitID.trackID()
                                      << "; overwriting it.";
 #endif
       hitID.setTrackID(recordTrackID);
     }
     // Check if this parent fits the boundary-crossing criteria
     if (recordTrackWithHistory->crossedBoundary() && recordTrackWithHistory->getIDAtBoundary() == (int)recordTrackID) {
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("DoFineCalo") << "Recording track " << recordTrackID << " as source of hit " << hit->getUnitID()
                                      << "; crossed boundary at pos=("
                                      << recordTrackWithHistory->getPositionAtBoundary().x() << ","
                                      << recordTrackWithHistory->getPositionAtBoundary().y() << ","
                                      << recordTrackWithHistory->getPositionAtBoundary().z() << ")"
                                      << " mom=(" << recordTrackWithHistory->getMomentumAtBoundary().x() << ","
                                      << recordTrackWithHistory->getMomentumAtBoundary().y() << ","
                                      << recordTrackWithHistory->getMomentumAtBoundary().z() << ","
                                      << recordTrackWithHistory->getMomentumAtBoundary().e() << ")"
                                      << " id@boundary=" << recordTrackWithHistory->getIDAtBoundary()
                                      << "; decayChain: " << printableDecayChain(decayChain);
 #endif
       break;
     }
     // This parent track did not fit criteria - go to the next parent
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("DoFineCalo") << "Track " << recordTrackID << " did not cross the boundary or fit other criteria";
 #endif
     recordTrackID = recordTrackWithHistory->parentID();
     if (recordTrackID <= 0) {
       // Track ID 0 is not a track;
       // This means that no parent of the currentTrack fitted the criteria
       throw cms::Exception("Unknown", "CaloSD")
           << "Hit " << hit->getUnitID() << " does not have any parent track that passes the criteria!"
           << " decayChain so far: " << printableDecayChain(decayChain);
     }
   }
   // Parentage traversal done - do the bookeeping for the found ancestor track
   recordTrackWithHistory->save();
   hitID.setFineTrackID(recordTrackID);
   hit->setID(hitID);  // Actually overwrite the ID for the hit
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("DoFineCalo") << "Stored the following bookeeping for hit " << hit->getUnitID()
                                  << " hitID.trackID()=" << hitID.trackID()
                                  << " hitID.fineTrackID()=" << hitID.fineTrackID()
                                  << " recordTrackWithHistory->trackID()=" << recordTrackWithHistory->trackID()
                                  << " recordTrackWithHistory->saved()=" << recordTrackWithHistory->saved();
 #endif
 }

References SensitiveDetector::cmsTrackInformation(), TrackWithHistory::crossedBoundary(), TrackInformation::crossedBoundary(), currentID, Exception, CaloHitID::fineTrackID(), TrackWithHistory::getIDAtBoundary(), TrackWithHistory::getMomentumAtBoundary(), TrackWithHistory::getPositionAtBoundary(), SimTrackManager::getTrackByID(), createfilelist::int, m_trackManager, TrackWithHistory::parentID(), printableDecayChain(), TrackWithHistory::save(), TrackWithHistory::saved(), CaloHitID::setFineTrackID(), CaloHitID::setTrackID(), TrackInformation::storeTrack(), CaloHitID::trackID(), and TrackWithHistory::trackID().

Referenced by createNewHit().

◆ hitExists()

bool CaloSD::hitExists ( const G4Step * aStep )

protected

Definition at line 393 of file CaloSD.cc.

                                           {
   // Update if in the same detector, time-slice and for same track
   if (currentID == previousID) {
     updateHit(currentHit);
     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
   posGlobal = aStep->GetPreStepPoint()->GetPosition();
   if (currentID.trackID() != previousID.trackID()) {
     resetForNewPrimary(aStep);
   }
   return checkHit();
 }

References checkHit(), currentHit, currentID, posGlobal, previousID, resetForNewPrimary(), CaloHitID::trackID(), and updateHit().

Referenced by ProcessHits().

◆ ignoreRejection()

void CaloSD::ignoreRejection ( )

inlineprotected

Definition at line 105 of file CaloSD.h.

105 { ignoreReject = true; }

References ignoreReject.

Referenced by HGCalSD::setDetUnitId().

◆ initEvent()

void CaloSD::initEvent ( const BeginOfEvent * )

protectedvirtual

Reimplemented in HGCSD.

Definition at line 817 of file CaloSD.cc.

817 {}

Referenced by update().

◆ Initialize()

void CaloSD::Initialize ( G4HCofThisEvent * HCE )

override

Definition at line 332 of file CaloSD.cc.

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

References bysipixelclustmulteventfilter_cfi::collectionName, hcID, theHC, and totalHits.

◆ initRun()

void CaloSD::initRun ( )

protectedvirtual

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

Definition at line 815 of file CaloSD.cc.

815 {}

Referenced by update().

◆ isItFineCalo()

bool CaloSD::isItFineCalo ( const G4VTouchable * touch )

Definition at line 312 of file CaloSD.cc.

                                                    {
   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 == 0) ? "Unknown" : lv->GetName();
       edm::LogVerbatim("CaloSim") << "CaloSD: volume " << name1 << ":" << detector.name << " at Level "
                                   << detector.level << " Flag " << ok;
 #endif
       if (ok)
         break;
     }
   }
   return ok;
 }

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

Referenced by createNewHit().

◆ printableDecayChain()

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

staticprotected

Definition at line 445 of file CaloSD.cc.

                                                                                {
   /*
   Takes a vector of ints (representing trackIDs), and returns a formatted string
   for debugging purposes
   */
   std::stringstream ss;
   for (long unsigned int i = 0; i < decayChain.size(); i++) {
     if (i > 0)
       ss << " <- ";
     ss << decayChain[i];
   }
   return ss.str();
 }

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

Referenced by hitBookkeepingFineCalo().

◆ PrintAll()

void CaloSD::PrintAll ( )

override

Definition at line 367 of file CaloSD.cc.

                       {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD: Collection " << theHC->GetName();
 #endif
   theHC->PrintAllHits();
 }

References theHC.

◆ printDetectorLevels()

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

protected

Definition at line 1102 of file CaloSD.cc.

                                                                 {
   //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();
 }

References mps_fire::i, cuy::ii, personalPlayback::level, Skims_PA_cff::name, MetAnalyzer::pv(), and AlCaHLTBitMon_QueryRunRegistry::string.

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

◆ processHit()

void CaloSD::processHit ( const G4Step * step )

inlineprotected

Definition at line 110 of file CaloSD.h.

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

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

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

◆ ProcessHits() [1/2]

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

override

Definition at line 241 of file CaloSD.cc.

                                                                  {
   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);
   if (edep <= 0.0) {
     return false;
   }
  
   if (G4TrackToParticleID::isGammaElectronPositron(track)) {
     edepositEM = edep;
   } else {
     edepositHAD = edep;
   }
  
   unsigned int unitID = setDetUnitId(&fFakeStep);
  
   if (unitID > 0) {
     double time = 0;
     int primaryID = getTrackID(track);
     uint16_t depth = getDepth(&fFakeStep);
     currentID.setID(unitID, time, primaryID, depth);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CaloSim") << "CaloSD:: GetSpotInfo for Unit 0x" << std::hex << currentID.unitID() << std::dec
                                 << " Edeposit = " << edepositEM << " " << edepositHAD;
 #endif
     // Update if in the same detector, time-slice and for same track
     if (currentID == previousID) {
       updateHit(currentHit);
     } else {
       posGlobal = aSpot->GetEnergySpot()->GetPosition();
       // Reset entry point for new primary
       if (currentID.trackID() != previousID.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 = createNewHit(&fFakeStep, track);
       }
     }
     return true;
   }
   return false;
 }

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

◆ ProcessHits() [2/2]

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

overridevirtual

Implements SensitiveDetector.

Definition at line 159 of file CaloSD.cc.

                                                              {
   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
   // apply shower library or parameterisation
   if (isParameterized) {
     if (getFromLibrary(aStep)) {
       // for parameterized showers the primary track should be killed
       aStep->GetTrack()->SetTrackStatus(fStopAndKill);
       auto tv = aStep->GetSecondary();
       auto vol = aStep->GetPreStepPoint()->GetPhysicalVolume();
       for (auto& tk : *tv) {
         if (tk->GetVolume() == vol) {
           tk->SetTrackStatus(fStopAndKill);
         }
       }
       return true;
     }
   }
  
   // ignore steps without energy deposit
   edepositEM = edepositHAD = 0.f;
   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.setID(unitID, time, primaryID, depth);
   } else {
     if (aStep->GetTotalEnergyDeposit() > 0.0 && (!ignoreReject)) {
       const G4TouchableHistory* touch = static_cast<const G4TouchableHistory*>(theTrack->GetTouchable());
       edm::LogVerbatim("CaloSim") << "CaloSD::ProcessHits: unitID= " << unitID << " currUnit=   " << currentID.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;
   }
  
   if (aStep->GetTotalEnergyDeposit() == 0.0) {
     return false;
   }
  
   double energy = getEnergyDeposit(aStep);
   if (energy > 0.0) {
     if (G4TrackToParticleID::isGammaElectronPositron(theTrack)) {
       edepositEM = energy;
     } else {
       edepositHAD = energy;
     }
 #ifdef EDM_ML_DEBUG
     G4TouchableHistory* touch = (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 << ") previousID= (" << previousID << ")";
 #endif
     if (!hitExists(aStep)) {
       currentHit = createNewHit(aStep, aStep->GetTrack());
     } else {
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("DoFineCalo") << "Not creating new hit, only updating currentHit " << currentHit->getUnitID();
 #endif
     }
     return true;
   }
   return false;
 }

References createNewHit(), currentHit, currentID, TauDecayModes::dec, LEDCalibrationChannels::depth, edepositEM, edepositHAD, HCALHighEnergyHPDFilter_cfi::energy, getDepth(), getEnergyDeposit(), getFromLibrary(), getTrackID(), CaloG4Hit::getUnitID(), hitExists(), ignoreReject, G4TrackToParticleID::isGammaElectronPositron(), isParameterized, SensitiveDetector::NaNTrap(), previousID, setDetUnitId(), CaloHitID::setID(), protons_cff::time, and CaloHitID::unitID().

◆ reset()

void CaloSD::reset ( void )

overridevirtual

Reimplemented from SensitiveCaloDetector.

Definition at line 809 of file CaloSD.cc.

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

References fpCaloG4HitAllocator.

◆ resetForNewPrimary()

void CaloSD::resetForNewPrimary ( const G4Step * aStep )

protected

Definition at line 670 of file CaloSD.cc.

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

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

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

◆ saveHit()

bool CaloSD::saveHit ( CaloG4Hit * aHit )

private

Definition at line 895 of file CaloSD.cc.

                                     {
   int tkID;
   bool ok = true;
  
   double time = aHit->getTimeSlice();
   if (corrTOFBeam)
     time += correctT;
  
   // Do track bookkeeping a little differently for fine tracking
   if (doFineCalo_ && aHit->getID().hasFineTrackID()) {
     tkID = aHit->getID().fineTrackID();
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("DoFineCalo") << "Saving hit " << aHit->getUnitID() << " with trackID=" << tkID;
 #endif
     // Check if the track is actually in the trackManager
     if (m_trackManager) {
       if (!m_trackManager->trackExists(tkID)) {
         ok = false;
         throw cms::Exception("Unknown", "CaloSD")
             << "aHit " << aHit->getUnitID() << " has fine trackID " << tkID << ", which is NOT IN THE TRACK MANAGER";
       }
     } else {
       ok = false;
       throw cms::Exception("Unknown", "CaloSD") << "m_trackManager not set, saveHit ok=false!";
     }
     // Take the aHit-information and move it to the actual PCaloHitContainer
     slave.get()->processHits(
         aHit->getUnitID(), aHit->getEM() / CLHEP::GeV, aHit->getHadr() / CLHEP::GeV, time, tkID, 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 " << aHit->getUnitID() << " with trackID=" << tkID
                                    << " (no fineTrackID)";
 #endif
     slave.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;
 }

References correctT, corrTOFBeam, TauDecayModes::dec, doFineCalo_, Exception, CaloHitID::fineTrackID(), CaloG4Hit::getDepth(), CaloG4Hit::getEM(), CaloG4Hit::getHadr(), CaloG4Hit::getID(), CaloG4Hit::getTimeSlice(), CaloG4Hit::getTrackID(), CaloG4Hit::getUnitID(), GeV, SimTrackManager::giveMotherNeeded(), CaloHitID::hasFineTrackID(), m_trackManager, convertSQLiteXML::ok, slave, protons_cff::time, and SimTrackManager::trackExists().

Referenced by update().

◆ setDetUnitId()

uint32_t CaloSD::setDetUnitId ( const G4Step * step )

overridepure virtual

Implements SensitiveDetector.

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

Referenced by ProcessHits().

◆ setNumberCheckedHits()

void CaloSD::setNumberCheckedHits ( int val )

inlineprotected

Definition at line 119 of file CaloSD.h.

119 { nCheckedHits = val; }

References nCheckedHits, and heppy_batch::val.

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

◆ setParameterized()

void CaloSD::setParameterized ( bool val )

inlineprotected

Definition at line 107 of file CaloSD.h.

107 { isParameterized = val; }

References isParameterized, and heppy_batch::val.

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

◆ setToGlobal()

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

protected

Definition at line 389 of file CaloSD.cc.

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

References DTRecHitClients_cfi::local.

◆ setToLocal()

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

protected

Definition at line 385 of file CaloSD.cc.

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

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

◆ setTrackID()

int CaloSD::setTrackID ( const G4Step * aStep )

protectedvirtual

Definition at line 839 of file CaloSD.cc.

                                           {
   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.trackID()) {
     resetForNewPrimary(aStep);
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CaloSim") << "CaloSD::setTrackID for " << GetName()
                               << " trackID= " << aStep->GetTrack()->GetTrackID() << " primaryID= " << primaryID;
 #endif
   return primaryID;
 }

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

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

◆ setUseMap()

void CaloSD::setUseMap ( bool val )

inlineprotected

Definition at line 108 of file CaloSD.h.

108 { useMap = val; }

References useMap, and heppy_batch::val.

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

◆ storeHit()

void CaloSD::storeHit ( CaloG4Hit * hit )

private

Definition at line 883 of file CaloSD.cc.

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

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

Referenced by createNewHit().

◆ update() [1/5]

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

overrideprotected

Definition at line 741 of file CaloSD.cc.

                                      {
   endEvent();
   slave.get()->ReserveMemory(theHC->entries());
  
   int count(0);
   int wrong(0);
   double eEM(0.0);
   double eHAD(0.0);
   double eEM2(0.0);
   double eHAD2(0.0);
   double tt(0.0);
   double zloc(0.0);
   double zglob(0.0);
   double ee(0.0);
   int hc_entries = theHC->entries();
   for (int i = 0; i < hc_entries; ++i) {
     if (!saveHit((*theHC)[i])) {
       ++wrong;
     }
     ++count;
     double x = (*theHC)[i]->getEM();
     eEM += x;
     eEM2 += x * x;
     x = (*theHC)[i]->getHadr();
     eHAD += x;
     eHAD2 += x * x;
     tt += (*theHC)[i]->getTimeSlice();
     ee += (*theHC)[i]->getIncidentEnergy();
     zglob += std::abs((*theHC)[i]->getEntry().z());
     zloc += std::abs((*theHC)[i]->getEntryLocal().z());
   }
  
   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);
   tt *= norm;
   ee *= norm;
   zglob *= norm;
   zloc *= norm;
  
 #ifdef EDM_ML_DEBUG
   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
   tkMap.erase(tkMap.begin(), tkMap.end());
   std::vector<std::unique_ptr<CaloG4Hit>>().swap(reusehit);
   if (useMap)
     hitMap.erase(hitMap.begin(), hitMap.end());
 }

References funct::abs(), submitPVResolutionJobs::count, endEvent(), edm::roottree::getEntry(), hitMap, mps_fire::i, 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().

◆ 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 707 of file CaloSD.cc.

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

References clearHits(), and initEvent().

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

◆ update() [3/5]

void CaloSD::update ( const BeginOfRun * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 705 of file CaloSD.cc.

705 { initRun(); }

References initRun().

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

◆ update() [4/5]

void CaloSD::update ( const BeginOfTrack * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack * >.

Definition at line 963 of file CaloSD.cc.

                                            {
   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
  
     if (theHC->entries() > 0)
       cleanHitCollection();
   }
 }

References cleanHitCollection(), SensitiveDetector::cmsTrackInformation(), TrackInformation::isPrimary(), primAncestor, and theHC.

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

◆ update() [5/5]

void CaloSD::update ( const EndOfTrack * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfTrack * >.

Definition at line 715 of file CaloSD.cc.

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

References SensitiveDetector::cmsTrackInformation(), TrackInformation::getIDonCaloSurface(), createfilelist::int, 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().

◆ updateHit()

void CaloSD::updateHit ( CaloG4Hit * aHit )

protected

Definition at line 659 of file CaloSD.cc.

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

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

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

Member Data Documentation

◆ cleanIndex

int CaloSD::cleanIndex

private

Definition at line 173 of file CaloSD.h.

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

◆ correctT

double CaloSD::correctT

private

Definition at line 180 of file CaloSD.h.

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

◆ corrTOFBeam

bool CaloSD::corrTOFBeam

private

Definition at line 169 of file CaloSD.h.

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

◆ currentHit

CaloG4Hit* CaloSD::currentHit

protected

Definition at line 143 of file CaloSD.h.

Referenced by checkHit(), ZdcSD::getFromLibrary(), hitExists(), processHit(), and ProcessHits().

◆ currentID

CaloHitID CaloSD::currentID

protected

Definition at line 139 of file CaloSD.h.

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

◆ doFineCalo_

bool CaloSD::doFineCalo_

private

Definition at line 181 of file CaloSD.h.

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

◆ edepositEM

float CaloSD::edepositEM

protected

Definition at line 137 of file CaloSD.h.

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

◆ edepositHAD

float CaloSD::edepositHAD

protected

Definition at line 137 of file CaloSD.h.

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

◆ eMinFine_

double CaloSD::eMinFine_

private

Definition at line 182 of file CaloSD.h.

Referenced by CaloSD().

◆ eminHit

double CaloSD::eminHit

protected

Definition at line 141 of file CaloSD.h.

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

◆ eminHitD

double CaloSD::eminHitD

private

Definition at line 179 of file CaloSD.h.

Referenced by CaloSD(), and filterHit().

◆ energyCut

double CaloSD::energyCut

protected

Definition at line 141 of file CaloSD.h.

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

◆ entranceLocal

G4ThreeVector CaloSD::entranceLocal

protected

Definition at line 134 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), ZdcSD::getFromLibrary(), ProcessHits(), and resetForNewPrimary().

◆ entrancePoint

G4ThreeVector CaloSD::entrancePoint

protected

Definition at line 133 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), ZdcSD::getFromLibrary(), ProcessHits(), and resetForNewPrimary().

◆ fineDetectors_

std::vector<Detector> CaloSD::fineDetectors_

private

Definition at line 187 of file CaloSD.h.

Referenced by CaloSD(), and isItFineCalo().

◆ forceSave

bool CaloSD::forceSave

protected

Definition at line 148 of file CaloSD.h.

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

◆ hcID

int CaloSD::hcID

private

Definition at line 171 of file CaloSD.h.

Referenced by Initialize().

◆ hitMap

std::map<CaloHitID, CaloG4Hit*> CaloSD::hitMap

private

Definition at line 184 of file CaloSD.h.

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

◆ ignoreReject

bool CaloSD::ignoreReject

private

Definition at line 167 of file CaloSD.h.

Referenced by ignoreRejection(), and ProcessHits().

◆ ignoreTrackID

bool CaloSD::ignoreTrackID

private

Definition at line 165 of file CaloSD.h.

Referenced by CaloSD().

◆ incidentEnergy

float CaloSD::incidentEnergy

protected

Definition at line 136 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), ZdcSD::getFromLibrary(), ProcessHits(), and resetForNewPrimary().

◆ isParameterized

bool CaloSD::isParameterized

private

Definition at line 166 of file CaloSD.h.

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

◆ kmaxIon

double CaloSD::kmaxIon

protected

Definition at line 146 of file CaloSD.h.

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

◆ kmaxNeutron

double CaloSD::kmaxNeutron

protected

Definition at line 146 of file CaloSD.h.

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

◆ kmaxProton

double CaloSD::kmaxProton

protected

Definition at line 146 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 158 of file CaloSD.h.

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

◆ meanResponse

std::unique_ptr<CaloMeanResponse> CaloSD::meanResponse

private

Definition at line 161 of file CaloSD.h.

Referenced by CaloSD(), and getResponseWt().

◆ nCheckedHits

int CaloSD::nCheckedHits

private

Definition at line 176 of file CaloSD.h.

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

◆ posGlobal

G4ThreeVector CaloSD::posGlobal

protected

Definition at line 135 of file CaloSD.h.

Referenced by CaloSD(), createNewHit(), ZdcSD::getFromLibrary(), hitExists(), and ProcessHits().

◆ previousID

CaloHitID CaloSD::previousID

protected

Definition at line 139 of file CaloSD.h.

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

◆ primAncestor

int CaloSD::primAncestor

private

Definition at line 172 of file CaloSD.h.

Referenced by CaloSD(), and update().

◆ primIDSaved

int CaloSD::primIDSaved

private

Definition at line 175 of file CaloSD.h.

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

◆ reusehit

std::vector<std::unique_ptr<CaloG4Hit> > CaloSD::reusehit

private

Definition at line 186 of file CaloSD.h.

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

◆ slave

std::unique_ptr<CaloSlaveSD> CaloSD::slave

private

Definition at line 160 of file CaloSD.h.

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

◆ suppressHeavy

bool CaloSD::suppressHeavy

protected

Definition at line 145 of file CaloSD.h.

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

◆ theHC

CaloG4HitCollection* CaloSD::theHC

private

Definition at line 163 of file CaloSD.h.

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

◆ timeSlice

float CaloSD::timeSlice

private

Definition at line 178 of file CaloSD.h.

Referenced by CaloSD().

◆ tkMap

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

private

Definition at line 185 of file CaloSD.h.

Referenced by createNewHit(), and update().

◆ tmaxHit

double CaloSD::tmaxHit

protected

Definition at line 141 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

private

Definition at line 174 of file CaloSD.h.

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

◆ useMap

bool CaloSD::useMap

private

Definition at line 168 of file CaloSD.h.

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

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

◆ getAttenuation()

◆ getDepth()

◆ getEnergyDeposit()

◆ getFromLibrary()

◆ getNumberOfHits()

◆ getResponseWt()

◆ getTrackID()

◆ hitBookkeepingFineCalo()

◆ hitExists()

◆ ignoreRejection()

◆ initEvent()

◆ Initialize()

◆ initRun()

◆ isItFineCalo()

◆ printableDecayChain()

◆ PrintAll()

◆ printDetectorLevels()

◆ processHit()

◆ ProcessHits() [1/2]

◆ ProcessHits() [2/2]

◆ reset()

◆ resetForNewPrimary()

◆ saveHit()

◆ setDetUnitId()

◆ setNumberCheckedHits()

◆ setParameterized()

◆ setToGlobal()

◆ setToLocal()

◆ setTrackID()

◆ setUseMap()

◆ storeHit()

◆ update() [1/5]

◆ update() [2/5]

◆ update() [3/5]

◆ update() [4/5]

◆ update() [5/5]

◆ updateHit()

Member Data Documentation

◆ cleanIndex

◆ correctT

◆ corrTOFBeam

◆ currentHit

◆ currentID

◆ doFineCalo_

◆ edepositEM

◆ edepositHAD

◆ eMinFine_

◆ eminHit

◆ eminHitD

◆ energyCut

◆ entranceLocal

◆ entrancePoint

◆ fineDetectors_

◆ forceSave

◆ hcID

◆ hitMap