Inheritance diagram for HcalTestAnalysis:

Public Member Functions
void	beginRun (edm::EventSetup const &) override

	HcalTestAnalysis (const edm::ParameterSet &p)

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

void	registerConsumes (edm::ConsumesCollector) override

	~HcalTestAnalysis () 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
bool	isMT () const

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

	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	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 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_

const int	centralTower_

unsigned int	count_

edm::ESGetToken< HcalDDDSimConstants, HcalSimNumberingRecord >	ddconsToken_

double	edepEB_

double	edepEE_

double	edepHB_

double	edepHE_

double	edepHO_

double	edepl_ [20]

double	eta0_

const std::string	fileName_

std::vector< int >	group_

const HcalDDDSimConstants *	hcons_

const int	laygroup_

const edm::ParameterSet	m_Anal

double	mudist_ [20]

std::unique_ptr< HcalQie >	myqie_

const std::vector< std::string >	names_

int	nGroup_

int	nTower_

std::unique_ptr< HcalNumberingFromDDD >	numberingFromDDD_

std::unique_ptr< HcalTestNumberingScheme >	org_

double	phi0_

std::vector< int >	tower_

std::unique_ptr< HcalTestHistoClass >	tuples_

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 42 of file HcalTestAnalysis.cc.

Constructor & Destructor Documentation

◆ HcalTestAnalysis()

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

Definition at line 104 of file HcalTestAnalysis.cc.

References addTower_, centralTower_, count_, group_, mps_fire::i, layerGrouping(), laygroup_, myqie_, nGroup_, nTower_, numberingFromDDD_, org_, AlCaHLTBitMon_ParallelJobs::p, tower_, towersToAdd(), and tuples_.

     : addTower_(3),
       m_Anal(p.getParameter<edm::ParameterSet>("HcalTestAnalysis")),
       eta0_(m_Anal.getParameter<double>("Eta0")),
       phi0_(m_Anal.getParameter<double>("Phi0")),
       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")),
       hcons_(nullptr) {
   org_.reset(nullptr);
 
   tuples_.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_ = std::make_unique<HcalQie>(p);
 
   produces<HcalTestHistoClass>();
 }

◆ ~HcalTestAnalysis()

HcalTestAnalysis::~HcalTestAnalysis ( )

override

Definition at line 139 of file HcalTestAnalysis.cc.

References count_, and org_.

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

Member Function Documentation

◆ beginRun()

void HcalTestAnalysis::beginRun ( edm::EventSetup const & es )

overridevirtual

Reimplemented from SimWatcher.

Definition at line 224 of file HcalTestAnalysis.cc.

References ddconsToken_, edm::EventSetup::getData(), hcons_, names_, numberingFromDDD_, and org_.

                                                        {
   // Numbering From DDD
   hcons_ = &es.getData(ddconsToken_);
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis:: Initialise HcalNumberingFromDDD for " << names_[0];
   numberingFromDDD_ = std::make_unique<HcalNumberingFromDDD>(hcons_);
 
   // Numbering scheme
   org_ = std::make_unique<HcalTestNumberingScheme>(false);
 }

◆ fill()

void HcalTestAnalysis::fill ( const EndOfEvent * ev )

private

Definition at line 391 of file HcalTestAnalysis.cc.

References caloHitCache_, MillePedeFileConverter_cfg::e, PVValHelper::eta, CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getPosition(), CaloG4Hit::getTimeSlice(), CaloG4Hit::getUnitID(), HcalEndcap, HcalForward, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, dqmiolumiharvest::j, nano_mu_digi_cff::layer, dqm-mbProfile::log, names_, numberingFromDDD_, org_, phi, AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), theta(), hcalRecHitTable_cff::time, timeOfFlight(), 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;
   int hchc_entries = theHCHC->entries();
   if (HCHCid >= 0 && theHCHC != nullptr) {
     for (j = 0; j < hchc_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;
   int ebhc_entries = theEBHC->entries();
   if (EBHCid >= 0 && theEBHC != nullptr) {
     for (j = 0; j < ebhc_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;
   int eehc_entries = theEEHC->entries();
   if (EEHCid >= 0 && theEEHC != nullptr) {
     for (j = 0; j < eehc_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;
 }

◆ layerAnalysis()

void HcalTestAnalysis::layerAnalysis ( )

private

Definition at line 645 of file HcalTestAnalysis.cc.

References edepEB_, edepEE_, edepHB_, edepHE_, edepHO_, edepl_, mps_fire::i, 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_);
 }

◆ layerGrouping()

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

private

Definition at line 154 of file HcalTestAnalysis.cc.

References watchdog::group, 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;
 }

◆ produce()

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

overridevirtual

Implements SimProducer.

Definition at line 149 of file HcalTestAnalysis.cc.

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

                                                                 {
   e.put(std::move(tuples_));
   tuples_.reset(nullptr);
 }

◆ qieAnalysis()

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

private

Definition at line 523 of file HcalTestAnalysis.cc.

References addTower_, caloHitCache_, hippyaddtobaddatafiles::cd(), centralTower_, eta0_, watchdog::group, group_, HcalBarrel, HcalEndcap, hfClusterShapes_cfi::hits, mps_fire::i, hit::id, heavyIonCSV_trainingSettings::idx, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, nano_mu_digi_cff::layer, myqie_, nGroup_, nt, nTower_, numberingFromDDD_, org_, phi0_, AlCaHLTBitMon_QueryRunRegistry::string, tower_, tuples_, 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);
     }
   }
 }

