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=default
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept 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 noexcept
bool	wantsStreamRuns () const noexcept
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

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
const std::string	cfile_
std::map< DetId, double >	corrFactor_
int	depth
const double	eHighHB_
const double	eHighHE_
const double	eHighHF_
const double	eLowHB_
const double	eLowHE_
const double	eLowHF_
const double	eMin_
const bool	extraHist_
const bool	fillHist_
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
const 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_
const bool	Noise_
double	rnnum_
double	rnnumber
const int	runMax_
const int	runMin_
const bool	runNZS_
bool	theRecalib_
const edm::EDGetTokenT< GenEventInfoProduct >	tok_ew_
const edm::EDGetTokenT< L1GlobalTriggerReadoutRecord >	tok_gtRec_
const edm::EDGetTokenT< HBHEDigiCollection >	tok_hbhedigi_
const edm::EDGetTokenT< HBHERecHitCollection >	tok_hbherecoMB_
const edm::EDGetTokenT< HFDigiCollection >	tok_hfdigi_
const edm::EDGetTokenT< HFRecHitCollection >	tok_hfrecoMB_
const edm::EDGetTokenT< L1GlobalTriggerObjectMapRecord >	tok_hltL1GtMap_
const edm::EDGetTokenT< HODigiCollection >	tok_hodigi_
const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	tok_htopo_
const edm::EDGetTokenT< QIE10DigiCollection >	tok_qie10digi_
const edm::EDGetTokenT< QIE11DigiCollection >	tok_qie11digi_
int	trigbit
const 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)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

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

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

Definition at line 67 of file RecAnalyzerMinbias.cc.

Constructor & Destructor Documentation

RecAnalyzerMinbias()

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

explicit

Definition at line 129 of file RecAnalyzerMinbias.cc.

References funct::abs(), cfile_, corrFactor_, depth, hcalRecHitTable_cff::detId, eHighHB_, eHighHE_, eHighHF_, eLowHB_, eLowHE_, eLowHF_, edm::ParameterSet::getUntrackedParameter(), HcalBarrel, HcalEndcap, HcalForward, hcalID_, HcalOuter, EcalPhiSymFlatTableProducers_cfi::id, ieta, ignoreL1_, timingPdfMaker::infile, iphi, isotrackApplyRegressor::k, nano_mu_digi_cff::rawId, runNZS_, theRecalib_, and trigbit_.

    : runNZS_(iConfig.getParameter<bool>("runNZS")),
      Noise_(iConfig.getParameter<bool>("noise")),
      ignoreL1_(iConfig.getUntrackedParameter<bool>("ignoreL1", false)),
      fillHist_(iConfig.getUntrackedParameter<bool>("fillHisto", false)),
      extraHist_(iConfig.getUntrackedParameter<bool>("extraHisto", false)),
      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)),
      runMin_(iConfig.getUntrackedParameter<int>("RunMin", 308327)),
      runMax_(iConfig.getUntrackedParameter<int>("RunMax", 315250)),
      trigbit_(iConfig.getUntrackedParameter<std::vector<int>>("triggerBits")),
      cfile_(iConfig.getUntrackedParameter<std::string>("corrFile")),
      init_(false),
      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"))),
      tok_htopo_(esConsumes<HcalTopology, HcalRecNumberingRecord, edm::Transition::BeginRun>()) {
  usesResource("TFileService");

  // get constants for DetID's
  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");

  // Read correction factors
  std::ifstream infile(cfile_.c_str());
  if (!infile.is_open()) {
    theRecalib_ = false;
    edm::LogWarning("RecAnalyzerMinbias") << "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::LogVerbatim("RecAnalyzerMinbias") << "Reads " << nrec << " correction factors for " << ndets << " detIds";
    theRecalib_ = (ndets > 0);
  }

  edm::LogVerbatim("RecAnalyzerMinbias") << " Flags (ReCalib): " << theRecalib_ << " (IgnoreL1): " << ignoreL1_
                                         << " (NZS) " << runNZS_ << " and with " << ieta.size()
                                         << " detId for full histogram";
  edm::LogVerbatim("RecAnalyzerMinbias") << "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::LogVerbatim("RecAnalyzerMinbias") << "DetId[" << k << "] " << HcalDetId(id);
  }
  edm::LogVerbatim("RecAnalyzerMinbias") << "Select on " << trigbit_.size() << " L1 Trigger selection";
  for (unsigned int k = 0; k < trigbit_.size(); ++k)
    edm::LogVerbatim("RecAnalyzerMinbias") << "Bit[" << k << "] " << trigbit_[k];
}

~RecAnalyzerMinbias()

