Inheritance diagram for HOSimHitStudy:

Public Member Functions
	HOSimHitStudy (const edm::ParameterSet &ps)

	~HOSimHitStudy () override=default

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	analyzeHits ()

void	beginJob () override

void	beginRun (edm::Run const &, edm::EventSetup const &) override

void	endJob () override

void	endRun (edm::Run const &, edm::EventSetup const &) override

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Private Attributes
std::vector< PCaloHit >	ecalHits

TH1F *	edep_ [3]

TH1F *	edepT_ [3]

TH1F *	edepTW_ [3]

TH1F *	edepTWT_ [3]

TH1F *	edepZon_ [3]

TH1F *	edepZonT_ [3]

TH1F *	eEB_

TH1F *	eEBHB_

TH1F *	eEBHBHO_

TH1F *	eEBHBHOT_

TH1F *	eEBHBT_

TH1F *	eEBT_

TProfile *	eHO17_

TProfile *	eHO17T_

TProfile *	eHO18_

TProfile *	eHO18T_

TProfile *	eHO1_

TProfile *	eHO1T_

TProfile2D *	eHO2_

TProfile2D *	eHO2T_

TH1F *	eHOE17_ [15]

TH1F *	eHOE17T_ [15]

TH1F *	eHOE18_ [15]

TH1F *	eHOE18T_ [15]

TH1F *	eHOE_ [15]

TH1F *	eHOET_ [15]

TH1F *	eHOEta17_ [15]

TH1F *	eHOEta17T_ [15]

TH1F *	eHOEta18_ [15]

TH1F *	eHOEta18T_ [15]

TH1F *	eHOEta_ [15]

TH1F *	eHOEtaT_ [15]

double	eInc

TH1F *	eneInc_

double	etaInc

TH1F *	etaInc_

const std::string	g4Label_

std::vector< PCaloHit >	hcalHits

TH1F *	hit_ [3]

const std::vector< std::string >	hitLab_

TH1F *	hitTow_ [3]

const double	maxEnergy_

TProfile *	nHO1_

TProfile *	nHO1T_

TProfile2D *	nHO2_

TProfile2D *	nHO2T_

TH1F *	nHOE1_ [15]

TH1F *	nHOE1T_ [15]

TProfile *	nHOE2_ [15]

TProfile *	nHOE2T_ [15]

TH1F *	nHOEta1_ [15]

TH1F *	nHOEta1T_ [15]

TProfile *	nHOEta2_ [15]

TProfile *	nHOEta2T_ [15]

double	phiInc

TH1F *	phiInc_

const bool	print_

const double	scaleEB_

const double	scaleHB_

const double	scaleHO_

const bool	scheme_

const double	tcut_

TH1F *	time_ [3]

const edm::EDGetTokenT< edm::HepMCProduct >	tok_evt_

const std::vector< edm::EDGetTokenT< edm::PCaloHitContainer > >	toks_calo_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Detailed Description

Definition at line 41 of file HOSimHitStudy.cc.

Constructor & Destructor Documentation

◆ HOSimHitStudy()

HOSimHitStudy::HOSimHitStudy ( const edm::ParameterSet & ps )

Definition at line 81 of file HOSimHitStudy.cc.

References g4Label_, hitLab_, Skims_PA_cff::name, AlCaHLTBitMon_QueryRunRegistry::string, and edm::vector_transform().

     : g4Label_(ps.getUntrackedParameter<std::string>("ModuleLabel", "g4SimHits")),
       hitLab_(ps.getParameter<std::vector<std::string>>("HitCollection")),
       maxEnergy_(ps.getUntrackedParameter<double>("MaxEnergy", 200.0)),
       scaleEB_(ps.getUntrackedParameter<double>("ScaleEB", 1.0)),
       scaleHB_(ps.getUntrackedParameter<double>("ScaleHB", 100.0)),
       scaleHO_(ps.getUntrackedParameter<double>("ScaleHO", 2.0)),
       tcut_(ps.getUntrackedParameter<double>("TimeCut", 100.0)),
       scheme_(ps.getUntrackedParameter<bool>("TestNumbering", false)),
       print_(ps.getUntrackedParameter<bool>("PrintExcessEnergy", true)),
       tok_evt_(consumes<edm::HepMCProduct>(
           edm::InputTag(ps.getUntrackedParameter<std::string>("SourceLabel", "VtxSmeared")))),
       toks_calo_{edm::vector_transform(hitLab_, [this](const std::string &name) {
         return consumes<edm::PCaloHitContainer>(edm::InputTag{g4Label_, name});
       })} {
   usesResource(TFileService::kSharedResource);
 
   edm::LogVerbatim("HitStudy") << "HOSimHitStudy::Module Label: " << g4Label_ << "   Hits: " << hitLab_[0] << ", "
                                << hitLab_[1] << "   MaxEnergy: " << maxEnergy_ << "   Scale factor for EB " << scaleEB_
                                << ", for HB " << scaleHB_ << " and for HO " << scaleHO_ << " time Cut " << tcut_;
 }

◆ ~HOSimHitStudy()

HOSimHitStudy::~HOSimHitStudy ( )

overridedefault

Member Function Documentation

◆ analyze()

void HOSimHitStudy::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 364 of file HOSimHitStudy.cc.

