#include <ZdcSD.h>

Inheritance diagram for ZdcSD:

Public Member Functions
uint32_t	setDetUnitId (const G4Step *step) override

void	setNumberingScheme (ZdcNumberingScheme *scheme)

	ZdcSD (const std::string &, const SensitiveDetectorCatalog &, edm::ParameterSet const &, const SimTrackManager *)

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

void	clear () override

void	clearHits () override

void	DrawAll () override

void	EndOfEvent (G4HCofThisEvent *eventHC) override

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

void	Initialize (G4HCofThisEvent *HCE) override

bool	isItFineCalo (const G4VTouchable *touch)

void	PrintAll () override

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)

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

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

void	Initialize (G4HCofThisEvent *eventHC) override

bool	isCaloSD () const

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

	~SensitiveDetector () override

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

void	slotForUpdate (const BeginOfRun * iT)

virtual	~Observer ()

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

void	slotForUpdate (const BeginOfEvent * iT)

virtual	~Observer ()

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

void	slotForUpdate (const BeginOfTrack * iT)

virtual	~Observer ()

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

void	slotForUpdate (const EndOfTrack * iT)

virtual	~Observer ()

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

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

Protected Member Functions
double	getEnergyDeposit (const G4Step *) override

bool	getFromLibrary (const G4Step *) override

void	initRun () override

Protected Member Functions inherited from CaloSD
bool	checkHit ()

CaloG4Hit *	createNewHit (const G4Step , const G4Track )

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

double	getResponseWt (const G4Track *)

virtual int	getTrackID (const G4Track *)

bool	hitExists (const G4Step *)

void	ignoreRejection ()

virtual void	initEvent (const BeginOfEvent *)

void	printDetectorLevels (const G4VTouchable *) const

void	processHit (const G4Step *step)

void	resetForNewPrimary (const G4Step *)

void	setNumberCheckedHits (int val)

void	setParameterized (bool val)

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

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

virtual int	setTrackID (const G4Step *)

void	setUseMap (bool val)

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

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 Attributes
std::vector< ZdcShowerLibrary::Hit >	hits

std::unique_ptr< ZdcNumberingScheme >	numberingScheme

std::unique_ptr< ZdcShowerLibrary >	showerLibrary

double	thFibDir

bool	useShowerHits

bool	useShowerLibrary

int	verbosity

double	zdcHitEnergyCut

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

CaloHitID	currentID

float	edepositEM

float	edepositHAD

double	eminHit

double	energyCut

G4ThreeVector	entranceLocal

G4ThreeVector	entrancePoint

bool	forceSave

float	incidentEnergy

double	kmaxIon

double	kmaxNeutron

double	kmaxProton

G4ThreeVector	posGlobal

CaloHitID	previousID

bool	suppressHeavy

double	tmaxHit

Detailed Description

Definition at line 15 of file ZdcSD.h.

Constructor & Destructor Documentation

◆ ZdcSD()

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

Definition at line 30 of file ZdcSD.cc.

References edm::ParameterSet::getParameter(), Skims_PA_cff::name, AlCaHLTBitMon_ParallelJobs::p, setNumberingScheme(), CaloSD::setParameterized(), showerLibrary, thFibDir, useShowerHits, useShowerLibrary, verbosity, and zdcHitEnergyCut.

     : CaloSD(name, clg, p, manager) {
   edm::ParameterSet m_ZdcSD = p.getParameter<edm::ParameterSet>("ZdcSD");
   useShowerLibrary = m_ZdcSD.getParameter<bool>("UseShowerLibrary");
   useShowerHits = m_ZdcSD.getParameter<bool>("UseShowerHits");
   zdcHitEnergyCut = m_ZdcSD.getParameter<double>("ZdcHitEnergyCut") * GeV;
   thFibDir = m_ZdcSD.getParameter<double>("FiberDirection");
   verbosity = m_ZdcSD.getParameter<int>("Verbosity");
   int verbn = verbosity / 10;
   verbosity %= 10;
   setNumberingScheme(new ZdcNumberingScheme(verbn));
 
   edm::LogVerbatim("ForwardSim") << "***************************************************\n"
                                  << "*                                                 *\n"
                                  << "* Constructing a ZdcSD  with name " << name << "   *\n"
                                  << "*                                                 *\n"
                                  << "***************************************************";
 
   edm::LogVerbatim("ForwardSim") << "\nUse of shower library is set to " << useShowerLibrary
                                  << "\nUse of Shower hits method is set to " << useShowerHits;
 
   edm::LogVerbatim("ForwardSim") << "\nEnergy Threshold Cut set to " << zdcHitEnergyCut / GeV << " (GeV)";
 
   if (useShowerLibrary) {
     showerLibrary = std::make_unique<ZdcShowerLibrary>(name, p);
     setParameterized(true);
   } else {
     showerLibrary.reset(nullptr);
   }
 }

