Inheritance diagram for HcalTB04Analysis:

Public Member Functions
	HcalTB04Analysis (const edm::ParameterSet &p)

	HcalTB04Analysis (const HcalTB04Analysis &)=delete

const HcalTB04Analysis &	operator= (const HcalTB04Analysis &)=delete

void	produce (edm::Event &, const edm::EventSetup &) override

	~HcalTB04Analysis () override

Public Member Functions inherited from SimProducer
const SimProducer &	operator= (const SimProducer &)=delete

void	registerProducts (edm::ProducesCollector producesCollector)

	SimProducer ()

	SimProducer (const SimProducer &)=delete

Public Member Functions inherited from SimWatcher
virtual void	beginRun (edm::EventSetup const &)

bool	isMT () const

const SimWatcher &	operator= (const SimWatcher &)=delete

virtual void	registerConsumes (edm::ConsumesCollector)

	SimWatcher ()

	SimWatcher (const SimWatcher &)=delete

virtual	~SimWatcher ()

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 EndOfEvent *>
	Observer ()

void	slotForUpdate (const EndOfEvent * iT)

virtual	~Observer ()

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

void	slotForUpdate (const G4Step * iT)

virtual	~Observer ()

Private Member Functions
void	clear ()

void	fillBuffer (const EndOfEvent *evt)

void	fillEvent (PHcalTB04Info &)

void	finalAnalysis ()

void	init ()

void	qieAnalysis (CLHEP::HepRandomEngine *)

double	scale (int det, int layer)

double	timeOfFlight (int det, int layer, double eta)

int	unitID (uint32_t id)

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

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

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

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

void	xtalAnalysis (CLHEP::HepRandomEngine *)

Private Attributes
G4RotationMatrix *	beamline_RM

const double	beamOffset

int	count

std::vector< CaloHit >	ecalHitCache

const double	ecalNoise

double	eecalq

double	eecals

double	ehcalq

double	ehcals

std::vector< double >	enois

std::vector< double >	eqeta

std::vector< double >	eqie

std::vector< double >	eqlay

std::vector< double >	eqphi

std::vector< double >	eseta

std::vector< double >	esime

std::vector< double >	esimh

std::vector< double >	eslay

std::vector< double >	esphi

double	etaInit

double	etotq

double	etots

int	evNum

std::vector< CaloHit >	hcalHitCache

std::vector< CaloHit >	hcalHitLayer

const bool	hcalOnly

HcalTB04Histo *	histo

int	iceta

int	icphi

std::vector< uint32_t >	idEcal

std::vector< int >	idHcal

std::vector< uint32_t >	idTower

std::vector< int >	idXtal

const edm::ParameterSet	m_Anal

const int	mode

HcalQie *	myQie

const std::vector< std::string >	names

int	nCrystal

int	nPrimary

int	nTower

int	particleType

double	phiInit

double	pInit

bool	pvFound

G4ThreeVector	pvMomentum

G4ThreeVector	pvPosition

int	pvType

G4ThreeVector	pvUVW

const double	scaleHB0

const double	scaleHB16

const double	scaleHE0

const double	scaleHO

std::vector< double >	secEkin

std::vector< G4ThreeVector >	secMomentum

std::vector< int >	secPartID

std::vector< int >	secTrackID

std::vector< int >	shortLivedSecondaries

const int	type

Additional Inherited Members
Protected Member Functions inherited from SimProducer
template<class T >
void	produces ()

template<class T >
void	produces (const std::string &instanceName)

Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 73 of file HcalTB04Analysis.cc.

Constructor & Destructor Documentation

◆ HcalTB04Analysis() [1/2]

HcalTB04Analysis::HcalTB04Analysis ( const edm::ParameterSet & p )

Definition at line 150 of file HcalTB04Analysis.cc.

References beamline_RM, beamOffset, ecalNoise, JetChargeProducer_cfi::exp, edm::ParameterSet::getParameter(), hcalOnly, histo, iceta, icphi, init(), m_Anal, mode, myQie, and AlCaHLTBitMon_ParallelJobs::p.

     : m_Anal(p.getParameter<edm::ParameterSet>("HcalTB04Analysis")),
       hcalOnly(m_Anal.getParameter<bool>("HcalOnly")),
       mode(m_Anal.getParameter<int>("Mode")),
       type(m_Anal.getParameter<int>("Type")),
       ecalNoise(m_Anal.getParameter<double>("EcalNoise")),
       beamOffset(-m_Anal.getParameter<double>("BeamPosition") * CLHEP::cm),
       scaleHB0(m_Anal.getParameter<double>("ScaleHB0")),
       scaleHB16(m_Anal.getParameter<double>("ScaleHB16")),
       scaleHO(m_Anal.getParameter<double>("ScaleHO")),
       scaleHE0(m_Anal.getParameter<double>("ScaleHE0")),
       names(m_Anal.getParameter<std::vector<std::string> >("Names")),
       myQie(nullptr),
       histo(nullptr) {
   double fMinEta = m_Anal.getParameter<double>("MinEta");
   double fMaxEta = m_Anal.getParameter<double>("MaxEta");
   double fMinPhi = m_Anal.getParameter<double>("MinPhi");
   double fMaxPhi = m_Anal.getParameter<double>("MaxPhi");
   double beamEta = (fMaxEta + fMinEta) / 2.;
   double beamPhi = (fMaxPhi + fMinPhi) / 2.;
   double beamThet = 2 * atan(exp(-beamEta));
   if (beamPhi < 0)
     beamPhi += twopi;
   iceta = static_cast<int>(beamEta / 0.087) + 1;
   icphi = static_cast<int>(std::fabs(beamPhi) / 0.087) + 5;
   if (icphi > 72)
     icphi -= 73;
 
   produces<PHcalTB04Info>();
 
   beamline_RM = new G4RotationMatrix;
   beamline_RM->rotateZ(-beamPhi);
   beamline_RM->rotateY(-beamThet);
 
   edm::LogVerbatim("HcalTBSim")
       << "HcalTB04:: Initialised as observer of BeginOf Job/BeginOfRun/BeginOfEvent/G4Step/EndOfEvent with Parameter "
          "values:\n \thcalOnly = "
       << hcalOnly << "\tecalNoise = " << ecalNoise << "\n\tMode = " << mode << " (0: HB2 Standard; 1:HB2 Segmented)"
       << "\tType = " << type << " (0: HB; 1 HE; 2 HB+HE)\n\tbeamOffset = " << beamOffset << "\ticeta = " << iceta
       << "\ticphi = " << icphi << "\n\tbeamline_RM = " << *beamline_RM;
 
   init();
 
   myQie = new HcalQie(p);
   histo = new HcalTB04Histo(m_Anal);
 }

◆ HcalTB04Analysis() [2/2]

HcalTB04Analysis::HcalTB04Analysis ( const HcalTB04Analysis & )

delete

◆ ~HcalTB04Analysis()

HcalTB04Analysis::~HcalTB04Analysis ( )

override

Definition at line 197 of file HcalTB04Analysis.cc.

References count, histo, and myQie.

                                     {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "\n -------->  Total number of selected entries : " << count << "\nPointers:: QIE "
                                 << myQie << " Histo " << histo;
 #endif
   if (myQie) {
     delete myQie;
     myQie = nullptr;
   }
   if (histo) {
     delete histo;
     histo = nullptr;
   }
 }

Member Function Documentation

◆ clear()

void HcalTB04Analysis::clear ( void )

private

Definition at line 1033 of file HcalTB04Analysis.cc.

References ecalHitCache, enois, eqie, esime, esimh, etaInit, hcalHitCache, hcalHitLayer, mps_fire::i, nCrystal, nPrimary, nTower, particleType, phiInit, pInit, pvFound, pvMomentum, pvPosition, pvType, pvUVW, secEkin, secMomentum, secPartID, secTrackID, and shortLivedSecondaries.

