#include <HcalTestAnalysis.h>

Inheritance diagram for HcalTestAnalysis:

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

	~HcalTestAnalysis () override

Public Member Functions inherited from SimWatcher
	SimWatcher ()

virtual	~SimWatcher ()

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

void	slotForUpdate (const BeginOfJob * iT)

virtual	~Observer ()

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	fill (const EndOfEvent *ev)

void	layerAnalysis ()

std::vector< int >	layerGrouping (int)

void	qieAnalysis (CLHEP::HepRandomEngine *)

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

std::vector< int >	towersToAdd (int centre, int nadd)

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

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 EndOfEvent *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...

Private Attributes
int	addTower_

std::vector< CaloHit >	caloHitCache_

int	centralTower_

unsigned int	count_

double	edepEB_

double	edepEE_

double	edepHB_

double	edepHE_

double	edepHO_

double	edepl_ [20]

double	eta0_

std::string	fileName_

std::vector< int >	group_

const HcalDDDSimConstants *	hcons_

double	mudist_ [20]

std::unique_ptr< HcalQie >	myqie_

std::vector< std::string >	names_

int	nGroup_

int	nTower_

std::unique_ptr< HcalNumberingFromDDD >	numberingFromDDD_

HcalTestNumberingScheme *	org_

double	phi0_

std::vector< int >	tower_

HcalTestHistoClass *	tuples_

std::unique_ptr< HcalTestHistoManager >	tuplesManager_

Additional Inherited Members

Detailed Description

Definition at line 34 of file HcalTestAnalysis.h.

Constructor & Destructor Documentation

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

Definition at line 32 of file HcalTestAnalysis.cc.

References addTower_, centralTower_, count_, eta0_, fileName_, edm::ParameterSet::getParameter(), group_, mps_fire::i, layerGrouping(), myqie_, names_, nGroup_, nTower_, numberingFromDDD_, phi0_, AlCaHLTBitMon_QueryRunRegistry::string, tower_, towersToAdd(), and tuplesManager_.

     : addTower_(3), tuples_(nullptr), hcons_(nullptr), org_(nullptr) {
   edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("HcalTestAnalysis");
   eta0_ = m_Anal.getParameter<double>("Eta0");
   phi0_ = m_Anal.getParameter<double>("Phi0");
   int laygroup = m_Anal.getParameter<int>("LayerGrouping");
   centralTower_ = m_Anal.getParameter<int>("CentralTower");
   names_ = m_Anal.getParameter<std::vector<std::string> >("Names");
   fileName_ = m_Anal.getParameter<std::string>("FileName");
 
   tuplesManager_.reset(nullptr);
   numberingFromDDD_.reset(nullptr);
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: Initialised as observer of begin/end events"
                               << " and of G4step";
 
   count_ = 0;
   group_ = layerGrouping(laygroup);
   nGroup_ = 0;
   for (unsigned int i = 0; i < group_.size(); i++)
     if (group_[i] > nGroup_)
       nGroup_ = group_[i];
   tower_ = towersToAdd(centralTower_, addTower_);
   nTower_ = tower_.size() / 2;
 
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: initialised for " << nGroup_ << " Longitudinal groups and "
                               << nTower_ << " towers";
 
   // qie
   myqie_.reset(new HcalQie(p));
 }

HcalTestAnalysis::~HcalTestAnalysis ( )

override

Definition at line 63 of file HcalTestAnalysis.cc.

References count_, org_, and tuples_.

                                     {
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis: -------->  Total number of selected entries : " << count_;
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis: Pointers:: HistoClass " << tuples_ << ", Numbering Scheme " << org_;
 }

Member Function Documentation

void HcalTestAnalysis::fill ( const EndOfEvent * ev )

private

Definition at line 315 of file HcalTestAnalysis.cc.