◆ ~ZdcSD()

ZdcSD::~ZdcSD ( )

overridedefault

Member Function Documentation

◆ getEnergyDeposit()

double ZdcSD::getEnergyDeposit ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 120 of file ZdcSD.cc.

References a, funct::abs(), HLT_2023v12_cff::beta, ALCARECOTkAlJpsiMuMu_cff::charge, ztail::d, PVValHelper::eta, f, ForwardName::getName(), dqm-mbProfile::log, SiStripPI::max, SiStripPI::min, PbPb_ZMuSkimMuonDPG_cff::particleType, phi, pi, alignCSCRings::r, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), theta(), and thFibDir.

                                                   {
   double NCherPhot = 0.;
 
   // preStepPoint information
   G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
 
   const G4ThreeVector& hitPoint = preStepPoint->GetPosition();
   const G4ThreeVector& hit_mom = preStepPoint->GetMomentumDirection();
   G4double stepL = aStep->GetStepLength() / cm;
   G4double beta = preStepPoint->GetBeta();
   G4double charge = preStepPoint->GetCharge();
   if (charge == 0.0) 
     return 0.0;
 
   // theTrack information
   G4Track* theTrack = aStep->GetTrack();
   G4ThreeVector localPoint = theTrack->GetTouchable()->GetHistory()->GetTopTransform().TransformPoint(hitPoint);
 
 #ifdef EDM_ML_DEBUG
   const G4ThreeVector& vert_mom = theTrack->GetVertexMomentumDirection();
 
   // calculations
   float costheta =
       vert_mom.z() / sqrt(vert_mom.x() * vert_mom.x() + vert_mom.y() * vert_mom.y() + vert_mom.z() * vert_mom.z());
   float theta = std::acos(std::min(std::max(costheta, -1.f), 1.f));
   float eta = -std::log(std::tan(theta * 0.5f));
   float phi = -100.;
   if (vert_mom.x() != 0)
     phi = std::atan2(vert_mom.y(), vert_mom.x());
   if (phi < 0.)
     phi += twopi;
 
   // Get the total energy deposit
   double stepE = aStep->GetTotalEnergyDeposit();
 
   // postStepPoint information
   G4StepPoint* postStepPoint = aStep->GetPostStepPoint();
   G4VPhysicalVolume* postPV = postStepPoint->GetPhysicalVolume();
   G4VPhysicalVolume* currentPV = preStepPoint->GetPhysicalVolume();
   std::string nameVolume = ForwardName::getName(currentPV->GetName());
   std::string postnameVolume = ForwardName::getName(postPV->GetName());
   G4String particleType = theTrack->GetDefinition()->GetParticleName();
   edm::LogVerbatim("ForwardSim") << "ZdcSD::  getEnergyDeposit: \n"
                                  << "  preStepPoint: " << nameVolume << "," << stepL << "," << stepE << "," << beta
                                  << "," << charge << "\n"
                                  << "  postStepPoint: " << postnameVolume << "," << costheta << "," << theta << ","
                                  << eta << "," << phi << "," << particleType << " id= " << theTrack->GetTrackID()
                                  << " Etot(GeV)= " << theTrack->GetTotalEnergy() / GeV;
 #endif
   const double bThreshold = 0.67;
   if (beta > bThreshold) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("ForwardSim") << "ZdcSD::  getEnergyDeposit:  pass ";
 #endif
     const float nMedium = 1.4925;
     // float photEnSpectrDL = 10714.285714;
     //       photEnSpectrDL = (1./400.nm-1./700.nm)*10000000.cm/nm; /* cm-1  */
 
     const float photEnSpectrDE = 1.24;
     // 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
 
     const float effPMTandTransport = 0.15;
 
     // Check these values
     const float thFullRefl = 23.;
     float thFullReflRad = thFullRefl * pi / 180.;
 
     float thFibDirRad = thFibDir * pi / 180.;
 
     // at which theta the point is located:
     //   float th1 = hitPoint.theta();
 
     // theta of charged particle in LabRF(hit momentum direction):
     float costh = hit_mom.z() / sqrt(hit_mom.x() * hit_mom.x() + hit_mom.y() * hit_mom.y() + hit_mom.z() * hit_mom.z());
     float th = acos(std::min(std::max(costh, -1.f), 1.f));
     // just in case (can do both standard ranges of phi):
     if (th < 0.)
       th += twopi;
 
     // theta of cone with Cherenkov photons w.r.t.direction of charged part.:
     float costhcher = 1. / (nMedium * beta);
     float thcher = acos(std::min(std::max(costhcher, -1.f), 1.f));
 
     // diff thetas of charged part. and quartz direction in LabRF:
     float DelFibPart = std::abs(th - thFibDirRad);
 
     // define real distances:
     float d = std::abs(std::tan(th) - std::tan(thFibDirRad));
 
     float a = std::tan(thFibDirRad) + std::tan(std::abs(thFibDirRad - thFullReflRad));
     float r = std::tan(th) + std::tan(std::abs(th - thcher));
 
     // define losses d_qz in cone of full reflection inside quartz direction
     float d_qz = -1;
 #ifdef EDM_ML_DEBUG
     float variant = -1;
 #endif
     // if (d > (r+a))
     if (DelFibPart > (thFullReflRad + thcher)) {
 #ifdef EDM_ML_DEBUG
       variant = 0.;
 #endif
       d_qz = 0.;
     } else {
       // if ((DelFibPart + thcher) < thFullReflRad )  [(d+r) < a]
       if ((th + thcher) < (thFibDirRad + thFullReflRad) && (th - thcher) > (thFibDirRad - thFullReflRad)) {
 #ifdef EDM_ML_DEBUG
         variant = 1.;
 #endif
         d_qz = 1.;
       } else {
         // if ((thcher - DelFibPart ) > thFullReflRad )  [(r-d) > a]
         if ((thFibDirRad + thFullReflRad) < (th + thcher) && (thFibDirRad - thFullReflRad) > (th - thcher)) {
 #ifdef EDM_ML_DEBUG
           variant = 2.;
 #endif
           d_qz = 0.;
         } else {
 #ifdef EDM_ML_DEBUG
           variant = 3.;  // d_qz is calculated below
 #endif
           // use crossed length of circles(cone projection) - dC1/dC2 :
           float arg_arcos = 0.;
           float 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);
           float th_arcos = acos(std::min(std::max(arg_arcos, -1.f), 1.f));
           d_qz = th_arcos / twopi;
           d_qz = std::abs(d_qz);
 #ifdef EDM_ML_DEBUG
           edm::LogVerbatim("ForwardSim") << "  d_qz: " << r << "," << a << "," << d << " " << tan_arcos << " "
                                          << arg_arcos;
           edm::LogVerbatim("ForwardSim") << "," << arg_arcos;
           edm::LogVerbatim("ForwardSim") << " " << d_qz;
           edm::LogVerbatim("ForwardSim") << " " << th_arcos;
           edm::LogVerbatim("ForwardSim") << "," << d_qz;
 #endif
         }
       }
     }
     double meanNCherPhot = 0.;
     int poissNCherPhot = 0;
     if (d_qz > 0) {
       meanNCherPhot = 370. * charge * charge * (1. - 1. / (nMedium * nMedium * beta * beta)) * photEnSpectrDE * stepL;
 
       poissNCherPhot = std::max((int)G4Poisson(meanNCherPhot), 0);
       NCherPhot = poissNCherPhot * effPMTandTransport * d_qz;
     }
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("ForwardSim") << "ZdcSD::  getEnergyDeposit:  gED: " << stepE << "," << costh << "," << th << ","
                                    << costhcher << "," << thcher << "," << DelFibPart << "," << d << "," << a << ","
                                    << r << "," << hitPoint << "," << hit_mom << "," << vert_mom << "," << localPoint
                                    << "," << charge << "," << beta << "," << stepL << "," << d_qz << "," << variant
                                    << "," << meanNCherPhot << "," << poissNCherPhot << "," << NCherPhot;
 #endif
     // --constants-----------------
     // << "," << photEnSpectrDE
     // << "," << nMedium
     // << "," << bThreshold
     // << "," << thFibDirRad
     // << "," << thFullReflRad
     // << "," << effPMTandTransport
     // --other variables-----------
     // << "," << curprocess
     // << "," << nameProcess
     // << "," << name
     // << "," << rad
     // << "," << mat
 
   } else {
     // determine failure mode: beta, charge, and/or nameVolume
     if (beta <= bThreshold)
       edm::LogVerbatim("ForwardSim") << "ZdcSD::  getEnergyDeposit: fail beta=" << beta;
 
 #ifdef EDM_ML_DEBUG
     std::string nameVolume = ForwardName::getName(currentPV->GetName());
     if (!(nameVolume == "ZDC_EMFiber" || nameVolume == "ZDC_HadFiber"))
       edm::LogVerbatim("ForwardSim") << "ZdcSD::  getEnergyDeposit: fail nv=" << nameVolume;
 #endif
   }
 
   return NCherPhot;
 }