References analyzeHits(), MillePedeFileConverter_cfg::e, ecalHits, eInc, etaInc, edm::HepMCProduct::GetEvent(), hcalHits, hitLab_, mps_fire::i, edm::HandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, phiInc, tok_evt_, and toks_calo_.

                                                                     {
   edm::LogVerbatim("HitStudy") << "HOSimHitStudy::Run = " << e.id().run() << " Event = " << e.id().event();
 
   const edm::Handle<edm::HepMCProduct> &EvtHandle = e.getHandle(tok_evt_);
   const HepMC::GenEvent *myGenEvent = EvtHandle->GetEvent();
 
   eInc = etaInc = phiInc = 0;
   HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
   if (p != myGenEvent->particles_end()) {
     eInc = (*p)->momentum().e();
     etaInc = (*p)->momentum().eta();
     phiInc = (*p)->momentum().phi();
   }
 
   edm::LogVerbatim("HitStudy") << "HOSimHitStudy::Energy = " << eInc << " Eta = " << etaInc
                                << " Phi = " << phiInc / CLHEP::deg;
 
   for (int i = 0; i < 2; i++) {
     bool getHits = false;
     if (i == 0)
       ecalHits.clear();
     else
       hcalHits.clear();
     const edm::Handle<edm::PCaloHitContainer> &hitsCalo = e.getHandle(toks_calo_[i]);
     if (hitsCalo.isValid())
       getHits = true;
     edm::LogVerbatim("HitStudy") << "HOSimHitStudy::Input flag " << hitLab_[i] << " getHits flag " << getHits;
 
     if (getHits) {
       unsigned int isiz;
       if (i == 0) {
         ecalHits.insert(ecalHits.end(), hitsCalo->begin(), hitsCalo->end());
         isiz = ecalHits.size();
       } else {
         hcalHits.insert(hcalHits.end(), hitsCalo->begin(), hitsCalo->end());
         isiz = hcalHits.size();
       }
       edm::LogVerbatim("HitStudy") << "HOSimHitStudy:: Hit buffer for " << hitLab_[i] << " has " << isiz << " hits";
     }
   }
   analyzeHits();
 }

◆ analyzeHits()

void HOSimHitStudy::analyzeHits ( )

protected

Definition at line 407 of file HOSimHitStudy.cc.

References funct::abs(), TauDecayModes::dec, hcalRecHitTable_cff::depth, ecalHits, edep_, edepT_, edepTW_, edepTWT_, edepZon_, edepZonT_, eEB_, eEBHB_, eEBHBHO_, eEBHBHOT_, eEBHBT_, eEBT_, eHO17_, eHO17T_, eHO18_, eHO18T_, eHO1_, eHO1T_, eHO2_, eHO2T_, eHOE17_, eHOE17T_, eHOE18_, eHOE18T_, eHOE_, eHOET_, eHOEta17_, eHOEta18_, eHOEta_, eInc, eneInc_, PVValHelper::eta, etaInc, etaInc_, HcalBarrel, hcalHits, HcalOuter, hit_, hitTow_, mps_fire::i, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, dqmiolumiharvest::j, dqmdumpme::k, maxEnergy_, nHO1_, nHO1T_, nHO2_, nHO2T_, nHOE1_, nHOE1T_, nHOE2_, nHOE2T_, nHOEta1_, nHOEta1T_, nHOEta2_, nHOEta2T_, convertSQLiteXML::ok, phi, phiInc, phiInc_, print_, scaleEB_, scaleHB_, scaleHO_, scheme_, AlCaHLTBitMon_QueryRunRegistry::string, tcut_, hcalRecHitTable_cff::time, time_, HcalTestNumberingScheme::unpackHcalIndex(), and ecaldqm::zside().

