#include <ZdcShowerLibrary.h>

Classes
struct	Hit

Public Member Functions
int	encodePartID (G4int parCode)

int	getEnergyFromLibrary (const G4ThreeVector &posHit, const G4ThreeVector &momDir, double energy, G4int parCode, HcalZDCDetId::Section section, bool side, int channel)

std::vector< Hit > &	getHits (const G4Step *aStep, bool &ok)

void	initRun (G4ParticleTable *theParticleTable)

int	photonFluctuation (double eav, double esig, double edis)

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

	~ZdcShowerLibrary ()=default

Private Attributes
G4int	anuePDG

G4int	anumuPDG

G4int	anutauPDG

G4int	emPDG

G4int	epPDG

G4int	etaPDG

G4int	gammaPDG

G4int	geantinoPDG

std::vector< ZdcShowerLibrary::Hit >	hits

int	npe

G4int	nuePDG

G4int	numuPDG

G4int	nutauPDG

G4int	pi0PDG

bool	verbose

Detailed Description

Definition at line 22 of file ZdcShowerLibrary.h.

Constructor & Destructor Documentation

◆ ZdcShowerLibrary()

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

Definition at line 19 of file ZdcShowerLibrary.cc.

References edm::ParameterSet::getUntrackedParameter(), hits, npe, AlCaHLTBitMon_ParallelJobs::p, and verbose.

                                                                                   {
   edm::ParameterSet m_HS = p.getParameter<edm::ParameterSet>("ZdcShowerLibrary");
   verbose = m_HS.getUntrackedParameter<int>("Verbosity", 0);
 
   npe = 9;  // number of channels or fibers where the energy will be deposited
   hits.reserve(npe);
 }

◆ ~ZdcShowerLibrary()

ZdcShowerLibrary::~ZdcShowerLibrary ( )

default

Member Function Documentation

◆ encodePartID()

int ZdcShowerLibrary::encodePartID ( G4int parCode )

Definition at line 272 of file ZdcShowerLibrary.cc.

References emPDG, epPDG, and gammaPDG.

Referenced by getEnergyFromLibrary(), and getHits().

                                                 {
   G4int iparCode = 1;
   if (parCode == emPDG || parCode == epPDG || parCode == gammaPDG) {
     iparCode = 0;
   } else {
     return iparCode;
   }
   return iparCode;
 }

◆ getEnergyFromLibrary()

int ZdcShowerLibrary::getEnergyFromLibrary	(	const G4ThreeVector &	posHit,
		const G4ThreeVector &	momDir,
		double	energy,
		G4int	parCode,
		HcalZDCDetId::Section	section,
		bool	side,
		int	channel
	)

Definition at line 148 of file ZdcShowerLibrary.cc.

References encodePartID(), hcalRecHitTable_cff::energy, fact, HMesonGammaMonitor_cff::nphotons, photonFluctuation(), funct::pow(), hgcalPlots::section, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), theXChannelBoundaries, and ecalPiZeroTask_cfi::X0.

