CaloTowersAnalyzer Class Reference

#include <CaloTowersAnalyzer.h>

Inheritance diagram for CaloTowersAnalyzer:

Public Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	CaloTowersAnalyzer (edm::ParameterSet const &conf)
void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
	~CaloTowersAnalyzer () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Types
typedef math::RhoEtaPhiVector	Vector

Private Member Functions
double	dR (double eta1, double phi1, double eta2, double phi2)

Private Attributes
MonitorElement *	emean_vs_ieta_E
MonitorElement *	emean_vs_ieta_E1
MonitorElement *	emean_vs_ieta_EH
MonitorElement *	emean_vs_ieta_EH1
MonitorElement *	emean_vs_ieta_H
MonitorElement *	emean_vs_ieta_H1
MonitorElement *	emEnergyTiming_HB
MonitorElement *	emEnergyTiming_HE
MonitorElement *	emEnergyTiming_HF
MonitorElement *	emEnergyTiming_High_HB
MonitorElement *	emEnergyTiming_Low_HB
MonitorElement *	emEnergyTiming_Low_HE
MonitorElement *	emEnergyTiming_profile_HB
MonitorElement *	emEnergyTiming_profile_HE
MonitorElement *	emEnergyTiming_profile_HF
MonitorElement *	emEnergyTiming_profile_High_HB
MonitorElement *	emEnergyTiming_profile_Low_HB
MonitorElement *	emEnergyTiming_profile_Low_HE
MonitorElement *	emTiming_HB
MonitorElement *	emTiming_HE
MonitorElement *	emTiming_HF
double	etaMax [3]
double	etaMin [3]
MonitorElement *	hadEnergyTiming_HB
MonitorElement *	hadEnergyTiming_HE
MonitorElement *	hadEnergyTiming_HF
MonitorElement *	hadEnergyTiming_High_HB
MonitorElement *	hadEnergyTiming_Low_HB
MonitorElement *	hadEnergyTiming_Low_HE
MonitorElement *	hadEnergyTiming_Low_HF
MonitorElement *	hadEnergyTiming_profile_HB
MonitorElement *	hadEnergyTiming_profile_HE
MonitorElement *	hadEnergyTiming_profile_HF
MonitorElement *	hadEnergyTiming_profile_High_HB
MonitorElement *	hadEnergyTiming_profile_Low_HB
MonitorElement *	hadEnergyTiming_profile_Low_HE
MonitorElement *	hadEnergyTiming_profile_Low_HF
MonitorElement *	hadTiming_HB
MonitorElement *	hadTiming_HE
MonitorElement *	hadTiming_HF
edm::ESGetToken < HcalDDDRecConstants, HcalRecNumberingRecord >	hcalDDDRecConstantsToken_
std::string	hcalselector_
int	ieta_bins_
float	ieta_max_
float	ieta_min_
int	imc
int	iphi_bins_
float	iphi_max_
float	iphi_min_
int	isub
MonitorElement *	mapEnergy_E
MonitorElement *	mapEnergy_EH
MonitorElement *	mapEnergy_H
MonitorElement *	mapEnergy_HB
MonitorElement *	mapEnergy_HE
MonitorElement *	mapEnergy_HF
MonitorElement *	mapEnergy_N
MonitorElement *	mapEnergyEcal_HB
MonitorElement *	mapEnergyEcal_HE
MonitorElement *	mapEnergyEcal_HF
MonitorElement *	mapEnergyHcal_HB
MonitorElement *	mapEnergyHcal_HE
MonitorElement *	mapEnergyHcal_HF
std::string	mc_
MonitorElement *	meEnergyEcal_HB
MonitorElement *	meEnergyEcal_HE
MonitorElement *	meEnergyEcal_HF
MonitorElement *	meEnergyEcalTower_HB
MonitorElement *	meEnergyEcalTower_HE
MonitorElement *	meEnergyEcalTower_HF
MonitorElement *	meEnergyHcal_HB
MonitorElement *	meEnergyHcal_HE
MonitorElement *	meEnergyHcal_HF
MonitorElement *	meEnergyHcalTower_HB
MonitorElement *	meEnergyHcalTower_HE
MonitorElement *	meEnergyHcalTower_HF
MonitorElement *	meEnergyHcalvsEcal_HB
MonitorElement *	meEnergyHcalvsEcal_HE
MonitorElement *	meEnergyHcalvsEcal_HF
MonitorElement *	meEnergyHO_HB
MonitorElement *	meEnergyHO_HE
MonitorElement *	meEnergyHO_HF
MonitorElement *	meIphiCaloTower_HFM
MonitorElement *	meIphiCaloTower_HFP
MonitorElement *	meIphiHcalTower_HBM
MonitorElement *	meIphiHcalTower_HBP
MonitorElement *	meIphiHcalTower_HEM
MonitorElement *	meIphiHcalTower_HEP
MonitorElement *	meNumFiredTowers_HB
MonitorElement *	meNumFiredTowers_HE
MonitorElement *	meNumFiredTowers_HF
MonitorElement *	MET_HB
MonitorElement *	MET_HE
MonitorElement *	MET_HF
MonitorElement *	meTotEnergy_HB
MonitorElement *	meTotEnergy_HE
MonitorElement *	meTotEnergy_HF
int	nevent
MonitorElement *	Ntowers_vs_ieta
MonitorElement *	numBadCellsEcal_EB
MonitorElement *	numBadCellsEcal_EE
MonitorElement *	numBadCellsHcal_HB
MonitorElement *	numBadCellsHcal_HE
MonitorElement *	numBadCellsHcal_HF
MonitorElement *	numPrbCellsEcal_EB
MonitorElement *	numPrbCellsEcal_EE
MonitorElement *	numPrbCellsHcal_HB
MonitorElement *	numPrbCellsHcal_HE
MonitorElement *	numPrbCellsHcal_HF
MonitorElement *	numRcvCellsEcal_EB
MonitorElement *	numRcvCellsEcal_EE
MonitorElement *	numRcvCellsHcal_HB
MonitorElement *	numRcvCellsHcal_HE
MonitorElement *	numRcvCellsHcal_HF
MonitorElement *	occupancy_map
MonitorElement *	occupancy_vs_ieta
std::string	outputFile_
MonitorElement *	phiMET_HB
MonitorElement *	phiMET_HE
MonitorElement *	phiMET_HF
MonitorElement *	SET_HB
MonitorElement *	SET_HE
MonitorElement *	SET_HF
edm::EDGetTokenT < CaloTowerCollection >	tok_towers_
bool	useAllHistos_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 41 of file CaloTowersAnalyzer.h.

Member Typedef Documentation

typedef math::RhoEtaPhiVector CaloTowersAnalyzer::Vector

private

Definition at line 58 of file CaloTowersAnalyzer.h.

Constructor & Destructor Documentation

CaloTowersAnalyzer::CaloTowersAnalyzer

(

edm::ParameterSet const &

conf

)

