#include <Calibration/HcalCalibAlgos/test/RecAnalyzerMinbias.cc>

Inheritance diagram for RecAnalyzerMinbias:

Classes
struct	myInfo

Public Member Functions
	RecAnalyzerMinbias (const edm::ParameterSet &)

	~RecAnalyzerMinbias () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () 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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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)

Private Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override

void	analyzeHcal (const HBHERecHitCollection &, const HFRecHitCollection &, int, bool, double)

void	beginJob () override

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

void	endJob () override

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

Private Attributes
int	cells

std::map< DetId, double >	corrFactor_

int	depth

double	eHighHB_

double	eHighHE_

double	eHighHF_

double	eLowHB_

double	eLowHE_

double	eLowHF_

double	eMin_

bool	extraHist_

bool	fillHist_

edm::Service< TFileService >	fs_

TH1D *	h_ [4]

TH1D *	h_AmplitudeHB_

TH1D *	h_AmplitudeHBtest_

TH1D *	h_AmplitudeHE_

TH1D *	h_AmplitudeHEtest_

TH1D *	h_AmplitudeHF_

TH1D *	h_AmplitudeHFtest_

TH2D *	hb_

TH2D *	hbhe_

TProfile *	hbherun_

int	HBHEsize

TProfile *	hbrun_

std::vector< unsigned int >	hcalID_

TH2D *	he_

TProfile *	herun_

TH2D *	hf_

TProfile *	hfrun_

int	HFsize

std::map< HcalDetId, TH1D * >	histHC_

std::vector< TH1D * >	histo_

int	ieta

bool	ignoreL1_

bool	init_

int	iphi

float	mom0_MB

float	mom1_MB

float	mom2_MB

float	mom3_MB

float	mom4_MB

std::map< std::pair< int, HcalDetId >, myInfo >	myMap_

int	mysubd

TTree *	myTree1_

TTree *	myTree_

bool	Noise_

double	rnnum_

double	rnnumber

int	runMax_

int	runMin_

bool	runNZS_

bool	theRecalib_

edm::EDGetTokenT< GenEventInfoProduct >	tok_ew_

edm::EDGetTokenT< L1GlobalTriggerReadoutRecord >	tok_gtRec_

edm::EDGetTokenT< HBHEDigiCollection >	tok_hbhedigi_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbherecoMB_

edm::EDGetTokenT< HFDigiCollection >	tok_hfdigi_

edm::EDGetTokenT< HFRecHitCollection >	tok_hfrecoMB_

edm::EDGetTokenT< L1GlobalTriggerObjectMapRecord >	tok_hltL1GtMap_

edm::EDGetTokenT< HODigiCollection >	tok_hodigi_

edm::EDGetTokenT< QIE10DigiCollection >	tok_qie10digi_

edm::EDGetTokenT< QIE11DigiCollection >	tok_qie11digi_

int	trigbit

std::vector< int >	trigbit_

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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)

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

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

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

Detailed Description

Description: Performs phi-symmetry studies of HB/HE/HF channels

Implementation: <Notes on="" implementation>="">

Definition at line 68 of file RecAnalyzerMinbias.cc.

Constructor & Destructor Documentation

◆ RecAnalyzerMinbias()

RecAnalyzerMinbias::RecAnalyzerMinbias ( const edm::ParameterSet & iConfig )

explicit

Definition at line 128 of file RecAnalyzerMinbias.cc.

                                                                      : init_(false) {
   usesResource("TFileService");
  
   // get name of output file with histogramms
   runNZS_ = iConfig.getParameter<bool>("runNZS");
   Noise_ = iConfig.getParameter<bool>("noise");
   eLowHB_ = iConfig.getParameter<double>("eLowHB");
   eHighHB_ = iConfig.getParameter<double>("eHighHB");
   eLowHE_ = iConfig.getParameter<double>("eLowHE");
   eHighHE_ = iConfig.getParameter<double>("eHighHE");
   eLowHF_ = iConfig.getParameter<double>("eLowHF");
   eHighHF_ = iConfig.getParameter<double>("eHighHF");
   eMin_ = iConfig.getUntrackedParameter<double>("eMin", 2.0);
   // The following run range is suited to study 2017 commissioning period
   runMin_ = iConfig.getUntrackedParameter<int>("RunMin", 308327);
   runMax_ = iConfig.getUntrackedParameter<int>("RunMax", 315250);
   trigbit_ = iConfig.getUntrackedParameter<std::vector<int>>("triggerBits");
   ignoreL1_ = iConfig.getUntrackedParameter<bool>("ignoreL1", false);
   std::string cfile = iConfig.getUntrackedParameter<std::string>("corrFile");
   fillHist_ = iConfig.getUntrackedParameter<bool>("fillHisto", false);
   extraHist_ = iConfig.getUntrackedParameter<bool>("extraHisto", false);
   std::vector<int> ieta = iConfig.getUntrackedParameter<std::vector<int>>("hcalIeta");
   std::vector<int> iphi = iConfig.getUntrackedParameter<std::vector<int>>("hcalIphi");
   std::vector<int> depth = iConfig.getUntrackedParameter<std::vector<int>>("hcalDepth");
  
   // get token names of modules, producing object collections
   tok_hbherecoMB_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInputMB"));
   tok_hfrecoMB_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("hfInputMB"));
   tok_hltL1GtMap_ = consumes<L1GlobalTriggerObjectMapRecord>(edm::InputTag("hltL1GtObjectMap"));
   tok_ew_ = consumes<GenEventInfoProduct>(edm::InputTag("generator"));
   tok_hbhedigi_ = consumes<HBHEDigiCollection>(iConfig.getParameter<edm::InputTag>("hcalDigiCollectionTag"));
   tok_qie11digi_ = consumes<QIE11DigiCollection>(iConfig.getParameter<edm::InputTag>("hcalDigiCollectionTag"));
   tok_hodigi_ = consumes<HODigiCollection>(iConfig.getParameter<edm::InputTag>("hcalDigiCollectionTag"));
   tok_hfdigi_ = consumes<HFDigiCollection>(iConfig.getParameter<edm::InputTag>("hcalDigiCollectionTag"));
   tok_qie10digi_ = consumes<QIE10DigiCollection>(iConfig.getParameter<edm::InputTag>("hcalDigiCollectionTag"));
   tok_gtRec_ = consumes<L1GlobalTriggerReadoutRecord>(edm::InputTag("gtDigisAlCaMB"));
  
   // Read correction factors
   std::ifstream infile(cfile.c_str());
   if (!infile.is_open()) {
     theRecalib_ = false;
     edm::LogWarning("RecAnalyzer") << "Cannot open '" << cfile << "' for the correction file";
   } else {
     unsigned int ndets(0), nrec(0);
     while (true) {
       unsigned int id;
       double cfac;
       infile >> id >> cfac;
       if (!infile.good())
         break;
       HcalDetId detId(id);
       nrec++;
       std::map<DetId, double>::iterator itr = corrFactor_.find(detId);
       if (itr == corrFactor_.end()) {
         corrFactor_[detId] = cfac;
         ndets++;
       }
     }
     infile.close();
     edm::LogInfo("RecAnalyzer") << "Reads " << nrec << " correction factors for " << ndets << " detIds";
     theRecalib_ = (ndets > 0);
   }
  
   edm::LogInfo("RecAnalyzer") << " Flags (ReCalib): " << theRecalib_ << " (IgnoreL1): " << ignoreL1_ << " (NZS) "
                               << runNZS_ << " and with " << ieta.size() << " detId for full histogram";
   edm::LogInfo("RecAnalyzer") << "Thresholds for HB " << eLowHB_ << ":" << eHighHB_ << "  for HE " << eLowHE_ << ":"
                               << eHighHE_ << "  for HF " << eLowHF_ << ":" << eHighHF_;
   for (unsigned int k = 0; k < ieta.size(); ++k) {
     HcalSubdetector subd =
         ((std::abs(ieta[k]) > 29) ? HcalForward
                                   : (std::abs(ieta[k]) > 16) ? HcalEndcap
                                                              : ((std::abs(ieta[k]) == 16) && (depth[k] == 3))
                                                                    ? HcalEndcap
                                                                    : (depth[k] == 4) ? HcalOuter : HcalBarrel);
     unsigned int id = (HcalDetId(subd, ieta[k], iphi[k], depth[k])).rawId();
     hcalID_.push_back(id);
     edm::LogInfo("RecAnalyzer") << "DetId[" << k << "] " << HcalDetId(id);
   }
   edm::LogInfo("RecAnalyzer") << "Select on " << trigbit_.size() << " L1 Trigger selection";
   for (unsigned int k = 0; k < trigbit_.size(); ++k)
     edm::LogInfo("RecAnalyzer") << "Bit[" << k << "] " << trigbit_[k];
 }