Referenced by getHits().

                                                         {
   int nphotons = -1;
 
   energy = energy / CLHEP::GeV;
 
   edm::LogVerbatim("ZdcShower") << "\n ZdcShowerLibrary::getEnergyFromLibrary input/output variables:"
                                 << " phi: " << momDir.phi() / CLHEP::deg << " theta: " << momDir.theta() / CLHEP::deg
                                 << " xin : " << hitPoint.x() << " yin : " << hitPoint.y() << " zin : " << hitPoint.z()
                                 << " track en: " << energy << "(GeV)"
                                 << " section: " << section << " side: " << side << " channel: " << channel
                                 << " partID: " << parCode;
 
   // these varables are not used for now
   //float phi   = momDir.phi() / CLHEP::deg;
   //float theta = momDir.theta() / CLHEP::deg;
   //float zin = hitPoint.z();
   //int isection = int(section);
   //int iside = (side)? 1 : 2;
 
   int iparCode = encodePartID(parCode);
 
   double eav = 0.;
   double esig = 0.;
   double edis = 0.;
 
   float xin = hitPoint.x();
   float yin = hitPoint.y();
   float fact = 0.;
 
   if (section == 1 && iparCode != 0) {
     if (channel < 5)
       if (((theXChannelBoundaries[channel - 1]) < (xin + X0)) && ((xin + X0) <= theXChannelBoundaries[channel]))
         fact = 0.18;
     if (channel == 5)
       if (theXChannelBoundaries[channel - 1] < xin + X0)
         fact = 0.18;
   }
 
   if (section == 2 && iparCode != 0) {
     if (channel == 1)
       fact = 0.34;
     if (channel == 2)
       fact = 0.24;
     if (channel == 3)
       fact = 0.17;
     if (channel == 4)
       fact = 0.07;
   }
   if (section == 1 && iparCode == 0) {
     if (channel < 5)
       if (((theXChannelBoundaries[channel - 1]) < (xin + X0)) && ((xin + X0) <= theXChannelBoundaries[channel]))
         fact = 1.;
     if (channel == 5)
       if (theXChannelBoundaries[channel - 1] < xin + X0)
         fact = 1.0;
   }
 
   //change to cm for parametrization
   yin = yin / CLHEP::cm;
   xin = xin / CLHEP::cm;
 
   if (iparCode == 0) {
     eav = ((((((-0.0002 * xin - 2.0e-13) * xin + 0.0022) * xin + 1.0e-11) * xin - 0.0217) * xin - 3.0e-10) * xin +
            1.0028) *
           (((0.0001 * yin + 0.0056) * yin + 0.0508) * yin + 1.0) * 300.0 * pow((energy / 300.0), 0.99);  // EM
     esig = ((((((0.0005 * xin - 1.0e-12) * xin - 0.0052) * xin + 5.0e-11) * xin + 0.032) * xin - 2.0e-10) * xin + 1.0) *
            (((0.0006 * yin + 0.0071) * yin - 0.031) * yin + 1.0) * 30.0 * pow((energy / 300.0), 0.54);  // EM
     edis = 1.0;
   } else {
     eav = ((((((-0.0002 * xin - 2.0e-13) * xin + 0.0022) * xin + 1.0e-11) * xin - 0.0217) * xin - 3.0e-10) * xin +
            1.0028) *
           (((0.0001 * yin + 0.0056) * yin + 0.0508) * yin + 1.0) * 300.0 * pow((energy / 300.0), 1.12);  // HD
     esig = ((((((0.0005 * xin - 1.0e-12) * xin - 0.0052) * xin + 5.0e-11) * xin + 0.032) * xin - 2.0e-10) * xin + 1.0) *
            (((0.0006 * yin + 0.0071) * yin - 0.031) * yin + 1.0) * 54.0 * pow((energy / 300.0), 0.93);  //HD
     edis = 3.0;
   }
 
   if (eav < 0. || edis < 0.) {
     edm::LogVerbatim("ZdcShower") << " Negative everage energy from parametrization \n"
                                   << " xin: " << xin << "(cm)"
                                   << " yin: " << yin << "(cm)"
                                   << " track en: " << energy << "(GeV)"
                                   << " eaverage: " << eav / CLHEP::GeV << " (GeV)"
                                   << " esigma: " << esig / CLHEP::GeV << "  (GeV)"
                                   << " edist: " << edis << " (GeV)"
                                   << " dE hit: " << nphotons / CLHEP::GeV << " (GeV)";
     return 0;
   }
 
   // Convert from GeV to MeV for the code
   eav = eav * CLHEP::GeV;
   esig = esig * CLHEP::GeV;
 
   while (nphotons == -1 || nphotons > int(eav + 5. * esig))
     nphotons = static_cast<int>(fact * photonFluctuation(eav, esig, edis));
 
   edm::LogVerbatim("ZdcShower") << " track en: " << energy << "(GeV)"
                                 << " eaverage: " << eav / CLHEP::GeV << " (GeV)"
                                 << " esigma: " << esig / CLHEP::GeV << "  (GeV)"
                                 << " edist: " << edis << " (GeV)"
                                 << " dE hit: " << nphotons / CLHEP::GeV << " (GeV)";
 
   return nphotons;
 }

