#include <SimG4CMS/Forward/interface/CastorSD.h>

Inheritance diagram for CastorSD:

Public Member Functions
	CastorSD (const std::string &, const SensitiveDetectorCatalog &clg, edm::ParameterSet const &, const SimTrackManager *)

uint32_t	setDetUnitId (const G4Step *step) override

void	setNumberingScheme (CastorNumberingScheme *scheme)

	~CastorSD () override

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

Protected Member Functions
double	getEnergyDeposit (const G4Step *) override

bool	getFromLibrary (const G4Step *) override

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

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

virtual void	endEvent ()

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

virtual bool	filterHit (CaloG4Hit *, double)

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

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

virtual uint16_t	getDepth (const G4Step *)

int	getNumberOfHits (int k=0)

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

virtual int	getTrackID (const G4Track *)

bool	hitExists (const G4Step *, int k)

void	ignoreRejection ()

virtual void	initEvent (const BeginOfEvent *)

virtual void	initRun ()

void	printDetectorLevels (const G4VTouchable *) const

void	processHit (const G4Step *step)

virtual void	processSecondHit (const G4Step , const G4Track )

void	resetForNewPrimary (const G4Step *)

void	setNumberCheckedHits (int val, int k=0)

void	setParameterized (bool val)

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

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

virtual int	setTrackID (const G4Step *)

void	setUseMap (bool val)

std::string	shortreprID (const CaloHitID &ID)

std::string	shortreprID (const CaloG4Hit *hit)

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

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

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

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

void	update (const ::EndOfEvent *) override

void	updateHit (CaloG4Hit *, int k)

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

Local3DPoint	ConvertToLocal3DPoint (const G4ThreeVector &point) const

Local3DPoint	FinalStepPosition (const G4Step *step, coordinates) const

Local3DPoint	InitialStepPosition (const G4Step *step, coordinates) const

Local3DPoint	LocalPostStepPosition (const G4Step *step) const

Local3DPoint	LocalPreStepPosition (const G4Step *step) const

void	NaNTrap (const G4Step *step) const

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

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

Private Member Functions
uint32_t	rotateUnitID (uint32_t, const G4Track *, const CastorShowerEvent &)

Private Attributes
double	energyThresholdSL

G4LogicalVolume *	lvC3EF

G4LogicalVolume *	lvC3HF

G4LogicalVolume *	lvC4EF

G4LogicalVolume *	lvC4HF

G4LogicalVolume *	lvCAST

double	non_compensation_factor

CastorNumberingScheme *	numberingScheme

CastorShowerLibrary *	showerLibrary

bool	useShowerLibrary

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

Description: Stores hits of Castor in appropriate container

Usage: Used in sensitive detector builder

Definition at line 33 of file CastorSD.h.

Constructor & Destructor Documentation

◆ CastorSD()

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

Definition at line 31 of file CastorSD.cc.

