RecAnalyzerMinbias Class Reference

#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 const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
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 const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
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
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
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 ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
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)

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 ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 131 of file RecAnalyzerMinbias.cc.

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_, iphi, gen::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_.

                                                                     :
  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];
}

RecAnalyzerMinbias::~RecAnalyzerMinbias ( )

override

Definition at line 221 of file RecAnalyzerMinbias.cc.

{}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 435 of file RecAnalyzerMinbias.cc.

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

                                                                                    {

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

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

private

Definition at line 610 of file RecAnalyzerMinbias.cc.

References edm::SortedCollection< T, SORT >::begin(), corrFactor_, KineDebug3::count(), DEFINE_FWK_MODULE, eHighHB_, eHighHE_, eHighHF_, eLowHB_, eLowHE_, eLowHF_, eMin_, edm::SortedCollection< T, SORT >::end(), CaloRecHit::energy(), 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(), and endRun().

                                                            {
  // 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
}

void RecAnalyzerMinbias::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 255 of file RecAnalyzerMinbias.cc.

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_, AnalysisDataFormats_SUSYBSMObjects::hc, 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_, dataset::name, rnnum_, rnnumber, runMax_, runMin_, AlCaHLTBitMon_QueryRunRegistry::string, fftjetcommon_cfi::title, trigbit, TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

                                   {

   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();
 }

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

overrideprivate

Definition at line 360 of file RecAnalyzerMinbias.cc.

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

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

void RecAnalyzerMinbias::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 323 of file RecAnalyzerMinbias.cc.

References cells, depth, fillHist_, ieta, info(), iphi, mom0_MB, mom1_MB, mom2_MB, mom3_MB, mom4_MB, myMap_, mysubd, myTree_, rnnumber, RecAnalyzerMinbias::myInfo::runcheck, RecAnalyzerMinbias::myInfo::theMB0, RecAnalyzerMinbias::myInfo::theMB1, RecAnalyzerMinbias::myInfo::theMB2, RecAnalyzerMinbias::myInfo::theMB3, RecAnalyzerMinbias::myInfo::theMB4, and trigbit.

Referenced by o2olib.O2ORunMgr::executeJob().

                                {

  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
}

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

inlineoverrideprivate

Definition at line 81 of file RecAnalyzerMinbias.cc.

References analyzeHcal().

{}

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

static

Definition at line 223 of file RecAnalyzerMinbias.cc.

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

                                                                                    {

  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);
}

Member Data Documentation

int RecAnalyzerMinbias::cells

private

Definition at line 114 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

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

private

Definition at line 93 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

int RecAnalyzerMinbias::depth

private

Definition at line 114 of file RecAnalyzerMinbias.cc.

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

double RecAnalyzerMinbias::eHighHB_

private

Definition at line 90 of file RecAnalyzerMinbias.cc.

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

double RecAnalyzerMinbias::eHighHE_

private

Definition at line 90 of file RecAnalyzerMinbias.cc.

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

double RecAnalyzerMinbias::eHighHF_

private

Definition at line 91 of file RecAnalyzerMinbias.cc.

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

double RecAnalyzerMinbias::eLowHB_

private

Definition at line 90 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

double RecAnalyzerMinbias::eLowHE_

private

Definition at line 90 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

double RecAnalyzerMinbias::eLowHF_

private

Definition at line 91 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

double RecAnalyzerMinbias::eMin_

private

Definition at line 91 of file RecAnalyzerMinbias.cc.

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

bool RecAnalyzerMinbias::extraHist_

private

Definition at line 89 of file RecAnalyzerMinbias.cc.

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

bool RecAnalyzerMinbias::fillHist_

private

Definition at line 89 of file RecAnalyzerMinbias.cc.

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

edm::Service<TFileService> RecAnalyzerMinbias::fs_

private

Definition at line 87 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and beginRun().

TH1D* RecAnalyzerMinbias::h_[4]

private

Definition at line 96 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

TH1D* RecAnalyzerMinbias::h_AmplitudeHB_

private

Definition at line 100 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

TH1D* RecAnalyzerMinbias::h_AmplitudeHBtest_

private

Definition at line 98 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

TH1D * RecAnalyzerMinbias::h_AmplitudeHE_

private

Definition at line 100 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

TH1D * RecAnalyzerMinbias::h_AmplitudeHEtest_

private

Definition at line 98 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

TH1D * RecAnalyzerMinbias::h_AmplitudeHF_

private

Definition at line 100 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

TH1D* RecAnalyzerMinbias::h_AmplitudeHFtest_

private

Definition at line 99 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

TH2D * RecAnalyzerMinbias::hb_

private

Definition at line 97 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