◆ getHits()

std::vector< ZdcShowerLibrary::Hit > & ZdcShowerLibrary::getHits	(	const G4Step *	aStep,
		bool &	ok
	)

Definition at line 48 of file ZdcShowerLibrary.cc.

References ZdcShowerLibrary::Hit::DeEM, ZdcShowerLibrary::Hit::DeHad, hcalRecHitTable_cff::depth, ZdcShowerLibrary::Hit::depth, ZdcShowerLibrary::Hit::detID, HcalZDCDetId::EM, encodePartID(), hcalRecHitTable_cff::energy, ZdcShowerLibrary::Hit::entryLocal, geantinoPDG, getEnergyFromLibrary(), HcalZDCDetId::HAD, hits, mps_fire::i, npe, convertSQLiteXML::ok, ZdcShowerLibrary::Hit::position, hgcalPlots::section, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), theXChannelBoundaries, theZHadChannelBoundaries, theZSectionBoundaries, ZdcShowerLibrary::Hit::time, HLT_2024v13_cff::track, funct::true, ecalPiZeroTask_cfi::X0, geometryCSVtoXML::xx, Y0, geometryCSVtoXML::yy, Z0, and geometryCSVtoXML::zz.

                                                                                        {
   G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
   G4StepPoint* postStepPoint = aStep->GetPostStepPoint();
   G4Track* track = aStep->GetTrack();
 
   const G4DynamicParticle* aParticle = track->GetDynamicParticle();
   const G4ThreeVector& momDir = aParticle->GetMomentumDirection();
   double energy = preStepPoint->GetKineticEnergy();
   G4ThreeVector hitPoint = preStepPoint->GetPosition();
   const G4ThreeVector& hitPointOrig = preStepPoint->GetPosition();
   G4int parCode = track->GetDefinition()->GetPDGEncoding();
 
   hits.clear();
 
   ok = false;
   if (parCode == geantinoPDG)
     return hits;
   ok = true;
 
   G4ThreeVector pos;
   G4ThreeVector posLocal;
   double tSlice = (postStepPoint->GetGlobalTime()) / CLHEP::nanosecond;
 
   int nHit = 0;
   HcalZDCDetId::Section section;
   bool side = false;
   int channel = 0;
   double xx, yy, zz;
   double xxlocal, yylocal, zzlocal;
 
   ZdcShowerLibrary::Hit oneHit;
   side = (hitPointOrig.z() > 0.) ? true : false;
 
   float xWidthEM = fabs(theXChannelBoundaries[0] - theXChannelBoundaries[1]);
   float zWidthEM = fabs(theZSectionBoundaries[0] - theZSectionBoundaries[1]);
   float zWidthHAD = fabs(theZHadChannelBoundaries[0] - theZHadChannelBoundaries[1]);
 
   for (int i = 0; i < npe; i++) {
     if (i < 5) {
       section = HcalZDCDetId::EM;
       channel = i + 1;
       xxlocal = theXChannelBoundaries[i] + (xWidthEM / 2.);
       xx = xxlocal + X0;
       yy = 0.0;
       yylocal = yy + Y0;
       zzlocal = theZSectionBoundaries[0] + (zWidthEM / 2.);
       zz = (hitPointOrig.z() > 0.) ? zzlocal + Z0 : zzlocal - Z0;
       pos = G4ThreeVector(xx, yy, zz);
       posLocal = G4ThreeVector(xxlocal, yylocal, zzlocal);
     }
     if (i > 4) {
       section = HcalZDCDetId::HAD;
       channel = i - 4;
       xxlocal = 0.0;
       xx = xxlocal + X0;
       yylocal = 0;
       yy = yylocal + Y0;
       zzlocal = (hitPointOrig.z() > 0.) ? theZHadChannelBoundaries[i - 5] + (zWidthHAD / 2.)
                                         : theZHadChannelBoundaries[i - 5] - (zWidthHAD / 2.);
       zz = (hitPointOrig.z() > 0.) ? zzlocal + Z0 : zzlocal - Z0;
       pos = G4ThreeVector(xx, yy, zz);
       posLocal = G4ThreeVector(xxlocal, yylocal, zzlocal);
     }
 
     oneHit.position = pos;
     oneHit.entryLocal = posLocal;
     oneHit.depth = channel;
     oneHit.time = tSlice;
     oneHit.detID = HcalZDCDetId(section, side, channel);
 
     // Note: coodinates of hit are relative to center of detector (X0,Y0,Z0)
     hitPoint.setX(hitPointOrig.x() - X0);
     hitPoint.setY(hitPointOrig.y() - Y0);
     double setZ = (hitPointOrig.z() > 0.) ? hitPointOrig.z() - Z0 : fabs(hitPointOrig.z()) - Z0;
     hitPoint.setZ(setZ);
 
     int dE = getEnergyFromLibrary(hitPoint, momDir, energy, parCode, section, side, channel);
 
     int iparCode = encodePartID(parCode);
     if (iparCode == 0) {
       oneHit.DeEM = dE;
       oneHit.DeHad = 0.;
     } else {
       oneHit.DeEM = 0;
       oneHit.DeHad = dE;
     }
 
     hits.push_back(oneHit);
 
     edm::LogVerbatim("ZdcShower") << "\nZdcShowerLibrary:Generated Hit " << nHit << " orig hit pos " << hitPointOrig
                                   << " orig hit pos local coord" << hitPoint << " new position "
                                   << (hits[nHit].position) << " Channel " << (hits[nHit].depth) << " side " << side
                                   << " Time " << (hits[nHit].time) << " DetectorID " << (hits[nHit].detID)
                                   << " Had Energy " << (hits[nHit].DeHad) << " EM Energy  " << (hits[nHit].DeEM)
                                   << "\n";
     nHit++;
   }
   return hits;
 }