References g4SimHits_cfi::CastorShowerLibrary, energyThresholdSL, edm::ParameterSet::getParameter(), lvC3EF, lvC3HF, lvC4EF, lvC4HF, lvCAST, Skims_PA_cff::name, non_compensation_factor, AlCaHLTBitMon_ParallelJobs::p, setNumberingScheme(), CaloSD::setParameterized(), showerLibrary, and useShowerLibrary.

     : CaloSD(name, clg, p, manager),
       numberingScheme(nullptr),
       lvC3EF(nullptr),
       lvC3HF(nullptr),
       lvC4EF(nullptr),
       lvC4HF(nullptr),
       lvCAST(nullptr) {
   edm::ParameterSet m_CastorSD = p.getParameter<edm::ParameterSet>("CastorSD");
   useShowerLibrary = m_CastorSD.getParameter<bool>("useShowerLibrary");
   energyThresholdSL = m_CastorSD.getParameter<double>("minEnergyInGeVforUsingSLibrary");
   energyThresholdSL = energyThresholdSL * CLHEP::GeV;  //  Convert GeV => MeV
 
   non_compensation_factor = m_CastorSD.getParameter<double>("nonCompensationFactor");
 
   if (useShowerLibrary) {
     showerLibrary = new CastorShowerLibrary(name, p);
     setParameterized(true);
   }
   setNumberingScheme(new CastorNumberingScheme());
 
   edm::LogVerbatim("ForwardSim") << "********************************************************\n"
                                  << "* Constructing a CastorSD  with name " << GetName() << "\n"
                                  << "* Using Castor Shower Library: " << useShowerLibrary << "\n"
                                  << "********************************************************";
 
   const G4LogicalVolumeStore* lvs = G4LogicalVolumeStore::GetInstance();
   for (auto lv : *lvs) {
     if (strcmp((lv->GetName()).c_str(), "C3EF") == 0) {
       lvC3EF = lv;
     }
     if (strcmp((lv->GetName()).c_str(), "C3HF") == 0) {
       lvC3HF = lv;
     }
     if (strcmp((lv->GetName()).c_str(), "C4EF") == 0) {
       lvC4EF = lv;
     }
     if (strcmp((lv->GetName()).c_str(), "C4HF") == 0) {
       lvC4HF = lv;
     }
     if (strcmp((lv->GetName()).c_str(), "CAST") == 0) {
       lvCAST = lv;
     }
     if (lvC3EF != nullptr && lvC3HF != nullptr && lvC4EF != nullptr && lvC4HF != nullptr && lvCAST != nullptr) {
       break;
     }
   }
   edm::LogVerbatim("ForwardSim") << "CastorSD:: LogicalVolume pointers\n"
                                  << lvC3EF << " for C3EF; " << lvC3HF << " for C3HF; " << lvC4EF << " for C4EF; "
                                  << lvC4HF << " for C4HF; " << lvCAST << " for CAST. ";
 }

◆ ~CastorSD()

CastorSD::~CastorSD ( )

override

Definition at line 88 of file CastorSD.cc.

References showerLibrary.

88 { delete showerLibrary; }

CastorSD::showerLibrary

CastorShowerLibrary * showerLibrary

Definition: CastorSD.h:47

Member Function Documentation

◆ getEnergyDeposit()

double CastorSD::getEnergyDeposit ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 92 of file CastorSD.cc.

