40 #include "CLHEP/Units/GlobalPhysicalConstants.h"
55 outputFile_ = conf.
getParameter<std::string>(
"outputFile");
56 histfile_ =
new TFile(outputFile_.c_str(),
"RECREATE");
68 etamax=pset.getParameter<
double>(
"Etamax");
69 phimin=pset.getParameter<
double>(
"Phimin");
70 phimax=pset.getParameter<
double>(
"Phimax");
71 ptmax=pset.getParameter<
double>(
"Ptmax");
72 pmax=pset.getParameter<
double>(
"Pmax");
73 eopmax=pset.getParameter<
double>(
"Eopmax");
74 eopmaxsht=pset.getParameter<
double>(
"Eopmaxsht");
75 detamin=pset.getParameter<
double>(
"Detamin");
76 detamax=pset.getParameter<
double>(
"Detamax");
77 dphimin=pset.getParameter<
double>(
"Dphimin");
78 dphimax=pset.getParameter<
double>(
"Dphimax");
79 detamatchmin=pset.getParameter<
double>(
"Detamatchmin");
80 detamatchmax=pset.getParameter<
double>(
"Detamatchmax");
81 dphimatchmin=pset.getParameter<
double>(
"Dphimatchmin");
82 dphimatchmax=pset.getParameter<
double>(
"Dphimatchmax");
83 fhitsmax=pset.getParameter<
double>(
"Fhitsmax");
84 lhitsmax=pset.getParameter<
double>(
"Lhitsmax");
85 nbineta=pset.getParameter<
int>(
"Nbineta");
86 nbineta2D=pset.getParameter<
int>(
"Nbineta2D");
87 nbinp=pset.getParameter<
int>(
"Nbinp");
88 nbinpt=pset.getParameter<
int>(
"Nbinpt");
89 nbinp2D=pset.getParameter<
int>(
"Nbinp2D");
90 nbinpt2D=pset.getParameter<
int>(
"Nbinpt2D");
91 nbinpteff=pset.getParameter<
int>(
"Nbinpteff");
92 nbinphi=pset.getParameter<
int>(
"Nbinphi");
93 nbinphi2D=pset.getParameter<
int>(
"Nbinphi2D");
94 nbineop=pset.getParameter<
int>(
"Nbineop");
95 nbineop2D=pset.getParameter<
int>(
"Nbineop2D");
96 nbinfhits=pset.getParameter<
int>(
"Nbinfhits");
97 nbinlhits=pset.getParameter<
int>(
"Nbinlhits");
98 nbinxyz=pset.getParameter<
int>(
"Nbinxyz");
99 nbindeta=pset.getParameter<
int>(
"Nbindeta");
100 nbindphi=pset.getParameter<
int>(
"Nbindphi");
101 nbindetamatch=pset.getParameter<
int>(
"Nbindetamatch");
102 nbindphimatch=pset.getParameter<
int>(
"Nbindphimatch");
103 nbindetamatch2D=pset.getParameter<
int>(
"Nbindetamatch2D");
104 nbindphimatch2D=pset.getParameter<
int>(
"Nbindphimatch2D");
105 nbinmee= pset.getParameter<
int>(
"Nbinmee");
106 meemin=pset.getParameter<
double>(
"Meemin");
107 meemax=pset.getParameter<
double>(
"Meemax");
108 nbinhoe= pset.getParameter<
int>(
"Nbinhoe");
109 hoemin=pset.getParameter<
double>(
"Hoemin");
110 hoemax=pset.getParameter<
double>(
"Hoemax");
131 if ( locJet != std::string::npos ) {
132 std::cout <<
"Matching objects are GenJets " << std::endl;
135 std::cout <<
"Didn't recognize input matching objects!! " << std::endl;
138 std::string htitle, hlabel;
139 hlabel=
"h_"+type_+
"Num"; htitle=
"# "+type_+
"s";
140 h_matchingObjectNum =
new TH1F( hlabel.c_str(), htitle.c_str(), nbinfhits,0.,fhitsmax );
143 histNum_=
new TH1F(
"h_recEleNum",
"# rec electrons",20, 0.,20.);
146 hlabel=
"h_matchingObject_eta"; htitle=type_+
" #eta";
147 h_matchingObjectEta =
new TH1F( hlabel.c_str(), htitle.c_str(), nbineta,etamin,etamax);
148 hlabel=
"h_matchingObject_abseta"; htitle=type_+
" |#eta|";
149 h_matchingObjectAbsEta =
new TH1F( hlabel.c_str(), htitle.c_str(), nbineta/2,0.,etamax);
150 hlabel=
"h_matchingObject_P"; htitle=type_+
" p";
151 h_matchingObjectP =
new TH1F( hlabel.c_str(), htitle.c_str(), nbinp,0.,pmax);
152 hlabel=
"h_matchingObject_Pt"; htitle=type_+
" pt";
153 h_matchingObjectPt =
new TH1F( hlabel.c_str(),htitle.c_str(), nbinpteff,5.,ptmax);
154 hlabel=
"h_matchingObject_phi"; htitle=type_+
" phi";
155 h_matchingObjectPhi =
new TH1F( hlabel.c_str(), htitle.c_str(), nbinphi,phimin,phimax);
156 hlabel=
"h_matchingObject_z"; htitle=type_+
" z";
157 h_matchingObjectZ =
new TH1F( hlabel.c_str(), htitle.c_str(), nbinxyz, -25, 25 );
160 h_ele_EoverP_all =
new TH1F(
"h_ele_EoverP_all",
"ele E/P_{vertex}, all reco electrons", nbineop,0.,eopmax);
161 h_ele_EseedOP_all =
new TH1F(
"h_ele_EseedOP_all",
"ele E_{seed}/P_{vertex}, all reco electrons", nbineop,0.,eopmax);
162 h_ele_EoPout_all =
new TH1F(
"h_ele_EoPout_all",
"ele E_{seed}/P_{out}, all reco electrons", nbineop,0.,eopmax);
163 h_ele_EeleOPout_all =
new TH1F(
"h_ele_EeleOPout_all",
"ele E_{ele}/P_{out}, all reco electrons", nbineop,0.,eopmax);
164 h_ele_dEtaSc_propVtx_all =
new TH1F(
"h_ele_dEtaSc_propVtx_all",
"ele #eta_{sc} - #eta_{tr}, prop from vertex, all reco electrons", nbindetamatch,detamatchmin,detamatchmax);
165 h_ele_dPhiSc_propVtx_all =
new TH1F(
"h_ele_dPhiSc_propVtx_all",
"ele #phi_{sc} - #phi_{tr}, prop from vertex, all reco electrons", nbindphimatch,dphimatchmin,dphimatchmax);
166 h_ele_dEtaCl_propOut_all =
new TH1F(
"h_ele_dEtaCl_propOut_all",
"ele #eta_{cl} - #eta_{tr}, prop from outermost, all reco electrons", nbindetamatch,detamatchmin,detamatchmax);
167 h_ele_dPhiCl_propOut_all =
new TH1F(
"h_ele_dPhiCl_propOut_all",
"ele #phi_{cl} - #phi_{tr}, prop from outermost, all reco electrons", nbindphimatch,dphimatchmin,dphimatchmax);
168 h_ele_TIP_all =
new TH1F(
"h_ele_TIP_all",
"ele vertex transverse radius, all reco electrons", 100,0.,0.2);
169 h_ele_HoE_all =
new TH1F(
"h_ele_HoE_all",
"ele hadronic energy / em energy, all reco electrons", nbinhoe, hoemin, hoemax) ;
170 h_ele_vertexEta_all =
new TH1F(
"h_ele_vertexEta_all",
"ele eta, all reco electrons", nbineta,etamin,etamax);
171 h_ele_vertexEta_all->Sumw2();
172 h_ele_vertexPt_all =
new TH1F(
"h_ele_vertexPt_all",
"ele p_{T}, all reco electrons", nbinpteff,5.,ptmax);
173 h_ele_vertexPt_all->Sumw2();
174 h_ele_mee_all =
new TH1F(
"h_ele_mee_all",
"ele pairs invariant mass, all reco electrons", nbinmee, meemin, meemax );
175 h_ele_mee_os =
new TH1F(
"h_ele_mee_os",
"ele pairs invariant mass, opp. sign", nbinmee, meemin, meemax );
178 h_ele_E2mnE1vsMee_all =
new TH2F(
"h_ele_E2mnE1vsMee_all",
"E2 - E1 vs ele pairs invariant mass, all electrons", nbinmee, meemin, meemax, 100, -50., 50. );
179 h_ele_E2mnE1vsMee_egeg_all =
new TH2F(
"h_ele_E2mnE1vsMee_egeg_all",
"E2 - E1 vs ele pairs invariant mass, ecal driven pairs, all electrons", nbinmee, meemin, meemax, 100, -50., 50. );
183 h_ele_matchingObjectPt_matched =
new TH1F(
"h_ele_matchingObjectPt_matched",
"Efficiency vs matching jet p_{T}", nbinpteff,5.,ptmax);
184 h_ele_matchingObjectPt_matched->Sumw2();
185 h_ele_matchingObjectAbsEta_matched =
new TH1F(
"h_ele_matchingObjectAbsEta_matched",
"Efficiency vs matching jet |#eta|", nbineta/2,0.,2.5);
186 h_ele_matchingObjectAbsEta_matched->Sumw2();
187 h_ele_matchingObjectEta_matched =
new TH1F(
"h_ele_matchingObjectEta_matched",
"Efficiency vs matching jet #eta", nbineta,etamin,etamax);
188 h_ele_matchingObjectEta_matched ->Sumw2();
189 h_ele_matchingObjectPhi_matched =
new TH1F(
"h_ele_matchingObjectPhi_matched",
"Efficiency vs matching jet phi", nbinphi,phimin,phimax);
190 h_ele_matchingObjectPhi_matched->Sumw2();
191 h_ele_matchingObjectZ_matched =
new TH1F(
"h_ele_matchingObjectZ_matched",
"Efficiency vs matching jet z", nbinxyz,-25,25);
192 h_ele_matchingObjectZ_matched->Sumw2();
194 h_ele_charge =
new TH1F(
"h_ele_charge",
"ele charge", 5,-2.,2.);
195 h_ele_chargeVsEta =
new TH2F(
"h_ele_chargeVsEta",
"ele charge vs eta", nbineta2D,etamin,etamax,5,-2.,2.);
196 h_ele_chargeVsPhi =
new TH2F(
"h_ele_chargeVsPhi",
"ele charge vs phi", nbinphi2D,phimin,phimax,5,-2.,2.);
197 h_ele_chargeVsPt =
new TH2F(
"h_ele_chargeVsPt",
"ele charge vs pt", nbinpt,0.,100.,5,-2.,2.);
198 h_ele_vertexP =
new TH1F(
"h_ele_vertexP",
"ele momentum", nbinp,0.,pmax);
199 h_ele_vertexPt =
new TH1F(
"h_ele_vertexPt",
"ele transverse momentum", nbinpt,0.,ptmax);
200 h_ele_vertexPtVsEta =
new TH2F(
"h_ele_vertexPtVsEta",
"ele transverse momentum vs eta",nbineta2D,etamin,etamax,nbinpt2D,0.,ptmax);
201 h_ele_vertexPtVsPhi =
new TH2F(
"h_ele_vertexPtVsPhi",
"ele transverse momentum vs phi",nbinphi2D,phimin,phimax,nbinpt2D,0.,ptmax);
202 h_ele_vertexEta =
new TH1F(
"h_ele_vertexEta",
"ele momentum eta", nbineta,etamin,etamax);
203 h_ele_vertexEtaVsPhi =
new TH2F(
"h_ele_vertexEtaVsPhi",
"ele momentum eta vs phi",nbineta2D,etamin,etamax,nbinphi2D,phimin,phimax );
204 h_ele_vertexPhi =
new TH1F(
"h_ele_vertexPhi",
"ele momentum #phi", nbinphi,phimin,phimax);
205 h_ele_vertexX =
new TH1F(
"h_ele_vertexX",
"ele vertex x", nbinxyz,-0.1,0.1 );
206 h_ele_vertexY =
new TH1F(
"h_ele_vertexY",
"ele vertex y", nbinxyz,-0.1,0.1 );
207 h_ele_vertexZ =
new TH1F(
"h_ele_vertexZ",
"ele vertex z", nbinxyz,-25, 25 );
208 h_ele_vertexTIP =
new TH1F(
"h_ele_vertexTIP",
"ele transverse impact parameter (wrt bs)", 90,0.,0.15);
209 h_ele_vertexTIPVsEta =
new TH2F(
"h_ele_vertexTIPVsEta",
"ele transverse impact parameter (wrt bs) vs eta", nbineta2D,etamin,etamax,45,0.,0.15);
210 h_ele_vertexTIPVsPhi =
new TH2F(
"h_ele_vertexTIPVsPhi",
"ele transverse impact parameter (wrt bs) vs phi", nbinphi2D,phimin,phimax,45,0.,0.15);
211 h_ele_vertexTIPVsPt =
new TH2F(
"h_ele_vertexTIPVsPt",
"ele transverse impact parameter (wrt bs) vs Pt", nbinpt2D,0.,ptmax,45,0.,0.15);
213 h_ele_PoPmatchingObject =
new TH1F(
"h_ele_PoPmatchingObject",
"ele momentum / matching jet momentum", 75,0.,1.5);
214 h_ele_PoPmatchingObjectVsEta =
new TH2F(
"h_ele_PoPmatchingObjectVsEta",
"ele momentum / matching jet momentum vs eta", nbineta2D,etamin,etamax,50,0.,1.5);
215 h_ele_PoPmatchingObjectVsPhi =
new TH2F(
"h_ele_PoPmatchingObjectVsPhi",
"ele momentum / matching jet momentum vs phi", nbinphi2D,phimin,phimax,50,0.,1.5);
216 h_ele_PoPmatchingObjectVsPt =
new TH2F(
"h_ele_PoPmatchingObjectVsPt",
"ele momentum / matching jet momentum vs eta", nbinpt2D,0.,ptmax,50,0.,1.5);
217 h_ele_PoPmatchingObject_barrel =
new TH1F(
"h_ele_PoPmatchingObject_barrel",
"ele momentum / matching jet momentum, barrel",75,0.,1.5);
218 h_ele_PoPmatchingObject_endcaps =
new TH1F(
"h_ele_PoPmatchingObject_endcaps",
"ele momentum / matching jet momentum, endcaps",75,0.,1.5);
220 h_ele_EtaMnEtamatchingObject =
new TH1F(
"h_ele_EtaMnEtamatchingObject",
"ele momentum eta - matching jet eta",nbindeta,detamin,detamax);
221 h_ele_EtaMnEtamatchingObjectVsEta =
new TH2F(
"h_ele_EtaMnEtamatchingObjectVsEta",
"ele momentum eta - matching jet eta vs eta",nbineta2D,etamin,etamax,nbindeta/2,detamin,detamax);
222 h_ele_EtaMnEtamatchingObjectVsPhi =
new TH2F(
"h_ele_EtaMnEtamatchingObjectVsPhi",
"ele momentum eta - matching jet eta vs phi",nbinphi2D,phimin,phimax,nbindeta/2,detamin,detamax);
223 h_ele_EtaMnEtamatchingObjectVsPt =
new TH2F(
"h_ele_EtaMnEtamatchingObjectVsPt",
"ele momentum eta - matching jet eta vs pt",nbinpt,0.,ptmax,nbindeta/2,detamin,detamax);
224 h_ele_PhiMnPhimatchingObject =
new TH1F(
"h_ele_PhiMnPhimatchingObject",
"ele momentum phi - matching jet phi",nbindphi,dphimin,dphimax);
225 h_ele_PhiMnPhimatchingObject2 =
new TH1F(
"h_ele_PhiMnPhimatchingObject2",
"ele momentum phi - matching jet phi",nbindphimatch2D,dphimatchmin,dphimatchmax);
226 h_ele_PhiMnPhimatchingObjectVsEta =
new TH2F(
"h_ele_PhiMnPhimatchingObjectVsEta",
"ele momentum phi - matching SC phi vs eta",nbineta2D,etamin,etamax,nbindphi/2,dphimin,dphimax);
227 h_ele_PhiMnPhimatchingObjectVsPhi =
new TH2F(
"h_ele_PhiMnPhimatchingObjectVsPhi",
"ele momentum phi - matching SC phi vs phi",nbinphi2D,phimin,phimax,nbindphi/2,dphimin,dphimax);
228 h_ele_PhiMnPhimatchingObjectVsPt =
new TH2F(
"h_ele_PhiMnPhimatchingObjectVsPt",
"ele momentum phi - matching SC phi vs pt",nbinpt2D,0.,ptmax,nbindphi/2,dphimin,dphimax);
231 histSclEn_ =
new TH1F(
"h_scl_energy",
"ele supercluster energy",nbinp,0.,pmax);
232 histSclEoEmatchingObject_barrel =
new TH1F(
"h_scl_EoEmatchingObject_barrel",
"ele supercluster energy / matching jet energy, barrel",50,0.2,1.2);
233 histSclEoEmatchingObject_endcaps =
new TH1F(
"h_scl_EoEmatchingObject_endcaps",
"ele supercluster energy / matching jet energy, endcaps",50,0.2,1.2);
234 histSclEt_ =
new TH1F(
"h_scl_et",
"ele supercluster transverse energy",nbinpt,0.,ptmax);
235 histSclEtVsEta_ =
new TH2F(
"h_scl_etVsEta",
"ele supercluster transverse energy vs eta",nbineta2D,etamin,etamax,nbinpt,0.,ptmax);
236 histSclEtVsPhi_ =
new TH2F(
"h_scl_etVsPhi",
"ele supercluster transverse energy vs phi",nbinphi2D,phimin,phimax,nbinpt,0.,ptmax);
237 histSclEtaVsPhi_ =
new TH2F(
"h_scl_etaVsPhi",
"ele supercluster eta vs phi",nbinphi2D,phimin,phimax,nbineta2D,etamin,etamax);
238 histSclEta_ =
new TH1F(
"h_scl_eta",
"ele supercluster eta",nbineta,etamin,etamax);
239 histSclPhi_ =
new TH1F(
"h_scl_phi",
"ele supercluster phi",nbinphi,phimin,phimax);
241 histSclSigIEtaIEta_ =
new TH1F(
"h_scl_sigietaieta",
"ele supercluster sigma ieta ieta",100,0.,0.05);
242 histSclSigIEtaIEta_barrel_ =
new TH1F(
"h_scl_sigietaieta_barrel",
"ele supercluster sigma ieta ieta, barrel",100,0.,0.05);
243 histSclSigIEtaIEta_endcaps_ =
new TH1F(
"h_scl_sigietaieta_endcaps",
"ele supercluster sigma ieta ieta, endcaps",100,0.,0.05);
244 histSclE1x5_ =
new TH1F(
"h_scl_E1x5",
"ele supercluster energy in 1x5",nbinp,0.,pmax);
245 histSclE1x5_barrel_ =
new TH1F(
"h_scl_E1x5_barrel",
"ele supercluster energy in 1x5 barrel",nbinp,0.,pmax);
246 histSclE1x5_endcaps_ =
new TH1F(
"h_scl_E1x5_endcaps",
"ele supercluster energy in 1x5 endcaps",nbinp,0.,pmax);
247 histSclE2x5max_ =
new TH1F(
"h_scl_E2x5max",
"ele supercluster energy in 2x5 max",nbinp,0.,pmax);
248 histSclE2x5max_barrel_ =
new TH1F(
"h_scl_E2x5max_barrel",
"ele supercluster energy in 2x5 max barrel",nbinp,0.,pmax);
249 histSclE2x5max_endcaps_ =
new TH1F(
"h_scl_E2x5max_endcaps",
"ele supercluster energy in 2x5 max endcaps",nbinp,0.,pmax);
250 histSclE5x5_ =
new TH1F(
"h_scl_E5x5",
"ele supercluster energy in 5x5",nbinp,0.,pmax);
251 histSclE5x5_barrel_ =
new TH1F(
"h_scl_E5x5_barrel",
"ele supercluster energy in 5x5 barrel",nbinp,0.,pmax);
252 histSclE5x5_endcaps_ =
new TH1F(
"h_scl_E5x5_endcaps",
"ele supercluster energy in 5x5 endcaps",nbinp,0.,pmax);
255 h_ele_ambiguousTracks =
new TH1F(
"h_ele_ambiguousTracks",
"ele # ambiguous tracks", 5,0.,5.);
256 h_ele_ambiguousTracksVsEta =
new TH2F(
"h_ele_ambiguousTracksVsEta",
"ele # ambiguous tracks vs eta", nbineta2D,etamin,etamax,5,0.,5.);
257 h_ele_ambiguousTracksVsPhi =
new TH2F(
"h_ele_ambiguousTracksVsPhi",
"ele # ambiguous tracks vs phi", nbinphi2D,phimin,phimax,5,0.,5.);
258 h_ele_ambiguousTracksVsPt =
new TH2F(
"h_ele_ambiguousTracksVsPt",
"ele # ambiguous tracks vs pt", nbinpt2D,0.,ptmax,5,0.,5.);
259 h_ele_foundHits =
new TH1F(
"h_ele_foundHits",
"ele track # found hits", nbinfhits,0.,fhitsmax);
260 h_ele_foundHitsVsEta =
new TH2F(
"h_ele_foundHitsVsEta",
"ele track # found hits vs eta", nbineta2D,etamin,etamax,nbinfhits,0.,fhitsmax);
261 h_ele_foundHitsVsPhi =
new TH2F(
"h_ele_foundHitsVsPhi",
"ele track # found hits vs phi", nbinphi2D,phimin,phimax,nbinfhits,0.,fhitsmax);
262 h_ele_foundHitsVsPt =
new TH2F(
"h_ele_foundHitsVsPt",
"ele track # found hits vs pt", nbinpt2D,0.,ptmax,nbinfhits,0.,fhitsmax);
263 h_ele_lostHits =
new TH1F(
"h_ele_lostHits",
"ele track # lost hits", 5,0.,5.);
264 h_ele_lostHitsVsEta =
new TH2F(
"h_ele_lostHitsVsEta",
"ele track # lost hits vs eta", nbineta2D,etamin,etamax,nbinlhits,0.,lhitsmax);
265 h_ele_lostHitsVsPhi =
new TH2F(
"h_ele_lostHitsVsPhi",
"ele track # lost hits vs eta", nbinphi2D,phimin,phimax,nbinlhits,0.,lhitsmax);
266 h_ele_lostHitsVsPt =
new TH2F(
"h_ele_lostHitsVsPt",
"ele track # lost hits vs eta", nbinpt2D,0.,ptmax,nbinlhits,0.,lhitsmax);
267 h_ele_chi2 =
new TH1F(
"h_ele_chi2",
"ele track #chi^{2}", 100,0.,15.);
268 h_ele_chi2VsEta =
new TH2F(
"h_ele_chi2VsEta",
"ele track #chi^{2} vs eta", nbineta2D,etamin,etamax,50,0.,15.);
269 h_ele_chi2VsPhi =
new TH2F(
"h_ele_chi2VsPhi",
"ele track #chi^{2} vs phi", nbinphi2D,phimin,phimax,50,0.,15.);
270 h_ele_chi2VsPt =
new TH2F(
"h_ele_chi2VsPt",
"ele track #chi^{2} vs pt", nbinpt2D,0.,ptmax,50,0.,15.);
271 h_ele_PinMnPout =
new TH1F(
"h_ele_PinMnPout",
"ele track inner p - outer p, mean" ,nbinp,0.,200.);
272 h_ele_PinMnPout_mode =
new TH1F(
"h_ele_PinMnPout_mode",
"ele track inner p - outer p, mode" ,nbinp,0.,100.);
273 h_ele_PinMnPoutVsEta_mode =
new TH2F(
"h_ele_PinMnPoutVsEta_mode",
"ele track inner p - outer p vs eta, mode" ,nbineta2D, etamin,etamax,nbinp2D,0.,100.);
274 h_ele_PinMnPoutVsPhi_mode =
new TH2F(
"h_ele_PinMnPoutVsPhi_mode",
"ele track inner p - outer p vs phi, mode" ,nbinphi2D, phimin,phimax,nbinp2D,0.,100.);
275 h_ele_PinMnPoutVsPt_mode =
new TH2F(
"h_ele_PinMnPoutVsPt_mode",
"ele track inner p - outer p vs pt, mode" ,nbinpt2D, 0.,ptmax,nbinp2D,0.,100.);
276 h_ele_PinMnPoutVsE_mode =
new TH2F(
"h_ele_PinMnPoutVsE_mode",
"ele track inner p - outer p vs E, mode" ,nbinp2D, 0.,200.,nbinp2D,0.,100.);
277 h_ele_PinMnPoutVsChi2_mode =
new TH2F(
"h_ele_PinMnPoutVsChi2_mode",
"ele track inner p - outer p vs track chi2, mode" ,50, 0.,20.,nbinp2D,0.,100.);
278 h_ele_outerP =
new TH1F(
"h_ele_outerP",
"ele track outer p, mean", nbinp,0.,pmax);
279 h_ele_outerP_mode =
new TH1F(
"h_ele_outerP_mode",
"ele track outer p, mode", nbinp,0.,pmax);
280 h_ele_outerPVsEta_mode =
new TH2F(
"h_ele_outerPVsEta_mode",
"ele track outer p vs eta mode", nbineta2D,etamin,etamax,50,0.,pmax);
281 h_ele_outerPt =
new TH1F(
"h_ele_outerPt",
"ele track outer p_{T}, mean", nbinpt,0.,ptmax);
282 h_ele_outerPt_mode =
new TH1F(
"h_ele_outerPt_mode",
"ele track outer p_{T}, mode", nbinpt,0.,ptmax);
283 h_ele_outerPtVsEta_mode =
new TH2F(
"h_ele_outerPtVsEta_mode",
"ele track outer p_{T} vs eta, mode", nbineta2D,etamin,etamax,nbinpt2D,0.,ptmax);
284 h_ele_outerPtVsPhi_mode =
new TH2F(
"h_ele_outerPtVsPhi_mode",
"ele track outer p_{T} vs phi, mode", nbinphi2D,phimin,phimax,nbinpt2D,0.,ptmax);
285 h_ele_outerPtVsPt_mode =
new TH2F(
"h_ele_outerPtVsPt_mode",
"ele track outer p_{T} vs pt, mode", nbinpt2D,0.,100.,nbinpt2D,0.,ptmax);
288 h_ele_EoP =
new TH1F(
"h_ele_EoP",
"ele E/P_{vertex}", nbineop,0.,eopmax); h_ele_EoP->Sumw2();
289 h_ele_EoP_eg =
new TH1F(
"h_ele_EoP_eg",
"ele E/P_{vertex}, ecal driven", nbineop,0.,eopmax);
290 h_ele_EoP_eg->Sumw2();
291 h_ele_EoP_barrel =
new TH1F(
"h_ele_EoP_barrel",
"ele E/P_{vertex} barrel", nbineop,0.,eopmax);
292 h_ele_EoP_barrel->Sumw2();
293 h_ele_EoP_eg_barrel =
new TH1F(
"h_ele_EoP_eg_barrel",
"ele E/P_{vertex}, ecal driven barrel", nbineop,0.,eopmax);
294 h_ele_EoP_eg_barrel->Sumw2();
295 h_ele_EoP_endcaps =
new TH1F(
"h_ele_EoP_endcaps",
"ele E/P_{vertex} endcaps", nbineop,0.,eopmax);
296 h_ele_EoP_endcaps->Sumw2();
297 h_ele_EoP_eg_endcaps =
new TH1F(
"h_ele_EoP_eg_endcaps",
"ele E/P_{vertex}, ecal driven endcaps", nbineop,0.,eopmax);
298 h_ele_EoP_eg_endcaps->Sumw2();
299 h_ele_EoPVsEta =
new TH2F(
"h_ele_EoPVsEta",
"ele E/P_{vertex} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
300 h_ele_EoPVsPhi =
new TH2F(
"h_ele_EoPVsPhi",
"ele E/P_{vertex} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
301 h_ele_EoPVsE =
new TH2F(
"h_ele_EoPVsE",
"ele E/P_{vertex} vs E", 50,0.,pmax ,50,0.,5.);
302 h_ele_EseedOP =
new TH1F(
"h_ele_EseedOP",
"ele E_{seed}/P_{vertex}", nbineop,0.,eopmax);
303 h_ele_EseedOP->Sumw2();
304 h_ele_EseedOP_eg =
new TH1F(
"h_ele_EseedOP_eg",
"ele E_{seed}/P_{vertex}, ecal driven", nbineop,0.,eopmax);
305 h_ele_EseedOP_eg->Sumw2();
306 h_ele_EseedOP_barrel =
new TH1F(
"h_ele_EseedOP_barrel",
"ele E_{seed}/P_{vertex} barrel", nbineop,0.,eopmax);
307 h_ele_EseedOP_barrel->Sumw2();
308 h_ele_EseedOP_eg_barrel =
new TH1F(
"h_ele_EseedOP_eg_barrel",
"ele E_{seed}/P_{vertex}, ecal driven barrel", nbineop,0.,eopmax);
309 h_ele_EseedOP_eg_barrel->Sumw2();
310 h_ele_EseedOP_endcaps =
new TH1F(
"h_ele_EseedOP_endcaps",
"ele E_{seed}/P_{vertex} endcaps", nbineop,0.,eopmax);
311 h_ele_EseedOP_endcaps->Sumw2();
312 h_ele_EseedOP_eg_endcaps =
new TH1F(
"h_ele_EseedOP_eg_endcaps",
"ele E_{seed}/P_{vertex}, ecal driven, endcaps", nbineop,0.,eopmax);
313 h_ele_EseedOP_eg_endcaps->Sumw2();
314 h_ele_EseedOPVsEta =
new TH2F(
"h_ele_EseedOPVsEta",
"ele E_{seed}/P_{vertex} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
315 h_ele_EseedOPVsPhi =
new TH2F(
"h_ele_EseedOPVsPhi",
"ele E_{seed}/P_{vertex} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
316 h_ele_EseedOPVsE =
new TH2F(
"h_ele_EseedOPVsE",
"ele E_{seed}/P_{vertex} vs E", 50,0.,pmax ,50,0.,5.);
317 h_ele_EoPout =
new TH1F(
"h_ele_EoPout",
"ele E_{seed}/P_{out}", nbineop,0.,eopmax);
318 h_ele_EoPout->Sumw2();
319 h_ele_EoPout_eg =
new TH1F(
"h_ele_EoPout_eg",
"ele E_{seed}/P_{out}, ecal driven", nbineop,0.,eopmax);
320 h_ele_EoPout_eg->Sumw2();
321 h_ele_EoPout_barrel =
new TH1F(
"h_ele_EoPout_barrel",
"ele E_{seed}/P_{out} barrel", nbineop,0.,eopmax);
322 h_ele_EoPout_barrel->Sumw2();
323 h_ele_EoPout_eg_barrel =
new TH1F(
"h_ele_EoPout_eg_barrel",
"ele E_{seed}/P_{out}, ecal driven, barrel", nbineop,0.,eopmax);
324 h_ele_EoPout_eg_barrel->Sumw2();
325 h_ele_EoPout_endcaps =
new TH1F(
"h_ele_EoPout_endcaps",
"ele E_{seed}/P_{out} endcaps", nbineop,0.,eopmax);
326 h_ele_EoPout_endcaps->Sumw2();
327 h_ele_EoPout_eg_endcaps =
new TH1F(
"h_ele_EoPout_eg_endcaps",
"ele E_{seed}/P_{out}, ecal driven, endcaps", nbineop,0.,eopmax);
328 h_ele_EoPout_eg_endcaps->Sumw2();
329 h_ele_EoPoutVsEta =
new TH2F(
"h_ele_EoPoutVsEta",
"ele E/P_{out} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
330 h_ele_EoPoutVsPhi =
new TH2F(
"h_ele_EoPoutVsPhi",
"ele E/P_{out} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
331 h_ele_EoPoutVsE =
new TH2F(
"h_ele_EoPoutVsE",
"ele E/P_{out} vs E", nbinp2D,0.,pmax,nbineop2D,0.,eopmaxsht);
332 h_ele_EeleOPout =
new TH1F(
"h_ele_EeleOPout",
"ele E_{ele}/P_{out}", nbineop,0.,eopmax);
333 h_ele_EeleOPout->Sumw2();
334 h_ele_EeleOPout_eg =
new TH1F(
"h_ele_EeleOPout_eg",
"ele E_{ele}/P_{out}, ecal driven", nbineop,0.,eopmax);
335 h_ele_EeleOPout_eg->Sumw2();
336 h_ele_EeleOPout_barrel =
new TH1F(
"h_ele_EeleOPout_barrel",
"ele E_{ele}/P_{out} barrel", nbineop,0.,eopmax);
337 h_ele_EeleOPout_barrel->Sumw2();
338 h_ele_EeleOPout_eg_barrel =
new TH1F(
"h_ele_EeleOPout_eg_barrel",
"ele E_{ele}/P_{out}, ecal driven, barrel", nbineop,0.,eopmax);
339 h_ele_EeleOPout_eg_barrel->Sumw2();
340 h_ele_EeleOPout_endcaps =
new TH1F(
"h_ele_EeleOPout_endcaps",
"ele E_{ele}/P_{out} endcaps", nbineop,0.,eopmax);
341 h_ele_EeleOPout_endcaps->Sumw2();
342 h_ele_EeleOPout_eg_endcaps =
new TH1F(
"h_ele_EeleOPout_eg_endcaps",
"ele E_{ele}/P_{out}, ecal driven, endcaps", nbineop,0.,eopmax);
343 h_ele_EeleOPout_eg_endcaps->Sumw2();
344 h_ele_EeleOPoutVsEta =
new TH2F(
"h_ele_EeleOPoutVsEta",
"ele E_{ele}/P_{out} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
345 h_ele_EeleOPoutVsPhi =
new TH2F(
"h_ele_EeleOPoutVsPhi",
"ele E_{ele}/P_{out} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
346 h_ele_EeleOPoutVsE =
new TH2F(
"h_ele_EeleOPoutVsE",
"ele E_{ele}/P_{out} vs E", nbinp2D,0.,pmax,nbineop2D,0.,eopmaxsht);
347 h_ele_dEtaSc_propVtx =
new TH1F(
"h_ele_dEtaSc_propVtx",
"ele #eta_{sc} - #eta_{tr}, prop from vertex", nbindetamatch,detamatchmin,detamatchmax);
348 h_ele_dEtaSc_propVtx->Sumw2();
349 h_ele_dEtaSc_propVtx_eg =
new TH1F(
"h_ele_dEtaSc_propVtx_eg",
"ele #eta_{sc} - #eta_{tr}, prop from vertex, ecal driven", nbindetamatch,detamatchmin,detamatchmax);
350 h_ele_dEtaSc_propVtx_eg->Sumw2();
351 h_ele_dEtaSc_propVtx_barrel =
new TH1F(
"h_ele_dEtaSc_propVtx_barrel",
"ele #eta_{sc} - #eta_{tr}, prop from vertex, barrel", nbindetamatch,detamatchmin,detamatchmax);
352 h_ele_dEtaSc_propVtx_barrel->Sumw2();
353 h_ele_dEtaSc_propVtx_eg_barrel =
new TH1F(
"h_ele_dEtaSc_propVtx_eg_barrel",
"ele #eta_{sc} - #eta_{tr}, prop from vertex, ecal driven, barrel", nbindetamatch,detamatchmin,detamatchmax);
354 h_ele_dEtaSc_propVtx_eg_barrel->Sumw2();
355 h_ele_dEtaSc_propVtx_endcaps =
new TH1F(
"h_ele_dEtaSc_propVtx_endcaps",
"ele #eta_{sc} - #eta_{tr}, prop from vertex, endcaps", nbindetamatch,detamatchmin,detamatchmax);
356 h_ele_dEtaSc_propVtx_endcaps->Sumw2();
357 h_ele_dEtaSc_propVtx_eg_endcaps =
new TH1F(
"h_ele_dEtaSc_propVtx_eg_endcaps",
"ele #eta_{sc} - #eta_{tr}, prop from vertex, ecal driven, endcaps", nbindetamatch,detamatchmin,detamatchmax);
358 h_ele_dEtaSc_propVtx_eg_endcaps->Sumw2();
359 h_ele_dEtaScVsEta_propVtx =
new TH2F(
"h_ele_dEtaScVsEta_propVtx",
"ele #eta_{sc} - #eta_{tr} vs eta, prop from vertex", nbineta2D,etamin,etamax,nbindetamatch2D,detamatchmin,detamatchmax);
360 h_ele_dEtaScVsPhi_propVtx =
new TH2F(
"h_ele_dEtaScVsPhi_propVtx",
"ele #eta_{sc} - #eta_{tr} vs phi, prop from vertex", nbinphi2D,phimin,phimax,nbindetamatch2D,detamatchmin,detamatchmax);
361 h_ele_dEtaScVsPt_propVtx =
new TH2F(
"h_ele_dEtaScVsPt_propVtx",
"ele #eta_{sc} - #eta_{tr} vs pt, prop from vertex", nbinpt2D,0.,ptmax,nbindetamatch2D,detamatchmin,detamatchmax);
362 h_ele_dPhiSc_propVtx =
new TH1F(
"h_ele_dPhiSc_propVtx",
"ele #phi_{sc} - #phi_{tr}, prop from vertex", nbindphimatch,dphimatchmin,dphimatchmax);
363 h_ele_dPhiSc_propVtx->Sumw2();
364 h_ele_dPhiSc_propVtx_eg =
new TH1F(
"h_ele_dPhiSc_propVtx_eg",
"ele #phi_{sc} - #phi_{tr}, prop from vertex, ecal driven", nbindphimatch,dphimatchmin,dphimatchmax);
365 h_ele_dPhiSc_propVtx_eg->Sumw2();
366 h_ele_dPhiSc_propVtx_barrel =
new TH1F(
"h_ele_dPhiSc_propVtx_barrel",
"ele #phi_{sc} - #phi_{tr}, prop from vertex, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
367 h_ele_dPhiSc_propVtx_barrel->Sumw2();
368 h_ele_dPhiSc_propVtx_eg_barrel =
new TH1F(
"h_ele_dPhiSc_propVtx_eg_barrel",
"ele #phi_{sc} - #phi_{tr}, prop from vertex, ecal driven, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
369 h_ele_dPhiSc_propVtx_eg_barrel->Sumw2();
370 h_ele_dPhiSc_propVtx_endcaps =
new TH1F(
"h_ele_dPhiSc_propVtx_endcaps",
"ele #phi_{sc} - #phi_{tr}, prop from vertex, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
371 h_ele_dPhiSc_propVtx_endcaps->Sumw2();
372 h_ele_dPhiSc_propVtx_eg_endcaps =
new TH1F(
"h_ele_dPhiSc_propVtx_eg_endcaps",
"ele #phi_{sc} - #phi_{tr}, prop from vertex, ecal driven, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
373 h_ele_dPhiSc_propVtx_eg_endcaps->Sumw2();
374 h_ele_dPhiScVsEta_propVtx =
new TH2F(
"h_ele_dPhiScVsEta_propVtx",
"ele #phi_{sc} - #phi_{tr} vs eta, prop from vertex", nbineta2D,etamin,etamax,nbindphimatch2D,dphimatchmin,dphimatchmax);
375 h_ele_dPhiScVsPhi_propVtx =
new TH2F(
"h_ele_dPhiScVsPhi_propVtx",
"ele #phi_{sc} - #phi_{tr} vs phi, prop from vertex", nbinphi2D,phimin,phimax,nbindphimatch2D,dphimatchmin,dphimatchmax);
376 h_ele_dPhiScVsPt_propVtx =
new TH2F(
"h_ele_dPhiScVsPt_propVtx",
"ele #phi_{sc} - #phi_{tr} vs pt, prop from vertex", nbinpt2D,0.,ptmax,nbindphimatch2D,dphimatchmin,dphimatchmax);
377 h_ele_dEtaCl_propOut =
new TH1F(
"h_ele_dEtaCl_propOut",
"ele #eta_{cl} - #eta_{tr}, prop from outermost", nbindetamatch,detamatchmin,detamatchmax);
378 h_ele_dEtaCl_propOut->Sumw2();
379 h_ele_dEtaCl_propOut_eg =
new TH1F(
"h_ele_dEtaCl_propOut_eg",
"ele #eta_{cl} - #eta_{tr}, prop from outermost, ecal driven", nbindetamatch,detamatchmin,detamatchmax);
380 h_ele_dEtaCl_propOut_eg->Sumw2();
381 h_ele_dEtaCl_propOut_barrel =
new TH1F(
"h_ele_dEtaCl_propOut_barrel",
"ele #eta_{cl} - #eta_{tr}, prop from outermost, barrel", nbindetamatch,detamatchmin,detamatchmax);
382 h_ele_dEtaCl_propOut_barrel->Sumw2();
383 h_ele_dEtaCl_propOut_eg_barrel =
new TH1F(
"h_ele_dEtaCl_propOut_eg_barrel",
"ele #eta_{cl} - #eta_{tr}, prop from outermost, ecal driven, barrel", nbindetamatch,detamatchmin,detamatchmax);
384 h_ele_dEtaCl_propOut_eg_barrel->Sumw2();
385 h_ele_dEtaCl_propOut_endcaps =
new TH1F(
"h_ele_dEtaCl_propOut_endcaps",
"ele #eta_{cl} - #eta_{tr}, prop from outermost, endcaps", nbindetamatch,detamatchmin,detamatchmax);
386 h_ele_dEtaCl_propOut_endcaps->Sumw2();
387 h_ele_dEtaCl_propOut_eg_endcaps =
new TH1F(
"h_ele_dEtaCl_propOut_eg_endcaps",
"ele #eta_{cl} - #eta_{tr}, prop from outermost, ecal driven, endcaps", nbindetamatch,detamatchmin,detamatchmax);
388 h_ele_dEtaCl_propOut_eg_endcaps->Sumw2();
389 h_ele_dEtaClVsEta_propOut =
new TH2F(
"h_ele_dEtaClVsEta_propOut",
"ele #eta_{cl} - #eta_{tr} vs eta, prop from out", nbineta2D,etamin,etamax,nbindetamatch2D,detamatchmin,detamatchmax);
390 h_ele_dEtaClVsPhi_propOut =
new TH2F(
"h_ele_dEtaClVsPhi_propOut",
"ele #eta_{cl} - #eta_{tr} vs phi, prop from out", nbinphi2D,phimin,phimax,nbindetamatch2D,detamatchmin,detamatchmax);
391 h_ele_dEtaClVsPt_propOut =
new TH2F(
"h_ele_dEtaScVsPt_propOut",
"ele #eta_{cl} - #eta_{tr} vs pt, prop from out", nbinpt2D,0.,ptmax,nbindetamatch2D,detamatchmin,detamatchmax);
392 h_ele_dPhiCl_propOut =
new TH1F(
"h_ele_dPhiCl_propOut",
"ele #phi_{cl} - #phi_{tr}, prop from outermost", nbindphimatch,dphimatchmin,dphimatchmax);
393 h_ele_dPhiCl_propOut->Sumw2();
394 h_ele_dPhiCl_propOut_eg =
new TH1F(
"h_ele_dPhiCl_propOut_eg",
"ele #phi_{cl} - #phi_{tr}, prop from outermost, ecal driven", nbindphimatch,dphimatchmin,dphimatchmax);
395 h_ele_dPhiCl_propOut_eg->Sumw2();
396 h_ele_dPhiCl_propOut_barrel =
new TH1F(
"h_ele_dPhiCl_propOut_barrel",
"ele #phi_{cl} - #phi_{tr}, prop from outermost, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
397 h_ele_dPhiCl_propOut_barrel->Sumw2();
398 h_ele_dPhiCl_propOut_eg_barrel =
new TH1F(
"h_ele_dPhiCl_propOut_eg_barrel",
"ele #phi_{cl} - #phi_{tr}, prop from outermost, ecal driven, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
399 h_ele_dPhiCl_propOut_eg_barrel->Sumw2();
400 h_ele_dPhiCl_propOut_endcaps =
new TH1F(
"h_ele_dPhiCl_propOut_endcaps",
"ele #phi_{cl} - #phi_{tr}, prop from outermost, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
401 h_ele_dPhiCl_propOut_endcaps->Sumw2();
402 h_ele_dPhiCl_propOut_eg_endcaps =
new TH1F(
"h_ele_dPhiCl_propOut_eg_endcaps",
"ele #phi_{cl} - #phi_{tr}, prop from outermost, ecal driven, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
403 h_ele_dPhiCl_propOut_eg_endcaps->Sumw2();
404 h_ele_dPhiClVsEta_propOut =
new TH2F(
"h_ele_dPhiClVsEta_propOut",
"ele #phi_{cl} - #phi_{tr} vs eta, prop from out", nbineta2D,etamin,etamax,nbindphimatch2D,dphimatchmin,dphimatchmax);
405 h_ele_dPhiClVsPhi_propOut =
new TH2F(
"h_ele_dPhiClVsPhi_propOut",
"ele #phi_{cl} - #phi_{tr} vs phi, prop from out", nbinphi2D,phimin,phimax,nbindphimatch2D,dphimatchmin,dphimatchmax);
406 h_ele_dPhiClVsPt_propOut =
new TH2F(
"h_ele_dPhiSClsPt_propOut",
"ele #phi_{cl} - #phi_{tr} vs pt, prop from out", nbinpt2D,0.,ptmax,nbindphimatch2D,dphimatchmin,dphimatchmax);
407 h_ele_dEtaEleCl_propOut =
new TH1F(
"h_ele_dEtaEleCl_propOut",
"ele #eta_{EleCl} - #eta_{tr}, prop from outermost", nbindetamatch,detamatchmin,detamatchmax);
408 h_ele_dEtaEleCl_propOut->Sumw2();
409 h_ele_dEtaEleCl_propOut_eg =
new TH1F(
"h_ele_dEtaEleCl_propOut_eg",
"ele #eta_{EleCl} - #eta_{tr}, prop from outermost, ecal driven", nbindetamatch,detamatchmin,detamatchmax);
410 h_ele_dEtaEleCl_propOut_eg->Sumw2();
411 h_ele_dEtaEleCl_propOut_barrel =
new TH1F(
"h_ele_dEtaEleCl_propOut_barrel",
"ele #eta_{EleCl} - #eta_{tr}, prop from outermost, barrel", nbindetamatch,detamatchmin,detamatchmax);
412 h_ele_dEtaEleCl_propOut_barrel->Sumw2();
413 h_ele_dEtaEleCl_propOut_eg_barrel =
new TH1F(
"h_ele_dEtaEleCl_propOut_eg_barrel",
"ele #eta_{EleCl} - #eta_{tr}, prop from outermost, ecal driven, barrel", nbindetamatch,detamatchmin,detamatchmax);
414 h_ele_dEtaEleCl_propOut_eg_barrel->Sumw2();
415 h_ele_dEtaEleCl_propOut_endcaps =
new TH1F(
"h_ele_dEtaEleCl_propOut_endcaps",
"ele #eta_{EleCl} - #eta_{tr}, prop from outermost, endcaps", nbindetamatch,detamatchmin,detamatchmax);
416 h_ele_dEtaEleCl_propOut_endcaps->Sumw2();
417 h_ele_dEtaEleCl_propOut_eg_endcaps =
new TH1F(
"h_ele_dEtaEleCl_propOut_eg_endcaps",
"ele #eta_{EleCl} - #eta_{tr}, prop from outermost, ecal driven, endcaps", nbindetamatch,detamatchmin,detamatchmax);
418 h_ele_dEtaEleCl_propOut_eg_endcaps->Sumw2();
419 h_ele_dEtaEleClVsEta_propOut =
new TH2F(
"h_ele_dEtaEleClVsEta_propOut",
"ele #eta_{EleCl} - #eta_{tr} vs eta, prop from out", nbineta2D,etamin,etamax,nbindetamatch2D,detamatchmin,detamatchmax);
420 h_ele_dEtaEleClVsPhi_propOut =
new TH2F(
"h_ele_dEtaEleClVsPhi_propOut",
"ele #eta_{EleCl} - #eta_{tr} vs phi, prop from out", nbinphi2D,phimin,phimax,nbindetamatch2D,detamatchmin,detamatchmax);
421 h_ele_dEtaEleClVsPt_propOut =
new TH2F(
"h_ele_dEtaScVsPt_propOut",
"ele #eta_{EleCl} - #eta_{tr} vs pt, prop from out", nbinpt2D,0.,ptmax,nbindetamatch2D,detamatchmin,detamatchmax);
422 h_ele_dPhiEleCl_propOut =
new TH1F(
"h_ele_dPhiEleCl_propOut",
"ele #phi_{EleCl} - #phi_{tr}, prop from outermost", nbindphimatch,dphimatchmin,dphimatchmax);
423 h_ele_dPhiEleCl_propOut->Sumw2();
424 h_ele_dPhiEleCl_propOut_eg =
new TH1F(
"h_ele_dPhiEleCl_propOut_eg",
"ele #phi_{EleCl} - #phi_{tr}, prop from outermost, ecal driven", nbindphimatch,dphimatchmin,dphimatchmax);
425 h_ele_dPhiEleCl_propOut_eg->Sumw2();
426 h_ele_dPhiEleCl_propOut_barrel =
new TH1F(
"h_ele_dPhiEleCl_propOut_barrel",
"ele #phi_{EleCl} - #phi_{tr}, prop from outermost, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
427 h_ele_dPhiEleCl_propOut_barrel->Sumw2();
428 h_ele_dPhiEleCl_propOut_eg_barrel =
new TH1F(
"h_ele_dPhiEleCl_propOut_eg_barrel",
"ele #phi_{EleCl} - #phi_{tr}, prop from outermost, ecal driven, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
429 h_ele_dPhiEleCl_propOut_eg_barrel->Sumw2();
430 h_ele_dPhiEleCl_propOut_endcaps =
new TH1F(
"h_ele_dPhiEleCl_propOut_endcaps",
"ele #phi_{EleCl} - #phi_{tr}, prop from outermost, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
431 h_ele_dPhiEleCl_propOut_endcaps->Sumw2();
432 h_ele_dPhiEleCl_propOut_eg_endcaps =
new TH1F(
"h_ele_dPhiEleCl_propOut_eg_endcaps",
"ele #phi_{EleCl} - #phi_{tr}, prop from outermost, ecal driven, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
433 h_ele_dPhiEleCl_propOut_eg_endcaps->Sumw2();
434 h_ele_dPhiEleClVsEta_propOut =
new TH2F(
"h_ele_dPhiEleClVsEta_propOut",
"ele #phi_{EleCl} - #phi_{tr} vs eta, prop from out", nbineta2D,etamin,etamax,nbindphimatch2D,dphimatchmin,dphimatchmax);
435 h_ele_dPhiEleClVsPhi_propOut =
new TH2F(
"h_ele_dPhiEleClVsPhi_propOut",
"ele #phi_{EleCl} - #phi_{tr} vs phi, prop from out", nbinphi2D,phimin,phimax,nbindphimatch2D,dphimatchmin,dphimatchmax);
436 h_ele_dPhiEleClVsPt_propOut =
new TH2F(
"h_ele_dPhiSEleClsPt_propOut",
"ele #phi_{EleCl} - #phi_{tr} vs pt, prop from out", nbinpt2D,0.,ptmax,nbindphimatch2D,dphimatchmin,dphimatchmax);
438 h_ele_HoE =
new TH1F(
"h_ele_HoE",
"ele hadronic energy / em energy", nbinhoe, hoemin, hoemax) ;
440 h_ele_HoE_eg =
new TH1F(
"h_ele_HoE_eg",
"ele hadronic energy / em energy, ecal driven", nbinhoe, hoemin, hoemax) ;
441 h_ele_HoE_eg->Sumw2();
442 h_ele_HoE_barrel =
new TH1F(
"h_ele_HoE_barrel",
"ele hadronic energy / em energy, barrel", nbinhoe, hoemin, hoemax) ;
443 h_ele_HoE_barrel->Sumw2();
444 h_ele_HoE_eg_barrel =
new TH1F(
"h_ele_HoE_eg_barrel",
"ele hadronic energy / em energy, ecal driven, barrel", nbinhoe, hoemin, hoemax) ;
445 h_ele_HoE_eg_barrel->Sumw2();
446 h_ele_HoE_endcaps =
new TH1F(
"h_ele_HoE_endcaps",
"ele hadronic energy / em energy, endcaps", nbinhoe, hoemin, hoemax) ;
447 h_ele_HoE_endcaps->Sumw2();
448 h_ele_HoE_eg_endcaps =
new TH1F(
"h_ele_HoE_eg_endcaps",
"ele hadronic energy / em energy, ecal driven, endcaps", nbinhoe, hoemin, hoemax) ;
449 h_ele_HoE_eg_endcaps->Sumw2();
450 h_ele_HoE_fiducial =
new TH1F(
"h_ele_HoE_fiducial",
"ele hadronic energy / em energy",nbinhoe, hoemin, hoemax) ;
451 h_ele_HoEVsEta =
new TH2F(
"h_ele_HoEVsEta",
"ele hadronic energy / em energy vs eta", nbineta,etamin,etamax,nbinhoe, hoemin, hoemax) ;
452 h_ele_HoEVsPhi =
new TH2F(
"h_ele_HoEVsPhi",
"ele hadronic energy / em energy vs phi", nbinphi2D,phimin,phimax,nbinhoe, hoemin, hoemax) ;
453 h_ele_HoEVsE =
new TH2F(
"h_ele_HoEVsE",
"ele hadronic energy / em energy vs E", nbinp, 0.,300.,nbinhoe, hoemin, hoemax) ;
455 h_ele_seed_dphi2_ =
new TH1F(
"h_ele_seedDphi2",
"ele seed dphi 2nd layer", 50,-0.003,+0.003) ;
456 h_ele_seed_dphi2VsEta_ =
new TH2F(
"h_ele_seedDphi2VsEta",
"ele seed dphi 2nd layer vs eta", nbineta2D,etamin,etamax,50,-0.003,+0.003) ;
457 h_ele_seed_dphi2VsPt_ =
new TH2F(
"h_ele_seedDphi2VsPt",
"ele seed dphi 2nd layer vs pt", nbinpt2D,0.,ptmax,50,-0.003,+0.003) ;
458 h_ele_seed_drz2_ =
new TH1F(
"h_ele_seedDrz2",
"ele seed dr/dz 2nd layer", 50,-0.03,+0.03) ;
459 h_ele_seed_drz2VsEta_ =
new TH2F(
"h_ele_seedDrz2VsEta",
"ele seed dr/dz 2nd layer vs eta", nbineta2D,etamin,etamax,50,-0.03,+0.03) ;
460 h_ele_seed_drz2VsPt_ =
new TH2F(
"h_ele_seedDrz2VsPt",
"ele seed dr/dz 2nd layer vs pt", nbinpt2D,0.,ptmax,50,-0.03,+0.03) ;
461 h_ele_seed_subdet2_ =
new TH1F(
"h_ele_seedSubdet2",
"ele seed subdet 2nd layer", 10,0.,10.) ;
464 h_ele_classes =
new TH1F(
"h_ele_classes",
"electron classes", 20,0.0,20.);
465 h_ele_eta =
new TH1F(
"h_ele_eta",
"ele electron eta", nbineta/2,0.0,etamax);
466 h_ele_eta_golden =
new TH1F(
"h_ele_eta_golden",
"ele electron eta golden", nbineta/2,0.0,etamax);
467 h_ele_eta_bbrem =
new TH1F(
"h_ele_eta_bbrem",
"ele electron eta bbrem", nbineta/2,0.0,etamax);
468 h_ele_eta_narrow =
new TH1F(
"h_ele_eta_narrow",
"ele electron eta narrow", nbineta/2,0.0,etamax);
469 h_ele_eta_shower =
new TH1F(
"h_ele_eta_show",
"el, electron eta showering", nbineta/2,0.0,etamax);
470 h_ele_PinVsPoutGolden_mode =
new TH2F(
"h_ele_PinVsPoutGolden_mode",
"ele track inner p vs outer p, golden, mode" ,nbinp2D,0.,pmax,50,0.,pmax);
471 h_ele_PinVsPoutShowering_mode =
new TH2F(
"h_ele_PinVsPoutShowering_mode",
"ele track inner p vs outer p vs eta, showering, mode" ,nbinp2D,0.,pmax,50,0.,pmax);
472 h_ele_PinVsPoutGolden_mean =
new TH2F(
"h_ele_PinVsPoutGolden_mean",
"ele track inner p vs outer p, golden, mean" ,nbinp2D,0.,pmax,50,0.,pmax);
473 h_ele_PinVsPoutShowering_mean =
new TH2F(
"h_ele_PinVsPoutShowering_mean",
"ele track inner p vs outer p, showering, mean" ,nbinp2D,0.,pmax,50,0.,pmax);
474 h_ele_PtinVsPtoutGolden_mode =
new TH2F(
"h_ele_PtinVsPtoutGolden_mode",
"ele track inner pt vs outer pt, golden, mode" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
475 h_ele_PtinVsPtoutShowering_mode =
new TH2F(
"h_ele_PtinVsPtoutShowering_mode",
"ele track inner pt vs outer pt, showering, mode" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
476 h_ele_PtinVsPtoutGolden_mean =
new TH2F(
"h_ele_PtinVsPtoutGolden_mean",
"ele track inner pt vs outer pt, golden, mean" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
477 h_ele_PtinVsPtoutShowering_mean =
new TH2F(
"h_ele_PtinVsPtoutShowering_mean",
"ele track inner pt vs outer pt, showering, mean" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
478 histSclEoEmatchingObjectGolden_barrel =
new TH1F(
"h_scl_EoEmatchingObject golden, barrel",
"ele SC energy over matching jet energy, golden, barrel",100,0.2,1.2);
479 histSclEoEmatchingObjectGolden_endcaps =
new TH1F(
"h_scl_EoEmatchingObject golden, endcaps",
"ele SC energy over matching jet energy, golden, endcaps",100,0.2,1.2);
480 histSclEoEmatchingObjectShowering_barrel =
new TH1F(
"h_scl_EoEmatchingObject showering, barrel",
"ele SC energy over matching jet energy, showering, barrel",100,0.2,1.2);
481 histSclEoEmatchingObjectShowering_endcaps =
new TH1F(
"h_scl_EoEmatchingObject showering, endcaps",
"ele SC energy over matching jet energy, showering, endcaps",100,0.2,1.2);
484 h_ele_tkSumPt_dr03 =
new TH1F(
"h_ele_tkSumPt_dr03",
"tk isolation sum, dR=0.3",100,0.0,20.);
485 h_ele_ecalRecHitSumEt_dr03=
new TH1F(
"h_ele_ecalRecHitSumEt_dr03",
"ecal isolation sum, dR=0.3",100,0.0,20.);
486 h_ele_hcalDepth1TowerSumEt_dr03=
new TH1F(
"h_ele_hcalDepth1TowerSumEt_dr03",
"hcal depth1 isolation sum, dR=0.3",100,0.0,20.);
487 h_ele_hcalDepth2TowerSumEt_dr03=
new TH1F(
"h_ele_hcalDepth2TowerSumEt_dr03",
"hcal depth2 isolation sum, dR=0.3",100,0.0,20.);
488 h_ele_tkSumPt_dr04=
new TH1F(
"h_ele_tkSumPt_dr04",
"hcal isolation sum",100,0.0,20.);
489 h_ele_ecalRecHitSumEt_dr04=
new TH1F(
"h_ele_ecalRecHitSumEt_dr04",
"ecal isolation sum, dR=0.4",100,0.0,20.);
490 h_ele_hcalDepth1TowerSumEt_dr04=
new TH1F(
"h_ele_hcalDepth1TowerSumEt_dr04",
"hcal depth1 isolation sum, dR=0.4",100,0.0,20.);
491 h_ele_hcalDepth2TowerSumEt_dr04=
new TH1F(
"h_ele_hcalDepth2TowerSumEt_dr04",
"hcal depth2 isolation sum, dR=0.4",100,0.0,20.);
494 h_ele_fbrem =
new TH1F(
"h_ele_fbrem",
"ele brem fraction, mode",100,0.,1.);
495 h_ele_fbremVsEta_mode =
new TProfile(
"h_ele_fbremvsEtamode",
"mean ele brem fraction vs eta, mode",nbineta2D,etamin,etamax,0.,1.);
496 h_ele_fbremVsEta_mean =
new TProfile(
"h_ele_fbremvsEtamean",
"mean ele brem fraction vs eta, mean",nbineta2D,etamin,etamax,0.,1.);
499 h_ele_mva =
new TH1F(
"h_ele_mva",
"ele mva",100,-1.,1.);
500 h_ele_provenance =
new TH1F(
"h_ele_provenance",
"ele provenance",5,-2.,3.);
503 h_matchingObjectNum -> GetXaxis()-> SetTitle(
"N_{GenJet}");
504 h_matchingObjectNum -> GetYaxis()-> SetTitle(
"Events");
505 h_matchingObjectEta -> GetXaxis()-> SetTitle(
"#eta_{GenJet}");
506 h_matchingObjectEta -> GetYaxis()-> SetTitle(
"Events");
507 h_matchingObjectP -> GetXaxis()-> SetTitle(
"E_{GenJet} (GeV)");
508 h_matchingObjectP -> GetYaxis()-> SetTitle(
"Events");
509 h_ele_foundHits -> GetXaxis()-> SetTitle(
"N_{hits}");
510 h_ele_foundHits -> GetYaxis()-> SetTitle(
"Events");
511 h_ele_ambiguousTracks -> GetXaxis()-> SetTitle(
"N_{ambiguous tracks}");
512 h_ele_ambiguousTracks -> GetYaxis()-> SetTitle(
"Events");
513 h_ele_lostHits -> GetXaxis()-> SetTitle(
"N_{lost hits}");
514 h_ele_lostHits -> GetYaxis()-> SetTitle(
"Events");
515 h_ele_chi2 -> GetXaxis()-> SetTitle(
"#Chi^{2}");
516 h_ele_chi2 -> GetYaxis()-> SetTitle(
"Events");
517 h_ele_charge -> GetXaxis()-> SetTitle(
"charge");
518 h_ele_charge -> GetYaxis()-> SetTitle(
"Events");
519 h_ele_vertexP -> GetXaxis()-> SetTitle(
"p_{vertex} (GeV/c)");
520 h_ele_vertexP -> GetYaxis()-> SetTitle(
"Events");
521 h_ele_vertexPt -> GetXaxis()-> SetTitle(
"p_{T vertex} (GeV/c)");
522 h_ele_vertexPt -> GetYaxis()-> SetTitle(
"Events");
523 h_ele_vertexEta -> GetXaxis()-> SetTitle(
"#eta");
524 h_ele_vertexEta -> GetYaxis()-> SetTitle(
"Events");
525 h_ele_vertexPhi -> GetXaxis()-> SetTitle(
"#phi (rad)");
526 h_ele_vertexPhi -> GetYaxis()-> SetTitle(
"Events");
527 h_ele_PoPmatchingObject -> GetXaxis()-> SetTitle(
"P/P_{gen jet}");
528 h_ele_PoPmatchingObject -> GetYaxis()-> SetTitle(
"Events");
529 h_ele_PoPmatchingObject_barrel -> GetXaxis()-> SetTitle(
"P/P_{gen jet}");
530 h_ele_PoPmatchingObject_barrel -> GetYaxis()-> SetTitle(
"Events");
531 h_ele_PoPmatchingObject_endcaps -> GetXaxis()-> SetTitle(
"P/P_{gen jet}");
532 h_ele_PoPmatchingObject_endcaps -> GetYaxis()-> SetTitle(
"Events");
533 histSclSigIEtaIEta_-> GetXaxis()-> SetTitle(
"#sigma_{i#eta i#eta}") ;
534 histSclSigIEtaIEta_-> GetYaxis()-> SetTitle(
"Events") ;
535 histSclSigIEtaIEta_barrel_-> GetXaxis()-> SetTitle(
"#sigma_{i#eta i#eta}") ;
536 histSclSigIEtaIEta_barrel_-> GetYaxis()-> SetTitle(
"Events") ;
537 histSclSigIEtaIEta_endcaps_-> GetXaxis()-> SetTitle(
"#sigma_{i#eta i#eta}") ;
538 histSclSigIEtaIEta_endcaps_-> GetYaxis()-> SetTitle(
"Events") ;
539 histSclE1x5_-> GetXaxis()-> SetTitle(
"E1x5 (GeV)") ;
540 histSclE1x5_-> GetYaxis()-> SetTitle(
"Events") ;
541 histSclE1x5_barrel_-> GetXaxis()-> SetTitle(
"E1x5 (GeV)") ;
542 histSclE1x5_barrel_-> GetYaxis()-> SetTitle(
"Events") ;
543 histSclE1x5_endcaps_-> GetXaxis()-> SetTitle(
"E1x5 (GeV)") ;
544 histSclE1x5_endcaps_-> GetYaxis()-> SetTitle(
"Events") ;
545 histSclE2x5max_-> GetXaxis()-> SetTitle(
"E2x5 (GeV)") ;
546 histSclE2x5max_-> GetYaxis()-> SetTitle(
"Events") ;
547 histSclE2x5max_barrel_-> GetXaxis()-> SetTitle(
"E2x5 (GeV)") ;
548 histSclE2x5max_barrel_-> GetYaxis()-> SetTitle(
"Events") ;
549 histSclE2x5max_endcaps_-> GetXaxis()-> SetTitle(
"E2x5 (GeV)") ;
550 histSclE2x5max_endcaps_-> GetYaxis()-> SetTitle(
"Events") ;
551 histSclE5x5_-> GetXaxis()-> SetTitle(
"E5x5 (GeV)") ;
552 histSclE5x5_-> GetYaxis()-> SetTitle(
"Events") ;
553 histSclE5x5_barrel_-> GetXaxis()-> SetTitle(
"E5x5 (GeV)") ;
554 histSclE5x5_barrel_-> GetYaxis()-> SetTitle(
"Events") ;
555 histSclE5x5_endcaps_-> GetXaxis()-> SetTitle(
"E5x5 (GeV)") ;
556 histSclE5x5_endcaps_-> GetYaxis()-> SetTitle(
"Events") ;
557 h_ele_EtaMnEtamatchingObject -> GetXaxis()-> SetTitle(
"#eta_{rec} - #eta_{jet}");
558 h_ele_EtaMnEtamatchingObject -> GetYaxis()-> SetTitle(
"Events");
559 h_ele_PhiMnPhimatchingObject -> GetXaxis()-> SetTitle(
"#phi_{rec} - #phi_{jet} (rad)");
560 h_ele_PhiMnPhimatchingObject -> GetYaxis()-> SetTitle(
"Events");
561 h_ele_PinMnPout -> GetXaxis()-> SetTitle(
"P_{vertex} - P_{out} (GeV/c)");
562 h_ele_PinMnPout -> GetYaxis()-> SetTitle(
"Events");
563 h_ele_PinMnPout_mode -> GetXaxis()-> SetTitle(
"P_{vertex} - P_{out}, mode (GeV/c)");
564 h_ele_PinMnPout_mode -> GetYaxis()-> SetTitle(
"Events");
565 h_ele_outerP -> GetXaxis()-> SetTitle(
"P_{out} (GeV/c)");
566 h_ele_outerP -> GetYaxis()-> SetTitle(
"Events");
567 h_ele_outerP_mode -> GetXaxis()-> SetTitle(
"P_{out} (GeV/c)");
568 h_ele_outerP_mode -> GetYaxis()-> SetTitle(
"Events");
569 h_ele_outerPt -> GetXaxis()-> SetTitle(
"P_{T out} (GeV/c)");
570 h_ele_outerPt -> GetYaxis()-> SetTitle(
"Events");
571 h_ele_outerPt_mode -> GetXaxis()-> SetTitle(
"P_{T out} (GeV/c)");
572 h_ele_outerPt_mode -> GetYaxis()-> SetTitle(
"Events");
573 h_ele_EoP -> GetXaxis()-> SetTitle(
"E/P_{vertex}");
574 h_ele_EoP -> GetYaxis()-> SetTitle(
"Events");
575 h_ele_EseedOP -> GetXaxis()-> SetTitle(
"E_{seed}/P_{vertex}");
576 h_ele_EseedOP -> GetYaxis()-> SetTitle(
"Events");
577 h_ele_EoPout -> GetXaxis()-> SetTitle(
"E_{seed}/P_{out}");
578 h_ele_EoPout -> GetYaxis()-> SetTitle(
"Events");
579 h_ele_EeleOPout -> GetXaxis()-> SetTitle(
"E_{ele}/P_{out}");
580 h_ele_EeleOPout -> GetYaxis()-> SetTitle(
"Events");
581 h_ele_EoP_barrel -> GetXaxis()-> SetTitle(
"E/P_{vertex}");
582 h_ele_EoP_barrel -> GetYaxis()-> SetTitle(
"Events");
583 h_ele_EseedOP_barrel -> GetXaxis()-> SetTitle(
"E_{seed}/P_{vertex}");
584 h_ele_EseedOP_barrel -> GetYaxis()-> SetTitle(
"Events");
585 h_ele_EoPout_barrel -> GetXaxis()-> SetTitle(
"E_{seed}/P_{out}");
586 h_ele_EoPout_barrel -> GetYaxis()-> SetTitle(
"Events");
587 h_ele_EeleOPout_barrel -> GetXaxis()-> SetTitle(
"E_{ele}/P_{out}");
588 h_ele_EeleOPout_barrel -> GetYaxis()-> SetTitle(
"Events");
589 h_ele_EoP_endcaps -> GetXaxis()-> SetTitle(
"E/P_{vertex}");
590 h_ele_EoP_endcaps -> GetYaxis()-> SetTitle(
"Events");
591 h_ele_EseedOP_endcaps -> GetXaxis()-> SetTitle(
"E_{seed}/P_{vertex}");
592 h_ele_EseedOP_endcaps -> GetYaxis()-> SetTitle(
"Events");
593 h_ele_EoPout_endcaps -> GetXaxis()-> SetTitle(
"E_{seed}/P_{out}");
594 h_ele_EoPout_endcaps -> GetYaxis()-> SetTitle(
"Events");
595 h_ele_EeleOPout_endcaps -> GetXaxis()-> SetTitle(
"E_{ele}/P_{out}");
596 h_ele_EeleOPout_endcaps -> GetYaxis()-> SetTitle(
"Events");
597 h_ele_vertexX-> GetXaxis()-> SetTitle(
"x (cm)");
598 h_ele_vertexX-> GetYaxis()-> SetTitle(
"Events");
599 h_ele_vertexY-> GetXaxis()-> SetTitle(
"y (cm)");
600 h_ele_vertexY-> GetYaxis()-> SetTitle(
"Events");
601 h_ele_vertexZ-> GetXaxis()-> SetTitle(
"z (cm)");
602 h_ele_vertexZ-> GetYaxis()-> SetTitle(
"Events");
603 h_ele_vertexTIP-> GetXaxis()-> SetTitle(
"TIP (cm)");
604 h_ele_vertexTIP-> GetYaxis()-> SetTitle(
"Events");
605 h_ele_TIP_all-> GetXaxis()-> SetTitle(
"r_{T} (cm)");
606 h_ele_TIP_all-> GetYaxis()-> SetTitle(
"Events");
607 h_ele_vertexTIPVsEta-> GetYaxis()-> SetTitle(
"TIP (cm)");
608 h_ele_vertexTIPVsEta-> GetXaxis()-> SetTitle(
"#eta");
609 h_ele_vertexTIPVsPhi-> GetYaxis()-> SetTitle(
"TIP (cm)");
610 h_ele_vertexTIPVsPhi-> GetXaxis()-> SetTitle(
"#phi (rad)");
611 h_ele_vertexTIPVsPt-> GetYaxis()-> SetTitle(
"TIP (cm)");
612 h_ele_vertexTIPVsPt-> GetXaxis()-> SetTitle(
"p_{T} (GeV/c)");
613 h_ele_dEtaSc_propVtx-> GetXaxis()-> SetTitle(
"#eta_{sc} - #eta_{tr}");
614 h_ele_dEtaSc_propVtx-> GetXaxis()-> SetTitle(
"#eta_{sc} - #eta_{tr}");
615 h_ele_dEtaSc_propVtx-> GetYaxis()-> SetTitle(
"Events");
616 h_ele_dEtaCl_propOut-> GetXaxis()-> SetTitle(
"#eta_{seedcl} - #eta_{tr}");
617 h_ele_dEtaCl_propOut-> GetYaxis()-> SetTitle(
"Events");
618 h_ele_dEtaEleCl_propOut-> GetXaxis()-> SetTitle(
"#eta_{elecl} - #eta_{tr}");
619 h_ele_dEtaEleCl_propOut-> GetYaxis()-> SetTitle(
"Events");
620 h_ele_dPhiSc_propVtx-> GetXaxis()-> SetTitle(
"#phi_{sc} - #phi_{tr} (rad)");
621 h_ele_dPhiSc_propVtx-> GetYaxis()-> SetTitle(
"Events");
622 h_ele_dPhiCl_propOut-> GetXaxis()-> SetTitle(
"#phi_{seedcl} - #phi_{tr} (rad)");
623 h_ele_dPhiCl_propOut-> GetYaxis()-> SetTitle(
"Events");
624 h_ele_dPhiEleCl_propOut-> GetXaxis()-> SetTitle(
"#phi_{elecl} - #phi_{tr} (rad)");
625 h_ele_dPhiEleCl_propOut-> GetYaxis()-> SetTitle(
"Events");
626 h_ele_dEtaSc_propVtx_barrel-> GetXaxis()-> SetTitle(
"#eta_{sc} - #eta_{tr}");
627 h_ele_dEtaSc_propVtx_barrel-> GetYaxis()-> SetTitle(
"Events");
628 h_ele_dEtaCl_propOut_barrel-> GetXaxis()-> SetTitle(
"#eta_{seedcl} - #eta_{tr}");
629 h_ele_dEtaCl_propOut_barrel-> GetYaxis()-> SetTitle(
"Events");
630 h_ele_dEtaEleCl_propOut_barrel-> GetXaxis()-> SetTitle(
"#eta_{elecl} - #eta_{tr}");
631 h_ele_dEtaEleCl_propOut_barrel-> GetYaxis()-> SetTitle(
"Events");
632 h_ele_dPhiSc_propVtx_barrel-> GetXaxis()-> SetTitle(
"#phi_{sc} - #phi_{tr} (rad)");
633 h_ele_dPhiSc_propVtx_barrel-> GetYaxis()-> SetTitle(
"Events");
634 h_ele_dPhiCl_propOut_barrel-> GetXaxis()-> SetTitle(
"#phi_{seedcl} - #phi_{tr} (rad)");
635 h_ele_dPhiCl_propOut_barrel-> GetYaxis()-> SetTitle(
"Events");
636 h_ele_dPhiEleCl_propOut_barrel-> GetXaxis()-> SetTitle(
"#phi_{elecl} - #phi_{tr} (rad)");
637 h_ele_dPhiEleCl_propOut_barrel-> GetYaxis()-> SetTitle(
"Events");
638 h_ele_dEtaSc_propVtx_endcaps-> GetXaxis()-> SetTitle(
"#eta_{sc} - #eta_{tr}");
639 h_ele_dEtaSc_propVtx_endcaps-> GetYaxis()-> SetTitle(
"Events");
640 h_ele_dEtaCl_propOut_endcaps-> GetXaxis()-> SetTitle(
"#eta_{seedcl} - #eta_{tr}");
641 h_ele_dEtaCl_propOut_endcaps-> GetYaxis()-> SetTitle(
"Events");
642 h_ele_dEtaEleCl_propOut_endcaps-> GetXaxis()-> SetTitle(
"#eta_{elecl} - #eta_{tr}");
643 h_ele_dEtaEleCl_propOut_endcaps-> GetYaxis()-> SetTitle(
"Events");
644 h_ele_dPhiSc_propVtx_endcaps-> GetXaxis()-> SetTitle(
"#phi_{sc} - #phi_{tr} (rad)");
645 h_ele_dPhiSc_propVtx_endcaps-> GetYaxis()-> SetTitle(
"Events");
646 h_ele_dPhiCl_propOut_endcaps-> GetXaxis()-> SetTitle(
"#phi_{seedcl} - #phi_{tr} (rad)");
647 h_ele_dPhiCl_propOut_endcaps-> GetYaxis()-> SetTitle(
"Events");
648 h_ele_dPhiEleCl_propOut_endcaps-> GetXaxis()-> SetTitle(
"#phi_{elecl} - #phi_{tr} (rad)");
649 h_ele_dPhiEleCl_propOut_endcaps-> GetYaxis()-> SetTitle(
"Events");
650 h_ele_HoE-> GetXaxis()-> SetTitle(
"H/E") ;
651 h_ele_HoE-> GetYaxis()-> SetTitle(
"Events") ;
652 h_ele_HoE_barrel-> GetXaxis()-> SetTitle(
"H/E") ;
653 h_ele_HoE_barrel-> GetYaxis()-> SetTitle(
"Events") ;
654 h_ele_HoE_endcaps-> GetXaxis()-> SetTitle(
"H/E") ;
655 h_ele_HoE_endcaps-> GetYaxis()-> SetTitle(
"Events") ;
656 h_ele_HoE_fiducial-> GetXaxis()-> SetTitle(
"H/E") ;
657 h_ele_HoE_fiducial-> GetYaxis()-> SetTitle(
"Events") ;
658 h_ele_fbrem-> GetXaxis()-> SetTitle(
"P_{in} - P_{out} / P_{in}");
659 h_ele_fbrem-> GetYaxis()-> SetTitle(
"Events");
660 h_ele_seed_dphi2_-> GetXaxis()-> SetTitle(
"#phi_{hit}-#phi_{pred} (rad)") ;
661 h_ele_seed_dphi2_-> GetYaxis()-> SetTitle(
"Events") ;
662 h_ele_seed_drz2_-> GetXaxis()-> SetTitle(
"r(z)_{hit}-r(z)_{pred} (cm)") ;
663 h_ele_seed_drz2_-> GetYaxis()-> SetTitle(
"Events") ;
664 h_ele_seed_subdet2_-> GetXaxis()-> SetTitle(
"2nd hit subdet Id") ;
665 h_ele_seed_subdet2_-> GetYaxis()-> SetTitle(
"Events") ;
666 h_ele_classes-> GetXaxis()-> SetTitle(
"class Id") ;
667 h_ele_classes-> GetYaxis()-> SetTitle(
"Events") ;
668 h_ele_EoverP_all-> GetXaxis()-> SetTitle(
"E/P_{vertex}");
669 h_ele_EoverP_all-> GetYaxis()-> SetTitle(
"Events");
670 h_ele_EseedOP_all-> GetXaxis()-> SetTitle(
"E_{seed}/P_{vertex}");
671 h_ele_EseedOP_all-> GetYaxis()-> SetTitle(
"Events");
672 h_ele_EoPout_all -> GetXaxis()-> SetTitle(
"E_{seed}/P_{out}");
673 h_ele_EoPout_all-> GetYaxis()-> SetTitle(
"Events");
674 h_ele_EeleOPout_all-> GetXaxis()-> SetTitle(
"E_{ele}/P_{out}");
675 h_ele_EeleOPout_all-> GetYaxis()-> SetTitle(
"Events");
676 h_ele_dEtaSc_propVtx_all-> GetXaxis()-> SetTitle(
"#eta_{sc} - #eta_{tr}");
677 h_ele_dEtaSc_propVtx_all-> GetYaxis()-> SetTitle(
"Events");
678 h_ele_dPhiSc_propVtx_all-> GetXaxis()-> SetTitle(
"#phi_{sc} - #phi_{tr} (rad)");
679 h_ele_dPhiSc_propVtx_all-> GetYaxis()-> SetTitle(
"Events");
680 h_ele_dEtaCl_propOut_all-> GetXaxis()-> SetTitle(
"#eta_{sc} - #eta_{tr}");
681 h_ele_dEtaCl_propOut_all-> GetYaxis()-> SetTitle(
"Events");
682 h_ele_dPhiCl_propOut_all-> GetXaxis()-> SetTitle(
"#phi_{sc} - #phi_{tr} (rad)");
683 h_ele_dPhiCl_propOut_all-> GetYaxis()-> SetTitle(
"Events");
684 h_ele_HoE_all-> GetXaxis()-> SetTitle(
"H/E") ;
685 h_ele_HoE_all-> GetYaxis()-> SetTitle(
"Events");
686 h_ele_mee_all-> GetXaxis()-> SetTitle(
"m_{ee} (GeV/c^{2})");
687 h_ele_mee_all-> GetYaxis()-> SetTitle(
"Events");
688 h_ele_mee_os-> GetXaxis()-> SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
689 h_ele_mee_os-> GetYaxis()-> SetTitle(
"Events");
690 h_ele_E2mnE1vsMee_all-> GetXaxis()-> SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
691 h_ele_E2mnE1vsMee_all-> GetYaxis()-> SetTitle(
"E2 - E1 (GeV)");
692 h_ele_E2mnE1vsMee_egeg_all-> GetXaxis()-> SetTitle(
"m_{e^{+}e^{-}} (GeV/c^{2})");
693 h_ele_E2mnE1vsMee_egeg_all-> GetYaxis()-> SetTitle(
"E2 - E1 (GeV)");
694 histNum_-> GetXaxis()-> SetTitle(
"N_{ele}");
695 histNum_-> GetYaxis()-> SetTitle(
"Events");
696 h_ele_fbremVsEta_mode-> GetXaxis()-> SetTitle(
"#eta");
697 h_ele_fbremVsEta_mean-> GetXaxis()-> SetTitle(
"#eta");
705 std::cout <<
"efficiency calculation " << std::endl;
707 TH1F *h_ele_etaEff = (TH1F*)h_ele_matchingObjectEta_matched->Clone(
"h_ele_etaEff");
708 h_ele_etaEff->Reset();
709 h_ele_etaEff->Divide(h_ele_matchingObjectEta_matched,h_matchingObjectEta,1,1,
"b");
710 h_ele_etaEff->Print();
711 h_ele_etaEff->GetXaxis()->SetTitle(
"#eta");
712 h_ele_etaEff->GetYaxis()->SetTitle(
"Efficiency");
715 TH1F *h_ele_zEff = (TH1F*)h_ele_matchingObjectZ_matched->Clone(
"h_ele_zEff");
717 h_ele_zEff->Divide(h_ele_matchingObjectZ_matched,h_matchingObjectZ,1,1,
"b");
719 h_ele_zEff->GetXaxis()->SetTitle(
"z (cm)");
720 h_ele_zEff->GetYaxis()->SetTitle(
"Efficiency");
723 TH1F *h_ele_absetaEff = (TH1F*)h_ele_matchingObjectAbsEta_matched->Clone(
"h_ele_absetaEff");
724 h_ele_absetaEff->Reset();
725 h_ele_absetaEff->Divide(h_ele_matchingObjectAbsEta_matched,h_matchingObjectAbsEta,1,1,
"b");
726 h_ele_absetaEff->GetXaxis()->SetTitle(
"|#eta|");
727 h_ele_absetaEff->GetYaxis()->SetTitle(
"Efficiency");
730 TH1F *h_ele_ptEff = (TH1F*)h_ele_matchingObjectPt_matched->Clone(
"h_ele_ptEff");
731 h_ele_ptEff->Reset();
732 h_ele_ptEff->Divide(h_ele_matchingObjectPt_matched,h_matchingObjectPt,1,1,
"b");
733 h_ele_ptEff->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
734 h_ele_ptEff->GetYaxis()->SetTitle(
"Efficiency");
737 TH1F *h_ele_phiEff = (TH1F*)h_ele_matchingObjectPhi_matched->Clone(
"h_ele_phiEff");
738 h_ele_phiEff->Reset();
739 h_ele_phiEff->Divide(h_ele_matchingObjectPhi_matched,h_matchingObjectPhi,1,1,
"b");
740 h_ele_phiEff->GetXaxis()->SetTitle(
"#phi (rad)");
741 h_ele_phiEff->GetYaxis()->SetTitle(
"Efficiency");
744 TH1F *h_ele_etaEff_all = (TH1F*)h_ele_vertexEta_all->Clone(
"h_ele_etaEff_all");
745 h_ele_etaEff_all->Reset();
746 h_ele_etaEff_all->Divide(h_ele_vertexEta_all,h_matchingObjectEta,1,1,
"b");
747 h_ele_etaEff_all->Print();
748 h_ele_etaEff_all->GetXaxis()->SetTitle(
"#eta");
749 h_ele_etaEff_all->GetYaxis()->SetTitle(
"N_{rec}/N_{matching jet}");
751 TH1F *h_ele_ptEff_all = (TH1F*)h_ele_vertexPt_all->Clone(
"h_ele_ptEff_all");
752 h_ele_ptEff_all->Reset();
753 h_ele_ptEff_all->Divide(h_ele_vertexPt_all,h_matchingObjectPt,1,1,
"b");
754 h_ele_ptEff_all->Print();
755 h_ele_ptEff_all->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
756 h_ele_ptEff_all->GetYaxis()->SetTitle(
"N_{rec}/N_{matching gen jet}");
759 TH1F *h_ele_eta_goldenFrac = (TH1F*)h_ele_eta_golden->Clone(
"h_ele_eta_goldenFrac");
760 h_ele_eta_goldenFrac->Reset();
761 h_ele_eta_goldenFrac->GetXaxis()->SetTitle(
"|#eta|");
762 h_ele_eta_goldenFrac->GetYaxis()->SetTitle(
"Fraction of electrons");
763 h_ele_eta_goldenFrac->Divide(h_ele_eta_golden,h_ele_eta,1,1);
764 h_ele_eta_goldenFrac->SetTitle(
"fraction of golden electrons vs eta");
765 TH1F *h_ele_eta_bbremFrac = (TH1F*)h_ele_eta_bbrem->Clone(
"h_ele_eta_bbremFrac");
766 h_ele_eta_bbremFrac->Reset();
767 h_ele_eta_bbremFrac->Divide(h_ele_eta_bbrem,h_ele_eta,1,1);
768 h_ele_eta_bbremFrac->GetXaxis()->SetTitle(
"|#eta|");
769 h_ele_eta_bbremFrac->GetYaxis()->SetTitle(
"Fraction of electrons");
770 h_ele_eta_bbremFrac->SetTitle(
"fraction of big brem electrons vs eta");
771 TH1F *h_ele_eta_narrowFrac = (TH1F*)h_ele_eta_narrow->Clone(
"h_ele_eta_narrowFrac");
772 h_ele_eta_narrowFrac->Reset();
773 h_ele_eta_narrowFrac->Divide(h_ele_eta_narrow,h_ele_eta,1,1);
774 h_ele_eta_narrowFrac->GetXaxis()->SetTitle(
"|#eta|");
775 h_ele_eta_narrowFrac->GetYaxis()->SetTitle(
"Fraction of electrons");
776 h_ele_eta_narrowFrac->SetTitle(
"fraction of narrow electrons vs eta");
777 TH1F *h_ele_eta_showerFrac = (TH1F*)h_ele_eta_shower->Clone(
"h_ele_eta_showerFrac");
778 h_ele_eta_showerFrac->Reset();
779 h_ele_eta_showerFrac->Divide(h_ele_eta_shower,h_ele_eta,1,1);
780 h_ele_eta_showerFrac->GetXaxis()->SetTitle(
"|#eta|");
781 h_ele_eta_showerFrac->GetYaxis()->SetTitle(
"Fraction of electrons");
782 h_ele_eta_showerFrac->SetTitle(
"fraction of showering electrons vs eta");
785 TH1F *h_ele_xOverX0VsEta =
new TH1F(
"h_ele_xOverx0VsEta",
"mean X/X_0 vs eta",nbineta/2,0.0,2.5);
786 for (
int ibin=1;ibin<h_ele_fbremVsEta_mean->GetNbinsX()+1;ibin++) {
788 if (h_ele_fbremVsEta_mean->GetBinContent(ibin)>0.) xOverX0 = -
log(h_ele_fbremVsEta_mean->GetBinContent(ibin));
789 h_ele_xOverX0VsEta->SetBinContent(ibin,xOverX0);
793 TProfile *p_ele_PoPmatchingObjectVsEta = h_ele_PoPmatchingObjectVsEta->ProfileX();
794 p_ele_PoPmatchingObjectVsEta->SetTitle(
"mean ele momentum / matching SC energy vs eta");
795 p_ele_PoPmatchingObjectVsEta->GetXaxis()->SetTitle(
"#eta");
796 p_ele_PoPmatchingObjectVsEta->GetYaxis()->SetTitle(
"<P/E_{matching gen jet}>");
797 p_ele_PoPmatchingObjectVsEta->Write();
798 TProfile *p_ele_PoPmatchingObjectVsPhi = h_ele_PoPmatchingObjectVsPhi->ProfileX();
799 p_ele_PoPmatchingObjectVsPhi->SetTitle(
"mean ele momentum / gen momentum vs phi");
800 p_ele_PoPmatchingObjectVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
801 p_ele_PoPmatchingObjectVsPhi->GetYaxis()->SetTitle(
"<P/E_{matching gen jet}>");
802 p_ele_PoPmatchingObjectVsPhi->Write();
803 TProfile *p_ele_EtaMnEtamatchingObjectVsEta = h_ele_EtaMnEtamatchingObjectVsEta->ProfileX();
804 p_ele_EtaMnEtamatchingObjectVsEta->GetXaxis()->SetTitle(
"#eta");
805 p_ele_EtaMnEtamatchingObjectVsEta->GetYaxis()->SetTitle(
"<#eta_{rec} - #eta_{matching gen jet}>");
806 p_ele_EtaMnEtamatchingObjectVsEta->Write();
807 TProfile *p_ele_EtaMnEtamatchingObjectVsPhi = h_ele_EtaMnEtamatchingObjectVsPhi->ProfileX();
808 p_ele_EtaMnEtamatchingObjectVsPhi-> GetXaxis()-> SetTitle(
"#phi");
809 p_ele_EtaMnEtamatchingObjectVsPhi-> GetYaxis()-> SetTitle(
"<#eta_{rec} - #eta_{matching gen jet}>");
810 p_ele_EtaMnEtamatchingObjectVsPhi->Write();
811 TProfile *p_ele_PhiMnPhimatchingObjectVsEta = h_ele_PhiMnPhimatchingObjectVsEta->ProfileX();
812 p_ele_PhiMnPhimatchingObjectVsEta-> GetXaxis()-> SetTitle(
"#eta");
813 p_ele_PhiMnPhimatchingObjectVsEta-> GetYaxis()-> SetTitle(
"<#phi_{rec} - #phi_{matching gen jet}> (rad)");
814 p_ele_PhiMnPhimatchingObjectVsEta->Write();
815 TProfile *p_ele_PhiMnPhimatchingObjectVsPhi = h_ele_PhiMnPhimatchingObjectVsPhi->ProfileX();
816 p_ele_PhiMnPhimatchingObjectVsPhi-> GetXaxis()-> SetTitle(
"#phi");
817 p_ele_PhiMnPhimatchingObjectVsPhi-> GetYaxis()-> SetTitle(
"<#phi_{rec} - #phi_{matching gen jet}> (rad)");
818 p_ele_PhiMnPhimatchingObjectVsPhi->Write();
819 TProfile *p_ele_vertexPtVsEta = h_ele_vertexPtVsEta->ProfileX();
820 p_ele_vertexPtVsEta->GetXaxis()->SetTitle(
"#eta");
821 p_ele_vertexPtVsEta->GetYaxis()->SetTitle(
"<p_{T}> (GeV/c)");
822 p_ele_vertexPtVsEta->Write();
823 TProfile *p_ele_vertexPtVsPhi = h_ele_vertexPtVsPhi->ProfileX();
824 p_ele_vertexPtVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
825 p_ele_vertexPtVsPhi->GetYaxis()->SetTitle(
"<p_{T}> (GeV/c)");
826 p_ele_vertexPtVsPhi->Write();
827 TProfile *p_ele_EoPVsEta = h_ele_EoPVsEta->ProfileX();
828 p_ele_EoPVsEta->GetXaxis()->SetTitle(
"#eta");
829 p_ele_EoPVsEta->GetYaxis()->SetTitle(
"<E/P_{vertex}>");
830 p_ele_EoPVsEta->Write();
831 TProfile *p_ele_EoPVsPhi = h_ele_EoPVsPhi->ProfileX();
832 p_ele_EoPVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
833 p_ele_EoPVsPhi->GetYaxis()->SetTitle(
"<E/P_{vertex}>");
834 p_ele_EoPVsPhi->Write();
835 TProfile *p_ele_EoPoutVsEta = h_ele_EoPoutVsEta->ProfileX();
836 p_ele_EoPoutVsEta->GetXaxis()->SetTitle(
"#eta");
837 p_ele_EoPoutVsEta->GetYaxis()->SetTitle(
"<E_{seed}/P_{out}>");
838 p_ele_EoPoutVsEta->Write();
839 TProfile *p_ele_EoPoutVsPhi = h_ele_EoPoutVsPhi->ProfileX();
840 p_ele_EoPoutVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
841 p_ele_EoPoutVsPhi->GetYaxis()->SetTitle(
"<E_{seed}/P_{out}>");
842 p_ele_EoPoutVsPhi->Write();
843 TProfile *p_ele_EeleOPoutVsEta = h_ele_EeleOPoutVsEta->ProfileX();
844 p_ele_EeleOPoutVsEta->SetTitle(
"mean ele Eele/pout vs eta");
845 p_ele_EeleOPoutVsEta->GetXaxis()->SetTitle(
"#eta");
846 p_ele_EeleOPoutVsEta->GetYaxis()->SetTitle(
"<E_{ele}/P_{out}>");
847 p_ele_EeleOPoutVsEta->Write();
848 TProfile *p_ele_EeleOPoutVsPhi = h_ele_EeleOPoutVsPhi->ProfileX();
849 p_ele_EeleOPoutVsPhi->SetTitle(
"mean ele Eele/pout vs phi");
850 p_ele_EeleOPoutVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
851 p_ele_EeleOPoutVsPhi->GetYaxis()->SetTitle(
"<E_{ele}/P_{out}>");
852 p_ele_EeleOPoutVsPhi->Write();
853 TProfile *p_ele_HoEVsEta = h_ele_HoEVsEta->ProfileX();
854 p_ele_HoEVsEta->GetXaxis()->SetTitle(
"#eta");
855 p_ele_HoEVsEta->GetYaxis()->SetTitle(
"<H/E>");
856 p_ele_HoEVsEta->Write();
857 TProfile *p_ele_HoEVsPhi = h_ele_HoEVsPhi->ProfileX();
858 p_ele_HoEVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
859 p_ele_HoEVsPhi->GetYaxis()->SetTitle(
"<H/E>");
860 p_ele_HoEVsPhi->Write();
861 TProfile *p_ele_chi2VsEta = h_ele_chi2VsEta->ProfileX();
862 p_ele_chi2VsEta->GetXaxis()->SetTitle(
"#eta");
863 p_ele_chi2VsEta->GetYaxis()->SetTitle(
"<#Chi^{2}>");
864 p_ele_chi2VsEta->Write();
865 TProfile *p_ele_chi2VsPhi = h_ele_chi2VsPhi->ProfileX();
866 p_ele_chi2VsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
867 p_ele_chi2VsPhi->GetYaxis()->SetTitle(
"<#Chi^{2}>");
868 p_ele_chi2VsPhi->Write();
869 TProfile *p_ele_foundHitsVsEta = h_ele_foundHitsVsEta->ProfileX();
870 p_ele_foundHitsVsEta->GetXaxis()->SetTitle(
"#eta");
871 p_ele_foundHitsVsEta->GetYaxis()->SetTitle(
"<N_{hits}>");
872 p_ele_foundHitsVsEta->Write();
873 TProfile *p_ele_foundHitsVsPhi = h_ele_foundHitsVsPhi->ProfileX();
874 p_ele_foundHitsVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
875 p_ele_foundHitsVsPhi->GetYaxis()->SetTitle(
"<N_{hits}>");
876 p_ele_foundHitsVsPhi->Write();
877 TProfile *p_ele_lostHitsVsEta = h_ele_lostHitsVsEta->ProfileX();
878 p_ele_lostHitsVsEta->GetXaxis()->SetTitle(
"#eta");
879 p_ele_lostHitsVsEta->GetYaxis()->SetTitle(
"<N_{hits}>");
880 p_ele_lostHitsVsEta->Write();
881 TProfile *p_ele_lostHitsVsPhi = h_ele_lostHitsVsPhi->ProfileX();
882 p_ele_lostHitsVsPhi->GetXaxis()->SetTitle(
"#phi (rad)");
883 p_ele_lostHitsVsPhi->GetYaxis()->SetTitle(
"<N_{hits}>");
884 p_ele_lostHitsVsPhi->Write();
885 TProfile *p_ele_vertexTIPVsEta = h_ele_vertexTIPVsEta->ProfileX();
886 p_ele_vertexTIPVsEta->SetTitle(
"mean tip (wrt gen vtx) vs eta");
887 p_ele_vertexTIPVsEta->GetXaxis()->SetTitle(
"#eta");
888 p_ele_vertexTIPVsEta->GetYaxis()->SetTitle(
"<TIP> (cm)");
889 p_ele_vertexTIPVsEta->Write();
890 TProfile *p_ele_vertexTIPVsPhi = h_ele_vertexTIPVsPhi->ProfileX();
891 p_ele_vertexTIPVsPhi->SetTitle(
"mean tip (wrt gen vtx) vs phi");
892 p_ele_vertexTIPVsPhi->GetXaxis()->SetTitle(
"#phi");
893 p_ele_vertexTIPVsPhi->GetYaxis()->SetTitle(
"<TIP> (cm)");
894 p_ele_vertexTIPVsPhi->Write();
895 TProfile *p_ele_vertexTIPVsPt = h_ele_vertexTIPVsPt->ProfileX();
896 p_ele_vertexTIPVsPt->SetTitle(
"mean tip (wrt gen vtx) vs phi");
897 p_ele_vertexTIPVsPt->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
898 p_ele_vertexTIPVsPt->GetYaxis()->SetTitle(
"<TIP> (cm)");
899 p_ele_vertexTIPVsPt->Write();
903 h_matchingObjectNum->Write();
910 h_matchingObjectEta->Write();
911 h_matchingObjectAbsEta->Write();
912 h_matchingObjectP->Write();
913 h_matchingObjectPt->Write();
914 h_matchingObjectPhi->Write();
915 h_matchingObjectZ->Write();
918 h_ele_EoverP_all->Write();
919 h_ele_EseedOP_all->Write();
920 h_ele_EoPout_all->Write();
921 h_ele_EeleOPout_all ->Write();
922 h_ele_dEtaSc_propVtx_all->Write();
923 h_ele_dPhiSc_propVtx_all->Write();
924 h_ele_dEtaCl_propOut_all ->Write();
925 h_ele_dPhiCl_propOut_all->Write();
926 h_ele_HoE_all->Write();
927 h_ele_TIP_all->Write();
928 h_ele_vertexPt_all->Write();
929 h_ele_vertexEta_all->Write();
930 h_ele_mee_all->Write();
931 h_ele_mee_os->Write();
932 h_ele_E2mnE1vsMee_all ->Write();
933 h_ele_E2mnE1vsMee_egeg_all->Write();
936 h_ele_charge->Write();
937 h_ele_chargeVsEta->Write();
938 h_ele_chargeVsPhi->Write();
939 h_ele_chargeVsPt->Write();
940 h_ele_vertexP->Write();
941 h_ele_vertexPt->Write();
942 h_ele_vertexPtVsEta->Write();
943 h_ele_vertexPtVsPhi->Write();
944 h_ele_matchingObjectPt_matched->Write();
945 h_ele_vertexEta->Write();
946 h_ele_vertexEtaVsPhi->Write();
947 h_ele_matchingObjectAbsEta_matched->Write();
948 h_ele_matchingObjectEta_matched->Write();
949 h_ele_vertexPhi->Write();
950 h_ele_matchingObjectPhi_matched->Write();
951 h_ele_vertexX->Write();
952 h_ele_vertexY ->Write();
953 h_ele_vertexZ->Write();
954 h_ele_vertexTIP->Write();
955 h_ele_matchingObjectZ_matched->Write();
956 h_ele_vertexTIPVsEta->Write();
957 h_ele_vertexTIPVsPhi->Write();
958 h_ele_vertexTIPVsPt->Write();
959 h_ele_PoPmatchingObject->Write();
960 h_ele_PoPmatchingObjectVsEta ->Write();
961 h_ele_PoPmatchingObjectVsPhi->Write();
962 h_ele_PoPmatchingObjectVsPt->Write();
963 h_ele_PoPmatchingObject_barrel ->Write();
964 h_ele_PoPmatchingObject_endcaps->Write();
965 h_ele_EtaMnEtamatchingObject->Write();
966 h_ele_EtaMnEtamatchingObjectVsEta ->Write();
967 h_ele_EtaMnEtamatchingObjectVsPhi->Write();
968 h_ele_EtaMnEtamatchingObjectVsPt->Write();
969 h_ele_PhiMnPhimatchingObject ->Write();
970 h_ele_PhiMnPhimatchingObject2 ->Write();
971 h_ele_PhiMnPhimatchingObjectVsEta->Write();
972 h_ele_PhiMnPhimatchingObjectVsPhi->Write();
973 h_ele_PhiMnPhimatchingObjectVsPt->Write();
977 histSclEoEmatchingObject_barrel->Write();
978 histSclEoEmatchingObject_endcaps->Write();
980 histSclEtVsEta_->Write();
981 histSclEtVsPhi_->Write();
982 histSclEtaVsPhi_ ->Write();
983 histSclEta_->Write();
984 histSclPhi_->Write();
985 histSclSigIEtaIEta_->Write();
986 histSclSigIEtaIEta_barrel_->Write();
987 histSclSigIEtaIEta_endcaps_->Write();
988 histSclE1x5_->Write();
989 histSclE1x5_barrel_->Write();
990 histSclE1x5_endcaps_->Write();
991 histSclE2x5max_->Write();
992 histSclE2x5max_barrel_->Write();
993 histSclE2x5max_endcaps_->Write();
994 histSclE5x5_->Write();
995 histSclE5x5_barrel_->Write();
996 histSclE5x5_endcaps_->Write();
999 h_ele_ambiguousTracks->Write();
1000 h_ele_ambiguousTracksVsEta->Write();
1001 h_ele_ambiguousTracksVsPhi->Write();
1002 h_ele_ambiguousTracksVsPt->Write();
1003 h_ele_foundHits->Write();
1004 h_ele_foundHitsVsEta->Write();
1005 h_ele_foundHitsVsPhi->Write();
1006 h_ele_foundHitsVsPt->Write();
1007 h_ele_lostHits->Write();
1008 h_ele_lostHitsVsEta->Write();
1009 h_ele_lostHitsVsPhi->Write();
1010 h_ele_lostHitsVsPt->Write();
1011 h_ele_chi2 ->Write();
1012 h_ele_chi2VsEta ->Write();
1013 h_ele_chi2VsPhi ->Write();
1014 h_ele_chi2VsPt->Write();
1015 h_ele_PinMnPout->Write();
1016 h_ele_PinMnPout_mode->Write();
1017 h_ele_PinMnPoutVsEta_mode->Write();
1018 h_ele_PinMnPoutVsPhi_mode->Write();
1019 h_ele_PinMnPoutVsPt_mode->Write();
1020 h_ele_PinMnPoutVsE_mode->Write();
1021 h_ele_PinMnPoutVsChi2_mode->Write();
1022 h_ele_outerP ->Write();
1023 h_ele_outerP_mode->Write();
1024 h_ele_outerPVsEta_mode->Write();
1025 h_ele_outerPt->Write();
1026 h_ele_outerPt_mode ->Write();
1027 h_ele_outerPtVsEta_mode->Write();
1028 h_ele_outerPtVsPhi_mode->Write();
1029 h_ele_outerPtVsPt_mode->Write();
1032 h_ele_EoP ->Write();
1033 h_ele_EoP_eg ->Write();
1034 h_ele_EoP_barrel ->Write();
1035 h_ele_EoP_eg_barrel ->Write();
1036 h_ele_EoP_endcaps ->Write();
1037 h_ele_EoP_eg_endcaps ->Write();
1038 h_ele_EoPVsEta ->Write();
1039 h_ele_EoPVsPhi->Write();
1040 h_ele_EoPVsE->Write();
1041 h_ele_EseedOP ->Write();
1042 h_ele_EseedOP_eg ->Write();
1043 h_ele_EseedOP_barrel ->Write();
1044 h_ele_EseedOP_eg_barrel ->Write();
1045 h_ele_EseedOP_endcaps ->Write();
1046 h_ele_EseedOP_eg_endcaps ->Write();
1047 h_ele_EseedOPVsEta ->Write();
1048 h_ele_EseedOPVsPhi->Write();
1049 h_ele_EseedOPVsE->Write();
1050 h_ele_EoPout->Write();
1051 h_ele_EoPout_eg->Write();
1052 h_ele_EoPout_barrel->Write();
1053 h_ele_EoPout_eg_barrel->Write();
1054 h_ele_EoPout_endcaps->Write();
1055 h_ele_EoPout_eg_endcaps->Write();
1056 h_ele_EoPoutVsEta->Write();
1057 h_ele_EoPoutVsPhi->Write();
1058 h_ele_EoPoutVsE ->Write();
1059 h_ele_EeleOPout->Write();
1060 h_ele_EeleOPout_eg->Write();
1061 h_ele_EeleOPout_barrel->Write();
1062 h_ele_EeleOPout_eg_barrel->Write();
1063 h_ele_EeleOPout_endcaps->Write();
1064 h_ele_EeleOPout_eg_endcaps->Write();
1065 h_ele_EeleOPoutVsEta->Write();
1066 h_ele_EeleOPoutVsPhi->Write();
1067 h_ele_EeleOPoutVsE ->Write();
1068 h_ele_dEtaSc_propVtx->Write();
1069 h_ele_dEtaSc_propVtx_eg->Write();
1070 h_ele_dEtaSc_propVtx_barrel->Write();
1071 h_ele_dEtaSc_propVtx_eg_barrel->Write();
1072 h_ele_dEtaSc_propVtx_endcaps->Write();
1073 h_ele_dEtaSc_propVtx_eg_endcaps->Write();
1074 h_ele_dEtaScVsEta_propVtx->Write();
1075 h_ele_dEtaScVsPhi_propVtx->Write();
1076 h_ele_dEtaScVsPt_propVtx ->Write();
1077 h_ele_dPhiSc_propVtx->Write();
1078 h_ele_dPhiSc_propVtx_eg->Write();
1079 h_ele_dPhiSc_propVtx_barrel->Write();
1080 h_ele_dPhiSc_propVtx_eg_barrel->Write();
1081 h_ele_dPhiSc_propVtx_endcaps->Write();
1082 h_ele_dPhiSc_propVtx_eg_endcaps->Write();
1083 h_ele_dPhiScVsEta_propVtx ->Write();
1084 h_ele_dPhiScVsPhi_propVtx->Write();
1085 h_ele_dPhiScVsPt_propVtx->Write();
1086 h_ele_dEtaCl_propOut->Write();
1087 h_ele_dEtaCl_propOut_eg->Write();
1088 h_ele_dEtaCl_propOut_barrel->Write();
1089 h_ele_dEtaCl_propOut_eg_barrel->Write();
1090 h_ele_dEtaCl_propOut_endcaps->Write();
1091 h_ele_dEtaCl_propOut_eg_endcaps->Write();
1092 h_ele_dEtaClVsEta_propOut->Write();
1093 h_ele_dEtaClVsPhi_propOut->Write();
1094 h_ele_dEtaClVsPt_propOut->Write();
1095 h_ele_dPhiCl_propOut->Write();
1096 h_ele_dPhiCl_propOut_eg->Write();
1097 h_ele_dPhiCl_propOut_barrel->Write();
1098 h_ele_dPhiCl_propOut_eg_barrel->Write();
1099 h_ele_dPhiCl_propOut_endcaps->Write();
1100 h_ele_dPhiCl_propOut_eg_endcaps->Write();
1101 h_ele_dPhiClVsEta_propOut->Write();
1102 h_ele_dPhiClVsPhi_propOut->Write();
1103 h_ele_dPhiClVsPt_propOut->Write();
1104 h_ele_dEtaEleCl_propOut->Write();
1105 h_ele_dEtaEleCl_propOut_eg->Write();
1106 h_ele_dEtaEleCl_propOut_barrel->Write();
1107 h_ele_dEtaEleCl_propOut_eg_barrel->Write();
1108 h_ele_dEtaEleCl_propOut_endcaps->Write();
1109 h_ele_dEtaEleCl_propOut_eg_endcaps->Write();
1110 h_ele_dEtaEleClVsEta_propOut->Write();
1111 h_ele_dEtaEleClVsPhi_propOut->Write();
1112 h_ele_dEtaEleClVsPt_propOut->Write();
1113 h_ele_dPhiEleCl_propOut->Write();
1114 h_ele_dPhiEleCl_propOut_eg->Write();
1115 h_ele_dPhiEleCl_propOut_barrel->Write();
1116 h_ele_dPhiEleCl_propOut_eg_barrel->Write();
1117 h_ele_dPhiEleCl_propOut_endcaps->Write();
1118 h_ele_dPhiEleCl_propOut_eg_endcaps->Write();
1119 h_ele_dPhiEleClVsEta_propOut->Write();
1120 h_ele_dPhiEleClVsPhi_propOut->Write();
1121 h_ele_dPhiEleClVsPt_propOut->Write();
1123 h_ele_HoE_eg->Write();
1124 h_ele_HoE_barrel->Write();
1125 h_ele_HoE_eg_barrel->Write();
1126 h_ele_HoE_endcaps->Write();
1127 h_ele_HoE_eg_endcaps->Write();
1128 h_ele_HoE_fiducial->Write();
1129 h_ele_HoEVsEta->Write();
1130 h_ele_HoEVsPhi->Write();
1131 h_ele_HoEVsE->Write();
1133 h_ele_seed_dphi2_->Write();
1134 h_ele_seed_subdet2_->Write();
1135 TProfile *p_ele_seed_dphi2VsEta_ = h_ele_seed_dphi2VsEta_->ProfileX();
1136 p_ele_seed_dphi2VsEta_->SetTitle(
"mean ele seed dphi 2nd layer vs eta");
1137 p_ele_seed_dphi2VsEta_->GetXaxis()->SetTitle(
"#eta");
1138 p_ele_seed_dphi2VsEta_->GetYaxis()->SetTitle(
"<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)");
1139 p_ele_seed_dphi2VsEta_->SetMinimum(-0.004);
1140 p_ele_seed_dphi2VsEta_->SetMaximum(0.004);
1141 p_ele_seed_dphi2VsEta_->Write();
1142 TProfile *p_ele_seed_dphi2VsPt_ = h_ele_seed_dphi2VsPt_->ProfileX();
1143 p_ele_seed_dphi2VsPt_->SetTitle(
"mean ele seed dphi 2nd layer vs pt");
1144 p_ele_seed_dphi2VsPt_->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
1145 p_ele_seed_dphi2VsPt_->GetYaxis()->SetTitle(
"<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)");
1146 p_ele_seed_dphi2VsPt_->SetMinimum(-0.004);
1147 p_ele_seed_dphi2VsPt_->SetMaximum(0.004);
1148 p_ele_seed_dphi2VsPt_->Write();
1149 h_ele_seed_drz2_->Write();
1150 TProfile *p_ele_seed_drz2VsEta_ = h_ele_seed_drz2VsEta_->ProfileX();
1151 p_ele_seed_drz2VsEta_->SetTitle(
"mean ele seed dr(dz) 2nd layer vs eta");
1152 p_ele_seed_drz2VsEta_->GetXaxis()->SetTitle(
"#eta");
1153 p_ele_seed_drz2VsEta_->GetYaxis()->SetTitle(
"<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)");
1154 p_ele_seed_drz2VsEta_->SetMinimum(-0.15);
1155 p_ele_seed_drz2VsEta_->SetMaximum(0.15);
1156 p_ele_seed_drz2VsEta_->Write();
1157 TProfile *p_ele_seed_drz2VsPt_ = h_ele_seed_drz2VsPt_->ProfileX();
1158 p_ele_seed_drz2VsPt_->SetTitle(
"mean ele seed dr(dz) 2nd layer vs pt");
1159 p_ele_seed_drz2VsPt_->GetXaxis()->SetTitle(
"p_{T} (GeV/c)");
1160 p_ele_seed_drz2VsPt_->GetYaxis()->SetTitle(
"<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)");
1161 p_ele_seed_drz2VsPt_->SetMinimum(-0.15);
1162 p_ele_seed_drz2VsPt_->SetMaximum(0.15);
1163 p_ele_seed_drz2VsPt_->Write();
1166 h_ele_classes->Write();
1168 h_ele_eta_golden->Write();
1169 h_ele_eta_bbrem->Write();
1170 h_ele_eta_narrow->Write();
1171 h_ele_eta_shower->Write();
1172 h_ele_PinVsPoutGolden_mode->Write();
1173 h_ele_PinVsPoutShowering_mode->Write();
1174 h_ele_PinVsPoutGolden_mean->Write();
1175 h_ele_PinVsPoutShowering_mean->Write();
1176 h_ele_PtinVsPtoutGolden_mode->Write();
1177 h_ele_PtinVsPtoutShowering_mode->Write();
1178 h_ele_PtinVsPtoutGolden_mean->Write();
1179 h_ele_PtinVsPtoutShowering_mean->Write();
1180 histSclEoEmatchingObjectGolden_barrel->Write();
1181 histSclEoEmatchingObjectGolden_endcaps->Write();
1182 histSclEoEmatchingObjectShowering_barrel->Write();
1183 histSclEoEmatchingObjectShowering_endcaps->Write();
1186 h_ele_fbrem->Write();
1187 h_ele_fbremVsEta_mode->Write();
1188 h_ele_fbremVsEta_mean->Write();
1189 h_ele_etaEff->Write();
1190 h_ele_zEff->Write();
1191 h_ele_phiEff->Write();
1192 h_ele_absetaEff->Write();
1193 h_ele_ptEff->Write();
1194 h_ele_etaEff_all->Write();
1195 h_ele_ptEff_all->Write();
1196 h_ele_eta_goldenFrac->Write();
1197 h_ele_eta_bbremFrac->Write();
1198 h_ele_eta_narrowFrac->Write();
1199 h_ele_eta_showerFrac->Write();
1200 h_ele_xOverX0VsEta->Write();
1204 h_ele_provenance->Write();
1207 h_ele_tkSumPt_dr03->GetXaxis()->SetTitle(
"TkIsoSum, cone 0.3 (GeV/c)");
1208 h_ele_tkSumPt_dr03->GetYaxis()->SetTitle(
"Events");
1209 h_ele_tkSumPt_dr03->Write();
1210 h_ele_ecalRecHitSumEt_dr03->GetXaxis()->SetTitle(
"EcalIsoSum, cone 0.3 (GeV)");
1211 h_ele_ecalRecHitSumEt_dr03->GetYaxis()->SetTitle(
"Events");
1212 h_ele_ecalRecHitSumEt_dr03->Write();
1213 h_ele_hcalDepth1TowerSumEt_dr03->GetXaxis()->SetTitle(
"Hcal1IsoSum, cone 0.3 (GeV)");
1214 h_ele_hcalDepth1TowerSumEt_dr03->GetYaxis()->SetTitle(
"Events");
1215 h_ele_hcalDepth1TowerSumEt_dr03->Write();
1216 h_ele_hcalDepth2TowerSumEt_dr03->GetXaxis()->SetTitle(
"Hcal2IsoSum, cone 0.3 (GeV)");
1217 h_ele_hcalDepth2TowerSumEt_dr03->GetYaxis()->SetTitle(
"Events");
1218 h_ele_hcalDepth2TowerSumEt_dr03->Write();
1219 h_ele_tkSumPt_dr04->GetXaxis()->SetTitle(
"TkIsoSum, cone 0.4 (GeV/c)");
1220 h_ele_tkSumPt_dr04->GetYaxis()->SetTitle(
"Events");
1221 h_ele_tkSumPt_dr04->Write();
1222 h_ele_ecalRecHitSumEt_dr04->GetXaxis()->SetTitle(
"EcalIsoSum, cone 0.4 (GeV)");
1223 h_ele_ecalRecHitSumEt_dr04->GetYaxis()->SetTitle(
"Events");
1224 h_ele_ecalRecHitSumEt_dr04->Write();
1225 h_ele_hcalDepth1TowerSumEt_dr04->GetXaxis()->SetTitle(
"Hcal1IsoSum, cone 0.4 (GeV)");
1226 h_ele_hcalDepth1TowerSumEt_dr04->GetYaxis()->SetTitle(
"Events");
1227 h_ele_hcalDepth1TowerSumEt_dr04->Write();
1228 h_ele_hcalDepth2TowerSumEt_dr04->GetXaxis()->SetTitle(
"Hcal2IsoSum, cone 0.4 (GeV)");
1229 h_ele_hcalDepth2TowerSumEt_dr04->GetYaxis()->SetTitle(
"Events");
1230 h_ele_hcalDepth2TowerSumEt_dr04->Write();
1237 std::cout <<
"analyzing new event " << std::endl;
1241 iEvent.
getByLabel(electronCollection_,gsfElectrons);
1242 edm::LogInfo(
"")<<
"\n\n =================> Treating event "<<iEvent.
id()<<
" Number of electrons "<<gsfElectrons.
product()->size();
1246 iEvent.
getByLabel(matchingObjectCollection_,genJets);
1251 const BeamSpot bs = *recoBeamSpotHandle;
1253 histNum_->Fill((*gsfElectrons).size());
1256 for (reco::GsfElectronCollection::const_iterator gsfIter=gsfElectrons->begin();
1257 gsfIter!=gsfElectrons->end(); gsfIter++){
1259 if (gsfIter->pt()>maxPt_ ||
std::abs(gsfIter->eta())>maxAbsEta_)
continue;
1260 h_ele_EoverP_all -> Fill( gsfIter->eSuperClusterOverP() );
1261 h_ele_EseedOP_all -> Fill( gsfIter->eSeedClusterOverP() );
1262 h_ele_EoPout_all -> Fill( gsfIter->eSeedClusterOverPout() );
1263 h_ele_EeleOPout_all -> Fill( gsfIter->eEleClusterOverPout() );
1264 h_ele_dEtaSc_propVtx_all -> Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
1265 h_ele_dPhiSc_propVtx_all -> Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
1266 h_ele_dEtaCl_propOut_all -> Fill(gsfIter->deltaEtaSeedClusterTrackAtCalo());
1267 h_ele_dPhiCl_propOut_all -> Fill(gsfIter->deltaPhiSeedClusterTrackAtCalo());
1268 h_ele_HoE_all -> Fill( gsfIter->hadronicOverEm() );
1269 double d = gsfIter->vertex().x()*gsfIter->vertex().x()
1270 +gsfIter->vertex().y()*gsfIter->vertex().y();
1271 h_ele_TIP_all -> Fill(
sqrt(d) );
1272 h_ele_vertexEta_all -> Fill( gsfIter->eta() );
1273 h_ele_vertexPt_all -> Fill( gsfIter->pt() );
1274 float enrj1=gsfIter->superCluster()->energy();
1276 for (reco::GsfElectronCollection::const_iterator gsfIter2=gsfIter+1;
1277 gsfIter2!=gsfElectrons->end(); gsfIter2++){
1279 float mee2 = p12.Dot(p12);
1280 h_ele_mee_all -> Fill(
sqrt(mee2));
1281 float enrj2=gsfIter2->superCluster()->energy();
1282 h_ele_E2mnE1vsMee_all->Fill(
sqrt(mee2),enrj2-enrj1);
1283 if (gsfIter->ecalDrivenSeed() && gsfIter2->ecalDrivenSeed()) h_ele_E2mnE1vsMee_egeg_all->Fill(
sqrt(mee2),enrj2-enrj1);
1284 if (gsfIter->charge()*gsfIter2->charge()<0.) h_ele_mee_os -> Fill(
sqrt(mee2));
1289 int matchingObjectNum=0;
1291 for (reco::GenJetCollection::const_iterator moIter = genJets->begin();
1292 moIter != genJets->end(); ++moIter) {
1295 matchingObjectNum++;
1297 if (moIter->energy()/cosh(moIter->eta())> maxPt_ ||
std::abs(moIter->eta())> maxAbsEta_)
continue;
1304 h_matchingObjectEta -> Fill( moIter->eta() );
1305 h_matchingObjectAbsEta -> Fill(
std::abs(moIter->eta()) );
1306 h_matchingObjectP -> Fill( moIter->energy() );
1307 h_matchingObjectPt -> Fill( moIter->energy()/cosh(moIter->eta()) );
1308 h_matchingObjectPhi -> Fill( moIter->phi() );
1309 h_matchingObjectZ -> Fill( moIter->vz() );
1312 bool okGsfFound =
false;
1313 double gsfOkRatio = 999999.;
1317 for (reco::GsfElectronCollection::const_iterator gsfIter=gsfElectrons->begin();
1318 gsfIter!=gsfElectrons->end(); gsfIter++){
1320 double dphi = gsfIter->
phi()-moIter->phi();
1322 dphi = dphi < 0? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
1324 if ( deltaR < deltaR_ ){
1327 double tmpGsfRatio = gsfIter->p()/moIter->energy();
1329 gsfOkRatio = tmpGsfRatio;
1330 bestGsfElectron=*gsfIter;
1341 h_ele_charge -> Fill( bestGsfElectron.
charge() );
1342 h_ele_chargeVsEta -> Fill( bestGsfElectron.
eta(),bestGsfElectron.
charge() );
1343 h_ele_chargeVsPhi -> Fill( bestGsfElectron.
phi(),bestGsfElectron.
charge() );
1344 h_ele_chargeVsPt -> Fill( bestGsfElectron.
pt(),bestGsfElectron.
charge() );
1345 h_ele_vertexP -> Fill( bestGsfElectron.
p() );
1346 h_ele_vertexPt -> Fill( bestGsfElectron.
pt() );
1347 h_ele_vertexPtVsEta -> Fill( bestGsfElectron.
eta(),bestGsfElectron.
pt() );
1348 h_ele_vertexPtVsPhi -> Fill( bestGsfElectron.
phi(),bestGsfElectron.
pt() );
1349 h_ele_vertexEta -> Fill( bestGsfElectron.
eta() );
1351 h_ele_matchingObjectPt_matched -> Fill( moIter->energy()/cosh(moIter->eta()) );
1352 h_ele_matchingObjectPhi_matched -> Fill( moIter->phi() );
1353 h_ele_matchingObjectAbsEta_matched -> Fill(
std::abs(moIter->eta()) );
1354 h_ele_matchingObjectEta_matched -> Fill( moIter->eta() );
1355 h_ele_vertexEtaVsPhi -> Fill( bestGsfElectron.
phi(),bestGsfElectron.
eta() );
1356 h_ele_vertexPhi -> Fill( bestGsfElectron.
phi() );
1357 h_ele_vertexX -> Fill( bestGsfElectron.
vertex().x() );
1358 h_ele_vertexY -> Fill( bestGsfElectron.
vertex().y() );
1359 h_ele_vertexZ -> Fill( bestGsfElectron.
vertex().z() );
1360 h_ele_matchingObjectZ_matched -> Fill( moIter->vz() );
1366 h_ele_vertexTIP -> Fill( d );
1367 h_ele_vertexTIPVsEta -> Fill( bestGsfElectron.
eta(), d );
1368 h_ele_vertexTIPVsPhi -> Fill( bestGsfElectron.
phi(), d );
1369 h_ele_vertexTIPVsPt -> Fill( bestGsfElectron.
pt(), d );
1370 h_ele_EtaMnEtamatchingObject -> Fill( bestGsfElectron.
eta()-moIter->eta());
1371 h_ele_EtaMnEtamatchingObjectVsEta -> Fill( bestGsfElectron.
eta(), bestGsfElectron.
eta()-moIter->eta());
1372 h_ele_EtaMnEtamatchingObjectVsPhi -> Fill( bestGsfElectron.
phi(), bestGsfElectron.
eta()-moIter->eta());
1373 h_ele_EtaMnEtamatchingObjectVsPt -> Fill( bestGsfElectron.
pt(), bestGsfElectron.
eta()-moIter->eta());
1374 h_ele_PhiMnPhimatchingObject -> Fill( bestGsfElectron.
phi()-moIter->phi());
1375 h_ele_PhiMnPhimatchingObject2 -> Fill( bestGsfElectron.
phi()-moIter->phi());
1376 h_ele_PhiMnPhimatchingObjectVsEta -> Fill( bestGsfElectron.
eta(), bestGsfElectron.
phi()-moIter->phi());
1377 h_ele_PhiMnPhimatchingObjectVsPhi -> Fill( bestGsfElectron.
phi(), bestGsfElectron.
phi()-moIter->phi());
1378 h_ele_PhiMnPhimatchingObjectVsPt -> Fill( bestGsfElectron.
pt(), bestGsfElectron.
phi()-moIter->phi());
1379 h_ele_PoPmatchingObject -> Fill( bestGsfElectron.
p()/moIter->energy());
1380 h_ele_PoPmatchingObjectVsEta -> Fill( bestGsfElectron.
eta(), bestGsfElectron.
p()/moIter->energy());
1381 h_ele_PoPmatchingObjectVsPhi -> Fill( bestGsfElectron.
phi(), bestGsfElectron.
p()/moIter->energy());
1382 h_ele_PoPmatchingObjectVsPt -> Fill( bestGsfElectron.
py(), bestGsfElectron.
p()/moIter->energy());
1383 if (bestGsfElectron.
isEB()) h_ele_PoPmatchingObject_barrel -> Fill( bestGsfElectron.
p()/moIter->energy());
1384 if (bestGsfElectron.
isEE()) h_ele_PoPmatchingObject_endcaps -> Fill( bestGsfElectron.
p()/moIter->energy());
1389 histSclEn_->Fill(sclRef->energy());
1390 double R=TMath::Sqrt(sclRef->x()*sclRef->x() + sclRef->y()*sclRef->y() +sclRef->z()*sclRef->z());
1391 double Rt=TMath::Sqrt(sclRef->x()*sclRef->x() + sclRef->y()*sclRef->y());
1392 histSclEt_->Fill(sclRef->energy()*(Rt/
R));
1393 histSclEtVsEta_->Fill(sclRef->eta(),sclRef->energy()*(Rt/
R));
1394 histSclEtVsPhi_->Fill(sclRef->phi(),sclRef->energy()*(Rt/
R));
1395 if (bestGsfElectron.
classification() < 100) histSclEoEmatchingObject_barrel->Fill(sclRef->energy()/moIter->energy());
1396 if (bestGsfElectron.
classification() >= 100) histSclEoEmatchingObject_endcaps->Fill(sclRef->energy()/moIter->energy());
1397 histSclEta_->Fill(sclRef->eta());
1398 histSclEtaVsPhi_->Fill(sclRef->phi(),sclRef->eta());
1399 histSclPhi_->Fill(sclRef->phi());
1401 if (bestGsfElectron.
isEB()) histSclSigIEtaIEta_barrel_->Fill(bestGsfElectron.
scSigmaIEtaIEta());
1402 if (bestGsfElectron.
isEE()) histSclSigIEtaIEta_endcaps_->Fill(bestGsfElectron.
scSigmaIEtaIEta());
1403 histSclE1x5_->Fill(bestGsfElectron.
scE1x5());
1404 if (bestGsfElectron.
isEB()) histSclE1x5_barrel_->Fill(bestGsfElectron.
scE1x5());
1405 if (bestGsfElectron.
isEE()) histSclE1x5_endcaps_->Fill(bestGsfElectron.
scE1x5());
1406 histSclE2x5max_->Fill(bestGsfElectron.
scE2x5Max());
1407 if (bestGsfElectron.
isEB()) histSclE2x5max_barrel_->Fill(bestGsfElectron.
scE2x5Max());
1408 if (bestGsfElectron.
isEE()) histSclE2x5max_endcaps_->Fill(bestGsfElectron.
scE2x5Max());
1409 histSclE5x5_->Fill(bestGsfElectron.
scE5x5());
1410 if (bestGsfElectron.
isEB()) histSclE5x5_barrel_->Fill(bestGsfElectron.
scE5x5());
1411 if (bestGsfElectron.
isEE()) histSclE5x5_endcaps_->Fill(bestGsfElectron.
scE5x5());
1419 h_ele_foundHits -> Fill( bestGsfElectron.
gsfTrack()->numberOfValidHits() );
1420 h_ele_foundHitsVsEta -> Fill( bestGsfElectron.
eta(), bestGsfElectron.
gsfTrack()->numberOfValidHits() );
1421 h_ele_foundHitsVsPhi -> Fill( bestGsfElectron.
phi(), bestGsfElectron.
gsfTrack()->numberOfValidHits() );
1422 h_ele_foundHitsVsPt -> Fill( bestGsfElectron.
pt(), bestGsfElectron.
gsfTrack()->numberOfValidHits() );
1423 h_ele_lostHits -> Fill( bestGsfElectron.
gsfTrack()->numberOfLostHits() );
1424 h_ele_lostHitsVsEta -> Fill( bestGsfElectron.
eta(), bestGsfElectron.
gsfTrack()->numberOfLostHits() );
1425 h_ele_lostHitsVsPhi -> Fill( bestGsfElectron.
phi(), bestGsfElectron.
gsfTrack()->numberOfLostHits() );
1426 h_ele_lostHitsVsPt -> Fill( bestGsfElectron.
pt(), bestGsfElectron.
gsfTrack()->numberOfLostHits() );
1427 h_ele_chi2 -> Fill( bestGsfElectron.
gsfTrack()->normalizedChi2() );
1428 h_ele_chi2VsEta -> Fill( bestGsfElectron.
eta(), bestGsfElectron.
gsfTrack()->normalizedChi2() );
1429 h_ele_chi2VsPhi -> Fill( bestGsfElectron.
phi(), bestGsfElectron.
gsfTrack()->normalizedChi2() );
1430 h_ele_chi2VsPt -> Fill( bestGsfElectron.
pt(), bestGsfElectron.
gsfTrack()->normalizedChi2() );
1434 h_ele_PinMnPout -> Fill( bestGsfElectron.
gsfTrack()->innerMomentum().R() - bestGsfElectron.
gsfTrack()->outerMomentum().R() );
1435 h_ele_outerP -> Fill( bestGsfElectron.
gsfTrack()->outerMomentum().R() );
1436 h_ele_outerPt -> Fill( bestGsfElectron.
gsfTrack()->outerMomentum().Rho() );
1447 h_ele_outerPVsEta_mode -> Fill(bestGsfElectron.
eta(), bestGsfElectron.
trackMomentumOut().R() );
1449 h_ele_outerPtVsEta_mode -> Fill(bestGsfElectron.
eta(), bestGsfElectron.
trackMomentumOut().Rho() );
1450 h_ele_outerPtVsPhi_mode -> Fill(bestGsfElectron.
phi(), bestGsfElectron.
trackMomentumOut().Rho() );
1451 h_ele_outerPtVsPt_mode -> Fill(bestGsfElectron.
pt(), bestGsfElectron.
trackMomentumOut().Rho() );
1456 h_ele_seed_dphi2_-> Fill(elseed->dPhi2());
1457 h_ele_seed_dphi2VsEta_-> Fill(bestGsfElectron.
eta(), elseed->dPhi2());
1458 h_ele_seed_dphi2VsPt_-> Fill(bestGsfElectron.
pt(), elseed->dPhi2()) ;
1459 h_ele_seed_drz2_-> Fill(elseed->dRz2());
1460 h_ele_seed_drz2VsEta_-> Fill(bestGsfElectron.
eta(), elseed->dRz2());
1461 h_ele_seed_drz2VsPt_-> Fill(bestGsfElectron.
pt(), elseed->dRz2());
1462 h_ele_seed_subdet2_-> Fill(elseed->subDet2());
1557 if (bestGsfElectron.
isEB()) h_ele_HoE_barrel -> Fill(bestGsfElectron.
hadronicOverEm());
1559 if (bestGsfElectron.
isEE()) h_ele_HoE_endcaps -> Fill(bestGsfElectron.
hadronicOverEm());
1563 h_ele_HoEVsEta -> Fill( bestGsfElectron.
eta(),bestGsfElectron.
hadronicOverEm());
1569 if (bestGsfElectron.
isEE()) eleClass+=10;
1570 h_ele_classes ->Fill(eleClass);
1579 double fbrem_mean=0.;
1581 fbrem_mean = 1. - bestGsfElectron.
gsfTrack()->outerMomentum().R()/bestGsfElectron.
gsfTrack()->innerMomentum().R();
1582 double fbrem_mode = bestGsfElectron.
fbrem();
1583 h_ele_fbrem->Fill(fbrem_mode);
1584 h_ele_fbremVsEta_mode->Fill(bestGsfElectron.
eta(),fbrem_mode);
1586 h_ele_fbremVsEta_mean->Fill(bestGsfElectron.
eta(),fbrem_mean);
1595 h_ele_PinVsPoutShowering_mean -> Fill(bestGsfElectron.
gsfTrack()->outerMomentum().R(), bestGsfElectron.
gsfTrack()->innerMomentum().R());
1603 h_ele_PtinVsPtoutShowering_mean -> Fill(bestGsfElectron.
gsfTrack()->outerMomentum().Rho(), bestGsfElectron.
gsfTrack()->innerMomentum().Rho());
1605 h_ele_mva->Fill(bestGsfElectron.
mva());
1606 if (bestGsfElectron.
ecalDrivenSeed()) h_ele_provenance->Fill(1.);
1612 h_ele_tkSumPt_dr03->Fill(bestGsfElectron.
dr03TkSumPt());
1616 h_ele_tkSumPt_dr04->Fill(bestGsfElectron.
dr04TkSumPt());
1629 h_matchingObjectNum->Fill(matchingObjectNum);
T getParameter(std::string const &) const
virtual double p() const
magnitude of momentum vector
float eSuperClusterOverP() const
float dr04HcalDepth1TowerSumEt() const
virtual const Point & vertex() const
vertex position
float dr04TkSumPt() const
bool getByType(Handle< PROD > &result) const
float dr03HcalDepth2TowerSumEt() const
SuperClusterRef superCluster() const
reference to a SuperCluster
virtual double eta() const
momentum pseudorapidity
double deltaR(double eta1, double phi1, double eta2, double phi2)
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
float dr04HcalDepth2TowerSumEt() const
float deltaEtaSuperClusterTrackAtVtx() const
float hadronicOverEm() const
float eSeedClusterOverP() const
float dr04EcalRecHitSumEt() const
math::XYZVector trackMomentumOut() const
math::XYZVector trackMomentumAtVtx() const
float deltaPhiSeedClusterTrackAtCalo() const
float deltaPhiSuperClusterTrackAtVtx() const
virtual int charge() const
electric charge
float dr03TkSumPt() const
float dr03HcalDepth1TowerSumEt() const
float eEleClusterOverPout() const
float eSeedClusterOverPout() const
float scSigmaIEtaIEta() const
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
float deltaPhiEleClusterTrackAtCalo() const
GsfTrackRef gsfTrack() const
reference to a GsfTrack
Log< T >::type log(const T &t)
float deltaEtaEleClusterTrackAtCalo() const
virtual double pt() const
transverse momentum
Classification classification() const
REF castTo() const
cast to a concrete type
GsfTrackRefVector::size_type ambiguousGsfTracksSize() const
T const * product() const
float dr03EcalRecHitSumEt() const
float deltaEtaSeedClusterTrackAtCalo() const
bool trackerDrivenSeed() const
const Point & position() const
position
virtual ~GsfElectronMCFakeAnalyzer()
virtual double phi() const
momentum azimuthal angle
SuperClusterRef pflowSuperCluster() const
virtual double py() const
y coordinate of momentum vector
Power< A, B >::type pow(const A &a, const B &b)
GsfElectronMCFakeAnalyzer(const edm::ParameterSet &conf)
bool ecalDrivenSeed() const