test
CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EfficiencyPlotter.cc
Go to the documentation of this file.
1 #include "DQMOffline/Muon/interface/EfficiencyPlotter.h"
2 
3 // Framework
4 #include <FWCore/Framework/interface/Event.h>
5 #include "DataFormats/Common/interface/Handle.h"
6 #include <FWCore/Framework/interface/ESHandle.h>
7 #include <FWCore/Framework/interface/MakerMacros.h>
8 #include <FWCore/Framework/interface/EventSetup.h>
9 #include <FWCore/ParameterSet/interface/ParameterSet.h>
10 
11 
12 #include "DQMServices/Core/interface/DQMStore.h"
13 #include "DQMServices/Core/interface/MonitorElement.h"
14 #include "FWCore/MessageLogger/interface/MessageLogger.h"
15 #include "FWCore/Framework/interface/Run.h"
16 
17 #include <iostream>
18 #include <stdio.h>
19 #include <string>
20 #include <math.h>
21 #include "TF1.h"
22 #include "TH1F.h"
23 
24 using namespace edm;
25 using namespace std;
26 
27 //#define DEBUG
28 
29 EfficiencyPlotter::EfficiencyPlotter(const edm::ParameterSet& ps){
30 #ifdef DEBUG
31  cout << "EfficiencyPlotter(): Constructor " << endl;
32 #endif
33  parameters = ps;
34 
35  etaBin = parameters.getParameter<int>("etaBin");
36  etaMin = parameters.getParameter<double>("etaMin");
37  etaMax = parameters.getParameter<double>("etaMax");
38 
39  phiBin = parameters.getParameter<int>("phiBin");
40  phiMin = parameters.getParameter<double>("phiMin");
41  phiMax = parameters.getParameter<double>("phiMax");
42 
43  ptBin = parameters.getParameter<int>("ptBin");
44  ptMin = parameters.getParameter<double>("ptMin");
45  ptMax = parameters.getParameter<double>("ptMax");
46 
47  vtxBin = parameters.getParameter<int>("vtxBin");
48  vtxMin = parameters.getParameter<double>("vtxMin");
49  vtxMax = parameters.getParameter<double>("vtxMax");
50 
51  ID_ = parameters.getParameter<string>("MuonID");
52 }
53 EfficiencyPlotter::~EfficiencyPlotter(){}
54 
55 void EfficiencyPlotter::dqmEndJob(DQMStore::IBooker & ibooker, DQMStore::IGetter & igetter) {
56 
57  ibooker.setCurrentFolder("Muons/EfficiencyAnalyzer");
58 
59  // efficiency plot
60  h_eff_eta_ID = ibooker.book1D("Eff_eta_"+ID_, ID_+" Eff. vs #eta", etaBin, etaMin, etaMax);
61  h_eff_hp_eta_ID = ibooker.book1D("Eff_hp_eta_"+ID_, "High Pt "+ID_+" Eff. vs #eta", etaBin, etaMin, etaMax);
62  h_eff_phi_ID = ibooker.book1D("Eff_phi_"+ID_, ID_+" Eff. vs #phi", phiBin, phiMin, phiMax);
63  h_eff_pt_ID = ibooker.book1D("Eff_pt_"+ID_, ID_+" Eff. vs Pt", ptBin, ptMin, ptMax);
64  h_eff_pt_EB_ID = ibooker.book1D("Eff_pt_EB_"+ID_, "Barrel: "+ID_+" Eff. vs Pt", ptBin, ptMin, ptMax);
65  h_eff_pt_EE_ID = ibooker.book1D("Eff_pt_EE_"+ID_, "Endcap: "+ID_+" Eff. vs Pt", ptBin, ptMin, ptMax);
66  h_eff_pt_detIsoID = ibooker.book1D("Eff_pt_detIso"+ID_, "detIso"+ID_+" Efficiency vs Pt", ptBin, ptMin, ptMax);
67  h_eff_pt_EB_detIsoID = ibooker.book1D("Eff_pt_EB_detIso"+ID_, "Barrel: detIso"+ID_+" Eff. vs Pt", ptBin, ptMin, ptMax);
68  h_eff_pt_EE_detIsoID = ibooker.book1D("Eff_pt_EE_detIso"+ID_, "Endcap: detIso"+ID_+" Eff. vs Pt", ptBin, ptMin, ptMax);
69  h_eff_pt_pfIsoID = ibooker.book1D("Eff_pt_pfIso"+ID_, "pfIso"+ID_+" Eff. vs Pt", ptBin, ptMin, ptMax);
70  h_eff_pt_EB_pfIsoID = ibooker.book1D("Eff_pt_EB_pfIso"+ID_, "Barrel: pfIso"+ID_+" Eff. vs Pt", ptBin, ptMin, ptMax);
71  h_eff_pt_EE_pfIsoID = ibooker.book1D("Eff_pt_EE_pfIso"+ID_, "Endcap: pfIso"+ID_+" Eff. vs Pt", ptBin, ptMin, ptMax);
72  h_eff_vtx_detIsoID = ibooker.book1D("Eff_vtx_detIso"+ID_, "detIso"+ID_+" Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
73  h_eff_vtx_pfIsoID = ibooker.book1D("Eff_vtx_pfIso"+ID_, "pfIso"+ID_+" Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
74  h_eff_vtx_EB_detIsoID = ibooker.book1D("Eff_vtx_EB_detIso"+ID_, "Barrel: detIso"+ID_+" Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
75  h_eff_vtx_EB_pfIsoID = ibooker.book1D("Eff_vtx_EB_pfIso"+ID_, "Barrel: pfIso"+ID_+" Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
76  h_eff_vtx_EE_detIsoID = ibooker.book1D("Eff_vtx_EE_detIso"+ID_, "Endcap: detIso"+ID_+" Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
77  h_eff_vtx_EE_pfIsoID = ibooker.book1D("Eff_vtx_EE_pfIso"+ID_, "Endcap: pfIso"+ID_+" Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
78 
79  h_eff_pt_pfIsodBID = ibooker.book1D("Eff_pt_pfIsodB"+ID_, "pfIso"+ID_+" (deltaBeta) Eff. vs Pt", ptBin, ptMin, ptMax);
80  h_eff_pt_EB_pfIsodBID = ibooker.book1D("Eff_pt_EB_pfIsodB"+ID_, "Barrel: pfIso"+ID_+" (deltaBeta) Eff. vs Pt", ptBin, ptMin, ptMax);
81  h_eff_pt_EE_pfIsodBID = ibooker.book1D("Eff_pt_EE_pfIsodB"+ID_, "Endcap: pfIso"+ID_+" (deltaBeta) Eff. vs Pt", ptBin, ptMin, ptMax);
82  h_eff_vtx_pfIsodBID = ibooker.book1D("Eff_vtx_pfIsodB"+ID_, "pfIso"+ID_+" (deltaBeta) Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
83  h_eff_vtx_EB_pfIsodBID = ibooker.book1D("Eff_vtx_EB_pfIsodB"+ID_, "Barrel: pfIso"+ID_+" (deltaBeta) Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
84  h_eff_vtx_EE_pfIsodBID = ibooker.book1D("Eff_vtx_EE_pfIsodB"+ID_, "Endcap: pfIso"+ID_+" (deltaBeta) Eff. vs nVtx", vtxBin, vtxMin, vtxMax);
85 
86  // This prevents this ME to be normalized when drawn into the GUI
87  h_eff_eta_ID ->setEfficiencyFlag();
88  h_eff_hp_eta_ID ->setEfficiencyFlag();
89  h_eff_phi_ID ->setEfficiencyFlag();
90  h_eff_pt_ID ->setEfficiencyFlag();
91  h_eff_pt_EB_ID ->setEfficiencyFlag();
92  h_eff_pt_EE_ID ->setEfficiencyFlag();
93  h_eff_pt_detIsoID ->setEfficiencyFlag();
94  h_eff_pt_EB_detIsoID ->setEfficiencyFlag();
95  h_eff_pt_EE_detIsoID ->setEfficiencyFlag();
96  h_eff_pt_pfIsoID ->setEfficiencyFlag();
97  h_eff_pt_EB_pfIsoID ->setEfficiencyFlag();
98  h_eff_pt_EE_pfIsoID ->setEfficiencyFlag();
99  h_eff_vtx_detIsoID ->setEfficiencyFlag();
100  h_eff_vtx_pfIsoID ->setEfficiencyFlag();
101  h_eff_vtx_EB_detIsoID->setEfficiencyFlag();
102  h_eff_vtx_EB_pfIsoID ->setEfficiencyFlag();
103  h_eff_vtx_EE_detIsoID->setEfficiencyFlag();
104  h_eff_vtx_EE_pfIsoID ->setEfficiencyFlag();
105 
106  h_eff_pt_pfIsodBID ->setEfficiencyFlag();
107  h_eff_pt_EB_pfIsodBID ->setEfficiencyFlag();
108  h_eff_pt_EE_pfIsodBID ->setEfficiencyFlag();
109  h_eff_vtx_pfIsodBID ->setEfficiencyFlag();
110  h_eff_vtx_EB_pfIsodBID ->setEfficiencyFlag();
111  h_eff_vtx_EE_pfIsodBID ->setEfficiencyFlag();
112 
113 
114  // AXIS TITLES....
115  h_eff_hp_eta_ID ->setAxisTitle("#eta", 1);
116  h_eff_eta_ID ->setAxisTitle("#eta", 1);
117  h_eff_phi_ID ->setAxisTitle("#phi", 1);
118  h_eff_pt_ID ->setAxisTitle("p_{T} (GeV)", 1);
119  h_eff_pt_EB_ID ->setAxisTitle("p_{T} (GeV)", 1);
120  h_eff_pt_EE_ID ->setAxisTitle("p_{T} (GeV)", 1);
121  h_eff_pt_detIsoID ->setAxisTitle("p_{T} (GeV)", 1);
122  h_eff_pt_EB_detIsoID ->setAxisTitle("p_{T} (GeV)", 1);
123  h_eff_pt_EE_detIsoID ->setAxisTitle("p_{T} (GeV)", 1);
124  h_eff_pt_pfIsoID ->setAxisTitle("p_{T} (GeV)", 1);
125  h_eff_pt_EB_pfIsoID ->setAxisTitle("p_{T} (GeV)", 1);
126  h_eff_pt_EE_pfIsoID ->setAxisTitle("p_{T} (GeV)", 1);
127  h_eff_vtx_detIsoID ->setAxisTitle("Number of PV", 1);
128  h_eff_vtx_pfIsoID ->setAxisTitle("Number of PV", 1);
129  h_eff_vtx_EB_detIsoID->setAxisTitle("Number of PV", 1);
130  h_eff_vtx_EB_pfIsoID ->setAxisTitle("Number of PV", 1);
131  h_eff_vtx_EE_detIsoID->setAxisTitle("Number of PV", 1);
132  h_eff_vtx_EE_pfIsoID ->setAxisTitle("Number of PV", 1);
133 
134  h_eff_pt_pfIsodBID ->setAxisTitle("p_{T} (GeV)", 1);
135  h_eff_pt_EB_pfIsodBID ->setAxisTitle("p_{T} (GeV)", 1);
136  h_eff_pt_EE_pfIsodBID ->setAxisTitle("p_{T} (GeV)", 1);
137  h_eff_vtx_pfIsodBID ->setAxisTitle("Number of PV", 1);
138  h_eff_vtx_EB_pfIsodBID ->setAxisTitle("Number of PV", 1);
139  h_eff_vtx_EE_pfIsodBID ->setAxisTitle("Number of PV", 1);
140 
141 
143  string inputdir = "Muons/EfficiencyAnalyzer/"+ID_;
144  string numpath_pt = inputdir+"/passProbes_ID_pt";
145  string denpath_pt = inputdir+"/allProbes_pt";
146 
147  MonitorElement *Numerator_pt = igetter.get(numpath_pt);
148  MonitorElement *Denominator_pt = igetter.get(denpath_pt);
149 
150  if (Numerator_pt && Denominator_pt){
151  TH1F *h_numerator_pt = Numerator_pt->getTH1F();
152  TH1F *h_denominator_pt = Denominator_pt->getTH1F();
153  TH1F *h_eff_pt = h_eff_pt_ID->getTH1F();
154 
155  if (h_eff_pt->GetSumw2N() == 0) h_eff_pt->Sumw2();
156  h_eff_pt->Divide(h_numerator_pt, h_denominator_pt, 1., 1., "B");
157  }
158 
160  string numpath_EB_pt = inputdir+"/passProbes_ID_EB_pt";
161  string denpath_EB_pt = inputdir+"/allProbes_EB_pt";
162 
163  MonitorElement *Numerator_EB_pt = igetter.get(numpath_EB_pt);
164  MonitorElement *Denominator_EB_pt = igetter.get(denpath_EB_pt);
165 
166  if (Numerator_EB_pt && Denominator_EB_pt){
167  TH1F *h_numerator_EB_pt = Numerator_EB_pt->getTH1F();
168  TH1F *h_denominator_EB_pt = Denominator_EB_pt->getTH1F();
169  TH1F *h_eff_EB_pt = h_eff_pt_EB_ID->getTH1F();
170 
171  if (h_eff_EB_pt->GetSumw2N() == 0) h_eff_EB_pt->Sumw2();
172  h_eff_EB_pt->Divide(h_numerator_EB_pt, h_denominator_EB_pt, 1., 1., "B");
173  }
174 
176  string numpath_EE_pt = inputdir+"/passProbes_ID_EE_pt";
177  string denpath_EE_pt = inputdir+"/allProbes_EE_pt";
178 
179  MonitorElement *Numerator_EE_pt = igetter.get(numpath_EE_pt);
180  MonitorElement *Denominator_EE_pt = igetter.get(denpath_EE_pt);
181 
182  if (Numerator_EE_pt && Denominator_EE_pt){
183  TH1F *h_numerator_EE_pt = Numerator_EE_pt->getTH1F();
184  TH1F *h_denominator_EE_pt = Denominator_EE_pt->getTH1F();
185  TH1F *h_eff_EE_pt = h_eff_pt_EE_ID->getTH1F();
186 
187  if (h_eff_EE_pt->GetSumw2N() == 0) h_eff_EE_pt->Sumw2();
188  h_eff_EE_pt->Divide(h_numerator_EE_pt, h_denominator_EE_pt, 1., 1., "B");
189  }
190 
192  string numpath_eta = inputdir+"/passProbes_ID_eta";
193  string denpath_eta = inputdir+"/allProbes_eta";
194 
195  MonitorElement *Numerator_eta = igetter.get(numpath_eta);
196  MonitorElement *Denominator_eta = igetter.get(denpath_eta);
197 
198  if (Numerator_eta && Denominator_eta){
199 
200  TH1F *h_numerator_eta = Numerator_eta->getTH1F();
201  TH1F *h_denominator_eta = Denominator_eta->getTH1F();
202 
203  TH1F *h_eff_eta = h_eff_eta_ID->getTH1F();
204 
205  if (h_eff_eta->GetSumw2N() == 0) h_eff_eta->Sumw2();
206 
207  h_eff_eta->Divide(h_numerator_eta, h_denominator_eta, 1., 1., "B");
208 
209  }
210 
212 
213  string numpath_hp_eta = inputdir+"/passProbes_ID_hp_eta";
214  string denpath_hp_eta = inputdir+"/allProbes_hp_eta";
215 
216  MonitorElement *Numerator_hp_eta = igetter.get(numpath_hp_eta);
217  MonitorElement *Denominator_hp_eta = igetter.get(denpath_hp_eta);
218 
219  if (Numerator_hp_eta && Denominator_hp_eta){
220 
221  TH1F *h_numerator_hp_eta = Numerator_hp_eta->getTH1F();
222  TH1F *h_denominator_hp_eta = Denominator_hp_eta->getTH1F();
223 
224  TH1F *h_eff_hp_eta = h_eff_hp_eta_ID->getTH1F();
225 
226  if (h_eff_hp_eta->GetSumw2N() == 0) h_eff_hp_eta->Sumw2();
227 
228  h_eff_hp_eta->Divide(h_numerator_hp_eta, h_denominator_hp_eta, 1., 1., "B");
229 
230  }
231 
233 
234  string numpath_phi = inputdir+"/passProbes_ID_phi";
235  string denpath_phi = inputdir+"/allProbes_phi";
236 
237  MonitorElement *Numerator_phi = igetter.get(numpath_phi);
238  MonitorElement *Denominator_phi = igetter.get(denpath_phi);
239 
240  if (Numerator_phi && Denominator_phi){
241 
242  TH1F *h_numerator_phi = Numerator_phi->getTH1F();
243  TH1F *h_denominator_phi = Denominator_phi->getTH1F();
244 
245  TH1F *h_eff_phi = h_eff_phi_ID->getTH1F();
246 
247  if (h_eff_phi->GetSumw2N() == 0) h_eff_phi->Sumw2();
248 
249  h_eff_phi->Divide(h_numerator_phi, h_denominator_phi, 1., 1., "B");
250 
251  }
252 
253 
255  string numpath_detIso_pt = inputdir+"/passProbes_detIsoID_pt";
256  string denpath_detIso_pt = inputdir+"/allProbes_ID_pt";
257 
258  MonitorElement *Numerator_detIso_pt = igetter.get(numpath_detIso_pt);
259  MonitorElement *Denominator_detIso_pt = igetter.get(denpath_detIso_pt);
260 
261  if (Numerator_detIso_pt && Denominator_detIso_pt){
262 
263  TH1F *h_numerator_detIso_pt = Numerator_detIso_pt->getTH1F();
264  TH1F *h_denominator_detIso_pt = Denominator_detIso_pt->getTH1F();
265 
266  TH1F *h_eff_detIso_pt = h_eff_pt_detIsoID->getTH1F();
267 
268  if (h_eff_detIso_pt->GetSumw2N() == 0) h_eff_detIso_pt->Sumw2();
269  h_eff_detIso_pt->Divide(h_numerator_detIso_pt, h_denominator_detIso_pt, 1., 1., "B");
270 
271  }
272 
273 
275  string numpath_detIso_EB_pt = inputdir+"/passProbes_EB_detIsoID_pt";
276  string denpath_detIso_EB_pt = inputdir+"/allProbes_EB_ID_pt";
277 
278  MonitorElement *Numerator_detIso_EB_pt = igetter.get(numpath_detIso_EB_pt);
279  MonitorElement *Denominator_detIso_EB_pt = igetter.get(denpath_detIso_EB_pt);
280 
281  if (Numerator_detIso_EB_pt && Denominator_detIso_EB_pt){
282 
283  TH1F *h_numerator_detIso_EB_pt = Numerator_detIso_EB_pt->getTH1F();
284  TH1F *h_denominator_detIso_EB_pt = Denominator_detIso_EB_pt->getTH1F();
285 
286  TH1F *h_eff_detIso_EB_pt = h_eff_pt_EB_detIsoID->getTH1F();
287 
288  if (h_eff_detIso_EB_pt->GetSumw2N() == 0) h_eff_detIso_EB_pt->Sumw2();
289 
290  h_eff_detIso_EB_pt->Divide(h_numerator_detIso_EB_pt, h_denominator_detIso_EB_pt, 1., 1., "B");
291 
292  }
293 
294 
296  string numpath_detIso_EE_pt = inputdir+"/passProbes_EE_detIsoID_pt";
297  string denpath_detIso_EE_pt = inputdir+"/allProbes_EE_ID_pt";
298 
299  MonitorElement *Numerator_detIso_EE_pt = igetter.get(numpath_detIso_EE_pt);
300  MonitorElement *Denominator_detIso_EE_pt = igetter.get(denpath_detIso_EE_pt);
301 
302  if (Numerator_detIso_EE_pt && Denominator_detIso_EE_pt){
303 
304  TH1F *h_numerator_detIso_EE_pt = Numerator_detIso_EE_pt->getTH1F();
305  TH1F *h_denominator_detIso_EE_pt = Denominator_detIso_EE_pt->getTH1F();
306 
307  TH1F *h_eff_detIso_EE_pt = h_eff_pt_EE_detIsoID->getTH1F();
308 
309  if (h_eff_detIso_EE_pt->GetSumw2N() == 0) h_eff_detIso_EE_pt->Sumw2();
310 
311  h_eff_detIso_EE_pt->Divide(h_numerator_detIso_EE_pt, h_denominator_detIso_EE_pt, 1., 1., "B");
312 
313  }
314 
315 
317  string numpath_pfIso_pt = inputdir+"/passProbes_pfIsoID_pt";
318  string denpath_pfIso_pt = inputdir+"/allProbes_ID_pt";
319 
320  MonitorElement *Numerator_pfIso_pt = igetter.get(numpath_pfIso_pt);
321  MonitorElement *Denominator_pfIso_pt = igetter.get(denpath_pfIso_pt);
322 
323  if (Numerator_pfIso_pt && Denominator_pfIso_pt){
324 
325  TH1F *h_numerator_pfIso_pt = Numerator_pfIso_pt->getTH1F();
326  TH1F *h_denominator_pfIso_pt = Denominator_pfIso_pt->getTH1F();
327 
328  TH1F *h_eff_pfIso_pt = h_eff_pt_pfIsoID->getTH1F();
329 
330  if (h_eff_pfIso_pt->GetSumw2N() == 0) h_eff_pfIso_pt->Sumw2();
331 
332  h_eff_pfIso_pt->Divide(h_numerator_pfIso_pt, h_denominator_pfIso_pt, 1., 1., "B");
333  }
334 
335 
337 
338  string numpath_pfIso_EB_pt = inputdir+"/passProbes_EB_pfIsoID_pt";
339  string denpath_pfIso_EB_pt = inputdir+"/allProbes_EB_ID_pt";
340 
341  MonitorElement *Numerator_pfIso_EB_pt = igetter.get(numpath_pfIso_EB_pt);
342  MonitorElement *Denominator_pfIso_EB_pt = igetter.get(denpath_pfIso_EB_pt);
343 
344  if (Numerator_pfIso_EB_pt && Denominator_pfIso_EB_pt){
345 
346  TH1F *h_numerator_pfIso_EB_pt = Numerator_pfIso_EB_pt->getTH1F();
347  TH1F *h_denominator_pfIso_EB_pt = Denominator_pfIso_EB_pt->getTH1F();
348 
349  TH1F *h_eff_pfIso_EB_pt = h_eff_pt_EB_pfIsoID->getTH1F();
350 
351  if (h_eff_pfIso_EB_pt->GetSumw2N() == 0) h_eff_pfIso_EB_pt->Sumw2();
352  h_eff_pfIso_EB_pt->Divide(h_numerator_pfIso_EB_pt, h_denominator_pfIso_EB_pt, 1., 1., "B");
353 
354  }
355 
356 
358  string numpath_pfIso_EE_pt = inputdir+"/passProbes_EE_pfIsoID_pt";
359  string denpath_pfIso_EE_pt = inputdir+"/allProbes_EE_ID_pt";
360 
361  MonitorElement *Numerator_pfIso_EE_pt = igetter.get(numpath_pfIso_EE_pt);
362  MonitorElement *Denominator_pfIso_EE_pt = igetter.get(denpath_pfIso_EE_pt);
363 
364  if (Numerator_pfIso_EE_pt && Denominator_pfIso_EE_pt){
365 
366  TH1F *h_numerator_pfIso_EE_pt = Numerator_pfIso_EE_pt->getTH1F();
367  TH1F *h_denominator_pfIso_EE_pt = Denominator_pfIso_EE_pt->getTH1F();
368 
369  TH1F *h_eff_pfIso_EE_pt = h_eff_pt_EE_pfIsoID->getTH1F();
370 
371  if (h_eff_pfIso_EE_pt->GetSumw2N() == 0) h_eff_pfIso_EE_pt->Sumw2();
372 
373  h_eff_pfIso_EE_pt->Divide(h_numerator_pfIso_EE_pt, h_denominator_pfIso_EE_pt, 1., 1., "B");
374 
375  }
376 
377 
379 
380  string numpath_pfIso_nvtx = inputdir+"/passProbes_pfIsoID_nVtx";
381  string denpath_pfIso_nvtx = inputdir+"/allProbes_ID_nVtx";
382 
383  MonitorElement *Numerator_pfIso_nvtx = igetter.get(numpath_pfIso_nvtx);
384  MonitorElement *Denominator_pfIso_nvtx = igetter.get(denpath_pfIso_nvtx);
385 
386  if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
387 
388  TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->getTH1F();
389  TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->getTH1F();
390 
391  TH1F *h_eff_pfIso_nvtx = h_eff_vtx_pfIsoID->getTH1F();
392 
393  if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
394 
395  h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1., "B");
396 
397  }
398 
399 
400 
402  string numpath_detIso_nvtx = inputdir+"/passProbes_detIsoID_nVtx";
403  string denpath_detIso_nvtx = inputdir+"/allProbes_ID_nVtx";
404 
405  MonitorElement *Numerator_detIso_nvtx = igetter.get(numpath_detIso_nvtx);
406  MonitorElement *Denominator_detIso_nvtx = igetter.get(denpath_detIso_nvtx);
407 
408  if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
409 
410  TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->getTH1F();
411  TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->getTH1F();
412 
413  TH1F *h_eff_detIso_nvtx = h_eff_vtx_detIsoID->getTH1F();
414 
415  if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
416 
417  h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1., "B");
418 
419  }
420 
421 
423  numpath_detIso_nvtx = inputdir+"/passProbes_EB_detIsoID_nVtx";
424  denpath_detIso_nvtx = inputdir+"/allProbes_EB_ID_nVtx";
425 
426  Numerator_detIso_nvtx = igetter.get(numpath_detIso_nvtx);
427  Denominator_detIso_nvtx = igetter.get(denpath_detIso_nvtx);
428 
429  if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
430  TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->getTH1F();
431  TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->getTH1F();
432 
433  TH1F *h_eff_detIso_nvtx = h_eff_vtx_EB_detIsoID->getTH1F();
434 
435  if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
436  h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1., "B");
437  }
438 
440  numpath_detIso_nvtx = inputdir+"/passProbes_EE_detIsoID_nVtx";
441  denpath_detIso_nvtx = inputdir+"/allProbes_EE_ID_nVtx";
442 
443  Numerator_detIso_nvtx = igetter.get(numpath_detIso_nvtx);
444  Denominator_detIso_nvtx = igetter.get(denpath_detIso_nvtx);
445 
446  if (Numerator_detIso_nvtx && Denominator_detIso_nvtx){
447 
448  TH1F *h_numerator_detIso_nvtx = Numerator_detIso_nvtx->getTH1F();
449  TH1F *h_denominator_detIso_nvtx = Denominator_detIso_nvtx->getTH1F();
450 
451  TH1F *h_eff_detIso_nvtx = h_eff_vtx_EE_detIsoID->getTH1F();
452 
453  if (h_eff_detIso_nvtx->GetSumw2N() == 0) h_eff_detIso_nvtx->Sumw2();
454 
455  h_eff_detIso_nvtx->Divide(h_numerator_detIso_nvtx, h_denominator_detIso_nvtx, 1., 1., "B");
456 
457  }
458 
459 
461  numpath_pfIso_nvtx = inputdir+"/passProbes_EB_pfIsoID_nVtx";
462  denpath_pfIso_nvtx = inputdir+"/allProbes_EB_ID_nVtx";
463 
464  Numerator_pfIso_nvtx = igetter.get(numpath_pfIso_nvtx);
465  Denominator_pfIso_nvtx = igetter.get(denpath_pfIso_nvtx);
466 
467  if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
468 
469  TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->getTH1F();
470  TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->getTH1F();
471 
472  TH1F *h_eff_pfIso_nvtx = h_eff_vtx_EB_pfIsoID->getTH1F();
473 
474  if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
475 
476  h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1., "B");
477  }
478 
480  numpath_pfIso_nvtx = inputdir+"/passProbes_EE_pfIsoID_nVtx";
481  denpath_pfIso_nvtx = inputdir+"/allProbes_EE_ID_nVtx";
482 
483  Numerator_pfIso_nvtx = igetter.get(numpath_pfIso_nvtx);
484  Denominator_pfIso_nvtx = igetter.get(denpath_pfIso_nvtx);
485 
486  if (Numerator_pfIso_nvtx && Denominator_pfIso_nvtx){
487 
488  TH1F *h_numerator_pfIso_nvtx = Numerator_pfIso_nvtx->getTH1F();
489  TH1F *h_denominator_pfIso_nvtx = Denominator_pfIso_nvtx->getTH1F();
490 
491  TH1F *h_eff_pfIso_nvtx = h_eff_vtx_EE_pfIsoID->getTH1F();
492 
493  if (h_eff_pfIso_nvtx->GetSumw2N() == 0) h_eff_pfIso_nvtx->Sumw2();
494 
495  h_eff_pfIso_nvtx->Divide(h_numerator_pfIso_nvtx, h_denominator_pfIso_nvtx, 1., 1., "B");
496 
497  }
498 
499 
500 
502  string numpath_pfIsodB_pt = inputdir+"/passProbes_pfIsodBID_pt";
503  string denpath_pfIsodB_pt = inputdir+"/allProbes_ID_pt";
504 
505  MonitorElement *Numerator_pfIsodB_pt = igetter.get(numpath_pfIsodB_pt);
506  MonitorElement *Denominator_pfIsodB_pt = igetter.get(denpath_pfIsodB_pt);
507 
508  if (Numerator_pfIsodB_pt && Denominator_pfIsodB_pt){
509 
510  TH1F *h_numerator_pfIsodB_pt = Numerator_pfIsodB_pt->getTH1F();
511  TH1F *h_denominator_pfIsodB_pt = Denominator_pfIsodB_pt->getTH1F();
512 
513  TH1F *h_eff_pfIsodB_pt = h_eff_pt_pfIsodBID->getTH1F();
514 
515  if (h_eff_pfIsodB_pt->GetSumw2N() == 0) h_eff_pfIsodB_pt->Sumw2();
516 
517  h_eff_pfIsodB_pt->Divide(h_numerator_pfIsodB_pt, h_denominator_pfIsodB_pt, 1., 1., "B");
518  }
519 
520 
522 
523  string numpath_pfIsodB_EB_pt = inputdir+"/passProbes_EB_pfIsodBID_pt";
524  string denpath_pfIsodB_EB_pt = inputdir+"/allProbes_EB_ID_pt";
525 
526  MonitorElement *Numerator_pfIsodB_EB_pt = igetter.get(numpath_pfIsodB_EB_pt);
527  MonitorElement *Denominator_pfIsodB_EB_pt = igetter.get(denpath_pfIsodB_EB_pt);
528 
529  if (Numerator_pfIsodB_EB_pt && Denominator_pfIsodB_EB_pt){
530 
531  TH1F *h_numerator_pfIsodB_EB_pt = Numerator_pfIsodB_EB_pt->getTH1F();
532  TH1F *h_denominator_pfIsodB_EB_pt = Denominator_pfIsodB_EB_pt->getTH1F();
533 
534  TH1F *h_eff_pfIsodB_EB_pt = h_eff_pt_EB_pfIsodBID->getTH1F();
535 
536  if (h_eff_pfIsodB_EB_pt->GetSumw2N() == 0) h_eff_pfIsodB_EB_pt->Sumw2();
537  h_eff_pfIsodB_EB_pt->Divide(h_numerator_pfIsodB_EB_pt, h_denominator_pfIsodB_EB_pt, 1., 1., "B");
538 
539  }
540 
541 
543  string numpath_pfIsodB_EE_pt = inputdir+"/passProbes_EE_pfIsodBID_pt";
544  string denpath_pfIsodB_EE_pt = inputdir+"/allProbes_EE_ID_pt";
545 
546  MonitorElement *Numerator_pfIsodB_EE_pt = igetter.get(numpath_pfIsodB_EE_pt);
547  MonitorElement *Denominator_pfIsodB_EE_pt = igetter.get(denpath_pfIsodB_EE_pt);
548 
549  if (Numerator_pfIsodB_EE_pt && Denominator_pfIsodB_EE_pt){
550 
551  TH1F *h_numerator_pfIsodB_EE_pt = Numerator_pfIsodB_EE_pt->getTH1F();
552  TH1F *h_denominator_pfIsodB_EE_pt = Denominator_pfIsodB_EE_pt->getTH1F();
553 
554  TH1F *h_eff_pfIsodB_EE_pt = h_eff_pt_EE_pfIsodBID->getTH1F();
555 
556  if (h_eff_pfIsodB_EE_pt->GetSumw2N() == 0) h_eff_pfIsodB_EE_pt->Sumw2();
557 
558  h_eff_pfIsodB_EE_pt->Divide(h_numerator_pfIsodB_EE_pt, h_denominator_pfIsodB_EE_pt, 1., 1., "B");
559 
560  }
561 
562 
564 
565  string numpath_pfIsodB_nvtx = inputdir+"/passProbes_pfIsodBID_nVtx";
566  string denpath_pfIsodB_nvtx = inputdir+"/allProbes_ID_nVtx";
567 
568  MonitorElement *Numerator_pfIsodB_nvtx = igetter.get(numpath_pfIsodB_nvtx);
569  MonitorElement *Denominator_pfIsodB_nvtx = igetter.get(denpath_pfIsodB_nvtx);
570 
571  if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx){
572 
573  TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->getTH1F();
574  TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->getTH1F();
575 
576  TH1F *h_eff_pfIsodB_nvtx = h_eff_vtx_pfIsodBID->getTH1F();
577 
578  if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0) h_eff_pfIsodB_nvtx->Sumw2();
579 
580  h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1., "B");
581 
582  }
583 
584 
586  numpath_pfIsodB_nvtx = inputdir+"/passProbes_EB_pfIsodBID_nVtx";
587  denpath_pfIsodB_nvtx = inputdir+"/allProbes_EB_ID_nVtx";
588 
589  Numerator_pfIsodB_nvtx = igetter.get(numpath_pfIsodB_nvtx);
590  Denominator_pfIsodB_nvtx = igetter.get(denpath_pfIsodB_nvtx);
591 
592  if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx){
593 
594  TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->getTH1F();
595  TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->getTH1F();
596 
597  TH1F *h_eff_pfIsodB_nvtx = h_eff_vtx_EB_pfIsodBID->getTH1F();
598 
599  if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0) h_eff_pfIsodB_nvtx->Sumw2();
600 
601  h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1., "B");
602  }
603 
604 
606  numpath_pfIsodB_nvtx = inputdir+"/passProbes_EE_pfIsodBID_nVtx";
607  denpath_pfIsodB_nvtx = inputdir+"/allProbes_EE_ID_nVtx";
608 
609  Numerator_pfIsodB_nvtx = igetter.get(numpath_pfIsodB_nvtx);
610  Denominator_pfIsodB_nvtx = igetter.get(denpath_pfIsodB_nvtx);
611 
612  if (Numerator_pfIsodB_nvtx && Denominator_pfIsodB_nvtx){
613 
614  TH1F *h_numerator_pfIsodB_nvtx = Numerator_pfIsodB_nvtx->getTH1F();
615  TH1F *h_denominator_pfIsodB_nvtx = Denominator_pfIsodB_nvtx->getTH1F();
616 
617  TH1F *h_eff_pfIsodB_nvtx = h_eff_vtx_EE_pfIsodBID->getTH1F();
618 
619  if (h_eff_pfIsodB_nvtx->GetSumw2N() == 0) h_eff_pfIsodB_nvtx->Sumw2();
620 
621  h_eff_pfIsodB_nvtx->Divide(h_numerator_pfIsodB_nvtx, h_denominator_pfIsodB_nvtx, 1., 1., "B");
622 
623  }
624 
625 
626 }
627 
DQMStore::IGetter::get
MonitorElement * get(const std::string &path)
Definition: DQMStore.cc:304
EfficiencyPlotter::EfficiencyPlotter
EfficiencyPlotter(const edm::ParameterSet &ps)
Constructor.
Definition: EfficiencyPlotter.cc:29
EfficiencyPlotter::dqmEndJob
void dqmEndJob(DQMStore::IBooker &, DQMStore::IGetter &) override
DQM Client Diagnostic.
Definition: EfficiencyPlotter.cc:55
diMuonHistograms_cfi.etaBin
int etaBin
Definition: diMuonHistograms_cfi.py:31
EfficiencyPlotter::~EfficiencyPlotter
virtual ~EfficiencyPlotter()
Destructor.
Definition: EfficiencyPlotter.cc:53
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
HLT_25ns10e33_v2_cff.parameters
tuple parameters
Definition: HLT_25ns10e33_v2_cff.py:3864
PtMinSelector_cfg.ptMin
tuple ptMin
Definition: PtMinSelector_cfg.py:44
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:276
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1457
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
MonitorElement::setEfficiencyFlag
void setEfficiencyFlag(void)
Definition: MonitorElement.h:159
HLT_25ns10e33_v2_cff.etaMax
tuple etaMax
Definition: HLT_25ns10e33_v2_cff.py:8668
HLT_25ns10e33_v2_cff.etaMin
tuple etaMin
Definition: HLT_25ns10e33_v2_cff.py:8667
Generated for CMSSW Reference Manual by doxygen 1.8.5