CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes
CaloTowersAnalyzer Class Reference

#include <CaloTowersAnalyzer.h>

Inheritance diagram for CaloTowersAnalyzer:
DQMEDAnalyzer edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > > edm::stream::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

virtual void analyze (edm::Event const &e, edm::EventSetup const &c)
 
virtual void beginJob ()
 
virtual void beginRun ()
 
virtual void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
 CaloTowersAnalyzer (edm::ParameterSet const &conf)
 
virtual void endJob ()
 
virtual void endRun ()
 
 ~CaloTowersAnalyzer ()
 
- Public Member Functions inherited from DQMEDAnalyzer
virtual void beginRun (edm::Run const &, edm::EventSetup const &) final
 
virtual void beginStream (edm::StreamID id) final
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer (void)
 
virtual void endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
virtual void endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
uint32_t streamId () const
 
- Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
 EDAnalyzer ()=default
 
- Public Member Functions inherited from edm::stream::EDAnalyzerBase
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDAnalyzerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 
itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Types

typedef math::RhoEtaPhiVector Vector
 

Private Member Functions

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

Private Attributes

MonitorElementemean_vs_ieta_E
 
MonitorElementemean_vs_ieta_E1
 
MonitorElementemean_vs_ieta_EH
 
MonitorElementemean_vs_ieta_EH1
 
MonitorElementemean_vs_ieta_H
 
MonitorElementemean_vs_ieta_H1
 
MonitorElementemEnergyTiming_HB
 
MonitorElementemEnergyTiming_HE
 
MonitorElementemEnergyTiming_HF
 
MonitorElementemEnergyTiming_High_HB
 
MonitorElementemEnergyTiming_Low_HB
 
MonitorElementemEnergyTiming_Low_HE
 
MonitorElementemEnergyTiming_profile_HB
 
MonitorElementemEnergyTiming_profile_HE
 
MonitorElementemEnergyTiming_profile_HF
 
MonitorElementemEnergyTiming_profile_High_HB
 
MonitorElementemEnergyTiming_profile_Low_HB
 
MonitorElementemEnergyTiming_profile_Low_HE
 
MonitorElementemTiming_HB
 
MonitorElementemTiming_HE
 
MonitorElementemTiming_HF
 
double etaMax [3]
 
double etaMin [3]
 
MonitorElementhadEnergyTiming_HB
 
MonitorElementhadEnergyTiming_HE
 
MonitorElementhadEnergyTiming_HF
 
MonitorElementhadEnergyTiming_High_HB
 
MonitorElementhadEnergyTiming_Low_HB
 
MonitorElementhadEnergyTiming_Low_HE
 
MonitorElementhadEnergyTiming_Low_HF
 
MonitorElementhadEnergyTiming_profile_HB
 
MonitorElementhadEnergyTiming_profile_HE
 
MonitorElementhadEnergyTiming_profile_HF
 
MonitorElementhadEnergyTiming_profile_High_HB
 
MonitorElementhadEnergyTiming_profile_Low_HB
 
MonitorElementhadEnergyTiming_profile_Low_HE
 
MonitorElementhadEnergyTiming_profile_Low_HF
 
MonitorElementhadTiming_HB
 
MonitorElementhadTiming_HE
 
MonitorElementhadTiming_HF
 
std::string hcalselector_
 
int imc
 
int isub
 
MonitorElementmapEnergy_E
 
MonitorElementmapEnergy_EH
 
MonitorElementmapEnergy_H
 
MonitorElementmapEnergy_HB
 
MonitorElementmapEnergy_HE
 
MonitorElementmapEnergy_HF
 
MonitorElementmapEnergy_N
 
MonitorElementmapEnergyEcal_HB
 
MonitorElementmapEnergyEcal_HE
 
MonitorElementmapEnergyEcal_HF
 
MonitorElementmapEnergyHcal_HB
 
MonitorElementmapEnergyHcal_HE
 
MonitorElementmapEnergyHcal_HF
 
std::string mc_
 
MonitorElementmeEnergyEcal_HB
 
MonitorElementmeEnergyEcal_HE
 
MonitorElementmeEnergyEcal_HF
 
MonitorElementmeEnergyEcalTower_HB
 
MonitorElementmeEnergyEcalTower_HE
 
MonitorElementmeEnergyEcalTower_HF
 
MonitorElementmeEnergyHcal_HB
 
MonitorElementmeEnergyHcal_HE
 
MonitorElementmeEnergyHcal_HF
 
MonitorElementmeEnergyHcalTower_HB
 
MonitorElementmeEnergyHcalTower_HE
 
MonitorElementmeEnergyHcalTower_HF
 
MonitorElementmeEnergyHcalvsEcal_HB
 
MonitorElementmeEnergyHcalvsEcal_HE
 
MonitorElementmeEnergyHcalvsEcal_HF
 
MonitorElementmeEnergyHO_HB
 
MonitorElementmeEnergyHO_HE
 
MonitorElementmeEnergyHO_HF
 
MonitorElementmeNumFiredTowers_HB
 
MonitorElementmeNumFiredTowers_HE
 
MonitorElementmeNumFiredTowers_HF
 
MonitorElementMET_HB
 
MonitorElementMET_HE
 
MonitorElementMET_HF
 
MonitorElementmeTotEnergy_HB
 
MonitorElementmeTotEnergy_HE
 
MonitorElementmeTotEnergy_HF
 
int nevent
 
MonitorElementNtowers_vs_ieta
 
MonitorElementnumBadCellsEcal_EB
 
MonitorElementnumBadCellsEcal_EE
 
MonitorElementnumBadCellsHcal_HB
 
MonitorElementnumBadCellsHcal_HE
 
MonitorElementnumBadCellsHcal_HF
 
MonitorElementnumPrbCellsEcal_EB
 
MonitorElementnumPrbCellsEcal_EE
 
MonitorElementnumPrbCellsHcal_HB
 
MonitorElementnumPrbCellsHcal_HE
 
MonitorElementnumPrbCellsHcal_HF
 
MonitorElementnumRcvCellsEcal_EB
 
MonitorElementnumRcvCellsEcal_EE
 
MonitorElementnumRcvCellsHcal_HB
 
MonitorElementnumRcvCellsHcal_HE
 
MonitorElementnumRcvCellsHcal_HF
 
MonitorElementoccupancy_map
 
MonitorElementoccupancy_vs_ieta
 
std::string outputFile_
 
MonitorElementphiMET_HB
 
MonitorElementphiMET_HE
 
MonitorElementphiMET_HF
 
MonitorElementSET_HB
 
MonitorElementSET_HE
 
MonitorElementSET_HF
 
edm::EDGetTokenT
< CaloTowerCollection
tok_towers_
 
bool useAllHistos_
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...> CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T...> HasAbility
 
typedef
CacheTypes::LuminosityBlockCache 
LuminosityBlockCache
 
typedef
LuminosityBlockContextT
< LuminosityBlockCache,
RunCache, GlobalCache
LuminosityBlockContext
 
typedef
CacheTypes::LuminosityBlockSummaryCache 
LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache,
GlobalCache
RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 
- Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase ModuleType
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr
< dqmDetails::NoCache
globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
 
static std::shared_ptr
< dqmDetails::NoCache
globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
 
static void globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
 
static void globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
 
- Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::stream::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 32 of file CaloTowersAnalyzer.h.

Member Typedef Documentation

Definition at line 52 of file CaloTowersAnalyzer.h.