Referenced by analyze().

                                 {
   //initialize
   int nhit[3];
   double etot[3], etotT[3];
   std::vector<unsigned int> ebID, hbID, hoID;
   std::vector<double> ebEtow, hbEtow, hoEtow;
   std::vector<double> ebEtowT, hbEtowT, hoEtowT;
   for (int k = 0; k < 3; k++) {
     nhit[k] = 0;
     etot[k] = 0;
     etotT[k] = 0;
   }
   eneInc_->Fill(eInc);
   etaInc_->Fill(etaInc);
   phiInc_->Fill(phiInc);
 
   double eHO17 = 0, eHO18 = 0, eHO17T = 0, eHO18T = 0;
   double eHOE17[15], eHOE18[15], eHOE17T[15], eHOE18T[15];
   for (int k = 0; k < 15; k++) {
     eHOE17[k] = eHOE18[k] = eHOE17T[k] = eHOE18T[k] = 0;
   }
   // Loop over containers
   for (int k = 0; k < 2; k++) {
     int nHit;
     if (k == 0) {
       nHit = ecalHits.size();
     } else {
       nHit = hcalHits.size();
     }
     for (int i = 0; i < nHit; i++) {
       double edep, time;
       int indx = -1;
       unsigned int id_;
       if (k == 0) {
         indx = 0;
         edep = ecalHits[i].energy();
         time = ecalHits[i].time();
         id_ = ecalHits[i].id();
       } else {
         edep = hcalHits[i].energy();
         time = hcalHits[i].time();
         id_ = hcalHits[i].id();
         int subdet, zside, depth, eta, phi, lay;
         if (scheme_) {
           HcalTestNumberingScheme::unpackHcalIndex(id_, subdet, zside, depth, eta, phi, lay);
         } else {
           HcalDetId id = HcalDetId(id_);
           subdet = id.subdet();
           zside = id.zside();
           depth = id.depth();
           eta = id.ietaAbs();
           phi = id.iphi();
           lay = -1;
         }
         edm::LogVerbatim("HitStudy") << "HOSimHitStudy:: Hit " << k << " Subdet:" << subdet << " zside:" << zside
                                      << " depth:" << depth << " layer:" << lay << " eta:" << eta << " phi:" << phi;
         if (subdet == static_cast<int>(HcalBarrel))
           indx = 1;
         else if (subdet == static_cast<int>(HcalOuter)) {
           indx = 2;
           if (lay == 18) {
             eHO17 += edep;
             if (time < tcut_)
               eHO17T += edep;
             if (eta >= 0 && eta < 15) {
               eHOE17[eta] += edep;
               if (time < tcut_)
                 eHOE17T[eta] += edep;
             }
           } else {
             eHO18 += edep;
             if (time < tcut_)
               eHO18T += edep;
             if (eta >= 0 && eta < 15) {
               eHOE18[eta] += edep;
               if (time < tcut_)
                 eHOE18T[eta] += edep;
             }
           }
         }
       }
       if (indx >= 0) {
         double edepT = edep;
         time_[indx]->Fill(time);
         edep_[indx]->Fill(edep);
         etot[indx] += edep;
         if (time < tcut_) {
           etotT[indx] += edep;
           edepT_[indx]->Fill(edep);
           edepT = 0;
         }
         nhit[indx]++;
         if (indx == 0) {
           bool ok = false;
           for (unsigned int j = 0; j < ebID.size(); j++) {
             if (id_ == ebID[j]) {
               ebEtow[j] += edep;
               ebEtowT[j] += edepT;
               ok = true;
               break;
             }
           }
           if (!ok) {
             ebID.push_back(id_);
             ebEtow.push_back(edep);
             ebEtowT.push_back(edepT);
           }
         } else if (indx == 1) {
           bool ok = false;
           for (unsigned int j = 0; j < hbID.size(); j++) {
             if (id_ == hbID[j]) {
               hbEtow[j] += edep;
               hbEtowT[j] += edepT;
               ok = true;
               break;
             }
           }
           if (!ok) {
             hbID.push_back(id_);
             hbEtow.push_back(edep);
             hbEtowT.push_back(edepT);
           }
         } else {
           bool ok = false;
           for (unsigned int j = 0; j < hoID.size(); j++) {
             if (id_ == hoID[j]) {
               hoEtow[j] += edep;
               hoEtowT[j] += edepT;
               ok = true;
               break;
             }
           }
           if (!ok) {
             hoID.push_back(id_);
             hoEtow.push_back(edep);
             hoEtowT.push_back(edepT);
           }
         }
       }
     }
   }
 
   // Now for towers and total energy deposits
   for (int k = 0; k < 3; k++) {
     hit_[k]->Fill(double(nhit[k]));
     edepZon_[k]->Fill(etot[k]);
     edepZonT_[k]->Fill(etotT[k]);
   }
   hitTow_[0]->Fill(double(ebEtow.size()));
   for (unsigned int i = 0; i < ebEtow.size(); i++) {
     edepTW_[0]->Fill(ebEtow[i]);
     edepTWT_[0]->Fill(ebEtowT[i]);
   }
   hitTow_[1]->Fill(double(hbEtow.size()));
   for (unsigned int i = 0; i < hbEtow.size(); i++) {
     edepTW_[1]->Fill(hbEtow[i]);
     edepTWT_[1]->Fill(hbEtowT[i]);
   }
   hitTow_[2]->Fill(double(hoEtow.size()));
   for (unsigned int i = 0; i < hoEtow.size(); i++) {
     edepTW_[2]->Fill(hoEtow[i]);
     edepTWT_[2]->Fill(hoEtowT[i]);
   }
   double eEB = scaleEB_ * etot[0];
   double eEBHB = eEB + scaleHB_ * etot[1];
   double eEBHBHO = eEBHB + scaleHB_ * scaleHO_ * etot[2];
   eEB_->Fill(eEB);
   eEBHB_->Fill(eEBHB);
   eEBHBHO_->Fill(eEBHBHO);
   double eEBT = scaleEB_ * etotT[0];
   double eEBHBT = eEBT + scaleHB_ * etotT[1];
   double eEBHBHOT = eEBHBT + scaleHB_ * scaleHO_ * etotT[2];
   eEBT_->Fill(eEBT);
   eEBHBT_->Fill(eEBHBT);
   eEBHBHOT_->Fill(eEBHBHOT);
   eHO1_->Fill(etaInc, eHO17 + eHO18);
   eHO2_->Fill(etaInc, phiInc, eHO17 + eHO18);
   eHO17_->Fill(etaInc, eHO17);
   eHO18_->Fill(etaInc, eHO18);
   eHO1T_->Fill(etaInc, eHO17T + eHO18T);
   eHO2T_->Fill(etaInc, phiInc, eHO17T + eHO18T);
   eHO17T_->Fill(etaInc, eHO17T);
   eHO18T_->Fill(etaInc, eHO18T);
   int nHO = 0, nHOT = 0;
   if (eHO17 > 0)
     nHO++;
   if (eHO17T > 0)
     nHOT++;
   if (eHO18 > 0)
     nHO++;
   if (eHO18T > 0)
     nHOT++;
   nHO1_->Fill(etaInc, (double)(nHO));
   nHO2_->Fill(etaInc, phiInc, (double)(nHO));
   nHO1T_->Fill(etaInc, (double)(nHOT));
   nHO2T_->Fill(etaInc, phiInc, (double)(nHOT));
   int ieta = 15;
   for (int k = 0; k < 15; ++k) {
     if (std::abs(etaInc) < 0.087 * (k + 1)) {
       ieta = k;
       break;
     }
   }
   if (ieta >= 0 && ieta < 15) {
     eHOE_[ieta]->Fill(eHO17 + eHO18);
     eHOE17_[ieta]->Fill(eHO17);
     eHOE18_[ieta]->Fill(eHO18);
     eHOET_[ieta]->Fill(eHO17T + eHO18T);
     eHOE17T_[ieta]->Fill(eHO17T);
     eHOE18T_[ieta]->Fill(eHO18T);
     eHOEta_[ieta]->Fill(eHOE17[ieta] + eHOE18[ieta]);
     eHOEta17_[ieta]->Fill(eHOE17[ieta]);
     eHOEta18_[ieta]->Fill(eHOE18[ieta]);
     nHOE1_[ieta]->Fill((double)(nHO));
     nHOE2_[ieta]->Fill(phiInc, (double)(nHO));
     nHOE1T_[ieta]->Fill((double)(nHOT));
     nHOE2T_[ieta]->Fill(phiInc, (double)(nHOT));
     int nHOE = 0, nHOET = 0;
     if (eHOE17[ieta] > 0)
       nHOE++;
     if (eHOE17T[ieta] > 0)
       nHOET++;
     if (eHOE18[ieta] > 0)
       nHOE++;
     if (eHOE18T[ieta] > 0)
       nHOET++;
     nHOEta1_[ieta]->Fill((double)(nHOE));
     nHOEta2_[ieta]->Fill(phiInc, (double)(nHOE));
     nHOEta1T_[ieta]->Fill((double)(nHOET));
     nHOEta2T_[ieta]->Fill(phiInc, (double)(nHOET));
   }
 
   edm::LogVerbatim("HitStudy") << "HOSimHitStudy::analyzeHits: Hits in EB " << nhit[0] << " in " << ebEtow.size()
                                << " towers with total E " << etot[0] << "|" << etotT[0]
                                << "\n                            Hits in HB " << nhit[1] << " in " << hbEtow.size()
                                << " towers with total E " << etot[1] << "|" << etotT[1]
                                << "\n                            Hits in HO " << nhit[2] << " in " << hoEtow.size()
                                << " towers with total E " << etot[2] << "|" << etotT[2]
                                << "\n                            Energy in EB " << eEB << "|" << eEBT << " with HB "
                                << eEBHB << "|" << eEBHBT << " and with HO " << eEBHBHO << "|" << eEBHBHOT
                                << "\n                            E in HO layers " << eHO17 << "|" << eHO17T << " "
                                << eHO18 << "|" << eHO18T << " number of HO hits " << nHO << "|" << nHOT;
 
   if (eEBHBHO > 0.75 * maxEnergy_ && print_) {
     edm::LogVerbatim("HitStudy") << "HOSimHitStudy::Event with excessive energy: EB = " << eEB << " EB+HB = " << eEBHB
                                  << " EB+HB+HO = " << eEBHBHO;
     const std::string Dets[3] = {"EB", "HB", "HO"};
     for (int k = 0; k < 2; k++) {
       int nHit;
       if (k == 0) {
         nHit = ecalHits.size();
       } else {
         nHit = hcalHits.size();
       }
       for (int i = 0; i < nHit; i++) {
         double edep, time;
         int indx = -1;
         unsigned int id_;
         int ieta, iphi, depth = 0;
         if (k == 0) {
           indx = 0;
           edep = ecalHits[i].energy();
           time = ecalHits[i].time();
           id_ = ecalHits[i].id();
           EBDetId id = EBDetId(id_);
           ieta = id.ieta();
           iphi = id.iphi();
         } else {
           indx = -1;
           edep = hcalHits[i].energy();
           time = hcalHits[i].time();
           id_ = hcalHits[i].id();
           HcalDetId id = HcalDetId(id_);
           int subdet = id.subdet();
           if (subdet == static_cast<int>(HcalBarrel))
             indx = 1;
           else if (subdet == static_cast<int>(HcalOuter))
             indx = 2;
           ieta = id.ieta();
           iphi = id.iphi();
           depth = id.depth();
         }
         if (indx >= 0) {
           edm::LogVerbatim("HitStudy") << "HOSimHitStudy::" << Dets[indx] << " " << i << std::hex << id_ << std::dec
                                        << " (" << ieta << "|" << iphi << "|" << depth << ") " << std::setw(8) << edep
                                        << " " << std::setw(8) << time;
         }
       }
     }
   }
 }