Referenced by init(), and update().

                              {
   pvFound = false;
   pvType = -2;
   pvPosition = G4ThreeVector();
   pvMomentum = G4ThreeVector();
   pvUVW = G4ThreeVector();
   secTrackID.clear();
   secPartID.clear();
   secMomentum.clear();
   secEkin.clear();
   shortLivedSecondaries.clear();
 
   ecalHitCache.erase(ecalHitCache.begin(), ecalHitCache.end());
   hcalHitCache.erase(hcalHitCache.begin(), hcalHitCache.end());
   hcalHitLayer.erase(hcalHitLayer.begin(), hcalHitLayer.end());
   nPrimary = particleType = 0;
   pInit = etaInit = phiInit = 0;
 
   esimh.clear();
   eqie.clear();
   esimh.reserve(nTower);
   eqie.reserve(nTower);
   for (int i = 0; i < nTower; i++) {
     esimh.push_back(0.);
     eqie.push_back(0.);
   }
   esime.clear();
   enois.clear();
   esime.reserve(nCrystal);
   enois.reserve(nCrystal);
   for (int i = 0; i < nCrystal; i++) {
     esime.push_back(0.);
     enois.push_back(0.);
   }
 }

◆ fillBuffer()

void HcalTB04Analysis::fillBuffer ( const EndOfEvent * evt )

private

Definition at line 481 of file HcalTB04Analysis.cc.

References TauDecayModes::dec, MillePedeFileConverter_cfg::e, ecalHitCache, PVValHelper::eta, etaInit, evNum, CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getEntry(), CaloG4Hit::getTimeSlice(), CaloG4Hit::getTrackID(), HcalTBNumberingScheme::getUnitID(), CaloG4Hit::getUnitID(), watchdog::group, hcalHitCache, hcalHitLayer, hfClusterShapes_cfi::hits, mps_fire::i, heavyIonCSV_trainingSettings::idx, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, dqmiolumiharvest::j, nano_mu_digi_cff::layer, dqm-mbProfile::log, SiStripPI::max, SiStripPI::min, mode, names, npart, nPrimary, AlCaHLTBitMon_ParallelJobs::p, particleType, phi, phiInit, pInit, conifer::pow(), multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, scale(), heavyFlavorDQMFirstStep_cff::sdName, jetUpdater_cfi::sort, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), theta(), hcalRecHitTable_cff::time, timeOfFlight(), unitID(), HcalTestNumbering::unpackHcalIndex(), and z.