References funct::abs(), corrFactor_, depth, eHighHB_, eHighHE_, eHighHF_, eLowHB_, eLowHE_, eLowHF_, eMin_, extraHist_, fillHist_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HcalBarrel, HcalEndcap, HcalForward, hcalID_, HcalOuter, triggerObjects_cff::id, ieta, ignoreL1_, timingPdfMaker::infile, iphi, dqmdumpme::k, Noise_, runMax_, runMin_, runNZS_, AlCaHLTBitMon_QueryRunRegistry::string, theRecalib_, tok_ew_, tok_gtRec_, tok_hbhedigi_, tok_hbherecoMB_, tok_hfdigi_, tok_hfrecoMB_, tok_hltL1GtMap_, tok_hodigi_, tok_qie10digi_, tok_qie11digi_, and trigbit_.

◆ ~RecAnalyzerMinbias()

RecAnalyzerMinbias::~RecAnalyzerMinbias ( )

override

Definition at line 211 of file RecAnalyzerMinbias.cc.

211 {}

Member Function Documentation

◆ analyze()

void RecAnalyzerMinbias::analyze	(	edm::Event const &	iEvent,
		edm::EventSetup const &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 423 of file RecAnalyzerMinbias.cc.

                                                                                     {
   rnnum_ = (double)iEvent.run();
  
   if (extraHist_) {
     double amplitudefullHB(0), amplitudefullHE(0), amplitudefullHF(0);
     edm::Handle<HBHEDigiCollection> hbhedigi;
     iEvent.getByToken(tok_hbhedigi_, hbhedigi);
     if (hbhedigi.isValid()) {
       for (auto const& digi : *(hbhedigi.product())) {
         int nTS = digi.size();
         double amplitudefullTSs = 0.;
         if (digi.id().subdet() == HcalBarrel) {
           if (nTS <= 10) {
             for (int i = 0; i < nTS; i++)
               amplitudefullTSs += digi.sample(i).adc();
             h_AmplitudeHBtest_->Fill(amplitudefullTSs);
             amplitudefullHB += amplitudefullTSs;
           }
         }
         if (digi.id().subdet() == HcalEndcap) {
           if (nTS <= 10) {
             for (int i = 0; i < nTS; i++)
               amplitudefullTSs += digi.sample(i).adc();
             h_AmplitudeHEtest_->Fill(amplitudefullTSs);
             amplitudefullHE += amplitudefullTSs;
           }
         }
       }
     }
  
     edm::Handle<QIE11DigiCollection> qie11digi;
     iEvent.getByToken(tok_qie11digi_, qie11digi);
     if (qie11digi.isValid()) {
       for (QIE11DataFrame const digi : *(qie11digi.product())) {
         double amplitudefullTSs = 0.;
         if (HcalDetId(digi.id()).subdet() == HcalBarrel) {
           for (int i = 0; i < digi.samples(); i++)
             amplitudefullTSs += digi[i].adc();
           h_AmplitudeHBtest_->Fill(amplitudefullTSs);
           amplitudefullHB += amplitudefullTSs;
         }
         if (HcalDetId(digi.id()).subdet() == HcalEndcap) {
           for (int i = 0; i < digi.samples(); i++)
             amplitudefullTSs += digi[i].adc();
           h_AmplitudeHEtest_->Fill(amplitudefullTSs);
           amplitudefullHE += amplitudefullTSs;
         }
       }
     }
  
     edm::Handle<HFDigiCollection> hfdigi;
     iEvent.getByToken(tok_hfdigi_, hfdigi);
     if (hfdigi.isValid()) {
       for (auto const& digi : *(hfdigi.product())) {
         int nTS = digi.size();
         double amplitudefullTSs = 0.;
         if (digi.id().subdet() == HcalForward) {
           if (nTS <= 10) {
             for (int i = 0; i < nTS; i++)
               amplitudefullTSs += digi.sample(i).adc();
             h_AmplitudeHFtest_->Fill(amplitudefullTSs);
             amplitudefullHF += amplitudefullTSs;
           }
         }
       }
     }
  
     edm::Handle<QIE10DigiCollection> qie10digi;
     iEvent.getByToken(tok_qie10digi_, qie10digi);
     if (qie10digi.isValid()) {
       for (QIE10DataFrame const digi : *(qie10digi.product())) {
         double amplitudefullTSs = 0.;
         if (HcalDetId(digi.id()).subdet() == HcalForward) {
           for (int i = 0; i < digi.samples(); i++)
             amplitudefullTSs += digi[i].adc();
           h_AmplitudeHFtest_->Fill(amplitudefullTSs);
           amplitudefullHF += amplitudefullTSs;
         }
       }
     }
  
     h_AmplitudeHB_->Fill(amplitudefullHB);
     h_AmplitudeHE_->Fill(amplitudefullHE);
     h_AmplitudeHF_->Fill(amplitudefullHF);
   }
  
   edm::Handle<HBHERecHitCollection> hbheMB;
   iEvent.getByToken(tok_hbherecoMB_, hbheMB);
   if (!hbheMB.isValid()) {
     edm::LogWarning("RecAnalyzer") << "HcalCalibAlgos: Error! can't get hbhe product!";
     return;
   }
   const HBHERecHitCollection HithbheMB = *(hbheMB.product());
   HBHEsize = HithbheMB.size();
   edm::LogVerbatim("RecAnalyzer") << "HBHE MB size of collection " << HithbheMB.size();
   if (HithbheMB.size() < 5100 && runNZS_) {
     edm::LogWarning("RecAnalyzer") << "HBHE problem " << rnnum_ << " size " << HBHEsize;
   }
  
   edm::Handle<HFRecHitCollection> hfMB;
   iEvent.getByToken(tok_hfrecoMB_, hfMB);
   if (!hfMB.isValid()) {
     edm::LogWarning("RecAnalyzer") << "HcalCalibAlgos: Error! can't get hf product!";
     return;
   }
   const HFRecHitCollection HithfMB = *(hfMB.product());
   edm::LogVerbatim("RecAnalyzer") << "HF MB size of collection " << HithfMB.size();
   HFsize = HithfMB.size();
   if (HithfMB.size() < 1700 && runNZS_) {
     edm::LogWarning("RecAnalyzer") << "HF problem " << rnnum_ << " size " << HFsize;
   }
  
   bool select(false);
   if (!trigbit_.empty()) {
     edm::Handle<L1GlobalTriggerObjectMapRecord> gtObjectMapRecord;
     iEvent.getByToken(tok_hltL1GtMap_, gtObjectMapRecord);
     if (gtObjectMapRecord.isValid()) {
       const std::vector<L1GlobalTriggerObjectMap>& objMapVec = gtObjectMapRecord->gtObjectMap();
       for (std::vector<L1GlobalTriggerObjectMap>::const_iterator itMap = objMapVec.begin(); itMap != objMapVec.end();
            ++itMap) {
         bool resultGt = (*itMap).algoGtlResult();
         if (resultGt) {
           int algoBit = (*itMap).algoBitNumber();
           if (std::find(trigbit_.begin(), trigbit_.end(), algoBit) != trigbit_.end()) {
             select = true;
             break;
           }
         }
       }
     }
   }
  
   if (!trigbit_.empty() || select)
     myTree1_->Fill();
  
   //event weight for FLAT sample and PU information
   double eventWeight = 1.0;
   edm::Handle<GenEventInfoProduct> genEventInfo;
   iEvent.getByToken(tok_ew_, genEventInfo);
   if (genEventInfo.isValid())
     eventWeight = genEventInfo->weight();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("RecAnalyzer") << "Test HB " << HBHEsize << " HF " << HFsize << " Trigger " << trigbit_.size() << ":"
                                   << select << ":" << ignoreL1_ << " Wt " << eventWeight;
 #endif
   if (ignoreL1_ || (!trigbit_.empty() && select)) {
     analyzeHcal(HithbheMB, HithfMB, 1, true, eventWeight);
   } else if ((!ignoreL1_) && (trigbit_.empty())) {
     edm::Handle<L1GlobalTriggerObjectMapRecord> gtObjectMapRecord;
     iEvent.getByToken(tok_hltL1GtMap_, gtObjectMapRecord);
     if (gtObjectMapRecord.isValid()) {
       const std::vector<L1GlobalTriggerObjectMap>& objMapVec = gtObjectMapRecord->gtObjectMap();
       bool ok(false);
       for (std::vector<L1GlobalTriggerObjectMap>::const_iterator itMap = objMapVec.begin(); itMap != objMapVec.end();
            ++itMap) {
         bool resultGt = (*itMap).algoGtlResult();
         if (resultGt) {
           int algoBit = (*itMap).algoBitNumber();
           analyzeHcal(HithbheMB, HithfMB, algoBit, (!ok), eventWeight);
           ok = true;
         }
       }
       if (!ok) {
         edm::LogInfo("RecAnalyzer") << "No passed L1 Trigger found";
       }
     }
   }
 }

References ecalMGPA::adc(), analyzeHcal(), extraHist_, spr::find(), genParticles2HepMC_cfi::genEventInfo, L1GlobalTriggerObjectMapRecord::gtObjectMap(), h_AmplitudeHB_, h_AmplitudeHBtest_, h_AmplitudeHE_, h_AmplitudeHEtest_, h_AmplitudeHF_, h_AmplitudeHFtest_, HBHEsize, HcalBarrel, HcalEndcap, HcalForward, HFsize, mps_fire::i, iEvent, ignoreL1_, edm::HandleBase::isValid(), myTree1_, convertSQLiteXML::ok, edm::Handle< T >::product(), rnnum_, runNZS_, singleTopDQM_cfi::select, edm::SortedCollection< T, SORT >::size(), tok_ew_, tok_hbhedigi_, tok_hbherecoMB_, tok_hfdigi_, tok_hfrecoMB_, tok_hltL1GtMap_, tok_qie10digi_, tok_qie11digi_, and trigbit_.

◆ analyzeHcal()

void RecAnalyzerMinbias::analyzeHcal	(	const HBHERecHitCollection &	HithbheMB,
		const HFRecHitCollection &	HithfMB,
		int	algoBit,
		bool	fill,
		double	weight
	)

private

Definition at line 592 of file RecAnalyzerMinbias.cc.

                                                                                                                      {
   // Signal part for HB HE
   int count(0), countHB(0), countHE(0), count2(0), count2HB(0), count2HE(0);
   for (HBHERecHitCollection::const_iterator hbheItr = HithbheMB.begin(); hbheItr != HithbheMB.end(); hbheItr++) {
     // Recalibration of energy
     DetId mydetid = hbheItr->id().rawId();
     double icalconst(1.);
     if (theRecalib_) {
       std::map<DetId, double>::iterator itr = corrFactor_.find(mydetid);
       if (itr != corrFactor_.end())
         icalconst = itr->second;
     }
     HBHERecHit aHit(hbheItr->id(), hbheItr->energy() * icalconst, hbheItr->time());
     double energyhit = aHit.energy();
     DetId id = (*hbheItr).detid();
     HcalDetId hid = HcalDetId(id);
     double eLow = (hid.subdet() == HcalEndcap) ? eLowHE_ : eLowHB_;
     double eHigh = (hid.subdet() == HcalEndcap) ? eHighHE_ : eHighHB_;
     ++count;
     if (id.subdetId() == HcalBarrel)
       ++countHB;
     else
       ++countHE;
     if (fill) {
       for (unsigned int i = 0; i < hcalID_.size(); i++) {
         if (hcalID_[i] == id.rawId()) {
           histo_[i]->Fill(energyhit);
           break;
         }
       }
       if (fillHist_) {
         std::map<HcalDetId, TH1D*>::iterator itr1 = histHC_.find(hid);
         if (itr1 != histHC_.end())
           itr1->second->Fill(energyhit);
       }
       h_[hid.subdet() - 1]->Fill(energyhit);
       if (energyhit > eMin_) {
         hbhe_->Fill(hid.ieta(), hid.iphi());
         ++count2;
         if (id.subdetId() == HcalBarrel) {
           ++count2HB;
           hb_->Fill(hid.ieta(), hid.iphi());
         } else {
           ++count2HE;
           he_->Fill(hid.ieta(), hid.iphi());
         }
       }
     }
     if (!fillHist_) {
       if (Noise_ || runNZS_ || (energyhit >= eLow && energyhit <= eHigh)) {
         std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr1 =
             myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
         if (itr1 == myMap_.end()) {
           myInfo info;
           myMap_[std::pair<int, HcalDetId>(algoBit, hid)] = info;
           itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
         }
         itr1->second.theMB0 += weight;
         itr1->second.theMB1 += (weight * energyhit);
         itr1->second.theMB2 += (weight * energyhit * energyhit);
         itr1->second.theMB3 += (weight * energyhit * energyhit * energyhit);
         itr1->second.theMB4 += (weight * energyhit * energyhit * energyhit * energyhit);
         itr1->second.runcheck = rnnum_;
       }
     }
   }  // HBHE_MB
   if (fill) {
     if (count > 0)
       hbherun_->Fill(rnnum_, (double)(count2) / count);
     if (countHB > 0)
       hbrun_->Fill(rnnum_, (double)(count2HB) / countHB);
     if (countHE > 0)
       herun_->Fill(rnnum_, (double)(count2HE) / countHE);
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("RecAnalyzer") << "HBHE " << count2 << ":" << count << ":" << (double)(count2) / count << "\t HB "
                                   << count2HB << ":" << countHB << ":" << (double)(count2HB) / countHB << "\t HE "
                                   << count2HE << ":" << countHE << ":" << (double)(count2HE) / countHE;
 #endif
   int countHF(0), count2HF(0);
   // Signal part for HF
   for (HFRecHitCollection::const_iterator hfItr = HithfMB.begin(); hfItr != HithfMB.end(); hfItr++) {
     // Recalibration of energy
     DetId mydetid = hfItr->id().rawId();
     double icalconst(1.);
     if (theRecalib_) {
       std::map<DetId, double>::iterator itr = corrFactor_.find(mydetid);
       if (itr != corrFactor_.end())
         icalconst = itr->second;
     }
     HFRecHit aHit(hfItr->id(), hfItr->energy() * icalconst, hfItr->time());
  
     double energyhit = aHit.energy();
     DetId id = (*hfItr).detid();
     HcalDetId hid = HcalDetId(id);
     ++countHF;
     if (fill) {
       for (unsigned int i = 0; i < hcalID_.size(); i++) {
         if (hcalID_[i] == id.rawId()) {
           histo_[i]->Fill(energyhit);
           break;
         }
       }
       if (fillHist_) {
         std::map<HcalDetId, TH1D*>::iterator itr1 = histHC_.find(hid);
         if (itr1 != histHC_.end())
           itr1->second->Fill(energyhit);
       }
       h_[hid.subdet() - 1]->Fill(energyhit);
       if (energyhit > eMin_) {
         hf_->Fill(hid.ieta(), hid.iphi());
         ++count2HF;
       }
     }
  
     //
     // Remove PMT hits
     //
     if (!fillHist_) {
       if (((Noise_ || runNZS_) && fabs(energyhit) <= 40.) || (energyhit >= eLowHF_ && energyhit <= eHighHF_)) {
         std::map<std::pair<int, HcalDetId>, myInfo>::iterator itr1 =
             myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
         if (itr1 == myMap_.end()) {
           myInfo info;
           myMap_[std::pair<int, HcalDetId>(algoBit, hid)] = info;
           itr1 = myMap_.find(std::pair<int, HcalDetId>(algoBit, hid));
         }
         itr1->second.theMB0 += weight;
         itr1->second.theMB1 += (weight * energyhit);
         itr1->second.theMB2 += (weight * energyhit * energyhit);
         itr1->second.theMB3 += (weight * energyhit * energyhit * energyhit);
         itr1->second.theMB4 += (weight * energyhit * energyhit * energyhit * energyhit);
         itr1->second.runcheck = rnnum_;
       }
     }
   }
   if (fill && countHF > 0)
     hfrun_->Fill(rnnum_, (double)(count2HF) / countHF);
 #ifdef EDM_ML_DEBUG
   if (count)
     edm::LogVerbatim("RecAnalyzer") << "HF " << count2HF << ":" << countHF << ":" << (double)(count2HF) / countHF;
 #endif
 }

References edm::SortedCollection< T, SORT >::begin(), corrFactor_, KineDebug3::count(), eHighHB_, eHighHE_, eHighHF_, eLowHB_, eLowHE_, eLowHF_, eMin_, edm::SortedCollection< T, SORT >::end(), CaloRecHit::energy(), ntuplemaker::fill, fillHist_, h_, hb_, hbhe_, hbherun_, hbrun_, HcalBarrel, HcalEndcap, hcalID_, he_, herun_, hf_, hfrun_, histHC_, histo_, mps_fire::i, HcalDetId::ieta(), info(), HcalDetId::iphi(), myMap_, Noise_, DetId::rawId(), rnnum_, runNZS_, HcalDetId::subdet(), theRecalib_, and mps_merge::weight.

Referenced by analyze().

◆ beginJob()

void RecAnalyzerMinbias::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 244 of file RecAnalyzerMinbias.cc.

                                   {
   std::string hc[5] = {"Empty", "HB", "HE", "HO", "HF"};
   char name[700], title[700];
   hbhe_ = fs_->make<TH2D>("hbhe", "Noise in HB/HE", 61, -30.5, 30.5, 72, 0.5, 72.5);
   hb_ = fs_->make<TH2D>("hb", "Noise in HB", 61, -16.5, 16.5, 72, 0.5, 72.5);
   he_ = fs_->make<TH2D>("he", "Noise in HE", 61, -30.5, 30.5, 72, 0.5, 72.5);
   hf_ = fs_->make<TH2D>("hf", "Noise in HF", 82, -41.5, 41.5, 72, 0.5, 72.5);
   int nbin = (runMax_ - runMin_ + 1);
   sprintf(title, "Fraction of channels in HB/HE with E > %4.1f GeV vs Run number", eMin_);
   hbherun_ = fs_->make<TProfile>("hbherun", title, nbin, runMin_ - 0.5, runMax_ + 0.5, 0.0, 1.0);
   sprintf(title, "Fraction of channels in HB with E > %4.1f GeV vs Run number", eMin_);
   hbrun_ = fs_->make<TProfile>("hbrun", title, nbin, runMin_ - 0.5, runMax_ + 0.5, 0.0, 1.0);
   sprintf(title, "Fraction of channels in HE with E > %4.1f GeV vs Run number", eMin_);
   herun_ = fs_->make<TProfile>("herun", title, nbin, runMin_ - 0.5, runMax_ + 0.5, 0.0, 1.0);
   sprintf(title, "Fraction of channels in HF with E > %4.1f GeV vs Run number", eMin_);
   hfrun_ = fs_->make<TProfile>("hfrun", title, nbin, runMin_ - 0.5, runMax_ + 0.5, 0.0, 1.0);
   for (int idet = 1; idet <= 4; idet++) {
     sprintf(name, "%s", hc[idet].c_str());
     sprintf(title, "Noise distribution for %s", hc[idet].c_str());
     h_[idet - 1] = fs_->make<TH1D>(name, title, 48, -6., 6.);
   }
  
   for (const auto& hcalid : hcalID_) {
     HcalDetId id = HcalDetId(hcalid);
     int subdet = id.subdetId();
     sprintf(name, "%s%d_%d_%d", hc[subdet].c_str(), id.ieta(), id.iphi(), id.depth());
     sprintf(title,
             "Energy Distribution for %s ieta %d iphi %d depth %d",
             hc[subdet].c_str(),
             id.ieta(),
             id.iphi(),
             id.depth());
     double xmin = (subdet == 4) ? -10 : -1;
     double xmax = (subdet == 4) ? 90 : 9;
     TH1D* hh = fs_->make<TH1D>(name, title, 50, xmin, xmax);
     histo_.push_back(hh);
   };
  
   if (extraHist_) {
     h_AmplitudeHBtest_ = fs_->make<TH1D>("h_AmplitudeHBtest", "", 5000, 0., 5000.);
     h_AmplitudeHEtest_ = fs_->make<TH1D>("h_AmplitudeHEtest", "", 3000, 0., 3000.);
     h_AmplitudeHFtest_ = fs_->make<TH1D>("h_AmplitudeHFtest", "", 10000, 0., 10000.);
     h_AmplitudeHB_ = fs_->make<TH1D>("h_AmplitudeHB", "", 100000, 0., 100000.);
     h_AmplitudeHE_ = fs_->make<TH1D>("h_AmplitudeHE", "", 300000, 0., 300000.);
     h_AmplitudeHF_ = fs_->make<TH1D>("h_AmplitudeHF", "", 100000, 0., 1000000.);
   }
  
   if (!fillHist_) {
     myTree_ = fs_->make<TTree>("RecJet", "RecJet Tree");
     myTree_->Branch("cells", &cells, "cells/I");
     myTree_->Branch("mysubd", &mysubd, "mysubd/I");
     myTree_->Branch("depth", &depth, "depth/I");
     myTree_->Branch("ieta", &ieta, "ieta/I");
     myTree_->Branch("iphi", &iphi, "iphi/I");
     myTree_->Branch("mom0_MB", &mom0_MB, "mom0_MB/F");
     myTree_->Branch("mom1_MB", &mom1_MB, "mom1_MB/F");
     myTree_->Branch("mom2_MB", &mom2_MB, "mom2_MB/F");
     myTree_->Branch("mom3_MB", &mom3_MB, "mom3_MB/F");
     myTree_->Branch("mom4_MB", &mom4_MB, "mom4_MB/F");
     myTree_->Branch("trigbit", &trigbit, "trigbit/I");
     myTree_->Branch("rnnumber", &rnnumber, "rnnumber/D");
   }
   myTree1_ = fs_->make<TTree>("RecJet1", "RecJet1 Tree");
   myTree1_->Branch("rnnum_", &rnnum_, "rnnum_/D");
   myTree1_->Branch("HBHEsize", &HBHEsize, "HBHEsize/I");
   myTree1_->Branch("HFsize", &HFsize, "HFsize/I");
  
   myMap_.clear();
 }

References cells, depth, eMin_, extraHist_, fillHist_, fs_, h_, h_AmplitudeHB_, h_AmplitudeHBtest_, h_AmplitudeHE_, h_AmplitudeHEtest_, h_AmplitudeHF_, h_AmplitudeHFtest_, hb_, hbhe_, hbherun_, HBHEsize, hbrun_, hcalID_, he_, herun_, hf_, hfrun_, HFsize, histo_, ieta, iphi, TFileService::make(), mom0_MB, mom1_MB, mom2_MB, mom3_MB, mom4_MB, myMap_, mysubd, myTree1_, myTree_, Skims_PA_cff::name, ZMuMuCategoriesSequences_cff::nbin, rnnum_, rnnumber, runMax_, runMin_, AlCaHLTBitMon_QueryRunRegistry::string, overlapproblemtsosanalyzer_cfi::title, trigbit, TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

◆ beginRun()

void RecAnalyzerMinbias::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	iS
	)

overrideprivate

Definition at line 349 of file RecAnalyzerMinbias.cc.

                                                                         {
   if (!init_) {
     init_ = true;
     if (fillHist_) {
       edm::ESHandle<HcalTopology> htopo;
       iS.get<IdealGeometryRecord>().get(htopo);
       if (htopo.isValid()) {
         const HcalTopology* hcaltopology = htopo.product();
  
         char name[700], title[700];
         // For HB
         int maxDepthHB = hcaltopology->maxDepthHB();
         int nbinHB = (Noise_) ? 18 : int(2000 * eHighHB_);
         double x_min = (Noise_) ? -3. : 0.;
         double x_max = (Noise_) ? 3. : 2. * eHighHB_;
         for (int eta = -50; eta < 50; eta++) {
           for (int phi = 0; phi < 100; phi++) {
             for (int depth = 1; depth <= maxDepthHB; depth++) {
               HcalDetId cell(HcalBarrel, eta, phi, depth);
               if (hcaltopology->valid(cell)) {
                 sprintf(name, "HBeta%dphi%ddep%d", eta, phi, depth);
                 sprintf(title, "HB #eta %d #phi %d depth %d", eta, phi, depth);
                 TH1D* h = fs_->make<TH1D>(name, title, nbinHB, x_min, x_max);
                 histHC_[cell] = h;
               }
             }
           }
         }
         // For HE
         int maxDepthHE = hcaltopology->maxDepthHE();
         int nbinHE = (Noise_) ? 18 : int(2000 * eHighHE_);
         x_min = (Noise_) ? -3. : 0.;
         x_max = (Noise_) ? 3. : 2. * eHighHE_;
         for (int eta = -50; eta < 50; eta++) {
           for (int phi = 0; phi < 100; phi++) {
             for (int depth = 1; depth <= maxDepthHE; depth++) {
               HcalDetId cell(HcalEndcap, eta, phi, depth);
               if (hcaltopology->valid(cell)) {
                 sprintf(name, "HEeta%dphi%ddep%d", eta, phi, depth);
                 sprintf(title, "HE #eta %d #phi %d depth %d", eta, phi, depth);
                 TH1D* h = fs_->make<TH1D>(name, title, nbinHE, x_min, x_max);
                 histHC_[cell] = h;
               }
             }
           }
         }
         // For HF
         int maxDepthHF = 4;
         int nbinHF = (Noise_) ? 200 : int(2000 * eHighHF_);
         x_min = (Noise_) ? -10. : 0.;
         x_max = (Noise_) ? 10. : 2. * eHighHF_;
         for (int eta = -50; eta < 50; eta++) {
           for (int phi = 0; phi < 100; phi++) {
             for (int depth = 1; depth <= maxDepthHF; depth++) {
               HcalDetId cell(HcalForward, eta, phi, depth);
               if (hcaltopology->valid(cell)) {
                 sprintf(name, "HFeta%dphi%ddep%d", eta, phi, depth);
                 sprintf(title, "Energy (HF #eta %d #phi %d depth %d)", eta, phi, depth);
                 TH1D* h = fs_->make<TH1D>(name, title, nbinHF, x_min, x_max);
                 histHC_[cell] = h;
               }
             }
           }
         }
       }
     }
   }
 }

References depth, eHighHB_, eHighHE_, eHighHF_, PVValHelper::eta, fillHist_, fs_, edm::EventSetup::get(), get, h, HcalBarrel, HcalEndcap, HcalForward, histHC_, init_, edm::ESHandleBase::isValid(), TFileService::make(), hcalSLHCTopologyConstants_cfi::maxDepthHB, HcalTopology::maxDepthHB(), hcalSLHCTopologyConstants_cfi::maxDepthHE, HcalTopology::maxDepthHE(), Skims_PA_cff::name, Noise_, phi, edm::ESHandle< T >::product(), overlapproblemtsosanalyzer_cfi::title, and HcalTopology::valid().

◆ endJob()

void RecAnalyzerMinbias::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 316 of file RecAnalyzerMinbias.cc.

                                 {
   if (!fillHist_) {
     cells = 0;
     for (const auto& itr : myMap_) {
       edm::LogVerbatim("RecAnalyzer") << "Fired trigger bit number " << itr.first.first;
       myInfo info = itr.second;
       if (info.theMB0 > 0) {
         mom0_MB = info.theMB0;
         mom1_MB = info.theMB1;
         mom2_MB = info.theMB2;
         mom3_MB = info.theMB3;
         mom4_MB = info.theMB4;
         rnnumber = info.runcheck;
         trigbit = itr.first.first;
         mysubd = itr.first.second.subdet();
         depth = itr.first.second.depth();
         iphi = itr.first.second.iphi();
         ieta = itr.first.second.ieta();
         edm::LogVerbatim("RecAnalyzer") << " Result=  " << trigbit << " " << mysubd << " " << ieta << " " << iphi
                                         << " mom0  " << mom0_MB << " mom1 " << mom1_MB << " mom2 " << mom2_MB
                                         << " mom3 " << mom3_MB << " mom4 " << mom4_MB;
         myTree_->Fill();
         cells++;
       }
     }
     edm::LogVerbatim("RecAnalyzer") << "cells"
                                     << " " << cells;
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("RecAnalyzer") << "Exiting from RecAnalyzerMinbias::endjob";
 #endif
 }

References cells, depth, fillHist_, ieta, info(), iphi, mom0_MB, mom1_MB, mom2_MB, mom3_MB, mom4_MB, myMap_, mysubd, myTree_, rnnumber, and trigbit.

Referenced by o2olib.O2ORunMgr::executeJob().

◆ endRun()

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

inlineoverrideprivate

Definition at line 80 of file RecAnalyzerMinbias.cc.

80 {}

◆ fillDescriptions()

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

static

Definition at line 213 of file RecAnalyzerMinbias.cc.

                                                                                     {
   std::vector<int> iarray;
   edm::ParameterSetDescription desc;
   desc.add<bool>("runNZS", true);
   desc.add<bool>("noise", false);
   desc.add<double>("eLowHB", 4);
   desc.add<double>("eHighHB", 100);
   desc.add<double>("eLowHE", 4);
   desc.add<double>("eHighHE", 150);
   desc.add<double>("eLowHF", 10);
   desc.add<double>("eHighHF", 150);
   // Suitable cutoff to remove fluctuation of pedestal
   desc.addUntracked<double>("eMin", 2.0);
   // The following run range is suited to study 2017 commissioning period
   desc.addUntracked<int>("runMin", 308327);
   desc.addUntracked<int>("runMax", 308347);
   desc.addUntracked<std::vector<int>>("triggerBits", iarray);
   desc.addUntracked<bool>("ignoreL1", false);
   desc.addUntracked<std::string>("corrFile", "CorFactor.txt");
   desc.addUntracked<bool>("fillHisto", false);
   desc.addUntracked<bool>("extraHisto", false);
   desc.addUntracked<std::vector<int>>("hcalIeta", iarray);
   desc.addUntracked<std::vector<int>>("hcalIphi", iarray);
   desc.addUntracked<std::vector<int>>("hcalDepth", iarray);
   desc.add<edm::InputTag>("hbheInputMB", edm::InputTag("hbherecoMB"));
   desc.add<edm::InputTag>("hfInputMB", edm::InputTag("hfrecoMB"));
   desc.add<edm::InputTag>("gtDigisAlCaMB", edm::InputTag("gtDigisAlCaMB"));
   desc.add<edm::InputTag>("hcalDigiCollectionTag", edm::InputTag("hcalDigis"));
   descriptions.add("recAnalyzerMinbias", desc);
 }

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), HLT_2018_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

Member Data Documentation

◆ cells

int RecAnalyzerMinbias::cells

private

Definition at line 111 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ corrFactor_

std::map<DetId, double> RecAnalyzerMinbias::corrFactor_

private

Definition at line 92 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

◆ depth

int RecAnalyzerMinbias::depth

private

Definition at line 111 of file RecAnalyzerMinbias.cc.

Referenced by dqm_interfaces.DirID::__eq__(), dqm_interfaces.DirID::__repr__(), beginJob(), beginRun(), endJob(), RecAnalyzerMinbias(), and dqm_interfaces.DirWalkerDB::run().

◆ eHighHB_

double RecAnalyzerMinbias::eHighHB_

private

Definition at line 89 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginRun(), and RecAnalyzerMinbias().

◆ eHighHE_

double RecAnalyzerMinbias::eHighHE_

private

Definition at line 89 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginRun(), and RecAnalyzerMinbias().

◆ eHighHF_

double RecAnalyzerMinbias::eHighHF_

private

Definition at line 90 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginRun(), and RecAnalyzerMinbias().

◆ eLowHB_

double RecAnalyzerMinbias::eLowHB_

private

Definition at line 89 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

◆ eLowHE_

double RecAnalyzerMinbias::eLowHE_

private

Definition at line 89 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

◆ eLowHF_

double RecAnalyzerMinbias::eLowHF_

private

Definition at line 90 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

◆ eMin_

double RecAnalyzerMinbias::eMin_

private

Definition at line 90 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginJob(), and RecAnalyzerMinbias().

◆ extraHist_

bool RecAnalyzerMinbias::extraHist_

private

Definition at line 88 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), beginJob(), and RecAnalyzerMinbias().