RecAnalyzerMinbias::~RecAnalyzerMinbias ( )

overridedefault

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 418 of file RecAnalyzerMinbias.cc.

References gpuClustering::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_.

                                                                                    {
  rnnum_ = (double)iEvent.run();

  if (extraHist_) {
    double amplitudefullHB(0), amplitudefullHE(0), amplitudefullHF(0);
    const edm::Handle<HBHEDigiCollection>& hbhedigi = iEvent.getHandle(tok_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;
          }
        }
      }
    }

    const edm::Handle<QIE11DigiCollection>& qie11digi = iEvent.getHandle(tok_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;
        }
      }
    }

    const edm::Handle<HFDigiCollection>& hfdigi = iEvent.getHandle(tok_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;
          }
        }
      }
    }

    const edm::Handle<QIE10DigiCollection>& qie10digi = iEvent.getHandle(tok_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);
  }

  const edm::Handle<HBHERecHitCollection>& hbheMB = iEvent.getHandle(tok_hbherecoMB_);
  if (!hbheMB.isValid()) {
    edm::LogWarning("RecAnalyzerMinbias") << "HcalCalibAlgos: Error! can't get hbhe product!";
    return;
  }
  const HBHERecHitCollection HithbheMB = *(hbheMB.product());
  HBHEsize = HithbheMB.size();
  edm::LogVerbatim("RecAnalyzerMinbias") << "HBHE MB size of collection " << HithbheMB.size();
  if (HithbheMB.size() < 5100 && runNZS_) {
    edm::LogWarning("RecAnalyzerMinbias") << "HBHE problem " << rnnum_ << " size " << HBHEsize;
  }

  const edm::Handle<HFRecHitCollection> hfMB = iEvent.getHandle(tok_hfrecoMB_);
  if (!hfMB.isValid()) {
    edm::LogWarning("RecAnalyzerMinbias") << "HcalCalibAlgos: Error! can't get hf product!";
    return;
  }
  const HFRecHitCollection HithfMB = *(hfMB.product());
  edm::LogVerbatim("RecAnalyzerMinbias") << "HF MB size of collection " << HithfMB.size();
  HFsize = HithfMB.size();
  if (HithfMB.size() < 1700 && runNZS_) {
    edm::LogWarning("RecAnalyzerMinbias") << "HF problem " << rnnum_ << " size " << HFsize;
  }

  bool select(false);
  if (!trigbit_.empty()) {
    const edm::Handle<L1GlobalTriggerObjectMapRecord>& gtObjectMapRecord = iEvent.getHandle(tok_hltL1GtMap_);
    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;
  const edm::Handle<GenEventInfoProduct>& genEventInfo = iEvent.getHandle(tok_ew_);
  if (genEventInfo.isValid())
    eventWeight = genEventInfo->weight();
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("RecAnalyzerMinbias") << "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())) {
    const edm::Handle<L1GlobalTriggerObjectMapRecord>& gtObjectMapRecord = iEvent.getHandle(tok_hltL1GtMap_);
    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::LogVerbatim("RecAnalyzerMinbias") << "No passed L1 Trigger found";
      }
    }
  }
}

analyzeHcal()

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

private

Definition at line 578 of file RecAnalyzerMinbias.cc.

References edm::SortedCollection< T, SORT >::begin(), corrFactor_, submitPVResolutionJobs::count, eHighHB_, eHighHE_, eHighHF_, eLowHB_, eLowHE_, eLowHF_, eMin_, edm::SortedCollection< T, SORT >::end(), CaloRecHit::energy(), EcalPhiSymFlatTableProducers_cfi::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(), nano_mu_digi_cff::rawId, rnnum_, runNZS_, HcalDetId::subdet(), theRecalib_, and mps_merge::weight.

Referenced by analyze().

                                                                                                                     {
  // 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("RecAnalyzerMinbias") << "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("RecAnalyzerMinbias")
        << "HF " << count2HF << ":" << countHF << ":" << (double)(count2HF) / countHF;
#endif
}

beginJob()

void RecAnalyzerMinbias::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 242 of file RecAnalyzerMinbias.cc.

References cells, depth, eMin_, extraHist_, fillHist_, make_classfiles::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, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::hh, histo_, ieta, iphi, mom0_MB, mom1_MB, mom2_MB, mom3_MB, mom4_MB, myMap_, mysubd, myTree1_, myTree_, mergeVDriftHistosByStation::name, rnnum_, rnnumber, runMax_, runMin_, AlCaHLTBitMon_QueryRunRegistry::string, runGCPTkAlMap::title, trigbit, TrackerOfflineValidation_Dqm_cff::xmax, and TrackerOfflineValidation_Dqm_cff::xmin.

                                  {
  edm::Service<TFileService> fs;
  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();
}