Referenced by update().

                                                        {
   std::vector<CaloHit> hhits, hhitl;
   int idHC, j;
   CaloG4HitCollection* theHC;
   std::map<int, float, std::less<int> > primaries;
 #ifdef EDM_ML_DEBUG
   double etot1 = 0, etot2 = 0;
 #endif
 
   // Look for the Hit Collection of HCal
   G4HCofThisEvent* allHC = (*evt)()->GetHCofThisEvent();
   std::string sdName = names[0];
   idHC = G4SDManager::GetSDMpointer()->GetCollectionID(sdName);
   theHC = (CaloG4HitCollection*)allHC->GetHC(idHC);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Hit Collection for " << sdName << " of ID " << idHC
                                 << " is obtained at " << theHC << " with " << theHC->entries() << " entries";
 #endif
   int thehc_entries = theHC->entries();
   if (idHC >= 0 && theHC != nullptr) {
     hhits.reserve(theHC->entries());
     hhitl.reserve(theHC->entries());
     for (j = 0; j < thehc_entries; j++) {
       CaloG4Hit* aHit = (*theHC)[j];
       double e = aHit->getEnergyDeposit() / CLHEP::GeV;
       double time = aHit->getTimeSlice();
       math::XYZPoint pos = aHit->getEntry();
       unsigned int id = aHit->getUnitID();
       double theta = pos.theta();
       double eta = -std::log(std::tan(theta * 0.5));
       double phi = pos.phi();
       int det, z, group, ieta, iphi, layer;
       HcalTestNumbering::unpackHcalIndex(id, det, z, group, ieta, iphi, layer);
       double jitter = time - timeOfFlight(det, layer, eta);
       if (jitter < 0)
         jitter = 0;
       if (e < 0 || e > 1.)
         e = 0;
       double escl = e * scale(det, layer);
       unsigned int idx = HcalTBNumberingScheme::getUnitID(id, mode);
       CaloHit hit(det, layer, escl, eta, phi, jitter, idx);
       hhits.push_back(hit);
       CaloHit hitl(det, layer, escl, eta, phi, jitter, id);
       hhitl.push_back(hitl);
       primaries[aHit->getTrackID()] += e;
 #ifdef EDM_ML_DEBUG
       etot1 += escl;
       edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Hcal Hit i/p " << j << "  ID 0x" << std::hex << id << " 0x"
                                     << idx << std::dec << " time " << std::setw(6) << time << " " << std::setw(6)
                                     << jitter << " theta " << std::setw(8) << theta << " eta " << std::setw(8) << eta
                                     << " phi " << std::setw(8) << phi << " e " << std::setw(8) << e << " "
                                     << std::setw(8) << escl;
 #endif
     }
   }
 
   // Add hits in the same channel within same time slice
   std::vector<CaloHit>::iterator itr;
   int nHit = hhits.size();
   std::vector<CaloHit*> hits(nHit);
   for (j = 0, itr = hhits.begin(); itr != hhits.end(); j++, itr++) {
     hits[j] = &hhits[j];
   }
   sort(hits.begin(), hits.end(), CaloHitIdMore());
   std::vector<CaloHit*>::iterator k1, k2;
 #ifdef EDM_ML_DEBUG
   int nhit = 0;
 #endif
   for (k1 = hits.begin(); k1 != hits.end(); k1++) {
     int det = (**k1).det();
     int layer = (**k1).layer();
     double ehit = (**k1).e();
     double eta = (**k1).eta();
     double phi = (**k1).phi();
     double jitter = (**k1).t();
     uint32_t unitID = (**k1).id();
     int jump = 0;
     for (k2 = k1 + 1; k2 != hits.end() && std::fabs(jitter - (**k2).t()) < 1 && unitID == (**k2).id(); k2++) {
       ehit += (**k2).e();
       jump++;
     }
     CaloHit hit(det, layer, ehit, eta, phi, jitter, unitID);
     hcalHitCache.push_back(hit);
     k1 += jump;
 #ifdef EDM_ML_DEBUG
     nhit++;
     etot2 += ehit;
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Hcal Hit store " << nhit << "  ID 0x" << std::hex << unitID
                                   << std::dec << " time " << std::setw(6) << jitter << " eta " << std::setw(8) << eta
                                   << " phi " << std::setw(8) << phi << " e " << std::setw(8) << ehit;
 #endif
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Stores " << nhit << " HCal hits from " << nHit
                                 << " input hits E(Hcal) " << etot1 << " " << etot2;
 #endif
   //Repeat for Hit in each layer (hhits and hhitl sizes are the same)
   for (j = 0, itr = hhitl.begin(); itr != hhitl.end(); j++, itr++) {
     hits[j] = &hhitl[j];
   }
   sort(hits.begin(), hits.end(), CaloHitIdMore());
 #ifdef EDM_ML_DEBUG
   int nhitl = 0;
   double etotl = 0;
 #endif
   for (k1 = hits.begin(); k1 != hits.end(); k1++) {
     int det = (**k1).det();
     int layer = (**k1).layer();
     double ehit = (**k1).e();
     double eta = (**k1).eta();
     double phi = (**k1).phi();
     double jitter = (**k1).t();
     uint32_t unitID = (**k1).id();
     int jump = 0;
     for (k2 = k1 + 1; k2 != hits.end() && std::fabs(jitter - (**k2).t()) < 1 && unitID == (**k2).id(); k2++) {
       ehit += (**k2).e();
       jump++;
     }
     CaloHit hit(det, layer, ehit, eta, phi, jitter, unitID);
     hcalHitLayer.push_back(hit);
     k1 += jump;
 #ifdef EDM_ML_DEBUG
     nhitl++;
     etotl += ehit;
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Hcal Hit store " << nhitl << "  ID 0x" << std::hex << unitID
                                   << std::dec << " time " << std::setw(6) << jitter << " eta " << std::setw(8) << eta
                                   << " phi " << std::setw(8) << phi << " e " << std::setw(8) << ehit;
 #endif
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Stores " << nhitl << " HCal hits from " << nHit
                                 << " input hits E(Hcal) " << etot1 << " " << etotl;
 #endif
   // Look for the Hit Collection of ECal
   std::vector<CaloHit> ehits;
   sdName = names[1];
   idHC = G4SDManager::GetSDMpointer()->GetCollectionID(sdName);
   theHC = (CaloG4HitCollection*)allHC->GetHC(idHC);
 #ifdef EDM_ML_DEBUG
   etot1 = etot2 = 0;
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Hit Collection for " << sdName << " of ID " << idHC
                                 << " is obtained at " << theHC << " with " << theHC->entries() << " entries";
 #endif
   if (idHC >= 0 && theHC != nullptr) {
     thehc_entries = theHC->entries();
     ehits.reserve(theHC->entries());
     for (j = 0; j < thehc_entries; j++) {
       CaloG4Hit* aHit = (*theHC)[j];
       double e = aHit->getEnergyDeposit() / CLHEP::GeV;
       double time = aHit->getTimeSlice();
       if (e < 0 || e > 100000.)
         e = 0;
       if (e > 0) {
         math::XYZPoint pos = aHit->getEntry();
         unsigned int id = aHit->getUnitID();
         double theta = pos.theta();
         double eta = -std::log(std::tan(theta * 0.5));
         double phi = pos.phi();
         int det, z, group, ieta, iphi, layer;
         HcalTestNumbering::unpackHcalIndex(id, det, z, group, ieta, iphi, layer);
         CaloHit hit(det, 0, e, eta, phi, time, id);
         ehits.push_back(hit);
         primaries[aHit->getTrackID()] += e;
 #ifdef EDM_ML_DEBUG
         etot1 += e;
         edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Ecal Hit i/p " << j << "  ID 0x" << std::hex << id
                                       << std::dec << " time " << std::setw(6) << time << " theta " << std::setw(8)
                                       << theta << " eta " << std::setw(8) << eta << " phi " << std::setw(8) << phi
                                       << " e " << std::setw(8) << e;
 #endif
       }
     }
   }
 
   // Add hits in the same channel within same time slice
   nHit = ehits.size();
   std::vector<CaloHit*> hite(nHit);
   for (j = 0, itr = ehits.begin(); itr != ehits.end(); j++, itr++) {
     hite[j] = &ehits[j];
   }
   sort(hite.begin(), hite.end(), CaloHitIdMore());
 #ifdef EDM_ML_DEBUG
   nhit = 0;
 #endif
   for (k1 = hite.begin(); k1 != hite.end(); k1++) {
     int det = (**k1).det();
     int layer = (**k1).layer();
     double ehit = (**k1).e();
     double eta = (**k1).eta();
     double phi = (**k1).phi();
     double jitter = (**k1).t();
     uint32_t unitID = (**k1).id();
     int jump = 0;
     for (k2 = k1 + 1; k2 != hite.end() && std::fabs(jitter - (**k2).t()) < 1 && unitID == (**k2).id(); k2++) {
       ehit += (**k2).e();
       jump++;
     }
     CaloHit hit(det, layer, ehit, eta, phi, jitter, unitID);
     ecalHitCache.push_back(hit);
     k1 += jump;
 #ifdef EDM_ML_DEBUG
     nhit++;
     etot2 += ehit;
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Ecal Hit store " << nhit << "  ID 0x" << std::hex << unitID
                                   << std::dec << " time " << std::setw(6) << jitter << " eta " << std::setw(8) << eta
                                   << " phi " << std::setw(8) << phi << " e " << std::setw(8) << ehit;
 #endif
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Stores " << nhit << " ECal hits from " << nHit
                                 << " input hits E(Ecal) " << etot1 << " " << etot2;
 #endif
   // Find Primary info:
   nPrimary = static_cast<int>(primaries.size());
   int trackID = 0;
   G4PrimaryParticle* thePrim = nullptr;
   int nvertex = (*evt)()->GetNumberOfPrimaryVertex();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Event has " << nvertex << " verteices";
 #endif
   if (nvertex <= 0)
     edm::LogWarning("HcalTBSim") << "HcalTB04Analysis::EndOfEvent ERROR: no vertex found for event " << evNum;
   for (int i = 0; i < nvertex; i++) {
     G4PrimaryVertex* avertex = (*evt)()->GetPrimaryVertex(i);
     if (avertex == nullptr) {
       edm::LogWarning("HcalTBSim") << "HcalTB04Analysis::EndOfEvent ERR: pointer to vertex = 0 for event " << evNum;
     } else {
       edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::Vertex number :" << i << " " << avertex->GetPosition();
       int npart = avertex->GetNumberOfParticle();
       if (npart == 0)
         edm::LogWarning("HcalTBSim") << "HcalTB04Analysis::End Of Event ERR: no primary!";
       if (thePrim == nullptr)
         thePrim = avertex->GetPrimary(trackID);
     }
   }
 
   if (thePrim != nullptr) {
     double px = thePrim->GetPx();
     double py = thePrim->GetPy();
     double pz = thePrim->GetPz();
     double p = std::sqrt(pow(px, 2.) + pow(py, 2.) + pow(pz, 2.));
     pInit = p / CLHEP::GeV;
     if (p == 0)
       edm::LogWarning("HcalTBSim") << "HcalTB04Analysis:: EndOfEvent ERR: primary has p=0 ";
     else {
       double costheta = pz / p;
       double theta = acos(std::min(std::max(costheta, -1.), 1.));
       etaInit = -std::log(std::tan(theta / 2));
       if (px != 0 || py != 0)
         phiInit = std::atan2(py, px);
     }
     particleType = thePrim->GetPDGcode();
   } else
     edm::LogWarning("HcalTBSim") << "HcalTB04Analysis::EndOfEvent ERR: could not find primary";
 }

◆ fillEvent()

void HcalTB04Analysis::fillEvent ( PHcalTB04Info & product )

private

Definition at line 958 of file HcalTB04Analysis.cc.

References TauDecayModes::dec, ecalHitCache, eecalq, eecals, ehcalq, ehcals, enois, eqeta, eqie, eqlay, eqphi, eseta, esime, esimh, eslay, esphi, etaInit, etotq, etots, evNum, watchdog::group, hcalHitCache, mps_fire::i, idHcal, idXtal, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, nPrimary, particleType, phiInit, pInit, pvMomentum, pvPosition, pvType, pvUVW, PHcalTB04Info::saveHit(), secEkin, secMomentum, secPartID, secTrackID, PHcalTB04Info::setEdep(), PHcalTB04Info::setEdepHcal(), PHcalTB04Info::setIDs(), PHcalTB04Info::setLongProf(), PHcalTB04Info::setPrimary(), PHcalTB04Info::setTrnsProf(), PHcalTB04Info::setVtxPrim(), PHcalTB04Info::setVtxSec(), HcalTestNumbering::unpackHcalIndex(), x, y, and z.

Referenced by produce().

                                                        {
   //Setup the ID's
   product.setIDs(idHcal, idXtal);
 
   //Beam Information
   product.setPrimary(nPrimary, particleType, pInit, etaInit, phiInit);
 
   //Energy deposits in the crystals and towers
   product.setEdepHcal(esimh, eqie);
   product.setEdepHcal(esime, enois);
 
   // Total Energy
   product.setEdep(etots, eecals, ehcals, etotq, eecalq, ehcalq);
 
   // Lateral Profile
   product.setTrnsProf(eseta, eqeta, esphi, eqphi);
 
   // Longitudianl profile
   product.setLongProf(eslay, eqlay);
 
   //Save Hits
   int nhit = 0;
   std::vector<CaloHit>::iterator itr;
   for (itr = ecalHitCache.begin(); itr != ecalHitCache.end(); itr++) {
     uint32_t id = itr->id();
     int det, z, group, ieta, iphi, lay;
     HcalTestNumbering::unpackHcalIndex(id, det, z, group, ieta, iphi, lay);
     product.saveHit(det, lay, ieta, iphi, itr->e(), itr->t());
     nhit++;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Save Hit " << std::setw(3) << nhit << " ID 0x" << std::hex
                                   << group << std::dec << " " << std::setw(2) << det << " " << std::setw(2) << lay
                                   << " " << std::setw(1) << z << " " << std::setw(3) << ieta << " " << std::setw(3)
                                   << iphi << " T " << std::setw(6) << itr->t() << " E " << std::setw(6) << itr->e();
 #endif
   }
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Saves " << nhit << " hits from Crystals";
 #ifdef EDM_ML_DEBUG
   int nhit0 = nhit;
 #endif
   nhit = 0;
   for (itr = hcalHitCache.begin(); itr != hcalHitCache.end(); itr++) {
     uint32_t id = itr->id();
     int det, z, group, ieta, iphi, lay;
     HcalTestNumbering::unpackHcalIndex(id, det, z, group, ieta, iphi, lay);
     product.saveHit(det, lay, ieta, iphi, itr->e(), itr->t());
     nhit++;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Save Hit " << std::setw(3) << nhit + nhit0 << " ID 0x"
                                   << std::hex << group << std::dec << " " << std::setw(2) << det << " " << std::setw(2)
                                   << lay << " " << std::setw(1) << z << " " << std::setw(3) << ieta << " "
                                   << std::setw(3) << iphi << " T " << std::setw(6) << itr->t() << " E " << std::setw(6)
                                   << itr->e();
 #endif
   }
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Saves " << nhit << " hits from HCal";
 
   //Vertex associated quantities
   product.setVtxPrim(evNum,
                      pvType,
                      pvPosition.x(),
                      pvPosition.y(),
                      pvPosition.z(),
                      pvUVW.x(),
                      pvUVW.y(),
                      pvUVW.z(),
                      pvMomentum.x(),
                      pvMomentum.y(),
                      pvMomentum.z());
   for (unsigned int i = 0; i < secTrackID.size(); i++) {
     product.setVtxSec(
         secTrackID[i], secPartID[i], secMomentum[i].x(), secMomentum[i].y(), secMomentum[i].z(), secEkin[i]);
   }
 }

◆ finalAnalysis()

void HcalTB04Analysis::finalAnalysis ( )

private

Definition at line 847 of file HcalTB04Analysis.cc.

References StorageManager_cfg::e1, eecalq, eecals, ehcalq, ehcals, enois, eqeta, eqie, eqlay, eqphi, eseta, esime, esimh, eslay, esphi, etaInit, etotq, etots, HcalTB04Histo::fillEdep(), HcalTB04Histo::fillLongProf(), HcalTB04Histo::fillPrimary(), HcalTB04Histo::fillTrnsProf(), watchdog::group, histo, mps_fire::i, iceta, icphi, l1ctLayer2EG_cff::id, idTower, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, nano_mu_digi_cff::layer, nCrystal, nTower, phiInit, pInit, HcalTestNumbering::unpackHcalIndex(), and z.

Referenced by update().

                                      {
   //Beam Information
   histo->fillPrimary(pInit, etaInit, phiInit);
 
   // Total Energy
   eecals = ehcals = eecalq = ehcalq = 0.;
   for (int i = 0; i < nTower; i++) {
     ehcals += esimh[i];
     ehcalq += eqie[i];
   }
   for (int i = 0; i < nCrystal; i++) {
     eecals += esime[i];
     eecalq += enois[i];
   }
   etots = eecals + ehcals;
   etotq = eecalq + ehcalq;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Energy deposit at Sim Level (Total) " << etots << " (ECal) "
                                 << eecals << " (HCal) " << ehcals
                                 << "\nHcalTB04Analysis:: Energy deposit at Qie Level (Total) " << etotq << " (ECal) "
                                 << eecalq << " (HCal) " << ehcalq;
 #endif
   histo->fillEdep(etots, eecals, ehcals, etotq, eecalq, ehcalq);
 
   // Lateral Profile
   for (int i = 0; i < 5; i++) {
     eseta[i] = 0.;
     eqeta[i] = 0.;
   }
   for (int i = 0; i < 3; i++) {
     esphi[i] = 0.;
     eqphi[i] = 0.;
   }
   double e1 = 0, e2 = 0;
   unsigned int id;
   for (int i = 0; i < nTower; i++) {
     int det, z, group, ieta, iphi, layer;
     id = idTower[i];
     HcalTestNumbering::unpackHcalIndex(id, det, z, group, ieta, iphi, layer);
     iphi -= (icphi - 1);
     if (icphi > 4) {
       if (ieta == 0)
         ieta = 2;
       else
         ieta = -1;
     } else {
       ieta = ieta - iceta + 2;
     }
     if (iphi >= 0 && iphi < 3 && ieta >= 0 && ieta < 5) {
       eseta[ieta] += esimh[i];
       esphi[iphi] += esimh[i];
       e1 += esimh[i];
       eqeta[ieta] += eqie[i];
       eqphi[iphi] += eqie[i];
       e2 += eqie[i];
     }
   }
   for (int i = 0; i < 3; i++) {
     if (e1 > 0)
       esphi[i] /= e1;
     if (e2 > 0)
       eqphi[i] /= e2;
   }
   for (int i = 0; i < 5; i++) {
     if (e1 > 0)
       eseta[i] /= e1;
     if (e2 > 0)
       eqeta[i] /= e2;
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Energy fraction along Eta and Phi (Sim/Qie)";
   for (int i = 0; i < 5; i++)
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: [" << i << "] Eta Sim = " << eseta[i] << " Qie = " << eqeta[i]
                                   << " Phi Sim = " << esphi[i] << " Qie = " << eqphi[i];
 #endif
   histo->fillTrnsProf(eseta, eqeta, esphi, eqphi);
 
   // Longitudianl profile
   for (int i = 0; i < 20; i++) {
     eslay[i] = 0.;
     eqlay[i] = 0.;
   }
   e1 = 0;
   e2 = 0;
   for (int i = 0; i < nTower; i++) {
     int det, z, group, ieta, iphi, layer;
     id = idTower[i];
     HcalTestNumbering::unpackHcalIndex(id, det, z, group, ieta, iphi, layer);
     iphi -= (icphi - 1);
     layer -= 1;
     if (iphi >= 0 && iphi < 3 && layer >= 0 && layer < 20) {
       eslay[layer] += esimh[i];
       e1 += esimh[i];
       eqlay[layer] += eqie[i];
       e2 += eqie[i];
     }
   }
   for (int i = 0; i < 20; i++) {
     if (e1 > 0)
       eslay[i] /= e1;
     if (e2 > 0)
       eqlay[i] /= e2;
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Energy fraction along Layer";
   for (int i = 0; i < 20; i++)
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: [" << i << "] Sim = " << eslay[i] << " Qie = " << eqlay[i];
 #endif
   histo->fillLongProf(eslay, eqlay);
 }

◆ init()

void HcalTB04Analysis::init ( void )

private

Definition at line 222 of file HcalTB04Analysis.cc.

References clear(), count, TauDecayModes::dec, eqeta, eqlay, eqphi, eseta, eslay, esphi, evNum, HcalTBNumberingScheme::getUnitIDs(), hcalOnly, mps_fire::i, globals_cff::id1, idEcal, idHcal, idTower, idXtal, mode, nCrystal, nTower, HcalTestNumbering::packHcalIndex(), and unitID().

Referenced by HcalTB04Analysis().

                             {
   idTower = HcalTBNumberingScheme::getUnitIDs(type, mode);
   nTower = idTower.size();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Save information from " << nTower << " HCal towers";
 #endif
   idHcal.reserve(nTower);
   for (int i = 0; i < nTower; i++) {
     int id = unitID(idTower[i]);
     idHcal.push_back(id);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalTBSim") << "\tTower[" << i << "] Original " << std::hex << idTower[i] << " Stored "
                                   << idHcal[i] << std::dec;
 #endif
   }
 
   if (!hcalOnly) {
     int det = 10;
     uint32_t id1;
     nCrystal = 0;
     for (int lay = 1; lay < 8; lay++) {
       for (int icr = 1; icr < 8; icr++) {
         id1 = HcalTestNumbering::packHcalIndex(det, 0, 1, icr, lay, 1);
         int id = unitID(id1);
         idEcal.push_back(id1);
         idXtal.push_back(id);
         nCrystal++;
       }
     }
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Save information from " << nCrystal << " ECal Crystals";
 #ifdef EDM_ML_DEBUG
     for (int i = 0; i < nCrystal; i++) {
       edm::LogVerbatim("HcalTBSim") << "\tCrystal[" << i << "] Original " << std::hex << idEcal[i] << " Stored "
                                     << idXtal[i] << std::dec;
     }
 #endif
   }
   // Profile vectors
   eseta.reserve(5);
   eqeta.reserve(5);
   esphi.reserve(3);
   eqphi.reserve(3);
   eslay.reserve(20);
   eqlay.reserve(20);
   for (int i = 0; i < 5; i++) {
     eseta.push_back(0.);
     eqeta.push_back(0.);
   }
   for (int i = 0; i < 3; i++) {
     esphi.push_back(0.);
     eqphi.push_back(0.);
   }
   for (int i = 0; i < 20; i++) {
     eslay.push_back(0.);
     eqlay.push_back(0.);
   }
 
   // counter
   count = 0;
   evNum = 0;
   clear();
 }

◆ operator=()

const HcalTB04Analysis& HcalTB04Analysis::operator= ( const HcalTB04Analysis & )

delete

◆ produce()

void HcalTB04Analysis::produce	(	edm::Event &	e,
		const edm::EventSetup &
	)

overridevirtual

Implements SimProducer.

Definition at line 216 of file HcalTB04Analysis.cc.

References MillePedeFileConverter_cfg::e, fillEvent(), and eostools::move().

                                                                 {
   std::unique_ptr<PHcalTB04Info> product(new PHcalTB04Info);
   fillEvent(*product);
   e.put(std::move(product));
 }

◆ qieAnalysis()

void HcalTB04Analysis::qieAnalysis ( CLHEP::HepRandomEngine * engine )

private

Definition at line 737 of file HcalTB04Analysis.cc.

References hippyaddtobaddatafiles::cd(), TauDecayModes::dec, eqie, esimh, HcalQie::getCode(), HcalQie::getEnergy(), hcalHitCache, hfClusterShapes_cfi::hits, l1ctLayer2EG_cff::id, hit::id, idTower, myQie, nTower, and DQMOfflineHeavyIons_cff::todo.

Referenced by update().

                                                                {
   int hittot = hcalHitCache.size();
   if (hittot <= 0)
     hittot = 1;
   std::vector<CaloHit> hits(hittot);
   std::vector<int> todo(nTower, 0);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::qieAnalysis: Size " << hits.size() << " " << todo.size() << " "
                                 << idTower.size() << " " << esimh.size() << " " << eqie.size();
 #endif
   // Loop over all HCal hits
   for (unsigned int k1 = 0; k1 < hcalHitCache.size(); k1++) {
     CaloHit hit = hcalHitCache[k1];
     uint32_t id = hit.id();
     int nhit = 0;
     double esim = hit.e();
     hits[nhit] = hit;
     for (unsigned int k2 = k1 + 1; k2 < hcalHitCache.size(); k2++) {
       hit = hcalHitCache[k2];
       if (hit.id() == id) {
         nhit++;
         hits[nhit] = hit;
         esim += hit.e();
       }
     }
     k1 += nhit;
     nhit++;
     std::vector<int> cd = myQie->getCode(nhit, hits, engine);
     double eq = myQie->getEnergy(cd);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::  ID 0x" << std::hex << id << std::dec << " registers " << esim
                                   << " energy from " << nhit << " hits starting with hit # " << k1
                                   << " energy with noise " << eq;
 #endif
     for (int k2 = 0; k2 < nTower; k2++) {
       if (id == idTower[k2]) {
         todo[k2] = 1;
         esimh[k2] = esim;
         eqie[k2] = eq;
       }
     }
   }
 
   // Towers with no hit
   for (int k2 = 0; k2 < nTower; k2++) {
     if (todo[k2] == 0) {
       std::vector<int> cd = myQie->getCode(0, hits, engine);
       double eq = myQie->getEnergy(cd);
       esimh[k2] = 0;
       eqie[k2] = eq;
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::  ID 0x" << std::hex << idTower[k2] << std::dec
                                     << " registers " << esimh[k2] << " energy from hits and energy after QIE analysis "
                                     << eqie[k2];
 #endif
     }
   }
 }

◆ scale()

double HcalTB04Analysis::scale	(	int	det,
		int	layer
	)

private

Definition at line 1080 of file HcalTB04Analysis.cc.

References HcalBarrel, nano_mu_digi_cff::layer, scaleHB0, scaleHB16, scaleHE0, scaleHO, and createJobs::tmp.

Referenced by fillBuffer().

                                                  {
   double tmp = 1.;
   if (det == static_cast<int>(HcalBarrel)) {
     if (layer == 1)
       tmp = scaleHB0;
     else if (layer == 17)
       tmp = scaleHB16;
     else if (layer > 17)
       tmp = scaleHO;
   } else {
     if (layer <= 2)
       tmp = scaleHE0;
   }
   return tmp;
 }

◆ timeOfFlight()

double HcalTB04Analysis::timeOfFlight	(	int	det,
		int	layer,
		double	eta
	)

private

Definition at line 1096 of file HcalTB04Analysis.cc.

References beamOffset, funct::cos(), PVValHelper::eta, JetChargeProducer_cfi::exp, HcalBarrel, nano_mu_digi_cff::layer, funct::sin(), theta(), and createJobs::tmp.

Referenced by fillBuffer().

                                                                     {
   double theta = 2.0 * std::atan(std::exp(-eta));
   double dist = beamOffset;
   if (det == static_cast<int>(HcalBarrel)) {
     const double rLay[19] = {1836.0,
                              1902.0,
                              1962.0,
                              2022.0,
                              2082.0,
                              2142.0,
                              2202.0,
                              2262.0,
                              2322.0,
                              2382.0,
                              2448.0,
                              2514.0,
                              2580.0,
                              2646.0,
                              2712.0,
                              2776.0,
                              2862.5,
                              3847.0,
                              4052.0};
     if (layer > 0 && layer <= 19)
       dist += rLay[layer - 1] * mm / sin(theta);
   } else {
     const double zLay[19] = {4034.0,
                              4032.0,
                              4123.0,
                              4210.0,
                              4297.0,
                              4384.0,
                              4471.0,
                              4558.0,
                              4645.0,
                              4732.0,
                              4819.0,
                              4906.0,
                              4993.0,
                              5080.0,
                              5167.0,
                              5254.0,
                              5341.0,
                              5428.0,
                              5515.0};
     if (layer > 0 && layer <= 19)
       dist += zLay[layer - 1] * mm / cos(theta);
   }
 
   double tmp = dist / c_light / ns;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::timeOfFlight " << tmp << " for det/lay " << det << " " << layer
                                 << " eta/theta " << eta << " " << theta / deg << " dist " << dist;
 #endif
   return tmp;
 }

◆ unitID()

int HcalTB04Analysis::unitID ( uint32_t id )

private

Definition at line 1069 of file HcalTB04Analysis.cc.

References watchdog::group, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, HcalTestNumbering::unpackHcalIndex(), and z.

Referenced by fillBuffer(), and init().

                                         {
   int det, z, group, ieta, iphi, lay;
   HcalTestNumbering::unpackHcalIndex(id, det, z, group, ieta, iphi, lay);
   group = (det & 15) << 20;
   group += ((lay - 1) & 31) << 15;
   group += (z & 1) << 14;
   group += (ieta & 127) << 7;
   group += (iphi & 127);
   return group;
 }

◆ update() [1/4]

void HcalTB04Analysis::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 285 of file HcalTB04Analysis.cc.

References hcalOnly, names, ECalSD::setNumberingScheme(), HCalSD::setNumberingScheme(), and AlCaHLTBitMon_QueryRunRegistry::string.

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

                                                    {
   int irun = (*run)()->GetRunID();
   edm::LogVerbatim("HcalTBSim") << " =====> Begin of Run = " << irun;
 
   G4SDManager* sd = G4SDManager::GetSDMpointerIfExist();
   if (sd != nullptr) {
     std::string sdname = names[0];
     G4VSensitiveDetector* aSD = sd->FindSensitiveDetector(sdname);
     if (aSD == nullptr) {
       edm::LogWarning("HcalTBSim") << "HcalTB04Analysis::beginOfRun: No SD"
                                    << " with name " << sdname << " in this "
                                    << "Setup";
     } else {
       HCalSD* theCaloSD = dynamic_cast<HCalSD*>(aSD);
       edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::beginOfRun: Finds SD with name " << theCaloSD->GetName()
                                     << " in this Setup";
       HcalNumberingScheme* org = new HcalTestNumberingScheme(false);
       theCaloSD->setNumberingScheme(org);
       edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::beginOfRun: set a new numbering scheme";
     }
     if (!hcalOnly) {
       sdname = names[1];
       aSD = sd->FindSensitiveDetector(sdname);
       if (aSD == nullptr) {
         edm::LogWarning("HcalTBSim") << "HcalTB04Analysis::beginOfRun: No SD"
                                      << " with name " << sdname << " in this "
                                      << "Setup";
       } else {
         ECalSD* theCaloSD = dynamic_cast<ECalSD*>(aSD);
         edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::beginOfRun: Finds SD with name " << theCaloSD->GetName()
                                       << " in this Setup";
         EcalNumberingScheme* org = new HcalTB04XtalNumberingScheme();
         theCaloSD->setNumberingScheme(org);
         edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::beginOfRun: set a new numbering scheme";
       }
     }
   } else {
     edm::LogWarning("HcalTBSim") << "HcalTB04Analysis::beginOfRun: Could "
                                  << "not get SD Manager!";
   }
 }

◆ update() [2/4]

void HcalTB04Analysis::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 327 of file HcalTB04Analysis.cc.

References clear(), and evNum.

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

                                                      {
   clear();
 #ifdef EDM_ML_DEBUG
   evNum = (*evt)()->GetEventID();
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis: =====> Begin of event = " << evNum;
 #endif
 }

◆ update() [3/4]

void HcalTB04Analysis::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 335 of file HcalTB04Analysis.cc.

References position, pvFound, pvMomentum, pvPosition, pvType, pvUVW, secEkin, secMomentum, secPartID, secTrackID, and shortLivedSecondaries.

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

                                                  {
   if (aStep != nullptr) {
     //Get Step properties
     G4ThreeVector thePreStepPoint = aStep->GetPreStepPoint()->GetPosition();
     G4ThreeVector thePostStepPoint;
 
     // Get Tracks properties
     G4Track* aTrack = aStep->GetTrack();
     int trackID = aTrack->GetTrackID();
     int parentID = aTrack->GetParentID();
     const G4ThreeVector& position = aTrack->GetPosition();
     G4ThreeVector momentum = aTrack->GetMomentum();
     G4String partType = aTrack->GetDefinition()->GetParticleType();
     G4String partSubType = aTrack->GetDefinition()->GetParticleSubType();
     int partPDGEncoding = aTrack->GetDefinition()->GetPDGEncoding();
 #ifdef EDM_ML_DEBUG
     bool isPDGStable = aTrack->GetDefinition()->GetPDGStable();
 #endif
     double pDGlifetime = aTrack->GetDefinition()->GetPDGLifeTime();
     double gammaFactor = aStep->GetPreStepPoint()->GetGamma();
 
     if (!pvFound) {  //search for v1
       double stepDeltaEnergy = aStep->GetDeltaEnergy();
       double kinEnergy = aTrack->GetKineticEnergy();
 
       // look for DeltaE > 10% kinEnergy of particle, or particle death - Ek=0
       if (trackID == 1 && parentID == 0 && ((kinEnergy == 0.) || (std::fabs(stepDeltaEnergy / kinEnergy) > 0.1))) {
         pvType = -1;
         if (kinEnergy == 0.) {
           pvType = 0;
         } else {
           if (std::fabs(stepDeltaEnergy / kinEnergy) > 0.1)
             pvType = 1;
         }
         pvFound = true;
         pvPosition = position;
         pvMomentum = momentum;
         // Rotated coord.system:
         pvUVW = (*beamline_RM) * (pvPosition);
 
         //Volume name requires some checks:
         G4String thePostPVname = "NoName";
         G4StepPoint* thePostPoint = aStep->GetPostStepPoint();
         if (thePostPoint) {
           thePostStepPoint = thePostPoint->GetPosition();
           G4VPhysicalVolume* thePostPV = thePostPoint->GetPhysicalVolume();
           if (thePostPV)
             thePostPVname = thePostPV->GetName();
         }
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: V1 found at: " << thePostStepPoint
                                       << " G4VPhysicalVolume: " << thePostPVname;
         edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::fill_v1Pos: Primary Track momentum: " << pvMomentum
                                       << " psoition " << pvPosition << " u/v/w " << pvUVW;
 #endif
       }
     } else {
       // watch for secondaries originating @v1, including the surviving primary
       if ((trackID != 1 && parentID == 1 && (aTrack->GetCurrentStepNumber() == 1) && (thePreStepPoint == pvPosition)) ||
           (trackID == 1 && thePreStepPoint == pvPosition)) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::A secondary...  PDG:" << partPDGEncoding
                                       << " TrackID:" << trackID << " ParentID:" << parentID
                                       << " stable: " << isPDGStable << " Tau: " << pDGlifetime
                                       << " cTauGamma=" << c_light * pDGlifetime * gammaFactor * 1000.
                                       << "um GammaFactor: " << gammaFactor;
 #endif
         secTrackID.push_back(trackID);
         secPartID.push_back(partPDGEncoding);
         secMomentum.push_back(momentum);
         secEkin.push_back(aTrack->GetKineticEnergy());
 
         // Check for short-lived secondaries: cTauGamma<100um
         double ctaugamma_um = CLHEP::c_light * pDGlifetime * gammaFactor * 1000.;
         if ((ctaugamma_um > 0.) && (ctaugamma_um < 100.)) {  //short-lived secondary
           shortLivedSecondaries.push_back(trackID);
         } else {  //normal secondary - enter into the V1-calorimetric tree
                   //          histos->fill_v1cSec (aTrack);
         }
       }
       // Also watch for tertiary particles coming from
       // short-lived secondaries from V1
       if (aTrack->GetCurrentStepNumber() == 1) {
         if (!shortLivedSecondaries.empty()) {
           int pid = parentID;
           std::vector<int>::iterator pos1 = shortLivedSecondaries.begin();
           std::vector<int>::iterator pos2 = shortLivedSecondaries.end();
           std::vector<int>::iterator pos;
           for (pos = pos1; pos != pos2; pos++) {
             if (*pos == pid) {  //ParentID is on the list of short-lived
                                 // secondary
 #ifdef EDM_ML_DEBUG
               edm::LogVerbatim("HcalTBSim")
                   << "HcalTB04Analysis:: A tertiary...  PDG:" << partPDGEncoding << " TrackID:" << trackID
                   << " ParentID:" << parentID << " stable: " << isPDGStable << " Tau: " << pDGlifetime
                   << " cTauGamma=" << c_light * pDGlifetime * gammaFactor * 1000. << "um GammaFactor: " << gammaFactor;
 #endif
             }
           }
         }
       }
     }
   }
 }

◆ update() [4/4]

void HcalTB04Analysis::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent *>.

Definition at line 440 of file HcalTB04Analysis.cc.

References count, ecalHitCache, fillBuffer(), finalAnalysis(), hcalHitCache, hcalOnly, qieAnalysis(), and xtalAnalysis().

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

                                                    {
   count++;
 
   //fill the buffer
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::Fill event " << (*evt)()->GetEventID();
 #endif
   fillBuffer(evt);
 
   //QIE analysis
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::Do QIE analysis with " << hcalHitCache.size() << " hits";
 #endif
   CLHEP::HepRandomEngine* engine = G4Random::getTheEngine();
   qieAnalysis(engine);
 
   //Energy in Crystal Matrix
   if (!hcalOnly) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::Do Xtal analysis with " << ecalHitCache.size() << " hits";
 #endif
     xtalAnalysis(engine);
   }
 
   //Final Analysis
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::Final analysis";
 #endif
   finalAnalysis();
 
   int iEvt = (*evt)()->GetEventID();
   if (iEvt < 10)
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Event " << iEvt;
   else if ((iEvt < 100) && (iEvt % 10 == 0))
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Event " << iEvt;
   else if ((iEvt < 1000) && (iEvt % 100 == 0))
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Event " << iEvt;
   else if ((iEvt < 10000) && (iEvt % 1000 == 0))
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis:: Event " << iEvt;
 }

◆ xtalAnalysis()

void HcalTB04Analysis::xtalAnalysis ( CLHEP::HepRandomEngine * engine )

private

Definition at line 797 of file HcalTB04Analysis.cc.

References TauDecayModes::dec, ecalHitCache, ecalNoise, enois, esime, l1ctLayer2EG_cff::id, idEcal, and nCrystal.

Referenced by update().

                                                                 {
   CLHEP::RandGaussQ randGauss(*engine);
 
   // Crystal Data
   std::vector<int> iok(nCrystal, 0);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::xtalAnalysis: Size " << iok.size() << " " << idEcal.size() << " "
                                 << esime.size() << " " << enois.size();
 #endif
   for (unsigned int k1 = 0; k1 < ecalHitCache.size(); k1++) {
     uint32_t id = ecalHitCache[k1].id();
     int nhit = 0;
     double esim = ecalHitCache[k1].e();
     for (unsigned int k2 = k1 + 1; k2 < ecalHitCache.size(); k2++) {
       if (ecalHitCache[k2].id() == id) {
         nhit++;
         esim += ecalHitCache[k2].e();
       }
     }
     k1 += nhit;
     nhit++;
     double eq = esim + randGauss.fire(0., ecalNoise);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::  ID 0x" << std::hex << id << std::dec << " registers " << esim
                                   << " energy from " << nhit << " hits starting with hit # " << k1
                                   << " energy with noise " << eq;
 #endif
     for (int k2 = 0; k2 < nCrystal; k2++) {
       if (id == idEcal[k2]) {
         iok[k2] = 1;
         esime[k2] = esim;
         enois[k2] = eq;
       }
     }
   }
 
   // Crystals with no hit
   for (int k2 = 0; k2 < nCrystal; k2++) {
     if (iok[k2] == 0) {
       esime[k2] = 0;
       enois[k2] = randGauss.fire(0., ecalNoise);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalTBSim") << "HcalTB04Analysis::  ID 0x" << std::hex << idEcal[k2] << std::dec
                                     << " registers " << esime[k2] << " energy from hits and energy from noise "
                                     << enois[k2];
 #endif
     }
   }
 }

Member Data Documentation

◆ beamline_RM

G4RotationMatrix* HcalTB04Analysis::beamline_RM

private

Definition at line 120 of file HcalTB04Analysis.cc.

Referenced by HcalTB04Analysis().

◆ beamOffset

const double HcalTB04Analysis::beamOffset

private

Definition at line 112 of file HcalTB04Analysis.cc.

Referenced by HcalTB04Analysis(), and timeOfFlight().

◆ count

int HcalTB04Analysis::count

private

Definition at line 123 of file HcalTB04Analysis.cc.

Referenced by init(), update(), and ~HcalTB04Analysis().

◆ ecalHitCache

std::vector<CaloHit> HcalTB04Analysis::ecalHitCache

private

Definition at line 131 of file HcalTB04Analysis.cc.

Referenced by clear(), fillBuffer(), fillEvent(), update(), and xtalAnalysis().

◆ ecalNoise

const double HcalTB04Analysis::ecalNoise

private

Definition at line 112 of file HcalTB04Analysis.cc.

Referenced by HcalTB04Analysis(), and xtalAnalysis().

◆ eecalq

double HcalTB04Analysis::eecalq

private

Definition at line 135 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), and finalAnalysis().

◆ eecals

double HcalTB04Analysis::eecals

private

Definition at line 135 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), and finalAnalysis().

◆ ehcalq

double HcalTB04Analysis::ehcalq

private

Definition at line 135 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), and finalAnalysis().