◆ fillHist_

bool RecAnalyzerMinbias::fillHist_

private

Definition at line 88 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginJob(), beginRun(), endJob(), and RecAnalyzerMinbias().

◆ fs_

edm::Service<TFileService> RecAnalyzerMinbias::fs_

private

Definition at line 86 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and beginRun().

◆ h_

TH1D* RecAnalyzerMinbias::h_[4]

private

Definition at line 95 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

◆ h_AmplitudeHB_

TH1D* RecAnalyzerMinbias::h_AmplitudeHB_

private

Definition at line 99 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

◆ h_AmplitudeHBtest_

TH1D* RecAnalyzerMinbias::h_AmplitudeHBtest_

private

Definition at line 97 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

◆ h_AmplitudeHE_

TH1D * RecAnalyzerMinbias::h_AmplitudeHE_

private

Definition at line 99 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

◆ h_AmplitudeHEtest_

TH1D * RecAnalyzerMinbias::h_AmplitudeHEtest_

private

Definition at line 97 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

◆ h_AmplitudeHF_

TH1D * RecAnalyzerMinbias::h_AmplitudeHF_

private

Definition at line 99 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

◆ h_AmplitudeHFtest_

TH1D* RecAnalyzerMinbias::h_AmplitudeHFtest_

private