◆ beginJob()

void HOSimHitStudy::beginJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 119 of file HOSimHitStudy.cc.

References edep_, edepT_, edepTW_, edepTWT_, edepZon_, edepZonT_, eEB_, eEBHB_, eEBHBHO_, eEBHBHOT_, eEBHBT_, eEBT_, eHO17_, eHO17T_, eHO18_, eHO18T_, eHO1_, eHO1T_, eHO2_, eHO2T_, eHOE17_, eHOE17T_, eHOE18_, eHOE18T_, eHOE_, eHOET_, eHOEta17_, eHOEta17T_, eHOEta18_, eHOEta18T_, eHOEta_, eHOEtaT_, eneInc_, etaInc_, hit_, hitTow_, mps_fire::i, createfilelist::int, maxEnergy_, Skims_PA_cff::name, nHO1_, nHO1T_, nHO2_, nHO2T_, nHOE1_, nHOE1T_, nHOE2_, nHOE2T_, nHOEta1_, nHOEta1T_, nHOEta2_, nHOEta2T_, phiInc_, AlCaHLTBitMon_QueryRunRegistry::string, tcut_, compare::tfile, time_, runGCPTkAlMap::title, and L1TOccupancyClient_cfi::ymax.

                              {
   edm::Service<TFileService> tfile;
 
   if (!tfile.isAvailable())
     throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                       << "please add it to config file";
   std::string dets[3] = {"EB", "HB", "HO"};
   char name[60], title[100];
   double ymax = maxEnergy_;
   sprintf(title, "Incident Energy (GeV)");
   eneInc_ = tfile->make<TH1F>("EneInc", title, 1000, 0., ymax);
   eneInc_->GetXaxis()->SetTitle(title);
   eneInc_->GetYaxis()->SetTitle("Events");
   sprintf(title, "Incident #eta");
   etaInc_ = tfile->make<TH1F>("EtaInc", title, 200, -5., 5.);
   etaInc_->GetXaxis()->SetTitle(title);
   etaInc_->GetYaxis()->SetTitle("Events");
   sprintf(title, "Incident #phi");
   phiInc_ = tfile->make<TH1F>("PhiInc", title, 200, -3.1415926, 3.1415926);
   phiInc_->GetXaxis()->SetTitle(title);
   phiInc_->GetYaxis()->SetTitle("Events");
   int itcut = (int)(tcut_);
   for (int i = 0; i < 3; i++) {
     sprintf(name, "Hit%d", i);
     sprintf(title, "Number of hits in %s", dets[i].c_str());
     hit_[i] = tfile->make<TH1F>(name, title, 1000, 0., 20000.);
     hit_[i]->GetXaxis()->SetTitle(title);
     hit_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "Time%d", i);
     sprintf(title, "Time of the hit (ns) in %s", dets[i].c_str());
     time_[i] = tfile->make<TH1F>(name, title, 1200, 0., 1200.);
     time_[i]->GetXaxis()->SetTitle(title);
     time_[i]->GetYaxis()->SetTitle("Hits");
     if (i > 0)
       ymax = 1.0;
     else
       ymax = 50.0;
     sprintf(name, "Edep%d", i);
     sprintf(title, "Energy deposit (GeV) in %s", dets[i].c_str());
     edep_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     edep_[i]->GetXaxis()->SetTitle(title);
     edep_[i]->GetYaxis()->SetTitle("Hits");
     sprintf(name, "EdepT%d", i);
     sprintf(title, "Energy deposit (GeV) in %s for t < %d ns", dets[i].c_str(), itcut);
     edepT_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     edepT_[i]->GetXaxis()->SetTitle(title);
     edepT_[i]->GetYaxis()->SetTitle("Hits");
     sprintf(name, "HitTow%d", i);
     sprintf(title, "Number of towers with hits in %s", dets[i].c_str());
     hitTow_[i] = tfile->make<TH1F>(name, title, 1000, 0., 20000.);
     hitTow_[i]->GetXaxis()->SetTitle(title);
     hitTow_[i]->GetYaxis()->SetTitle("Events");
     if (i > 0)
       ymax = 0.05 * maxEnergy_;
     else
       ymax = maxEnergy_;
     sprintf(name, "EdepTW%d", i);
     sprintf(title, "Energy deposit (GeV) in %s Tower", dets[i].c_str());
     edepTW_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     edepTW_[i]->GetXaxis()->SetTitle(title);
     edepTW_[i]->GetYaxis()->SetTitle("Towers");
     sprintf(name, "EdepTWT%d", i);
     sprintf(title, "Energy deposit (GeV) in %s Tower for t < %d ns", dets[i].c_str(), itcut);
     edepTWT_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     edepTWT_[i]->GetXaxis()->SetTitle(title);
     edepTWT_[i]->GetYaxis()->SetTitle("Towers");
     sprintf(name, "EdepZone%d", i);
     sprintf(title, "Energy deposit (GeV) in %s", dets[i].c_str());
     edepZon_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     edepZon_[i]->GetXaxis()->SetTitle(title);
     edepZon_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EdepZoneT%d", i);
     sprintf(title, "Energy deposit (GeV) in %s for t < %d ns", dets[i].c_str(), itcut);
     edepZonT_[i] = tfile->make<TH1F>(name, title, 5000, 0., ymax);
     edepZonT_[i]->GetXaxis()->SetTitle(title);
     edepZonT_[i]->GetYaxis()->SetTitle("Events");
   }
   sprintf(title, "Energy Measured in EB (GeV)");
   eEB_ = tfile->make<TH1F>("EEB", title, 5000, 0., maxEnergy_);
   eEB_->GetXaxis()->SetTitle(title);
   eEB_->GetYaxis()->SetTitle("Events");
   sprintf(title, "Energy Measured in EB+HB (GeV)");
   eEBHB_ = tfile->make<TH1F>("EEBHB", title, 5000, 0., maxEnergy_);
   eEBHB_->GetXaxis()->SetTitle(title);
   eEBHB_->GetYaxis()->SetTitle("Events");
   sprintf(title, "Energy Measured in EB+HB+HO (GeV)");
   eEBHBHO_ = tfile->make<TH1F>("EEBHBHO", title, 5000, 0., maxEnergy_);
   eEBHBHO_->GetXaxis()->SetTitle(title);
   eEBHBHO_->GetYaxis()->SetTitle("Events");
   sprintf(title, "Energy Measured in EB (GeV) for t < %d ns", itcut);
   eEBT_ = tfile->make<TH1F>("EEBT", title, 5000, 0., maxEnergy_);
   eEBT_->GetXaxis()->SetTitle(title);
   eEBT_->GetYaxis()->SetTitle("Events");
   sprintf(title, "Energy Measured in EB+HB (GeV) for t < %d ns", itcut);
   eEBHBT_ = tfile->make<TH1F>("EEBHBT", title, 5000, 0., maxEnergy_);
   eEBHBT_->GetXaxis()->SetTitle(title);
   eEBHBT_->GetYaxis()->SetTitle("Events");
   sprintf(title, "Energy Measured in EB+HB+HO (GeV) for t < %d ns", itcut);
   eEBHBHOT_ = tfile->make<TH1F>("EEBHBHOT", title, 5000, 0., maxEnergy_);
   eEBHBHOT_->GetXaxis()->SetTitle(title);
   eEBHBHOT_->GetYaxis()->SetTitle("Events");
   sprintf(title, "SimHit energy in HO");
   eHO1_ = tfile->make<TProfile>("EHO1", title, 30, -1.305, 1.305);
   eHO1_->GetXaxis()->SetTitle(title);
   eHO1_->GetYaxis()->SetTitle("Events");
   eHO2_ = tfile->make<TProfile2D>("EHO2", title, 30, -1.305, 1.305, 72, -3.1415926, 3.1415926);
   eHO2_->GetXaxis()->SetTitle(title);
   eHO2_->GetYaxis()->SetTitle("Events");
   sprintf(title, "SimHit energy in HO Layer 17");
   eHO17_ = tfile->make<TProfile>("EHO17", title, 30, -1.305, 1.305);
   eHO17_->GetXaxis()->SetTitle(title);
   eHO17_->GetYaxis()->SetTitle("Events");
   sprintf(title, "SimHit energy in HO Layer 18");
   eHO18_ = tfile->make<TProfile>("EHO18", title, 30, -1.305, 1.305);
   eHO18_->GetXaxis()->SetTitle(title);
   eHO18_->GetYaxis()->SetTitle("Events");
   sprintf(title, "SimHit energy in HO for t < %d ns", itcut);
   eHO1T_ = tfile->make<TProfile>("EHO1T", title, 30, -1.305, 1.305);
   eHO1T_->GetXaxis()->SetTitle(title);
   eHO1T_->GetYaxis()->SetTitle("Events");
   eHO2T_ = tfile->make<TProfile2D>("EHO2T", title, 30, -1.305, 1.305, 72, -3.1415926, 3.1415926);
   eHO2T_->GetXaxis()->SetTitle(title);
   eHO2T_->GetYaxis()->SetTitle("Events");
   sprintf(title, "SimHit energy in HO Layer 17 for t < %d ns", itcut);
   eHO17T_ = tfile->make<TProfile>("EHO17T", title, 30, -1.305, 1.305);
   eHO17T_->GetXaxis()->SetTitle(title);
   eHO17T_->GetYaxis()->SetTitle("Events");
   sprintf(title, "SimHit energy in HO Layer 18 for t < %d ns", itcut);
   eHO18T_ = tfile->make<TProfile>("EHO18T", title, 30, -1.305, 1.305);
   eHO18T_->GetXaxis()->SetTitle(title);
   eHO18T_->GetYaxis()->SetTitle("Events");
   sprintf(title, "Number of layers hit in HO");
   nHO1_ = tfile->make<TProfile>("NHO1", title, 30, -1.305, 1.305);
   nHO1_->GetXaxis()->SetTitle(title);
   nHO1_->GetYaxis()->SetTitle("Events");
   nHO2_ = tfile->make<TProfile2D>("NHO2", title, 30, -1.305, 1.305, 72, -3.1415926, 3.1415926);
   nHO2_->GetXaxis()->SetTitle(title);
   nHO2_->GetYaxis()->SetTitle("Events");
   sprintf(title, "Number of layers hit in HO for t < %d ns", itcut);
   nHO1T_ = tfile->make<TProfile>("NHO1T", title, 30, -1.305, 1.305);
   nHO1T_->GetXaxis()->SetTitle(title);
   nHO1T_->GetYaxis()->SetTitle("Events");
   nHO2T_ = tfile->make<TProfile2D>("NHO2T", title, 30, -1.305, 1.305, 72, -3.1415926, 3.1415926);
   nHO2T_->GetXaxis()->SetTitle(title);
   nHO2T_->GetYaxis()->SetTitle("Events");
   for (int i = 0; i < 15; i++) {
     sprintf(name, "EHOE%d", i + 1);
     sprintf(title, "SimHit energy in HO (Beam in #eta=%d bin)", i + 1);
     eHOE_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOE_[i]->GetXaxis()->SetTitle(title);
     eHOE_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOE17%d", i + 1);
     sprintf(title, "SimHit energy in Layer 17 (Beam in #eta=%d bin)", i + 1);
     eHOE17_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOE17_[i]->GetXaxis()->SetTitle(title);
     eHOE17_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOE18%d", i + 1);
     sprintf(title, "SimHit energy in Layer 18 (Beam in #eta=%d bin)", i + 1);
     eHOE18_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOE18_[i]->GetXaxis()->SetTitle(title);
     eHOE18_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOE%dT", i + 1);
     sprintf(title, "SimHit energy in HO (Beam in #eta=%d bin, t < %d ns)", i + 1, itcut);
     eHOET_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOET_[i]->GetXaxis()->SetTitle(title);
     eHOET_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOE17%dT", i + 1);
     sprintf(title, "SimHit energy in Layer 17 (Beam in #eta=%d bin, t < %d ns)", i + 1, itcut);
     eHOE17T_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOE17T_[i]->GetXaxis()->SetTitle(title);
     eHOE17T_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOE18%dT", i + 1);
     sprintf(title, "SimHit energy in Layer 18 (Beam in #eta=%d bin, t < %d ns)", i + 1, itcut);
     eHOE18T_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOE18T_[i]->GetXaxis()->SetTitle(title);
     eHOE18T_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOEta%d", i + 1);
     sprintf(title, "SimHit energy in HO #eta bin %d (Beam in #eta=%d bin)", i + 1, i + 1);
     eHOEta_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOEta_[i]->GetXaxis()->SetTitle(title);
     eHOEta_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOEta17%d", i + 1);
     sprintf(title, "SimHit energy in Layer 17 #eta bin %d (Beam in #eta=%d bin)", i + 1, i + 1);
     eHOEta17_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOEta17_[i]->GetXaxis()->SetTitle(title);
     eHOEta17_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOEta18%d", i + 1);
     sprintf(title, "SimHit energy in Layer 18 #eta bin %d (Beam in #eta=%d bin)", i + 1, i + 1);
     eHOEta18_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOEta18_[i]->GetXaxis()->SetTitle(title);
     eHOEta18_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOEta%dT", i + 1);
     sprintf(title, "SimHit energy in HO #eta bin %d (Beam in #eta=%d bin, t < %d ns)", i + 1, i + 1, itcut);
     eHOEtaT_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOEtaT_[i]->GetXaxis()->SetTitle(title);
     eHOEtaT_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOEta17%dT", i + 1);
     sprintf(title, "SimHit energy in Layer 17 #eta bin %d (Beam in #eta=%d bin, t < %d ns)", i + 1, i + 1, itcut);
     eHOEta17T_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOEta17T_[i]->GetXaxis()->SetTitle(title);
     eHOEta17T_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "EHOEta18%dT", i + 1);
     sprintf(title, "SimHit energy in Layer 18 #eta bin %d (Beam in #eta=%d bin, t < %d ns)", i + 1, i + 1, itcut);
     eHOEta18T_[i] = tfile->make<TH1F>(name, title, 1000, 0., 0.25);
     eHOEta18T_[i]->GetXaxis()->SetTitle(title);
     eHOEta18T_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "NHOE1%d", i + 1);
     sprintf(title, "Number of layers hit in HO  (Beam in #eta=%d bin)", i + 1);
     nHOE1_[i] = tfile->make<TH1F>(name, title, 20, 0, 20.);
     nHOE1_[i]->GetXaxis()->SetTitle(title);
     nHOE1_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "NHOE2%d", i + 1);
     nHOE2_[i] = tfile->make<TProfile>(name, title, 72, -3.1415926, 3.1415926);
     nHOE2_[i]->GetXaxis()->SetTitle(title);
     nHOE2_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "NHOE1%dT", i + 1);
     sprintf(title, "Number of layers hit in HO (Beam in #eta=%d bin, t < %d ns)", i + 1, itcut);
     nHOE1T_[i] = tfile->make<TH1F>(name, title, 20, 0, 20.);
     nHOE1T_[i]->GetXaxis()->SetTitle(title);
     nHOE1T_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "NHOE2%dT", i + 1);
     nHOE2T_[i] = tfile->make<TProfile>(name, title, 72, -3.1415926, 3.1415926);
     nHOE2T_[i]->GetXaxis()->SetTitle(title);
     nHOE2T_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "NHOEta1%d", i + 1);
     sprintf(title, "Number of layers hit in HO #eta bin %d (Beam in #eta=%d bin)", i + 1, i + 1);
     nHOEta1_[i] = tfile->make<TH1F>(name, title, 20, 0, 20.);
     nHOEta1_[i]->GetXaxis()->SetTitle(title);
     nHOEta1_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "NHOEta2%d", i + 1);
     nHOEta2_[i] = tfile->make<TProfile>(name, title, 72, -3.1415926, 3.1415926);
     nHOEta2_[i]->GetXaxis()->SetTitle(title);
     nHOEta2_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "NHOEta1%dT", i + 1);
     sprintf(title, "Number of layers hit in HO #eta bin %d (Beam in #eta=%d bin, t < %d ns)", i + 1, i + 1, itcut);
     nHOEta1T_[i] = tfile->make<TH1F>(name, title, 20, 0, 20.);
     nHOEta1T_[i]->GetXaxis()->SetTitle(title);
     nHOEta1T_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "NHOEta2%dT", i + 1);
     nHOEta2T_[i] = tfile->make<TProfile>(name, title, 72, -3.1415926, 3.1415926);
     nHOEta2T_[i]->GetXaxis()->SetTitle(title);
     nHOEta2T_[i]->GetYaxis()->SetTitle("Events");
   }
 }