Constructor & Destructor Documentation

CaloTowersAnalyzer::CaloTowersAnalyzer ( edm::ParameterSet const &  conf)

Definition at line 6 of file CaloTowersAnalyzer.cc.

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

6  {
7 
8  tok_towers_ = consumes<CaloTowerCollection>(conf.getUntrackedParameter<edm::InputTag>("CaloTowerCollectionLabel"));
9  // DQM ROOT output
10 
11  outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");
12 
13  hcalselector_ = conf.getUntrackedParameter<std::string>("hcalselector", "all");
14 
15  useAllHistos_ = conf.getUntrackedParameter<bool>("useAllHistos", false);
16 }
edm::EDGetTokenT< CaloTowerCollection > tok_towers_
tuple conf
Definition: dbtoconf.py:185
CaloTowersAnalyzer::~CaloTowersAnalyzer ( )

Definition at line 397 of file CaloTowersAnalyzer.cc.

397  {
398 
399 }

Member Function Documentation

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

Implements edm::stream::EDAnalyzerBase.

Definition at line 406 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, MonitorElement::Fill(), HcalObjRepresent::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, meNumFiredTowers_HB, meNumFiredTowers_HE, meNumFiredTowers_HF, caloMETBenchmarkGeneric_cfi::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_, and useAllHistos_.