References caloHitCache_, MillePedeFileConverter_cfg::e, PVValHelper::eta, CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getPosition(), CaloG4Hit::getTimeSlice(), HcalTestNumberingScheme::getUnitID(), CaloG4Hit::getUnitID(), GeV, HcalEndcap, HcalForward, cmsBatch::log, names_, numberingFromDDD_, org_, HcalTestNumberingScheme::packHcalIndex(), phi, AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), theta(), protons_cff::time, timeOfFlight(), HcalTestNumberingScheme::unpackHcalIndex(), and ecaldqm::zside().

Referenced by update().

                                                  {
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis: Fill event " << (*evt)()->GetEventID();
 
   // access to the G4 hit collections
   G4HCofThisEvent* allHC = (*evt)()->GetHCofThisEvent();
 
   int nhc = 0, neb = 0, nef = 0, j = 0;
   caloHitCache_.erase(caloHitCache_.begin(), caloHitCache_.end());
 
   // Hcal
   int HCHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names_[0]);
   CaloG4HitCollection* theHCHC = (CaloG4HitCollection*)allHC->GetHC(HCHCid);
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis :: Hit Collection for " << names_[0] << " of ID " << HCHCid
                               << " is obtained at " << theHCHC;
   if (HCHCid >= 0 && theHCHC != nullptr) {
     for (j = 0; j < theHCHC->entries(); j++) {
       CaloG4Hit* aHit = (*theHCHC)[j];
 
       double e = aHit->getEnergyDeposit() / GeV;
       double time = aHit->getTimeSlice();
 
       math::XYZPoint pos = aHit->getPosition();
       double theta = pos.theta();
       double eta = -log(tan(theta * 0.5));
       double phi = pos.phi();
 
       uint32_t unitID = aHit->getUnitID();
       int subdet, zside, layer, etaIndex, phiIndex, lay;
       org_->unpackHcalIndex(unitID, subdet, zside, layer, etaIndex, phiIndex, lay);
       double jitter = time - timeOfFlight(subdet, lay, eta);
       if (jitter < 0)
         jitter = 0;
       CaloHit hit(subdet, lay, e, eta, phi, jitter, unitID);
       caloHitCache_.push_back(hit);
       nhc++;
 
       std::string det = "HB";
       if (subdet == static_cast<int>(HcalForward)) {
         det = "HF";
       } else if (subdet == static_cast<int>(HcalEndcap)) {
         if (etaIndex <= 20) {
           det = "HES";
         } else {
           det = "HED";
         }
       }
 
       edm::LogVerbatim("HcalSim") << "HcalTest: " << det << "  layer " << std::setw(2) << layer << " 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;
     }
   }
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::HCAL hits : " << nhc;
 
   // EB
   int EBHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names_[1]);
   CaloG4HitCollection* theEBHC = (CaloG4HitCollection*)allHC->GetHC(EBHCid);
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis :: Hit Collection for " << names_[1] << " of ID " << EBHCid
                               << " is obtained at " << theEBHC;
   if (EBHCid >= 0 && theEBHC != nullptr) {
     for (j = 0; j < theEBHC->entries(); j++) {
       CaloG4Hit* aHit = (*theEBHC)[j];
 
       double e = aHit->getEnergyDeposit() / GeV;
       double time = aHit->getTimeSlice();
       std::string det = "EB";
 
       math::XYZPoint pos = aHit->getPosition();
       double theta = pos.theta();
       double eta = -log(tan(theta / 2.));
       double phi = pos.phi();
 
       HcalNumberingFromDDD::HcalID id = numberingFromDDD_->unitID(eta, phi, 1, 1);
       uint32_t unitID = org_->getUnitID(id);
       int subdet, zside, layer, ieta, iphi, lay;
       org_->unpackHcalIndex(unitID, subdet, zside, layer, ieta, iphi, lay);
       subdet = 10;
       layer = 0;
       unitID = org_->packHcalIndex(subdet, zside, layer, ieta, iphi, lay);
       CaloHit hit(subdet, lay, e, eta, phi, time, unitID);
       caloHitCache_.push_back(hit);
       neb++;
       edm::LogVerbatim("HcalSim") << "HcalTest: " << det << "  layer " << std::setw(2) << layer << " 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;
     }
   }
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::EB hits : " << neb;
 
   // EE
   int EEHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names_[2]);
   CaloG4HitCollection* theEEHC = (CaloG4HitCollection*)allHC->GetHC(EEHCid);
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis :: Hit Collection for " << names_[2] << " of ID " << EEHCid
                               << " is obtained at " << theEEHC;
   if (EEHCid >= 0 && theEEHC != nullptr) {
     for (j = 0; j < theEEHC->entries(); j++) {
       CaloG4Hit* aHit = (*theEEHC)[j];
 
       double e = aHit->getEnergyDeposit() / GeV;
       double time = aHit->getTimeSlice();
       std::string det = "EE";
 
       math::XYZPoint pos = aHit->getPosition();
       double theta = pos.theta();
       double eta = -log(tan(theta / 2.));
       double phi = pos.phi();
 
       HcalNumberingFromDDD::HcalID id = numberingFromDDD_->unitID(eta, phi, 1, 1);
       uint32_t unitID = org_->getUnitID(id);
       int subdet, zside, layer, ieta, iphi, lay;
       org_->unpackHcalIndex(unitID, subdet, zside, layer, ieta, iphi, lay);
       subdet = 11;
       layer = 0;
       unitID = org_->packHcalIndex(subdet, zside, layer, ieta, iphi, lay);
       CaloHit hit(subdet, lay, e, eta, phi, time, unitID);
       caloHitCache_.push_back(hit);
       nef++;
       edm::LogVerbatim("HcalSim") << "HcalTest: " << det << "  layer " << std::setw(2) << layer << " 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;
     }
   }
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::EE hits : " << nef;
 }