◆ registerConsumes()

void HcalTestAnalysis::registerConsumes ( edm::ConsumesCollector cc )

overridevirtual

Reimplemented from SimWatcher.

Definition at line 144 of file HcalTestAnalysis.cc.

References edm::BeginRun, gpuPixelDoublets::cc, and ddconsToken_.

                                                                {
   ddconsToken_ = cc.esConsumes<HcalDDDSimConstants, HcalSimNumberingRecord, edm::Transition::BeginRun>();
   edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::Initialize ESGetToken for HcalDDDSimConstants";
 }

◆ timeOfFlight()

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

private

Definition at line 659 of file HcalTestAnalysis.cc.

References funct::cos(), PVValHelper::eta, JetChargeProducer_cfi::exp, HcalBarrel, nano_mu_digi_cff::layer, funct::sin(), theta(), and createJobs::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;
 }

◆ towersToAdd()

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

private

Definition at line 184 of file HcalTestAnalysis.cc.

References funct::abs(), PVValHelper::eta, TrackingDataMCValidation_Standalone_cff::etamax, TrackingDataMCValidation_Standalone_cff::etamin, 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;
 }

◆ update() [1/4]

void HcalTestAnalysis::update ( const BeginOfRun * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 235 of file HcalTestAnalysis.cc.

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

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

                                                    {
   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_.get()) {
         theCaloSD->setNumberingScheme(org_.get());
         edm::LogVerbatim("HcalSim") << "HcalTestAnalysis::beginOfRun: set a new numbering scheme";
       }
     }
   } else {
     edm::LogWarning("HcalSim") << "HcalTestAnalysis::beginOfRun: Could not get"
                                << " SD Manager!";
   }
 }

◆ update() [2/4]

