CMS 3D CMS Logo
ElectronMcFakePostValidator.cc
Go to the documentation of this file.
1 
2 #include "Validation/RecoEgamma/plugins/ElectronMcFakePostValidator.h"
3 #include "DQMServices/Core/interface/DQMStore.h"
4 #include "FWCore/MessageLogger/interface/MessageLogger.h"
5 
6 ElectronMcFakePostValidator::ElectronMcFakePostValidator(const edm::ParameterSet& conf)
7  : ElectronDqmHarvesterBase(conf) {
8  // histos bining and limits
9 
10  edm::ParameterSet histosSet = conf.getParameter<edm::ParameterSet>("histosCfg");
11 
12  set_EfficiencyFlag = histosSet.getParameter<bool>("EfficiencyFlag");
13  set_StatOverflowFlag = histosSet.getParameter<bool>("StatOverflowFlag");
14 }
15 
16 ElectronMcFakePostValidator::~ElectronMcFakePostValidator() {}
17 
18 void ElectronMcFakePostValidator::finalize(DQMStore::IBooker& iBooker, DQMStore::IGetter& iGetter) {
19  setBookIndex(-1);
20  setBookPrefix("h_ele");
21  setBookEfficiencyFlag(set_EfficiencyFlag);
22  setBookStatOverflowFlag(set_StatOverflowFlag);
23 
24  edm::LogInfo("ElectronMcFakePostValidator::finalize") << "efficiency calculation ";
25  bookH1andDivide(
26  iBooker, iGetter, "etaEff", "matchingObjectEta_matched", "matchingObject_eta", "#eta", "Efficiency", "");
27  bookH1andDivide(iBooker, iGetter, "zEff", "matchingObjectZ_matched", "matchingObject_z", "z (cm)", "Efficiency", "");
28  bookH1andDivide(iBooker,
29  iGetter,
30  "absetaEff",
31  "matchingObjectAbsEta_matched",
32  "matchingObject_abseta",
33  "|#eta|",
34  "Efficiency",
35  "");
36  bookH1andDivide(
37  iBooker, iGetter, "ptEff", "matchingObjectPt_matched", "matchingObject_Pt", "p_{T} (GeV/c)", "Efficiency", "");
38  bookH1andDivide(
39  iBooker, iGetter, "phiEff", "matchingObjectPhi_matched", "matchingObject_phi", "#phi (rad)", "Efficiency", "");
40  // bookH2andDivide(iBooker,iGetter, "ptEtaEff","matchingObjectPtEta_matched","matchingObjectPtEta","#eta","p_{T} (GeV/c)","");
41  //
42  // bookH1andDivide(iBooker,iGetter, "etaQmisid","matchingObjectEta_matched_qmisid","h_simEta","#eta","q misId","");
43  // bookH1andDivide(iBooker,iGetter, "zQmisid","matchingObjectZ_matched_qmisid","h_simZ","z (cm)","q misId","");
44  // bookH1andDivide(iBooker,iGetter, "absetaQmisid","matchingObjectAbsEta_matched_qmisid","h_simAbsEta","|#eta|","q misId","");
45  // bookH1andDivide(iBooker,iGetter, "ptQmisid","matchingObjectPt_matched_qmisid","h_simPt","p_{T} (GeV/c)","q misId","");
46 
47  edm::LogInfo("ElectronMcFakePostValidator::finalize") << "all reco electrons ";
48  bookH1andDivide(iBooker, iGetter, "etaEff_all", "vertexEta_all", "matchingObject_eta", "#eta", "N_{rec}/N_{gen}", "");
49  bookH1andDivide(
50  iBooker, iGetter, "ptEff_all", "vertexPt_all", "matchingObject_Pt", "p_{T} (GeV/c)", "N_{rec}/N_{gen}", "");
51 
52  edm::LogInfo("ElectronMcFakePostValidator::finalize") << "classes";
53  bookH1andDivide(iBooker,
54  iGetter,
55  "eta_goldenFrac",
56  "eta_golden",
57  "h_ele_eta",
58  "|#eta|",
59  "Fraction of electrons",
60  "fraction of golden electrons vs eta");
61  bookH1andDivide(iBooker,
62  iGetter,
63  "eta_bbremFrac",
64  "eta_bbrem",
65  "h_ele_eta",
66  "|#eta|",
67  "Fraction of electrons",
68  "fraction of big brem electrons vs eta");
69  // bookH1andDivide(iBooker,iGetter, "eta_narrowFrac","eta_narrow","h_ele_eta","|#eta|","Fraction of electrons","fraction of narrow electrons vs eta");
70  bookH1andDivide(iBooker,
71  iGetter,
72  "eta_showerFrac",
73  "eta_shower",
74  "h_ele_eta",
75  "|#eta|",
76  "Fraction of electrons",
77  "fraction of showering electrons vs eta");
78 
79  // fbrem
80  MonitorElement* p1_ele_fbremVsEta_mean = get(iGetter, "fbremvsEtamean");
81  TAxis* etaAxis = p1_ele_fbremVsEta_mean->getTProfile()->GetXaxis();
82  MonitorElement* h1_ele_xOverX0VsEta = bookH1withSumw2(
83  iBooker, "xOverx0VsEta", "mean X/X_0 vs eta", etaAxis->GetNbins(), etaAxis->GetXmin(), etaAxis->GetXmax());
84  for (int ibin = 1; ibin < etaAxis->GetNbins() + 1; ibin++) {
85  double xOverX0 = 0.;
86  if (p1_ele_fbremVsEta_mean->getBinContent(ibin) > 0.) {
87  xOverX0 = -log(p1_ele_fbremVsEta_mean->getBinContent(ibin));
88  }
89  h1_ele_xOverX0VsEta->setBinContent(ibin, xOverX0);
90  }
91 
92 
93  MonitorElement* h1_ele_provenance = get(iGetter, "provenance");
94  if (h1_ele_provenance->getBinContent(3) > 0) {
95  h1_ele_provenance->getTH1F()->Scale(1. / h1_ele_provenance->getBinContent(3));
96  }
97  MonitorElement* h1_ele_provenance_barrel = get(iGetter, "provenance_barrel");
98  if (h1_ele_provenance_barrel->getBinContent(3) > 0) {
99  h1_ele_provenance_barrel->getTH1F()->Scale(1. / h1_ele_provenance_barrel->getBinContent(3));
100  }
101  MonitorElement* h1_ele_provenance_endcaps = get(iGetter, "provenance_endcaps");
102  if (h1_ele_provenance_endcaps->getBinContent(3) > 0) {
103  h1_ele_provenance_endcaps->getTH1F()->Scale(1. / h1_ele_provenance_endcaps->getBinContent(3));
104  }
105 
106  // profiles from 2D histos
107  profileX(iBooker, iGetter, "PoPmatchingObjectVsEta", "", "#eta", "<P/P_{gen}>");
108  profileX(iBooker, iGetter, "PoPmatchingObjectVsPhi", "", "#phi (rad)", "<P/P_{gen}>");
109  profileX(iBooker, iGetter, "EtaMnEtamatchingObjectVsEta", "", "#eta", "<#eta_{rec} - #eta_{gen}>");
110  profileX(iBooker, iGetter, "EtaMnEtamatchingObjectVsPhi", "", "#phi (rad)", "<#eta_{rec} - #eta_{gen}>");
111  profileX(iBooker, iGetter, "PhiMnPhimatchingObjectVsEta", "", "#eta", "<#phi_{rec} - #phi_{gen}> (rad)");
112  profileX(iBooker, iGetter, "PhiMnPhimatchingObjectVsPhi", "", "#phi (rad)", "");
113  profileX(iBooker, iGetter, "vertexPtVsEta", "mean ele transverse momentum vs eta", "#eta", "<p_{T}> (GeV/c)");
114  profileX(iBooker, iGetter, "vertexPtVsPhi", "mean ele transverse momentum vs phi", "#phi (rad)", "<p_{T}> (GeV/c)");
115  profileX(iBooker, iGetter, "EoPVsEta", "mean ele E/p vs eta", "#eta", "<E/P_{vertex}>");
116  profileX(iBooker, iGetter, "EoPVsPhi", "mean ele E/p vs phi", "#phi (rad)", "<E/P_{vertex}>");
117  profileX(iBooker, iGetter, "EoPoutVsEta", "mean ele E/pout vs eta", "#eta", "<E_{seed}/P_{out}>");
118  profileX(iBooker, iGetter, "EoPoutVsPhi", "mean ele E/pout vs phi", "#phi (rad)", "<E_{seed}/P_{out}>");
119  profileX(iBooker, iGetter, "EeleOPoutVsEta", "mean ele Eele/pout vs eta", "#eta", "<E_{ele}/P_{out}>");
120  profileX(iBooker, iGetter, "EeleOPoutVsPhi", "mean ele Eele/pout vs phi", "#phi (rad)", "<E_{ele}/P_{out}>");
121  profileX(iBooker, iGetter, "HoEVsEta", "mean ele H/E vs eta", "#eta", "<H/E>");
122  profileX(iBooker, iGetter, "HoEVsPhi", "mean ele H/E vs phi", "#phi (rad)", "<H/E>");
123  profileX(iBooker, iGetter, "chi2VsEta", "mean ele track chi2 vs eta", "#eta", "<#Chi^{2}>");
124  profileX(iBooker, iGetter, "chi2VsPhi", "mean ele track chi2 vs phi", "#phi (rad)", "<#Chi^{2}>");
125  profileX(iBooker, iGetter, "foundHitsVsEta", "mean ele track # found hits vs eta", "#eta", "<N_{hits}>");
126  profileX(iBooker, iGetter, "foundHitsVsPhi", "mean ele track # found hits vs phi", "#phi (rad)", "<N_{hits}>");
127  profileX(iBooker, iGetter, "lostHitsVsEta", "mean ele track # lost hits vs eta", "#eta", "<N_{hits}>");
128  profileX(iBooker, iGetter, "lostHitsVsPhi", "mean ele track # lost hits vs phi", "#phi (rad)", "<N_{hits}>");
129  profileX(iBooker, iGetter, "vertexTIPVsEta", "mean tip (wrt gen vtx) vs eta", "#eta", "<TIP> (cm)");
130  profileX(iBooker, iGetter, "vertexTIPVsPhi", "mean tip (wrt gen vtx) vs phi", "#phi", "<TIP> (cm)");
131  profileX(iBooker, iGetter, "vertexTIPVsPt", "mean tip (wrt gen vtx) vs phi", "p_{T} (GeV/c)", "<TIP> (cm)");
132  profileX(iBooker,
133  iGetter,
134  "seedDphi2_VsEta",
135  "mean ele seed dphi 2nd layer vs eta",
136  "#eta",
137  "<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)",
138  -0.004,
139  0.004);
140  profileX(iBooker,
141  iGetter,
142  "seedDphi2_VsPt",
143  "mean ele seed dphi 2nd layer vs pt",
144  "p_{T} (GeV/c)",
145  "<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)",
146  -0.004,
147  0.004);
148  profileX(iBooker,
149  iGetter,
150  "seedDrz2_VsEta",
151  "mean ele seed dr(dz) 2nd layer vs eta",
152  "#eta",
153  "<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)",
154  -0.15,
155  0.15);
156  profileX(iBooker,
157  iGetter,
158  "seedDrz2_VsPt",
159  "mean ele seed dr(dz) 2nd layer vs pt",
160  "p_{T} (GeV/c)",
161  "<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)",
162  -0.15,
163  0.15);
164  profileX(iBooker,
165  iGetter,
166  "seedDphi2Pos_VsEta",
167  "mean ele seed dphi 2nd layer positron vs eta",
168  "#eta",
169  "<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)",
170  -0.004,
171  0.004);
172  profileX(iBooker,
173  iGetter,
174  "seedDphi2Pos_VsPt",
175  "mean ele seed dphi 2nd layer positron vs pt",
176  "p_{T} (GeV/c)",
177  "<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)",
178  -0.004,
179  0.004);
180  profileX(iBooker,
181  iGetter,
182  "seedDrz2Pos_VsEta",
183  "mean ele seed dr(dz) 2nd layer positron vs eta",
184  "#eta",
185  "<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)",
186  -0.15,
187  0.15);
188  profileX(iBooker,
189  iGetter,
190  "seedDrz2Pos_VsPt",
191  "mean ele seed dr(dz) 2nd layer positron vs pt",
192  "p_{T} (GeV/c)",
193  "<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)",
194  -0.15,
195  0.15);
196 
197 }
ElectronDqmHarvesterBase::setBookEfficiencyFlag
void setBookEfficiencyFlag(const bool &)
Definition: ElectronDqmHarvesterBase.cc:34
ElectronMcFakePostValidator::ElectronMcFakePostValidator
ElectronMcFakePostValidator(const edm::ParameterSet &conf)
Definition: ElectronMcFakePostValidator.cc:6
ElectronDqmHarvesterBase::setBookPrefix
void setBookPrefix(const std::string &)
Definition: ElectronDqmHarvesterBase.cc:30
MessageLogger.h
ElectronDqmHarvesterBase::bookH1andDivide
MonitorElement * bookH1andDivide(DQMStore::IBooker &iBooker, DQMStore::IGetter &, const std::string &name, MonitorElement *num, MonitorElement *denom, const std::string &titleX, const std::string &titleY, const std::string &title="")
Definition: ElectronDqmHarvesterBase.cc:326
ElectronDqmHarvesterBase::bookH1withSumw2
MonitorElement * bookH1withSumw2(DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")
Definition: ElectronDqmHarvesterBase.cc:212
DQMStore.h
dqm::legacy::MonitorElement
Definition: MonitorElement.h:462
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
dqm::legacy::MonitorElement::getTH1F
virtual TH1F * getTH1F() const
Definition: MonitorElement.h:479
ElectronDqmHarvesterBase
Definition: ElectronDqmHarvesterBase.h:16
ElectronMcFakePostValidator::~ElectronMcFakePostValidator
~ElectronMcFakePostValidator() override
Definition: ElectronMcFakePostValidator.cc:16
ElectronMcFakePostValidator::set_StatOverflowFlag
bool set_StatOverflowFlag
Definition: ElectronMcFakePostValidator.h:22
ElectronMcFakePostValidator.h
ElectronDqmHarvesterBase::profileX
MonitorElement * profileX(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter, MonitorElement *me2d, const std::string &title="", const std::string &titleX="", const std::string &titleY="", Double_t minimum=-1111, Double_t maximum=-1111)
Definition: ElectronDqmHarvesterBase.cc:407
dqm::legacy::MonitorElement::getTProfile
virtual TProfile * getTProfile() const
Definition: MonitorElement.h:507
ElectronMcFakePostValidator::finalize
void finalize(DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter) override
Definition: ElectronMcFakePostValidator.cc:18
edm::ParameterSet
Definition: ParameterSet.h:47
ElectronMcFakePostValidator::h1_ele_xOverX0VsEta
MonitorElement * h1_ele_xOverX0VsEta
Definition: ElectronMcFakePostValidator.h:25
dqm::impl::MonitorElement::setBinContent
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:691
ElectronDqmHarvesterBase::get
MonitorElement * get(DQMStore::IGetter &iGetter, const std::string &name)
Definition: ElectronDqmHarvesterBase.cc:120
ElectronDqmHarvesterBase::setBookStatOverflowFlag
void setBookStatOverflowFlag(const bool &)
Definition: ElectronDqmHarvesterBase.cc:36
dqm::implementation::IGetter
Definition: DQMStore.h:484
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
ElectronDqmHarvesterBase::setBookIndex
void setBookIndex(short)
Definition: ElectronDqmHarvesterBase.cc:32
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
dqm::implementation::IBooker
Definition: DQMStore.h:43
dqm::impl::MonitorElement::getBinContent
virtual double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:592
ElectronMcFakePostValidator::set_EfficiencyFlag
bool set_EfficiencyFlag
Definition: ElectronMcFakePostValidator.h:21
Generated for CMSSW Reference Manual by doxygen 1.8.16