Definition at line 98 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

◆ hb_

TH2D * RecAnalyzerMinbias::hb_

private

Definition at line 96 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

◆ hbhe_

TH2D* RecAnalyzerMinbias::hbhe_

private

Definition at line 96 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

◆ hbherun_

TProfile* RecAnalyzerMinbias::hbherun_

private

Definition at line 100 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

◆ HBHEsize

int RecAnalyzerMinbias::HBHEsize

private

Definition at line 113 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

◆ hbrun_

TProfile * RecAnalyzerMinbias::hbrun_

private

Definition at line 100 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

◆ hcalID_

std::vector<unsigned int> RecAnalyzerMinbias::hcalID_

private

Definition at line 93 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginJob(), and RecAnalyzerMinbias().

◆ he_

TH2D * RecAnalyzerMinbias::he_

private

Definition at line 96 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

◆ herun_

TProfile * RecAnalyzerMinbias::herun_

private

Definition at line 100 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

◆ hf_

TH2D * RecAnalyzerMinbias::hf_

private

Definition at line 96 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

◆ hfrun_

TProfile * RecAnalyzerMinbias::hfrun_

private

Definition at line 100 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

◆ HFsize

int RecAnalyzerMinbias::HFsize

private

Definition at line 113 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

◆ histHC_