◆ ehcals

double HcalTB04Analysis::ehcals

private

Definition at line 135 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), and finalAnalysis().

◆ enois

std::vector<double> HcalTB04Analysis::enois

private

Definition at line 133 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), finalAnalysis(), and xtalAnalysis().

◆ eqeta

std::vector<double> HcalTB04Analysis::eqeta

private

Definition at line 134 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), finalAnalysis(), and init().

◆ eqie

std::vector<double> HcalTB04Analysis::eqie

private

Definition at line 133 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), finalAnalysis(), and qieAnalysis().

◆ eqlay

std::vector<double> HcalTB04Analysis::eqlay

private

Definition at line 134 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), finalAnalysis(), and init().

◆ eqphi

std::vector<double> HcalTB04Analysis::eqphi

private

Definition at line 134 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), finalAnalysis(), and init().

◆ eseta

std::vector<double> HcalTB04Analysis::eseta

private

Definition at line 134 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), finalAnalysis(), and init().

◆ esime

std::vector<double> HcalTB04Analysis::esime

private

Definition at line 133 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), finalAnalysis(), and xtalAnalysis().

◆ esimh

std::vector<double> HcalTB04Analysis::esimh

private

Definition at line 133 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), finalAnalysis(), and qieAnalysis().

◆ eslay