void HcalTestAnalysis::layerAnalysis ( )

private

Definition at line 566 of file HcalTestAnalysis.cc.

References edepEB_, edepEE_, edepHB_, edepHE_, edepHO_, edepl_, HcalTestHistoClass::fillLayers(), mps_fire::i, MeV, mudist_, and tuples_.

Referenced by update().

                                      {
   int i = 0;
   edm::LogVerbatim("HcalSim") << "\n ===>>> HcalTestAnalysis: Energy deposit "
                               << "\n at EB : " << std::setw(6) << edepEB_ / MeV << "\n at EE : " << std::setw(6)
                               << edepEE_ / MeV << "\n at HB : " << std::setw(6) << edepHB_ / MeV
                               << "\n at HE : " << std::setw(6) << edepHE_ / MeV << "\n at HO : " << std::setw(6)
                               << edepHO_ / MeV << "\n ---- HcalTestAnalysis: Energy deposit in Layers";
   for (i = 0; i < 20; i++)
     edm::LogVerbatim("HcalSim") << " Layer " << std::setw(2) << i << " E " << std::setw(8) << edepl_[i] / MeV << " MeV";
 
   tuples_->fillLayers(edepl_, edepHO_, edepHB_ + edepHE_, mudist_);
 }

std::vector< int > HcalTestAnalysis::layerGrouping ( int group )

private

Definition at line 68 of file HcalTestAnalysis.cc.

References mps_fire::i, and groupFilesInBlocks::temp.

Referenced by HcalTestAnalysis().

                                                         {
   std::vector<int> temp(19);
   if (group <= 1) {
     int grp[19] = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2};
     for (int i = 0; i < 19; i++)
       temp[i] = grp[i];
   } else if (group == 2) {
     int grp[19] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19};
     for (int i = 0; i < 19; i++)
       temp[i] = grp[i];
   } else if (group == 3) {
     int grp[19] = {1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7};
     for (int i = 0; i < 19; i++)
       temp[i] = grp[i];
   } else if (group == 4) {
     int grp[19] = {1, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 5, 6, 6, 6, 6, 6, 7, 7};
     for (int i = 0; i < 19; i++)
       temp[i] = grp[i];
   } else {
     int grp[19] = {1, 1, 2, 3, 3, 4, 4, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 7, 7};
     for (int i = 0; i < 19; i++)
       temp[i] = grp[i];
   }
 
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: Layer Grouping ";
   for (int i = 0; i < 19; i++)
     edm::LogVerbatim("HcalSim") << "HcalTestAnalysis: Group[" << i << "] = " << temp[i];
   return temp;
 }

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