◆ initRun()

void ZdcShowerLibrary::initRun ( G4ParticleTable * theParticleTable )

Definition at line 27 of file ZdcShowerLibrary.cc.

References anuePDG, anumuPDG, anutauPDG, emPDG, epPDG, etaPDG, gammaPDG, geantinoPDG, nuePDG, numuPDG, nutauPDG, and pi0PDG.

                                                                 {
   G4String parName;
   emPDG = theParticleTable->FindParticle(parName = "e-")->GetPDGEncoding();
   epPDG = theParticleTable->FindParticle(parName = "e+")->GetPDGEncoding();
   gammaPDG = theParticleTable->FindParticle(parName = "gamma")->GetPDGEncoding();
   pi0PDG = theParticleTable->FindParticle(parName = "pi0")->GetPDGEncoding();
   etaPDG = theParticleTable->FindParticle(parName = "eta")->GetPDGEncoding();
   nuePDG = theParticleTable->FindParticle(parName = "nu_e")->GetPDGEncoding();
   numuPDG = theParticleTable->FindParticle(parName = "nu_mu")->GetPDGEncoding();
   nutauPDG = theParticleTable->FindParticle(parName = "nu_tau")->GetPDGEncoding();
   anuePDG = theParticleTable->FindParticle(parName = "anti_nu_e")->GetPDGEncoding();
   anumuPDG = theParticleTable->FindParticle(parName = "anti_nu_mu")->GetPDGEncoding();
   anutauPDG = theParticleTable->FindParticle(parName = "anti_nu_tau")->GetPDGEncoding();
   geantinoPDG = theParticleTable->FindParticle(parName = "geantino")->GetPDGEncoding();
   edm::LogVerbatim("ZdcShower") << "ZdcShowerLibrary: Particle codes for e- = " << emPDG << ", e+ = " << epPDG
                                 << ", gamma = " << gammaPDG << ", pi0 = " << pi0PDG << ", eta = " << etaPDG
                                 << ", geantino = " << geantinoPDG << "\n        nu_e = " << nuePDG
                                 << ", nu_mu = " << numuPDG << ", nu_tau = " << nutauPDG << ", anti_nu_e = " << anuePDG
                                 << ", anti_nu_mu = " << anumuPDG << ", anti_nu_tau = " << anutauPDG;
 }