std::vector<double> HcalTB04Analysis::eslay

private

Definition at line 134 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), finalAnalysis(), and init().

◆ esphi

std::vector<double> HcalTB04Analysis::esphi

private

Definition at line 134 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), finalAnalysis(), and init().

◆ etaInit

double HcalTB04Analysis::etaInit

private

Definition at line 130 of file HcalTB04Analysis.cc.

Referenced by clear(), fillBuffer(), fillEvent(), and finalAnalysis().

◆ etotq

double HcalTB04Analysis::etotq

private

Definition at line 135 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), and finalAnalysis().

◆ etots

double HcalTB04Analysis::etots

private

Definition at line 135 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), and finalAnalysis().

◆ evNum

int HcalTB04Analysis::evNum

private

Definition at line 138 of file HcalTB04Analysis.cc.

Referenced by fillBuffer(), fillEvent(), init(), and update().

◆ hcalHitCache

std::vector<CaloHit> HcalTB04Analysis::hcalHitCache

private

Definition at line 132 of file HcalTB04Analysis.cc.

Referenced by clear(), fillBuffer(), fillEvent(), qieAnalysis(), and update().

◆ hcalHitLayer

std::vector<CaloHit> HcalTB04Analysis::hcalHitLayer

private

Definition at line 132 of file HcalTB04Analysis.cc.

