#include <HCalSD.h>

Inheritance diagram for HCalSD:

Public Member Functions
	HCalSD (const std::string &, const HcalDDDSimConstants , const HcalDDDRecConstants , const HcalSimulationConstants , const HBHEDarkening , const HBHEDarkening , const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager )

uint32_t	setDetUnitId (const G4Step *step) override

void	setNumberingScheme (HcalNumberingScheme *)

	~HCalSD () override=default

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

void	clear () override

void	clearHits () override

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

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

void	Initialize (G4HCofThisEvent *HCE) override

bool	isItFineCalo (const G4VTouchable *touch)

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

void	PrintAll () override

G4bool	ProcessHits (G4Step step, G4TouchableHistory ) override

bool	ProcessHits (G4GFlashSpot aSpot, G4TouchableHistory ) override

void	reset () override

	~CaloSD () override

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

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

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

void	Initialize (G4HCofThisEvent *eventHC) override

bool	isCaloSD () const

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

	~SensitiveDetector () override

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

void	slotForUpdate (const BeginOfRun * iT)

virtual	~Observer ()

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

void	slotForUpdate (const BeginOfEvent * iT)

virtual	~Observer ()

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

void	slotForUpdate (const BeginOfTrack * iT)

virtual	~Observer ()

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

void	slotForUpdate (const EndOfTrack * iT)

virtual	~Observer ()

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

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

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

void	slotForUpdate (const BeginOfJob * iT)

virtual	~Observer ()

Protected Member Functions
void	endEvent () override

bool	filterHit (CaloG4Hit *, double) override

double	getEnergyDeposit (const G4Step *) override

bool	getFromLibrary (const G4Step *) override

void	initEvent (const BeginOfEvent *) override

void	initRun () override

void	update (const BeginOfRun *) override

void	update (const ::EndOfEvent *) override

void	update (const EndOfTrack *trk) override

void	update (const BeginOfTrack *trk) override

void	update (const BeginOfEvent *) override

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

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

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

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

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

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

virtual uint16_t	getDepth (const G4Step *)

int	getNumberOfHits (int k=0)

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

virtual int	getTrackID (const G4Track *)

bool	hitExists (const G4Step *, int k)

void	ignoreRejection ()

void	printDetectorLevels (const G4VTouchable *) const

void	processHit (const G4Step *step)

virtual void	processSecondHit (const G4Step , const G4Track )

void	resetForNewPrimary (const G4Step *)

void	setNumberCheckedHits (int val, int k=0)

void	setParameterized (bool val)

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

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

virtual int	setTrackID (const G4Step *)

void	setUseMap (bool val)

std::string	shortreprID (const CaloHitID &ID)

std::string	shortreprID (const CaloG4Hit *hit)

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

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

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

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

void	update (const ::EndOfEvent *) override

void	updateHit (CaloG4Hit *, int k)

Protected Member Functions inherited from SensitiveDetector
TrackInformation *	cmsTrackInformation (const G4Track *aTrack)

Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point) const

Local3DPoint	FinalStepPosition (const G4Step *step, coordinates) const

Local3DPoint	InitialStepPosition (const G4Step *step, coordinates) const

Local3DPoint	LocalPostStepPosition (const G4Step *step) const

Local3DPoint	LocalPreStepPosition (const G4Step *step) const

void	NaNTrap (const G4Step *step) const

void	setNames (const std::vector< std::string > &)

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

Private Member Functions
void	fillLogVolumeVector (const std::string &, const std::vector< std::string > &, std::vector< const G4LogicalVolume *> &)

void	getFromHFLibrary (const G4Step *step, bool &isKilled)

void	getFromParam (const G4Step *step, bool &isKilled)

void	getHitFibreBundle (const G4Step *step, bool type)

void	getHitPMT (const G4Step *step)

void	hitForFibre (const G4Step *step)

bool	isItConicalBundle (const G4LogicalVolume *)

bool	isItFibre (const G4LogicalVolume *)

bool	isItFibre (const G4String &)

bool	isItHF (const G4Step *)

bool	isItHF (const G4String &)

bool	isItinFidVolume (const G4ThreeVector &)

bool	isItPMT (const G4LogicalVolume *)

bool	isItScintillator (const G4Material *)

bool	isItStraightBundle (const G4LogicalVolume *)

double	layerWeight (int, const G4ThreeVector &, int, int)

void	modifyDepth (HcalNumberingFromDDD::HcalID &id)

void	plotHF (const G4ThreeVector &pos, bool emType)

void	plotProfile (const G4Step *step, const G4ThreeVector &pos, double edep, double time, int id)

void	printVolume (const G4VTouchable *touch) const

void	readWeightFromFile (const std::string &)

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

uint32_t	setDetUnitId (HcalNumberingFromDDD::HcalID &tmp)

Private Attributes
bool	agingFlagHB

bool	agingFlagHE

bool	applyFidCut

double	betaThr

double	birk1

double	birk2

double	birk3

double	deliveredLumi

int	depth_

std::vector< int >	detNull_

TH1F *	dist_ [9]

double	eminHitHB

double	eminHitHE

double	eminHitHF

double	eminHitHO

std::vector< const G4LogicalVolume * >	fibre1LV

std::vector< const G4LogicalVolume * >	fibre2LV

std::vector< const G4LogicalVolume * >	fibreLV

std::vector< std::string >	fibreNames

std::vector< double >	gpar

const HcalDDDSimConstants *	hcalConstants_

const HcalSimulationConstants *	hcalSimConstants_

std::vector< int >	hfLevels

std::vector< const G4LogicalVolume * >	hfLV

std::vector< std::string >	hfNames

std::unique_ptr< HFShower >	hfshower

TH1F *	hit_ [9]

TH1F *	hzvem

TH1F *	hzvhad

bool	isHF

std::map< uint32_t, double >	layerWeights

const HBHEDarkening *	m_HBDarkening

std::unique_ptr< HcalTestNS >	m_HcalTestNS

const HBHEDarkening *	m_HEDarkening

std::unique_ptr< HFDarkening >	m_HFDarkening

std::vector< const G4Material * >	materials

std::vector< std::string >	matNames

bool	neutralDensity

std::unique_ptr< HcalNumberingFromDDD >	numberingFromDDD

std::unique_ptr< HcalNumberingScheme >	numberingScheme

std::vector< const G4LogicalVolume * >	pmtLV

std::unique_ptr< HFShowerFibreBundle >	showerBundle

std::unique_ptr< HFShowerLibrary >	showerLibrary

std::unique_ptr< HFShowerParam >	showerParam

std::unique_ptr< HFShowerPMT >	showerPMT

bool	testNS_

bool	testNumber

TH1F *	time_ [9]

bool	useBirk

bool	useFibreBundle

bool	useHF

bool	useLayerWt

bool	useParam

bool	usePMTHit

bool	useShowerLibrary

double	weight_

Static Private Attributes
static constexpr double	maxRoff_ = 450.0

static constexpr double	maxZ_ = 10000.0

static constexpr double	minRoff_ = -1500.0

static constexpr double	slopeHE_ = 0.4

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

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

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

CaloG4Hit *	currentHit [2]

CaloHitID	currentID [2]

std::string	detName_

float	edepositEM

float	edepositHAD

double	eminHit

double	energyCut

G4ThreeVector	entranceLocal

G4ThreeVector	entrancePoint

bool	forceSave

std::vector< std::string >	hcn_

float	incidentEnergy

double	kmaxIon

double	kmaxNeutron

double	kmaxProton

int	nHC_

G4ThreeVector	posGlobal

CaloHitID	previousID [2]

bool	suppressHeavy

double	tmaxHit

std::vector< int >	useResMap_

Detailed Description

Definition at line 38 of file HCalSD.h.

Constructor & Destructor Documentation

◆ HCalSD()

HCalSD::HCalSD	(	const std::string &	name,
		const HcalDDDSimConstants *	hcns,
		const HcalDDDRecConstants *	hcnr,
		const HcalSimulationConstants *	hscs,
		const HBHEDarkening *	hbd,
		const HBHEDarkening *	hed,
		const SensitiveDetectorCatalog &	clg,
		edm::ParameterSet const &	p,
		const SimTrackManager *	manager
	)

Definition at line 43 of file HCalSD.cc.