◆ photonFluctuation()

int ZdcShowerLibrary::photonFluctuation	(	double	eav,
		double	esig,
		double	edis
	)

Definition at line 259 of file ZdcShowerLibrary.cc.

References createfilelist::int.

Referenced by getEnergyFromLibrary().

                                                                             {
   int nphot = 0;
   double efluct = 0.;
   if (edis == 1.0)
     efluct = eav + esig * CLHEP::RandGaussQ::shoot();
   if (edis == 3.0)
     efluct = eav + esig * CLHEP::RandLandau::shoot();
   nphot = int(efluct);
   if (nphot < 0)
     nphot = 0;
   return nphot;
 }

Member Data Documentation

◆ anuePDG

G4int ZdcShowerLibrary::anuePDG

private

Definition at line 55 of file ZdcShowerLibrary.h.

Referenced by initRun().

◆ anumuPDG

G4int ZdcShowerLibrary::anumuPDG

private

Definition at line 55 of file ZdcShowerLibrary.h.

Referenced by initRun().

◆ anutauPDG

G4int ZdcShowerLibrary::anutauPDG

private

Definition at line 55 of file ZdcShowerLibrary.h.

Referenced by initRun().

◆ emPDG

G4int ZdcShowerLibrary::emPDG

private

Definition at line 53 of file ZdcShowerLibrary.h.

Referenced by encodePartID(), and initRun().

◆ epPDG

G4int ZdcShowerLibrary::epPDG

private

Definition at line 53 of file ZdcShowerLibrary.h.

Referenced by encodePartID(), and initRun().

◆ etaPDG

G4int ZdcShowerLibrary::etaPDG

private

Definition at line 54 of file ZdcShowerLibrary.h.

Referenced by initRun().

◆ gammaPDG

G4int ZdcShowerLibrary::gammaPDG

private

Definition at line 53 of file ZdcShowerLibrary.h.

Referenced by encodePartID(), and initRun().

◆ geantinoPDG

G4int ZdcShowerLibrary::geantinoPDG

private

Definition at line 55 of file ZdcShowerLibrary.h.

Referenced by getHits(), and initRun().

◆ hits

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

private

Definition at line 57 of file ZdcShowerLibrary.h.

Referenced by getHits(), and ZdcShowerLibrary().

◆ npe

int ZdcShowerLibrary::npe

private

Definition at line 56 of file ZdcShowerLibrary.h.

Referenced by getHits(), and ZdcShowerLibrary().

◆ nuePDG

G4int ZdcShowerLibrary::nuePDG

private

Definition at line 54 of file ZdcShowerLibrary.h.

Referenced by initRun().

◆ numuPDG

G4int ZdcShowerLibrary::numuPDG

private

Definition at line 54 of file ZdcShowerLibrary.h.

Referenced by initRun().

◆ nutauPDG

G4int ZdcShowerLibrary::nutauPDG

private

Definition at line 54 of file ZdcShowerLibrary.h.

Referenced by initRun().

◆ pi0PDG

G4int ZdcShowerLibrary::pi0PDG

private

Definition at line 54 of file ZdcShowerLibrary.h.

Referenced by initRun().

◆ verbose

bool ZdcShowerLibrary::verbose

private

Definition at line 52 of file ZdcShowerLibrary.h.

Referenced by ZdcShowerLibrary().

Classes

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ZdcShowerLibrary()

◆ ~ZdcShowerLibrary()

Member Function Documentation

◆ encodePartID()

◆ getEnergyFromLibrary()

◆ getHits()

◆ initRun()

◆ photonFluctuation()

Member Data Documentation

◆ anuePDG

◆ anumuPDG

◆ anutauPDG

◆ emPDG

◆ epPDG

◆ etaPDG

◆ gammaPDG

◆ geantinoPDG

◆ hits

◆ npe

◆ nuePDG

◆ numuPDG

◆ nutauPDG

◆ pi0PDG

◆ verbose