CMS 3D CMS Logo
PhotonValidatorMiniAOD.cc
Go to the documentation of this file.
1 #include <memory>
2 #include "Validation/RecoEgamma/plugins/PhotonValidatorMiniAOD.h"
3 //
4 #include "FWCore/Framework/interface/Frameworkfwd.h"
5 #include "FWCore/Framework/interface/EDAnalyzer.h"
6 #include "FWCore/MessageLogger/interface/MessageLogger.h"
7 #include "FWCore/Common/interface/TriggerNames.h"
8 #include "FWCore/Framework/interface/Event.h"
9 #include "FWCore/Framework/interface/MakerMacros.h"
10 #include "FWCore/Utilities/interface/InputTag.h"
11 #include "FWCore/ParameterSet/interface/ParameterSet.h"
12 //#include "CommonTools/UtilAlgos/interface/TFileService.h"
13 #include "FWCore/ServiceRegistry/interface/Service.h"
14 #include "DataFormats/Common/interface/Handle.h"
15 #include "DataFormats/Common/interface/Ptr.h"
16 //#include "DataFormats/Common/interface/PtrVector.h"
17 #include "DataFormats/VertexReco/interface/Vertex.h"
18 #include "DataFormats/Common/interface/TriggerResults.h"
19 //#include "DataFormats/EgammaCandidates/interface/PhotonCore.h"
20 //#include "DataFormats/EgammaCandidates/interface/Photon.h"
21 #include "DataFormats/PatCandidates/interface/Photon.h"
22 #include "DataFormats/PatCandidates/interface/PackedGenParticle.h"
23 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
24 
25 #include "CLHEP/Units/GlobalPhysicalConstants.h"
26 #include "TTree.h"
27 #include "TFile.h"
28 #include "TH1F.h"
29 #include "TH2F.h"
30 
31 // **********************************************************************
32 
33 PhotonValidatorMiniAOD::PhotonValidatorMiniAOD(const edm::ParameterSet &iConfig)
34  : outputFileName_(iConfig.getParameter<std::string>("outputFileName")),
35  photonToken_(consumes<edm::View<pat::Photon> >(
36  iConfig.getUntrackedParameter<edm::InputTag>("PhotonTag", edm::InputTag("slimmedPhotons")))),
37  genpartToken_(consumes<reco::GenParticleCollection>(
38  iConfig.getUntrackedParameter<edm::InputTag>("genpartTag", edm::InputTag("prunedGenParticles"))))
39 
40 {
41  parameters_ = iConfig;
42 }
43 
44 PhotonValidatorMiniAOD::~PhotonValidatorMiniAOD() {}
45 
46 void PhotonValidatorMiniAOD::bookHistograms(DQMStore::IBooker &iBooker,
47  edm::Run const &run,
48  edm::EventSetup const &es) {
49  double eMin = parameters_.getParameter<double>("eMin");
50  double eMax = parameters_.getParameter<double>("eMax");
51  int eBin = parameters_.getParameter<int>("eBin");
52 
53  double etMin = parameters_.getParameter<double>("etMin");
54  double etMax = parameters_.getParameter<double>("etMax");
55  int etBin = parameters_.getParameter<int>("etBin");
56 
57  double resMin = parameters_.getParameter<double>("resMin");
58  double resMax = parameters_.getParameter<double>("resMax");
59  int resBin = parameters_.getParameter<int>("resBin");
60 
61  double etaMin = parameters_.getParameter<double>("etaMin");
62  double etaMax = parameters_.getParameter<double>("etaMax");
63  int etaBin = parameters_.getParameter<int>("etaBin");
64 
65  double phiMin = parameters_.getParameter<double>("phiMin");
66  double phiMax = parameters_.getParameter<double>("phiMax");
67  int phiBin = parameters_.getParameter<int>("phiBin");
68 
69  double r9Min = parameters_.getParameter<double>("r9Min");
70  double r9Max = parameters_.getParameter<double>("r9Max");
71  int r9Bin = parameters_.getParameter<int>("r9Bin");
72 
73  iBooker.setCurrentFolder("EgammaV/PhotonValidatorMiniAOD/Photons");
74 
75  h_scEta_[0] = iBooker.book1D("scEta_miniAOD", " SC Eta ", etaBin, etaMin, etaMax);
76  h_scPhi_[0] = iBooker.book1D("scPhi_miniAOD", " SC Phi ", phiBin, phiMin, phiMax);
77 
78  std::string histname = " ";
79 
80  histname = "phoE";
81  h_phoE_[0][0] = iBooker.book1D(histname + "All_miniAOD", " Photon Energy: All ecal ", eBin, eMin, eMax);
82  h_phoE_[0][1] = iBooker.book1D(histname + "Barrel_miniAOD", " Photon Energy: barrel ", eBin, eMin, eMax);
83  h_phoE_[0][2] = iBooker.book1D(histname + "Endcap_miniAOD", " Photon Energy: Endcap ", eBin, eMin, eMax);
84 
85  histname = "phoEt";
86  h_phoEt_[0][0] =
87  iBooker.book1D(histname + "All_miniAOD", " Photon Transverse Energy: All ecal ", etBin, etMin, etMax);
88  h_phoEt_[0][1] =
89  iBooker.book1D(histname + "Barrel_miniAOD", " Photon Transverse Energy: Barrel ", etBin, etMin, etMax);
90  h_phoEt_[0][2] =
91  iBooker.book1D(histname + "Endcap_miniAOD", " Photon Transverse Energy: Endcap ", etBin, etMin, etMax);
92 
93  histname = "eRes";
94  h_phoERes_[0][0] = iBooker.book1D(
95  histname + "All_miniAOD", " Photon E/E_{true}: All ecal; E/E_{true} (GeV)", resBin, resMin, resMax);
96  h_phoERes_[0][1] = iBooker.book1D(
97  histname + "Barrel_miniAOD", "Photon E/E_{true}: Barrel; E/E_{true} (GeV)", resBin, resMin, resMax);
98  h_phoERes_[0][2] = iBooker.book1D(
99  histname + "Endcap_miniAOD", " Photon E/E_{true}: Endcap; E/E_{true} (GeV)", resBin, resMin, resMax);
100 
101  histname = "sigmaEoE";
102  h_phoSigmaEoE_[0][0] =
103  iBooker.book1D(histname + "All_miniAOD", "#sigma_{E}/E: All ecal; #sigma_{E}/E", 100, 0., 0.08);
104  h_phoSigmaEoE_[0][1] =
105  iBooker.book1D(histname + "Barrel_miniAOD", "#sigma_{E}/E: Barrel; #sigma_{E}/E", 100, 0., 0.08);
106  h_phoSigmaEoE_[0][2] =
107  iBooker.book1D(histname + "Endcap_miniAOD", "#sigma_{E}/E: Endcap, #sigma_{E}/E", 100, 0., 0.08);
108 
109  histname = "r9";
110  h_r9_[0][0] = iBooker.book1D(histname + "All_miniAOD", " r9: All Ecal", r9Bin, r9Min, r9Max);
111  h_r9_[0][1] = iBooker.book1D(histname + "Barrel_miniAOD", " r9: Barrel ", r9Bin, r9Min, r9Max);
112  h_r9_[0][2] = iBooker.book1D(histname + "Endcap_miniAOD", " r9: Endcap ", r9Bin, r9Min, r9Max);
113  histname = "full5x5_r9";
114  h_full5x5_r9_[0][0] = iBooker.book1D(histname + "All_miniAOD", " r9: All Ecal", r9Bin, r9Min, r9Max);
115  h_full5x5_r9_[0][1] = iBooker.book1D(histname + "Barrel_miniAOD", " r9: Barrel ", r9Bin, r9Min, r9Max);
116  h_full5x5_r9_[0][2] = iBooker.book1D(histname + "Endcap_miniAOD", " r9: Endcap ", r9Bin, r9Min, r9Max);
117  histname = "r1";
118  h_r1_[0][0] = iBooker.book1D(histname + "All_miniAOD", " e1x5/e5x5: All Ecal", r9Bin, r9Min, r9Max);
119  h_r1_[0][1] = iBooker.book1D(histname + "Barrel_miniAOD", " e1x5/e5x5: Barrel ", r9Bin, r9Min, r9Max);
120  h_r1_[0][2] = iBooker.book1D(histname + "Endcap_miniAOD", " e1x5/e5x5: Endcap ", r9Bin, r9Min, r9Max);
121  histname = "r2";
122  h_r2_[0][0] = iBooker.book1D(histname + "All_miniAOD", " e2x5/e5x5: All Ecal", r9Bin, r9Min, r9Max);
123  h_r2_[0][1] = iBooker.book1D(histname + "Barrel_miniAOD", " e2x5/e5x5: Barrel ", r9Bin, r9Min, r9Max);
124  h_r2_[0][2] = iBooker.book1D(histname + "Endcap_miniAOD", " e2x5/e5x5: Endcap ", r9Bin, r9Min, r9Max);
125  histname = "hOverE";
126  h_hOverE_[0][0] = iBooker.book1D(histname + "All_miniAOD", "H/E: All Ecal", 100, 0., 0.2);
127  h_hOverE_[0][1] = iBooker.book1D(histname + "Barrel_miniAOD", "H/E: Barrel ", 100, 0., 0.2);
128  h_hOverE_[0][2] = iBooker.book1D(histname + "Endcap_miniAOD", "H/E: Endcap ", 100, 0., 0.2);
129  //
130  histname = "newhOverE";
131  h_newhOverE_[0][0] = iBooker.book1D(histname + "All_miniAOD", "new H/E: All Ecal", 100, 0., 0.2);
132  h_newhOverE_[0][1] = iBooker.book1D(histname + "Barrel_miniAOD", "new H/E: Barrel ", 100, 0., 0.2);
133  h_newhOverE_[0][2] = iBooker.book1D(histname + "Endcap_miniAOD", "new H/E: Endcap ", 100, 0., 0.2);
134  //
135  histname = "sigmaIetaIeta";
136  h_sigmaIetaIeta_[0][0] = iBooker.book1D(histname + "All_miniAOD", "sigmaIetaIeta: All Ecal", 100, 0., 0.1);
137  h_sigmaIetaIeta_[0][1] = iBooker.book1D(histname + "Barrel_miniAOD", "sigmaIetaIeta: Barrel ", 100, 0., 0.05);
138  h_sigmaIetaIeta_[0][2] = iBooker.book1D(histname + "Endcap_miniAOD", "sigmaIetaIeta: Endcap ", 100, 0., 0.1);
139  histname = "full5x5_sigmaIetaIeta";
140  h_full5x5_sigmaIetaIeta_[0][0] =
141  iBooker.book1D(histname + "All_miniAOD", "Full5x5 sigmaIetaIeta: All Ecal", 100, 0., 0.1);
142  h_full5x5_sigmaIetaIeta_[0][1] =
143  iBooker.book1D(histname + "Barrel_miniAOD", "Full5x5 sigmaIetaIeta: Barrel ", 100, 0., 0.05);
144  h_full5x5_sigmaIetaIeta_[0][2] =
145  iBooker.book1D(histname + "Endcap_miniAOD", "Full5x5 sigmaIetaIeta: Endcap ", 100, 0., 0.1);
146  //
147  histname = "ecalRecHitSumEtConeDR04";
148  h_ecalRecHitSumEtConeDR04_[0][0] =
149  iBooker.book1D(histname + "All_miniAOD", "ecalRecHitSumEtDR04: All Ecal", etBin, etMin, 20.);
150  h_ecalRecHitSumEtConeDR04_[0][1] =
151  iBooker.book1D(histname + "Barrel_miniAOD", "ecalRecHitSumEtDR04: Barrel ", etBin, etMin, 20.);
152  h_ecalRecHitSumEtConeDR04_[0][2] =
153  iBooker.book1D(histname + "Endcap_miniAOD", "ecalRecHitSumEtDR04: Endcap ", etBin, etMin, 20.);
154  histname = "hcalTowerSumEtConeDR04";
155  h_hcalTowerSumEtConeDR04_[0][0] =
156  iBooker.book1D(histname + "All_miniAOD", "hcalTowerSumEtConeDR04: All Ecal", etBin, etMin, 20.);
157  h_hcalTowerSumEtConeDR04_[0][1] =
158  iBooker.book1D(histname + "Barrel_miniAOD", "hcalTowerSumEtConeDR04: Barrel ", etBin, etMin, 20.);
159  h_hcalTowerSumEtConeDR04_[0][2] =
160  iBooker.book1D(histname + "Endcap_miniAOD", "hcalTowerSumEtConeDR04: Endcap ", etBin, etMin, 20.);
161  //
162  histname = "hcalTowerBcSumEtConeDR04";
163  h_hcalTowerBcSumEtConeDR04_[0][0] =
164  iBooker.book1D(histname + "All_miniAOD", "hcalTowerBcSumEtConeDR04: All Ecal", etBin, etMin, 20.);
165  h_hcalTowerBcSumEtConeDR04_[0][1] =
166  iBooker.book1D(histname + "Barrel_miniAOD", "hcalTowerBcSumEtConeDR04: Barrel ", etBin, etMin, 20.);
167  h_hcalTowerBcSumEtConeDR04_[0][2] =
168  iBooker.book1D(histname + "Endcap_miniAOD", "hcalTowerBcSumEtConeDR04: Endcap ", etBin, etMin, 20.);
169  histname = "isoTrkSolidConeDR04";
170  h_isoTrkSolidConeDR04_[0][0] =
171  iBooker.book1D(histname + "All_miniAOD", "isoTrkSolidConeDR04: All Ecal", etBin, etMin, etMax * 0.1);
172  h_isoTrkSolidConeDR04_[0][1] =
173  iBooker.book1D(histname + "Barrel_miniAOD", "isoTrkSolidConeDR04: Barrel ", etBin, etMin, etMax * 0.1);
174  h_isoTrkSolidConeDR04_[0][2] =
175  iBooker.book1D(histname + "Endcap_miniAOD", "isoTrkSolidConeDR04: Endcap ", etBin, etMin, etMax * 0.1);
176  histname = "nTrkSolidConeDR04";
177  h_nTrkSolidConeDR04_[0][0] = iBooker.book1D(histname + "All_miniAOD", "nTrkSolidConeDR04: All Ecal", 20, 0., 20);
178  h_nTrkSolidConeDR04_[0][1] = iBooker.book1D(histname + "Barrel_miniAOD", "nTrkSolidConeDR04: Barrel ", 20, 0., 20);
179  h_nTrkSolidConeDR04_[0][2] = iBooker.book1D(histname + "Endcap_miniAOD", "nTrkSolidConeDR04: Endcap ", 20, 0., 20);
180 
181  // Infos from Particle Flow - isolation and ID
182  histname = "chargedHadIso";
183  h_chHadIso_[0] = iBooker.book1D(histname + "All_miniAOD", "PF chargedHadIso: All Ecal", etBin, etMin, 20.);
184  h_chHadIso_[1] = iBooker.book1D(histname + "Barrel_miniAOD", "PF chargedHadIso: Barrel", etBin, etMin, 20.);
185  h_chHadIso_[2] = iBooker.book1D(histname + "Endcap_miniAOD", "PF chargedHadIso: Endcap", etBin, etMin, 20.);
186  histname = "neutralHadIso";
187  h_nHadIso_[0] = iBooker.book1D(histname + "All_miniAOD", "PF neutralHadIso: All Ecal", etBin, etMin, 20.);
188  h_nHadIso_[1] = iBooker.book1D(histname + "Barrel_miniAOD", "PF neutralHadIso: Barrel", etBin, etMin, 20.);
189  h_nHadIso_[2] = iBooker.book1D(histname + "Endcap_miniAOD", "PF neutralHadIso: Endcap", etBin, etMin, 20.);
190  histname = "photonIso";
191  h_phoIso_[0] = iBooker.book1D(histname + "All_miniAOD", "PF photonIso: All Ecal", etBin, etMin, 20.);
192  h_phoIso_[1] = iBooker.book1D(histname + "Barrel_miniAOD", "PF photonIso: Barrel", etBin, etMin, 20.);
193  h_phoIso_[2] = iBooker.book1D(histname + "Endcap_miniAOD", "PF photonIso: Endcap", etBin, etMin, 20.);
194 }
195 
196 void PhotonValidatorMiniAOD::analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup) {
197  // ********************************************************************************
198  using namespace edm;
199 
200  // access edm objects
201  Handle<View<pat::Photon> > photonsHandle;
202  iEvent.getByToken(photonToken_, photonsHandle);
203  const auto &photons = *photonsHandle;
204  //
205  Handle<reco::GenParticleCollection> genParticles;
206  iEvent.getByToken(genpartToken_, genParticles);
207 
208  for (reco::GenParticleCollection::const_iterator mcIter = genParticles->begin(); mcIter != genParticles->end();
209  mcIter++) {
210  if (!(mcIter->pdgId() == 22))
211  continue;
212  if (mcIter->mother() != nullptr && !(mcIter->mother()->pdgId() == 25))
213  continue;
214  if (fabs(mcIter->eta()) > 2.5)
215  continue;
216  // if ( mcIter->pt() < 5) continue;
217  // if ( fabs(mcIter->eta()) > 1.4442 && fabs(mcIter->eta()) < 1.566) continue;
218 
219  float mcPhi = mcIter->phi();
220  float mcEta = mcIter->eta();
221  //mcEta = etaTransformation(mcEta, (*mcPho).primaryVertex().z() );
222  float mcEnergy = mcIter->energy();
223 
224  double dR = 9999999.;
225  float minDr = 10000.;
226  int iMatch = -1;
227  bool matched = false;
228  for (size_t ipho = 0; ipho < photons.size(); ipho++) {
229  Ptr<pat::Photon> pho = photons.ptrAt(ipho);
230 
231  double dphi = pho->phi() - mcPhi;
232  if (std::fabs(dphi) > CLHEP::pi) {
233  dphi = dphi < 0 ? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
234  }
235  double deta = pho->superCluster()->position().eta() - mcEta;
236 
237  dR = sqrt(pow((deta), 2) + pow(dphi, 2));
238  if (dR < 0.1 && dR < minDr) {
239  minDr = dR;
240  iMatch = ipho;
241  }
242  }
243 
244  if (iMatch > -1)
245  matched = true;
246  if (!matched)
247  continue;
248 
249  Ptr<pat::Photon> matchingPho = photons.ptrAt(iMatch);
250 
251  bool phoIsInBarrel = false;
252  bool phoIsInEndcap = false;
253 
254  float phoEta = matchingPho->superCluster()->position().eta();
255  if (fabs(phoEta) < 1.479) {
256  phoIsInBarrel = true;
257  } else {
258  phoIsInEndcap = true;
259  }
260 
261  float photonE = matchingPho->energy();
262  float sigmaEoE = matchingPho->getCorrectedEnergyError(matchingPho->getCandidateP4type()) / matchingPho->energy();
263  float photonEt = matchingPho->energy() / cosh(matchingPho->eta());
264  // float photonERegr1 = matchingPho->getCorrectedEnergy(reco::Photon::regression1);
265  //float photonERegr2 = matchingPho->getCorrectedEnergy(reco::Photon::regression2);
266  float r9 = matchingPho->r9();
267  float full5x5_r9 = matchingPho->full5x5_r9();
268  float r1 = matchingPho->r1x5();
269  float r2 = matchingPho->r2x5();
270  float sieie = matchingPho->sigmaIetaIeta();
271  float full5x5_sieie = matchingPho->full5x5_sigmaIetaIeta();
272  float hOverE = matchingPho->hadronicOverEm();
273  float newhOverE = matchingPho->hadTowOverEm();
274  float ecalIso = matchingPho->ecalRecHitSumEtConeDR04();
275  float hcalIso = matchingPho->hcalTowerSumEtConeDR04();
276  float newhcalIso = matchingPho->hcalTowerSumEtBcConeDR04();
277  float trkIso = matchingPho->trkSumPtSolidConeDR04();
278  float nIsoTrk = matchingPho->nTrkSolidConeDR04();
279  // PF related quantities
280  float chargedHadIso = matchingPho->chargedHadronIso();
281  float neutralHadIso = matchingPho->neutralHadronIso();
282  float photonIso = matchingPho->photonIso();
283  // float etOutsideMustache = matchingPho->etOutsideMustache();
284  // int nClusterOutsideMustache = matchingPho->nClusterOutsideMustache();
285  //float pfMVA = matchingPho->pfMVA();
286 
287  h_scEta_[0]->Fill(matchingPho->superCluster()->eta());
288  h_scPhi_[0]->Fill(matchingPho->superCluster()->phi());
289 
290  h_phoE_[0][0]->Fill(photonE);
291  h_phoEt_[0][0]->Fill(photonEt);
292 
293  h_phoERes_[0][0]->Fill(photonE / mcEnergy);
294  h_phoSigmaEoE_[0][0]->Fill(sigmaEoE);
295 
296  h_r9_[0][0]->Fill(r9);
297  h_full5x5_r9_[0][0]->Fill(full5x5_r9);
298  h_r1_[0][0]->Fill(r1);
299  h_r2_[0][0]->Fill(r2);
300 
301  h_sigmaIetaIeta_[0][0]->Fill(sieie);
302  h_full5x5_sigmaIetaIeta_[0][0]->Fill(full5x5_sieie);
303  h_hOverE_[0][0]->Fill(hOverE);
304  h_newhOverE_[0][0]->Fill(newhOverE);
305 
306  h_ecalRecHitSumEtConeDR04_[0][0]->Fill(ecalIso);
307  h_hcalTowerSumEtConeDR04_[0][0]->Fill(hcalIso);
308  h_hcalTowerBcSumEtConeDR04_[0][0]->Fill(newhcalIso);
309  h_isoTrkSolidConeDR04_[0][0]->Fill(trkIso);
310  h_nTrkSolidConeDR04_[0][0]->Fill(nIsoTrk);
311 
312  //
313  h_chHadIso_[0]->Fill(chargedHadIso);
314  h_nHadIso_[0]->Fill(neutralHadIso);
315  h_phoIso_[0]->Fill(photonIso);
316 
317  //
318  if (phoIsInBarrel) {
319  h_phoE_[0][1]->Fill(photonE);
320  h_phoEt_[0][1]->Fill(photonEt);
321  h_phoERes_[0][1]->Fill(photonE / mcEnergy);
322  h_phoSigmaEoE_[0][1]->Fill(sigmaEoE);
323 
324  h_r9_[0][1]->Fill(r9);
325  h_full5x5_r9_[0][1]->Fill(full5x5_r9);
326  h_r1_[0][1]->Fill(r1);
327  h_r2_[0][1]->Fill(r2);
328  h_sigmaIetaIeta_[0][1]->Fill(sieie);
329  h_full5x5_sigmaIetaIeta_[0][1]->Fill(full5x5_sieie);
330  h_hOverE_[0][1]->Fill(hOverE);
331  h_newhOverE_[0][1]->Fill(newhOverE);
332  h_ecalRecHitSumEtConeDR04_[0][1]->Fill(ecalIso);
333  h_hcalTowerSumEtConeDR04_[0][1]->Fill(hcalIso);
334  h_hcalTowerBcSumEtConeDR04_[0][1]->Fill(newhcalIso);
335  h_isoTrkSolidConeDR04_[0][1]->Fill(trkIso);
336  h_nTrkSolidConeDR04_[0][1]->Fill(nIsoTrk);
337  h_chHadIso_[1]->Fill(chargedHadIso);
338  h_nHadIso_[1]->Fill(neutralHadIso);
339  h_phoIso_[1]->Fill(photonIso);
340  }
341  if (phoIsInEndcap) {
342  h_phoE_[0][2]->Fill(photonE);
343  h_phoEt_[0][2]->Fill(photonEt);
344 
345  h_phoERes_[0][2]->Fill(photonE / mcEnergy);
346  h_phoSigmaEoE_[0][2]->Fill(sigmaEoE);
347  h_r9_[0][2]->Fill(r9);
348  h_full5x5_r9_[0][2]->Fill(full5x5_r9);
349  h_r1_[0][2]->Fill(r1);
350  h_r2_[0][2]->Fill(r2);
351  h_sigmaIetaIeta_[0][2]->Fill(sieie);
352  h_full5x5_sigmaIetaIeta_[0][2]->Fill(full5x5_sieie);
353  h_hOverE_[0][2]->Fill(hOverE);
354  h_newhOverE_[0][2]->Fill(newhOverE);
355  h_ecalRecHitSumEtConeDR04_[0][2]->Fill(ecalIso);
356  h_hcalTowerSumEtConeDR04_[0][2]->Fill(hcalIso);
357  h_hcalTowerBcSumEtConeDR04_[0][2]->Fill(newhcalIso);
358  h_isoTrkSolidConeDR04_[0][2]->Fill(trkIso);
359  h_nTrkSolidConeDR04_[0][2]->Fill(nIsoTrk);
360  h_chHadIso_[2]->Fill(chargedHadIso);
361  h_nHadIso_[2]->Fill(neutralHadIso);
362  h_phoIso_[2]->Fill(photonIso);
363  }
364 
365  } // end loop over gen photons
366 }
PhotonValidatorMiniAOD::h_newhOverE_
MonitorElement * h_newhOverE_[3][3]
Definition: PhotonValidatorMiniAOD.h:77
muonTagProbeFilters_cff.matched
matched
Definition: muonTagProbeFilters_cff.py:62
electrons_cff.r9
r9
Definition: electrons_cff.py:366
pat::Photon::photonIso
float photonIso() const
Definition: Photon.h:133
PhotonValidatorMiniAOD::h_phoIso_
MonitorElement * h_phoIso_[3]
Definition: PhotonValidatorMiniAOD.h:93
photonAnalyzer_cfi.etMin
etMin
Definition: photonAnalyzer_cfi.py:54
Handle.h
ALCARECOTkAlBeamHalo_cff.etaMin
etaMin
GeV.
Definition: ALCARECOTkAlBeamHalo_cff.py:32
PhotonValidatorMiniAOD::h_scPhi_
MonitorElement * h_scPhi_[2]
Definition: PhotonValidatorMiniAOD.h:68
PhotonValidatorMiniAOD::h_nTrkSolidConeDR04_
MonitorElement * h_nTrkSolidConeDR04_[3][3]
Definition: PhotonValidatorMiniAOD.h:82
MessageLogger.h
genParticles2HepMC_cfi.genParticles
genParticles
Definition: genParticles2HepMC_cfi.py:4
PhotonValidatorMiniAOD::h_phoERes_
MonitorElement * h_phoERes_[3][3]
Definition: PhotonValidatorMiniAOD.h:86
TriggerResults.h
etaBin
int etaBin(const l1t::HGCalMulticluster *cl)
Definition: L1EGammaEEProducer.cc:19
reco::Photon::sigmaIetaIeta
float sigmaIetaIeta() const
Definition: Photon.h:237
sistrip::View
View
Definition: ConstantsForView.h:26
PhotonValidatorMiniAOD::~PhotonValidatorMiniAOD
~PhotonValidatorMiniAOD() override
Definition: PhotonValidatorMiniAOD.cc:44
edm::Run
Definition: Run.h:45
BeamMonitor_cff.phiBin
phiBin
Definition: BeamMonitor_cff.py:75
edm
HLT enums.
Definition: AlignableModifier.h:19
ZMuMuCategoriesSequences_cff.trkIso
trkIso
Definition: ZMuMuCategoriesSequences_cff.py:140
PhotonValidatorMiniAOD::h_hOverE_
MonitorElement * h_hOverE_[3][3]
Definition: PhotonValidatorMiniAOD.h:76
Photon.h
PhotonValidatorMiniAOD::h_phoE_
MonitorElement * h_phoE_[2][3]
Definition: PhotonValidatorMiniAOD.h:84
PackedGenParticle.h
reco::Photon::hadTowOverEm
float hadTowOverEm() const
the ration of hadronic energy in towers behind the BCs in the SC and the SC energy
Definition: Photon.h:217
reco::GenParticleCollection
std::vector< GenParticle > GenParticleCollection
collection of GenParticles
Definition: GenParticleFwd.h:13
PhotonValidatorMiniAOD::h_r9_
MonitorElement * h_r9_[3][3]
Definition: PhotonValidatorMiniAOD.h:70
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
EDAnalyzer.h
PhotonValidatorMiniAOD::h_phoEt_
MonitorElement * h_phoEt_[2][3]
Definition: PhotonValidatorMiniAOD.h:85
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
edm::Handle
Definition: AssociativeIterator.h:50
PhotonValidatorMiniAOD::h_hcalTowerSumEtConeDR04_
MonitorElement * h_hcalTowerSumEtConeDR04_[3][3]
Definition: PhotonValidatorMiniAOD.h:79
PhotonValidatorMiniAOD::photonToken_
edm::EDGetTokenT< edm::View< pat::Photon > > photonToken_
Definition: PhotonValidatorMiniAOD.h:58
GenParticle.h
gedPhotonSequence_cff.photonIso
photonIso
Definition: gedPhotonSequence_cff.py:30
pat::Photon::neutralHadronIso
float neutralHadronIso() const
Definition: Photon.h:130
MakerMacros.h
photonValidator_cfi.resMin
resMin
Definition: photonValidator_cfi.py:91
PhotonValidatorMiniAOD.h
Service.h
PhotonValidatorMiniAOD::h_r2_
MonitorElement * h_r2_[3][3]
Definition: PhotonValidatorMiniAOD.h:75
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
photonAnalyzer_cfi.eMin
eMin
Definition: photonAnalyzer_cfi.py:50
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
AlignmentTrackSelector_cfi.phiMin
phiMin
Definition: AlignmentTrackSelector_cfi.py:18
reco::Photon::nTrkSolidConeDR04
int nTrkSolidConeDR04() const
Definition: Photon.h:428
AlignmentTrackSelector_cfi.phiMax
phiMax
Definition: AlignmentTrackSelector_cfi.py:17
PhotonValidatorMiniAOD::genpartToken_
edm::EDGetTokenT< reco::GenParticleCollection > genpartToken_
Definition: PhotonValidatorMiniAOD.h:59
photonAnalyzer_cfi.r9Min
r9Min
Definition: photonAnalyzer_cfi.py:70
PhotonValidatorMiniAOD::PhotonValidatorMiniAOD
PhotonValidatorMiniAOD(const edm::ParameterSet &)
Definition: PhotonValidatorMiniAOD.cc:33
PhotonValidatorMiniAOD::h_scEta_
MonitorElement * h_scEta_[2]
Definition: PhotonValidatorMiniAOD.h:67
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
Vertex.h
reco::Photon::full5x5_r9
float full5x5_r9() const
Definition: Photon.h:252
photonAnalyzer_cfi.r9Bin
r9Bin
Definition: photonAnalyzer_cfi.py:69
HLT_2018_cff.InputTag
InputTag
Definition: HLT_2018_cff.py:79016
EgHLTOffHistBins_cfi.hOverE
hOverE
Definition: EgHLTOffHistBins_cfi.py:17
edm::ParameterSet
Definition: ParameterSet.h:36
reco::Photon::hcalTowerSumEtBcConeDR04
float hcalTowerSumEtBcConeDR04() const
Hcal isolation sum subtracting the hadronic energy in towers behind the BCs in the SC.
Definition: Photon.h:418
reco::Photon::r9
float r9() const
Definition: Photon.h:240
Event.h
reco::LeafCandidate::eta
double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:152
reco::Photon::getCorrectedEnergyError
float getCorrectedEnergyError(P4type type) const
cosmicPhotonAnalyzer_cfi.eMax
eMax
Definition: cosmicPhotonAnalyzer_cfi.py:10
diffTwoXMLs.r2
r2
Definition: diffTwoXMLs.py:73
Photon
Definition: Photon.py:1
Ptr.h
TriggerNames.h
iEvent
int iEvent
Definition: GenABIO.cc:224
PhotonValidatorMiniAOD::h_full5x5_sigmaIetaIeta_
MonitorElement * h_full5x5_sigmaIetaIeta_[3][3]
Definition: PhotonValidatorMiniAOD.h:73
photonValidator_cfi.resBin
resBin
Definition: photonValidator_cfi.py:90
BPHMonitor_cfi.photons
photons
Definition: BPHMonitor_cfi.py:91
PhotonValidatorMiniAOD::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: PhotonValidatorMiniAOD.cc:46
edm::EventSetup
Definition: EventSetup.h:57
pat
Definition: HeavyIon.h:7
PhotonValidatorMiniAOD::h_full5x5_r9_
MonitorElement * h_full5x5_r9_[3][3]
Definition: PhotonValidatorMiniAOD.h:71
InputTag.h
edm::Ptr
Definition: AssociationVector.h:31
PhotonValidatorMiniAOD::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: PhotonValidatorMiniAOD.cc:196
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
cosmicPhotonAnalyzer_cfi.etMax
etMax
Definition: cosmicPhotonAnalyzer_cfi.py:11
std
Definition: JetResolutionObject.h:76
PhotonValidatorMiniAOD::h_phoSigmaEoE_
MonitorElement * h_phoSigmaEoE_[3][3]
Definition: PhotonValidatorMiniAOD.h:87
writedatasetfile.run
run
Definition: writedatasetfile.py:27
PhotonValidatorMiniAOD::h_nHadIso_
MonitorElement * h_nHadIso_[3]
Definition: PhotonValidatorMiniAOD.h:92
reco::Photon::hadronicOverEm
float hadronicOverEm() const
the total hadronic over electromagnetic fraction
Definition: Photon.h:208
PhotonValidatorMiniAOD::h_r1_
MonitorElement * h_r1_[3][3]
Definition: PhotonValidatorMiniAOD.h:74
reco::LeafCandidate::phi
double phi() const final
momentum azimuthal angle
Definition: LeafCandidate.h:148
PhotonValidatorMiniAOD::h_ecalRecHitSumEtConeDR04_
MonitorElement * h_ecalRecHitSumEtConeDR04_[3][3]
Definition: PhotonValidatorMiniAOD.h:78
Frameworkfwd.h
diffTwoXMLs.r1
r1
Definition: diffTwoXMLs.py:53
reco::Photon::trkSumPtSolidConeDR04
float trkSumPtSolidConeDR04() const
Definition: Photon.h:424
cosmicPhotonAnalyzer_cfi.r9Max
r9Max
Definition: cosmicPhotonAnalyzer_cfi.py:12
reco::Photon::getCandidateP4type
P4type getCandidateP4type() const
Definition: Photon.h:316
ALCARECOTkAlBeamHalo_cff.etaMax
etaMax
Definition: ALCARECOTkAlBeamHalo_cff.py:33
reco::Photon::ecalRecHitSumEtConeDR04
float ecalRecHitSumEtConeDR04() const
Definition: Photon.h:410
reco::Photon::r2x5
float r2x5() const
Definition: Photon.h:239
photonValidator_cfi.resMax
resMax
Definition: photonValidator_cfi.py:92
pat::Photon::chargedHadronIso
float chargedHadronIso() const
Definition: Photon.h:127
PhotonValidatorMiniAOD::parameters_
edm::ParameterSet parameters_
Definition: PhotonValidatorMiniAOD.h:64
dqm::implementation::IBooker
Definition: DQMStore.h:43
reco::LeafCandidate::energy
double energy() const final
energy
Definition: LeafCandidate.h:125
reco::Photon::full5x5_sigmaIetaIeta
float full5x5_sigmaIetaIeta() const
Definition: Photon.h:249
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:30
PhotonValidatorMiniAOD::h_hcalTowerBcSumEtConeDR04_
MonitorElement * h_hcalTowerBcSumEtConeDR04_[3][3]
Definition: PhotonValidatorMiniAOD.h:80
pi
const Double_t pi
Definition: trackSplitPlot.h:36
HGC3DClusterGenMatchSelector_cfi.dR
dR
Definition: HGC3DClusterGenMatchSelector_cfi.py:7
ParameterSet.h
reco::Photon::r1x5
float r1x5() const
Definition: Photon.h:238
PhotonValidatorMiniAOD::h_isoTrkSolidConeDR04_
MonitorElement * h_isoTrkSolidConeDR04_[3][3]
Definition: PhotonValidatorMiniAOD.h:81
photonAnalyzer_cfi.etBin
etBin
Definition: photonAnalyzer_cfi.py:53
edm::Event
Definition: Event.h:73
PhotonValidatorMiniAOD::h_sigmaIetaIeta_
MonitorElement * h_sigmaIetaIeta_[3][3]
Definition: PhotonValidatorMiniAOD.h:72
electrons_cff.sieie
sieie
Definition: electrons_cff.py:367
reco::Photon::hcalTowerSumEtConeDR04
float hcalTowerSumEtConeDR04() const
Hcal isolation sum.
Definition: Photon.h:412
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
PhotonValidatorMiniAOD::h_chHadIso_
MonitorElement * h_chHadIso_[3]
Definition: PhotonValidatorMiniAOD.h:91
photonAnalyzer_cfi.eBin
eBin
CAREFULL for the time being only Loos exists 1 => Tight Photon = GoodCandidate.
Definition: photonAnalyzer_cfi.py:49
pat::Photon::superCluster
reco::SuperClusterRef superCluster() const override
override the superCluster method from CaloJet, to access the internal storage of the supercluster
Generated for CMSSW Reference Manual by doxygen 1.8.16