private

Definition at line 444 of file HcalTestAnalysis.cc.

References addTower_, caloHitCache_, hippyaddtobaddatafiles::cd(), centralTower_, CaloHit::det(), CaloHit::e(), eta0_, HcalTestHistoClass::fillHits(), HcalTestHistoClass::fillQie(), HcalTestNumberingScheme::getUnitID(), edmIntegrityCheck::group, group_, HcalBarrel, HcalEndcap, hfClusterShapes_cfi::hits, mps_fire::i, CaloHit::id(), training_settings::idx, CaloHit::layer(), myqie_, nGroup_, nt, nTower_, numberingFromDDD_, org_, phi0_, AlCaHLTBitMon_QueryRunRegistry::string, tower_, tuples_, HcalTestNumberingScheme::unpackHcalIndex(), and ecaldqm::zside().

Referenced by update().

                                                                {
   //Fill tuple with hit information
   int hittot = caloHitCache_.size();
   tuples_->fillHits(caloHitCache_);
 
   //Get the index of the central tower
   HcalNumberingFromDDD::HcalID id = numberingFromDDD_->unitID(eta0_, phi0_, 1, 1);
   uint32_t unitID = org_->getUnitID(id);
   int subdet, zside, layer, ieta, iphi, lay;
   org_->unpackHcalIndex(unitID, subdet, zside, layer, ieta, iphi, lay);
   int laymax = 0;
   std::string det = "Unknown";
   if (subdet == static_cast<int>(HcalBarrel)) {
     laymax = 4;
     det = "HB";
   } else if (subdet == static_cast<int>(HcalEndcap)) {
     laymax = 2;
     det = "HES";
   }
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::Qie: " << det << " Eta " << ieta << " Phi " << iphi << " Laymax "
                               << laymax << " Hits " << hittot;
 
   if (laymax > 0 && hittot > 0) {
     std::vector<CaloHit> hits(hittot);
     std::vector<double> eqielay(80, 0.0), esimlay(80, 0.0), esimtot(4, 0.0);
     std::vector<double> eqietow(200, 0.0), esimtow(200, 0.0), eqietot(4, 0.0);
     int etac = (centralTower_ / 100) % 100;
     int phic = (centralTower_ % 100);
 
     for (int layr = 0; layr < nGroup_; layr++) {
       /*
       int layx, layy=20;
       for (int i=0; i<20; i++) 
     if (group_[i] == layr+1 && i < layy) layy = i+1;
       if (subdet == static_cast<int>(HcalBarrel)) {
     if      (layy < 2)    layx = 0;
     else if (layy < 17)   layx = 1;
     else if (layy == 17)  layx = 2;
     else                  layx = 3;
       } else {
     if      (layy < 2)    layx = 0;
     else                  layx = 1;
       }
       */
       for (int it = 0; it < nTower_; it++) {
         int nhit = 0;
         double esim = 0;
         for (int k1 = 0; k1 < hittot; k1++) {
           CaloHit hit = caloHitCache_[k1];
           int subdetc = hit.det();
           int layer = hit.layer();
           int group = 0;
           if (layer > 0 && layer < 20)
             group = group_[layer];
           if (subdetc == subdet && group == layr + 1) {
             int zsidec, ietac, iphic, idx;
             unitID = hit.id();
             org_->unpackHcalIndex(unitID, subdetc, zsidec, layer, ietac, iphic, lay);
             if (etac > 0 && phic > 0) {
               idx = ietac * 100 + iphic;
             } else if (etac > 0) {
               idx = ietac * 100;
             } else if (phic > 0) {
               idx = iphic;
             } else {
               idx = 0;
             }
             if (zsidec == zside && idx == tower_[it]) {
               hits[nhit] = hit;
               edm::LogVerbatim("HcalSim") << "HcalTest: Hit " << nhit << " " << hit;
               nhit++;
               esim += hit.e();
             }
           }
         }
 
         std::vector<int> cd = myqie_->getCode(nhit, hits, engine);
         double eqie = myqie_->getEnergy(cd);
 
         edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::Qie: Energy in layer " << layr << " Sim " << esim
                                     << " After QIE " << eqie;
         for (int i = 0; i < 4; i++) {
           if (tower_[nTower_ + it] <= i) {
             esimtot[i] += esim;
             eqietot[i] += eqie;
             esimlay[20 * i + layr] += esim;
             eqielay[20 * i + layr] += eqie;
             esimtow[50 * i + it] += esim;
             eqietow[50 * i + it] += eqie;
           }
         }
       }
     }
     edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::Qie: Total energy " << esimtot[3] << " (SimHit) " << eqietot[3]
                                 << " (After QIE)";
 
     std::vector<double> latphi(10);
     int nt = 2 * addTower_ + 1;
     for (int it = 0; it < nt; it++)
       latphi[it] = it - addTower_;
     for (int i = 0; i < 4; i++) {
       double scals = 1, scalq = 1;
       std::vector<double> latfs(10, 0.), latfq(10, 0.), longs(20), longq(20);
       if (esimtot[i] > 0)
         scals = 1. / esimtot[i];
       if (eqietot[i] > 0)
         scalq = 1. / eqietot[i];
       for (int it = 0; it < nTower_; it++) {
         int phib = it % nt;
         latfs[phib] += scals * esimtow[50 * i + it];
         latfq[phib] += scalq * eqietow[50 * i + it];
       }
       for (int layr = 0; layr <= nGroup_; layr++) {
         longs[layr] = scals * esimlay[20 * i + layr];
         longq[layr] = scalq * eqielay[20 * i + layr];
       }
       tuples_->fillQie(i, esimtot[i], eqietot[i], nGroup_, longs, longq, nt, latphi, latfs, latfq);
     }
   }
 }

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