Referenced by clear(), and fillBuffer().

◆ hcalOnly

const bool HcalTB04Analysis::hcalOnly

private

Definition at line 110 of file HcalTB04Analysis.cc.

Referenced by HcalTB04Analysis(), init(), and update().

◆ histo

HcalTB04Histo* HcalTB04Analysis::histo

private

Definition at line 117 of file HcalTB04Analysis.cc.

Referenced by finalAnalysis(), HcalTB04Analysis(), and ~HcalTB04Analysis().

◆ iceta

int HcalTB04Analysis::iceta

private

Definition at line 119 of file HcalTB04Analysis.cc.

Referenced by finalAnalysis(), and HcalTB04Analysis().

◆ icphi

int HcalTB04Analysis::icphi

private

Definition at line 119 of file HcalTB04Analysis.cc.

Referenced by finalAnalysis(), and HcalTB04Analysis().

◆ idEcal

std::vector<uint32_t> HcalTB04Analysis::idEcal

private

Definition at line 126 of file HcalTB04Analysis.cc.

Referenced by init(), and xtalAnalysis().

◆ idHcal

std::vector<int> HcalTB04Analysis::idHcal

private

Definition at line 125 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), and init().

◆ idTower

std::vector<uint32_t> HcalTB04Analysis::idTower