std::map<HcalDetId, TH1D*> RecAnalyzerMinbias::histHC_

private

Definition at line 102 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginRun().

◆ histo_

std::vector<TH1D*> RecAnalyzerMinbias::histo_

private

Definition at line 101 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

◆ ieta

int RecAnalyzerMinbias::ieta

private

Definition at line 111 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), endJob(), and RecAnalyzerMinbias().

◆ ignoreL1_

bool RecAnalyzerMinbias::ignoreL1_

private

Definition at line 87 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

◆ init_

bool RecAnalyzerMinbias::init_

private

Definition at line 88 of file RecAnalyzerMinbias.cc.

Referenced by beginRun().

◆ iphi

int RecAnalyzerMinbias::iphi

private

Definition at line 111 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), endJob(), and RecAnalyzerMinbias().

◆ mom0_MB

float RecAnalyzerMinbias::mom0_MB

private

Definition at line 112 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom1_MB

float RecAnalyzerMinbias::mom1_MB

private

Definition at line 112 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom2_MB

float RecAnalyzerMinbias::mom2_MB

private

Definition at line 112 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom3_MB

float RecAnalyzerMinbias::mom3_MB

private

Definition at line 112 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ mom4_MB

float RecAnalyzerMinbias::mom4_MB

private

Definition at line 112 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ myMap_