◆ getFromLibrary()

bool ZdcSD::getFromLibrary ( const G4Step * aStep )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 66 of file ZdcSD.cc.

References CaloSD::currentHit, CaloSD::currentID, CaloSD::edepositEM, CaloSD::edepositHAD, CaloSD::entranceLocal, CaloSD::entrancePoint, CaloG4Hit::getEM(), CaloG4Hit::getEntryLocal(), CaloG4Hit::getHadr(), CaloG4Hit::getIncidentEnergy(), CaloG4Hit::getTrackID(), CaloG4Hit::getUnitID(), hits, mps_fire::i, CaloSD::incidentEnergy, convertSQLiteXML::ok, CaloSD::posGlobal, CaloSD::processHit(), CaloSD::resetForNewPrimary(), CaloHitID::setID(), CaloSD::setTrackID(), showerLibrary, protons_cff::time, and zdcHitEnergyCut.

                                               {
   bool ok = true;
 
   auto const preStepPoint = aStep->GetPreStepPoint();
 
   double etrack = preStepPoint->GetKineticEnergy();
   int primaryID = setTrackID(aStep);
 
   hits.clear();
 
   // Reset entry point for new primary
   resetForNewPrimary(aStep);
 
   if (etrack >= zdcHitEnergyCut) {
     // create hits only if above threshold
 
 #ifdef EDM_ML_DEBUG
     auto const theTrack = aStep->GetTrack();
     edm::LogVerbatim("ForwardSim") << "----------------New track------------------------------\n"
                                    << "Incident EnergyTrack: " << etrack << " MeV \n"
                                    << "Zdc Cut Energy for Hits: " << zdcHitEnergyCut << " MeV \n"
                                    << "ZdcSD::getFromLibrary " << hits.size() << " hits for " << GetName() << " of "
                                    << primaryID << " with " << theTrack->GetDefinition()->GetParticleName() << " of "
                                    << etrack << " MeV\n";
 #endif
     hits.swap(showerLibrary.get()->getHits(aStep, ok));
   }
 
   incidentEnergy = etrack;
   entrancePoint = preStepPoint->GetPosition();
   for (unsigned int i = 0; i < hits.size(); i++) {
     posGlobal = hits[i].position;
     entranceLocal = hits[i].entryLocal;
     double time = hits[i].time;
     unsigned int unitID = hits[i].detID;
     edepositHAD = hits[i].DeHad;
     edepositEM = hits[i].DeEM;
     currentID.setID(unitID, time, primaryID, 0);
     processHit(aStep);
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("ForwardSim") << "ZdcSD: Final Hit number:" << i << "-->"
                                    << "New HitID: " << currentHit->getUnitID()
                                    << " New Hit trackID: " << currentHit->getTrackID()
                                    << " New EM Energy: " << currentHit->getEM() / GeV
                                    << " New HAD Energy: " << currentHit->getHadr() / GeV
                                    << " New HitEntryPoint: " << currentHit->getEntryLocal()
                                    << " New IncidentEnergy: " << currentHit->getIncidentEnergy() / GeV
                                    << " New HitPosition: " << posGlobal;
 #endif
   }
   return ok;
 }

