CMS 3D CMS Logo
PFClusterValidation.cc
Go to the documentation of this file.
1 #include "DQMServices/Core/interface/DQMEDAnalyzer.h"
2 #include "DQMServices/Core/interface/DQMStore.h"
3 #include "DataFormats/CaloRecHit/interface/CaloCluster.h"
4 #include "DataFormats/CaloRecHit/interface/CaloClusterFwd.h"
5 #include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"
6 #include "DataFormats/CaloTowers/interface/CaloTowerDetId.h"
7 #include "DataFormats/Common/interface/Handle.h"
8 #include "DataFormats/DetId/interface/DetId.h"
9 #include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
10 #include "DataFormats/HcalDetId/interface/HcalDetId.h"
11 #include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
12 #include "DataFormats/Math/interface/Vector3D.h"
13 #include "DataFormats/Math/interface/deltaR.h"
14 #include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
15 #include "DataFormats/ParticleFlowReco/interface/PFBlock.h"
16 #include "DataFormats/ParticleFlowReco/interface/PFBlockElementCluster.h"
17 #include "DataFormats/ParticleFlowReco/interface/PFBlockElementTrack.h"
18 #include "DataFormats/ParticleFlowReco/interface/PFCluster.h"
19 #include "DataFormats/ParticleFlowReco/interface/PFClusterFwd.h"
20 #include "DataFormats/ParticleFlowReco/interface/PFLayer.h"
21 #include "DataFormats/ParticleFlowReco/interface/PFRecHit.h"
22 #include "DataFormats/ParticleFlowReco/interface/PFRecHitFraction.h"
23 #include "DataFormats/PatCandidates/interface/PackedCandidate.h"
24 #include "FWCore/Framework/interface/EDAnalyzer.h"
25 #include "FWCore/Framework/interface/ESHandle.h"
26 #include "FWCore/Framework/interface/Event.h"
27 #include "FWCore/Framework/interface/EventSetup.h"
28 #include "FWCore/Framework/interface/Frameworkfwd.h"
29 #include "FWCore/MessageLogger/interface/MessageLogger.h"
30 #include "FWCore/ParameterSet/interface/ParameterSet.h"
31 #include "FWCore/PluginManager/interface/ModuleDef.h"
32 #include "FWCore/ServiceRegistry/interface/Service.h"
33 #include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
34 
35 #include <algorithm>
36 #include <cmath>
37 #include <iostream>
38 #include <string>
39 #include <utility>
40 #include <vector>
41 
42 class PFClusterValidation : public DQMEDAnalyzer {
43 public:
44  PFClusterValidation(edm::ParameterSet const& conf);
45  ~PFClusterValidation() override;
46  void analyze(edm::Event const& e, edm::EventSetup const& c) override;
47  void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
48 
49 private:
50  static constexpr double partR2 = 0.3 * 0.3; // dr cutoff (squared)
51  static double sumEnergy(edm::Handle<reco::PFClusterCollection> const& pfClusters, double eta, double phi);
52 
53  edm::EDGetTokenT<edm::HepMCProduct> hepMCTok_;
54  edm::EDGetTokenT<reco::PFClusterCollection> pfClusterECALTok_;
55  edm::EDGetTokenT<reco::PFClusterCollection> pfClusterHCALTok_;
56  edm::EDGetTokenT<reco::PFClusterCollection> pfClusterHOTok_;
57  edm::EDGetTokenT<reco::PFClusterCollection> pfClusterHFTok_;
58 
59  MonitorElement* emean_vs_eta_E_;
60  MonitorElement* emean_vs_eta_H_;
61  MonitorElement* emean_vs_eta_EH_;
62 
63  MonitorElement* emean_vs_eta_HF_;
64  MonitorElement* emean_vs_eta_HO_;
65  MonitorElement* emean_vs_eta_EHF_;
66  MonitorElement* emean_vs_eta_EHFO_;
67 
68  MonitorElement* ratio_Esummed_ECAL_0_;
69  MonitorElement* ratio_Esummed_HCAL_0_;
70  MonitorElement* ratio_Esummed_HO_0_;
71 
72  MonitorElement* ratio_Esummed_ECAL_1_;
73  MonitorElement* ratio_Esummed_HCAL_1_;
74  MonitorElement* ratio_Esummed_HO_1_;
75 
76  MonitorElement* ratio_Esummed_ECAL_2_;
77  MonitorElement* ratio_Esummed_HCAL_2_;
78  MonitorElement* ratio_Esummed_HO_2_;
79 
80  MonitorElement* ratio_Esummed_ECAL_3_;
81  MonitorElement* ratio_Esummed_HCAL_3_;
82  MonitorElement* ratio_Esummed_HO_3_;
83 
84  MonitorElement* ratio_Esummed_ECAL_4_;
85  MonitorElement* ratio_Esummed_HCAL_4_;
86  MonitorElement* ratio_Esummed_HO_4_;
87 
88  MonitorElement* ratio_Esummed_HF_5_;
89  MonitorElement* ratio_Esummed_HF_6_;
90 
91  MonitorElement* ratio_Esummed_ECAL_HCAL_0_;
92  MonitorElement* ratio_Esummed_ECAL_HCAL_HO_0_;
93  MonitorElement* ratio_Esummed_ECAL_HCAL_1_;
94  MonitorElement* ratio_Esummed_ECAL_HCAL_HO_1_;
95  MonitorElement* ratio_Esummed_ECAL_HCAL_2_;
96  MonitorElement* ratio_Esummed_ECAL_HCAL_HO_2_;
97  MonitorElement* ratio_Esummed_ECAL_HCAL_3_;
98  MonitorElement* ratio_Esummed_ECAL_HCAL_HO_3_;
99  MonitorElement* ratio_Esummed_ECAL_HCAL_4_;
100  MonitorElement* ratio_Esummed_ECAL_HCAL_HO_4_;
101 
102  MonitorElement* egen_MC_;
103 };
104 
105 PFClusterValidation::PFClusterValidation(const edm::ParameterSet& conf) {
106  hepMCTok_ = consumes<edm::HepMCProduct>(edm::InputTag("generatorSmeared"));
107 
108  pfClusterECALTok_ =
109  consumes<reco::PFClusterCollection>(conf.getUntrackedParameter<edm::InputTag>("pflowClusterECAL"));
110  pfClusterHCALTok_ =
111  consumes<reco::PFClusterCollection>(conf.getUntrackedParameter<edm::InputTag>("pflowClusterHCAL"));
112  pfClusterHOTok_ = consumes<reco::PFClusterCollection>(conf.getUntrackedParameter<edm::InputTag>("pflowClusterHO"));
113  pfClusterHFTok_ = consumes<reco::PFClusterCollection>(conf.getUntrackedParameter<edm::InputTag>("pflowClusterHF"));
114 }
115 
116 PFClusterValidation::~PFClusterValidation() {}
117 
118 void PFClusterValidation::bookHistograms(DQMStore::IBooker& ibooker,
119  edm::Run const& irun,
120  edm::EventSetup const& isetup) {
121  constexpr auto size = 100;
122  char histo[size];
123 
124  constexpr double etaBinsOffset[] = {
125  -5.191, -4.889, -4.716, -4.538, -4.363, -4.191, -4.013, -3.839, -3.664, -3.489, -3.314, -3.139, -2.964, -2.853,
126  -2.65, -2.5, -2.322, -2.172, -2.043, -1.93, -1.83, -1.74, -1.653, -1.566, -1.479, -1.392, -1.305, -1.218,
127  -1.131, -1.044, -0.957, -0.879, -0.783, -0.696, -0.609, -0.522, -0.435, -0.348, -0.261, -0.174, -0.087, 0,
128  0.087, 0.174, 0.261, 0.348, 0.435, 0.522, 0.609, 0.696, 0.783, 0.879, 0.957, 1.044, 1.131, 1.218,
129  1.305, 1.392, 1.479, 1.566, 1.653, 1.74, 1.83, 1.93, 2.043, 2.172, 2.322, 2.5, 2.65, 2.853,
130  2.964, 3.139, 3.314, 3.489, 3.664, 3.839, 4.013, 4.191, 4.363, 4.538, 4.716, 4.889, 5.191};
131  constexpr int etaBins = std::size(etaBinsOffset) - 1;
132 
133  ibooker.setCurrentFolder("ParticleFlow/PFClusterV");
134 
135  // These are the single pion scan histos
136 
137  strncpy(histo, "emean_vs_eta_E", size);
138  emean_vs_eta_E_ = ibooker.bookProfile(histo, histo, etaBins, etaBinsOffset, -100., 2000., " ");
139  strncpy(histo, "emean_vs_eta_H", size);
140  emean_vs_eta_H_ = ibooker.bookProfile(histo, histo, etaBins, etaBinsOffset, -100., 2000., " ");
141  strncpy(histo, "emean_vs_eta_EH", size);
142  emean_vs_eta_EH_ = ibooker.bookProfile(histo, histo, etaBins, etaBinsOffset, -100., 2000., " ");
143 
144  strncpy(histo, "emean_vs_eta_HF", size);
145  emean_vs_eta_HF_ = ibooker.bookProfile(histo, histo, etaBins, etaBinsOffset, -100., 2000., " ");
146  strncpy(histo, "emean_vs_eta_HO", size);
147  emean_vs_eta_HO_ = ibooker.bookProfile(histo, histo, etaBins, etaBinsOffset, -100., 2000., " ");
148 
149  strncpy(histo, "emean_vs_eta_EHF", size);
150  emean_vs_eta_EHF_ = ibooker.bookProfile(histo, histo, etaBins, etaBinsOffset, -100., 2000., " ");
151  strncpy(histo, "emean_vs_eta_EHFO", size);
152  emean_vs_eta_EHFO_ = ibooker.bookProfile(histo, histo, etaBins, etaBinsOffset, -100., 2000., " ");
153 
154  // 1D histos
155 
156  strncpy(histo, "Ratio_Esummed_ECAL_0", size);
157  ratio_Esummed_ECAL_0_ = ibooker.book1D(histo, histo, 50, 0., 5.);
158  strncpy(histo, "Ratio_Esummed_HCAL_0", size);
159  ratio_Esummed_HCAL_0_ = ibooker.book1D(histo, histo, 50, 0., 5.);
160  strncpy(histo, "Ratio_Esummed_HO_0", size);
161  ratio_Esummed_HO_0_ = ibooker.book1D(histo, histo, 50, 0., 5.);
162  strncpy(histo, "Ratio_Esummed_ECAL_HCAL_0", size);
163  ratio_Esummed_ECAL_HCAL_0_ = ibooker.book1D(histo, histo, 50, 0., 5.);
164  strncpy(histo, "Ratio_Esummed_ECAL_HCAL_HO_0", size);
165  ratio_Esummed_ECAL_HCAL_HO_0_ = ibooker.book1D(histo, histo, 50, 0., 5.);
166 
167  strncpy(histo, "Ratio_Esummed_ECAL_1", size);
168  ratio_Esummed_ECAL_1_ = ibooker.book1D(histo, histo, 50, 0., 5.);
169  strncpy(histo, "Ratio_Esummed_HCAL_1", size);
170  ratio_Esummed_HCAL_1_ = ibooker.book1D(histo, histo, 50, 0., 5.);
171  strncpy(histo, "Ratio_Esummed_HO_1", size);
172  ratio_Esummed_HO_1_ = ibooker.book1D(histo, histo, 50, 0., 5.);
173  strncpy(histo, "Ratio_Esummed_ECAL_HCAL_1", size);
174  ratio_Esummed_ECAL_HCAL_1_ = ibooker.book1D(histo, histo, 50, 0., 5.);
175  strncpy(histo, "Ratio_Esummed_ECAL_HCAL_HO_1", size);
176  ratio_Esummed_ECAL_HCAL_HO_1_ = ibooker.book1D(histo, histo, 50, 0., 5.);
177 
178  strncpy(histo, "Ratio_Esummed_ECAL_2", size);
179  ratio_Esummed_ECAL_2_ = ibooker.book1D(histo, histo, 50, 0., 5.);
180  strncpy(histo, "Ratio_Esummed_HCAL_2", size);
181  ratio_Esummed_HCAL_2_ = ibooker.book1D(histo, histo, 50, 0., 5.);
182  strncpy(histo, "Ratio_Esummed_HO_2", size);
183  ratio_Esummed_HO_2_ = ibooker.book1D(histo, histo, 50, 0., 5.);
184  strncpy(histo, "Ratio_Esummed_ECAL_HCAL_2", size);
185  ratio_Esummed_ECAL_HCAL_2_ = ibooker.book1D(histo, histo, 50, 0., 5.);
186  strncpy(histo, "Ratio_Esummed_ECAL_HCAL_HO_2", size);
187  ratio_Esummed_ECAL_HCAL_HO_2_ = ibooker.book1D(histo, histo, 50, 0., 5.);
188 
189  strncpy(histo, "Ratio_Esummed_ECAL_3", size);
190  ratio_Esummed_ECAL_3_ = ibooker.book1D(histo, histo, 50, 0., 5.);
191  strncpy(histo, "Ratio_Esummed_HCAL_3", size);
192  ratio_Esummed_HCAL_3_ = ibooker.book1D(histo, histo, 50, 0., 5.);
193  strncpy(histo, "Ratio_Esummed_HO_3", size);
194  ratio_Esummed_HO_3_ = ibooker.book1D(histo, histo, 50, 0., 5.);
195  strncpy(histo, "Ratio_Esummed_ECAL_HCAL_3", size);
196  ratio_Esummed_ECAL_HCAL_3_ = ibooker.book1D(histo, histo, 50, 0., 5.);
197  strncpy(histo, "Ratio_Esummed_ECAL_HCAL_HO_3", size);
198  ratio_Esummed_ECAL_HCAL_HO_3_ = ibooker.book1D(histo, histo, 50, 0., 5.);
199 
200  strncpy(histo, "Ratio_Esummed_ECAL_4", size);
201  ratio_Esummed_ECAL_4_ = ibooker.book1D(histo, histo, 50, 0., 5.);
202  strncpy(histo, "Ratio_Esummed_HCAL_4", size);
203  ratio_Esummed_HCAL_4_ = ibooker.book1D(histo, histo, 50, 0., 5.);
204  strncpy(histo, "Ratio_Esummed_HO_4", size);
205  ratio_Esummed_HO_4_ = ibooker.book1D(histo, histo, 50, 0., 5.);
206  strncpy(histo, "Ratio_Esummed_ECAL_HCAL_4", size);
207  ratio_Esummed_ECAL_HCAL_4_ = ibooker.book1D(histo, histo, 50, 0., 5.);
208  strncpy(histo, "Ratio_Esummed_ECAL_HCAL_HO_4", size);
209  ratio_Esummed_ECAL_HCAL_HO_4_ = ibooker.book1D(histo, histo, 50, 0., 5.);
210 
211  strncpy(histo, "Ratio_Esummed_HF_5", size);
212  ratio_Esummed_HF_5_ = ibooker.book1D(histo, histo, 50, 0., 5.);
213  strncpy(histo, "Ratio_Esummed_HF_6", size);
214  ratio_Esummed_HF_6_ = ibooker.book1D(histo, histo, 50, 0., 5.);
215 
216  strncpy(histo, "Egen_MC", size);
217  egen_MC_ = ibooker.book1D(histo, histo, 50, 0, 50);
218 }
219 
220 void PFClusterValidation::analyze(edm::Event const& event, edm::EventSetup const& c) {
221  double eta_MC = 0.;
222  double phi_MC = 0.;
223  double energy_MC = 0.;
224 
225  edm::Handle<edm::HepMCProduct> hepMC;
226  event.getByToken(hepMCTok_, hepMC);
227  if (not hepMC.isValid()) {
228  edm::LogWarning("PFClusterValidation") << "HepMCProduct not found";
229  return;
230  }
231 
232  // MC particle with highest pt is taken as a direction reference
233  double maxPt = -99999.;
234  const HepMC::GenEvent* myGenEvent = hepMC->GetEvent();
235  for (auto p = myGenEvent->particles_begin(); p != myGenEvent->particles_end(); ++p) {
236  double phip = (*p)->momentum().phi();
237  double etap = (*p)->momentum().eta();
238  double pt = (*p)->momentum().perp();
239  double energy = (*p)->momentum().e();
240  if (pt > maxPt) {
241  maxPt = pt;
242  energy_MC = energy;
243  phi_MC = phip;
244  eta_MC = etap;
245  }
246  }
247 
248  egen_MC_->Fill(energy_MC);
249 
250  edm::Handle<reco::PFClusterCollection> pfClusterECAL;
251  event.getByToken(pfClusterECALTok_, pfClusterECAL);
252 
253  edm::Handle<reco::PFClusterCollection> pfClusterHCAL;
254  event.getByToken(pfClusterHCALTok_, pfClusterHCAL);
255 
256  edm::Handle<reco::PFClusterCollection> pfClusterHO;
257  event.getByToken(pfClusterHOTok_, pfClusterHO);
258 
259  edm::Handle<reco::PFClusterCollection> pfClusterHF;
260  event.getByToken(pfClusterHFTok_, pfClusterHF);
261 
262  // sum the energy in a dR cone for each subsystem
263  const double Econe = sumEnergy(pfClusterECAL, eta_MC, phi_MC);
264  const double Hcone = sumEnergy(pfClusterHCAL, eta_MC, phi_MC);
265  const double HOcone = sumEnergy(pfClusterHO, eta_MC, phi_MC);
266  const double HFcone = sumEnergy(pfClusterHF, eta_MC, phi_MC);
267 
268  if (energy_MC > 0.) {
269  if (std::abs(eta_MC) < 0.5) {
270  ratio_Esummed_ECAL_0_->Fill(Econe / energy_MC);
271  ratio_Esummed_HCAL_0_->Fill(Hcone / energy_MC);
272  ratio_Esummed_HO_0_->Fill(HOcone / energy_MC);
273  ratio_Esummed_ECAL_HCAL_0_->Fill((Econe + Hcone) / energy_MC);
274  ratio_Esummed_ECAL_HCAL_HO_0_->Fill((Econe + Hcone + HOcone) / energy_MC);
275  } else if (std::abs(eta_MC) < 1.3 && std::abs(eta_MC) > 0.5) {
276  ratio_Esummed_ECAL_1_->Fill(Econe / energy_MC);
277  ratio_Esummed_HCAL_1_->Fill(Hcone / energy_MC);
278  ratio_Esummed_HO_1_->Fill(HOcone / energy_MC);
279  ratio_Esummed_ECAL_HCAL_1_->Fill((Econe + Hcone) / energy_MC);
280  ratio_Esummed_ECAL_HCAL_HO_1_->Fill((Econe + Hcone + HOcone) / energy_MC);
281  } else if (std::abs(eta_MC) < 2.1 && std::abs(eta_MC) > 1.3) {
282  ratio_Esummed_ECAL_2_->Fill(Econe / energy_MC);
283  ratio_Esummed_HCAL_2_->Fill(Hcone / energy_MC);
284  ratio_Esummed_HO_2_->Fill(HOcone / energy_MC);
285  ratio_Esummed_ECAL_HCAL_2_->Fill((Econe + Hcone) / energy_MC);
286  ratio_Esummed_ECAL_HCAL_HO_2_->Fill((Econe + Hcone + HOcone) / energy_MC);
287  } else if (std::abs(eta_MC) < 2.5 && std::abs(eta_MC) > 2.1) {
288  ratio_Esummed_ECAL_3_->Fill(Econe / energy_MC);
289  ratio_Esummed_HCAL_3_->Fill(Hcone / energy_MC);
290  ratio_Esummed_HO_3_->Fill(HOcone / energy_MC);
291  ratio_Esummed_ECAL_HCAL_3_->Fill((Econe + Hcone) / energy_MC);
292  ratio_Esummed_ECAL_HCAL_HO_3_->Fill((Econe + Hcone + HOcone) / energy_MC);
293  } else if (2.5 < std::abs(eta_MC) && std::abs(eta_MC) < 3.0) {
294  ratio_Esummed_ECAL_4_->Fill(Econe / energy_MC);
295  ratio_Esummed_HCAL_4_->Fill(Hcone / energy_MC);
296  ratio_Esummed_HO_4_->Fill(HOcone / energy_MC);
297  ratio_Esummed_ECAL_HCAL_4_->Fill((Econe + Hcone) / energy_MC);
298  ratio_Esummed_ECAL_HCAL_HO_4_->Fill((Econe + Hcone + HOcone) / energy_MC);
299  } else if (3.0 < std::abs(eta_MC) && std::abs(eta_MC) < 4.0) {
300  ratio_Esummed_HF_5_->Fill(HFcone / energy_MC);
301  } else if (4.0 < std::abs(eta_MC) && std::abs(eta_MC) < 5.0) {
302  ratio_Esummed_HF_6_->Fill(HFcone / energy_MC);
303  }
304  }
305 
306  emean_vs_eta_E_->Fill(eta_MC, Econe);
307  emean_vs_eta_H_->Fill(eta_MC, Hcone);
308  emean_vs_eta_EH_->Fill(eta_MC, Econe + Hcone);
309  emean_vs_eta_HF_->Fill(eta_MC, HFcone);
310  emean_vs_eta_HO_->Fill(eta_MC, HOcone);
311  emean_vs_eta_EHF_->Fill(eta_MC, Econe + Hcone + HFcone);
312  emean_vs_eta_EHFO_->Fill(eta_MC, Econe + Hcone + HFcone + HOcone);
313 }
314 
315 double PFClusterValidation::sumEnergy(edm::Handle<reco::PFClusterCollection> const& pfClusters,
316  double eta,
317  double phi) {
318  if (not pfClusters.isValid())
319  return 0.;
320 
321  double sum = 0.;
322  for (auto pf = pfClusters->begin(); pf != pfClusters->end(); ++pf) {
323  if (reco::deltaR2(eta, phi, pf->eta(), pf->phi()) < partR2) {
324  sum += pf->energy();
325  }
326  }
327 
328  return sum;
329 }
330 
331 #include "FWCore/Framework/interface/MakerMacros.h"
332 DEFINE_FWK_MODULE(PFClusterValidation);
dqm::impl::MonitorElement
Definition: MonitorElement.h:98
Handle.h
PFClusterValidation::ratio_Esummed_HCAL_1_
MonitorElement * ratio_Esummed_HCAL_1_
Definition: PFClusterValidation.cc:73
PFClusterValidation::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: PFClusterValidation.cc:118
MessageLogger.h
PFClusterValidation::partR2
static constexpr double partR2
Definition: PFClusterValidation.cc:50
ESHandle.h
PFBlockElementCluster.h
PFCandidate.h
PFClusterValidation::hepMCTok_
edm::EDGetTokenT< edm::HepMCProduct > hepMCTok_
Definition: PFClusterValidation.cc:53
DiDispStaMuonMonitor_cfi.pt
pt
Definition: DiDispStaMuonMonitor_cfi.py:39
edm::Run
Definition: Run.h:45
edm::EDGetTokenT< edm::HepMCProduct >
PFClusterValidation::ratio_Esummed_HCAL_4_
MonitorElement * ratio_Esummed_HCAL_4_
Definition: PFClusterValidation.cc:85
PFClusterValidation::ratio_Esummed_HO_2_
MonitorElement * ratio_Esummed_HO_2_
Definition: PFClusterValidation.cc:78
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
HLT_FULL_cff.etaBins
etaBins
Definition: HLT_FULL_cff.py:15288
PFClusterValidation::ratio_Esummed_ECAL_HCAL_HO_1_
MonitorElement * ratio_Esummed_ECAL_HCAL_HO_1_
Definition: PFClusterValidation.cc:94
timingPdfMaker.histo
histo
Definition: timingPdfMaker.py:278
DQMStore.h
PFClusterValidation::ratio_Esummed_ECAL_HCAL_1_
MonitorElement * ratio_Esummed_ECAL_HCAL_1_
Definition: PFClusterValidation.cc:93
L1TrackerEtMissProducer_cfi.maxPt
maxPt
Definition: L1TrackerEtMissProducer_cfi.py:14
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
EDAnalyzer.h
PFClusterValidation::emean_vs_eta_HF_
MonitorElement * emean_vs_eta_HF_
Definition: PFClusterValidation.cc:63
PFClusterValidation::ratio_Esummed_HO_1_
MonitorElement * ratio_Esummed_HO_1_
Definition: PFClusterValidation.cc:74
PFClusterValidation::ratio_Esummed_HCAL_3_
MonitorElement * ratio_Esummed_HCAL_3_
Definition: PFClusterValidation.cc:81
edm::Handle
Definition: AssociativeIterator.h:50
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
PFRecHit.h
HepMC::GenEvent
Definition: hepmc_rootio.cc:9
PFClusterValidation::ratio_Esummed_ECAL_HCAL_4_
MonitorElement * ratio_Esummed_ECAL_HCAL_4_
Definition: PFClusterValidation.cc:99
defaults_cfi.etaBinsOffset
etaBinsOffset
Definition: defaults_cfi.py:23
PFLayer.h
MakerMacros.h
PFClusterValidation::analyze
void analyze(edm::Event const &e, edm::EventSetup const &c) override
Definition: PFClusterValidation.cc:220
PFClusterValidation::ratio_Esummed_ECAL_HCAL_0_
MonitorElement * ratio_Esummed_ECAL_HCAL_0_
Definition: PFClusterValidation.cc:91
PFClusterValidation::emean_vs_eta_EHFO_
MonitorElement * emean_vs_eta_EHFO_
Definition: PFClusterValidation.cc:66
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
PFClusterValidation::ratio_Esummed_HCAL_0_
MonitorElement * ratio_Esummed_HCAL_0_
Definition: PFClusterValidation.cc:69
PFClusterValidation::ratio_Esummed_ECAL_HCAL_2_
MonitorElement * ratio_Esummed_ECAL_HCAL_2_
Definition: PFClusterValidation.cc:95
Service.h
PVValHelper::eta
Definition: PVValidationHelpers.h:70
PFClusterValidation::ratio_Esummed_ECAL_4_
MonitorElement * ratio_Esummed_ECAL_4_
Definition: PFClusterValidation.cc:84
PFClusterValidation::ratio_Esummed_ECAL_HCAL_HO_2_
MonitorElement * ratio_Esummed_ECAL_HCAL_HO_2_
Definition: PFClusterValidation.cc:96
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
dqm::implementation::IBooker::bookProfile
MonitorElement * bookProfile(TString const &name, TString const &title, int nchX, double lowX, double highX, int, double lowY, double highY, char const *option="s", FUNC onbooking=NOOP())
Definition: DQMStore.h:322
PFCluster.h
PFClusterValidation::pfClusterHFTok_
edm::EDGetTokenT< reco::PFClusterCollection > pfClusterHFTok_
Definition: PFClusterValidation.cc:57
CaloClusterFwd.h
HCALHighEnergyHPDFilter_cfi.energy
energy
Definition: HCALHighEnergyHPDFilter_cfi.py:5
PFClusterValidation::emean_vs_eta_E_
MonitorElement * emean_vs_eta_E_
Definition: PFClusterValidation.cc:59
PFClusterValidation::egen_MC_
MonitorElement * egen_MC_
Definition: PFClusterValidation.cc:102
PFClusterValidation::ratio_Esummed_ECAL_HCAL_HO_0_
MonitorElement * ratio_Esummed_ECAL_HCAL_HO_0_
Definition: PFClusterValidation.cc:92
DQMEDAnalyzer.h
EcalSubdetector.h
DQMEDAnalyzer
Definition: DQMEDAnalyzer.py:1
PFClusterValidation::ratio_Esummed_HF_6_
MonitorElement * ratio_Esummed_HF_6_
Definition: PFClusterValidation.cc:89
edm::ParameterSet
Definition: ParameterSet.h:47
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
Event.h
PVValHelper::phi
Definition: PVValidationHelpers.h:69
deltaR.h
PackedCandidate.h
HcalDetId.h
reco::deltaR2
constexpr auto deltaR2(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:16
PFClusterValidation::ratio_Esummed_HCAL_2_
MonitorElement * ratio_Esummed_HCAL_2_
Definition: PFClusterValidation.cc:77
ModuleDef.h
PFClusterFwd.h
PFClusterValidation::ratio_Esummed_ECAL_HCAL_3_
MonitorElement * ratio_Esummed_ECAL_HCAL_3_
Definition: PFClusterValidation.cc:97
PFClusterValidation::PFClusterValidation
PFClusterValidation(edm::ParameterSet const &conf)
Definition: PFClusterValidation.cc:105
PFClusterValidation::pfClusterHOTok_
edm::EDGetTokenT< reco::PFClusterCollection > pfClusterHOTok_
Definition: PFClusterValidation.cc:56
PFClusterValidation::ratio_Esummed_HF_5_
MonitorElement * ratio_Esummed_HF_5_
Definition: PFClusterValidation.cc:88
CaloTowerDetId.h
edm::EventSetup
Definition: EventSetup.h:58
PFClusterValidation::ratio_Esummed_ECAL_HCAL_HO_4_
MonitorElement * ratio_Esummed_ECAL_HCAL_HO_4_
Definition: PFClusterValidation.cc:100
HcalSubdetector.h
PFClusterValidation::pfClusterECALTok_
edm::EDGetTokenT< reco::PFClusterCollection > pfClusterECALTok_
Definition: PFClusterValidation.cc:54
PFClusterValidation
Definition: PFClusterValidation.cc:42
PFClusterValidation::pfClusterHCALTok_
edm::EDGetTokenT< reco::PFClusterCollection > pfClusterHCALTok_
Definition: PFClusterValidation.cc:55
PFClusterValidation::emean_vs_eta_EHF_
MonitorElement * emean_vs_eta_EHF_
Definition: PFClusterValidation.cc:65
PFRecHitFraction.h
PFClusterValidation::~PFClusterValidation
~PFClusterValidation() override
Definition: PFClusterValidation.cc:116
CaloTowerCollection.h
PFClusterValidation::emean_vs_eta_HO_
MonitorElement * emean_vs_eta_HO_
Definition: PFClusterValidation.cc:64
PFBlock.h
DDAxes::phi
PFClusterValidation::ratio_Esummed_ECAL_1_
MonitorElement * ratio_Esummed_ECAL_1_
Definition: PFClusterValidation.cc:72
PFClusterValidation::emean_vs_eta_H_
MonitorElement * emean_vs_eta_H_
Definition: PFClusterValidation.cc:60
DetId.h
Frameworkfwd.h
packedPFCandidateRefMixer_cfi.pf
pf
Definition: packedPFCandidateRefMixer_cfi.py:4
PFClusterValidation::ratio_Esummed_ECAL_3_
MonitorElement * ratio_Esummed_ECAL_3_
Definition: PFClusterValidation.cc:80
PFClusterValidation::ratio_Esummed_HO_4_
MonitorElement * ratio_Esummed_HO_4_
Definition: PFClusterValidation.cc:86
EventSetup.h
dqm::implementation::IBooker
Definition: DQMStore.h:43
PFClusterValidation::ratio_Esummed_HO_3_
MonitorElement * ratio_Esummed_HO_3_
Definition: PFClusterValidation.cc:82
PFClusterValidation::ratio_Esummed_HO_0_
MonitorElement * ratio_Esummed_HO_0_
Definition: PFClusterValidation.cc:70
PFClusterValidation::emean_vs_eta_EH_
MonitorElement * emean_vs_eta_EH_
Definition: PFClusterValidation.cc:61
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
ParameterSet.h
HepMCProduct.h
c
auto & c
Definition: CAHitNtupletGeneratorKernelsImpl.h:56
PFClusterValidation::ratio_Esummed_ECAL_2_
MonitorElement * ratio_Esummed_ECAL_2_
Definition: PFClusterValidation.cc:76
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
event
Definition: event.py:1
edm::Event
Definition: Event.h:73
PFClusterValidation::sumEnergy
static double sumEnergy(edm::Handle< reco::PFClusterCollection > const &pfClusters, double eta, double phi)
Definition: PFClusterValidation.cc:315
Vector3D.h
edm::InputTag
Definition: InputTag.h:15
PFBlockElementTrack.h
MuonClassifier_cff.hepMC
hepMC
Definition: MuonClassifier_cff.py:13
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
PFClusterValidation::ratio_Esummed_ECAL_0_
MonitorElement * ratio_Esummed_ECAL_0_
Definition: PFClusterValidation.cc:68
PFClusterValidation::ratio_Esummed_ECAL_HCAL_HO_3_
MonitorElement * ratio_Esummed_ECAL_HCAL_HO_3_
Definition: PFClusterValidation.cc:98
CaloCluster.h
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
Generated for CMSSW Reference Manual by doxygen 1.8.16