std::map<std::pair<int, HcalDetId>, myInfo> RecAnalyzerMinbias::myMap_

private

Definition at line 114 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginJob(), and endJob().

◆ mysubd

int RecAnalyzerMinbias::mysubd

private

Definition at line 111 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ myTree1_

TTree * RecAnalyzerMinbias::myTree1_

private

Definition at line 94 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

◆ myTree_

TTree* RecAnalyzerMinbias::myTree_

private

Definition at line 94 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ Noise_

bool RecAnalyzerMinbias::Noise_

private

Definition at line 87 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), beginRun(), and RecAnalyzerMinbias().

◆ rnnum_

double RecAnalyzerMinbias::rnnum_

private

Definition at line 104 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), analyzeHcal(), and beginJob().

◆ rnnumber

double RecAnalyzerMinbias::rnnumber

private

Definition at line 110 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ runMax_

int RecAnalyzerMinbias::runMax_

private

Definition at line 91 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and RecAnalyzerMinbias().

◆ runMin_

int RecAnalyzerMinbias::runMin_

private

Definition at line 91 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and RecAnalyzerMinbias().

◆ runNZS_

bool RecAnalyzerMinbias::runNZS_

private

Definition at line 87 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), analyzeHcal(), and RecAnalyzerMinbias().