◆ beginRun()

void HOSimHitStudy::beginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

inlineoverrideprotected

Definition at line 49 of file HOSimHitStudy.cc.

49 {}

◆ endJob()

void HOSimHitStudy::endJob ( void )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 50 of file HOSimHitStudy.cc.

50 {}

◆ endRun()

void HOSimHitStudy::endRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

inlineoverrideprotected

Definition at line 51 of file HOSimHitStudy.cc.

51 {}

◆ fillDescriptions()

void HOSimHitStudy::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 103 of file HOSimHitStudy.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, SummaryClient_cfi::labels, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                {
   edm::ParameterSetDescription desc;
   std::vector<std::string> labels = {"EcalHitsEB", "HcalHits"};
   desc.addUntracked<std::string>("SourceLabel", "generatorSmeared");
   desc.addUntracked<std::string>("ModuleLabel", "g4SimHits");
   desc.addUntracked<std::vector<std::string>>("HitCollection", labels);
   desc.addUntracked<double>("MaxEnergy", 50.0);
   desc.addUntracked<double>("ScaleEB", 1.02);
   desc.addUntracked<double>("ScaleHB", 104.4);
   desc.addUntracked<double>("ScaleHO", 2.33);
   desc.addUntracked<double>("TimeCut", 100.0);
   desc.addUntracked<bool>("PrintExcessEnergy", true);
   desc.addUntracked<bool>("TestNumbering", false);
   descriptions.add("HOSimHitStudy", desc);
 }

