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:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

virtual void analyze (edm::Event const &e, edm::EventSetup const &c)
 
virtual void beginJob ()
 
virtual void beginRun ()
 
 CaloTowersAnalyzer (edm::ParameterSet const &conf)
 
virtual void endJob ()
 
virtual void endRun ()
 
 ~CaloTowersAnalyzer ()
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- 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
 
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

DQMStoredbe_
 
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::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
- 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 31 of file CaloTowersAnalyzer.h.

Member Typedef Documentation

Definition at line 50 of file CaloTowersAnalyzer.h.

Constructor & Destructor Documentation

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

Definition at line 6 of file CaloTowersAnalyzer.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), dbe_, 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, edm::ParameterSet::getUntrackedParameter(), 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, cppFunctionSkipper::operator, outputFile_, phiMET_HB, phiMET_HE, phiMET_HF, SET_HB, SET_HE, SET_HF, DQMStore::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, tok_towers_, and useAllHistos_.

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

Definition at line 398 of file CaloTowersAnalyzer.cc.

398  {
399 
400 }

Member Function Documentation

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

Implements edm::EDAnalyzer.

Definition at line 408 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, 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, CaloMET_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_.

408  {
409 
410  nevent++;
411 
413  event.getByToken(tok_towers_,towers);
415 
416  double met;
417  double phimet;
418 
419  // ieta scan
420  //double partR = 0.3;
421  //double Rmin = 9999.;
422  //double Econe = 0.;
423  //double Hcone = 0.;
424  //double Ee1 = 0.;
425  //double Eh1 = 0.;
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  if(isub == 0 || isub == 1) {
673  met = sqrt(metx_HB*metx_HB + mety_HB*mety_HB);
674  Vector metv(metx_HB,mety_HB,metz_HB);
675  phimet = metv.phi();
676 
677  //Five oldest drawn histos first; the rest are not used
678  meEnergyHcal_HB-> Fill(sumEnergyHcal_HB);
679  meEnergyEcal_HB-> Fill(sumEnergyEcal_HB);
680  meNumFiredTowers_HB-> Fill(numFiredTowers_HB);
681  MET_HB -> Fill (met);
682  SET_HB -> Fill (sEt_HB);
683 
684  if (useAllHistos_){
685  meEnergyHcalvsEcal_HB->Fill(sumEnergyEcal_HB, sumEnergyHcal_HB);
686  meEnergyHO_HB-> Fill(sumEnergyHO_HB);
687  meTotEnergy_HB-> Fill(sumEnergyHcal_HB+sumEnergyEcal_HB
688  +sumEnergyHO_HB);
689  phiMET_HB -> Fill (phimet);
690  }
691  }
692 
693 
694  if(isub == 0 || isub == 2) {
695  met = sqrt(metx_HE*metx_HE + mety_HE*mety_HE);
696  Vector metv(metx_HE,mety_HE,metz_HE);
697  phimet = metv.phi();
698 
699  //Five oldest drawn histos first; the rest are not used
700  meEnergyHcal_HE-> Fill(sumEnergyHcal_HE);
701  meEnergyEcal_HE-> Fill(sumEnergyEcal_HE);
702  meNumFiredTowers_HE-> Fill(numFiredTowers_HE);
703  MET_HE -> Fill (met);
704  SET_HE -> Fill (sEt_HE);
705 
706  if (useAllHistos_){
707  meEnergyHcalvsEcal_HE->Fill(sumEnergyEcal_HE, sumEnergyHcal_HE);
708  meEnergyHO_HE-> Fill(sumEnergyHO_HE);
709  meTotEnergy_HE-> Fill(sumEnergyHcal_HE+sumEnergyEcal_HE
710  +sumEnergyHO_HE);
711  phiMET_HE -> Fill (phimet);
712  }
713  }
714 
715  if(isub == 0 || isub == 3) {
716  met = sqrt(metx_HF*metx_HF + mety_HF*mety_HF);
717  Vector metv(metx_HF,mety_HF,metz_HF);
718  phimet = metv.phi();
719 
720  //Five oldest drawn histos first; the rest are not used
721  meEnergyHcal_HF-> Fill(sumEnergyHcal_HF);
722  meEnergyEcal_HF-> Fill(sumEnergyEcal_HF);
723  meNumFiredTowers_HF-> Fill(numFiredTowers_HF);
724  MET_HF -> Fill (met);
725  SET_HF -> Fill (sEt_HF);
726 
727  if (useAllHistos_){
728  meEnergyHcalvsEcal_HF->Fill(sumEnergyEcal_HF, sumEnergyHcal_HF);
729  meEnergyHO_HF-> Fill(sumEnergyHO_HF);
730  meTotEnergy_HF-> Fill(sumEnergyHcal_HF+sumEnergyEcal_HF
731  +sumEnergyHO_HF);
732  phiMET_HF -> Fill (phimet);
733  }
734  }
735 
736 }
MonitorElement * mapEnergyEcal_HF
MonitorElement * meNumFiredTowers_HF
MonitorElement * hadEnergyTiming_profile_High_HB
MonitorElement * hadEnergyTiming_HF
tuple met
____________________________________________________________________________||
Definition: CaloMET_cfi.py:7
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
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