◆ initRun()

void ZdcSD::initRun ( )

overrideprotectedvirtual

Reimplemented from CaloSD.

Definition at line 64 of file ZdcSD.cc.

References hits.

64 { hits.clear(); }

ZdcSD::hits

std::vector< ZdcShowerLibrary::Hit > hits

Definition: ZdcSD.h:38

◆ setDetUnitId()

uint32_t ZdcSD::setDetUnitId ( const G4Step * step )

overridevirtual

Implements CaloSD.

Definition at line 309 of file ZdcSD.cc.

References numberingScheme.

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

◆ setNumberingScheme()

void ZdcSD::setNumberingScheme ( ZdcNumberingScheme * scheme )

Definition at line 313 of file ZdcSD.cc.

References numberingScheme, and generator_cfi::scheme.

Referenced by ZdcSD().

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

Member Data Documentation

◆ hits

std::vector<ZdcShowerLibrary::Hit> ZdcSD::hits

private

Definition at line 38 of file ZdcSD.h.

Referenced by getFromLibrary(), and initRun().

◆ numberingScheme

std::unique_ptr<ZdcNumberingScheme> ZdcSD::numberingScheme

private

Definition at line 37 of file ZdcSD.h.

Referenced by setDetUnitId(), and setNumberingScheme().