Member Data Documentation

◆ ecalHits

std::vector<PCaloHit> HOSimHitStudy::ecalHits

private

Definition at line 64 of file HOSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

◆ edep_

TH1F* HOSimHitStudy::edep_[3]

private

Definition at line 66 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ edepT_

TH1F* HOSimHitStudy::edepT_[3]

private

Definition at line 67 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ edepTW_

TH1F * HOSimHitStudy::edepTW_[3]

private

Definition at line 65 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ edepTWT_

TH1F * HOSimHitStudy::edepTWT_[3]

private

Definition at line 65 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ edepZon_

TH1F* HOSimHitStudy::edepZon_[3]

private

Definition at line 68 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ edepZonT_

TH1F * HOSimHitStudy::edepZonT_[3]

private

Definition at line 68 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eEB_

TH1F * HOSimHitStudy::eEB_

private

Definition at line 67 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eEBHB_

TH1F * HOSimHitStudy::eEBHB_

private

Definition at line 67 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eEBHBHO_

TH1F * HOSimHitStudy::eEBHBHO_

private

Definition at line 67 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eEBHBHOT_

TH1F * HOSimHitStudy::eEBHBHOT_

private

Definition at line 67 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eEBHBT_

TH1F * HOSimHitStudy::eEBHBT_