Reimplemented from edm::EDAnalyzer.

Definition at line 406 of file CaloTowersAnalyzer.cc.

References nevent.

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

Definition at line 394 of file CaloTowersAnalyzer.cc.

394 {}
double CaloTowersAnalyzer::dR ( double  eta1,
double  phi1,
double  eta2,
double  phi2 
)
private

Definition at line 738 of file CaloTowersAnalyzer.cc.

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

738  {
739  double PI = 3.1415926535898;
740  double deltaphi= phi1 - phi2;
741  if( phi2 > phi1 ) { deltaphi= phi2 - phi1;}
742  if(deltaphi > PI) { deltaphi = 2.*PI - deltaphi;}
743  double deltaeta = eta2 - eta1;
744  double tmp = sqrt(deltaeta* deltaeta + deltaphi*deltaphi);
745  return tmp;
746 }
#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

Reimplemented from edm::EDAnalyzer.

Definition at line 402 of file CaloTowersAnalyzer.cc.

References dbe_, outputFile_, and DQMStore::save().

402  {
403  if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
404 }
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", const uint32_t run=0, SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:2296
void CaloTowersAnalyzer::endRun ( void  )
virtual

Definition at line 396 of file CaloTowersAnalyzer.cc.

396 {}

Member Data Documentation

DQMStore* CaloTowersAnalyzer::dbe_
private

Definition at line 44 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer(), and endJob().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_E
private

Definition at line 64 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_E1
private

Definition at line 68 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_EH
private

Definition at line 66 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_EH1
private

Definition at line 70 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_H
private

Definition at line 65 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_H1
private

Definition at line 69 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HB
private

Definition at line 122 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HE
private

Definition at line 157 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HF
private

Definition at line 186 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_High_HB
private

Definition at line 123 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_Low_HB
private

Definition at line 121 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_Low_HE
private

Definition at line 156 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HB
private

Definition at line 125 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HE
private

Definition at line 159 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HF
private

Definition at line 187 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_High_HB
private

Definition at line 126 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_Low_HB
private

Definition at line 124 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_Low_HE
private

Definition at line 158 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emTiming_HB
private

Definition at line 118 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emTiming_HE
private

Definition at line 153 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::emTiming_HF
private

Definition at line 184 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

double CaloTowersAnalyzer::etaMax[3]
private

Definition at line 60 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

double CaloTowersAnalyzer::etaMin[3]
private

Definition at line 61 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HB
private

Definition at line 129 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HE
private

Definition at line 162 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HF
private

Definition at line 190 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_High_HB
private

Definition at line 130 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HB
private

Definition at line 128 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HE
private

Definition at line 161 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HF
private