◆ showerLibrary

std::unique_ptr<ZdcShowerLibrary> ZdcSD::showerLibrary

private

Definition at line 36 of file ZdcSD.h.

Referenced by getFromLibrary(), and ZdcSD().

◆ thFibDir

double ZdcSD::thFibDir

private

Definition at line 33 of file ZdcSD.h.

Referenced by getEnergyDeposit(), and ZdcSD().

◆ useShowerHits

bool ZdcSD::useShowerHits

private

Definition at line 32 of file ZdcSD.h.

Referenced by ZdcSD().

◆ useShowerLibrary

bool ZdcSD::useShowerLibrary

private

Definition at line 32 of file ZdcSD.h.

Referenced by ZdcSD().

◆ verbosity

int ZdcSD::verbosity

private

Definition at line 31 of file ZdcSD.h.

Referenced by ZdcSD().

◆ zdcHitEnergyCut

double ZdcSD::zdcHitEnergyCut

private

Definition at line 34 of file ZdcSD.h.

Referenced by getFromLibrary(), and ZdcSD().

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ZdcSD()

◆ ~ZdcSD()

Member Function Documentation

◆ getEnergyDeposit()

◆ getFromLibrary()

◆ initRun()

◆ setDetUnitId()

◆ setNumberingScheme()

Member Data Documentation

◆ hits

◆ numberingScheme

◆ showerLibrary

◆ thFibDir

◆ useShowerHits

◆ useShowerLibrary

◆ verbosity

◆ zdcHitEnergyCut