References agingFlagHB, agingFlagHE, applyFidCut, betaThr, birk1, birk2, birk3, deliveredLumi, dist_, eminHitHB, eminHitHE, eminHitHF, eminHitHO, fibre1LV, fibre2LV, fibreLV, fibreNames, geometryDiff::file, fillLogVolumeVector(), g, relativeConstraints::geom, HcalDDDSimConstants::getGparHF(), HcalDDDSimConstants::getMaxDepth(), gpar, hcalConstants_, HcalSimulationParameters::hcalMaterialNames_, hcalSimConstants_, HcalSimulationConstants::hcalsimpar(), HcalSimulationParameters::hfFibreConicalNames_, HcalSimulationParameters::hfFibreNames_, HcalSimulationParameters::hfFibreStraightNames_, hfLevels, HcalSimulationParameters::hfLevels_, hfLV, hfNames, HcalSimulationParameters::hfNames_, HcalSimulationParameters::hfPMTNames_, hfshower, hit_, hzvem, hzvhad, mps_fire::i, CaloSD::kmaxIon, CaloSD::kmaxNeutron, CaloSD::kmaxProton, SummaryClient_cfi::labels, SensitiveDetectorCatalog::logicalNames(), materialBudgetVolume_cfi::lvNames, m_HcalTestNS, m_HFDarkening, TFileDirectory::make(), materials, matNames, Skims_PA_cff::name, neutralDensity, numberingFromDDD, numberingScheme, AlCaHLTBitMon_ParallelJobs::p, pmtLV, readWeightFromFile(), generator_cfi::scheme, setNumberingScheme(), CaloSD::setParameterized(), showerBundle, showerLibrary, showerParam, showerPMT, AlCaHLTBitMon_QueryRunRegistry::string, CaloSD::suppressHeavy, groupFilesInBlocks::temp, testNS_, testNumber, compare::tfile, time_, runGCPTkAlMap::title, useBirk, useFibreBundle, useHF, useLayerWt, useParam, usePMTHit, useShowerLibrary, and relativeConstraints::value.

     : CaloSD(name,
              clg,
              p,
              manager,
              (float)(p.getParameter<edm::ParameterSet>("HCalSD").getParameter<double>("TimeSliceUnit")),
              p.getParameter<edm::ParameterSet>("HCalSD").getParameter<bool>("IgnoreTrackID")),
       hcalConstants_(hcns),
       hcalSimConstants_(hscs),
       m_HBDarkening(hbd),
       m_HEDarkening(hed),
       isHF(false),
       weight_(1.0),
       depth_(1) {
   numberingFromDDD.reset(nullptr);
   numberingScheme.reset(nullptr);
   showerLibrary.reset(nullptr);
   hfshower.reset(nullptr);
   showerParam.reset(nullptr);
   showerPMT.reset(nullptr);
   showerBundle.reset(nullptr);
   m_HFDarkening.reset(nullptr);
   m_HcalTestNS.reset(nullptr);
 
   //static SimpleConfigurable<double> bk1(0.013, "HCalSD:BirkC1");
   //static SimpleConfigurable<double> bk2(0.0568,"HCalSD:BirkC2");
   //static SimpleConfigurable<double> bk3(1.75,  "HCalSD:BirkC3");
   // Values from NIM 80 (1970) 239-244: as implemented in Geant3
 
   bool dd4hep = p.getParameter<bool>("g4GeometryDD4hepSource");
   edm::ParameterSet m_HC = p.getParameter<edm::ParameterSet>("HCalSD");
   useBirk = m_HC.getParameter<bool>("UseBirkLaw");
   double bunit = (CLHEP::g / (CLHEP::MeV * CLHEP::cm2));
   birk1 = m_HC.getParameter<double>("BirkC1") * bunit;
   birk2 = m_HC.getParameter<double>("BirkC2");
   birk3 = m_HC.getParameter<double>("BirkC3");
   useShowerLibrary = m_HC.getParameter<bool>("UseShowerLibrary");
   useParam = m_HC.getParameter<bool>("UseParametrize");
   testNumber = m_HC.getParameter<bool>("TestNumberingScheme");
   neutralDensity = m_HC.getParameter<bool>("doNeutralDensityFilter");
   usePMTHit = m_HC.getParameter<bool>("UsePMTHits");
   betaThr = m_HC.getParameter<double>("BetaThreshold");
   eminHitHB = m_HC.getParameter<double>("EminHitHB") * MeV;
   eminHitHE = m_HC.getParameter<double>("EminHitHE") * MeV;
   eminHitHO = m_HC.getParameter<double>("EminHitHO") * MeV;
   eminHitHF = m_HC.getParameter<double>("EminHitHF") * MeV;
   useFibreBundle = m_HC.getParameter<bool>("UseFibreBundleHits");
   deliveredLumi = m_HC.getParameter<double>("DelivLuminosity");
   agingFlagHB = m_HC.getParameter<bool>("HBDarkening");
   agingFlagHE = m_HC.getParameter<bool>("HEDarkening");
   bool agingFlagHF = m_HC.getParameter<bool>("HFDarkening");
   useHF = m_HC.getUntrackedParameter<bool>("UseHF", true);
   bool forTBHC = m_HC.getUntrackedParameter<bool>("ForTBHCAL", false);
   bool forTBH2 = m_HC.getUntrackedParameter<bool>("ForTBH2", false);
   useLayerWt = m_HC.getUntrackedParameter<bool>("UseLayerWt", false);
   std::string file = m_HC.getUntrackedParameter<std::string>("WtFile", "None");
   testNS_ = m_HC.getUntrackedParameter<bool>("TestNS", false);
   edm::ParameterSet m_HF = p.getParameter<edm::ParameterSet>("HFShower");
   applyFidCut = m_HF.getParameter<bool>("ApplyFiducialCut");
   bool dumpGeom = m_HC.getUntrackedParameter<bool>("DumpGeometry", false);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "***************************************************"
                               << "\n"
                               << "* Constructing a HCalSD  with name " << name << "\n"
                               << "\n"
                               << "***************************************************";
 #endif
   edm::LogVerbatim("HcalSim") << "HCalSD:: Use of HF code is set to " << useHF
                               << "\n         Use of shower parametrization set to " << useParam
                               << "\n         Use of shower library is set to " << useShowerLibrary
                               << "\n         Use PMT Hit is set to " << usePMTHit << " with beta Threshold " << betaThr
                               << "\n         USe of FibreBundle Hit set to " << useFibreBundle
                               << "\n         Use of Birks law is set to " << useBirk
                               << " with three constants kB = " << birk1 / bunit << ", C1 = " << birk2
                               << ", C2 = " << birk3;
   edm::LogVerbatim("HcalSim") << "HCalSD:: Suppression Flag " << suppressHeavy << " protons below " << kmaxProton
                               << " MeV,"
                               << " neutrons below " << kmaxNeutron << " MeV and"
                               << " ions below " << kmaxIon << " MeV\n"
                               << "         Threshold for storing hits in HB: " << eminHitHB << " HE: " << eminHitHE
                               << " HO: " << eminHitHO << " HF: " << eminHitHF << "\n"
                               << "Delivered luminosity for Darkening " << deliveredLumi << " Flag (HE) " << agingFlagHE
                               << " Flag (HB) " << agingFlagHB << " Flag (HF) " << agingFlagHF << "\n"
                               << "Application of Fiducial Cut " << applyFidCut
                               << "Flag for test number|neutral density filter " << testNumber << " " << neutralDensity;
 
   if (forTBHC) {
     useHF = false;
     matNames.emplace_back("Scintillator");
   } else {
     matNames = hcalSimConstants_->hcalsimpar()->hcalMaterialNames_;
   }
 
   HcalNumberingScheme* scheme;
   if (testNumber || forTBH2) {
     scheme = dynamic_cast<HcalNumberingScheme*>(new HcalTestNumberingScheme(forTBH2));
   } else {
     scheme = new HcalNumberingScheme();
   }
   setNumberingScheme(scheme);
 
   // always call getFromLibrary() method to identify HF region
   setParameterized(true);
 
   const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
   //  std::vector<G4LogicalVolume *>::const_iterator lvcite;
   const G4LogicalVolume* lv;
   std::string attribute, value;
 
   if (useHF) {
     if (useParam) {
       showerParam = std::make_unique<HFShowerParam>(name, hcalConstants_, hcalSimConstants_->hcalsimpar(), p);
     } else {
       if (useShowerLibrary) {
         showerLibrary = std::make_unique<HFShowerLibrary>(name, hcalConstants_, hcalSimConstants_->hcalsimpar(), p);
       }
       hfshower = std::make_unique<HFShower>(name, hcalConstants_, hcalSimConstants_->hcalsimpar(), p, 0);
     }
 
     // HF volume names
     hfNames = hcalSimConstants_->hcalsimpar()->hfNames_;
     const std::vector<int>& temp = hcalSimConstants_->hcalsimpar()->hfLevels_;
 #ifdef EDM_ML_DEBUG
     std::stringstream ss0;
     ss0 << "HCalSD: Names to be tested for Volume = HF has " << hfNames.size() << " elements";
 #endif
     int addlevel = dd4hep ? 1 : 0;
     for (unsigned int i = 0; i < hfNames.size(); ++i) {
       G4String namv(static_cast<std::string>(dd4hep::dd::noNamespace(hfNames[i])));
       lv = nullptr;
       for (auto lvol : *lvs) {
         if (dd4hep::dd::noNamespace(lvol->GetName()) == namv) {
           lv = lvol;
           break;
         }
       }
       hfLV.emplace_back(lv);
       hfLevels.emplace_back(temp[i] + addlevel);
 #ifdef EDM_ML_DEBUG
       ss0 << "\n        HF[" << i << "] = " << namv << " LV " << lv << " at level " << (temp[i] + addlevel);
 #endif
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << ss0.str();
 #endif
     // HF Fibre volume names
     fibreNames = hcalSimConstants_->hcalsimpar()->hfFibreNames_;
     fillLogVolumeVector("HFFibre", fibreNames, fibreLV);
     const std::vector<std::string>& pmtNames = hcalSimConstants_->hcalsimpar()->hfPMTNames_;
     fillLogVolumeVector("HFPMT", pmtNames, pmtLV);
     const std::vector<std::string>& straightNames = hcalSimConstants_->hcalsimpar()->hfFibreStraightNames_;
     fillLogVolumeVector("HFFibreBundleStraight", straightNames, fibre1LV);
     const std::vector<std::string>& conicalNames = hcalSimConstants_->hcalsimpar()->hfFibreConicalNames_;
     fillLogVolumeVector("HFFibreBundleConical", conicalNames, fibre2LV);
   }
 
   //Material list for HB/HE/HO sensitive detectors
   const G4MaterialTable* matTab = G4Material::GetMaterialTable();
   std::vector<G4Material*>::const_iterator matite;
   for (auto const& namx : matNames) {
     const G4Material* mat = nullptr;
     for (matite = matTab->begin(); matite != matTab->end(); ++matite) {
       if (static_cast<std::string>(dd4hep::dd::noNamespace((*matite)->GetName())) == namx) {
         mat = (*matite);
         break;
       }
     }
     materials.emplace_back(mat);
   }
 #ifdef EDM_ML_DEBUG
   std::stringstream ss1;
   for (unsigned int i = 0; i < matNames.size(); ++i) {
     if (i / 10 * 10 == i) {
       ss1 << "\n";
     }
     ss1 << "  " << matNames[i];
   }
   edm::LogVerbatim("HcalSim") << "HCalSD: Material names for HCAL: " << ss1.str();
 #endif
   if (useLayerWt) {
     readWeightFromFile(file);
   }
   numberingFromDDD = std::make_unique<HcalNumberingFromDDD>(hcalConstants_);
 
   //Special Geometry parameters
   gpar = hcalConstants_->getGparHF();
 #ifdef EDM_ML_DEBUG
   std::stringstream ss2;
   for (unsigned int ig = 0; ig < gpar.size(); ig++) {
     ss2 << "\n         gpar[" << ig << "] = " << gpar[ig] / cm << " cm";
   }
   edm::LogVerbatim("HcalSim") << "Maximum depth for HF " << hcalConstants_->getMaxDepth(2) << gpar.size()
                               << " gpar (cm)" << ss2.str();
 #endif
 
   //Test Hcal Numbering Scheme
   if (testNS_)
     m_HcalTestNS = std::make_unique<HcalTestNS>(hcnr);
 
   for (int i = 0; i < 9; ++i) {
     hit_[i] = time_[i] = dist_[i] = nullptr;
   }
   hzvem = hzvhad = nullptr;
 
   if (agingFlagHF) {
     m_HFDarkening = std::make_unique<HFDarkening>(m_HC.getParameter<edm::ParameterSet>("HFDarkeningParameterBlock"));
   }
 #ifdef plotDebug
   edm::Service<TFileService> tfile;
 
   if (tfile.isAvailable()) {
     static const char* const labels[] = {"HB",
                                          "HE",
                                          "HO",
                                          "HF Absorber",
                                          "HF PMT",
                                          "HF Absorber Long",
                                          "HF Absorber Short",
                                          "HF PMT Long",
                                          "HF PMT Short"};
     TFileDirectory hcDir = tfile->mkdir("ProfileFromHCalSD");
     char name[20], title[60];
     for (int i = 0; i < 9; ++i) {
       sprintf(title, "Hit energy in %s", labels[i]);
       sprintf(name, "HCalSDHit%d", i);
       hit_[i] = hcDir.make<TH1F>(name, title, 2000, 0., 2000.);
       sprintf(title, "Energy (MeV)");
       hit_[i]->GetXaxis()->SetTitle(title);
       hit_[i]->GetYaxis()->SetTitle("Hits");
       sprintf(title, "Time of the hit in %s", labels[i]);
       sprintf(name, "HCalSDTime%d", i);
       time_[i] = hcDir.make<TH1F>(name, title, 2000, 0., 2000.);
       sprintf(title, "Time (ns)");
       time_[i]->GetXaxis()->SetTitle(title);
       time_[i]->GetYaxis()->SetTitle("Hits");
       sprintf(title, "Longitudinal profile in %s", labels[i]);
       sprintf(name, "HCalSDDist%d", i);
       dist_[i] = hcDir.make<TH1F>(name, title, 2000, 0., 2000.);
       sprintf(title, "Distance (mm)");
       dist_[i]->GetXaxis()->SetTitle(title);
       dist_[i]->GetYaxis()->SetTitle("Hits");
     }
     if (useHF && (!useParam)) {
       hzvem = hcDir.make<TH1F>("hzvem", "Longitudinal Profile (EM Part)", 330, 0.0, 1650.0);
       hzvem->GetXaxis()->SetTitle("Longitudinal Profile (EM Part)");
       hzvhad = hcDir.make<TH1F>("hzvhad", "Longitudinal Profile (Had Part)", 330, 0.0, 1650.0);
       hzvhad->GetXaxis()->SetTitle("Longitudinal Profile (Hadronic Part)");
     }
   }
 #endif
   if (dumpGeom) {
     std::unique_ptr<HcalNumberingFromDDD> hcn = std::make_unique<HcalNumberingFromDDD>(hcalConstants_);
     const auto& lvNames = clg.logicalNames(name);
     HcalDumpGeometry geom(lvNames, hcn.get(), testNumber, false);
     geom.update();
   }
 }

◆ ~HCalSD()

HCalSD::~HCalSD ( )

overridedefault

Member Function Documentation

◆ endEvent()

void HCalSD::endEvent ( )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 1115 of file HCalSD.cc.

References detNull_.

                       {
 #ifdef printDebug
   int sum = detNull_[0] + detNull_[1] + detNull_[2];
   if (sum > 0)
     edm::LogVerbatim("HcalSim") << "NullDets " << detNull_[0] << " " << detNull_[1] << " " << detNull_[2] << " "
                                 << detNull_[3] << " " << (static_cast<float>(sum) / (sum + detNull_[3]));
 #endif
 }

◆ fillLogVolumeVector()

void HCalSD::fillLogVolumeVector	(	const std::string &	value,
		const std::vector< std::string > &	lvnames,
		std::vector< const G4LogicalVolume *> &	lvvec
	)

private

Definition at line 310 of file HCalSD.cc.

References mps_fire::i.

Referenced by HCalSD().

                                                                            {
   const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
   const G4LogicalVolume* lv;
   std::stringstream ss3;
   ss3 << "HCalSD: " << lvnames.size() << " names to be tested for Volume <" << value << ">:";
   for (unsigned int i = 0; i < lvnames.size(); ++i) {
     G4String namv(static_cast<std::string>(dd4hep::dd::noNamespace(lvnames[i])));
     lv = nullptr;
     for (auto lvol : *lvs) {
       if (dd4hep::dd::noNamespace(lvol->GetName()) == namv) {
         lv = lvol;
         break;
       }
     }
     lvvec.emplace_back(lv);
     if (i / 10 * 10 == i) {
       ss3 << "\n";
     }
     ss3 << "  " << namv;
   }
   edm::LogVerbatim("HcalSim") << ss3.str();
 }

◆ filterHit()

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

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 580 of file HCalSD.cc.

References eminHitHB, eminHitHE, eminHitHF, eminHitHO, CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getUnitID(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, remoteMonitoring_LASER_era2018_cfg::threshold, hcalRecHitTable_cff::time, and CaloSD::tmaxHit.

                                                    {
   double threshold = 0;
   DetId theId(aHit->getUnitID());
   switch (theId.subdetId()) {
     case HcalBarrel:
       threshold = eminHitHB;
       break;
     case HcalEndcap:
       threshold = eminHitHE;
       break;
     case HcalOuter:
       threshold = eminHitHO;
       break;
     case HcalForward:
       threshold = eminHitHF;
       break;
     default:
       break;
   }
   return ((time <= tmaxHit) && (aHit->getEnergyDeposit() > threshold));
 }

◆ getEnergyDeposit()

double HCalSD::getEnergyDeposit ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 402 of file HCalSD.cc.

References agingFlagHB, agingFlagHE, birk1, birk2, birk3, HBHEDarkening::degradation(), deliveredLumi, hcalRecHitTable_cff::depth, depth_, CaloSD::getAttenuation(), getHitFibreBundle(), getHitPMT(), HcalDDDSimConstants::getLayer0Wt(), CaloSD::getResponseWt(), hcalConstants_, hitForFibre(), hcalRecHitTable_cff::ieta, HcalDetId::ietaAbs(), HcalDetId::iphi(), isHF, isItConicalBundle(), isItFibre(), isItPMT(), isItScintillator(), isItStraightBundle(), G4TrackToParticleID::isMuon(), TrackInformation::isPrimary(), CaloSD::kmaxIon, CaloSD::kmaxNeutron, CaloSD::kmaxProton, layerWeight(), m_HBDarkening, m_HEDarkening, neutralDensity, phi, setDetUnitId(), showerBundle, showerPMT, DetId::subdetId(), CaloSD::suppressHeavy, testNumber, HcalTestNumbering::unpackHcalIndex(), useBirk, useFibreBundle, useLayerWt, usePMTHit, useShowerLibrary, weight_, z, and HcalDetId::zside().

                                                    {
   double destep(0.0);
   auto const lv = aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume();
   auto const theTrack = aStep->GetTrack();
 
   if (isHF) {
     if (useShowerLibrary && G4TrackToParticleID::isMuon(theTrack) && isItFibre(lv)) {
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalSim") << "HCalSD: Hit at Fibre in LV " << lv->GetName() << " for track "
                                   << aStep->GetTrack()->GetTrackID() << " ("
                                   << aStep->GetTrack()->GetDefinition()->GetParticleName() << ")";
 #endif
       hitForFibre(aStep);
     }
     return destep;
   }
 
   if (isItPMT(lv)) {
     if (usePMTHit && showerPMT) {
       getHitPMT(aStep);
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << "HCalSD: Hit from PMT parametrization in LV " << lv->GetName() << " for Track "
                                 << aStep->GetTrack()->GetTrackID() << " ("
                                 << aStep->GetTrack()->GetDefinition()->GetParticleName() << ")";
 #endif
     return destep;
 
   } else if (isItStraightBundle(lv)) {
     if (useFibreBundle && showerBundle) {
       getHitFibreBundle(aStep, false);
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << "HCalSD: Hit from straight FibreBundle in LV: " << lv->GetName() << " for track "
                                 << aStep->GetTrack()->GetTrackID() << " ("
                                 << aStep->GetTrack()->GetDefinition()->GetParticleName() << ")";
 #endif
     return destep;
 
   } else if (isItConicalBundle(lv)) {
     if (useFibreBundle && showerBundle) {
       getHitFibreBundle(aStep, true);
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << "HCalSD: Hit from conical FibreBundle PV: " << lv->GetName() << " for track "
                                 << aStep->GetTrack()->GetTrackID() << " ("
                                 << aStep->GetTrack()->GetDefinition()->GetParticleName() << ")";
 #endif
     return destep;
   }
 
   // normal hit
   destep = aStep->GetTotalEnergyDeposit();
 
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   uint32_t detid = setDetUnitId(aStep);
   int det(0), ieta(0), phi(0), z(0), lay, depth(-1);
   if (testNumber) {
     HcalTestNumbering::unpackHcalIndex(detid, det, z, depth, ieta, phi, lay);
     if (z == 0) {
       z = -1;
     }
   } else {
     HcalDetId hcid(detid);
     det = hcid.subdetId();
     ieta = hcid.ietaAbs();
     phi = hcid.iphi();
     z = hcid.zside();
   }
   lay = (touch->GetReplicaNumber(0) / 10) % 100 + 1;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "HCalSD: det: " << det << " ieta: " << ieta << " iphi: " << phi << " zside " << z
                               << "  lay: " << lay - 2;
 #endif
   if (depth_ == 0 && (det == 1 || det == 2) && ((!testNumber) || neutralDensity))
     weight_ = hcalConstants_->getLayer0Wt(det, phi, z);
   if (useLayerWt) {
     G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
     weight_ = layerWeight(det + 2, hitPoint, depth_, lay);
   }
 
   if (agingFlagHB && m_HBDarkening && det == 1) {
     double dweight = m_HBDarkening->degradation(deliveredLumi, ieta, lay);
     weight_ *= dweight;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << "HCalSD: HB Lumi: " << deliveredLumi << " coefficient = " << dweight
                                 << " Weight= " << weight_;
 #endif
   }
 
   if (agingFlagHE && m_HEDarkening && det == 2) {
     double dweight = m_HEDarkening->degradation(deliveredLumi, ieta, lay);
     weight_ *= dweight;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << "HCalSD: HB Lumi: " << deliveredLumi << " coefficient = " << dweight
                                 << " Weight= " << weight_;
 #endif
   }
 
   if (suppressHeavy) {
     TrackInformation* trkInfo = (TrackInformation*)(theTrack->GetUserInformation());
     if (trkInfo) {
       int pdg = theTrack->GetDefinition()->GetPDGEncoding();
       if (!(trkInfo->isPrimary())) {  // Only secondary particles
         double ke = theTrack->GetKineticEnergy();
         if (pdg / 1000000000 == 1 && (pdg / 10000) % 100 > 0 && (pdg / 10) % 100 > 0 && ke < kmaxIon)
           weight_ = 0;
         if ((pdg == 2212) && (ke < kmaxProton))
           weight_ = 0;
         if ((pdg == 2112) && (ke < kmaxNeutron))
           weight_ = 0;
       }
     }
   }
   double wt0(1.0);
   if (useBirk) {
     const G4Material* mat = aStep->GetPreStepPoint()->GetMaterial();
     if (isItScintillator(mat))
       wt0 = getAttenuation(aStep, birk1, birk2, birk3);
   }
   weight_ *= wt0;
   double wt1 = getResponseWt(theTrack);
   double wt2 = theTrack->GetWeight();
   double edep = weight_ * wt1 * destep;
   if (wt2 > 0.0) {
     edep *= wt2;
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "HCalSD: edep= " << edep << " Det: " << det + 2 << " depth= " << depth_
                               << " weight= " << weight_ << " wt0= " << wt0 << " wt1= " << wt1 << " wt2= " << wt2;
 #endif
   return edep;
 }

◆ getFromHFLibrary()

void HCalSD::getFromHFLibrary	(	const G4Step *	step,
		bool &	isKilled
	)

private

Definition at line 742 of file HCalSD.cc.

References CaloSD::currentID, hcalRecHitTable_cff::depth, CaloSD::edepositEM, CaloSD::edepositHAD, hfClusterShapes_cfi::hits, mps_fire::i, G4TrackToParticleID::isGammaElectronPositron(), isItinFidVolume(), plotHF(), plotProfile(), CaloSD::processHit(), CaloSD::resetForNewPrimary(), setDetUnitId(), CaloHitID::setID(), CaloSD::setTrackID(), showerLibrary, hcalRecHitTable_cff::time, and weight_.

Referenced by getFromLibrary().

                                                                  {
   std::vector<HFShowerLibrary::Hit> hits = showerLibrary->getHits(aStep, isKilled, weight_, false);
   if (!isKilled || hits.empty()) {
     return;
   }
 
   int primaryID = setTrackID(aStep);
 
   // Reset entry point for new primary
   resetForNewPrimary(aStep);
 
   auto const theTrack = aStep->GetTrack();
   int det = 5;
 
   if (G4TrackToParticleID::isGammaElectronPositron(theTrack)) {
     edepositEM = 1. * GeV;
     edepositHAD = 0.;
   } else {
     edepositEM = 0.;
     edepositHAD = 1. * GeV;
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "HCalSD::getFromLibrary " << hits.size() << " hits for " << GetName() << " of "
                               << primaryID << " with " << theTrack->GetDefinition()->GetParticleName() << " of "
                               << aStep->GetPreStepPoint()->GetKineticEnergy() / GeV << " GeV";
 #endif
   for (unsigned int i = 0; i < hits.size(); ++i) {
     G4ThreeVector hitPoint = hits[i].position;
     if (isItinFidVolume(hitPoint)) {
       int depth = hits[i].depth;
       double time = hits[i].time;
       unsigned int unitID = setDetUnitId(det, hitPoint, depth);
       currentID[0].setID(unitID, time, primaryID, 0);
 #ifdef plotDebug
       plotProfile(aStep, hitPoint, 1.0 * GeV, time, depth);
       bool emType = G4TrackToParticleID::isGammaElectronPositron(theTrack->GetDefinition()->GetPDGEncoding());
       plotHF(hitPoint, emType);
 #endif
       processHit(aStep);
     }
   }
 }

◆ getFromLibrary()

bool HCalSD::getFromLibrary ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 335 of file HCalSD.cc.

References funct::abs(), deliveredLumi, depth_, getFromHFLibrary(), getFromParam(), CaloSD::getNumberOfHits(), mps_fire::i, createfilelist::int, G4TrackToParticleID::isGammaElectronPositron(), isHF, isItHF(), G4TrackToParticleID::isMuon(), G4TrackToParticleID::isStableHadronIon(), HFDarkening::lowZLimit, m_HFDarkening, HFDarkening::numberOfZLayers, alignCSCRings::r, HLT_2023v12_cff::track, HFDarkening::upperZLimit, useParam, useShowerLibrary, weight_, and z.

                                                {
   auto const track = aStep->GetTrack();
   depth_ = (aStep->GetPreStepPoint()->GetTouchable()->GetReplicaNumber(0)) % 10;
   weight_ = 1.0;
   bool kill(false);
   isHF = isItHF(aStep);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "GetFromLibrary: "
                               << (aStep->GetPreStepPoint()->GetTouchable()->GetVolume(0)->GetName())
 
                               << " isHF " << isHF << " darken " << (m_HFDarkening != nullptr) << " useParam "
                               << useParam << " useShowerLibrary " << useShowerLibrary << " Muon? "
                               << G4TrackToParticleID::isMuon(track) << " electron? "
                               << G4TrackToParticleID::isGammaElectronPositron(track) << " Stable Hadron? "
                               << G4TrackToParticleID::isStableHadronIon(track);
 #endif
   if (isHF) {
     if (m_HFDarkening) {
       G4ThreeVector hitPoint = aStep->GetPreStepPoint()->GetPosition();
       const double invcm = 1. / CLHEP::cm;
       double r = hitPoint.perp() * invcm;
       double z = std::abs(hitPoint.z()) * invcm;
       double dose_acquired = 0.;
       if (z >= HFDarkening::lowZLimit && z <= HFDarkening::upperZLimit) {
         unsigned int hfZLayer = (unsigned int)((z - HFDarkening::lowZLimit) / 5);
         if (hfZLayer >= HFDarkening::upperZLimit)
           hfZLayer = (HFDarkening::upperZLimit - 1);
         float normalized_lumi = m_HFDarkening->int_lumi(deliveredLumi);
         for (int i = hfZLayer; i != HFDarkening::numberOfZLayers; ++i) {
           dose_acquired = m_HFDarkening->dose(i, r);
           weight_ *= m_HFDarkening->degradation(normalized_lumi * dose_acquired);
         }
       }
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalSim") << "HCalSD::getFromLibrary: HFLumiDarkening at "
                                   << "r= " << r << ", z= " << z << " Dose= " << dose_acquired << " weight= " << weight_;
 #endif
     }
 
     if (useParam) {
       getFromParam(aStep, kill);
 #ifdef EDM_ML_DEBUG
       G4String nameVolume = aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName();
       edm::LogVerbatim("HcalSim") << "HCalSD: " << getNumberOfHits() << " hits from parametrization in " << nameVolume
                                   << " for Track " << track->GetTrackID() << " ("
                                   << track->GetDefinition()->GetParticleName() << ")";
 #endif
     } else if (useShowerLibrary && !G4TrackToParticleID::isMuon(track)) {
       if (G4TrackToParticleID::isGammaElectronPositron(track) || G4TrackToParticleID::isStableHadronIon(track)) {
 #ifdef EDM_ML_DEBUG
         auto nameVolume = aStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetName();
         edm::LogVerbatim("HcalSim") << "HCalSD: Starts shower library from " << nameVolume << " for Track "
                                     << track->GetTrackID() << " (" << track->GetDefinition()->GetParticleName() << ")";
 
 #endif
         getFromHFLibrary(aStep, kill);
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "HCalSD::getFromLibrary ID= " << track->GetTrackID() << " ("
                               << track->GetDefinition()->GetParticleName() << ") kill= " << kill
                               << " weight= " << weight_ << " depth= " << depth_ << " isHF: " << isHF;
 #endif
   return kill;
 }

◆ getFromParam()

void HCalSD::getFromParam	(	const G4Step *	step,
		bool &	isKilled
	)

private

Definition at line 827 of file HCalSD.cc.

References CaloSD::currentID, hcalRecHitTable_cff::depth, CaloSD::edepositEM, CaloSD::edepositHAD, hfClusterShapes_cfi::hits, mps_fire::i, plotProfile(), CaloSD::processHit(), setDetUnitId(), CaloHitID::setID(), CaloSD::setTrackID(), showerParam, hcalRecHitTable_cff::time, and weight_.

Referenced by getFromLibrary().

                                                              {
   std::vector<HFShowerParam::Hit> hits = showerParam->getHits(aStep, weight_, isKilled);
   if (!isKilled || hits.empty()) {
     return;
   }
 
   int primaryID = setTrackID(aStep);
   int det = 5;
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "HCalSD::getFromParam " << hits.size() << " hits for " << GetName() << " of "
                               << primaryID << " with " << aStep->GetTrack()->GetDefinition()->GetParticleName()
                               << " of " << aStep->GetPreStepPoint()->GetKineticEnergy() / GeV
                               << " GeV in detector type " << det;
 #endif
   for (unsigned int i = 0; i < hits.size(); ++i) {
     G4ThreeVector hitPoint = hits[i].position;
     int depth = hits[i].depth;
     double time = hits[i].time;
     unsigned int unitID = setDetUnitId(det, hitPoint, depth);
     currentID[0].setID(unitID, time, primaryID, 0);
     edepositEM = hits[i].edep * GeV;
     edepositHAD = 0.;
 #ifdef plotDebug
     plotProfile(aStep, hitPoint, edepositEM, time, depth);
 #endif
     processHit(aStep);
   }
 }

◆ getHitFibreBundle()

void HCalSD::getHitFibreBundle	(	const G4Step *	step,
		bool	type
	)

private

Definition at line 916 of file HCalSD.cc.

References HLT_2023v12_cff::beta, CaloSD::currentID, TauDecayModes::dec, hcalRecHitTable_cff::depth, CaloSD::edepositEM, CaloSD::edepositHAD, CaloSD::energyCut, HcalForward, numberingFromDDD, phi, plotProfile(), CaloSD::processHit(), CaloSD::resetForNewPrimary(), findQualityFiles::rr, setDetUnitId(), CaloHitID::setID(), showerBundle, hcalRecHitTable_cff::time, and createJobs::tmp.

Referenced by getEnergyDeposit().

                                                              {
   auto const preStepPoint = aStep->GetPreStepPoint();
   auto const theTrack = aStep->GetTrack();
   double edep = showerBundle->getHits(aStep, type);
 
   if (edep >= 0.0) {
     edep *= GeV;
     double etrack = preStepPoint->GetKineticEnergy();
     int primaryID = 0;
     if (etrack >= energyCut) {
       primaryID = theTrack->GetTrackID();
     } else {
       primaryID = theTrack->GetParentID();
       if (primaryID == 0)
         primaryID = theTrack->GetTrackID();
     }
     // Reset entry point for new primary
     resetForNewPrimary(aStep);
     //
     int det = static_cast<int>(HcalForward);
     const G4ThreeVector& hitPoint = preStepPoint->GetPosition();
     double rr = hitPoint.x() * hitPoint.x() + hitPoint.y() * hitPoint.y();
     double phi = rr == 0. ? 0. : atan2(hitPoint.y(), hitPoint.x());
     double etaR = showerBundle->getRadius();
     int depth = 1;
     if (etaR < 0.) {
       depth = 2;
       etaR = -etaR;
     }
     if (hitPoint.z() < 0.)
       etaR = -etaR;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << "HCalSD::Hit for Detector " << det << " etaR " << etaR << " phi " << phi / deg
                                 << " depth " << depth;
 #endif
     double time = (aStep->GetPostStepPoint()->GetGlobalTime());
     uint32_t unitID = 0;
     if (numberingFromDDD) {
       HcalNumberingFromDDD::HcalID tmp = numberingFromDDD->unitID(det, etaR, phi, depth, 1);
       unitID = setDetUnitId(tmp);
     }
     if (type)
       currentID[0].setID(unitID, time, primaryID, 3);
     else
       currentID[0].setID(unitID, time, primaryID, 2);
 
     edepositHAD = aStep->GetTotalEnergyDeposit();
     edepositEM = -edepositHAD + edep;
 #ifdef plotDebug
     plotProfile(aStep, hitPoint, edep, time, depth);
 #endif
 #ifdef EDM_ML_DEBUG
     double beta = preStepPoint->GetBeta();
     edm::LogVerbatim("HcalSim") << "HCalSD::getHitFibreBundle 1 hit for " << GetName() << " of " << primaryID
                                 << " with " << theTrack->GetDefinition()->GetParticleName() << " of "
                                 << preStepPoint->GetKineticEnergy() / GeV << " GeV with velocity " << beta << " UnitID "
                                 << std::hex << unitID << std::dec;
 #endif
     processHit(aStep);
   }  // non-zero energy deposit
 }

◆ getHitPMT()

void HCalSD::getHitPMT ( const G4Step * step )

private

Definition at line 857 of file HCalSD.cc.

References HLT_2023v12_cff::beta, CaloSD::currentID, TauDecayModes::dec, hcalRecHitTable_cff::depth, CaloSD::edepositEM, CaloSD::edepositHAD, CaloSD::energyCut, HcalForward, numberingFromDDD, phi, plotProfile(), CaloSD::processHit(), CaloSD::resetForNewPrimary(), findQualityFiles::rr, setDetUnitId(), CaloHitID::setID(), showerPMT, hcalRecHitTable_cff::time, and createJobs::tmp.

Referenced by getEnergyDeposit().

                                           {
   auto const preStepPoint = aStep->GetPreStepPoint();
   auto const theTrack = aStep->GetTrack();
   double edep = showerPMT->getHits(aStep);
 
   if (edep >= 0.) {
     edep *= GeV;
     double etrack = preStepPoint->GetKineticEnergy();
     int primaryID = 0;
     if (etrack >= energyCut) {
       primaryID = theTrack->GetTrackID();
     } else {
       primaryID = theTrack->GetParentID();
       if (primaryID == 0)
         primaryID = theTrack->GetTrackID();
     }
     // Reset entry point for new primary
     resetForNewPrimary(aStep);
     //
     int det = static_cast<int>(HcalForward);
     const G4ThreeVector& hitPoint = preStepPoint->GetPosition();
     double rr = (hitPoint.x() * hitPoint.x() + hitPoint.y() * hitPoint.y());
     double phi = (rr == 0. ? 0. : atan2(hitPoint.y(), hitPoint.x()));
     double etaR = showerPMT->getRadius();
     int depth = 1;
     if (etaR < 0) {
       depth = 2;
       etaR = -etaR;
     }
     if (hitPoint.z() < 0)
       etaR = -etaR;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << "HCalSD::Hit for Detector " << det << " etaR " << etaR << " phi " << phi / deg
                                 << " depth " << depth;
 #endif
     double time = (aStep->GetPostStepPoint()->GetGlobalTime());
     uint32_t unitID = 0;
     if (numberingFromDDD) {
       HcalNumberingFromDDD::HcalID tmp = numberingFromDDD->unitID(det, etaR, phi, depth, 1);
       unitID = setDetUnitId(tmp);
     }
     currentID[0].setID(unitID, time, primaryID, 1);
 
     edepositHAD = aStep->GetTotalEnergyDeposit();
     edepositEM = -edepositHAD + edep;
 #ifdef plotDebug
     plotProfile(aStep, hitPoint, edep, time, depth);
 #endif
 #ifdef EDM_ML_DEBUG
     double beta = preStepPoint->GetBeta();
     edm::LogVerbatim("HcalSim") << "HCalSD::getHitPMT 1 hit for " << GetName() << " of " << primaryID << " with "
                                 << theTrack->GetDefinition()->GetParticleName() << " of "
                                 << preStepPoint->GetKineticEnergy() / GeV << " GeV with velocity " << beta << " UnitID "
                                 << std::hex << unitID << std::dec;
 #endif
     processHit(aStep);
   }
 }

◆ hitForFibre()

void HCalSD::hitForFibre ( const G4Step * step )

private

Definition at line 785 of file HCalSD.cc.

References CaloSD::currentID, hcalRecHitTable_cff::depth, CaloSD::edepositEM, CaloSD::edepositHAD, hfshower, hfClusterShapes_cfi::hits, mps_fire::i, G4TrackToParticleID::isGammaElectronPositron(), isItinFidVolume(), plotHF(), plotProfile(), CaloSD::processHit(), setDetUnitId(), CaloHitID::setID(), CaloSD::setTrackID(), hcalRecHitTable_cff::time, funct::true, and weight_.

Referenced by getEnergyDeposit().

                                             {  // if not ParamShower
 
   std::vector<HFShower::Hit> hits = hfshower->getHits(aStep, weight_);
   if (hits.empty()) {
     return;
   }
 
   auto const theTrack = aStep->GetTrack();
   int primaryID = setTrackID(aStep);
   int det = 5;
 
   if (G4TrackToParticleID::isGammaElectronPositron(theTrack)) {
     edepositEM = 1. * GeV;
     edepositHAD = 0.;
   } else {
     edepositEM = 0.;
     edepositHAD = 1. * GeV;
   }
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "HCalSD::hitForFibre " << hits.size() << " hits for " << GetName() << " of "
                               << primaryID << " with " << theTrack->GetDefinition()->GetParticleName() << " of "
                               << aStep->GetPreStepPoint()->GetKineticEnergy() / GeV << " GeV in detector type " << det;
 #endif
 
   for (unsigned int i = 0; i < hits.size(); ++i) {
     G4ThreeVector hitPoint = hits[i].position;
     if (isItinFidVolume(hitPoint)) {
       int depth = hits[i].depth;
       double time = hits[i].time;
       unsigned int unitID = setDetUnitId(det, hitPoint, depth);
       currentID[0].setID(unitID, time, primaryID, 0);
 #ifdef plotDebug
       plotProfile(aStep, hitPoint, edepositEM, time, depth);
       bool emType = (edepositEM > 0.) ? true : false;
       plotHF(hitPoint, emType);
 #endif
       processHit(aStep);
     }
   }
 }

◆ initEvent()

void HCalSD::initEvent ( const BeginOfEvent * )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 1109 of file HCalSD.cc.

References detNull_.

                                           {
 #ifdef printDebug
   detNull_ = {0, 0, 0, 0};
 #endif
 }

◆ initRun()

void HCalSD::initRun ( )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 578 of file HCalSD.cc.

578 {}

◆ isItConicalBundle()

bool HCalSD::isItConicalBundle ( const G4LogicalVolume * lv )

private

Definition at line 710 of file HCalSD.cc.

References fibre2LV.

Referenced by getEnergyDeposit().

                                                         {
   for (auto lvol : fibre2LV)
     if (lv == lvol) {
       return true;
     }
   return false;
 }

◆ isItFibre() [1/2]

bool HCalSD::isItFibre ( const G4LogicalVolume * lv )

private

Definition at line 678 of file HCalSD.cc.

References fibreLV.

Referenced by getEnergyDeposit().

                                                 {
   for (auto lvol : fibreLV)
     if (lv == lvol) {
       return true;
     }
   return false;
 }

◆ isItFibre() [2/2]

bool HCalSD::isItFibre ( const G4String & name )

private

Definition at line 686 of file HCalSD.cc.

References fibreNames, and Skims_PA_cff::name.

                                            {
   for (const auto& nam : fibreNames)
     if (name == static_cast<G4String>(nam)) {
       return true;
     }
   return false;
 }

◆ isItHF() [1/2]

bool HCalSD::isItHF ( const G4Step * aStep )

private

Definition at line 657 of file HCalSD.cc.

References hfLevels, hfLV, hfNames, and GCP_tree_cfg::levels.

Referenced by getFromLibrary().

                                        {
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   int levels = (touch->GetHistoryDepth()) + 1;
   for (unsigned int it = 0; it < hfNames.size(); ++it) {
     if (levels >= hfLevels[it]) {
       const G4LogicalVolume* lv = touch->GetVolume(levels - hfLevels[it])->GetLogicalVolume();
       if (lv == hfLV[it])
         return true;
     }
   }
   return false;
 }

◆ isItHF() [2/2]

bool HCalSD::isItHF ( const G4String & name )

private

Definition at line 670 of file HCalSD.cc.

References hfNames, and Skims_PA_cff::name.

                                         {
   for (const auto& nam : hfNames)
     if (name == static_cast<G4String>(nam)) {
       return true;
     }
   return false;
 }

◆ isItinFidVolume()

bool HCalSD::isItinFidVolume ( const G4ThreeVector & hitPoint )

private

Definition at line 726 of file HCalSD.cc.

References applyFidCut, RemoveAddSevLevel::flag, and HFFibreFiducial::PMTNumber().

Referenced by getFromHFLibrary(), and hitForFibre().

                                                           {
   bool flag = true;
   if (applyFidCut) {
     int npmt = HFFibreFiducial::PMTNumber(hitPoint);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << "HCalSD::isItinFidVolume:#PMT= " << npmt << " for hit point " << hitPoint;
 #endif
     if (npmt <= 0)
       flag = false;
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "HCalSD::isItinFidVolume: point " << hitPoint << " return flag " << flag;
 #endif
   return flag;
 }

◆ isItPMT()

bool HCalSD::isItPMT ( const G4LogicalVolume * lv )

private

Definition at line 694 of file HCalSD.cc.

References pmtLV.

Referenced by getEnergyDeposit().

                                               {
   for (auto lvol : pmtLV)
     if (lv == lvol) {
       return true;
     }
   return false;
 }

◆ isItScintillator()

bool HCalSD::isItScintillator ( const G4Material * mat )

private

Definition at line 718 of file HCalSD.cc.

References materials.

Referenced by getEnergyDeposit().

                                                    {
   for (auto amat : materials)
     if (amat == mat) {
       return true;
     }
   return false;
 }

◆ isItStraightBundle()

bool HCalSD::isItStraightBundle ( const G4LogicalVolume * lv )

private

Definition at line 702 of file HCalSD.cc.

References fibre1LV.

Referenced by getEnergyDeposit().

                                                          {
   for (auto lvol : fibre1LV)
     if (lv == lvol) {
       return true;
     }
   return false;
 }

◆ layerWeight()

double HCalSD::layerWeight	(	int	det,
		const G4ThreeVector &	pos,
		int	depth,
		int	lay
	)

private

Definition at line 1002 of file HCalSD.cc.

References TauDecayModes::dec, hcalRecHitTable_cff::depth, EDM_ML_DEBUG, layerWeights, modifyDepth(), numberingFromDDD, HcalTestNumbering::packHcalIndex(), and createJobs::tmp.

Referenced by getEnergyDeposit().

                                                                                 {
   double wt = 1.;
   if (numberingFromDDD) {
     //get the ID's as eta, phi, depth, ... indices
     HcalNumberingFromDDD::HcalID tmp =
         numberingFromDDD->unitID(det, math::XYZVectorD(pos.x(), pos.y(), pos.z()), depth, lay);
     modifyDepth(tmp);
     uint32_t id = HcalTestNumbering::packHcalIndex(tmp.subdet, tmp.zside, 1, tmp.etaR, tmp.phis, tmp.lay);
     std::map<uint32_t, double>::const_iterator ite = layerWeights.find(id);
     if (ite != layerWeights.end())
       wt = ite->second;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << "HCalSD::layerWeight: ID " << std::hex << id << std::dec << " (" << tmp.subdet << "/"
                                 << tmp.zside << "/1/" << tmp.etaR << "/" << tmp.phis << "/" << tmp.lay << ") Weight "
                                 << wt;
 #endif
   }
   return wt;
 }

◆ modifyDepth()

void HCalSD::modifyDepth ( HcalNumberingFromDDD::HcalID & id )

private

Definition at line 1095 of file HCalSD.cc.

References hcalRecHitTable_cff::depth, depth_, hcalConstants_, hcalRecHitTable_cff::ieta, HcalDDDSimConstants::maxHFDepth(), and testNumber.

Referenced by layerWeight(), and setDetUnitId().

                                                        {
   if (id.subdet == 4) {
     int ieta = (id.zside == 0) ? -id.etaR : id.etaR;
     if (hcalConstants_->maxHFDepth(ieta, id.phis) > 2) {
       if (id.depth <= 2) {
         if (G4UniformRand() > 0.5)
           id.depth += 2;
       }
     }
   } else if ((id.subdet == 1 || id.subdet == 2) && testNumber) {
     id.depth = (depth_ == 0) ? 1 : 2;
   }
 }

◆ plotHF()

void HCalSD::plotHF	(	const G4ThreeVector &	pos,
		bool	emType
	)

private

Definition at line 1084 of file HCalSD.cc.

References funct::abs(), gpar, hzvem, hzvhad, and gpuVertexFinder::zv.

Referenced by getFromHFLibrary(), and hitForFibre().

                                                               {
   double zv = std::abs(hitPoint.z()) - gpar[4];
   if (emType) {
     if (hzvem != nullptr)
       hzvem->Fill(zv);
   } else {
     if (hzvhad != nullptr)
       hzvhad->Fill(zv);
   }
 }

◆ plotProfile()

void HCalSD::plotProfile	(	const G4Step *	step,
		const G4ThreeVector &	pos,
		double	edep,
		double	time,
		int	id
	)

private

Definition at line 1022 of file HCalSD.cc.

References funct::abs(), hcalRecHitTable_cff::depth, dist_, newFWLiteAna::found, hit_, heavyIonCSV_trainingSettings::idx, cuy::ii, DTRecHitClients_cfi::local, dqmiodumpmetadata::n, Skims_PA_cff::name, names, hcalRecHitTable_cff::time, and time_.

Referenced by getFromHFLibrary(), getFromParam(), getHitFibreBundle(), getHitPMT(), and hitForFibre().

                                                                                                            {
   const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
   static const unsigned int names = 10;
   static const G4String modName[names] = {
       "HEModule", "HVQF", "HBModule", "MBAT", "MBBT", "MBBTC", "MBBT_R1P", "MBBT_R1M", "MBBT_R1PX", "MBBT_R1MX"};
   G4ThreeVector local;
   bool found = false;
   double depth = -2000;
   int idx = 4;
   for (int n = 0; n < touch->GetHistoryDepth(); ++n) {
     G4String name(static_cast<std::string>(dd4hep::dd::noNamespace(touch->GetVolume(n)->GetName())));
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalSim") << "plotProfile Depth " << n << " Name " << name;
 #endif
     for (unsigned int ii = 0; ii < names; ++ii) {
       if (name == modName[ii]) {
         found = true;
         int dn = touch->GetHistoryDepth() - n;
         local = touch->GetHistory()->GetTransform(dn).TransformPoint(global);
         if (ii == 0) {
           depth = local.z() - 4006.5;
           idx = 1;
         } else if (ii == 1) {
           depth = local.z() + 825.0;
           idx = 3;
         } else if (ii == 2) {
           depth = local.x() - 1775.;
           idx = 0;
         } else {
           depth = local.y() + 15.;
           idx = 2;
         }
         break;
       }
     }
     if (found)
       break;
   }
   if (!found)
     depth = std::abs(global.z()) - 11500;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "plotProfile Found " << found << " Global " << global << " Local " << local
                               << " depth " << depth << " ID " << id << " EDEP " << edep << " Time " << time;
 #endif
   if (hit_[idx] != nullptr)
     hit_[idx]->Fill(edep);
   if (time_[idx] != nullptr)
     time_[idx]->Fill(time, edep);
   if (dist_[idx] != nullptr)
     dist_[idx]->Fill(depth, edep);
   int jd = 2 * idx + id - 7;
   if (jd >= 0 && jd < 4) {
     jd += 5;
     if (hit_[jd] != nullptr)
       hit_[jd]->Fill(edep);
     if (time_[jd] != nullptr)
       time_[jd]->Fill(time, edep);
     if (dist_[jd] != nullptr)
       dist_[jd]->Fill(depth, edep);
   }
 }

◆ printVolume()

void HCalSD::printVolume ( const G4VTouchable * touch ) const

private

Definition at line 1124 of file HCalSD.cc.

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

Referenced by setDetUnitId().

                                                         {
   if (touch) {
 #ifdef EDM_ML_DEBUG
     int level = ((touch->GetHistoryDepth()) + 1);
     edm::LogVerbatim("CaloSimX") << "HCalSD::printVolume with " << level << " levels";
     static const std::string unknown("Unknown");
     //Get name and copy numbers
     for (int ii = 0; ii < level; ii++) {
       int i = level - ii - 1;
       G4VPhysicalVolume* pv = touch->GetVolume(i);
       G4String name = (pv != nullptr) ? pv->GetName() : unknown;
       G4int copyno = touch->GetReplicaNumber(i);
       edm::LogVerbatim("HcalSim") << "[" << ii << "] " << name << ":" << copyno;
     }
 #endif
   }
 }

◆ readWeightFromFile()

void HCalSD::readWeightFromFile ( const std::string & fName )

private

Definition at line 978 of file HCalSD.cc.

References TauDecayModes::dec, mps_splice::entry, fileCollector2::fName, recoMuon::in, timingPdfMaker::infile, layerWeights, HcalTestNumbering::packHcalIndex(), phi, useLayerWt, and ecaldqm::zside().

Referenced by HCalSD().

                                                       {
   std::ifstream infile;
   int entry = 0;
   infile.open(fName.c_str(), std::ios::in);
   if (infile) {
     int det, zside, etaR, phi, lay;
     double wt;
     while (infile >> det >> zside >> etaR >> phi >> lay >> wt) {
       uint32_t id = HcalTestNumbering::packHcalIndex(det, zside, 1, etaR, phi, lay);
       layerWeights.insert(std::pair<uint32_t, double>(id, wt));
       ++entry;
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalSim") << "HCalSD::readWeightFromFile:Entry " << entry << " ID " << std::hex << id
                                   << std::dec << " (" << det << "/" << zside << "/1/" << etaR << "/" << phi << "/"
                                   << lay << ") Weight " << wt;
 #endif
     }
     infile.close();
   }
   edm::LogVerbatim("HcalSim") << "HCalSD::readWeightFromFile: reads " << entry << " weights from " << fName;
   if (entry <= 0)
     useLayerWt = false;
 }

◆ setDetUnitId() [1/3]

uint32_t HCalSD::setDetUnitId ( const G4Step * step )

overridevirtual

Implements CaloSD.

Definition at line 536 of file HCalSD.cc.

References TauDecayModes::dec, hcalRecHitTable_cff::depth, HcalDetId::depth(), heavyIonCSV_trainingSettings::idx, HcalDetId::ietaAbs(), HcalDetId::iphi(), printVolume(), DetId::subdetId(), testNumber, HcalTestNumbering::unpackHcalIndex(), HLTMuonOfflineAnalyzer_cfi::z0, and HcalDetId::zside().

Referenced by getEnergyDeposit(), getFromHFLibrary(), getFromParam(), getHitFibreBundle(), getHitPMT(), hitForFibre(), and setDetUnitId().

                                                  {
   auto const prePoint = aStep->GetPreStepPoint();
   auto const touch = prePoint->GetTouchable();
   const G4ThreeVector& hitPoint = prePoint->GetPosition();
 
   int depth = (touch->GetReplicaNumber(0)) % 10 + 1;
   int lay = (touch->GetReplicaNumber(0) / 10) % 100 + 1;
   int det = (touch->GetReplicaNumber(1)) / 1000;
 
   uint32_t idx = setDetUnitId(det, hitPoint, depth, lay);
 #ifdef EDM_ML_DEBUG
   if (depth == 1) {
     edm::LogVerbatim("HcalSim") << "HCalSD: Check for " << det << ":" << depth << ":" << lay << " ID " << std::hex
                                 << idx << std::dec;
     int det0, z0, depth0, eta0, phi0, lay0(-1);
     if (testNumber) {
       HcalTestNumbering::unpackHcalIndex(idx, det0, z0, depth0, eta0, phi0, lay0);
     } else {
       HcalDetId hcid0(idx);
       det0 = hcid0.subdetId();
       eta0 = hcid0.ietaAbs();
       phi0 = hcid0.iphi();
       z0 = hcid0.zside();
       depth0 = hcid0.depth();
     }
     edm::LogVerbatim("HcalSim") << "HCalSD: det|z|depth|eta|phi|lay " << det0 << ":" << z0 << ":" << depth0 << ":"
                                 << eta0 << ":" << phi0 << ":" << lay0;
     printVolume(touch);
   }
 #endif
   return idx;
 }

◆ setDetUnitId() [2/3]

uint32_t HCalSD::setDetUnitId	(	int	det,
		const G4ThreeVector &	pos,
		int	depth,
		int	lay = `1`
	)

private

Definition at line 602 of file HCalSD.cc.

References funct::abs(), hcalRecHitTable_cff::depth, detNull_, PVValHelper::eta, hcalConstants_, l1ctLayer2EG_cff::id, HcalDDDSimConstants::isHE(), maxRoff_, maxZ_, minRoff_, numberingFromDDD, setDetUnitId(), slopeHE_, and createJobs::tmp.

                                                                                        {
   uint32_t id = 0;
   if (det == 0) {
 #ifdef printDebug
     double eta = std::abs(pos.eta());
 #endif
     if (std::abs(pos.z()) > maxZ_) {
       det = 5;
 #ifdef printDebug
       if (eta < 2.868)
         ++detNull_[2];
 #endif
     } else if (!(hcalConstants_->isHE())) {
       det = 3;
 #ifdef printDebug
       ++detNull_[0];
 #endif
     } else {
       double minR = minRoff_ + slopeHE_ * std::abs(pos.z());
       double maxR = maxRoff_ + slopeHE_ * std::abs(pos.z());
       det = ((pos.perp() > minR) && (pos.perp() < maxR)) ? 4 : 3;
 #ifdef printDebug
       ++detNull_[det - 3];
 #endif
     }
 #ifdef printDEBUG
     edm::LogVerbatim("HcalSim") << "Position " << pos.perp() << ":" << std::abs(pos.z()) << " Limits "
                                 << !(hcalConstants_->isHE()) << ":" << maxZ_ << " det " << det;
   } else {
     ++detNull_[3];
 #endif
   }
 
   if (numberingFromDDD.get()) {
     //get the ID's as eta, phi, depth, ... indices
     HcalNumberingFromDDD::HcalID tmp =
         numberingFromDDD->unitID(det, math::XYZVectorD(pos.x(), pos.y(), pos.z()), depth, lay);
     id = setDetUnitId(tmp);
   }
   return id;
 }

◆ setDetUnitId() [3/3]

uint32_t HCalSD::setDetUnitId ( HcalNumberingFromDDD::HcalID & tmp )

private

Definition at line 644 of file HCalSD.cc.

References TauDecayModes::dec, l1ctLayer2EG_cff::id, m_HcalTestNS, modifyDepth(), numberingScheme, convertSQLiteXML::ok, testNumber, and createJobs::tmp.

                                                              {
   modifyDepth(tmp);
   uint32_t id = (numberingScheme.get()) ? numberingScheme->getUnitID(tmp) : 0;
   if ((!testNumber) && m_HcalTestNS.get()) {
     bool ok = m_HcalTestNS->compare(tmp, id);
     if (!ok)
       edm::LogWarning("HcalSim") << "Det ID from HCalSD " << HcalDetId(id) << " " << std::hex << id << std::dec
                                  << " does not match one from relabller for " << tmp.subdet << ":" << tmp.etaR << ":"
                                  << tmp.phi << ":" << tmp.phis << ":" << tmp.depth << ":" << tmp.lay << std::endl;
   }
   return id;
 }

◆ setNumberingScheme()

void HCalSD::setNumberingScheme ( HcalNumberingScheme * scheme )

Definition at line 569 of file HCalSD.cc.

References numberingScheme, and generator_cfi::scheme.

Referenced by HCalSD(), HcalTestAnalysis::update(), SimG4HcalValidation::update(), and HcalTB04Analysis::update().

                                                            {
   if (scheme != nullptr) {
     edm::LogVerbatim("HcalSim") << "HCalSD: updates numbering scheme for " << GetName();
     numberingScheme.reset(scheme);
   }
 }

◆ update() [1/6]

void CaloSD::update

overrideprotected

Definition at line 750 of file CaloSD.cc.

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

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

◆ update() [2/6]

void CaloSD::update

overrideprotected

Definition at line 779 of file CaloSD.cc.

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

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

◆ update() [3/6]

void CaloSD::update

overrideprotected

Definition at line 1012 of file CaloSD.cc.

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

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

◆ update() [4/6]

void CaloSD::update

overrideprotected

Definition at line 758 of file CaloSD.cc.

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

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

◆ update() [5/6]

void CaloSD::update

overrideprotected

Definition at line 748 of file CaloSD.cc.

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

748 { initRun(); }

HCalSD::initRun

void initRun() override

Definition: HCalSD.cc:578

◆ update() [6/6]

void HCalSD::update ( const BeginOfJob * )

overrideprotectedvirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob *>.

Definition at line 576 of file HCalSD.cc.

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

576 {}

Member Data Documentation

◆ agingFlagHB

bool HCalSD::agingFlagHB

private

Definition at line 108 of file HCalSD.h.

Referenced by getEnergyDeposit(), and HCalSD().

◆ agingFlagHE

bool HCalSD::agingFlagHE

private

Definition at line 108 of file HCalSD.h.

Referenced by getEnergyDeposit(), and HCalSD().

◆ applyFidCut

bool HCalSD::applyFidCut

private

Definition at line 112 of file HCalSD.h.

Referenced by HCalSD(), and isItinFidVolume().

◆ betaThr

double HCalSD::betaThr

private

Definition at line 111 of file HCalSD.h.

Referenced by HCalSD().

◆ birk1

double HCalSD::birk1

private

Definition at line 111 of file HCalSD.h.

Referenced by getEnergyDeposit(), and HCalSD().

◆ birk2

double HCalSD::birk2

private

Definition at line 111 of file HCalSD.h.

Referenced by getEnergyDeposit(), and HCalSD().

◆ birk3

double HCalSD::birk3

private

Definition at line 111 of file HCalSD.h.

Referenced by getEnergyDeposit(), and HCalSD().

◆ deliveredLumi

double HCalSD::deliveredLumi

private

Definition at line 114 of file HCalSD.h.

Referenced by getEnergyDeposit(), getFromLibrary(), and HCalSD().

◆ depth_

int HCalSD::depth_

private

Definition at line 116 of file HCalSD.h.

Referenced by getEnergyDeposit(), getFromLibrary(), and modifyDepth().

◆ detNull_

std::vector<int> HCalSD::detNull_

private

Definition at line 126 of file HCalSD.h.

Referenced by endEvent(), initEvent(), and setDetUnitId().

◆ dist_

TH1F * HCalSD::dist_[9]

private

Definition at line 125 of file HCalSD.h.

Referenced by HCalSD(), and plotProfile().

◆ eminHitHB

double HCalSD::eminHitHB

private

Definition at line 113 of file HCalSD.h.

Referenced by filterHit(), and HCalSD().

◆ eminHitHE

double HCalSD::eminHitHE

private

Definition at line 113 of file HCalSD.h.

Referenced by filterHit(), and HCalSD().

◆ eminHitHF

double HCalSD::eminHitHF

private

Definition at line 113 of file HCalSD.h.

Referenced by filterHit(), and HCalSD().

◆ eminHitHO

double HCalSD::eminHitHO

private

Definition at line 113 of file HCalSD.h.

Referenced by filterHit(), and HCalSD().

◆ fibre1LV

std::vector<const G4LogicalVolume*> HCalSD::fibre1LV

private

Definition at line 123 of file HCalSD.h.

Referenced by HCalSD(), and isItStraightBundle().

◆ fibre2LV

std::vector<const G4LogicalVolume*> HCalSD::fibre2LV

private

Definition at line 123 of file HCalSD.h.

Referenced by HCalSD(), and isItConicalBundle().

◆ fibreLV

std::vector<const G4LogicalVolume*> HCalSD::fibreLV

private

Definition at line 123 of file HCalSD.h.

Referenced by HCalSD(), and isItFibre().

◆ fibreNames

std::vector<std::string> HCalSD::fibreNames

private

Definition at line 120 of file HCalSD.h.

Referenced by HCalSD(), and isItFibre().

◆ gpar

std::vector<double> HCalSD::gpar

private

Definition at line 117 of file HCalSD.h.

Referenced by HCalSD(), and plotHF().

◆ hcalConstants_

const HcalDDDSimConstants* HCalSD::hcalConstants_

private

Definition at line 96 of file HCalSD.h.

Referenced by getEnergyDeposit(), HCalSD(), modifyDepth(), and setDetUnitId().

◆ hcalSimConstants_

const HcalSimulationConstants* HCalSD::hcalSimConstants_

private

Definition at line 97 of file HCalSD.h.

Referenced by HCalSD().

◆ hfLevels

std::vector<int> HCalSD::hfLevels

private

Definition at line 118 of file HCalSD.h.

Referenced by HCalSD(), and isItHF().

◆ hfLV

std::vector<const G4LogicalVolume*> HCalSD::hfLV

private

Definition at line 123 of file HCalSD.h.

Referenced by HCalSD(), and isItHF().

◆ hfNames

std::vector<std::string> HCalSD::hfNames

private

Definition at line 119 of file HCalSD.h.

Referenced by HCalSD(), and isItHF().

◆ hfshower

std::unique_ptr<HFShower> HCalSD::hfshower

private

Definition at line 91 of file HCalSD.h.

Referenced by HCalSD(), and hitForFibre().

◆ hit_

TH1F* HCalSD::hit_[9]

private

Definition at line 125 of file HCalSD.h.

Referenced by HCalSD(), and plotProfile().

◆ hzvem

TH1F * HCalSD::hzvem

private

Definition at line 125 of file HCalSD.h.

Referenced by HCalSD(), and plotHF().

◆ hzvhad

TH1F * HCalSD::hzvhad

private

Definition at line 125 of file HCalSD.h.

Referenced by HCalSD(), and plotHF().

◆ isHF

bool HCalSD::isHF

private

Definition at line 107 of file HCalSD.h.

Referenced by getEnergyDeposit(), and getFromLibrary().

◆ layerWeights

std::map<uint32_t, double> HCalSD::layerWeights

private

Definition at line 124 of file HCalSD.h.

Referenced by layerWeight(), and readWeightFromFile().

◆ m_HBDarkening

const HBHEDarkening* HCalSD::m_HBDarkening

private

Definition at line 98 of file HCalSD.h.

Referenced by getEnergyDeposit().

◆ m_HcalTestNS

std::unique_ptr<HcalTestNS> HCalSD::m_HcalTestNS

private

Definition at line 101 of file HCalSD.h.

Referenced by HCalSD(), and setDetUnitId().

◆ m_HEDarkening

const HBHEDarkening* HCalSD::m_HEDarkening

private

Definition at line 99 of file HCalSD.h.

Referenced by getEnergyDeposit().

◆ m_HFDarkening

std::unique_ptr<HFDarkening> HCalSD::m_HFDarkening

private

Definition at line 100 of file HCalSD.h.

Referenced by getFromLibrary(), and HCalSD().

◆ materials

std::vector<const G4Material*> HCalSD::materials

private

Definition at line 122 of file HCalSD.h.

Referenced by HCalSD(), and isItScintillator().

◆ matNames

std::vector<std::string> HCalSD::matNames

private

Definition at line 121 of file HCalSD.h.

Referenced by HCalSD().

◆ maxRoff_

constexpr double HCalSD::maxRoff_ = 450.0

staticprivate

Definition at line 105 of file HCalSD.h.

Referenced by setDetUnitId().

◆ maxZ_

constexpr double HCalSD::maxZ_ = 10000.0

staticprivate

Definition at line 103 of file HCalSD.h.

Referenced by setDetUnitId().

◆ minRoff_

constexpr double HCalSD::minRoff_ = -1500.0

staticprivate

Definition at line 104 of file HCalSD.h.

Referenced by setDetUnitId().

◆ neutralDensity

bool HCalSD::neutralDensity

private

Definition at line 110 of file HCalSD.h.

Referenced by getEnergyDeposit(), and HCalSD().

◆ numberingFromDDD

std::unique_ptr<HcalNumberingFromDDD> HCalSD::numberingFromDDD

private

Definition at line 88 of file HCalSD.h.

Referenced by getHitFibreBundle(), getHitPMT(), HCalSD(), layerWeight(), and setDetUnitId().

◆ numberingScheme

std::unique_ptr<HcalNumberingScheme> HCalSD::numberingScheme

private

Definition at line 89 of file HCalSD.h.

Referenced by HCalSD(), setDetUnitId(), and setNumberingScheme().

◆ pmtLV

std::vector<const G4LogicalVolume*> HCalSD::pmtLV

private

Definition at line 123 of file HCalSD.h.

Referenced by HCalSD(), and isItPMT().

◆ showerBundle

std::unique_ptr<HFShowerFibreBundle> HCalSD::showerBundle

private

Definition at line 94 of file HCalSD.h.

Referenced by getEnergyDeposit(), getHitFibreBundle(), and HCalSD().

◆ showerLibrary

std::unique_ptr<HFShowerLibrary> HCalSD::showerLibrary

private

Definition at line 90 of file HCalSD.h.

Referenced by getFromHFLibrary(), and HCalSD().

◆ showerParam

std::unique_ptr<HFShowerParam> HCalSD::showerParam

private

Definition at line 92 of file HCalSD.h.

Referenced by getFromParam(), and HCalSD().

◆ showerPMT

std::unique_ptr<HFShowerPMT> HCalSD::showerPMT

private

Definition at line 93 of file HCalSD.h.

Referenced by getEnergyDeposit(), getHitPMT(), and HCalSD().

◆ slopeHE_

constexpr double HCalSD::slopeHE_ = 0.4

staticprivate

Definition at line 106 of file HCalSD.h.

Referenced by setDetUnitId().

◆ testNS_

bool HCalSD::testNS_

private

Definition at line 110 of file HCalSD.h.

Referenced by HCalSD().

◆ testNumber

bool HCalSD::testNumber

private

Definition at line 110 of file HCalSD.h.

Referenced by getEnergyDeposit(), HCalSD(), modifyDepth(), and setDetUnitId().

◆ time_

TH1F * HCalSD::time_[9]

private

Definition at line 125 of file HCalSD.h.

Referenced by HCalSD(), and plotProfile().

◆ useBirk

bool HCalSD::useBirk

private

Definition at line 109 of file HCalSD.h.

Referenced by getEnergyDeposit(), and HCalSD().

◆ useFibreBundle

bool HCalSD::useFibreBundle

private

Definition at line 109 of file HCalSD.h.

Referenced by getEnergyDeposit(), and HCalSD().

◆ useHF

bool HCalSD::useHF

private

Definition at line 112 of file HCalSD.h.

Referenced by HCalSD().

◆ useLayerWt

bool HCalSD::useLayerWt

private

Definition at line 109 of file HCalSD.h.

Referenced by getEnergyDeposit(), HCalSD(), and readWeightFromFile().

◆ useParam

bool HCalSD::useParam

private

Definition at line 112 of file HCalSD.h.

Referenced by getFromLibrary(), and HCalSD().

◆ usePMTHit

bool HCalSD::usePMTHit

private

Definition at line 109 of file HCalSD.h.

Referenced by getEnergyDeposit(), and HCalSD().

◆ useShowerLibrary

bool HCalSD::useShowerLibrary

private

Definition at line 112 of file HCalSD.h.

Referenced by getEnergyDeposit(), getFromLibrary(), and HCalSD().

◆ weight_

double HCalSD::weight_

private

Definition at line 115 of file HCalSD.h.

Referenced by getEnergyDeposit(), getFromHFLibrary(), getFromLibrary(), getFromParam(), and hitForFibre().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HCalSD()

◆ ~HCalSD()

Member Function Documentation

◆ endEvent()

◆ fillLogVolumeVector()

◆ filterHit()

◆ getEnergyDeposit()

◆ getFromHFLibrary()

◆ getFromLibrary()

◆ getFromParam()

◆ getHitFibreBundle()

◆ getHitPMT()

◆ hitForFibre()

◆ initEvent()

◆ initRun()

◆ isItConicalBundle()

◆ isItFibre() [1/2]

◆ isItFibre() [2/2]

◆ isItHF() [1/2]

◆ isItHF() [2/2]

◆ isItinFidVolume()

◆ isItPMT()

◆ isItScintillator()

◆ isItStraightBundle()

◆ layerWeight()

◆ modifyDepth()

◆ plotHF()

◆ plotProfile()

◆ printVolume()

◆ readWeightFromFile()

◆ setDetUnitId() [1/3]

◆ setDetUnitId() [2/3]

◆ setDetUnitId() [3/3]

◆ setNumberingScheme()

◆ update() [1/6]

◆ update() [2/6]

◆ update() [3/6]

◆ update() [4/6]

◆ update() [5/6]

◆ update() [6/6]

Member Data Documentation

◆ agingFlagHB

◆ agingFlagHE

◆ applyFidCut

◆ betaThr

◆ birk1

◆ birk2

◆ birk3

◆ deliveredLumi

◆ depth_

◆ detNull_

◆ dist_

◆ eminHitHB

◆ eminHitHE

◆ eminHitHF

◆ eminHitHO

◆ fibre1LV

◆ fibre2LV

◆ fibreLV

◆ fibreNames

◆ gpar

◆ hcalConstants_

◆ hcalSimConstants_

◆ hfLevels

◆ hfLV

◆ hfNames

◆ hfshower

◆ hit_

◆ hzvem

◆ hzvhad

◆ isHF

◆ layerWeights