private

Definition at line 126 of file HcalTB04Analysis.cc.

Referenced by finalAnalysis(), init(), and qieAnalysis().

◆ idXtal

std::vector<int> HcalTB04Analysis::idXtal

private

Definition at line 125 of file HcalTB04Analysis.cc.

Referenced by fillEvent(), and init().

◆ m_Anal

const edm::ParameterSet HcalTB04Analysis::m_Anal

private

Definition at line 109 of file HcalTB04Analysis.cc.

Referenced by HcalTB04Analysis().

◆ mode

const int HcalTB04Analysis::mode

private

Definition at line 111 of file HcalTB04Analysis.cc.

Referenced by fillBuffer(), HcalTB04Analysis(), and init().

◆ myQie

HcalQie* HcalTB04Analysis::myQie

private

Definition at line 116 of file HcalTB04Analysis.cc.

Referenced by HcalTB04Analysis(), qieAnalysis(), and ~HcalTB04Analysis().

◆ names

const std::vector<std::string> HcalTB04Analysis::names

private

Definition at line 114 of file HcalTB04Analysis.cc.

Referenced by fillBuffer(), and update().

◆ nCrystal

int HcalTB04Analysis::nCrystal

private

Definition at line 124 of file HcalTB04Analysis.cc.

Referenced by clear(), finalAnalysis(), init(), and xtalAnalysis().