beginRun()

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

overrideprivate

Definition at line 347 of file RecAnalyzerMinbias.cc.

References depth, eHighHB_, eHighHE_, eHighHF_, PVValHelper::eta, fillHist_, make_classfiles::fs, edm::EventSetup::getData(), h, HcalBarrel, HcalEndcap, HcalForward, histHC_, init_, hcalSLHCTopologyConstants_cfi::maxDepthHB, HcalTopology::maxDepthHB(), hcalSLHCTopologyConstants_cfi::maxDepthHE, HcalTopology::maxDepthHE(), mergeVDriftHistosByStation::name, Noise_, phi, runGCPTkAlMap::title, tok_htopo_, HcalTopology::valid(), ppsFastLocalSimulation_cfi::x_max, and ppsFastLocalSimulation_cfi::x_min.

                                                                        {
  if (!init_) {
    init_ = true;
    if (fillHist_) {
      edm::Service<TFileService> fs;
      const HcalTopology* hcaltopology = &iS.getData(tok_htopo_);

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

endJob()

void RecAnalyzerMinbias::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 315 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, and trigbit.

Referenced by o2olib.O2ORunMgr::executeJob().

                                {
  if (!fillHist_) {
    cells = 0;
    for (const auto& itr : myMap_) {
      edm::LogVerbatim("RecAnalyzerMinbias") << "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("RecAnalyzerMinbias")
            << " 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("RecAnalyzerMinbias") << "cells " << cells;
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("RecAnalyzerMinbias") << "Exiting from RecAnalyzerMinbias::endjob";
#endif
}

endRun()

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

inlineoverrideprivate

Definition at line 79 of file RecAnalyzerMinbias.cc.

{}

fillDescriptions()

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

static

Definition at line 211 of file RecAnalyzerMinbias.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, 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

cells

int RecAnalyzerMinbias::cells

private

Definition at line 111 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

cfile_

const std::string RecAnalyzerMinbias::cfile_

private

Definition at line 91 of file RecAnalyzerMinbias.cc.

Referenced by RecAnalyzerMinbias().

corrFactor_

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

private

Definition at line 93 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_

const double RecAnalyzerMinbias::eHighHB_

private

Definition at line 87 of file RecAnalyzerMinbias.cc.

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

eHighHE_

const double RecAnalyzerMinbias::eHighHE_

private

Definition at line 87 of file RecAnalyzerMinbias.cc.

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

eHighHF_

const double RecAnalyzerMinbias::eHighHF_

private

Definition at line 88 of file RecAnalyzerMinbias.cc.

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

eLowHB_

const double RecAnalyzerMinbias::eLowHB_

private

Definition at line 87 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

eLowHE_

const double RecAnalyzerMinbias::eLowHE_

private

Definition at line 87 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

eLowHF_

const double RecAnalyzerMinbias::eLowHF_

private

Definition at line 88 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

eMin_

const double RecAnalyzerMinbias::eMin_

private

Definition at line 88 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

extraHist_

const bool RecAnalyzerMinbias::extraHist_

private

Definition at line 86 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

fillHist_

const bool RecAnalyzerMinbias::fillHist_

private

Definition at line 86 of file RecAnalyzerMinbias.cc.

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

h_

TH1D* RecAnalyzerMinbias::h_[4]

private

Definition at line 96 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

h_AmplitudeHB_

TH1D* RecAnalyzerMinbias::h_AmplitudeHB_

private

Definition at line 100 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

h_AmplitudeHBtest_

TH1D* RecAnalyzerMinbias::h_AmplitudeHBtest_

private

Definition at line 98 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

h_AmplitudeHE_

TH1D * RecAnalyzerMinbias::h_AmplitudeHE_

private

Definition at line 100 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

h_AmplitudeHEtest_

TH1D * RecAnalyzerMinbias::h_AmplitudeHEtest_

private

Definition at line 98 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

h_AmplitudeHF_

TH1D * RecAnalyzerMinbias::h_AmplitudeHF_

private

Definition at line 100 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

h_AmplitudeHFtest_

TH1D* RecAnalyzerMinbias::h_AmplitudeHFtest_

private

Definition at line 99 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

hb_

TH2D * RecAnalyzerMinbias::hb_

private

Definition at line 97 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

hbhe_

TH2D* RecAnalyzerMinbias::hbhe_

private

Definition at line 97 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

hbherun_

TProfile* RecAnalyzerMinbias::hbherun_

private

Definition at line 101 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 101 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

hcalID_

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

private

Definition at line 94 of file RecAnalyzerMinbias.cc.

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

he_

TH2D * RecAnalyzerMinbias::he_

private

Definition at line 97 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

herun_

TProfile * RecAnalyzerMinbias::herun_

private

Definition at line 101 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

hf_

TH2D * RecAnalyzerMinbias::hf_

private

Definition at line 97 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginJob().

hfrun_

TProfile * RecAnalyzerMinbias::hfrun_

private

Definition at line 101 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 103 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginRun().

histo_

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

private

Definition at line 102 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_

const bool RecAnalyzerMinbias::ignoreL1_

private

Definition at line 85 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().

init_

bool RecAnalyzerMinbias::init_

private

Definition at line 92 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 95 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and beginJob().

myTree_

TTree* RecAnalyzerMinbias::myTree_

private

Definition at line 95 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

Noise_

const bool RecAnalyzerMinbias::Noise_

private

Definition at line 85 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and beginRun().

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_

const int RecAnalyzerMinbias::runMax_

private

Definition at line 89 of file RecAnalyzerMinbias.cc.

Referenced by beginJob().

runMin_

const int RecAnalyzerMinbias::runMin_

private

Definition at line 89 of file RecAnalyzerMinbias.cc.

Referenced by beginJob().

runNZS_

const bool RecAnalyzerMinbias::runNZS_

private

Definition at line 85 of file RecAnalyzerMinbias.cc.

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

theRecalib_

bool RecAnalyzerMinbias::theRecalib_

private

Definition at line 92 of file RecAnalyzerMinbias.cc.

Referenced by analyzeHcal(), and RecAnalyzerMinbias().

tok_ew_

const edm::EDGetTokenT<GenEventInfoProduct> RecAnalyzerMinbias::tok_ew_

private

Definition at line 118 of file RecAnalyzerMinbias.cc.

Referenced by analyze().

tok_gtRec_

const edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> RecAnalyzerMinbias::tok_gtRec_

private

Definition at line 124 of file RecAnalyzerMinbias.cc.

tok_hbhedigi_

const edm::EDGetTokenT<HBHEDigiCollection> RecAnalyzerMinbias::tok_hbhedigi_

private

Definition at line 119 of file RecAnalyzerMinbias.cc.

Referenced by analyze().

tok_hbherecoMB_

const edm::EDGetTokenT<HBHERecHitCollection> RecAnalyzerMinbias::tok_hbherecoMB_

private

Definition at line 115 of file RecAnalyzerMinbias.cc.

Referenced by analyze().

tok_hfdigi_

const edm::EDGetTokenT<HFDigiCollection> RecAnalyzerMinbias::tok_hfdigi_

private

Definition at line 122 of file RecAnalyzerMinbias.cc.

Referenced by analyze().

tok_hfrecoMB_

const edm::EDGetTokenT<HFRecHitCollection> RecAnalyzerMinbias::tok_hfrecoMB_

private

Definition at line 116 of file RecAnalyzerMinbias.cc.

Referenced by analyze().

tok_hltL1GtMap_

const edm::EDGetTokenT<L1GlobalTriggerObjectMapRecord> RecAnalyzerMinbias::tok_hltL1GtMap_

private

Definition at line 117 of file RecAnalyzerMinbias.cc.

Referenced by analyze().

tok_hodigi_

const edm::EDGetTokenT<HODigiCollection> RecAnalyzerMinbias::tok_hodigi_

private

Definition at line 121 of file RecAnalyzerMinbias.cc.

tok_htopo_

const edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> RecAnalyzerMinbias::tok_htopo_

private

Definition at line 125 of file RecAnalyzerMinbias.cc.

Referenced by beginRun().

tok_qie10digi_

const edm::EDGetTokenT<QIE10DigiCollection> RecAnalyzerMinbias::tok_qie10digi_

private

Definition at line 123 of file RecAnalyzerMinbias.cc.

Referenced by analyze().

tok_qie11digi_

const edm::EDGetTokenT<QIE11DigiCollection> RecAnalyzerMinbias::tok_qie11digi_

private

Definition at line 120 of file RecAnalyzerMinbias.cc.

Referenced by analyze().

trigbit

int RecAnalyzerMinbias::trigbit

private

Definition at line 111 of file RecAnalyzerMinbias.cc.

Referenced by beginJob(), and endJob().

trigbit_

const std::vector<int> RecAnalyzerMinbias::trigbit_

private

Definition at line 90 of file RecAnalyzerMinbias.cc.

Referenced by analyze(), and RecAnalyzerMinbias().