406  {
407 
408  nevent++;
409 
411  event.getByToken(tok_towers_,towers);
413 
414  double met;
415  double phimet;
416 
417  // ieta scan
418  //double partR = 0.3;
419  //double Rmin = 9999.;
420  //double Econe = 0.;
421  //double Hcone = 0.;
422  //double Ee1 = 0.;
423  //double Eh1 = 0.;
424  double ieta_MC = 9999;
425  double iphi_MC = 9999;
426 
427  // HB
428  double sumEnergyHcal_HB = 0.;
429  double sumEnergyEcal_HB = 0.;
430  double sumEnergyHO_HB = 0.;
431  Int_t numFiredTowers_HB = 0;
432  double metx_HB = 0.;
433  double mety_HB = 0.;
434  double metz_HB = 0.;
435  double sEt_HB = 0.;
436  // HE
437  double sumEnergyHcal_HE = 0.;
438  double sumEnergyEcal_HE = 0.;
439  double sumEnergyHO_HE = 0.;
440  Int_t numFiredTowers_HE = 0;
441  double metx_HE = 0.;
442  double mety_HE = 0.;
443  double metz_HE = 0.;
444  double sEt_HE = 0.;
445  // HF
446  double sumEnergyHcal_HF = 0.;
447  double sumEnergyEcal_HF = 0.;
448  double sumEnergyHO_HF = 0.;
449  Int_t numFiredTowers_HF = 0;
450  double metx_HF = 0.;
451  double mety_HF = 0.;
452  double metz_HF = 0.;
453  double sEt_HF = 0.;
454 
455  for ( cal = towers->begin(); cal != towers->end(); ++cal ) {
456 
457  double eE = cal->emEnergy();
458  double eH = cal->hadEnergy();
459  double eHO = cal->outerEnergy();
460  double etaT = cal->eta();
461  // double phiT = cal->phi();
462  double en = cal->energy();
463  double etT = cal->et();
464  double had_tm = cal->hcalTime();
465  double em_tm = cal->ecalTime();
466 
467  int numBadEcalCells = cal->numBadEcalCells();
468  int numRcvEcalCells = cal->numRecoveredEcalCells();
469  int numPrbEcalCells = cal->numProblematicEcalCells();
470 
471  int numBadHcalCells = cal->numBadHcalCells();
472  int numRcvHcalCells = cal->numRecoveredHcalCells();
473  int numPrbHcalCells = cal->numProblematicHcalCells();
474 
475  math::RhoEtaPhiVector mom(cal->et(), cal->eta(), cal->phi());
476  // Vector mom = cal->momentum();
477 
478  // cell properties
479  CaloTowerDetId idT = cal->id();
480  int ieta = idT.ieta();
481  if(ieta > 0) ieta -= 1;
482  int iphi = idT.iphi();
483 
484  // ecal: 0 EcalBarrel 1 EcalEndcap
485  // hcal: 0 hcalBarrel 1 HcalEndcap 2 HcalForward
486  std::vector<int> inEcals(2), inHcals(3);
487  unsigned int constitSize = cal->constituentsSize();
488  for(unsigned int ic=0; ic<constitSize; ic++){
489  DetId detId = cal->constituent(ic);
490  if(detId.det() == DetId::Ecal){
491  if( detId.subdetId()==EcalBarrel ) inEcals[0] =1;
492  else if( detId.subdetId()==EcalEndcap ) inEcals[1] =1;
493  }
494  if(detId.det() == DetId::Hcal){
495  if( HcalDetId(detId).subdet()==HcalBarrel ) inHcals[0] =1;
496  else if( HcalDetId(detId).subdet()==HcalEndcap ) inHcals[1] =1;
497  else if( HcalDetId(detId).subdet()==HcalForward ) inHcals[2] =1;
498  }
499  }
500  //All cell histos are used
501  if( inEcals[0] ) {
502  numBadCellsEcal_EB->Fill(numBadEcalCells);
503  numRcvCellsEcal_EB->Fill(numRcvEcalCells);
504  numPrbCellsEcal_EB->Fill(numPrbEcalCells);
505  }
506  if( inEcals[1] ) {
507  numBadCellsEcal_EE->Fill(numBadEcalCells);
508  numRcvCellsEcal_EE->Fill(numRcvEcalCells);
509  numPrbCellsEcal_EE->Fill(numPrbEcalCells);
510  }
511 
512  //Ntowers is used in EndJob, occupancy_map is used for occupancy vs ieta
513  Ntowers_vs_ieta -> Fill(double(ieta),1.);
514  occupancy_map -> Fill(double(ieta),double(iphi));
515 
516 
517  if((isub == 0 || isub == 1)
518  && (fabs(etaT) < etaMax[0] && fabs(etaT) >= etaMin[0] )) {
519 
520  //All cell histos are used
521  numBadCellsHcal_HB->Fill(numBadHcalCells);
522  numRcvCellsHcal_HB->Fill(numRcvHcalCells);
523  numPrbCellsHcal_HB->Fill(numPrbHcalCells);
524 
525  //Map energy histos are not used
526  if (useAllHistos_){
527  mapEnergy_HB -> Fill(double(ieta), double(iphi), en);
528  mapEnergyHcal_HB -> Fill(double(ieta), double(iphi), eH);
529  mapEnergyEcal_HB -> Fill(double(ieta), double(iphi), eE);
530  }
531  // std::cout << " e_ecal = " << eE << std::endl;
532 
533  // simple sums
534  sumEnergyHcal_HB += eH;
535  sumEnergyEcal_HB += eE;
536  sumEnergyHO_HB += eHO;
537 
538  numFiredTowers_HB++;
539 
540  //Not used
541  if (useAllHistos_){
544  }
545 
546  // MET, SET & phimet
547  // double etT = cal->et();
548  metx_HB += mom.x();
549  mety_HB += mom.y(); //etT * sin(phiT);
550  sEt_HB += etT;
551 
552  //Timing (all histos are used)
553  emTiming_HB->Fill(em_tm);
554  hadTiming_HB->Fill(had_tm);
555 
556  emEnergyTiming_Low_HB->Fill(eE, em_tm);
557  emEnergyTiming_HB->Fill(eE, em_tm);
558  emEnergyTiming_High_HB->Fill(eE, em_tm);
560  emEnergyTiming_profile_HB->Fill(eE, em_tm);
562 
563  hadEnergyTiming_Low_HB->Fill(eH, had_tm);
564  hadEnergyTiming_HB->Fill(eH, had_tm);
565  hadEnergyTiming_High_HB->Fill(eH, had_tm);
566  hadEnergyTiming_profile_Low_HB->Fill(eH, had_tm);
567  hadEnergyTiming_profile_HB->Fill(eH, had_tm);
569  }
570 
571  if((isub == 0 || isub == 2)
572  && (fabs(etaT) < etaMax[1] && fabs(etaT) >= etaMin[1] )) {
573 
574  //All cell histos are used
575  numBadCellsHcal_HE->Fill(numBadHcalCells);
576  numRcvCellsHcal_HE->Fill(numRcvHcalCells);
577  numPrbCellsHcal_HE->Fill(numPrbHcalCells);
578 
579  //Map energy histos are not used
580  if (useAllHistos_){
581  mapEnergy_HE -> Fill(double(ieta), double(iphi), en);
582  mapEnergyHcal_HE -> Fill(double(ieta), double(iphi), eH);
583  mapEnergyEcal_HE -> Fill(double(ieta), double(iphi), eE);
584  }
585  // std::cout << " e_ecal = " << eE << std::endl;
586 
587  // simple sums
588  sumEnergyHcal_HE += eH;
589  sumEnergyEcal_HE += eE;
590  sumEnergyHO_HE += eHO;
591 
592  numFiredTowers_HE++;
593 
594  //Not used
595  if (useAllHistos_){
598  }
599  // MET, SET & phimet
600  // double etT = cal->et();
601  metx_HE += mom.x();
602  mety_HE += mom.y(); //etT * sin(phiT);
603  sEt_HE += etT;
604 
605  //Timing (all histos are used)
606  emTiming_HE->Fill(em_tm);
607  hadTiming_HE->Fill(had_tm);
608 
609  emEnergyTiming_Low_HE->Fill(eE, em_tm);
610  emEnergyTiming_HE->Fill(eE, em_tm);
612  emEnergyTiming_profile_HE->Fill(eE, em_tm);
613 
614  hadEnergyTiming_Low_HE->Fill(eH, had_tm);
615  hadEnergyTiming_HE->Fill(eH, had_tm);
616  hadEnergyTiming_profile_Low_HE->Fill(eH, had_tm);
617  hadEnergyTiming_profile_HE->Fill(eH, had_tm);
618  }
619 
620  if((isub == 0 || isub == 3)
621  && (fabs(etaT) < etaMax[2] && fabs(etaT) >= etaMin[2] )) {
622 
623  //All cell histos are used
624  numBadCellsHcal_HF->Fill(numBadHcalCells);
625  numRcvCellsHcal_HF->Fill(numRcvHcalCells);
626  numPrbCellsHcal_HF->Fill(numPrbHcalCells);
627 
628  //Map energy histos are not used
629  if (useAllHistos_){
630  mapEnergy_HF -> Fill(double(ieta), double(iphi), en);
631  mapEnergyHcal_HF -> Fill(double(ieta), double(iphi), eH);
632  mapEnergyEcal_HF -> Fill(double(ieta), double(iphi), eE);
633  }
634  // std::cout << " e_ecal = " << eE << std::endl;
635 
636  // simple sums
637  sumEnergyHcal_HF += eH;
638  sumEnergyEcal_HF += eE;
639  sumEnergyHO_HF += eHO;
640 
641  numFiredTowers_HF++;
642 
643  //Not used
644  if (useAllHistos_){
647  }
648  // MET, SET & phimet
649  // double etT = cal->et();
650  metx_HF += mom.x();
651  mety_HF += mom.y(); //etT * sin(phiT);
652  sEt_HF += etT;
653 
654  //Timing (all histos are used)
655  emTiming_HF->Fill(em_tm);
656  hadTiming_HF->Fill(had_tm);
657  emEnergyTiming_HF->Fill(eE, em_tm);
658  emEnergyTiming_profile_HF->Fill(eE, em_tm);
659 
660  hadEnergyTiming_Low_HF->Fill(eH, had_tm);
661  hadEnergyTiming_HF->Fill(eH, had_tm);
662  hadEnergyTiming_profile_Low_HF->Fill(eH, had_tm);
663  hadEnergyTiming_profile_HF->Fill(eH, had_tm);
664  }
665 
666 
667 
668  } // end of Towers cycle
669 
670  //These are the six single pion histos; only the second set is used
671 
672  mapEnergy_N -> Fill(double(ieta_MC), double(iphi_MC), 1.);
673 
674  if(isub == 0 || isub == 1) {
675  met = sqrt(metx_HB*metx_HB + mety_HB*mety_HB);
676  Vector metv(metx_HB,mety_HB,metz_HB);
677  phimet = metv.phi();
678 
679  //Five oldest drawn histos first; the rest are not used
680  meEnergyHcal_HB-> Fill(sumEnergyHcal_HB);
681  meEnergyEcal_HB-> Fill(sumEnergyEcal_HB);
682  meNumFiredTowers_HB-> Fill(numFiredTowers_HB);
683  MET_HB -> Fill (met);
684  SET_HB -> Fill (sEt_HB);
685 
686  if (useAllHistos_){
687  meEnergyHcalvsEcal_HB->Fill(sumEnergyEcal_HB, sumEnergyHcal_HB);
688  meEnergyHO_HB-> Fill(sumEnergyHO_HB);
689  meTotEnergy_HB-> Fill(sumEnergyHcal_HB+sumEnergyEcal_HB
690  +sumEnergyHO_HB);
691  phiMET_HB -> Fill (phimet);
692  }
693  }
694 
695 
696  if(isub == 0 || isub == 2) {
697  met = sqrt(metx_HE*metx_HE + mety_HE*mety_HE);
698  Vector metv(metx_HE,mety_HE,metz_HE);
699  phimet = metv.phi();
700 
701  //Five oldest drawn histos first; the rest are not used
702  meEnergyHcal_HE-> Fill(sumEnergyHcal_HE);
703  meEnergyEcal_HE-> Fill(sumEnergyEcal_HE);
704  meNumFiredTowers_HE-> Fill(numFiredTowers_HE);
705  MET_HE -> Fill (met);
706  SET_HE -> Fill (sEt_HE);
707 
708  if (useAllHistos_){
709  meEnergyHcalvsEcal_HE->Fill(sumEnergyEcal_HE, sumEnergyHcal_HE);
710  meEnergyHO_HE-> Fill(sumEnergyHO_HE);
711  meTotEnergy_HE-> Fill(sumEnergyHcal_HE+sumEnergyEcal_HE
712  +sumEnergyHO_HE);
713  phiMET_HE -> Fill (phimet);
714  }
715  }
716 
717  if(isub == 0 || isub == 3) {
718  met = sqrt(metx_HF*metx_HF + mety_HF*mety_HF);
719  Vector metv(metx_HF,mety_HF,metz_HF);
720  phimet = metv.phi();
721 
722  //Five oldest drawn histos first; the rest are not used
723  meEnergyHcal_HF-> Fill(sumEnergyHcal_HF);
724  meEnergyEcal_HF-> Fill(sumEnergyEcal_HF);
725  meNumFiredTowers_HF-> Fill(numFiredTowers_HF);
726  MET_HF -> Fill (met);
727  SET_HF -> Fill (sEt_HF);
728 
729  if (useAllHistos_){
730  meEnergyHcalvsEcal_HF->Fill(sumEnergyEcal_HF, sumEnergyHcal_HF);
731  meEnergyHO_HF-> Fill(sumEnergyHO_HF);
732  meTotEnergy_HF-> Fill(sumEnergyHcal_HF+sumEnergyEcal_HF
733  +sumEnergyHO_HF);
734  phiMET_HF -> Fill (phimet);
735  }
736  }
737 
738 }
MonitorElement * mapEnergyEcal_HF
MonitorElement * meNumFiredTowers_HF
MonitorElement * hadEnergyTiming_profile_High_HB
MonitorElement * hadEnergyTiming_HF
MonitorElement * meEnergyHcalvsEcal_HE
MonitorElement * meEnergyHcalTower_HF
HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:30
MonitorElement * hadEnergyTiming_Low_HE
MonitorElement * meEnergyHcalTower_HE
MonitorElement * numRcvCellsHcal_HF
MonitorElement * MET_HE
MonitorElement * phiMET_HF
MonitorElement * emEnergyTiming_profile_Low_HB
MonitorElement * phiMET_HE
ROOT::Math::Plane3D::Vector Vector
Definition: EcalHitMaker.cc:29
MonitorElement * meEnergyEcalTower_HF
MonitorElement * hadEnergyTiming_profile_Low_HE
MonitorElement * emEnergyTiming_Low_HB
MonitorElement * emTiming_HE
MonitorElement * meEnergyEcalTower_HE
MonitorElement * numRcvCellsEcal_EB
std::vector< CaloTower >::const_iterator const_iterator
MonitorElement * meEnergyHcalTower_HB
MonitorElement * emEnergyTiming_profile_Low_HE
MonitorElement * numRcvCellsEcal_EE
MonitorElement * SET_HF
MonitorElement * meEnergyHcal_HF
MonitorElement * meEnergyHcal_HB
MonitorElement * meEnergyHO_HB
MonitorElement * emTiming_HB
MonitorElement * MET_HF
MonitorElement * occupancy_map
MonitorElement * emEnergyTiming_HE
MonitorElement * hadEnergyTiming_Low_HB
MonitorElement * hadEnergyTiming_profile_Low_HB
RhoEtaPhiVectorD RhoEtaPhiVector
spatial vector with cylindrical internal representation using pseudorapidity
Definition: Vector3D.h:32
MonitorElement * hadTiming_HE
MonitorElement * SET_HB
void Fill(long long x)
MonitorElement * mapEnergyEcal_HB
MonitorElement * meEnergyHO_HF
MonitorElement * phiMET_HB
MonitorElement * meEnergyEcal_HE
MonitorElement * hadEnergyTiming_profile_HB
MonitorElement * meEnergyEcal_HB
MonitorElement * numBadCellsHcal_HE
MonitorElement * meTotEnergy_HE
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
MonitorElement * hadEnergyTiming_Low_HF
MonitorElement * SET_HE
T sqrt(T t)
Definition: SSEVec.h:48
MonitorElement * numRcvCellsHcal_HB
MonitorElement * meNumFiredTowers_HB
MonitorElement * MET_HB
MonitorElement * numPrbCellsHcal_HB
MonitorElement * mapEnergyHcal_HF
MonitorElement * meEnergyHcal_HE
MonitorElement * hadTiming_HB
edm::EDGetTokenT< CaloTowerCollection > tok_towers_
MonitorElement * numBadCellsEcal_EE
MonitorElement * numPrbCellsEcal_EE
MonitorElement * mapEnergy_HE
MonitorElement * mapEnergy_HF
MonitorElement * mapEnergy_HB
MonitorElement * meEnergyEcalTower_HB
MonitorElement * emEnergyTiming_High_HB
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:37
MonitorElement * numRcvCellsHcal_HE
MonitorElement * emEnergyTiming_profile_HF
MonitorElement * hadEnergyTiming_HE
MonitorElement * numPrbCellsHcal_HE
MonitorElement * emEnergyTiming_profile_High_HB
MonitorElement * meTotEnergy_HB
MonitorElement * meTotEnergy_HF
MonitorElement * hadEnergyTiming_High_HB
MonitorElement * mapEnergyHcal_HE
MonitorElement * hadEnergyTiming_HB
MonitorElement * numPrbCellsHcal_HF
Definition: DetId.h:18
MonitorElement * hadEnergyTiming_profile_HF
MonitorElement * numBadCellsHcal_HB
MonitorElement * Ntowers_vs_ieta
MonitorElement * emEnergyTiming_HB
MonitorElement * meEnergyHO_HE
MonitorElement * hadEnergyTiming_profile_HE
MonitorElement * emEnergyTiming_Low_HE
MonitorElement * numPrbCellsEcal_EB
MonitorElement * numBadCellsHcal_HF
MonitorElement * hadTiming_HF
MonitorElement * numBadCellsEcal_EB
MonitorElement * mapEnergyEcal_HE
MonitorElement * emTiming_HF
MonitorElement * emEnergyTiming_profile_HE
MonitorElement * mapEnergy_N
int ieta() const
get the tower ieta
MonitorElement * meEnergyHcalvsEcal_HB
Detector det() const
get the detector field from this detid
Definition: DetId.h:35
MonitorElement * mapEnergyHcal_HB
MonitorElement * emEnergyTiming_profile_HB
MonitorElement * meEnergyHcalvsEcal_HF
MonitorElement * meNumFiredTowers_HE
MonitorElement * emEnergyTiming_HF
MonitorElement * hadEnergyTiming_profile_Low_HF
MonitorElement * meEnergyEcal_HF
void CaloTowersAnalyzer::beginJob ( void  )
virtual