◆ theRecalib_

bool RecAnalyzerMinbias::theRecalib_

private

Definition at line 87 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

◆ tok_ew_

edm::EDGetTokenT<GenEventInfoProduct> RecAnalyzerMinbias::tok_ew_

private

Definition at line 118 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

◆ tok_gtRec_

edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> RecAnalyzerMinbias::tok_gtRec_

private

Definition at line 124 of file RecAnalyzerMinbias.cc.

Referenced by RecAnalyzerMinbias().

◆ tok_hbhedigi_

edm::EDGetTokenT<HBHEDigiCollection> RecAnalyzerMinbias::tok_hbhedigi_

private

Definition at line 119 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

◆ tok_hbherecoMB_

edm::EDGetTokenT<HBHERecHitCollection> RecAnalyzerMinbias::tok_hbherecoMB_

private

Definition at line 115 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

◆ tok_hfdigi_

edm::EDGetTokenT<HFDigiCollection> RecAnalyzerMinbias::tok_hfdigi_

private

Definition at line 122 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

◆ tok_hfrecoMB_

edm::EDGetTokenT<HFRecHitCollection> RecAnalyzerMinbias::tok_hfrecoMB_

private

Definition at line 116 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

◆ tok_hltL1GtMap_

edm::EDGetTokenT<L1GlobalTriggerObjectMapRecord> RecAnalyzerMinbias::tok_hltL1GtMap_