private

Definition at line 580 of file HcalTestAnalysis.cc.

References funct::cos(), JetChargeProducer_cfi::exp, HcalBarrel, funct::sin(), theta(), and tmp.

Referenced by fill().

                                                                     {
   double theta = 2.0 * atan(exp(-eta));
   double dist = 0.;
   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 < 20)
       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 < 20)
       dist += zLay[layer - 1] * mm / cos(theta);
   }
   double tmp = dist / c_light / ns;
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::timeOfFlight " << tmp << " for det/lay " << det << " " << layer
                               << " eta/theta " << eta << " " << theta / deg << " dist " << dist;
   return tmp;
 }

std::vector< int > HcalTestAnalysis::towersToAdd	(	int	centre,
		int	nadd
	)

private

Definition at line 98 of file HcalTestAnalysis.cc.

References funct::abs(), PVValHelper::eta, mps_fire::i, SiStripPI::max, phi, phimax, phimin, and groupFilesInBlocks::temp.

Referenced by HcalTestAnalysis().

                                                                  {
   int etac = (centre / 100) % 100;
   int phic = (centre % 100);
   int etamin, etamax, phimin, phimax;
   if (etac > 0) {
     etamin = etac - nadd;
     etamax = etac + nadd;
   } else {
     etamin = etac;
     etamax = etac;
   }
   if (phic > 0) {
     phimin = phic - nadd;
     phimax = phic + nadd;
   } else {
     phimin = phic;
     phimax = phic;
   }
 
   int nbuf, kount = 0;
   nbuf = (etamax - etamin + 1) * (phimax - phimin + 1);
   std::vector<int> temp(2 * nbuf);
   for (int eta = etamin; eta <= etamax; eta++) {
     for (int phi = phimin; phi <= phimax; phi++) {
       temp[kount] = (eta * 100 + phi);
       temp[kount + nbuf] = std::max(abs(eta - etac), abs(phi - phic));
       kount++;
     }
   }
 
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: Towers to be considered for Central " << centre << " and " << nadd
                               << " on either side";
   for (int i = 0; i < nbuf; i++)
     edm::LogVerbatim("HcalSim") << "HcalTestAnalysis: Tower[" << std::setw(3) << i << "] " << temp[i] << " "
                                 << temp[nbuf + i];
   return temp;
 }