References a, funct::abs(), HLT_2023v12_cff::beta, ALCARECOTkAlJpsiMuMu_cff::charge, ztail::d, TrackInformation::getCastorHitPID(), TrackInformation::hasCastorHit(), G4TrackToParticleID::isStableHadronIon(), lvC3EF, lvC3HF, lvC4EF, lvC4HF, SiStripPI::max, SiStripPI::min, non_compensation_factor, pi, alignCSCRings::r, pfClustersFromCombinedCaloHF_cfi::scale, TrackInformation::setCastorHitPID(), mathSSE::sqrt(), and funct::tan().

                                                      {
   double NCherPhot = 0.;
 
   // Get theTrack
   auto const theTrack = aStep->GetTrack();
 
   // preStepPoint information *********************************************
   auto const preStepPoint = aStep->GetPreStepPoint();
   auto const currentPV = preStepPoint->GetPhysicalVolume();
   auto const currentLV = currentPV->GetLogicalVolume();
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("ForwardSim") << "CastorSD::getEnergyDeposit:"
                                  << "\n CurrentStepNumber , TrackID , ParentID, Particle , VertexPosition ,"
                                  << " LogicalVolumeAtVertex , PV, Time"
                                  << "\n  TRACKINFO: " << theTrack->GetCurrentStepNumber() << " , "
                                  << theTrack->GetTrackID() << " , " << theTrack->GetParentID() << " , "
                                  << theTrack->GetDefinition()->GetParticleName() << " , "
                                  << theTrack->GetVertexPosition() << " , "
                                  << theTrack->GetLogicalVolumeAtVertex()->GetName() << " , " << currentPV->GetName()
                                  << " , " << theTrack->GetGlobalTime();
 #endif
 
   // if particle moves from interaction point or "backwards (halo)
   const G4ThreeVector& hit_mom = preStepPoint->GetMomentumDirection();
   const G4ThreeVector& hitPoint = preStepPoint->GetPosition();
   double zint = hitPoint.z();
 
   // Check if theTrack is a muon (if so, DO NOT use Shower Library)
   G4int parCode = theTrack->GetDefinition()->GetPDGEncoding();
   bool isHad = G4TrackToParticleID::isStableHadronIon(theTrack);
 
   double meanNCherPhot = 0;
   G4double charge = preStepPoint->GetCharge();
   // VI: no Cerenkov light from neutrals
   if (0.0 == charge) {
     return meanNCherPhot;
   }
 
   G4double beta = preStepPoint->GetBeta();
   const double bThreshold = 0.67;
   // VI: no Cerenkov light from non-relativistic particles
   if (beta < bThreshold) {
     return meanNCherPhot;
   }
 
   // remember primary particle hitting the CASTOR detector
   TrackInformationExtractor TIextractor;
   TrackInformation& trkInfo = TIextractor(theTrack);
   if (!trkInfo.hasCastorHit()) {
     trkInfo.setCastorHitPID(parCode);
   }
   G4double stepl = aStep->GetStepLength() / CLHEP::cm;
 
   //int castorHitPID = trkInfo.hasCastorHit() ? std::abs(trkInfo.getCastorHitPID())
   //  : std::abs(parCode);
 
   // *************************************************
   // take into account light collection curve for plate
   //      double weight = curve_Castor(nameVolume, preStepPoint);
   //      double edep   = aStep->GetTotalEnergyDeposit() * weight;
   // *************************************************
 
   // *************************************************
   /*    comments for sensitive volumes:      
         C001 ...-... CP06,CBU1 ...-...CALM --- > fibres and bundle 
         for first release of CASTOR
         CASF  --- > quartz plate  for first and second releases of CASTOR  
         GF2Q, GFNQ, GR2Q, GRNQ 
         for tests with my own test geometry of HF (on ask of Gavrilov)
         C3TF, C4TF - for third release of CASTOR
   */
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("ForwardSim") << "CastorSD::getEnergyDeposit: for ID=" << theTrack->GetTrackID()
                                  << " LV: " << currentLV->GetName() << " isHad:" << isHad << " pdg=" << parCode
                                  << " castorPID=" << trkInfo.getCastorHitPID() << " sl= " << stepl
                                  << " Edep= " << aStep->GetTotalEnergyDeposit();
 #endif
   if (currentLV == lvC3EF || currentLV == lvC4EF || currentLV == lvC3HF || currentLV == lvC4HF) {
     const double nMedium = 1.4925;
     //     double photEnSpectrDL = (1./400.nm-1./700.nm)*10000000.cm/nm; /* cm-1  */
     //     double photEnSpectrDL = 10714.285714;
 
     const double photEnSpectrDE = 1.24; /* see below   */
     /*     E = 2pi*(1./137.)*(eV*cm/370.)/lambda  =     */
     /*       = 12.389184*(eV*cm)/lambda                 */
     /*     Emax = 12.389184*(eV*cm)/400nm*10-7cm/nm  = 3.01 eV     */
     /*     Emin = 12.389184*(eV*cm)/700nm*10-7cm/nm  = 1.77 eV     */
     /*     delE = Emax - Emin = 1.24 eV  --> */
     /*   */
     /* default for Castor nameVolume  == "CASF" or (C3TF & C4TF)  */
 
     const double thFullRefl = 23.; /* 23.dergee */
     const double thFullReflRad = thFullRefl * pi / 180.;
 
     /* default for Castor nameVolume  == "CASF" or (C3TF & C4TF)  */
     const double thFibDir = 45.; /* .dergee */
     /* for test HF geometry volumes:   
        if(nameVolume == "GF2Q" || nameVolume == "GFNQ" ||
        nameVolume == "GR2Q" || nameVolume == "GRNQ")
        thFibDir = 0.0; // .dergee
     */
     const double thFibDirRad = thFibDir * pi / 180.;
 
     // theta of charged particle in LabRF(hit momentum direction):
     double costh =
         hit_mom.z() / sqrt(hit_mom.x() * hit_mom.x() + hit_mom.y() * hit_mom.y() + hit_mom.z() * hit_mom.z());
     if (zint < 0)
       costh = -costh;
     double th = acos(std::min(std::max(costh, double(-1.)), double(1.)));
 
     // just in case (can do bot use):
     if (th < 0.)
       th += twopi;
 
     // theta of cone with Cherenkov photons w.r.t.direction of charged part.:
     double costhcher = 1. / (nMedium * beta);
     double thcher = acos(std::min(std::max(costhcher, double(-1.)), double(1.)));
 
     // diff thetas of charged part. and quartz direction in LabRF:
     double DelFibPart = std::abs(th - thFibDirRad);
 
     // define real distances:
     double d = std::abs(tan(th) - tan(thFibDirRad));
 
     double a = tan(thFibDirRad) + tan(std::abs(thFibDirRad - thFullReflRad));
     double r = tan(th) + tan(std::abs(th - thcher));
 
     // define losses d_qz in cone of full reflection inside quartz direction
     double d_qz;
 #ifdef EDM_ML_DEBUG
     int variant(0);
 #endif
     if (DelFibPart > (thFullReflRad + thcher)) {
       d_qz = 0.;
     } else {
       if ((th + thcher) < (thFibDirRad + thFullReflRad) && (th - thcher) > (thFibDirRad - thFullReflRad)) {
         d_qz = 1.;
 #ifdef EDM_ML_DEBUG
         variant = 1;
 #endif
       } else {
         if ((thFibDirRad + thFullReflRad) < (th + thcher) && (thFibDirRad - thFullReflRad) > (th - thcher)) {
           d_qz = 0.;
 #ifdef EDM_ML_DEBUG
           variant = 2;
 #endif
         } else {
 #ifdef EDM_ML_DEBUG
           variant = 3;
 #endif
           // use crossed length of circles(cone projection)
           // dC1/dC2 :
           double arg_arcos = 0.;
           double tan_arcos = 2. * a * d;
           if (tan_arcos != 0.)
             arg_arcos = (r * r - a * a - d * d) / tan_arcos;
           arg_arcos = std::abs(arg_arcos);
           double th_arcos = acos(std::min(std::max(arg_arcos, -1.), 1.));
           d_qz = std::abs(th_arcos / CLHEP::twopi);
         }
       }
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("ForwardSim") << "variant" << variant;
 #endif
 
     // !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
     meanNCherPhot = 370. * charge * charge * (1. - 1. / (nMedium * nMedium * beta * beta)) * photEnSpectrDE * stepl;
 
     double scale = isHad ? non_compensation_factor : 1.0;
 
     int poissNCherPhot = std::max(0, (int)G4Poisson(meanNCherPhot * scale));
 
     const double effPMTandTransport = 0.19;
     const double ReflPower = 0.1;
     double proba = d_qz + (1 - d_qz) * ReflPower;
     NCherPhot = poissNCherPhot * effPMTandTransport * proba * 0.307;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("ForwardSim") << " Nph= " << NCherPhot << " Np= " << poissNCherPhot
                                    << " eff= " << effPMTandTransport << " pb= " << proba << " Nmean= " << meanNCherPhot
                                    << " q=" << charge << " beta=" << beta << " nMedium= " << nMedium << " sl= " << stepl
                                    << " Nde=" << photEnSpectrDE;
 #endif
   }
   return NCherPhot;
 }

◆ getFromLibrary()

bool CastorSD::getFromLibrary ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 361 of file CastorSD.cc.

References funct::abs(), CaloSD::currentID, dot(), CaloSD::edepositEM, CaloSD::edepositHAD, energyThresholdSL, funct::false, CastorShowerLibrary::getShowerHits(), TrackInformation::hasCastorHit(), hfClusterShapes_cfi::hits, mps_fire::i, G4TrackToParticleID::isGammaElectronPositron(), G4TrackToParticleID::isStableHadronIon(), lvCAST, M_PI, non_compensation_factor, CaloSD::processHit(), dttmaxenums::R, CaloSD::resetForNewPrimary(), rotateUnitID(), pfClustersFromCombinedCaloHF_cfi::scale, TrackInformation::setCastorHitPID(), CaloHitID::setID(), CaloSD::setTrackID(), showerLibrary, mathSSE::sqrt(), theta(), hcalRecHitTable_cff::time, funct::true, and useShowerLibrary.

                                                  {
   //
   //   Method to get hits from the Shower Library
   //
   //   "updateHit" save the Hits to a CaloG4Hit container
   //
 
   auto const theTrack = aStep->GetTrack();
   auto parCode = theTrack->GetDefinition()->GetPDGEncoding();
 
   // remember primary particle hitting the CASTOR detector
   TrackInformationExtractor TIextractor;
   TrackInformation& trkInfo = TIextractor(theTrack);
   if (!trkInfo.hasCastorHit()) {
     trkInfo.setCastorHitPID(parCode);
   }
 
   if (!useShowerLibrary) {
     return false;
   }
 
   // preStepPoint information *********************************************
 
   auto const preStepPoint = aStep->GetPreStepPoint();
   auto const currentPV = preStepPoint->GetPhysicalVolume();
   auto const currentLV = currentPV->GetLogicalVolume();
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("ForwardSim") << "CastorSD::getFromLibrary: for ID=" << theTrack->GetTrackID()
                                  << " parentID= " << theTrack->GetParentID() << " "
                                  << theTrack->GetDefinition()->GetParticleName() << " LV: " << currentLV->GetName()
                                  << " PV: " << currentPV->GetName() << "\n eta= " << theTrack->GetPosition().eta()
                                  << " phi= " << theTrack->GetPosition().phi()
                                  << " z(cm)= " << theTrack->GetPosition().z() / CLHEP::cm
                                  << " time(ns)= " << theTrack->GetGlobalTime()
                                  << " E(GeV)= " << theTrack->GetTotalEnergy() / CLHEP::GeV;
 
 #endif
 
   const G4ThreeVector& hitPoint = preStepPoint->GetPosition();
   const G4ThreeVector& hit_mom = preStepPoint->GetMomentumDirection();
   double zint = hitPoint.z();
   double pz = hit_mom.z();
 
   // Check if theTrack moves backward
   bool backward = (pz * zint < 0.) ? true : false;
 
   // Check that theTrack is above the energy threshold to use Shower Library
   bool aboveThreshold = (theTrack->GetKineticEnergy() > energyThresholdSL) ? true : false;
 
   // Check if theTrack is a muon (if so, DO NOT use Shower Library)
   bool isEM = G4TrackToParticleID::isGammaElectronPositron(parCode);
   bool isHad = G4TrackToParticleID::isStableHadronIon(theTrack);
 
   // angle condition
   double theta_max = M_PI - 3.1305;  // angle in radians corresponding to -5.2 eta
   double R_mom = sqrt(hit_mom.x() * hit_mom.x() + hit_mom.y() * hit_mom.y());
   double theta = atan2(R_mom, std::abs(pz));
   bool angleok = (theta < theta_max) ? true : false;
 
   // OkToUse
   double R = sqrt(hitPoint.x() * hitPoint.x() + hitPoint.y() * hitPoint.y());
   bool dot = (zint < -14450. && R < 45.) ? true : false;
   bool inRange = (zint < -14700. || R > 193.) ? false : true;
   //bool OkToUse = (inRange && !dot) ? true : false;
 
   bool particleWithinShowerLibrary =
       aboveThreshold && (isEM || isHad) && (!backward) && inRange && !dot && angleok && currentLV == lvCAST;
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("ForwardSim") << "CastorSD::getFromLibrary: ID= " << theTrack->GetTrackID() << " E>E0 "
                                  << aboveThreshold << " isEM " << isEM << " isHad " << isHad << " backword " << backward
                                  << " Ok " << (inRange && !dot) << " angle " << angleok
                                  << " LV: " << currentLV->GetName() << "  " << (currentLV == lvCAST) << " "
                                  << particleWithinShowerLibrary << " Edep= " << aStep->GetTotalEnergyDeposit();
 #endif
 
   // Use Castor shower library if energy is above threshold, is not a muon
   // and is not moving backward
   if (!particleWithinShowerLibrary) {
     return false;
   }
 
   // ****    Call method to retrieve hits from the ShowerLibrary   ****
   // always kill primary
   bool isKilled(true);
   CastorShowerEvent hits = showerLibrary->getShowerHits(aStep, isKilled);
 
   int primaryID = setTrackID(aStep);
 
   // Reset entry point for new primary
   resetForNewPrimary(aStep);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("ForwardSim") << "CastorSD::getFromLibrary:  " << hits.getNhit() << " hits for " << GetName()
                                  << " from " << theTrack->GetDefinition()->GetParticleName() << " of "
                                  << preStepPoint->GetKineticEnergy() / CLHEP::GeV << " GeV and trackID "
                                  << theTrack->GetTrackID() << " isHad: " << isHad;
 #endif
 
   // Scale to correct energy
   double E_track = preStepPoint->GetTotalEnergy();
   double E_SLhit = hits.getPrimE() * CLHEP::GeV;
   double scale = E_track / E_SLhit;
 
   //Non compensation
   if (isHad) {
     scale *= non_compensation_factor;  // if hadronic extend the scale with the non-compensation factor
   }
   //  Loop over hits retrieved from the library
   for (unsigned int i = 0; i < hits.getNhit(); ++i) {
     // Get nPhotoElectrons and set edepositEM / edepositHAD accordingly
     double nPhotoElectrons = hits.getNphotons(i) * scale;
 
     if (isEM) {
       edepositEM = nPhotoElectrons;
       edepositHAD = 0.;
     } else {
       edepositEM = 0.;
       edepositHAD = nPhotoElectrons;
     }
 
     // Get hit position and time
     double time = hits.getTime(i);
 
     // Get hit detID
     unsigned int unitID = hits.getDetID(i);
 
     // Make the detID "rotation" from one sector to another taking into account the
     // sectors of the impinging particle (theTrack) and of the particle that produced
     // the 'hits' retrieved from shower library
     unsigned int rotatedUnitID = rotateUnitID(unitID, theTrack, hits);
     currentID[0].setID(rotatedUnitID, time, primaryID, 0);
     processHit(aStep);
   }  //  End of loop over hits
   return isKilled;
 }

◆ rotateUnitID()

uint32_t CastorSD::rotateUnitID	(	uint32_t	unitID,
		const G4Track *	track,
		const CastorShowerEvent &	shower
	)

private

Definition at line 299 of file CastorSD.cc.

References printConversionInfo::aux, CastorShowerEvent::getPrimPhi(), createfilelist::int, M_PI, fileinputsource_cfi::sec, HLT_2023v12_cff::track, and reco::btau::trackPhi.

Referenced by getFromLibrary().

                                                                                                       {
   // ==============================================================
   //
   //   o   Exploit Castor phi symmetry to return newUnitID for
   //       shower hits based on track phi coordinate
   //
   // ==============================================================
 
   // Get 'track' phi:
   double trackPhi = track->GetPosition().phi();
   if (trackPhi < 0)
     trackPhi += 2 * M_PI;
   // Get phi from primary that gave rise to SL 'shower':
   double showerPhi = shower.getPrimPhi();
   if (showerPhi < 0)
     showerPhi += 2 * M_PI;
   // Delta phi:
 
   //  Find the OctSector for which 'track' and 'shower' belong
   int trackOctSector = (int)(trackPhi / (M_PI / 4));
   int showerOctSector = (int)(showerPhi / (M_PI / 4));
 
   uint32_t newUnitID;
   uint32_t sec = ((unitID >> 4) & 0xF);
   uint32_t complement = (unitID & 0xFFFFFF0F);
 
   // Get 'track' z:
   double trackZ = track->GetPosition().z();
 
   int aux;
   int dSec = 2 * (trackOctSector - showerOctSector);
   if (trackZ > 0)  // Good for revision 1.9 of CastorNumberingScheme
   {
     int sec1 = sec - dSec;
     //    sec -= dSec ;
     if (sec1 < 0)
       sec1 += 16;
     if (sec1 > 15)
       sec1 -= 16;
     sec = (uint32_t)(sec1);
   } else {
     if (dSec < 0)
       sec += 16;
     sec += dSec;
     aux = (int)(sec / 16);
     sec -= aux * 16;
   }
   sec = sec << 4;
   newUnitID = complement | sec;
 
 #ifdef EDM_ML_DEBUG
   if (newUnitID != unitID) {
     edm::LogVerbatim("ForwardSim") << "\n CastorSD::rotateUnitID:  "
                                    << "\n     unitID = " << unitID << "\n  newUnitID = " << newUnitID;
   }
 #endif
 
   return newUnitID;
 }

◆ setDetUnitId()

uint32_t CastorSD::setDetUnitId ( const G4Step * step )

overridevirtual

Implements CaloSD.

Definition at line 283 of file CastorSD.cc.

References CastorNumberingScheme::getUnitID(), and numberingScheme.

                                                    {
   return (numberingScheme == nullptr ? 0 : numberingScheme->getUnitID(aStep));
 }

◆ setNumberingScheme()

void CastorSD::setNumberingScheme ( CastorNumberingScheme * scheme )

Definition at line 289 of file CastorSD.cc.

References numberingScheme, and generator_cfi::scheme.

Referenced by CastorSD().

                                                                {
   if (scheme != nullptr) {
     edm::LogVerbatim("ForwardSim") << "CastorSD: updates numbering scheme for " << GetName();
     delete numberingScheme;
     numberingScheme = scheme;
   }
 }

Member Data Documentation

◆ energyThresholdSL

double CastorSD::energyThresholdSL

private

Definition at line 52 of file CastorSD.h.

Referenced by CastorSD(), and getFromLibrary().

◆ lvC3EF

G4LogicalVolume* CastorSD::lvC3EF

private

Definition at line 48 of file CastorSD.h.

Referenced by CastorSD(), and getEnergyDeposit().

◆ lvC3HF

G4LogicalVolume * CastorSD::lvC3HF

private

Definition at line 48 of file CastorSD.h.

Referenced by CastorSD(), and getEnergyDeposit().

◆ lvC4EF

G4LogicalVolume * CastorSD::lvC4EF

private

Definition at line 48 of file CastorSD.h.

Referenced by CastorSD(), and getEnergyDeposit().

◆ lvC4HF

G4LogicalVolume * CastorSD::lvC4HF

private

Definition at line 48 of file CastorSD.h.

Referenced by CastorSD(), and getEnergyDeposit().

◆ lvCAST

G4LogicalVolume* CastorSD::lvCAST

private

Definition at line 49 of file CastorSD.h.

Referenced by CastorSD(), and getFromLibrary().

◆ non_compensation_factor

double CastorSD::non_compensation_factor

private

Definition at line 53 of file CastorSD.h.

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

◆ numberingScheme

CastorNumberingScheme* CastorSD::numberingScheme

private

Definition at line 46 of file CastorSD.h.

Referenced by setDetUnitId(), and setNumberingScheme().

◆ showerLibrary

CastorShowerLibrary* CastorSD::showerLibrary

private

Definition at line 47 of file CastorSD.h.

Referenced by CastorSD(), getFromLibrary(), and ~CastorSD().

◆ useShowerLibrary

bool CastorSD::useShowerLibrary

private

Definition at line 51 of file CastorSD.h.

Referenced by CastorSD(), and getFromLibrary().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ CastorSD()

◆ ~CastorSD()

Member Function Documentation

◆ getEnergyDeposit()

◆ getFromLibrary()

◆ rotateUnitID()

◆ setDetUnitId()

◆ setNumberingScheme()

Member Data Documentation

◆ energyThresholdSL

◆ lvC3EF

◆ lvC3HF

◆ lvC4EF

◆ lvC4HF

◆ lvCAST

◆ non_compensation_factor

◆ numberingScheme

◆ showerLibrary

◆ useShowerLibrary