Definition at line 7 of file CaloTowersAnalyzer.cc.

References edm::ParameterSet::getUntrackedParameter(), hcalselector_, outputFile_, AlCaHLTBitMon_QueryRunRegistry::string, tok_towers_, and useAllHistos_.

    : hcalDDDRecConstantsToken_{esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>()} {
  tok_towers_ = consumes<CaloTowerCollection>(conf.getUntrackedParameter<edm::InputTag>("CaloTowerCollectionLabel"));
  // DQM ROOT output

  outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");

  hcalselector_ = conf.getUntrackedParameter<std::string>("hcalselector", "all");

  useAllHistos_ = conf.getUntrackedParameter<bool>("useAllHistos", false);
}

CaloTowersAnalyzer::~CaloTowersAnalyzer ( )

override

Definition at line 452 of file CaloTowersAnalyzer.cc.

{}

Member Function Documentation

void CaloTowersAnalyzer::analyze ( edm::Event const &  event, edm::EventSetup const &    )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 454 of file CaloTowersAnalyzer.cc.

References DetId::det(), DetId::Ecal, EcalBarrel, EcalEndcap, emEnergyTiming_HB, emEnergyTiming_HE, emEnergyTiming_HF, emEnergyTiming_High_HB, emEnergyTiming_Low_HB, emEnergyTiming_Low_HE, emEnergyTiming_profile_HB, emEnergyTiming_profile_HE, emEnergyTiming_profile_HF, emEnergyTiming_profile_High_HB, emEnergyTiming_profile_Low_HB, emEnergyTiming_profile_Low_HE, emTiming_HB, emTiming_HE, emTiming_HF, etaMax, etaMin, dqm::impl::MonitorElement::Fill(), hadEnergyTiming_HB, hadEnergyTiming_HE, hadEnergyTiming_HF, hadEnergyTiming_High_HB, hadEnergyTiming_Low_HB, hadEnergyTiming_Low_HE, hadEnergyTiming_Low_HF, hadEnergyTiming_profile_HB, hadEnergyTiming_profile_HE, hadEnergyTiming_profile_HF, hadEnergyTiming_profile_High_HB, hadEnergyTiming_profile_Low_HB, hadEnergyTiming_profile_Low_HE, hadEnergyTiming_profile_Low_HF, hadTiming_HB, hadTiming_HE, hadTiming_HF, DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, CaloTowerDetId::ieta(), isub, mapEnergy_HB, mapEnergy_HE, mapEnergy_HF, mapEnergy_N, mapEnergyEcal_HB, mapEnergyEcal_HE, mapEnergyEcal_HF, mapEnergyHcal_HB, mapEnergyHcal_HE, mapEnergyHcal_HF, meEnergyEcal_HB, meEnergyEcal_HE, meEnergyEcal_HF, meEnergyEcalTower_HB, meEnergyEcalTower_HE, meEnergyEcalTower_HF, meEnergyHcal_HB, meEnergyHcal_HE, meEnergyHcal_HF, meEnergyHcalTower_HB, meEnergyHcalTower_HE, meEnergyHcalTower_HF, meEnergyHcalvsEcal_HB, meEnergyHcalvsEcal_HE, meEnergyHcalvsEcal_HF, meEnergyHO_HB, meEnergyHO_HE, meEnergyHO_HF, meIphiCaloTower_HFM, meIphiCaloTower_HFP, meIphiHcalTower_HBM, meIphiHcalTower_HBP, meIphiHcalTower_HEM, meIphiHcalTower_HEP, meNumFiredTowers_HB, meNumFiredTowers_HE, meNumFiredTowers_HF, objects.METAnalyzer::met, MET_HB, MET_HE, MET_HF, meTotEnergy_HB, meTotEnergy_HE, meTotEnergy_HF, nevent, Ntowers_vs_ieta, numBadCellsEcal_EB, numBadCellsEcal_EE, numBadCellsHcal_HB, numBadCellsHcal_HE, numBadCellsHcal_HF, numPrbCellsEcal_EB, numPrbCellsEcal_EE, numPrbCellsHcal_HB, numPrbCellsHcal_HE, numPrbCellsHcal_HF, numRcvCellsEcal_EB, numRcvCellsEcal_EE, numRcvCellsHcal_HB, numRcvCellsHcal_HE, numRcvCellsHcal_HF, occupancy_map, phiMET_HB, phiMET_HE, phiMET_HF, SET_HB, SET_HE, SET_HF, mathSSE::sqrt(), HcalDetId::subdet(), DetId::subdetId(), tok_towers_, HLT_FULL_cff::towers, and useAllHistos_.

                                                                             {
  nevent++;

  edm::Handle<CaloTowerCollection> towers;
  event.getByToken(tok_towers_, towers);
  CaloTowerCollection::const_iterator cal;

  double met;
  double phimet;

  // ieta scan
  // double partR  = 0.3;
  // double Rmin   = 9999.;
  // double Econe  = 0.;
  // double Hcone  = 0.;
  // double Ee1    = 0.;
  // double Eh1    = 0.;
  double ieta_MC = 9999;
  double iphi_MC = 9999;

  // HB
  double sumEnergyHcal_HB = 0.;
  double sumEnergyEcal_HB = 0.;
  double sumEnergyHO_HB = 0.;
  Int_t numFiredTowers_HB = 0;
  double metx_HB = 0.;
  double mety_HB = 0.;
  double metz_HB = 0.;
  double sEt_HB = 0.;
  // HE
  double sumEnergyHcal_HE = 0.;
  double sumEnergyEcal_HE = 0.;
  double sumEnergyHO_HE = 0.;
  Int_t numFiredTowers_HE = 0;
  double metx_HE = 0.;
  double mety_HE = 0.;
  double metz_HE = 0.;
  double sEt_HE = 0.;
  // HF
  double sumEnergyHcal_HF = 0.;
  double sumEnergyEcal_HF = 0.;
  double sumEnergyHO_HF = 0.;
  Int_t numFiredTowers_HF = 0;
  double metx_HF = 0.;
  double mety_HF = 0.;
  double metz_HF = 0.;
  double sEt_HF = 0.;

  for (cal = towers->begin(); cal != towers->end(); ++cal) {
    double eE = cal->emEnergy();
    double eH = cal->hadEnergy();
    double eHO = cal->outerEnergy();
    double etaT = cal->eta();
    //    double phiT   = cal->phi();
    double en = cal->energy();
    double etT = cal->et();
    double had_tm = cal->hcalTime();
    double em_tm = cal->ecalTime();

    int numBadEcalCells = cal->numBadEcalCells();
    int numRcvEcalCells = cal->numRecoveredEcalCells();
    int numPrbEcalCells = cal->numProblematicEcalCells();

    int numBadHcalCells = cal->numBadHcalCells();
    int numRcvHcalCells = cal->numRecoveredHcalCells();
    int numPrbHcalCells = cal->numProblematicHcalCells();

    math::RhoEtaPhiVector mom(cal->et(), cal->eta(), cal->phi());
    //  Vector mom  = cal->momentum();

    // cell properties
    CaloTowerDetId idT = cal->id();
    int ieta = idT.ieta();
    int iphi = idT.iphi();

    // ecal:  0 EcalBarrel  1 EcalEndcap
    // hcal:  0 hcalBarrel  1 HcalEndcap  2 HcalForward
    std::vector<int> inEcals(2), inHcals(3);
    unsigned int constitSize = cal->constituentsSize();
    for (unsigned int ic = 0; ic < constitSize; ic++) {
      DetId detId = cal->constituent(ic);
      if (detId.det() == DetId::Ecal) {
        if (detId.subdetId() == EcalBarrel)
          inEcals[0] = 1;
        else if (detId.subdetId() == EcalEndcap)
          inEcals[1] = 1;
      }
      if (detId.det() == DetId::Hcal) {
        if (HcalDetId(detId).subdet() == HcalBarrel)
          inHcals[0] = 1;
        else if (HcalDetId(detId).subdet() == HcalEndcap)
          inHcals[1] = 1;
        else if (HcalDetId(detId).subdet() == HcalForward)
          inHcals[2] = 1;
      }
    }
    // All cell histos are used
    if (inEcals[0]) {
      numBadCellsEcal_EB->Fill(numBadEcalCells);
      numRcvCellsEcal_EB->Fill(numRcvEcalCells);
      numPrbCellsEcal_EB->Fill(numPrbEcalCells);
    }
    if (inEcals[1]) {
      numBadCellsEcal_EE->Fill(numBadEcalCells);
      numRcvCellsEcal_EE->Fill(numRcvEcalCells);
      numPrbCellsEcal_EE->Fill(numPrbEcalCells);
    }

    // Ntowers is used in EndJob, occupancy_map is used for occupancy vs ieta
    Ntowers_vs_ieta->Fill(double(ieta), 1.);
    occupancy_map->Fill(double(ieta), double(iphi));

    if ((isub == 0 || isub == 1) && (fabs(etaT) < etaMax[0] && fabs(etaT) >= etaMin[0])) {
      // All cell histos are used
      numBadCellsHcal_HB->Fill(numBadHcalCells);
      numRcvCellsHcal_HB->Fill(numRcvHcalCells);
      numPrbCellsHcal_HB->Fill(numPrbHcalCells);

      // Map energy histos are not used
      if (useAllHistos_) {
        mapEnergy_HB->Fill(double(ieta), double(iphi), en);
        mapEnergyHcal_HB->Fill(double(ieta), double(iphi), eH);
        mapEnergyEcal_HB->Fill(double(ieta), double(iphi), eE);
      }
      //      std::cout << " e_ecal = " << eE << std::endl;

      if (eH > 0.) {  // iphi of HCAL component of calotower
        if (ieta > 0)
          meIphiHcalTower_HBP->Fill(double(iphi));
        else
          meIphiHcalTower_HBM->Fill(double(iphi));
      }

      //  simple sums
      sumEnergyHcal_HB += eH;
      sumEnergyEcal_HB += eE;
      sumEnergyHO_HB += eHO;

      numFiredTowers_HB++;

      // Not used
      if (useAllHistos_) {
        meEnergyEcalTower_HB->Fill(eE);
        meEnergyHcalTower_HB->Fill(eH);
      }

      // MET, SET & phimet
      //  double  etT = cal->et();
      metx_HB += mom.x();
      mety_HB += mom.y();  // etT * sin(phiT);
      sEt_HB += etT;

      // Timing (all histos are used)
      emTiming_HB->Fill(em_tm);
      hadTiming_HB->Fill(had_tm);

      emEnergyTiming_Low_HB->Fill(eE, em_tm);
      emEnergyTiming_HB->Fill(eE, em_tm);
      emEnergyTiming_High_HB->Fill(eE, em_tm);
      emEnergyTiming_profile_Low_HB->Fill(eE, em_tm);
      emEnergyTiming_profile_HB->Fill(eE, em_tm);
      emEnergyTiming_profile_High_HB->Fill(eE, em_tm);

      hadEnergyTiming_Low_HB->Fill(eH, had_tm);
      hadEnergyTiming_HB->Fill(eH, had_tm);
      hadEnergyTiming_High_HB->Fill(eH, had_tm);
      hadEnergyTiming_profile_Low_HB->Fill(eH, had_tm);
      hadEnergyTiming_profile_HB->Fill(eH, had_tm);
      hadEnergyTiming_profile_High_HB->Fill(eH, had_tm);
    }

    if ((isub == 0 || isub == 2) && (fabs(etaT) < etaMax[1] && fabs(etaT) >= etaMin[1])) {
      // All cell histos are used
      numBadCellsHcal_HE->Fill(numBadHcalCells);
      numRcvCellsHcal_HE->Fill(numRcvHcalCells);
      numPrbCellsHcal_HE->Fill(numPrbHcalCells);

      // Map energy histos are not used
      if (useAllHistos_) {
        mapEnergy_HE->Fill(double(ieta), double(iphi), en);
        mapEnergyHcal_HE->Fill(double(ieta), double(iphi), eH);
        mapEnergyEcal_HE->Fill(double(ieta), double(iphi), eE);
      }
      //      std::cout << " e_ecal = " << eE << std::endl;

      if (eH > 0.) {  // iphi of HCAL component of calotower
        if (ieta > 0)
          meIphiHcalTower_HEP->Fill(double(iphi));
        else
          meIphiHcalTower_HEM->Fill(double(iphi));
      }

      //  simple sums
      sumEnergyHcal_HE += eH;
      sumEnergyEcal_HE += eE;
      sumEnergyHO_HE += eHO;

      numFiredTowers_HE++;

      // Not used
      if (useAllHistos_) {
        meEnergyEcalTower_HE->Fill(eE);
        meEnergyHcalTower_HE->Fill(eH);
      }
      // MET, SET & phimet
      //  double  etT = cal->et();
      metx_HE += mom.x();
      mety_HE += mom.y();  // etT * sin(phiT);
      sEt_HE += etT;

      // Timing (all histos are used)
      emTiming_HE->Fill(em_tm);
      hadTiming_HE->Fill(had_tm);

      emEnergyTiming_Low_HE->Fill(eE, em_tm);
      emEnergyTiming_HE->Fill(eE, em_tm);
      emEnergyTiming_profile_Low_HE->Fill(eE, em_tm);
      emEnergyTiming_profile_HE->Fill(eE, em_tm);

      hadEnergyTiming_Low_HE->Fill(eH, had_tm);
      hadEnergyTiming_HE->Fill(eH, had_tm);
      hadEnergyTiming_profile_Low_HE->Fill(eH, had_tm);
      hadEnergyTiming_profile_HE->Fill(eH, had_tm);
    }

    if ((isub == 0 || isub == 3) && (fabs(etaT) < etaMax[2] && fabs(etaT) >= etaMin[2])) {
      // All cell histos are used
      numBadCellsHcal_HF->Fill(numBadHcalCells);
      numRcvCellsHcal_HF->Fill(numRcvHcalCells);
      numPrbCellsHcal_HF->Fill(numPrbHcalCells);

      // Map energy histos are not used
      if (useAllHistos_) {
        mapEnergy_HF->Fill(double(ieta), double(iphi), en);
        mapEnergyHcal_HF->Fill(double(ieta), double(iphi), eH);
        mapEnergyEcal_HF->Fill(double(ieta), double(iphi), eE);
      }
      //      std::cout << " e_ecal = " << eE << std::endl;

      if (ieta > 0)
        meIphiCaloTower_HFP->Fill(double(iphi));
      else
        meIphiCaloTower_HFM->Fill(double(iphi));

      //  simple sums
      sumEnergyHcal_HF += eH;
      sumEnergyEcal_HF += eE;
      sumEnergyHO_HF += eHO;

      numFiredTowers_HF++;

      // Not used
      if (useAllHistos_) {
        meEnergyEcalTower_HF->Fill(eE);
        meEnergyHcalTower_HF->Fill(eH);
      }
      // MET, SET & phimet
      //  double  etT = cal->et();
      metx_HF += mom.x();
      mety_HF += mom.y();  // etT * sin(phiT);
      sEt_HF += etT;

      // Timing (all histos are used)
      emTiming_HF->Fill(em_tm);
      hadTiming_HF->Fill(had_tm);
      emEnergyTiming_HF->Fill(eE, em_tm);
      emEnergyTiming_profile_HF->Fill(eE, em_tm);

      hadEnergyTiming_Low_HF->Fill(eH, had_tm);
      hadEnergyTiming_HF->Fill(eH, had_tm);
      hadEnergyTiming_profile_Low_HF->Fill(eH, had_tm);
      hadEnergyTiming_profile_HF->Fill(eH, had_tm);
    }

  }  // end of Towers cycle

  // These are the six single pion histos; only the second set is used

  mapEnergy_N->Fill(double(ieta_MC), double(iphi_MC), 1.);

  if (isub == 0 || isub == 1) {
    met = sqrt(metx_HB * metx_HB + mety_HB * mety_HB);
    Vector metv(metx_HB, mety_HB, metz_HB);
    phimet = metv.phi();

    // Five oldest drawn histos first; the rest are not used
    meEnergyHcal_HB->Fill(sumEnergyHcal_HB);
    meEnergyEcal_HB->Fill(sumEnergyEcal_HB);
    meNumFiredTowers_HB->Fill(numFiredTowers_HB);
    MET_HB->Fill(met);
    SET_HB->Fill(sEt_HB);

    if (useAllHistos_) {
      meEnergyHcalvsEcal_HB->Fill(sumEnergyEcal_HB, sumEnergyHcal_HB);
      meEnergyHO_HB->Fill(sumEnergyHO_HB);
      meTotEnergy_HB->Fill(sumEnergyHcal_HB + sumEnergyEcal_HB + sumEnergyHO_HB);
      phiMET_HB->Fill(phimet);
    }
  }

  if (isub == 0 || isub == 2) {
    met = sqrt(metx_HE * metx_HE + mety_HE * mety_HE);
    Vector metv(metx_HE, mety_HE, metz_HE);
    phimet = metv.phi();

    // Five oldest drawn histos first; the rest are not used
    meEnergyHcal_HE->Fill(sumEnergyHcal_HE);
    meEnergyEcal_HE->Fill(sumEnergyEcal_HE);
    meNumFiredTowers_HE->Fill(numFiredTowers_HE);
    MET_HE->Fill(met);
    SET_HE->Fill(sEt_HE);

    if (useAllHistos_) {
      meEnergyHcalvsEcal_HE->Fill(sumEnergyEcal_HE, sumEnergyHcal_HE);
      meEnergyHO_HE->Fill(sumEnergyHO_HE);
      meTotEnergy_HE->Fill(sumEnergyHcal_HE + sumEnergyEcal_HE + sumEnergyHO_HE);
      phiMET_HE->Fill(phimet);
    }
  }

  if (isub == 0 || isub == 3) {
    met = sqrt(metx_HF * metx_HF + mety_HF * mety_HF);
    Vector metv(metx_HF, mety_HF, metz_HF);
    phimet = metv.phi();

    // Five oldest drawn histos first; the rest are not used
    meEnergyHcal_HF->Fill(sumEnergyHcal_HF);
    meEnergyEcal_HF->Fill(sumEnergyEcal_HF);
    meNumFiredTowers_HF->Fill(numFiredTowers_HF);
    MET_HF->Fill(met);
    SET_HF->Fill(sEt_HF);

    if (useAllHistos_) {
      meEnergyHcalvsEcal_HF->Fill(sumEnergyEcal_HF, sumEnergyHcal_HF);
      meEnergyHO_HF->Fill(sumEnergyHO_HF);
      meTotEnergy_HF->Fill(sumEnergyHcal_HF + sumEnergyEcal_HF + sumEnergyHO_HF);
      phiMET_HF->Fill(phimet);
    }
  }
}