◆ nPrimary

int HcalTB04Analysis::nPrimary

private

Definition at line 129 of file HcalTB04Analysis.cc.

Referenced by clear(), fillBuffer(), and fillEvent().

◆ nTower

int HcalTB04Analysis::nTower

private

Definition at line 124 of file HcalTB04Analysis.cc.

Referenced by clear(), finalAnalysis(), init(), and qieAnalysis().

◆ particleType

int HcalTB04Analysis::particleType

private

Definition at line 129 of file HcalTB04Analysis.cc.

Referenced by clear(), fillBuffer(), and fillEvent().

◆ phiInit

double HcalTB04Analysis::phiInit

private

Definition at line 130 of file HcalTB04Analysis.cc.

Referenced by clear(), fillBuffer(), fillEvent(), and finalAnalysis().

◆ pInit

double HcalTB04Analysis::pInit

private

Definition at line 130 of file HcalTB04Analysis.cc.

Referenced by clear(), fillBuffer(), fillEvent(), and finalAnalysis().

◆ pvFound

bool HcalTB04Analysis::pvFound

private

Definition at line 137 of file HcalTB04Analysis.cc.

Referenced by clear(), and update().

◆ pvMomentum

G4ThreeVector HcalTB04Analysis::pvMomentum

private

Definition at line 139 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), and update().

◆ pvPosition

G4ThreeVector HcalTB04Analysis::pvPosition

private

Definition at line 139 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), and update().

◆ pvType

int HcalTB04Analysis::pvType

private

Definition at line 138 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), and update().

◆ pvUVW

G4ThreeVector HcalTB04Analysis::pvUVW

private

Definition at line 139 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), and update().

◆ scaleHB0

const double HcalTB04Analysis::scaleHB0

private

Definition at line 113 of file HcalTB04Analysis.cc.

Referenced by scale().

◆ scaleHB16

const double HcalTB04Analysis::scaleHB16

private

Definition at line 113 of file HcalTB04Analysis.cc.

Referenced by scale().

◆ scaleHE0

const double HcalTB04Analysis::scaleHE0

private

Definition at line 113 of file HcalTB04Analysis.cc.

Referenced by scale().

◆ scaleHO

const double HcalTB04Analysis::scaleHO

private

Definition at line 113 of file HcalTB04Analysis.cc.

Referenced by scale().

◆ secEkin

std::vector<double> HcalTB04Analysis::secEkin

private

Definition at line 142 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), and update().

◆ secMomentum

std::vector<G4ThreeVector> HcalTB04Analysis::secMomentum

private

Definition at line 141 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), and update().

◆ secPartID

std::vector<int> HcalTB04Analysis::secPartID

private

Definition at line 140 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), and update().

◆ secTrackID

std::vector<int> HcalTB04Analysis::secTrackID

private

Definition at line 140 of file HcalTB04Analysis.cc.

Referenced by clear(), fillEvent(), and update().

◆ shortLivedSecondaries

std::vector<int> HcalTB04Analysis::shortLivedSecondaries

private

Definition at line 143 of file HcalTB04Analysis.cc.

Referenced by clear(), and update().

◆ type

const int HcalTB04Analysis::type

private

Definition at line 111 of file HcalTB04Analysis.cc.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HcalTB04Analysis() [1/2]

◆ HcalTB04Analysis() [2/2]

◆ ~HcalTB04Analysis()

Member Function Documentation

◆ clear()

◆ fillBuffer()

◆ fillEvent()

◆ finalAnalysis()

◆ init()

◆ operator=()

◆ produce()

◆ qieAnalysis()

◆ scale()

◆ timeOfFlight()

◆ unitID()

◆ update() [1/4]

◆ update() [2/4]

◆ update() [3/4]

◆ update() [4/4]

◆ xtalAnalysis()

Member Data Documentation

◆ beamline_RM

◆ beamOffset

◆ count

◆ ecalHitCache

◆ ecalNoise

◆ eecalq

◆ eecals

◆ ehcalq

◆ ehcals

◆ enois

◆ eqeta

◆ eqie

◆ eqlay

◆ eqphi

◆ eseta

◆ esime

◆ esimh

◆ eslay

◆ esphi

◆ etaInit

◆ etotq

◆ etots

◆ evNum

◆ hcalHitCache

◆ hcalHitLayer

◆ hcalOnly

◆ histo

◆ iceta

◆ icphi

◆ idEcal

◆ idHcal

◆ idTower

◆ idXtal

◆ m_Anal

◆ mode

◆ myQie

◆ names

◆ nCrystal

◆ nPrimary

◆ nTower

◆ particleType

◆ phiInit

◆ pInit

◆ pvFound

◆ pvMomentum

◆ pvPosition

◆ pvType

◆ pvUVW

◆ scaleHB0

◆ scaleHB16

◆ scaleHE0

◆ scaleHO

◆ secEkin