private

Definition at line 68 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eEBT_

TH1F * HOSimHitStudy::eEBT_

private

Definition at line 68 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHO17_

TProfile * HOSimHitStudy::eHO17_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHO17T_

TProfile * HOSimHitStudy::eHO17T_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHO18_

TProfile * HOSimHitStudy::eHO18_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHO18T_

TProfile * HOSimHitStudy::eHO18T_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHO1_

TProfile* HOSimHitStudy::eHO1_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHO1T_

TProfile * HOSimHitStudy::eHO1T_

private

Definition at line 74 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHO2_

TProfile2D* HOSimHitStudy::eHO2_

private

Definition at line 77 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHO2T_

TProfile2D * HOSimHitStudy::eHO2T_

private

Definition at line 77 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHOE17_

TH1F* HOSimHitStudy::eHOE17_[15]

private

Definition at line 69 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHOE17T_

TH1F* HOSimHitStudy::eHOE17T_[15]

private

Definition at line 70 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHOE18_

TH1F * HOSimHitStudy::eHOE18_[15]

private

Definition at line 69 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHOE18T_

TH1F * HOSimHitStudy::eHOE18T_[15]

private

Definition at line 70 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHOE_

TH1F * HOSimHitStudy::eHOE_[15]

private

Definition at line 69 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHOET_

TH1F * HOSimHitStudy::eHOET_[15]

private

Definition at line 70 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHOEta17_

TH1F* HOSimHitStudy::eHOEta17_[15]

private

Definition at line 71 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHOEta17T_

TH1F* HOSimHitStudy::eHOEta17T_[15]

private

Definition at line 72 of file HOSimHitStudy.cc.

Referenced by beginJob().

◆ eHOEta18_

TH1F * HOSimHitStudy::eHOEta18_[15]

private

Definition at line 71 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHOEta18T_

TH1F * HOSimHitStudy::eHOEta18T_[15]

private

Definition at line 72 of file HOSimHitStudy.cc.

Referenced by beginJob().

◆ eHOEta_

TH1F * HOSimHitStudy::eHOEta_[15]

private

Definition at line 71 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ eHOEtaT_