private

Definition at line 117 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

◆ tok_hodigi_

edm::EDGetTokenT<HODigiCollection> RecAnalyzerMinbias::tok_hodigi_

private

Definition at line 121 of file RecAnalyzerMinbias.cc.

Referenced by RecAnalyzerMinbias().

◆ tok_qie10digi_

edm::EDGetTokenT<QIE10DigiCollection> RecAnalyzerMinbias::tok_qie10digi_

private

Definition at line 123 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

◆ tok_qie11digi_

edm::EDGetTokenT<QIE11DigiCollection> RecAnalyzerMinbias::tok_qie11digi_

private

Definition at line 120 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

◆ trigbit

int RecAnalyzerMinbias::trigbit

private

Definition at line 111 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

◆ trigbit_

std::vector<int> RecAnalyzerMinbias::trigbit_

private

Definition at line 103 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

Classes

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ RecAnalyzerMinbias()

◆ ~RecAnalyzerMinbias()

Member Function Documentation

◆ analyze()

◆ analyzeHcal()

◆ beginJob()

◆ beginRun()

◆ endJob()

◆ endRun()

◆ fillDescriptions()

Member Data Documentation

◆ cells

◆ corrFactor_

◆ depth

◆ eHighHB_

◆ eHighHE_

◆ eHighHF_

◆ eLowHB_

◆ eLowHE_

◆ eLowHF_

◆ eMin_

◆ extraHist_

◆ fillHist_

◆ fs_

◆ h_

◆ h_AmplitudeHB_

◆ h_AmplitudeHBtest_

◆ h_AmplitudeHE_

◆ h_AmplitudeHEtest_

◆ h_AmplitudeHF_

◆ h_AmplitudeHFtest_

◆ hb_

◆ hbhe_

◆ hbherun_

◆ HBHEsize

◆ hbrun_

◆ hcalID_

◆ he_

◆ herun_

◆ hf_

◆ hfrun_

◆ HFsize

◆ histHC_

◆ histo_

◆ ieta

◆ ignoreL1_

◆ init_

◆ iphi

◆ mom0_MB

◆ mom1_MB

◆ mom2_MB

◆ mom3_MB

◆ mom4_MB

◆ myMap_

◆ mysubd

◆ myTree1_

◆ myTree_

◆ Noise_

◆ rnnum_

◆ rnnumber

◆ runMax_

◆ runMin_

◆ runNZS_

◆ theRecalib_

◆ tok_ew_

◆ tok_gtRec_

◆ tok_hbhedigi_

◆ tok_hbherecoMB_

◆ tok_hfdigi_

◆ tok_hfrecoMB_

◆ tok_hltL1GtMap_

◆ tok_hodigi_