void HcalTestAnalysis::update ( const BeginOfJob * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 138 of file HcalTestAnalysis.cc.

References fileName_, hcons_, names_, numberingFromDDD_, org_, edm::ESHandle< T >::product(), and tuplesManager_.

                                                    {
   // Numbering From DDD
   edm::ESHandle<HcalDDDSimConstants> hdc;
   (*job)()->get<HcalSimNumberingRecord>().get(hdc);
   hcons_ = hdc.product();
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: Initialise HcalNumberingFromDDD for " << names_[0];
   numberingFromDDD_.reset(new HcalNumberingFromDDD(hcons_));
 
   // Ntuples
   tuplesManager_.reset(new HcalTestHistoManager(fileName_));
 
   // Numbering scheme
   org_ = new HcalTestNumberingScheme(false);
 }

void HcalTestAnalysis::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 154 of file HcalTestAnalysis.cc.

References HcalDDDSimConstants::cell(), centralTower_, HcalCellType::HcalCell::eta, PVValHelper::eta, eta0_, HcalBarrel, HcalEndcap, HcalForward, hcons_, heppy_loop::loop, names_, HcalCellType::HcalCell::ok, org_, phi, HcalCellType::HcalCell::phi, phi0_, sd, HCalSD::setNumberingScheme(), AlCaHLTBitMon_QueryRunRegistry::string, and tmp.

                                                    {
   int irun = (*run)()->GetRunID();
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: Begin of Run = " << irun;
 
   bool loop = true, eta = true, phi = true;
   int etac = (centralTower_ / 100) % 100;
   if (etac == 0) {
     etac = 1;
     eta = false;
   }
   int phic = (centralTower_ % 100);
   if (phic == 0) {
     phic = 1;
     phi = false;
   }
   int idet = static_cast<int>(HcalBarrel);
   while (loop) {
     HcalCellType::HcalCell tmp = hcons_->cell(idet, 1, 1, etac, phic);
     if (tmp.ok) {
       if (eta)
         eta0_ = tmp.eta;
       if (phi)
         phi0_ = tmp.phi;
       loop = false;
     } else if (idet == static_cast<int>(HcalBarrel)) {
       idet = static_cast<int>(HcalEndcap);
     } else if (idet == static_cast<int>(HcalEndcap)) {
       idet = static_cast<int>(HcalForward);
     } else {
       loop = false;
     }
   }
 
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: Central Tower " << centralTower_
                               << " corresponds to eta0 = " << eta0_ << " phi0 = " << phi0_;
 
   std::string sdname = names_[0];
   G4SDManager* sd = G4SDManager::GetSDMpointerIfExist();
   if (sd != nullptr) {
     G4VSensitiveDetector* aSD = sd->FindSensitiveDetector(sdname);
     if (aSD == nullptr) {
       edm::LogWarning("HcalSim") << "HcalTestAnalysis::beginOfRun: No SD with "
                                  << "name " << sdname << " in this Setup";
     } else {
       HCalSD* theCaloSD = dynamic_cast<HCalSD*>(aSD);
       edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::beginOfRun: Finds SD with name " << theCaloSD->GetName()
                                   << " in this Setup";
       if (org_) {
         theCaloSD->setNumberingScheme(org_);
         edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::beginOfRun: set a new numbering scheme";
       }
     }
   } else {
     edm::LogWarning("HcalSim") << "HcalTestAnalysis::beginOfRun: Could not get"
                                << " SD Manager!";
   }
 }

void HcalTestAnalysis::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 213 of file HcalTestAnalysis.cc.

References edepEB_, edepEE_, edepHB_, edepHE_, edepHO_, edepl_, mps_fire::i, mudist_, HcalTestHistoClass::setCounters(), and tuples_.

                                                      {
   // create tuple object
   tuples_ = new HcalTestHistoClass();
   // Reset counters
   tuples_->setCounters();
 
   int i = 0;
   edepEB_ = edepEE_ = edepHB_ = edepHE_ = edepHO_ = 0.;
   for (i = 0; i < 20; i++)
     edepl_[i] = 0.;
   for (i = 0; i < 20; i++)
     mudist_[i] = -1.;
 
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis: Begin of event = " << (*evt)()->GetEventID();
 }

void HcalTestAnalysis::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 293 of file HcalTestAnalysis.cc.

References count_, fill(), layerAnalysis(), qieAnalysis(), tuples_, and tuplesManager_.

                                                    {
   ++count_;
   // Fill event input
   fill(evt);
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: ---  after Fill";
 
   // Qie analysis
   CLHEP::HepRandomEngine* engine = G4Random::getTheEngine();
   qieAnalysis(engine);
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: ---  after QieAnalysis";
 
   // Layers tuples filling
   layerAnalysis();
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: ---  after LayerAnalysis";
 
   // Writing the data to the Tree
   tuplesManager_.get()->fillTree(tuples_);  // (no need to delete it...)
   tuples_ = nullptr;                        // but avoid to reuse it...
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: --- after fillTree";
 }

void HcalTestAnalysis::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 230 of file HcalTestAnalysis.cc.

References edepEB_, edepEE_, edepHB_, edepHE_, edepHO_, edepl_, PVValHelper::eta, eta0_, cmsBatch::log, MeV, mudist_, dataset::name, phi, phi0_, mathSSE::sqrt(), funct::tan(), and theta().

                                                  {
   if (aStep != nullptr) {
     G4VPhysicalVolume* curPV = aStep->GetPreStepPoint()->GetPhysicalVolume();
     G4String name = curPV->GetName();
     name.assign(name, 0, 3);
     double edeposit = aStep->GetTotalEnergyDeposit();
     int layer = -1;
     if (name == "EBR") {
       edepEB_ += edeposit;
     } else if (name == "EFR") {
       edepEE_ += edeposit;
     } else if (name == "HBS") {
       layer = (curPV->GetCopyNo() / 10) % 100;
       if (layer >= 0 && layer < 17) {
         edepHB_ += edeposit;
       } else {
         edm::LogWarning("HcalSim") << "HcalTestAnalysis::Error in HB " << curPV->GetName() << curPV->GetCopyNo();
         layer = -1;
       }
     } else if (name == "HES") {
       layer = (curPV->GetCopyNo() / 10) % 100;
       if (layer >= 0 && layer < 19) {
         edepHE_ += edeposit;
       } else {
         edm::LogWarning("HcalSim") << "HcalTestAnalysis::Error in HE " << curPV->GetName() << curPV->GetCopyNo();
         layer = -1;
       }
     } else if (name == "HTS") {
       layer = (curPV->GetCopyNo() / 10) % 100;
       if (layer >= 17 && layer < 20) {
         edepHO_ += edeposit;
       } else {
         edm::LogWarning("HcalSim") << "HcalTestAnalysis::Error in HO " << curPV->GetName() << curPV->GetCopyNo();
         layer = -1;
       }
     }
     if (layer >= 0 && layer < 20) {
       edepl_[layer] += edeposit;
 
       // Calculate the distance if it is a muon
       G4String part = aStep->GetTrack()->GetDefinition()->GetParticleName();
       if ((part == "mu-" || part == "mu+") && mudist_[layer] < 0) {
         math::XYZPoint pos(aStep->GetPreStepPoint()->GetPosition().x(),
                            aStep->GetPreStepPoint()->GetPosition().y(),
                            aStep->GetPreStepPoint()->GetPosition().z());
         double theta = pos.theta();
         double eta = -log(tan(theta * 0.5));
         double phi = pos.phi();
         double dist = sqrt((eta - eta0_) * (eta - eta0_) + (phi - phi0_) * (phi - phi0_));
         mudist_[layer] = dist * std::sqrt(pos.perp2());
       }
     }
 
     if (layer >= 0 && layer < 20) {
       edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: G4Step: " << name << " Layer " << std::setw(3) << layer
                                   << " Edep " << std::setw(6) << edeposit / MeV << " MeV";
     }
   } else {
     edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: G4Step: Null Step";
   }
 }