void CaloTowersAnalyzer::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 50 of file CaloTowersAnalyzer.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookProfile(), emean_vs_ieta_E, emean_vs_ieta_E1, emean_vs_ieta_EH, emean_vs_ieta_EH1, emean_vs_ieta_H, emean_vs_ieta_H1, emEnergyTiming_HB, emEnergyTiming_HE, emEnergyTiming_HF, emEnergyTiming_High_HB, emEnergyTiming_Low_HB, emEnergyTiming_Low_HE, emEnergyTiming_profile_HB, emEnergyTiming_profile_HE, emEnergyTiming_profile_HF, emEnergyTiming_profile_High_HB, emEnergyTiming_profile_Low_HB, emEnergyTiming_profile_Low_HE, emTiming_HB, emTiming_HE, emTiming_HF, etaMax, etaMin, hadEnergyTiming_HB, hadEnergyTiming_HE, hadEnergyTiming_HF, hadEnergyTiming_High_HB, hadEnergyTiming_Low_HB, hadEnergyTiming_Low_HE, hadEnergyTiming_Low_HF, hadEnergyTiming_profile_HB, hadEnergyTiming_profile_HE, hadEnergyTiming_profile_HF, hadEnergyTiming_profile_High_HB, hadEnergyTiming_profile_Low_HB, hadEnergyTiming_profile_Low_HE, hadEnergyTiming_profile_Low_HF, hadTiming_HB, hadTiming_HE, hadTiming_HF, hcalselector_, timingPdfMaker::histo, ieta_bins_, ieta_max_, ieta_min_, iphi_bins_, iphi_max_, iphi_min_, isub, mapEnergy_E, mapEnergy_EH, mapEnergy_H, mapEnergy_HB, mapEnergy_HE, mapEnergy_HF, mapEnergy_N, mapEnergyEcal_HB, mapEnergyEcal_HE, mapEnergyEcal_HF, mapEnergyHcal_HB, mapEnergyHcal_HE, mapEnergyHcal_HF, meEnergyEcal_HB, meEnergyEcal_HE, meEnergyEcal_HF, meEnergyEcalTower_HB, meEnergyEcalTower_HE, meEnergyEcalTower_HF, meEnergyHcal_HB, meEnergyHcal_HE, meEnergyHcal_HF, meEnergyHcalTower_HB, meEnergyHcalTower_HE, meEnergyHcalTower_HF, meEnergyHcalvsEcal_HB, meEnergyHcalvsEcal_HE, meEnergyHcalvsEcal_HF, meEnergyHO_HB, meEnergyHO_HE, meEnergyHO_HF, meIphiCaloTower_HFM, meIphiCaloTower_HFP, meIphiHcalTower_HBM, meIphiHcalTower_HBP, meIphiHcalTower_HEM, meIphiHcalTower_HEP, meNumFiredTowers_HB, meNumFiredTowers_HE, meNumFiredTowers_HF, MET_HB, MET_HE, MET_HF, meTotEnergy_HB, meTotEnergy_HE, meTotEnergy_HF, Ntowers_vs_ieta, numBadCellsEcal_EB, numBadCellsEcal_EE, numBadCellsHcal_HB, numBadCellsHcal_HE, numBadCellsHcal_HF, numPrbCellsEcal_EB, numPrbCellsEcal_EE, numPrbCellsHcal_HB, numPrbCellsHcal_HE, numPrbCellsHcal_HF, numRcvCellsEcal_EB, numRcvCellsEcal_EE, numRcvCellsHcal_HB, numRcvCellsHcal_HE, numRcvCellsHcal_HF, occupancy_map, occupancy_vs_ieta, outputFile_, phiMET_HB, phiMET_HE, phiMET_HF, SET_HB, SET_HE, SET_HF, dqm::implementation::NavigatorBase::setCurrentFolder(), and useAllHistos_.

                                                                 {
  etaMin[0] = 0.;
  etaMax[0] = 1.4;
  etaMin[1] = 1.4;
  etaMax[1] = 2.9;
  etaMin[2] = 2.9;
  etaMax[2] = 5.2;

  isub = 0;
  if (hcalselector_ == "HB")
    isub = 1;
  if (hcalselector_ == "HE")
    isub = 2;
  if (hcalselector_ == "HF")
    isub = 3;

  if (!outputFile_.empty()) {
    edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will be saved to '" << outputFile_.c_str() << "'";
  } else {
    edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will NOT be saved";
  }

  Char_t histo[100];

  ibooker.setCurrentFolder("CaloTowersD/CaloTowersTask");

  // These two histos are not drawn by our macros, but they are used
  // in the EndJob for norms and such so I am leaving them alone for now
  //-------------------------------------------------------------------------------------------
  sprintf(histo, "Ntowers_per_event_vs_ieta");
  Ntowers_vs_ieta = ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_);

  sprintf(histo, "CaloTowersTask_map_Nentries");
  mapEnergy_N = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
  //-------------------------------------------------------------------------------------------

  // These the single pion scan histos
  //-------------------------------------------------------------------------------------------
  // The first three are not used
  if (useAllHistos_) {
    sprintf(histo, "emean_vs_ieta_E");
    emean_vs_ieta_E = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
    sprintf(histo, "emean_vs_ieta_H");
    emean_vs_ieta_H = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
    sprintf(histo, "emean_vs_ieta_EH");
    emean_vs_ieta_EH = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
  }
  // These are drawn
  sprintf(histo, "emean_vs_ieta_E1");
  emean_vs_ieta_E1 = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
  sprintf(histo, "emean_vs_ieta_H1");
  emean_vs_ieta_H1 = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
  sprintf(histo, "emean_vs_ieta_EH1");
  emean_vs_ieta_EH1 = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
  //-------------------------------------------------------------------------------------------

  // Map energy histos are not drawn
  if (useAllHistos_) {
    sprintf(histo, "CaloTowersTask_map_energy_E");
    mapEnergy_E = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
    sprintf(histo, "CaloTowersTask_map_energy_H");
    mapEnergy_H = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
    sprintf(histo, "CaloTowersTask_map_energy_EH");
    mapEnergy_EH = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
  }

  // All ECAL cell histos are used
  // XXX: ECAL 0-25 [0-26, 26 bins]   HCAL 0-4 [0-5, 5 bins]
  sprintf(histo, "number_of_bad_cells_Ecal_EB");
  numBadCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
  sprintf(histo, "number_of_bad_cells_Ecal_EE");
  numBadCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26);
  sprintf(histo, "number_of_recovered_cells_Ecal_EB");
  numRcvCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
  sprintf(histo, "number_of_recovered_cells_Ecal_EE");
  numRcvCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26);
  sprintf(histo, "number_of_problematic_cells_Ecal_EB");
  numPrbCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
  sprintf(histo, "number_of_problematic_cells_Ecal_EE");
  numPrbCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26);

  // Occupancy vs. ieta is drawn, occupancy map is needed to draw it
  sprintf(histo, "CaloTowersTask_map_occupancy");
  occupancy_map = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);

  sprintf(histo, "CaloTowersTask_occupancy_vs_ieta");
  occupancy_vs_ieta = ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_);

  if (isub == 1 || isub == 0) {
    // All cell histos are used
    sprintf(histo, "number_of_bad_cells_Hcal_HB");
    numBadCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf(histo, "number_of_recovered_cells_Hcal_HB");
    numRcvCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf(histo, "number_of_problematic_cells_Hcal_HB");
    numPrbCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5);

    // These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL,
    // MET and SET
    //-------------------------------------------------------------------------------------------
    sprintf(histo, "CaloTowersTask_energy_HCAL_HB");
    meEnergyHcal_HB = ibooker.book1D(histo, histo, 4100, -200, 8000);

    sprintf(histo, "CaloTowersTask_energy_ECAL_HB");
    meEnergyEcal_HB = ibooker.book1D(histo, histo, 3100, -200, 6000);

    sprintf(histo, "CaloTowersTask_number_of_fired_towers_HB");
    meNumFiredTowers_HB = ibooker.book1D(histo, histo, 1000, 0, 2000);

    sprintf(histo, "CaloTowersTask_MET_HB");
    MET_HB = ibooker.book1D(histo, histo, 3000, 0., 3000.);

    sprintf(histo, "CaloTowersTask_SET_HB");
    SET_HB = ibooker.book1D(histo, histo, 8000, 0., 8000.);
    //-------------------------------------------------------------------------------------------

    // Timing histos and profiles -- all six are necessary
    //-------------------------------------------------------------------------------------------
    sprintf(histo, "CaloTowersTask_EM_Timing_HB");
    emTiming_HB = ibooker.book1D(histo, histo, 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_HAD_Timing_HB");
    hadTiming_HB = ibooker.book1D(histo, histo, 70, -48., 92.);

    // Energy-Timing histos are divided into low, medium and high to reduce
    // memory usage EM
    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_Low_HB");
    emEnergyTiming_Low_HB = ibooker.book2D(histo, histo, 40, 0., 40., 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_HB");
    emEnergyTiming_HB = ibooker.book2D(histo, histo, 200, 0., 400., 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_High_HB");
    emEnergyTiming_High_HB = ibooker.book2D(histo, histo, 200, 0., 3000., 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_profile_Low_HB");
    emEnergyTiming_profile_Low_HB = ibooker.bookProfile(histo, histo, 40, 0., 40., 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_profile_HB");
    emEnergyTiming_profile_HB = ibooker.bookProfile(histo, histo, 200, 0., 400., 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_profile_High_HB");
    emEnergyTiming_profile_High_HB = ibooker.bookProfile(histo, histo, 200, 0., 3000., 110, -120., 100.);

    // HAD
    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_Low_HB");
    hadEnergyTiming_Low_HB = ibooker.book2D(histo, histo, 40, 0., 40., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_HB");
    hadEnergyTiming_HB = ibooker.book2D(histo, histo, 100, 0., 200., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_High_HB");
    hadEnergyTiming_High_HB = ibooker.book2D(histo, histo, 300, 0., 3000., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HB");
    hadEnergyTiming_profile_Low_HB = ibooker.bookProfile(histo, histo, 40, 0., 40., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_profile_HB");
    hadEnergyTiming_profile_HB = ibooker.bookProfile(histo, histo, 100, 0., 200., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_profile_High_HB");
    hadEnergyTiming_profile_High_HB = ibooker.bookProfile(histo, histo, 300, 0., 3000., 70, -48., 92.);
    //-------------------------------------------------------------------------------------------

    sprintf(histo, "CaloTowersTask_Iphi_HCAL_component_of_tower_HBP");
    meIphiHcalTower_HBP = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_);
    sprintf(histo, "CaloTowersTask_Iphi_HCAL_component_of_tower_HBM");
    meIphiHcalTower_HBM = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_);

    //-------------------------------------------------------------------------------------------

    // Everything else is not drawn
    if (useAllHistos_) {
      sprintf(histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HB");
      meEnergyHcalvsEcal_HB = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);

      sprintf(histo, "CaloTowersTask_energy_OUTER_HB");
      meEnergyHO_HB = ibooker.book1D(histo, histo, 1640, -200, 8000);

      sprintf(histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HB");
      meEnergyEcalTower_HB = ibooker.book1D(histo, histo, 440, -200, 2000);

      sprintf(histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HB");
      meEnergyHcalTower_HB = ibooker.book1D(histo, histo, 440, -200, 2000);

      sprintf(histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HB");
      meTotEnergy_HB = ibooker.book1D(histo, histo, 400, 0., 2000.);

      sprintf(histo, "CaloTowersTask_map_energy_HB");
      mapEnergy_HB = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf(histo, "CaloTowersTask_map_energy_HCAL_HB");
      mapEnergyHcal_HB =
          ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf(histo, "CaloTowersTask_map_energy_ECAL_HB");
      mapEnergyEcal_HB =
          ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);

      sprintf(histo, "CaloTowersTask_phi_MET_HB");
      phiMET_HB = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898);
    }
  }

  if (isub == 2 || isub == 0) {
    // All cell histos are used
    sprintf(histo, "number_of_bad_cells_Hcal_HE");
    numBadCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf(histo, "number_of_recovered_cells_Hcal_HE");
    numRcvCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf(histo, "number_of_problematic_cells_Hcal_HE");
    numPrbCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5);

    // These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL,
    // MET and SET
    //-------------------------------------------------------------------------------------------
    sprintf(histo, "CaloTowersTask_energy_HCAL_HE");
    meEnergyHcal_HE = ibooker.book1D(histo, histo, 1240, -200, 6000);

    sprintf(histo, "CaloTowersTask_energy_ECAL_HE");
    meEnergyEcal_HE = ibooker.book1D(histo, histo, 1240, -200, 6000);

    sprintf(histo, "CaloTowersTask_number_of_fired_towers_HE");
    meNumFiredTowers_HE = ibooker.book1D(histo, histo, 1000, 0, 2000);

    sprintf(histo, "CaloTowersTask_MET_HE");
    MET_HE = ibooker.book1D(histo, histo, 1000, 0., 1000.);

    sprintf(histo, "CaloTowersTask_SET_HE");
    SET_HE = ibooker.book1D(histo, histo, 2000, 0., 2000.);
    //-------------------------------------------------------------------------------------------

    // Timing histos and profiles -- all six are necessary
    //-------------------------------------------------------------------------------------------
    sprintf(histo, "CaloTowersTask_EM_Timing_HE");
    emTiming_HE = ibooker.book1D(histo, histo, 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_HAD_Timing_HE");
    hadTiming_HE = ibooker.book1D(histo, histo, 70, -48., 92.);

    // Energy-Timing histos are divided into low and normal to reduce memory
    // usage EM
    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_Low_HE");
    emEnergyTiming_Low_HE = ibooker.book2D(histo, histo, 160, 0., 160., 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_HE");
    emEnergyTiming_HE = ibooker.book2D(histo, histo, 200, 0., 800., 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_profile_Low_HE");
    emEnergyTiming_profile_Low_HE = ibooker.bookProfile(histo, histo, 160, 0., 160., 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_profile_HE");
    emEnergyTiming_profile_HE = ibooker.bookProfile(histo, histo, 200, 0., 800., 110, -120., 100.);

    // HAD
    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_Low_HE");
    hadEnergyTiming_Low_HE = ibooker.book2D(histo, histo, 160, 0., 160., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_HE");
    hadEnergyTiming_HE = ibooker.book2D(histo, histo, 200, 0., 800., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HE");
    hadEnergyTiming_profile_Low_HE = ibooker.bookProfile(histo, histo, 160, 0., 160., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_profile_HE");
    hadEnergyTiming_profile_HE = ibooker.bookProfile(histo, histo, 200, 0., 800., 70, -48., 92.);
    //-------------------------------------------------------------------------------------------

    sprintf(histo, "CaloTowersTask_Iphi_HCAL_component_of_tower_HEP");
    meIphiHcalTower_HEP = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_);
    sprintf(histo, "CaloTowersTask_Iphi_HCAL_component_of_tower_HEM");
    meIphiHcalTower_HEM = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_);

    //-------------------------------------------------------------------------------------------

    // Everything else is not drawn
    if (useAllHistos_) {
      sprintf(histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HE");
      meEnergyHcalvsEcal_HE = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);

      sprintf(histo, "CaloTowersTask_energy_OUTER_HE");
      meEnergyHO_HE = ibooker.book1D(histo, histo, 440, -200, 2000);

      sprintf(histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HE");
      meEnergyEcalTower_HE = ibooker.book1D(histo, histo, 1100, -200, 2000);

      sprintf(histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HE");
      meEnergyHcalTower_HE = ibooker.book1D(histo, histo, 1100, -200, 2000);

      sprintf(histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HE");
      meTotEnergy_HE = ibooker.book1D(histo, histo, 400, 0., 2000.);

      sprintf(histo, "CaloTowersTask_map_energy_HE");
      mapEnergy_HE = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf(histo, "CaloTowersTask_map_energy_HCAL_HE");
      mapEnergyHcal_HE =
          ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf(histo, "CaloTowersTask_map_energy_ECAL_HE");
      mapEnergyEcal_HE =
          ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);

      sprintf(histo, "CaloTowersTask_phi_MET_HE");
      phiMET_HE = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898);
    }
  }

  if (isub == 3 || isub == 0) {
    // All cell histos are used
    sprintf(histo, "number_of_bad_cells_Hcal_HF");
    numBadCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf(histo, "number_of_recovered_cells_Hcal_HF");
    numRcvCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf(histo, "number_of_problematic_cells_Hcal_HF");
    numPrbCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5);

    // These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL,
    // MET and SET
    //-------------------------------------------------------------------------------------------
    sprintf(histo, "CaloTowersTask_energy_HCAL_HF");
    meEnergyHcal_HF = ibooker.book1D(histo, histo, 5001, -20., 1000000.);

    sprintf(histo, "CaloTowersTask_energy_ECAL_HF");
    meEnergyEcal_HF = ibooker.book1D(histo, histo, 3501, -20., 70000.);

    sprintf(histo, "CaloTowersTask_number_of_fired_towers_HF");
    meNumFiredTowers_HF = ibooker.book1D(histo, histo, 1000, 0., 2000.);

    sprintf(histo, "CaloTowersTask_MET_HF");
    MET_HF = ibooker.book1D(histo, histo, 500, 0., 500.);

    sprintf(histo, "CaloTowersTask_SET_HF");
    SET_HF = ibooker.book1D(histo, histo, 500, 0., 5000.);
    //-------------------------------------------------------------------------------------------

    // Timing histos and profiles -- all six are necessary
    //-------------------------------------------------------------------------------------------
    sprintf(histo, "CaloTowersTask_EM_Timing_HF");
    emTiming_HF = ibooker.book1D(histo, histo, 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_HAD_Timing_HF");
    hadTiming_HF = ibooker.book1D(histo, histo, 70, -48., 92.);

    // EM
    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_HF");
    emEnergyTiming_HF = ibooker.book2D(histo, histo, 150, 0., 300., 110, -120., 100.);

    sprintf(histo, "CaloTowersTask_EM_Energy_Timing_profile_HF");
    emEnergyTiming_profile_HF = ibooker.bookProfile(histo, histo, 150, 0., 300., 110, -120., 100.);

    // HAD (requires two different sets of histograms to lower RAM usage)
    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_Low_HF");
    hadEnergyTiming_Low_HF = ibooker.book2D(histo, histo, 40, 0., 40., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_HF");
    hadEnergyTiming_HF = ibooker.book2D(histo, histo, 200, 0., 600., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HF");
    hadEnergyTiming_profile_Low_HF = ibooker.bookProfile(histo, histo, 40, 0., 40., 70, -48., 92.);

    sprintf(histo, "CaloTowersTask_HAD_Energy_Timing_profile_HF");
    hadEnergyTiming_profile_HF = ibooker.bookProfile(histo, histo, 200, 0., 600., 70, -48., 92.);
    //-------------------------------------------------------------------------------------------

    sprintf(histo, "CaloTowersTask_Iphi_tower_HFP");
    meIphiCaloTower_HFP = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_);
    sprintf(histo, "CaloTowersTask_Iphi_tower_HFM");
    meIphiCaloTower_HFM = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_);

    //-------------------------------------------------------------------------------------------

    // Everything else is not drawn
    if (useAllHistos_) {
      sprintf(histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HF");
      meEnergyHcalvsEcal_HF = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);

      sprintf(histo, "CaloTowersTask_energy_OUTER_HF");
      meEnergyHO_HF = ibooker.book1D(histo, histo, 440, -200, 2000);

      sprintf(histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HF");
      meEnergyEcalTower_HF = ibooker.book1D(histo, histo, 440, -200, 2000);

      sprintf(histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HF");
      meEnergyHcalTower_HF = ibooker.book1D(histo, histo, 440, -200, 2000);

      sprintf(histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HF");
      meTotEnergy_HF = ibooker.book1D(histo, histo, 400, 0., 2000.);

      sprintf(histo, "CaloTowersTask_map_energy_HF");
      mapEnergy_HF = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf(histo, "CaloTowersTask_map_energy_HCAL_HF");
      mapEnergyHcal_HF =
          ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf(histo, "CaloTowersTask_map_energy_ECAL_HF");
      mapEnergyEcal_HF =
          ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);

      sprintf(histo, "CaloTowersTask_phi_MET_HF");
      phiMET_HF = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898);
    }
  }
}

void CaloTowersAnalyzer::dqmBeginRun ( const edm::Run &  run, const edm::EventSetup &  es  )

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 19 of file CaloTowersAnalyzer.cc.

References edm::EventSetup::getData(), HcalDDDRecConstants::getEtaRange(), HcalDDDRecConstants::getNPhi(), hcalDDDRecConstantsToken_, ieta_bins_, ieta_max_, ieta_min_, iphi_bins_, iphi_max_, and iphi_min_.

                                                                               {
  HcalDDDRecConstants const &hcons = es.getData(hcalDDDRecConstantsToken_);

  // Get Phi segmentation from geometry, use the max phi number so that all iphi
  // values are included.

  int NphiMax = hcons.getNPhi(0);

  NphiMax = (hcons.getNPhi(1) > NphiMax ? hcons.getNPhi(1) : NphiMax);
  NphiMax = (hcons.getNPhi(2) > NphiMax ? hcons.getNPhi(2) : NphiMax);
  NphiMax = (hcons.getNPhi(3) > NphiMax ? hcons.getNPhi(3) : NphiMax);

  // Center the iphi bins on the integers
  iphi_min_ = 0.5;
  iphi_max_ = NphiMax + 0.5;
  iphi_bins_ = (int)(iphi_max_ - iphi_min_);

  // Retain classic behavior, all plots have same ieta range.

  int iEtaMax = (hcons.getEtaRange(0).second > hcons.getEtaRange(1).second ? hcons.getEtaRange(0).second
                                                                           : hcons.getEtaRange(1).second);
  iEtaMax = (iEtaMax > hcons.getEtaRange(2).second ? iEtaMax : hcons.getEtaRange(2).second);
  iEtaMax = (iEtaMax > hcons.getEtaRange(3).second ? iEtaMax : hcons.getEtaRange(3).second);

  // Give an empty bin around the subdet ieta range to make it clear that all
  // ieta rings have been included
  ieta_min_ = -iEtaMax - 1.5;
  ieta_max_ = iEtaMax + 1.5;
  ieta_bins_ = (int)(ieta_max_ - ieta_min_);
}

double CaloTowersAnalyzer::dR ( double  eta1, double  phi1, double  eta2, double  phi2  )

private

Definition at line 795 of file CaloTowersAnalyzer.cc.

References HLT_FULL_cff::eta1, PI, mathSSE::sqrt(), and createJobs::tmp.

                                                                                {
  double PI = 3.1415926535898;
  double deltaphi = phi1 - phi2;
  if (phi2 > phi1) {
    deltaphi = phi2 - phi1;
  }
  if (deltaphi > PI) {
    deltaphi = 2. * PI - deltaphi;
  }
  double deltaeta = eta2 - eta1;
  double tmp = sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
  return tmp;
}

Member Data Documentation

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_E

private

Definition at line 79 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_E1

private

Definition at line 83 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_EH

private

Definition at line 81 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_EH1

private

Definition at line 85 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_H

private

Definition at line 80 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_H1

private

Definition at line 84 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HB

private

Definition at line 140 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HE

private

Definition at line 178 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HF

private

Definition at line 210 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_High_HB

private

Definition at line 141 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_Low_HB

private

Definition at line 139 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_Low_HE

private

Definition at line 177 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HB

private

Definition at line 143 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HE

private

Definition at line 180 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HF

private

Definition at line 211 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_High_HB

private

Definition at line 144 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_Low_HB

private

Definition at line 142 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_Low_HE

private

Definition at line 179 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emTiming_HB

private

Definition at line 136 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emTiming_HE

private

Definition at line 174 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emTiming_HF

private

Definition at line 208 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

double CaloTowersAnalyzer::etaMax[3]

private

Definition at line 69 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

double CaloTowersAnalyzer::etaMin[3]

private

Definition at line 70 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HB

private

Definition at line 147 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HE

private

Definition at line 183 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HF

private

Definition at line 214 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_High_HB

private

Definition at line 148 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HB

private

Definition at line 146 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HE

private

Definition at line 182 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HF

private

Definition at line 213 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HB

private

Definition at line 150 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HE

private

Definition at line 185 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HF

private

Definition at line 216 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_High_HB

private

Definition at line 151 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HB

private

Definition at line 149 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HE

private

Definition at line 184 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HF

private

Definition at line 215 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadTiming_HB

private

Definition at line 137 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadTiming_HE

private

Definition at line 175 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadTiming_HF

private

Definition at line 209 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> CaloTowersAnalyzer::hcalDDDRecConstantsToken_

private

Definition at line 61 of file CaloTowersAnalyzer.h.

Referenced by dqmBeginRun().

std::string CaloTowersAnalyzer::hcalselector_

private

Definition at line 54 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and CaloTowersAnalyzer().

int CaloTowersAnalyzer::ieta_bins_

private

Definition at line 75 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float CaloTowersAnalyzer::ieta_max_

private

Definition at line 76 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float CaloTowersAnalyzer::ieta_min_

private

Definition at line 76 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int CaloTowersAnalyzer::imc

private

Definition at line 66 of file CaloTowersAnalyzer.h.

int CaloTowersAnalyzer::iphi_bins_

private

Definition at line 72 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float CaloTowersAnalyzer::iphi_max_

private

Definition at line 73 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float CaloTowersAnalyzer::iphi_min_

private

Definition at line 73 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int CaloTowersAnalyzer::isub

private

Definition at line 63 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_E

private

Definition at line 92 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_EH

private

Definition at line 94 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_H

private

Definition at line 93 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HB

private

Definition at line 129 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HE

private

Definition at line 167 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HF

private

Definition at line 201 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_N

private

Definition at line 95 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HB

private

Definition at line 130 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HE

private

Definition at line 168 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HF

private

Definition at line 202 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HB

private

Definition at line 131 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HE

private

Definition at line 169 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HF

private

Definition at line 203 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::string CaloTowersAnalyzer::mc_

private

Definition at line 55 of file CaloTowersAnalyzer.h.

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HB

private

Definition at line 118 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HE

private

Definition at line 156 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HF

private

Definition at line 190 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HB

private

Definition at line 122 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HE

private

Definition at line 160 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HF

private

Definition at line 194 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HB

private

Definition at line 119 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HE

private

Definition at line 157 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HF

private

Definition at line 191 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HB

private

Definition at line 123 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HE

private

Definition at line 161 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HF

private

Definition at line 195 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HB

private

Definition at line 116 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HE

private

Definition at line 154 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HF

private

Definition at line 188 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HB

private

Definition at line 117 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HE

private

Definition at line 155 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HF

private

Definition at line 189 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiCaloTower_HFM

private

Definition at line 199 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiCaloTower_HFP

private

Definition at line 198 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiHcalTower_HBM

private

Definition at line 127 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiHcalTower_HBP

private

Definition at line 126 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiHcalTower_HEM

private

Definition at line 165 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiHcalTower_HEP

private

Definition at line 164 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HB

private

Definition at line 120 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HE

private

Definition at line 158 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HF

private

Definition at line 192 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::MET_HB

private

Definition at line 132 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::MET_HE

private

Definition at line 170 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::MET_HF

private

Definition at line 204 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HB

private

Definition at line 124 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HE

private

Definition at line 162 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HF

private

Definition at line 196 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

int CaloTowersAnalyzer::nevent

private

Definition at line 64 of file CaloTowersAnalyzer.h.

Referenced by analyze().

MonitorElement* CaloTowersAnalyzer::Ntowers_vs_ieta

private

Definition at line 87 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsEcal_EB

private

Definition at line 98 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsEcal_EE

private

Definition at line 99 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HB

private

Definition at line 105 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HE

private

Definition at line 106 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HF

private

Definition at line 107 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsEcal_EB

private

Definition at line 102 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsEcal_EE

private

Definition at line 103 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HB

private

Definition at line 111 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HE

private

Definition at line 112 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HF

private

Definition at line 113 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsEcal_EB

private

Definition at line 100 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsEcal_EE

private

Definition at line 101 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HB

private

Definition at line 108 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HE

private

Definition at line 109 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HF

private

Definition at line 110 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::occupancy_map

private

Definition at line 88 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::occupancy_vs_ieta

private

Definition at line 89 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

std::string CaloTowersAnalyzer::outputFile_

private

Definition at line 53 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::phiMET_HB

private

Definition at line 134 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::phiMET_HE

private

Definition at line 172 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::phiMET_HF

private

Definition at line 206 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::SET_HB

private

Definition at line 133 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::SET_HE

private

Definition at line 171 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::SET_HF

private

Definition at line 205 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<CaloTowerCollection> CaloTowersAnalyzer::tok_towers_

private

Definition at line 60 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

bool CaloTowersAnalyzer::useAllHistos_

private

Definition at line 56 of file CaloTowersAnalyzer.h.

Referenced by analyze(), bookHistograms(), and CaloTowersAnalyzer().