TH1F * HOSimHitStudy::eHOEtaT_[15]

private

Definition at line 72 of file HOSimHitStudy.cc.

Referenced by beginJob().

◆ eInc

double HOSimHitStudy::eInc

private

Definition at line 78 of file HOSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

◆ eneInc_

TH1F * HOSimHitStudy::eneInc_

private

Definition at line 66 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ etaInc

double HOSimHitStudy::etaInc

private

Definition at line 78 of file HOSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

◆ etaInc_

TH1F * HOSimHitStudy::etaInc_

private

Definition at line 66 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ g4Label_

const std::string HOSimHitStudy::g4Label_

private

Definition at line 57 of file HOSimHitStudy.cc.

Referenced by HOSimHitStudy().

◆ hcalHits

std::vector<PCaloHit> HOSimHitStudy::hcalHits

private

Definition at line 64 of file HOSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

◆ hit_

TH1F* HOSimHitStudy::hit_[3]

private

Definition at line 65 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ hitLab_

const std::vector<std::string> HOSimHitStudy::hitLab_

private

Definition at line 58 of file HOSimHitStudy.cc.

Referenced by analyze(), and HOSimHitStudy().

◆ hitTow_

TH1F * HOSimHitStudy::hitTow_[3]

private

Definition at line 66 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ maxEnergy_

const double HOSimHitStudy::maxEnergy_

private

Definition at line 59 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHO1_

TProfile* HOSimHitStudy::nHO1_

private

Definition at line 75 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHO1T_

TProfile * HOSimHitStudy::nHO1T_

private

Definition at line 75 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHO2_

TProfile2D * HOSimHitStudy::nHO2_

private

Definition at line 77 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHO2T_

TProfile2D * HOSimHitStudy::nHO2T_

private

Definition at line 77 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHOE1_

TH1F* HOSimHitStudy::nHOE1_[15]

private

Definition at line 73 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHOE1T_

TH1F * HOSimHitStudy::nHOE1T_[15]

private

Definition at line 73 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHOE2_

TProfile * HOSimHitStudy::nHOE2_[15]

private

Definition at line 75 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHOE2T_

TProfile * HOSimHitStudy::nHOE2T_[15]

private

Definition at line 75 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHOEta1_

TH1F * HOSimHitStudy::nHOEta1_[15]

private

Definition at line 73 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHOEta1T_

TH1F * HOSimHitStudy::nHOEta1T_[15]

private

Definition at line 73 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHOEta2_

TProfile* HOSimHitStudy::nHOEta2_[15]

private

Definition at line 76 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ nHOEta2T_

TProfile * HOSimHitStudy::nHOEta2T_[15]

private

Definition at line 76 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ phiInc

double HOSimHitStudy::phiInc

private

Definition at line 78 of file HOSimHitStudy.cc.

Referenced by analyze(), and analyzeHits().

◆ phiInc_

TH1F * HOSimHitStudy::phiInc_

private

Definition at line 66 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ print_

const bool HOSimHitStudy::print_

private

Definition at line 61 of file HOSimHitStudy.cc.

Referenced by analyzeHits().

◆ scaleEB_

const double HOSimHitStudy::scaleEB_

private

Definition at line 59 of file HOSimHitStudy.cc.

Referenced by analyzeHits().

◆ scaleHB_

const double HOSimHitStudy::scaleHB_

private

Definition at line 59 of file HOSimHitStudy.cc.

Referenced by analyzeHits().

◆ scaleHO_

const double HOSimHitStudy::scaleHO_

private

Definition at line 59 of file HOSimHitStudy.cc.

Referenced by analyzeHits().

◆ scheme_

const bool HOSimHitStudy::scheme_

private

Definition at line 61 of file HOSimHitStudy.cc.

Referenced by analyzeHits().

◆ tcut_

const double HOSimHitStudy::tcut_

private

Definition at line 60 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ time_

TH1F * HOSimHitStudy::time_[3]

private

Definition at line 65 of file HOSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ tok_evt_

const edm::EDGetTokenT<edm::HepMCProduct> HOSimHitStudy::tok_evt_

private

Definition at line 62 of file HOSimHitStudy.cc.

Referenced by analyze().

◆ toks_calo_

const std::vector<edm::EDGetTokenT<edm::PCaloHitContainer> > HOSimHitStudy::toks_calo_

private

Definition at line 63 of file HOSimHitStudy.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HOSimHitStudy()

◆ ~HOSimHitStudy()

Member Function Documentation

◆ analyze()

◆ analyzeHits()

◆ beginJob()

◆ beginRun()

◆ endJob()

◆ endRun()

◆ fillDescriptions()

Member Data Documentation

◆ ecalHits

◆ edep_

◆ edepT_

◆ edepTW_

◆ edepTWT_

◆ edepZon_

◆ edepZonT_

◆ eEB_

◆ eEBHB_

◆ eEBHBHO_

◆ eEBHBHOT_

◆ eEBHBT_

◆ eEBT_

◆ eHO17_

◆ eHO17T_

◆ eHO18_

◆ eHO18T_

◆ eHO1_

◆ eHO1T_

◆ eHO2_

◆ eHO2T_

◆ eHOE17_

◆ eHOE17T_

◆ eHOE18_

◆ eHOE18T_

◆ eHOE_

◆ eHOET_

◆ eHOEta17_

◆ eHOEta17T_

◆ eHOEta18_

◆ eHOEta18T_

◆ eHOEta_

◆ eHOEtaT_

◆ eInc

◆ eneInc_

◆ etaInc

◆ etaInc_

◆ g4Label_

◆ hcalHits

◆ hit_

◆ hitLab_

◆ hitTow_

◆ maxEnergy_

◆ nHO1_

◆ nHO1T_

◆ nHO2_

◆ nHO2T_

◆ nHOE1_

◆ nHOE1T_

◆ nHOE2_

◆ nHOE2T_

◆ nHOEta1_

◆ nHOEta1T_

◆ nHOEta2_

◆ nHOEta2T_

◆ phiInc

◆ phiInc_

◆ print_

◆ scaleEB_

◆ scaleHB_

◆ scaleHO_

◆ scheme_