Definition at line 404 of file CaloTowersAnalyzer.cc.

References nevent.

404 { nevent = 0; }
void CaloTowersAnalyzer::beginRun ( void  )
virtual

Definition at line 393 of file CaloTowersAnalyzer.cc.

393 {}
void CaloTowersAnalyzer::bookHistograms ( DQMStore::IBooker ibooker,
edm::Run const &  ,
edm::EventSetup const &   
)
overridevirtual

Implements DQMEDAnalyzer.

Definition at line 18 of file CaloTowersAnalyzer.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::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, 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, 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, DQMStore::IBooker::setCurrentFolder(), and useAllHistos_.

19 {
20 
21  etaMin[0] = 0.;
22  etaMax[0] = 1.4;
23  etaMin[1] = 1.4;
24  etaMax[1] = 2.9;
25  etaMin[2] = 2.9;
26  etaMax[2] = 5.2;
27 
28  isub = 0;
29  if(hcalselector_ == "HB") isub = 1;
30  if(hcalselector_ == "HE") isub = 2;
31  if(hcalselector_ == "HF") isub = 3;
32 
33  if ( outputFile_.size() != 0 ) {
34  edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will be saved to '" << outputFile_.c_str() << "'";
35  } else {
36  edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will NOT be saved";
37  }
38 
39  Char_t histo[100];
40 
41 
42  ibooker.setCurrentFolder("CaloTowersD/CaloTowersTask") ;
43 
44 
45 
46  //These two histos are not drawn by our macros, but they are used
47  //in the EndJob for norms and such so I am leaving them alone for now
48  //-------------------------------------------------------------------------------------------
49  sprintf (histo, "Ntowers_per_event_vs_ieta" );
50  Ntowers_vs_ieta = ibooker.book1D(histo, histo, 82, -41., 41.);
51 
52  sprintf (histo, "CaloTowersTask_map_Nentries" );
53  mapEnergy_N = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
54  //-------------------------------------------------------------------------------------------
55 
56  //These the single pion scan histos
57  //-------------------------------------------------------------------------------------------
58  //The first three are not used
59  if (useAllHistos_){
60  sprintf (histo, "emean_vs_ieta_E" );
61  emean_vs_ieta_E = ibooker.bookProfile(histo, histo, 82, -41., 41., 2100, -100., 2000.);
62  sprintf (histo, "emean_vs_ieta_H" );
63  emean_vs_ieta_H = ibooker.bookProfile(histo, histo, 82, -41., 41., 2100, -100., 2000.);
64  sprintf (histo, "emean_vs_ieta_EH" );
65  emean_vs_ieta_EH = ibooker.bookProfile(histo, histo, 82, -41., 41., 2100, -100., 2000.);
66  }
67  //These are drawn
68  sprintf (histo, "emean_vs_ieta_E1" );
69  emean_vs_ieta_E1 = ibooker.bookProfile(histo, histo, 82, -41., 41., 2100, -100., 2000.);
70  sprintf (histo, "emean_vs_ieta_H1" );
71  emean_vs_ieta_H1 = ibooker.bookProfile(histo, histo, 82, -41., 41., 2100, -100., 2000.);
72  sprintf (histo, "emean_vs_ieta_EH1" );
73  emean_vs_ieta_EH1 = ibooker.bookProfile(histo, histo, 82, -41., 41., 2100, -100., 2000.);
74  //-------------------------------------------------------------------------------------------
75 
76  //Map energy histos are not drawn
77  if (useAllHistos_){
78  sprintf (histo, "CaloTowersTask_map_energy_E" );
79  mapEnergy_E = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
80  sprintf (histo, "CaloTowersTask_map_energy_H");
81  mapEnergy_H = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
82  sprintf (histo, "CaloTowersTask_map_energy_EH" );
83  mapEnergy_EH = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
84  }
85 
86  //All ECAL cell histos are used
87  // XXX: ECAL 0-25 [0-26, 26 bins] HCAL 0-4 [0-5, 5 bins]
88  sprintf (histo, "number_of_bad_cells_Ecal_EB");
89  numBadCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
90  sprintf (histo, "number_of_bad_cells_Ecal_EE");
91  numBadCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26);
92  sprintf (histo, "number_of_recovered_cells_Ecal_EB");
93  numRcvCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
94  sprintf (histo, "number_of_recovered_cells_Ecal_EE");
95  numRcvCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26);
96  sprintf (histo, "number_of_problematic_cells_Ecal_EB");
97  numPrbCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
98  sprintf (histo, "number_of_problematic_cells_Ecal_EE");
99  numPrbCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26);
100 
101  //Occupancy vs. ieta is drawn, occupancy map is needed to draw it
102  sprintf (histo, "CaloTowersTask_map_occupancy" );
103  occupancy_map = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
104 
105  sprintf (histo, "CaloTowersTask_occupancy_vs_ieta" );
106  occupancy_vs_ieta = ibooker.book1D(histo, histo, 82, -41, 41);
107 
108  if( isub == 1 || isub == 0) {
109  //All cell histos are used
110  sprintf (histo, "number_of_bad_cells_Hcal_HB");
111  numBadCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5);
112  sprintf (histo, "number_of_recovered_cells_Hcal_HB");
113  numRcvCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5);
114  sprintf (histo, "number_of_problematic_cells_Hcal_HB");
115  numPrbCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5);
116 
117  //These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL, MET and SET
118  //-------------------------------------------------------------------------------------------
119  sprintf (histo, "CaloTowersTask_energy_HCAL_HB" ) ;
120  meEnergyHcal_HB = ibooker.book1D(histo, histo, 4100, -200, 8000);
121 
122  sprintf (histo, "CaloTowersTask_energy_ECAL_HB" ) ;
123  meEnergyEcal_HB = ibooker.book1D(histo, histo, 3100, -200, 6000);
124 
125  sprintf (histo, "CaloTowersTask_number_of_fired_towers_HB" ) ;
126  meNumFiredTowers_HB = ibooker.book1D(histo, histo, 1000, 0, 2000);
127 
128  sprintf (histo, "CaloTowersTask_MET_HB" ) ;
129  MET_HB = ibooker.book1D(histo, histo, 3000, 0. , 3000. ) ;
130 
131  sprintf (histo, "CaloTowersTask_SET_HB" ) ;
132  SET_HB = ibooker.book1D(histo, histo, 8000, 0. , 8000. ) ;
133  //-------------------------------------------------------------------------------------------
134 
135  //Timing histos and profiles -- all six are necessary
136  //-------------------------------------------------------------------------------------------
137  sprintf (histo, "CaloTowersTask_EM_Timing_HB" ) ;
138  emTiming_HB = ibooker.book1D(histo, histo, 110, -120., 100. ) ;
139 
140  sprintf (histo, "CaloTowersTask_HAD_Timing_HB" ) ;
141  hadTiming_HB = ibooker.book1D(histo, histo, 70, -48., 92. ) ;
142 
143  //Energy-Timing histos are divided into low, medium and high to reduce memory usage
144  //EM
145  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_Low_HB" ) ;
146  emEnergyTiming_Low_HB = ibooker.book2D(histo, histo, 40, 0. , 40., 110, -120., 100. ) ;
147 
148  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_HB" ) ;
149  emEnergyTiming_HB = ibooker.book2D(histo, histo, 200, 0. , 400., 110, -120., 100. ) ;
150 
151  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_High_HB" ) ;
152  emEnergyTiming_High_HB = ibooker.book2D(histo, histo, 200, 0. , 3000., 110, -120., 100. ) ;
153 
154  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_profile_Low_HB" ) ;
155  emEnergyTiming_profile_Low_HB = ibooker.bookProfile(histo, histo, 40, 0. , 40., 110, -120., 100. ) ;
156 
157  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_profile_HB" ) ;
158  emEnergyTiming_profile_HB = ibooker.bookProfile(histo, histo, 200, 0. , 400., 110, -120., 100. ) ;
159 
160  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_profile_High_HB" ) ;
161  emEnergyTiming_profile_High_HB = ibooker.bookProfile(histo, histo, 200, 0. , 3000., 110, -120., 100. ) ;
162 
163  //HAD
164  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_Low_HB" ) ;
165  hadEnergyTiming_Low_HB = ibooker.book2D(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;
166 
167  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_HB" ) ;
168  hadEnergyTiming_HB = ibooker.book2D(histo, histo, 100, 0. , 200., 70, -48., 92. ) ;
169 
170  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_High_HB" ) ;
171  hadEnergyTiming_High_HB = ibooker.book2D(histo, histo, 300, 0. , 3000., 70, -48., 92. ) ;
172 
173  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HB" ) ;
174  hadEnergyTiming_profile_Low_HB = ibooker.bookProfile(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;
175 
176  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_profile_HB" ) ;
177  hadEnergyTiming_profile_HB = ibooker.bookProfile(histo, histo, 100, 0. , 200., 70, -48., 92. ) ;
178 
179  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_profile_High_HB" ) ;
180  hadEnergyTiming_profile_High_HB = ibooker.bookProfile(histo, histo, 300, 0. , 3000., 70, -48., 92. ) ;
181  //-------------------------------------------------------------------------------------------
182 
183  //Everything else is not drawn
184  if (useAllHistos_){
185  sprintf (histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HB") ;
186  meEnergyHcalvsEcal_HB = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);
187 
188  sprintf (histo, "CaloTowersTask_energy_OUTER_HB" ) ;
189  meEnergyHO_HB = ibooker.book1D(histo, histo, 1640, -200, 8000);
190 
191  sprintf (histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HB" ) ;
192  meEnergyEcalTower_HB = ibooker.book1D(histo, histo, 440, -200, 2000);
193 
194  sprintf (histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HB" ) ;
195  meEnergyHcalTower_HB = ibooker.book1D(histo, histo, 440 , -200 , 2000);
196 
197  sprintf (histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HB" ) ;
198  meTotEnergy_HB = ibooker.book1D(histo, histo,400, 0., 2000.) ;
199 
200  sprintf (histo, "CaloTowersTask_map_energy_HB" );
201  mapEnergy_HB = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
202  sprintf (histo, "CaloTowersTask_map_energy_HCAL_HB");
203  mapEnergyHcal_HB = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
204  sprintf (histo, "CaloTowersTask_map_energy_ECAL_HB" );
205  mapEnergyEcal_HB = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
206 
207  sprintf (histo, "CaloTowersTask_phi_MET_HB" ) ;
208  phiMET_HB = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898 ) ;
209  }
210  }
211 
212  if( isub == 2 || isub == 0) {
213  //All cell histos are used
214  sprintf (histo, "number_of_bad_cells_Hcal_HE");
215  numBadCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5);
216  sprintf (histo, "number_of_recovered_cells_Hcal_HE");
217  numRcvCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5);
218  sprintf (histo, "number_of_problematic_cells_Hcal_HE");
219  numPrbCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5);
220 
221  //These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL, MET and SET
222  //-------------------------------------------------------------------------------------------
223  sprintf (histo, "CaloTowersTask_energy_HCAL_HE" ) ;
224  meEnergyHcal_HE = ibooker.book1D(histo, histo, 1240, -200, 6000);
225 
226  sprintf (histo, "CaloTowersTask_energy_ECAL_HE" ) ;
227  meEnergyEcal_HE = ibooker.book1D(histo, histo, 1240, -200, 6000);
228 
229  sprintf (histo, "CaloTowersTask_number_of_fired_towers_HE" ) ;
230  meNumFiredTowers_HE = ibooker.book1D(histo, histo, 1000, 0, 2000);
231 
232  sprintf (histo, "CaloTowersTask_MET_HE" ) ;
233  MET_HE = ibooker.book1D(histo, histo, 1000, 0. , 1000. ) ;
234 
235  sprintf (histo, "CaloTowersTask_SET_HE" ) ;
236  SET_HE = ibooker.book1D(histo, histo, 2000, 0. , 2000. ) ;
237  //-------------------------------------------------------------------------------------------
238 
239  //Timing histos and profiles -- all six are necessary
240  //-------------------------------------------------------------------------------------------
241  sprintf (histo, "CaloTowersTask_EM_Timing_HE" ) ;
242  emTiming_HE = ibooker.book1D(histo, histo, 110, -120., 100. ) ;
243 
244  sprintf (histo, "CaloTowersTask_HAD_Timing_HE" ) ;
245  hadTiming_HE = ibooker.book1D(histo, histo, 70, -48., 92. ) ;
246 
247  //Energy-Timing histos are divided into low and normal to reduce memory usage
248  //EM
249  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_Low_HE" ) ;
250  emEnergyTiming_Low_HE = ibooker.book2D(histo, histo, 160, 0. , 160., 110, -120., 100. ) ;
251 
252  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_HE" ) ;
253  emEnergyTiming_HE = ibooker.book2D(histo, histo, 200, 0. , 800., 110, -120., 100. ) ;
254 
255  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_profile_Low_HE" ) ;
256  emEnergyTiming_profile_Low_HE = ibooker.bookProfile(histo, histo, 160, 0. , 160., 110, -120., 100. ) ;
257 
258  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_profile_HE" ) ;
259  emEnergyTiming_profile_HE = ibooker.bookProfile(histo, histo, 200, 0. , 800., 110, -120., 100. ) ;
260 
261  //HAD
262  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_Low_HE" ) ;
263  hadEnergyTiming_Low_HE = ibooker.book2D(histo, histo, 160, 0. , 160., 70, -48., 92. ) ;
264 
265  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_HE" ) ;
266  hadEnergyTiming_HE = ibooker.book2D(histo, histo, 200, 0. , 800., 70, -48., 92. ) ;
267 
268  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HE" ) ;
269  hadEnergyTiming_profile_Low_HE = ibooker.bookProfile(histo, histo, 160, 0. , 160., 70, -48., 92. ) ;
270 
271  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_profile_HE" ) ;
272  hadEnergyTiming_profile_HE = ibooker.bookProfile(histo, histo, 200, 0. , 800., 70, -48., 92. ) ;
273  //-------------------------------------------------------------------------------------------
274 
275  //Everything else is not drawn
276  if (useAllHistos_){
277  sprintf (histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HE") ;
278  meEnergyHcalvsEcal_HE = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);
279 
280  sprintf (histo, "CaloTowersTask_energy_OUTER_HE" ) ;
281  meEnergyHO_HE = ibooker.book1D(histo, histo, 440, -200, 2000);
282 
283  sprintf (histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HE" ) ;
284  meEnergyEcalTower_HE = ibooker.book1D(histo, histo, 1100, -200, 2000);
285 
286  sprintf (histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HE" ) ;
287  meEnergyHcalTower_HE = ibooker.book1D(histo, histo, 1100 , -200 , 2000);
288 
289  sprintf (histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HE" ) ;
290  meTotEnergy_HE = ibooker.book1D(histo, histo,400, 0., 2000.) ;
291 
292  sprintf (histo, "CaloTowersTask_map_energy_HE" );
293  mapEnergy_HE = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
294  sprintf (histo, "CaloTowersTask_map_energy_HCAL_HE");
295  mapEnergyHcal_HE = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
296  sprintf (histo, "CaloTowersTask_map_energy_ECAL_HE" );
297  mapEnergyEcal_HE = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
298 
299  sprintf (histo, "CaloTowersTask_phi_MET_HE" ) ;
300  phiMET_HE = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898 ) ;
301  }
302  }
303 
304 
305  if( isub == 3 || isub == 0) {
306  //All cell histos are used
307  sprintf (histo, "number_of_bad_cells_Hcal_HF");
308  numBadCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5);
309  sprintf (histo, "number_of_recovered_cells_Hcal_HF");
310  numRcvCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5);
311  sprintf (histo, "number_of_problematic_cells_Hcal_HF");
312  numPrbCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5);
313 
314  //These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL, MET and SET
315  //-------------------------------------------------------------------------------------------
316  sprintf (histo, "CaloTowersTask_energy_HCAL_HF" ) ;
317  meEnergyHcal_HF = ibooker.book1D(histo, histo, 4040, -200, 20000);
318 
319  sprintf (histo, "CaloTowersTask_energy_ECAL_HF" ) ;
320  meEnergyEcal_HF = ibooker.book1D(histo, histo, 2440, -200, 12000);
321 
322  sprintf (histo, "CaloTowersTask_number_of_fired_towers_HF" ) ;
323  meNumFiredTowers_HF = ibooker.book1D(histo, histo, 1000, 0, 2000);
324 
325  sprintf (histo, "CaloTowersTask_MET_HF" ) ;
326  MET_HF = ibooker.book1D(histo, histo, 500, 0. , 500. ) ;
327 
328  sprintf (histo, "CaloTowersTask_SET_HF" ) ;
329  SET_HF = ibooker.book1D(histo, histo, 2000, 0. , 2000. ) ;
330  //-------------------------------------------------------------------------------------------
331 
332  //Timing histos and profiles -- all six are necessary
333  //-------------------------------------------------------------------------------------------
334  sprintf (histo, "CaloTowersTask_EM_Timing_HF" ) ;
335  emTiming_HF = ibooker.book1D(histo, histo, 110, -120., 100. ) ;
336 
337  sprintf (histo, "CaloTowersTask_HAD_Timing_HF" ) ;
338  hadTiming_HF = ibooker.book1D(histo, histo, 70, -48., 92. ) ;
339 
340  //EM
341  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_HF" ) ;
342  emEnergyTiming_HF = ibooker.book2D(histo, histo, 150, 0. , 300., 110, -120., 100. ) ;
343 
344  sprintf (histo, "CaloTowersTask_EM_Energy_Timing_profile_HF" ) ;
345  emEnergyTiming_profile_HF = ibooker.bookProfile(histo, histo, 150, 0. , 300., 110, -120., 100. ) ;
346 
347  //HAD (requires two different sets of histograms to lower RAM usage)
348  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_Low_HF" ) ;
349  hadEnergyTiming_Low_HF = ibooker.book2D(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;
350 
351  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_HF" ) ;
352  hadEnergyTiming_HF = ibooker.book2D(histo, histo, 200, 0. , 600., 70, -48., 92. ) ;
353 
354  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HF" ) ;
355  hadEnergyTiming_profile_Low_HF = ibooker.bookProfile(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;
356 
357  sprintf (histo, "CaloTowersTask_HAD_Energy_Timing_profile_HF" ) ;
358  hadEnergyTiming_profile_HF = ibooker.bookProfile(histo, histo, 200, 0. , 600., 70, -48., 92. ) ;
359  //-------------------------------------------------------------------------------------------
360 
361  //Everything else is not drawn
362  if (useAllHistos_){
363  sprintf (histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HF") ;
364  meEnergyHcalvsEcal_HF = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);
365 
366  sprintf (histo, "CaloTowersTask_energy_OUTER_HF" ) ;
367  meEnergyHO_HF = ibooker.book1D(histo, histo, 440, -200, 2000);
368 
369  sprintf (histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HF" ) ;
370  meEnergyEcalTower_HF = ibooker.book1D(histo, histo, 440, -200, 2000);
371 
372  sprintf (histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HF" ) ;
373  meEnergyHcalTower_HF = ibooker.book1D(histo, histo, 440 , -200 , 2000);
374 
375  sprintf (histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HF" ) ;
376  meTotEnergy_HF = ibooker.book1D(histo, histo, 400, 0., 2000.) ;
377 
378  sprintf (histo, "CaloTowersTask_map_energy_HF" );
379  mapEnergy_HF = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
380  sprintf (histo, "CaloTowersTask_map_energy_HCAL_HF");
381  mapEnergyHcal_HF = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
382  sprintf (histo, "CaloTowersTask_map_energy_ECAL_HF" );
383  mapEnergyEcal_HF = ibooker.book2D(histo, histo, 82, -41., 41., 72, 0., 72.);
384 
385  sprintf (histo, "CaloTowersTask_phi_MET_HF" ) ;
386  phiMET_HF = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898 ) ;
387  }
388  }
389 
390 }
MonitorElement * emean_vs_ieta_EH1
MonitorElement * mapEnergyEcal_HF
MonitorElement * meNumFiredTowers_HF
MonitorElement * hadEnergyTiming_profile_High_HB
MonitorElement * hadEnergyTiming_HF
MonitorElement * meEnergyHcalTower_HF
MonitorElement * meEnergyHcalvsEcal_HE
MonitorElement * hadEnergyTiming_Low_HE
MonitorElement * meEnergyHcalTower_HE
MonitorElement * numRcvCellsHcal_HF
MonitorElement * MET_HE
MonitorElement * phiMET_HF
MonitorElement * emEnergyTiming_profile_Low_HB
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
MonitorElement * phiMET_HE
MonitorElement * meEnergyEcalTower_HF
MonitorElement * hadEnergyTiming_profile_Low_HE
MonitorElement * emEnergyTiming_Low_HB
MonitorElement * emTiming_HE
MonitorElement * meEnergyEcalTower_HE
MonitorElement * numRcvCellsEcal_EB
MonitorElement * meEnergyHcalTower_HB
MonitorElement * emEnergyTiming_profile_Low_HE
MonitorElement * numRcvCellsEcal_EE
MonitorElement * SET_HF
MonitorElement * meEnergyHcal_HF
MonitorElement * meEnergyHcal_HB
MonitorElement * meEnergyHO_HB
MonitorElement * emTiming_HB
MonitorElement * MET_HF
MonitorElement * occupancy_map
MonitorElement * emEnergyTiming_HE
MonitorElement * hadEnergyTiming_Low_HB
MonitorElement * hadEnergyTiming_profile_Low_HB
MonitorElement * emean_vs_ieta_H
MonitorElement * hadTiming_HE
MonitorElement * SET_HB
MonitorElement * mapEnergyEcal_HB
MonitorElement * meEnergyHO_HF
MonitorElement * occupancy_vs_ieta
MonitorElement * mapEnergy_EH
MonitorElement * phiMET_HB
MonitorElement * meEnergyEcal_HE
MonitorElement * hadEnergyTiming_profile_HB
MonitorElement * meEnergyEcal_HB
MonitorElement * numBadCellsHcal_HE
MonitorElement * meTotEnergy_HE
MonitorElement * mapEnergy_E
MonitorElement * hadEnergyTiming_Low_HF
MonitorElement * SET_HE
MonitorElement * numRcvCellsHcal_HB
MonitorElement * meNumFiredTowers_HB
MonitorElement * MET_HB
MonitorElement * numPrbCellsHcal_HB
MonitorElement * mapEnergyHcal_HF
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
MonitorElement * meEnergyHcal_HE
MonitorElement * hadTiming_HB
MonitorElement * numBadCellsEcal_EE
MonitorElement * numPrbCellsEcal_EE
MonitorElement * mapEnergy_HE
MonitorElement * mapEnergy_HF
MonitorElement * mapEnergy_HB
MonitorElement * meEnergyEcalTower_HB
MonitorElement * emean_vs_ieta_H1
MonitorElement * emEnergyTiming_High_HB
MonitorElement * numRcvCellsHcal_HE
MonitorElement * emEnergyTiming_profile_HF
MonitorElement * hadEnergyTiming_HE
MonitorElement * numPrbCellsHcal_HE
MonitorElement * emEnergyTiming_profile_High_HB
MonitorElement * meTotEnergy_HB
MonitorElement * meTotEnergy_HF
MonitorElement * hadEnergyTiming_High_HB
MonitorElement * mapEnergyHcal_HE
MonitorElement * hadEnergyTiming_HB
MonitorElement * numPrbCellsHcal_HF
MonitorElement * hadEnergyTiming_profile_HF
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
MonitorElement * numBadCellsHcal_HB
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
MonitorElement * Ntowers_vs_ieta
MonitorElement * emEnergyTiming_HB
MonitorElement * meEnergyHO_HE
MonitorElement * hadEnergyTiming_profile_HE
MonitorElement * emEnergyTiming_Low_HE
MonitorElement * emean_vs_ieta_E
MonitorElement * numPrbCellsEcal_EB
MonitorElement * emean_vs_ieta_EH
MonitorElement * numBadCellsHcal_HF
MonitorElement * hadTiming_HF
MonitorElement * numBadCellsEcal_EB
MonitorElement * mapEnergyEcal_HE
MonitorElement * emTiming_HF
MonitorElement * emEnergyTiming_profile_HE
MonitorElement * mapEnergy_N
MonitorElement * meEnergyHcalvsEcal_HB
MonitorElement * mapEnergyHcal_HB
MonitorElement * emEnergyTiming_profile_HB
MonitorElement * meEnergyHcalvsEcal_HF
MonitorElement * emean_vs_ieta_E1
MonitorElement * meNumFiredTowers_HE
MonitorElement * emEnergyTiming_HF
MonitorElement * hadEnergyTiming_profile_Low_HF
MonitorElement * mapEnergy_H
MonitorElement * meEnergyEcal_HF
double CaloTowersAnalyzer::dR ( double  eta1,
double  phi1,
double  eta2,
double  phi2 
)
private

Definition at line 740 of file CaloTowersAnalyzer.cc.

References PI, mathSSE::sqrt(), and tmp.

740  {
741  double PI = 3.1415926535898;
742  double deltaphi= phi1 - phi2;
743  if( phi2 > phi1 ) { deltaphi= phi2 - phi1;}
744  if(deltaphi > PI) { deltaphi = 2.*PI - deltaphi;}
745  double deltaeta = eta2 - eta1;
746  double tmp = sqrt(deltaeta* deltaeta + deltaphi*deltaphi);
747  return tmp;
748 }
#define PI
T sqrt(T t)
Definition: SSEVec.h:48
std::vector< std::vector< double > > tmp
Definition: MVATrainer.cc:100
void CaloTowersAnalyzer::endJob ( void  )
virtual

Definition at line 401 of file CaloTowersAnalyzer.cc.

401  {
402 }
void CaloTowersAnalyzer::endRun ( void  )
virtual

Definition at line 395 of file CaloTowersAnalyzer.cc.

395 {}

Member Data Documentation

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_E
private

Definition at line 66 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_E1
private

Definition at line 70 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_EH
private

Definition at line 68 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_EH1
private

Definition at line 72 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_H
private

Definition at line 67 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_H1
private

Definition at line 71 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HB
private

Definition at line 124 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HE
private

Definition at line 159 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HF
private

Definition at line 188 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_High_HB
private

Definition at line 125 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_Low_HB
private

Definition at line 123 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_Low_HE
private

Definition at line 158 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HB
private

Definition at line 127 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HE
private

Definition at line 161 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HF
private

Definition at line 189 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_High_HB
private

Definition at line 128 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_Low_HB
private

Definition at line 126 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_Low_HE
private

Definition at line 160 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emTiming_HB
private

Definition at line 120 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emTiming_HE
private

Definition at line 155 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emTiming_HF
private

Definition at line 186 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

double CaloTowersAnalyzer::etaMax[3]
private

Definition at line 62 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

double CaloTowersAnalyzer::etaMin[3]
private

Definition at line 63 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HB
private

Definition at line 131 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HE
private

Definition at line 164 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HF
private

Definition at line 192 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_High_HB
private

Definition at line 132 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HB
private

Definition at line 130 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HE
private

Definition at line 163 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HF
private

Definition at line 191 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HB
private

Definition at line 134 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HE
private

Definition at line 166 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HF
private

Definition at line 194 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_High_HB
private

Definition at line 135 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HB
private

Definition at line 133 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HE
private

Definition at line 165 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HF
private

Definition at line 193 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadTiming_HB
private

Definition at line 121 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadTiming_HE
private

Definition at line 156 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadTiming_HF
private

Definition at line 187 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::string CaloTowersAnalyzer::hcalselector_
private

Definition at line 48 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and CaloTowersAnalyzer().

int CaloTowersAnalyzer::imc
private

Definition at line 59 of file CaloTowersAnalyzer.h.

int CaloTowersAnalyzer::isub
private

Definition at line 56 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_E
private

Definition at line 79 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_EH
private

Definition at line 81 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_H
private

Definition at line 80 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HB
private

Definition at line 113 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HE
private

Definition at line 148 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HF
private

Definition at line 179 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_N
private

Definition at line 82 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HB
private

Definition at line 114 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HE
private

Definition at line 149 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HF
private

Definition at line 180 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HB
private

Definition at line 115 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HE
private

Definition at line 150 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HF
private

Definition at line 181 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::string CaloTowersAnalyzer::mc_
private

Definition at line 49 of file CaloTowersAnalyzer.h.

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HB
private

Definition at line 105 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HE
private

Definition at line 140 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HF
private

Definition at line 171 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HB
private

Definition at line 109 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HE
private

Definition at line 144 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HF
private

Definition at line 175 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HB
private

Definition at line 106 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HE
private

Definition at line 141 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HF
private

Definition at line 172 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HB
private

Definition at line 110 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HE
private

Definition at line 145 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HF
private

Definition at line 176 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HB
private

Definition at line 103 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HE
private

Definition at line 138 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HF
private

Definition at line 169 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HB
private

Definition at line 104 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HE
private

Definition at line 139 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HF
private

Definition at line 170 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HB
private

Definition at line 107 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HE
private

Definition at line 142 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HF
private

Definition at line 173 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::MET_HB
private

Definition at line 116 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::MET_HE
private

Definition at line 151 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::MET_HF
private

Definition at line 182 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HB
private

Definition at line 111 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HE
private

Definition at line 146 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HF
private

Definition at line 177 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

int CaloTowersAnalyzer::nevent
private

Definition at line 57 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* CaloTowersAnalyzer::Ntowers_vs_ieta
private

Definition at line 74 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsEcal_EB
private

Definition at line 85 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsEcal_EE
private

Definition at line 86 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HB
private

Definition at line 92 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HE
private

Definition at line 93 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HF
private

Definition at line 94 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsEcal_EB
private

Definition at line 89 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsEcal_EE
private

Definition at line 90 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HB
private

Definition at line 98 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HE
private

Definition at line 99 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HF
private

Definition at line 100 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsEcal_EB
private

Definition at line 87 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsEcal_EE
private

Definition at line 88 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HB
private

Definition at line 95 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HE
private

Definition at line 96 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HF
private

Definition at line 97 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::occupancy_map
private

Definition at line 75 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::occupancy_vs_ieta
private

Definition at line 76 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

std::string CaloTowersAnalyzer::outputFile_
private

Definition at line 47 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::phiMET_HB
private

Definition at line 118 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::phiMET_HE
private

Definition at line 153 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::phiMET_HF
private

Definition at line 184 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::SET_HB
private

Definition at line 117 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::SET_HE
private

Definition at line 152 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::SET_HF
private

Definition at line 183 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<CaloTowerCollection> CaloTowersAnalyzer::tok_towers_
private

Definition at line 54 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

bool CaloTowersAnalyzer::useAllHistos_
private

Definition at line 50 of file CaloTowersAnalyzer.h.

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