void HcalTestAnalysis::update ( const BeginOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 294 of file HcalTestAnalysis.cc.

References edepEB_, edepEE_, edepHB_, edepHE_, edepHO_, edepl_, mps_fire::i, mudist_, and tuples_.

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

                                                      {
   // create tuple object
   tuples_ = std::make_unique<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();
 }

◆ update() [3/4]

void HcalTestAnalysis::update ( const EndOfEvent * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent *>.

Definition at line 374 of file HcalTestAnalysis.cc.

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

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

                                                    {
   ++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";
 }

◆ update() [4/4]

void HcalTestAnalysis::update ( const G4Step * )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 311 of file HcalTestAnalysis.cc.

References edepEB_, edepEE_, edepHB_, edepHE_, edepHO_, edepl_, PVValHelper::eta, eta0_, nano_mu_digi_cff::layer, dqm-mbProfile::log, mudist_, Skims_PA_cff::name, phi, phi0_, mathSSE::sqrt(), funct::tan(), and theta().

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

                                                  {
   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";
   }
 }

Member Data Documentation

◆ addTower_

int HcalTestAnalysis::addTower_

private

Definition at line 73 of file HcalTestAnalysis.cc.

Referenced by HcalTestAnalysis(), and qieAnalysis().

◆ caloHitCache_

std::vector<CaloHit> HcalTestAnalysis::caloHitCache_

private

Definition at line 93 of file HcalTestAnalysis.cc.

Referenced by fill(), and qieAnalysis().

◆ centralTower_

const int HcalTestAnalysis::centralTower_

private

Definition at line 81 of file HcalTestAnalysis.cc.

Referenced by HcalTestAnalysis(), qieAnalysis(), and update().

◆ count_

unsigned int HcalTestAnalysis::count_

private

Definition at line 98 of file HcalTestAnalysis.cc.

Referenced by HcalTestAnalysis(), update(), and ~HcalTestAnalysis().

◆ ddconsToken_

edm::ESGetToken<HcalDDDSimConstants, HcalSimNumberingRecord> HcalTestAnalysis::ddconsToken_

private

Definition at line 87 of file HcalTestAnalysis.cc.

Referenced by beginRun(), and registerConsumes().

◆ edepEB_

double HcalTestAnalysis::edepEB_

private

Definition at line 99 of file HcalTestAnalysis.cc.

Referenced by layerAnalysis(), and update().

◆ edepEE_

double HcalTestAnalysis::edepEE_

private

Definition at line 99 of file HcalTestAnalysis.cc.

Referenced by layerAnalysis(), and update().

◆ edepHB_

double HcalTestAnalysis::edepHB_

private

Definition at line 99 of file HcalTestAnalysis.cc.

Referenced by layerAnalysis(), and update().

◆ edepHE_

double HcalTestAnalysis::edepHE_

private

Definition at line 99 of file HcalTestAnalysis.cc.

Referenced by layerAnalysis(), and update().

◆ edepHO_

double HcalTestAnalysis::edepHO_

private

Definition at line 100 of file HcalTestAnalysis.cc.

Referenced by layerAnalysis(), and update().

◆ edepl_

double HcalTestAnalysis::edepl_[20]

private

Definition at line 100 of file HcalTestAnalysis.cc.

Referenced by layerAnalysis(), and update().

◆ eta0_

double HcalTestAnalysis::eta0_

private

Definition at line 80 of file HcalTestAnalysis.cc.

Referenced by qieAnalysis(), and update().

◆ fileName_

const std::string HcalTestAnalysis::fileName_

private

Definition at line 84 of file HcalTestAnalysis.cc.

◆ group_

std::vector<int> HcalTestAnalysis::group_

private

Definition at line 94 of file HcalTestAnalysis.cc.

Referenced by HcalTestAnalysis(), and qieAnalysis().

◆ hcons_

const HcalDDDSimConstants* HcalTestAnalysis::hcons_

private

Definition at line 89 of file HcalTestAnalysis.cc.

Referenced by beginRun(), and update().

◆ laygroup_

const int HcalTestAnalysis::laygroup_

private

Definition at line 81 of file HcalTestAnalysis.cc.

Referenced by HcalTestAnalysis().

◆ m_Anal

const edm::ParameterSet HcalTestAnalysis::m_Anal

private

Definition at line 79 of file HcalTestAnalysis.cc.

◆ mudist_

double HcalTestAnalysis::mudist_[20]

private

Definition at line 101 of file HcalTestAnalysis.cc.

Referenced by layerAnalysis(), and update().

◆ myqie_

std::unique_ptr<HcalQie> HcalTestAnalysis::myqie_

private

Definition at line 72 of file HcalTestAnalysis.cc.

Referenced by HcalTestAnalysis(), and qieAnalysis().

◆ names_

const std::vector<std::string> HcalTestAnalysis::names_

private

Definition at line 82 of file HcalTestAnalysis.cc.

Referenced by beginRun(), fill(), and update().

◆ nGroup_

int HcalTestAnalysis::nGroup_

private

Definition at line 95 of file HcalTestAnalysis.cc.

Referenced by HcalTestAnalysis(), and qieAnalysis().

◆ nTower_

int HcalTestAnalysis::nTower_

private

Definition at line 95 of file HcalTestAnalysis.cc.

Referenced by HcalTestAnalysis(), and qieAnalysis().

◆ numberingFromDDD_

std::unique_ptr<HcalNumberingFromDDD> HcalTestAnalysis::numberingFromDDD_

private

Definition at line 88 of file HcalTestAnalysis.cc.

Referenced by beginRun(), fill(), HcalTestAnalysis(), and qieAnalysis().

◆ org_

std::unique_ptr<HcalTestNumberingScheme> HcalTestAnalysis::org_

private

Definition at line 90 of file HcalTestAnalysis.cc.

Referenced by beginRun(), fill(), HcalTestAnalysis(), qieAnalysis(), update(), and ~HcalTestAnalysis().

◆ phi0_

double HcalTestAnalysis::phi0_

private

Definition at line 80 of file HcalTestAnalysis.cc.

Referenced by qieAnalysis(), and update().

◆ tower_

std::vector<int> HcalTestAnalysis::tower_

private

Definition at line 94 of file HcalTestAnalysis.cc.

Referenced by HcalTestAnalysis(), and qieAnalysis().

◆ tuples_

std::unique_ptr<HcalTestHistoClass> HcalTestAnalysis::tuples_

private

Definition at line 76 of file HcalTestAnalysis.cc.

Referenced by HcalTestAnalysis(), layerAnalysis(), produce(), qieAnalysis(), and update().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HcalTestAnalysis()

◆ ~HcalTestAnalysis()

Member Function Documentation

◆ beginRun()

◆ fill()

◆ layerAnalysis()

◆ layerGrouping()

◆ produce()

◆ qieAnalysis()

◆ registerConsumes()

◆ timeOfFlight()

◆ towersToAdd()

◆ update() [1/4]

◆ update() [2/4]

◆ update() [3/4]

◆ update() [4/4]

Member Data Documentation

◆ addTower_

◆ caloHitCache_

◆ centralTower_

◆ count_

◆ ddconsToken_

◆ edepEB_

◆ edepEE_

◆ edepHB_

◆ edepHE_

◆ edepHO_

◆ edepl_

◆ eta0_

◆ fileName_

◆ group_

◆ hcons_

◆ laygroup_

◆ m_Anal

◆ mudist_

◆ myqie_

◆ names_

◆ nGroup_

◆ nTower_

◆ numberingFromDDD_

◆ org_

◆ phi0_

◆ tower_

◆ tuples_