TH2D* RecAnalyzerMinbias::hbhe_

private

Definition at line 97 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

TProfile* RecAnalyzerMinbias::hbherun_

private

Definition at line 101 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

int RecAnalyzerMinbias::HBHEsize

private

Definition at line 116 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

TProfile * RecAnalyzerMinbias::hbrun_

private

Definition at line 101 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

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

private

Definition at line 94 of file RecAnalyzerMinbias.cc.

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

TH2D * RecAnalyzerMinbias::he_

private

Definition at line 97 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

TProfile * RecAnalyzerMinbias::herun_

private

Definition at line 101 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

TH2D * RecAnalyzerMinbias::hf_

private

Definition at line 97 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

TProfile * RecAnalyzerMinbias::hfrun_

private

Definition at line 101 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

int RecAnalyzerMinbias::HFsize

private

Definition at line 116 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 103 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginRun().

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

private

Definition at line 102 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

int RecAnalyzerMinbias::ieta

private

Definition at line 114 of file RecAnalyzerMinbias.cc.

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

bool RecAnalyzerMinbias::ignoreL1_

private

Definition at line 88 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

bool RecAnalyzerMinbias::init_

private

Definition at line 89 of file RecAnalyzerMinbias.cc.

Referenced by beginRun().

int RecAnalyzerMinbias::iphi

private

Definition at line 114 of file RecAnalyzerMinbias.cc.

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

float RecAnalyzerMinbias::mom0_MB

private

Definition at line 115 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float RecAnalyzerMinbias::mom1_MB

private

Definition at line 115 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float RecAnalyzerMinbias::mom2_MB

private

Definition at line 115 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float RecAnalyzerMinbias::mom3_MB

private

Definition at line 115 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

float RecAnalyzerMinbias::mom4_MB

private

Definition at line 115 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

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

private

Definition at line 117 of file RecAnalyzerMinbias.cc.

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

int RecAnalyzerMinbias::mysubd

private

Definition at line 114 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

TTree * RecAnalyzerMinbias::myTree1_

private

Definition at line 95 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

TTree* RecAnalyzerMinbias::myTree_

private

Definition at line 95 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

bool RecAnalyzerMinbias::Noise_

private

Definition at line 88 of file RecAnalyzerMinbias.cc.

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

double RecAnalyzerMinbias::rnnum_

private

Definition at line 105 of file RecAnalyzerMinbias.cc.

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

double RecAnalyzerMinbias::rnnumber

private

Definition at line 113 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

int RecAnalyzerMinbias::runMax_

private

Definition at line 92 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and RecAnalyzerMinbias().

int RecAnalyzerMinbias::runMin_

private

Definition at line 92 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and RecAnalyzerMinbias().

bool RecAnalyzerMinbias::runNZS_

private

Definition at line 88 of file RecAnalyzerMinbias.cc.

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

bool RecAnalyzerMinbias::theRecalib_

private

Definition at line 88 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

edm::EDGetTokenT<GenEventInfoProduct> RecAnalyzerMinbias::tok_ew_

private

Definition at line 121 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> RecAnalyzerMinbias::tok_gtRec_

private

Definition at line 127 of file RecAnalyzerMinbias.cc.

Referenced by RecAnalyzerMinbias().

edm::EDGetTokenT<HBHEDigiCollection> RecAnalyzerMinbias::tok_hbhedigi_

private

Definition at line 122 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

edm::EDGetTokenT<HBHERecHitCollection> RecAnalyzerMinbias::tok_hbherecoMB_

private

Definition at line 118 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

edm::EDGetTokenT<HFDigiCollection> RecAnalyzerMinbias::tok_hfdigi_

private

Definition at line 125 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

edm::EDGetTokenT<HFRecHitCollection> RecAnalyzerMinbias::tok_hfrecoMB_

private

Definition at line 119 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

edm::EDGetTokenT<L1GlobalTriggerObjectMapRecord> RecAnalyzerMinbias::tok_hltL1GtMap_

private

Definition at line 120 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

edm::EDGetTokenT<HODigiCollection> RecAnalyzerMinbias::tok_hodigi_

private

Definition at line 124 of file RecAnalyzerMinbias.cc.

Referenced by RecAnalyzerMinbias().

edm::EDGetTokenT<QIE10DigiCollection> RecAnalyzerMinbias::tok_qie10digi_

private

Definition at line 126 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

edm::EDGetTokenT<QIE11DigiCollection> RecAnalyzerMinbias::tok_qie11digi_

private

Definition at line 123 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

int RecAnalyzerMinbias::trigbit

private

Definition at line 114 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

std::vector<int> RecAnalyzerMinbias::trigbit_

private

Definition at line 104 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().