Definition at line 189 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HB
private

Definition at line 132 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HE
private

Definition at line 164 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HF
private

Definition at line 192 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_High_HB
private

Definition at line 133 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HB
private

Definition at line 131 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HE
private

Definition at line 163 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HF
private

Definition at line 191 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadTiming_HB
private

Definition at line 119 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadTiming_HE
private

Definition at line 154 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::hadTiming_HF
private

Definition at line 185 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

std::string CaloTowersAnalyzer::hcalselector_
private

Definition at line 46 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

int CaloTowersAnalyzer::imc
private

Definition at line 57 of file CaloTowersAnalyzer.h.

int CaloTowersAnalyzer::isub
private

Definition at line 54 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergy_E
private

Definition at line 77 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergy_EH
private

Definition at line 79 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergy_H
private

Definition at line 78 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HB
private

Definition at line 111 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HE
private

Definition at line 146 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HF
private

Definition at line 177 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergy_N
private

Definition at line 80 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HB
private

Definition at line 112 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HE
private

Definition at line 147 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HF
private

Definition at line 178 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HB
private

Definition at line 113 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HE
private

Definition at line 148 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HF
private

Definition at line 179 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

std::string CaloTowersAnalyzer::mc_
private

Definition at line 47 of file CaloTowersAnalyzer.h.

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HB
private

Definition at line 103 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HE
private

Definition at line 138 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HF
private

Definition at line 169 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HB
private

Definition at line 107 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HE
private

Definition at line 142 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HF
private

Definition at line 173 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HB
private

Definition at line 104 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HE
private

Definition at line 139 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HF
private

Definition at line 170 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HB
private

Definition at line 108 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HE
private

Definition at line 143 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HF
private

Definition at line 174 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HB
private

Definition at line 101 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HE
private

Definition at line 136 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HF
private

Definition at line 167 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HB
private

Definition at line 102 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HE
private

Definition at line 137 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HF
private

Definition at line 168 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HB
private

Definition at line 105 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HE
private

Definition at line 140 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HF
private

Definition at line 171 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::MET_HB
private

Definition at line 114 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::MET_HE
private

Definition at line 149 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::MET_HF
private

Definition at line 180 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HB
private

Definition at line 109 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HE
private

Definition at line 144 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HF
private

Definition at line 175 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

int CaloTowersAnalyzer::nevent
private

Definition at line 55 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* CaloTowersAnalyzer::Ntowers_vs_ieta
private

Definition at line 72 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numBadCellsEcal_EB
private

Definition at line 83 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numBadCellsEcal_EE
private

Definition at line 84 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HB
private

Definition at line 90 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HE
private

Definition at line 91 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HF
private

Definition at line 92 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numPrbCellsEcal_EB
private

Definition at line 87 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numPrbCellsEcal_EE
private

Definition at line 88 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HB
private

Definition at line 96 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HE
private

Definition at line 97 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HF
private

Definition at line 98 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numRcvCellsEcal_EB
private

Definition at line 85 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numRcvCellsEcal_EE
private

Definition at line 86 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HB
private

Definition at line 93 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HE
private

Definition at line 94 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HF
private

Definition at line 95 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::occupancy_map
private

Definition at line 73 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::occupancy_vs_ieta
private

Definition at line 74 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer().

std::string CaloTowersAnalyzer::outputFile_
private

Definition at line 45 of file CaloTowersAnalyzer.h.

Referenced by CaloTowersAnalyzer(), and endJob().

MonitorElement* CaloTowersAnalyzer::phiMET_HB
private

Definition at line 116 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::phiMET_HE
private

Definition at line 151 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::phiMET_HF
private

Definition at line 182 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::SET_HB
private

Definition at line 115 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::SET_HE
private

Definition at line 150 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::SET_HF
private

Definition at line 181 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

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

Definition at line 52 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

bool CaloTowersAnalyzer::useAllHistos_
private

Definition at line 48 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().