CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
GsfElectronMCAnalyzer.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: RecoEgamma/Examples
4 // Class: GsfElectronMCAnalyzer
5 //
13 //
14 // Original Author: Ursula Berthon
15 // Created: Mon Mar 27 13:22:06 CEST 2006
16 //
17 //
18 
19 // user include files
20 #include "RecoEgamma/Examples/plugins/GsfElectronMCAnalyzer.h"
21 #include "FWCore/ParameterSet/interface/ParameterSet.h"
22 #include "FWCore/Framework/interface/EDAnalyzer.h"
23 #include "FWCore/Framework/interface/Event.h"
24 #include "FWCore/Framework/interface/MakerMacros.h"
25 #include "FWCore/MessageLogger/interface/MessageLogger.h"
26 #include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
27 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
28 #include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
29 #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
30 #include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
31 #include "DataFormats/EgammaReco/interface/ElectronSeed.h"
32 #include "DataFormats/EgammaReco/interface/ElectronSeedFwd.h"
33 #include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
34 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
35 #include "DataFormats/HepMCCandidate/interface/GenParticleFwd.h"
36 
37 #include "CLHEP/Units/GlobalPhysicalConstants.h"
38 #include <iostream>
39 #include <vector>
40 #include "TMath.h"
41 #include "TFile.h"
42 #include "TH1F.h"
43 #include "TH1I.h"
44 #include "TH2F.h"
45 #include "TProfile.h"
46 #include "TTree.h"
47 #include <iostream>
48 
49 using namespace reco;
50 
51 GsfElectronMCAnalyzer::GsfElectronMCAnalyzer(const edm::ParameterSet& conf)
52  {
53  outputFile_ = conf.getParameter<std::string>("outputFile");
54  histfile_ = new TFile(outputFile_.c_str(),"RECREATE");
55  electronCollection_=conf.getParameter<edm::InputTag>("electronCollection");
56  mcTruthCollection_ = conf.getParameter<edm::InputTag>("mcTruthCollection");
57  readAOD_ = conf.getParameter<bool>("readAOD");
58  maxPt_ = conf.getParameter<double>("MaxPt");
59  maxAbsEta_ = conf.getParameter<double>("MaxAbsEta");
60  deltaR_ = conf.getParameter<double>("DeltaR");
61  matchingIDs_ = conf.getParameter<std::vector<int> >("MatchingID");
62  matchingMotherIDs_ = conf.getParameter<std::vector<int> >("MatchingMotherID");
63  edm::ParameterSet pset =
64  conf.getParameter<edm::ParameterSet>("HistosConfigurationMC") ;
65 
66  etamin=pset.getParameter<double>("Etamin");
67  etamax=pset.getParameter<double>("Etamax");
68  phimin=pset.getParameter<double>("Phimin");
69  phimax=pset.getParameter<double>("Phimax");
70  ptmax=pset.getParameter<double>("Ptmax");
71  pmax=pset.getParameter<double>("Pmax");
72  eopmax=pset.getParameter<double>("Eopmax");
73  eopmaxsht=pset.getParameter<double>("Eopmaxsht");
74  detamin=pset.getParameter<double>("Detamin");
75  detamax=pset.getParameter<double>("Detamax");
76  dphimin=pset.getParameter<double>("Dphimin");
77  dphimax=pset.getParameter<double>("Dphimax");
78  detamatchmin=pset.getParameter<double>("Detamatchmin");
79  detamatchmax=pset.getParameter<double>("Detamatchmax");
80  dphimatchmin=pset.getParameter<double>("Dphimatchmin");
81  dphimatchmax=pset.getParameter<double>("Dphimatchmax");
82  fhitsmax=pset.getParameter<double>("Fhitsmax");
83  lhitsmax=pset.getParameter<double>("Lhitsmax");
84  nbineta=pset.getParameter<int>("Nbineta");
85  nbineta2D=pset.getParameter<int>("Nbineta2D");
86  nbinp=pset.getParameter<int>("Nbinp");
87  nbinpt=pset.getParameter<int>("Nbinpt");
88  nbinp2D=pset.getParameter<int>("Nbinp2D");
89  nbinpt2D=pset.getParameter<int>("Nbinpt2D");
90  nbinpteff=pset.getParameter<int>("Nbinpteff");
91  nbinphi=pset.getParameter<int>("Nbinphi");
92  nbinphi2D=pset.getParameter<int>("Nbinphi2D");
93  nbineop=pset.getParameter<int>("Nbineop");
94  nbineop2D=pset.getParameter<int>("Nbineop2D");
95  nbinfhits=pset.getParameter<int>("Nbinfhits");
96  nbinlhits=pset.getParameter<int>("Nbinlhits");
97  nbinxyz=pset.getParameter<int>("Nbinxyz");
98  nbindeta=pset.getParameter<int>("Nbindeta");
99  nbindphi=pset.getParameter<int>("Nbindphi");
100  nbindetamatch=pset.getParameter<int>("Nbindetamatch");
101  nbindphimatch=pset.getParameter<int>("Nbindphimatch");
102  nbindetamatch2D=pset.getParameter<int>("Nbindetamatch2D");
103  nbindphimatch2D=pset.getParameter<int>("Nbindphimatch2D");
104  nbinpoptrue= pset.getParameter<int>("Nbinpoptrue");
105  poptruemin=pset.getParameter<double>("Poptruemin");
106  poptruemax=pset.getParameter<double>("Poptruemax");
107  nbinmee= pset.getParameter<int>("Nbinmee");
108  meemin=pset.getParameter<double>("Meemin");
109  meemax=pset.getParameter<double>("Meemax");
110  nbinhoe= pset.getParameter<int>("Nbinhoe");
111  hoemin=pset.getParameter<double>("Hoemin");
112  hoemax=pset.getParameter<double>("Hoemax");
113  }
114 
115 void GsfElectronMCAnalyzer::beginJob(){
116 
117  histfile_->cd();
118 
119  // mc truth
120  h_mcNum = new TH1F( "h_mcNum", "# mc particles", nbinfhits,0.,fhitsmax );
121  h_mcNum->Sumw2();
122  h_eleNum = new TH1F( "h_mcNum_ele", "# mc electrons", nbinfhits,0.,fhitsmax);
123  h_eleNum->Sumw2();
124  h_gamNum = new TH1F( "h_mcNum_gam", "# mc gammas", nbinfhits,0.,fhitsmax);
125  h_gamNum->Sumw2();
126 
127  // rec event
128  histNum_= new TH1F("h_recEleNum","# rec electrons",20, 0.,20.);
129 
130  // mc
131  h_simEta = new TH1F( "h_mc_eta", "gen #eta", nbineta,etamin,etamax);
132  h_simEta->Sumw2();
133  h_simAbsEta = new TH1F( "h_mc_abseta", "gen |#eta|", nbineta/2,0.,etamax);
134  h_simAbsEta->Sumw2();
135  h_simP = new TH1F( "h_mc_P", "gen p", nbinp,0.,pmax);
136  h_simP->Sumw2();
137  h_simPt = new TH1F( "h_mc_Pt", "gen pt", nbinpteff,5.,ptmax);
138  h_simPt->Sumw2();
139  h_simPhi = new TH1F( "h_mc_phi", "gen phi", nbinphi,phimin,phimax);
140  h_simPhi->Sumw2();
141  h_simZ = new TH1F( "h_mc_z", "gen z ", nbinxyz, -25, 25 );
142  h_simZ->Sumw2();
143  h_simPtEta = new TH2F( "h_mc_pteta", "gen pt vs #eta", nbineta2D,etamin,etamax, nbinpt2D,5.,ptmax );
144  h_simPtEta->Sumw2();
145 
146  // all electrons
147  h_ele_EoverP_all = new TH1F( "h_ele_EoverP_all", "ele E/P_{vertex}, all reco electrons", nbineop,0.,eopmax);
148  h_ele_EoverP_all->Sumw2();
149  h_ele_EoverP_all_barrel = new TH1F( "h_ele_EoverP_all_barrel", "ele E/P_{vertex}, all reco electrons, barrel", nbineop,0.,eopmax);
150  h_ele_EoverP_all_barrel->Sumw2();
151  h_ele_EoverP_all_endcaps = new TH1F( "h_ele_EoverP_all_endcaps", "ele E/P_{vertex}, all reco electrons, endcaps", nbineop,0.,eopmax);
152  h_ele_EoverP_all_endcaps->Sumw2();
153  h_ele_EseedOP_all = new TH1F( "h_ele_EseedOP_all", "ele E_{seed}/P_{vertex}, all reco electrons", nbineop,0.,eopmax);
154  h_ele_EseedOP_all->Sumw2();
155  h_ele_EseedOP_all_barrel = new TH1F( "h_ele_EseedOP_all_barrel", "ele E_{seed}/P_{vertex}, all reco electrons, barrel", nbineop,0.,eopmax);
156  h_ele_EseedOP_all_barrel->Sumw2();
157  h_ele_EseedOP_all_endcaps = new TH1F( "h_ele_EseedOP_all_endcaps", "ele E_{seed}/P_{vertex}, all reco electrons, endcaps", nbineop,0.,eopmax);
158  h_ele_EseedOP_all_endcaps->Sumw2();
159  h_ele_EoPout_all = new TH1F( "h_ele_EoPout_all", "ele E_{seed}/P_{out}, all reco electrons", nbineop,0.,eopmax);
160  h_ele_EoPout_all->Sumw2();
161  h_ele_EoPout_all_barrel = new TH1F( "h_ele_EoPout_all_barrel", "ele E_{seed}/P_{out}, all reco electrons barrel", nbineop,0.,eopmax);
162  h_ele_EoPout_all_barrel->Sumw2();
163  h_ele_EoPout_all_endcaps = new TH1F( "h_ele_EoPout_all_endcaps", "ele E_{seed}/P_{out}, all reco electrons endcaps", nbineop,0.,eopmax);
164  h_ele_EoPout_all_endcaps->Sumw2();
165  h_ele_EeleOPout_all = new TH1F( "h_ele_EeleOPout_all", "ele E_{ele}/P_{out}, all reco electrons", nbineop,0.,eopmax);
166  h_ele_EeleOPout_all->Sumw2();
167  h_ele_EeleOPout_all_barrel = new TH1F( "h_ele_EeleOPout_all_barrel", "ele E_{ele}/P_{out}, all reco electrons barrel", nbineop,0.,eopmax);
168  h_ele_EeleOPout_all_barrel->Sumw2();
169  h_ele_EeleOPout_all_endcaps = new TH1F( "h_ele_EeleOPout_all_endcaps", "ele E_{ele}/P_{out}, all reco electrons endcaps", nbineop,0.,eopmax);
170  h_ele_EeleOPout_all_endcaps->Sumw2();
171  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);
172  h_ele_dEtaSc_propVtx_all->Sumw2();
173  h_ele_dEtaSc_propVtx_all_barrel = new TH1F( "h_ele_dEtaSc_propVtx_all_barrel", "ele #eta_{sc} - #eta_{tr}, prop from vertex, all reco electrons barrel", nbindetamatch,detamatchmin,detamatchmax);
174  h_ele_dEtaSc_propVtx_all_barrel->Sumw2();
175  h_ele_dEtaSc_propVtx_all_endcaps = new TH1F( "h_ele_dEtaSc_propVtx_all_endcaps", "ele #eta_{sc} - #eta_{tr}, prop from vertex, all reco electrons endcaps", nbindetamatch,detamatchmin,detamatchmax);
176  h_ele_dEtaSc_propVtx_all_endcaps->Sumw2();
177  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);
178  h_ele_dPhiSc_propVtx_all->Sumw2();
179  h_ele_dPhiSc_propVtx_all_barrel = new TH1F( "h_ele_dPhiSc_propVtx_all_barrel", "ele #phi_{sc} - #phi_{tr}, prop from vertex, all reco electrons barrel", nbindphimatch,dphimatchmin,dphimatchmax);
180  h_ele_dPhiSc_propVtx_all_barrel->Sumw2();
181  h_ele_dPhiSc_propVtx_all_endcaps = new TH1F( "h_ele_dPhiSc_propVtx_all_endcaps", "ele #phi_{sc} - #phi_{tr}, prop from vertex, all reco electrons endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
182  h_ele_dPhiSc_propVtx_all_endcaps->Sumw2();
183  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);
184  h_ele_dEtaCl_propOut_all->Sumw2();
185  h_ele_dEtaCl_propOut_all_barrel = new TH1F( "h_ele_dEtaCl_propOut_all_barrel", "ele #eta_{cl} - #eta_{tr}, prop from outermost, all reco electrons barrel", nbindetamatch,detamatchmin,detamatchmax);
186  h_ele_dEtaCl_propOut_all_barrel->Sumw2();
187  h_ele_dEtaCl_propOut_all_endcaps = new TH1F( "h_ele_dEtaCl_propOut_all_endcaps", "ele #eta_{cl} - #eta_{tr}, prop from outermost, all reco electrons endcaps", nbindetamatch,detamatchmin,detamatchmax);
188  h_ele_dEtaCl_propOut_all_endcaps->Sumw2();
189  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);
190  h_ele_dPhiCl_propOut_all->Sumw2();
191  h_ele_dPhiCl_propOut_all_barrel = new TH1F( "h_ele_dPhiCl_propOut_all_barrel", "ele #phi_{cl} - #phi_{tr}, prop from outermost, all reco electrons barrel", nbindphimatch,dphimatchmin,dphimatchmax);
192  h_ele_dPhiCl_propOut_all_endcaps = new TH1F( "h_ele_dPhiCl_propOut_all_endcaps", "ele #phi_{cl} - #phi_{tr}, prop from outermost, all reco electrons endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
193  h_ele_dPhiCl_propOut_all_barrel->Sumw2();
194  h_ele_dPhiCl_propOut_all_endcaps->Sumw2();
195  h_ele_HoE_all = new TH1F("h_ele_HoE_all", "ele hadronic energy / em energy, all reco electrons", nbinhoe, hoemin, hoemax) ;
196  h_ele_HoE_all->Sumw2();
197  h_ele_HoE_all_barrel = new TH1F("h_ele_HoE_all_barrel", "ele hadronic energy / em energy, all reco electrons barrel", nbinhoe, hoemin, hoemax) ;
198  h_ele_HoE_all_barrel->Sumw2();
199  h_ele_HoE_all_endcaps = new TH1F("h_ele_HoE_all_endcaps", "ele hadronic energy / em energy, all reco electrons endcaps", nbinhoe, hoemin, hoemax) ;
200  h_ele_HoE_all_endcaps->Sumw2();
201  h_ele_vertexPt_all = new TH1F( "h_ele_vertexPt_all", "ele p_{T}, all reco electrons", nbinpteff,5.,ptmax);
202  h_ele_vertexPt_all->Sumw2();
203  h_ele_Et_all = new TH1F( "h_ele_Et_all", "ele SC E_{T}, all reco electrons", nbinpteff,5.,ptmax);
204  h_ele_Et_all->Sumw2();
205  h_ele_vertexEta_all = new TH1F( "h_ele_vertexEta_all", "ele eta, all reco electrons", nbineta,etamin,etamax);
206  h_ele_vertexEta_all->Sumw2();
207  h_ele_TIP_all = new TH1F( "h_ele_TIP_all", "ele vertex transverse radius, all reco electrons", 100,0.,0.2);
208  h_ele_TIP_all->Sumw2();
209  h_ele_TIP_all_barrel = new TH1F( "h_ele_TIP_all_barrel", "ele vertex transverse radius, all reco electrons barrel", 100,0.,0.2);
210  h_ele_TIP_all_barrel->Sumw2();
211  h_ele_TIP_all_endcaps = new TH1F( "h_ele_TIP_all_endcaps", "ele vertex transverse radius, all reco electrons endcaps", 100,0.,0.2);
212  h_ele_TIP_all_endcaps->Sumw2();
213  h_ele_mee_all = new TH1F( "h_ele_mee_all", "ele pairs invariant mass, all reco electrons", nbinmee, meemin, meemax );
214  h_ele_mee_all->Sumw2();
215  h_ele_mee_os = new TH1F( "h_ele_mee_os", "ele pairs invariant mass, opp. sign", nbinmee, meemin, meemax );
216  h_ele_mee_os->Sumw2();
217  h_ele_mee_os_ebeb = new TH1F( "h_ele_mee_os_ebeb", "ele pairs invariant mass, opp. sign, EB-EB", nbinmee, meemin, meemax );
218  h_ele_mee_os_ebeb->Sumw2();
219  h_ele_mee_os_ebee = new TH1F( "h_ele_mee_os_ebee", "ele pairs invariant mass, opp. sign, EB-EE", nbinmee, meemin, meemax );
220  h_ele_mee_os_ebee->Sumw2();
221  h_ele_mee_os_eeee = new TH1F( "h_ele_mee_os_eeee", "ele pairs invariant mass, opp. sign, EE-EE", nbinmee, meemin, meemax );
222  h_ele_mee_os_eeee->Sumw2();
223  h_ele_mee_os_gg = new TH1F( "h_ele_mee_os_gg", "ele pairs invariant mass, opp. sign, good-good", nbinmee, meemin, meemax );
224  h_ele_mee_os_gg->Sumw2();
225  h_ele_mee_os_gb = new TH1F( "h_ele_mee_os_gb", "ele pairs invariant mass, opp. sign, good-bad", nbinmee, meemin, meemax );
226  h_ele_mee_os_gb->Sumw2();
227  h_ele_mee_os_bb = new TH1F( "h_ele_mee_os_bb", "ele pairs invariant mass, opp. sign, bad-bad", nbinmee, meemin, meemax );
228  h_ele_mee_os_bb->Sumw2();
229 
230  // duplicates
231  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. );
232  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. );
233 
234  // charge ID
235  h_ele_ChargeMnChargeTrue = new TH1F( "h_ele_ChargeMnChargeTrue", "ele charge - gen charge ",5,-1.,4.);
236  h_ele_ChargeMnChargeTrue->Sumw2();
237  h_ele_simEta_matched_qmisid = new TH1F( "h_ele_eta_matched_qmisid","charge misid vs gen eta", nbineta,etamin,etamax);
238  h_ele_simEta_matched_qmisid->Sumw2();
239  h_ele_simAbsEta_matched_qmisid = new TH1F( "h_ele_abseta_matched_qmisid", "charge misid vs gen |eta|", nbineta/2,0.,etamax);
240  h_ele_simAbsEta_matched_qmisid->Sumw2();
241  h_ele_simPt_matched_qmisid = new TH1F( "h_ele_Pt_matched_qmisid", "charge misid vs gen transverse momentum", nbinpteff,5.,ptmax);
242  h_ele_simPt_matched_qmisid->Sumw2();
243  h_ele_simPhi_matched_qmisid = new TH1F( "h_ele_phi_matched_qmisid","charge misid vs gen phi", nbinphi,phimin,phimax);
244  h_ele_simPhi_matched_qmisid->Sumw2();
245  h_ele_simZ_matched_qmisid = new TH1F( "h_ele_z_matched_qmisid","charge misid vs gen z",nbinxyz, -25, 25 );
246  h_ele_simZ_matched_qmisid->Sumw2();
247 
248  // matched electrons
249  h_ele_charge = new TH1F( "h_ele_charge", "ele charge", 5,-2.,2.);
250  h_ele_charge->Sumw2();
251  h_ele_chargeVsEta = new TH2F( "h_ele_chargeVsEta", "ele charge vs eta", nbineta2D,etamin,etamax,5,-2.,2.);
252  h_ele_chargeVsPhi = new TH2F( "h_ele_chargeVsPhi", "ele charge vs phi", nbinphi2D,phimin,phimax,5,-2.,2.);
253  h_ele_chargeVsPt = new TH2F( "h_ele_chargeVsPt", "ele charge vs pt", nbinpt,0.,100.,5,-2.,2.);
254  h_ele_vertexP = new TH1F( "h_ele_vertexP", "ele momentum", nbinp,0.,pmax);
255  h_ele_vertexP->Sumw2();
256  h_ele_vertexPt = new TH1F( "h_ele_vertexPt", "ele transverse momentum", nbinpt,0.,ptmax);
257  h_ele_vertexPt->Sumw2();
258  h_ele_Et = new TH1F( "h_ele_Et", "ele transverse energy", nbinpt,0.,ptmax);
259  h_ele_Et->Sumw2();
260  h_ele_vertexPtVsEta = new TH2F( "h_ele_vertexPtVsEta", "ele transverse momentum vs eta",nbineta2D,etamin,etamax,nbinpt2D,0.,ptmax);
261  h_ele_vertexPtVsPhi = new TH2F( "h_ele_vertexPtVsPhi", "ele transverse momentum vs phi",nbinphi2D,phimin,phimax,nbinpt2D,0.,ptmax);
262  h_ele_simPt_matched = new TH1F( "h_ele_simPt_matched", "Efficiency vs gen transverse momentum", nbinpteff,5.,ptmax);
263  h_ele_vertexEta = new TH1F( "h_ele_vertexEta", "ele momentum eta", nbineta,etamin,etamax);
264  h_ele_vertexEta->Sumw2();
265  h_ele_vertexEtaVsPhi = new TH2F( "h_ele_vertexEtaVsPhi", "ele momentum eta vs phi",nbineta2D,etamin,etamax,nbinphi2D,phimin,phimax );
266  h_ele_simAbsEta_matched = new TH1F( "h_ele_simAbsEta_matched", "Efficiency vs gen |eta|", nbineta/2,0.,2.5);
267  h_ele_simAbsEta_matched->Sumw2();
268  h_ele_simEta_matched = new TH1F( "h_ele_simEta_matched", "Efficiency vs gen eta", nbineta,etamin,etamax);
269  h_ele_simEta_matched->Sumw2();
270  h_ele_simPtEta_matched = new TH2F( "h_ele_simPtEta_matched", "Efficiency vs pt #eta", nbineta2D,etamin,etamax, nbinpt2D,5.,ptmax );
271  h_ele_simPtEta_matched->Sumw2();
272  h_ele_simPhi_matched = new TH1F( "h_ele_simPhi_matched", "Efficiency vs gen phi", nbinphi,phimin,phimax);
273  h_ele_simPhi_matched->Sumw2();
274  h_ele_vertexPhi = new TH1F( "h_ele_vertexPhi", "ele momentum #phi", nbinphi,phimin,phimax);
275  h_ele_vertexPhi->Sumw2();
276  h_ele_vertexX = new TH1F( "h_ele_vertexX", "ele vertex x", nbinxyz,-0.1,0.1 );
277  h_ele_vertexX->Sumw2();
278  h_ele_vertexY = new TH1F( "h_ele_vertexY", "ele vertex y", nbinxyz,-0.1,0.1 );
279  h_ele_vertexY->Sumw2();
280  h_ele_vertexZ = new TH1F( "h_ele_vertexZ", "ele vertex z", nbinxyz,-25, 25 );
281  h_ele_vertexZ->Sumw2();
282  h_ele_simZ_matched = new TH1F( "h_ele_simZ_matched", "Efficiency vs gen vertex z", nbinxyz,-25,25);
283  h_ele_simZ_matched->Sumw2();
284  h_ele_vertexTIP = new TH1F( "h_ele_vertexTIP", "ele transverse impact parameter (wrt gen vtx)", 90,0.,0.15);
285  h_ele_vertexTIP->Sumw2();
286  h_ele_vertexTIPVsEta = new TH2F( "h_ele_vertexTIPVsEta", "ele transverse impact parameter (wrt gen vtx) vs eta", nbineta2D,etamin,etamax,45,0.,0.15);
287  h_ele_vertexTIPVsPhi = new TH2F( "h_ele_vertexTIPVsPhi", "ele transverse impact parameter (wrt gen vtx) vs phi", nbinphi2D,phimin,phimax,45,0.,0.15);
288  h_ele_vertexTIPVsPt = new TH2F( "h_ele_vertexTIPVsPt", "ele transverse impact parameter (wrt gen vtx) vs transverse momentum", nbinpt2D,0.,ptmax,45,0.,0.15);
289  h_ele_PoPtrue = new TH1F( "h_ele_PoPtrue", "ele momentum / gen momentum", nbinpoptrue,poptruemin,poptruemax);
290  h_ele_PoPtrue->Sumw2();
291  h_ele_PtoPttrue = new TH1F( "h_ele_PtoPttrue", "ele transverse momentum / gen transverse momentum", nbinpoptrue,poptruemin,poptruemax);
292  h_ele_PtoPttrue->Sumw2();
293  h_ele_PoPtrueVsEta = new TH2F( "h_ele_PoPtrueVsEta", "ele momentum / gen momentum vs eta", nbineta2D,etamin,etamax,50,poptruemin,poptruemax);
294  h_ele_PoPtrueVsPhi = new TH2F( "h_ele_PoPtrueVsPhi", "ele momentum / gen momentum vs phi", nbinphi2D,phimin,phimax,50,poptruemin,poptruemax);
295  h_ele_PoPtrueVsPt = new TH2F( "h_ele_PoPtrueVsPt", "ele momentum / gen momentum vs eta", nbinpt2D,0.,ptmax,50,poptruemin,poptruemax);
296  h_ele_PoPtrue_barrel = new TH1F( "h_ele_PoPtrue_barrel", "ele momentum / gen momentum, barrel",nbinpoptrue,poptruemin,poptruemax);
297  h_ele_PoPtrue_barrel->Sumw2();
298  h_ele_PoPtrue_endcaps = new TH1F( "h_ele_PoPtrue_endcaps", "ele momentum / gen momentum, endcaps",nbinpoptrue,poptruemin,poptruemax);
299  h_ele_PoPtrue_endcaps->Sumw2();
300  h_ele_PoPtrue_golden_barrel = new TH1F( "h_ele_PoPtrue_golden_barrel", "ele momentum / gen momentum, golden, barrel",nbinpoptrue,poptruemin,poptruemax);
301  h_ele_PoPtrue_golden_barrel->Sumw2();
302  h_ele_PoPtrue_golden_endcaps = new TH1F( "h_ele_PoPtrue_golden_endcaps", "ele momentum / gen momentum, golden, endcaps",nbinpoptrue,poptruemin,poptruemax);
303  h_ele_PoPtrue_golden_endcaps->Sumw2();
304  h_ele_PoPtrue_showering_barrel = new TH1F( "h_ele_PoPtrue_showering_barrel", "ele momentum / gen momentum, showering, barrel",nbinpoptrue,poptruemin,poptruemax);
305  h_ele_PoPtrue_showering_barrel->Sumw2();
306  h_ele_PoPtrue_showering_endcaps = new TH1F( "h_ele_PoPtrue_showering_endcaps", "ele momentum / gen momentum, showering, endcaps",nbinpoptrue,poptruemin,poptruemax);
307  h_ele_PoPtrue_showering_endcaps->Sumw2();
308  h_ele_PtoPttrue_barrel = new TH1F( "h_ele_PtoPttrue_barrel", "ele transverse momentum / gen transverse momentum, barrel",nbinpoptrue,poptruemin,poptruemax);
309  h_ele_PtoPttrue_barrel->Sumw2();
310  h_ele_PtoPttrue_endcaps = new TH1F( "h_ele_PtoPttrue_endcaps", "ele transverse momentum / gen transverse momentum, endcaps",nbinpoptrue,poptruemin,poptruemax);
311  h_ele_PtoPttrue_endcaps->Sumw2();
312  h_ele_EtaMnEtaTrue = new TH1F( "h_ele_EtaMnEtaTrue", "ele momentum eta - gen eta",nbindeta,detamin,detamax);
313  h_ele_EtaMnEtaTrue->Sumw2();
314  h_ele_EtaMnEtaTrue_barrel = new TH1F( "h_ele_EtaMnEtaTrue_barrel", "ele momentum eta - gen eta barrel",nbindeta,detamin,detamax);
315  h_ele_EtaMnEtaTrue_barrel->Sumw2();
316  h_ele_EtaMnEtaTrue_endcaps = new TH1F( "h_ele_EtaMnEtaTrue_endcaps", "ele momentum eta - gen eta endcaps",nbindeta,detamin,detamax);
317  h_ele_EtaMnEtaTrue_endcaps->Sumw2();
318  h_ele_EtaMnEtaTrueVsEta = new TH2F( "h_ele_EtaMnEtaTrueVsEta", "ele momentum eta - gen eta vs eta",nbineta2D,etamin,etamax,nbindeta/2,detamin,detamax);
319  h_ele_EtaMnEtaTrueVsPhi = new TH2F( "h_ele_EtaMnEtaTrueVsPhi", "ele momentum eta - gen eta vs phi",nbinphi2D,phimin,phimax,nbindeta/2,detamin,detamax);
320  h_ele_EtaMnEtaTrueVsPt = new TH2F( "h_ele_EtaMnEtaTrueVsPt", "ele momentum eta - gen eta vs pt",nbinpt,0.,ptmax,nbindeta/2,detamin,detamax);
321  h_ele_PhiMnPhiTrue = new TH1F( "h_ele_PhiMnPhiTrue", "ele momentum phi - gen phi",nbindphi,dphimin,dphimax);
322  h_ele_PhiMnPhiTrue->Sumw2();
323  h_ele_PhiMnPhiTrue_barrel = new TH1F( "h_ele_PhiMnPhiTrue_barrel", "ele momentum phi - gen phi barrel",nbindphi,dphimin,dphimax);
324  h_ele_PhiMnPhiTrue_barrel->Sumw2();
325  h_ele_PhiMnPhiTrue_endcaps = new TH1F( "h_ele_PhiMnPhiTrue_endcaps", "ele momentum phi - gen phi endcaps",nbindphi,dphimin,dphimax);
326  h_ele_PhiMnPhiTrue_endcaps->Sumw2();
327  h_ele_PhiMnPhiTrue2 = new TH1F( "h_ele_PhiMnPhiTrue2", "ele momentum phi - gen phi",nbindphimatch2D,dphimatchmin,dphimatchmax);
328  h_ele_PhiMnPhiTrueVsEta = new TH2F( "h_ele_PhiMnPhiTrueVsEta", "ele momentum phi - gen phi vs eta",nbineta2D,etamin,etamax,nbindphi/2,dphimin,dphimax);
329  h_ele_PhiMnPhiTrueVsPhi = new TH2F( "h_ele_PhiMnPhiTrueVsPhi", "ele momentum phi - gen phi vs phi",nbinphi2D,phimin,phimax,nbindphi/2,dphimin,dphimax);
330  h_ele_PhiMnPhiTrueVsPt = new TH2F( "h_ele_PhiMnPhiTrueVsPt", "ele momentum phi - gen phi vs pt",nbinpt2D,0.,ptmax,nbindphi/2,dphimin,dphimax);
331 
332  // matched electron, superclusters
333  histSclEn_ = new TH1F("h_scl_energy","ele supercluster energy",nbinp,0.,pmax);
334  histSclEn_->Sumw2();
335  histSclEoEtrue_barrel = new TH1F("h_scl_EoEtrue_barrel","ele supercluster energy / gen energy, barrel",50,0.2,1.2);
336  histSclEoEtrue_barrel->Sumw2();
337  histSclEoEtrue_barrel_eg = new TH1F("h_scl_EoEtrue_barrel_eg","ele supercluster energy / gen energy, barrel, ecal driven",50,0.2,1.2);
338  histSclEoEtrue_barrel_eg->Sumw2();
339  histSclEoEtrue_barrel_etagap = new TH1F("h_scl_EoEtrue_barrel_etagap","ele supercluster energy / gen energy, barrel, etagap",50,0.2,1.2);
340  histSclEoEtrue_barrel_etagap->Sumw2();
341  histSclEoEtrue_barrel_phigap = new TH1F("h_scl_EoEtrue_barrel_phigap","ele supercluster energy / gen energy, barrel, phigap",50,0.2,1.2);
342  histSclEoEtrue_barrel_phigap->Sumw2();
343  histSclEoEtrue_ebeegap = new TH1F("h_scl_EoEtrue_ebeegap","ele supercluster energy / gen energy, ebeegap",50,0.2,1.2);
344  histSclEoEtrue_ebeegap->Sumw2();
345  histSclEoEtrue_endcaps = new TH1F("h_scl_EoEtrue_endcaps","ele supercluster energy / gen energy, endcaps",50,0.2,1.2);
346  histSclEoEtrue_endcaps->Sumw2();
347  histSclEoEtrue_endcaps_eg = new TH1F("h_scl_EoEtrue_endcaps_eg","ele supercluster energy / gen energy, endcaps, ecal driven",50,0.2,1.2);
348  histSclEoEtrue_endcaps_eg->Sumw2();
349  histSclEoEtrue_endcaps_deegap = new TH1F("h_scl_EoEtrue_endcaps_deegap","ele supercluster energy / gen energy, endcaps, deegap",50,0.2,1.2);
350  histSclEoEtrue_endcaps_deegap->Sumw2();
351  histSclEoEtrue_endcaps_ringgap = new TH1F("h_scl_EoEtrue_endcaps_ringgap","ele supercluster energy / gen energy, endcaps, ringgap",50,0.2,1.2);
352  histSclEoEtrue_endcaps_ringgap->Sumw2();
353  histSclEoEtrue_barrel_new = new TH1F("h_scl_EoEtrue_barrel_new","ele supercluster energy / gen energy, barrel",nbinpoptrue,poptruemin,poptruemax);
354  histSclEoEtrue_barrel_new->Sumw2();
355  histSclEoEtrue_barrel_eg_new = new TH1F("h_scl_EoEtrue_barrel_eg_new","ele supercluster energy / gen energy, barrel, ecal driven",nbinpoptrue,poptruemin,poptruemax);
356  histSclEoEtrue_barrel_eg_new->Sumw2();
357  histSclEoEtrue_barrel_etagap_new = new TH1F("h_scl_EoEtrue_barrel_etagap_new","ele supercluster energy / gen energy, barrel, etagap",nbinpoptrue,poptruemin,poptruemax);
358  histSclEoEtrue_barrel_etagap_new->Sumw2();
359  histSclEoEtrue_barrel_phigap_new = new TH1F("h_scl_EoEtrue_barrel_phigap_new","ele supercluster energy / gen energy, barrel, phigap",nbinpoptrue,poptruemin,poptruemax);
360  histSclEoEtrue_barrel_phigap_new->Sumw2();
361  histSclEoEtrue_ebeegap_new = new TH1F("h_scl_EoEtrue_ebeegap_new","ele supercluster energy / gen energy, ebeegap",nbinpoptrue,poptruemin,poptruemax);
362  histSclEoEtrue_ebeegap_new->Sumw2();
363  histSclEoEtrue_endcaps_new = new TH1F("h_scl_EoEtrue_endcaps_new","ele supercluster energy / gen energy, endcaps",nbinpoptrue,poptruemin,poptruemax);
364  histSclEoEtrue_endcaps_new->Sumw2();
365  histSclEoEtrue_endcaps_eg_new = new TH1F("h_scl_EoEtrue_endcaps_eg_new","ele supercluster energy / gen energy, endcaps, ecal driven",nbinpoptrue,poptruemin,poptruemax);
366  histSclEoEtrue_endcaps_eg_new->Sumw2();
367  histSclEoEtrue_endcaps_deegap_new = new TH1F("h_scl_EoEtrue_endcaps_deegap_new","ele supercluster energy / gen energy, endcaps, deegap",nbinpoptrue,poptruemin,poptruemax);
368  histSclEoEtrue_endcaps_deegap_new->Sumw2();
369  histSclEoEtrue_endcaps_ringgap_new = new TH1F("h_scl_EoEtrue_endcaps_ringgap_new","ele supercluster energy / gen energy, endcaps, ringgap",nbinpoptrue,poptruemin,poptruemax);
370  histSclEoEtrue_endcaps_ringgap_new->Sumw2();
371  histSclEt_ = new TH1F("h_scl_et","ele supercluster transverse energy",nbinpt,0.,ptmax);
372  histSclEt_->Sumw2();
373  histSclEtVsEta_ = new TH2F("h_scl_etVsEta","ele supercluster transverse energy vs eta",nbineta2D,etamin,etamax,nbinpt,0.,ptmax);
374  histSclEtVsPhi_ = new TH2F("h_scl_etVsPhi","ele supercluster transverse energy vs phi",nbinphi2D,phimin,phimax,nbinpt,0.,ptmax);
375  histSclEtaVsPhi_ = new TH2F("h_scl_etaVsPhi","ele supercluster eta vs phi",nbinphi2D,phimin,phimax,nbineta2D,etamin,etamax);
376  histSclEta_ = new TH1F("h_scl_eta","ele supercluster eta",nbineta,etamin,etamax);
377  histSclEta_->Sumw2();
378  histSclPhi_ = new TH1F("h_scl_phi","ele supercluster phi",nbinphi,phimin,phimax);
379  histSclPhi_->Sumw2();
380 
381  histSclSigEtaEta_ = new TH1F("h_scl_sigetaeta","ele supercluster sigma eta eta",100,0.,0.05);
382  histSclSigEtaEta_->Sumw2();
383  histSclSigEtaEta_barrel_ = new TH1F("h_scl_sigetaeta_barrel","ele supercluster sigma eta eta barrel",100,0.,0.05);
384  histSclSigEtaEta_barrel_->Sumw2();
385  histSclSigEtaEta_endcaps_ = new TH1F("h_scl_sigetaeta_endcaps","ele supercluster sigma eta eta endcaps",100,0.,0.05);
386  histSclSigEtaEta_endcaps_->Sumw2();
387  histSclSigIEtaIEta_ = new TH1F("h_scl_sigietaieta","ele supercluster sigma ieta ieta",100,0.,0.05);
388  histSclSigIEtaIEta_->Sumw2();
389  histSclSigIEtaIEta_barrel_ = new TH1F("h_scl_sigietaieta_barrel","ele supercluster sigma ieta ieta, barrel",100,0.,0.05);
390  histSclSigIEtaIEta_barrel_->Sumw2();
391  histSclSigIEtaIEta_endcaps_ = new TH1F("h_scl_sigietaieta_endcaps","ele supercluster sigma ieta ieta, endcaps",100,0.,0.05);
392  histSclSigIEtaIEta_endcaps_->Sumw2();
393  histSclE1x5_ = new TH1F("h_scl_E1x5","ele supercluster energy in 1x5",nbinp,0., pmax);
394  histSclE1x5_->Sumw2();
395  histSclE1x5_barrel_ = new TH1F("h_scl_E1x5_barrel","ele supercluster energy in 1x5 barrel",nbinp,0., pmax);
396  histSclE1x5_barrel_->Sumw2();
397  histSclE1x5_endcaps_ = new TH1F("h_scl_E1x5_endcaps","ele supercluster energy in 1x5 endcaps",nbinp,0., pmax);
398  histSclE1x5_endcaps_->Sumw2();
399  histSclE2x5max_ = new TH1F("h_scl_E2x5max","ele supercluster energy in 2x5 max",nbinp,0.,pmax);
400  histSclE2x5max_->Sumw2();
401  histSclE2x5max_barrel_ = new TH1F("h_scl_E2x5max_barrel","ele supercluster energy in 2x5 max barrel",nbinp,0.,pmax);
402  histSclE2x5max_barrel_->Sumw2();
403  histSclE2x5max_endcaps_ = new TH1F("h_scl_E2x5max_endcaps","ele supercluster energy in 2x5 max endcaps",nbinp,0.,pmax);
404  histSclE2x5max_endcaps_->Sumw2();
405  histSclE5x5_ = new TH1F("h_scl_E5x5","ele supercluster energy in 5x5",nbinp,0.,pmax);
406  histSclE5x5_->Sumw2();
407  histSclE5x5_barrel_ = new TH1F("h_scl_E5x5_barrel","ele supercluster energy in 5x5 barrel",nbinp,0.,pmax);
408  histSclE5x5_barrel_->Sumw2();
409  histSclE5x5_endcaps_ = new TH1F("h_scl_E5x5_endcaps","ele supercluster energy in 5x5 endcaps",nbinp,0.,pmax);
410  histSclE5x5_endcaps_->Sumw2();
411  histSclSigEtaEta_eg_ = new TH1F("h_scl_sigetaeta_eg","ele supercluster sigma eta eta, ecal driven",100,0.,0.05);
412  histSclSigEtaEta_eg_->Sumw2();
413  histSclSigEtaEta_eg_barrel_ = new TH1F("h_scl_sigetaeta_eg_barrel","ele supercluster sigma eta eta, ecal driven barrel",100,0.,0.05);
414  histSclSigEtaEta_eg_barrel_->Sumw2();
415  histSclSigEtaEta_eg_endcaps_ = new TH1F("h_scl_sigetaeta_eg_endcaps","ele supercluster sigma eta eta, ecal driven endcaps",100,0.,0.05);
416  histSclSigEtaEta_eg_endcaps_->Sumw2();
417  histSclSigIEtaIEta_eg_ = new TH1F("h_scl_sigietaieta_eg","ele supercluster sigma ieta ieta, ecal driven",100,0.,0.05);
418  histSclSigIEtaIEta_eg_->Sumw2();
419  histSclSigIEtaIEta_eg_barrel_ = new TH1F("h_scl_sigietaieta_barrel_eg","ele supercluster sigma ieta ieta, barrel, ecal driven",100,0.,0.05);
420  histSclSigIEtaIEta_eg_barrel_->Sumw2();
421  histSclSigIEtaIEta_eg_endcaps_ = new TH1F("h_scl_sigietaieta_endcaps_eg","ele supercluster sigma ieta ieta, endcaps, ecal driven",100,0.,0.05);
422  histSclSigIEtaIEta_eg_endcaps_->Sumw2();
423  histSclE1x5_eg_ = new TH1F("h_scl_E1x5_eg","ele supercluster energy in 1x5, ecal driven",nbinp,0., pmax);
424  histSclE1x5_eg_->Sumw2();
425  histSclE1x5_eg_barrel_ = new TH1F("h_scl_E1x5_eg_barrel","ele supercluster energy in 1x5, ecal driven barrel",nbinp,0., pmax);
426  histSclE1x5_eg_barrel_->Sumw2();
427  histSclE1x5_eg_endcaps_ = new TH1F("h_scl_E1x5_eg_endcaps","ele supercluster energy in 1x5, ecal driven endcaps",nbinp,0., pmax);
428  histSclE1x5_eg_endcaps_->Sumw2();
429  histSclE2x5max_eg_ = new TH1F("h_scl_E2x5max_eg","ele supercluster energy in 2x5 max, ecal driven",nbinp,0.,pmax);
430  histSclE2x5max_eg_->Sumw2();
431  histSclE2x5max_eg_barrel_ = new TH1F("h_scl_E2x5max_eg_barrel","ele supercluster energy in 2x5 max, ecal driven barrel",nbinp,0.,pmax);
432  histSclE2x5max_eg_barrel_->Sumw2();
433  histSclE2x5max_eg_endcaps_ = new TH1F("h_scl_E2x5max_eg_endcaps","ele supercluster energy in 2x5 max, ecal driven endcaps",nbinp,0.,pmax);
434  histSclE2x5max_eg_endcaps_->Sumw2();
435  histSclE5x5_eg_ = new TH1F("h_scl_E5x5_eg","ele supercluster energy in 5x5, ecal driven",nbinp,0.,pmax);
436  histSclE5x5_eg_->Sumw2();
437  histSclE5x5_eg_barrel_ = new TH1F("h_scl_E5x5_eg_barrel","ele supercluster energy in 5x5, ecal driven barrel",nbinp,0.,pmax);
438  histSclE5x5_eg_barrel_->Sumw2();
439  histSclE5x5_eg_endcaps_ = new TH1F("h_scl_E5x5_eg_endcaps","ele supercluster energy in 5x5, ecal driven endcaps",nbinp,0.,pmax);
440  histSclE5x5_eg_endcaps_->Sumw2();
441 
442  histSclEoEtruePfVsEg = new TH2F("h_scl_EoEtruePfVsEg","ele supercluster energy / gen energy pflow vs eg",75,-0.1,1.4, 75, -0.1, 1.4);
443 
444  // matched electron, gsf tracks
445  h_ele_ambiguousTracks = new TH1F( "h_ele_ambiguousTracks", "ele # ambiguous tracks", 5,0.,5.);
446  h_ele_ambiguousTracks->Sumw2();
447  h_ele_ambiguousTracksVsEta = new TH2F( "h_ele_ambiguousTracksVsEta","ele # ambiguous tracks vs eta", nbineta2D,etamin,etamax,5,0.,5.);
448  h_ele_ambiguousTracksVsPhi = new TH2F( "h_ele_ambiguousTracksVsPhi", "ele # ambiguous tracks vs phi", nbinphi2D,phimin,phimax,5,0.,5.);
449  h_ele_ambiguousTracksVsPt = new TH2F( "h_ele_ambiguousTracksVsPt", "ele # ambiguous tracks vs pt", nbinpt2D,0.,ptmax,5,0.,5.);
450  h_ele_foundHits = new TH1F( "h_ele_foundHits", "ele track # found hits", nbinfhits,0.,fhitsmax);
451  h_ele_foundHits->Sumw2();
452  h_ele_foundHits_barrel = new TH1F( "h_ele_foundHits_barrel", "ele track # found hits, barrel", nbinfhits,0.,fhitsmax);
453  h_ele_foundHits_barrel->Sumw2();
454  h_ele_foundHits_endcaps = new TH1F( "h_ele_foundHits_endcaps", "ele track # found hits, endcaps", nbinfhits,0.,fhitsmax);
455  h_ele_foundHits_endcaps->Sumw2();
456  h_ele_foundHitsVsEta = new TH2F( "h_ele_foundHitsVsEta", "ele track # found hits vs eta", nbineta2D,etamin,etamax,nbinfhits,0.,fhitsmax);
457  h_ele_foundHitsVsPhi = new TH2F( "h_ele_foundHitsVsPhi", "ele track # found hits vs phi", nbinphi2D,phimin,phimax,nbinfhits,0.,fhitsmax);
458  h_ele_foundHitsVsPt = new TH2F( "h_ele_foundHitsVsPt", "ele track # found hits vs pt", nbinpt2D,0.,ptmax,nbinfhits,0.,fhitsmax);
459  h_ele_lostHits = new TH1F( "h_ele_lostHits", "ele track # lost hits", 5,0.,5.);
460  h_ele_lostHits->Sumw2();
461  h_ele_lostHits_barrel = new TH1F( "h_ele_lostHits_barrel", "ele track # lost hits, barrel", 5,0.,5.);
462  h_ele_lostHits_barrel->Sumw2();
463  h_ele_lostHits_endcaps = new TH1F( "h_ele_lostHits_endcaps", "ele track # lost hits, endcaps", 5,0.,5.);
464  h_ele_lostHits_endcaps->Sumw2();
465  h_ele_lostHitsVsEta = new TH2F( "h_ele_lostHitsVsEta", "ele track # lost hits vs eta", nbineta2D,etamin,etamax,nbinlhits,0.,lhitsmax);
466  h_ele_lostHitsVsPhi = new TH2F( "h_ele_lostHitsVsPhi", "ele track # lost hits vs eta", nbinphi2D,phimin,phimax,nbinlhits,0.,lhitsmax);
467  h_ele_lostHitsVsPt = new TH2F( "h_ele_lostHitsVsPt", "ele track # lost hits vs eta", nbinpt2D,0.,ptmax,nbinlhits,0.,lhitsmax);
468  h_ele_chi2 = new TH1F( "h_ele_chi2", "ele track #chi^{2}", 100,0.,15.);
469  h_ele_chi2->Sumw2();
470  h_ele_chi2_barrel = new TH1F( "h_ele_chi2_barrel", "ele track #chi^{2}, barrel", 100,0.,15.);
471  h_ele_chi2_barrel->Sumw2();
472  h_ele_chi2_endcaps = new TH1F( "h_ele_chi2_endcaps", "ele track #chi^{2}, endcaps", 100,0.,15.);
473  h_ele_chi2_endcaps->Sumw2();
474  h_ele_chi2VsEta = new TH2F( "h_ele_chi2VsEta", "ele track #chi^{2} vs eta", nbineta2D,etamin,etamax,50,0.,15.);
475  h_ele_chi2VsPhi = new TH2F( "h_ele_chi2VsPhi", "ele track #chi^{2} vs phi", nbinphi2D,phimin,phimax,50,0.,15.);
476  h_ele_chi2VsPt = new TH2F( "h_ele_chi2VsPt", "ele track #chi^{2} vs pt", nbinpt2D,0.,ptmax,50,0.,15.);
477  h_ele_PinMnPout = new TH1F( "h_ele_PinMnPout", "ele track inner p - outer p, mean of GSF components" ,nbinp,0.,200.);
478  h_ele_PinMnPout->Sumw2();
479  h_ele_PinMnPout_mode = new TH1F( "h_ele_PinMnPout_mode", "ele track inner p - outer p, mode of GSF components" ,nbinp,0.,100.);
480  h_ele_PinMnPout_mode->Sumw2();
481  h_ele_PinMnPoutVsEta_mode = new TH2F( "h_ele_PinMnPoutVsEta_mode", "ele track inner p - outer p vs eta, mode of GSF components" ,nbineta2D, etamin,etamax,nbinp2D,0.,100.);
482  h_ele_PinMnPoutVsPhi_mode = new TH2F( "h_ele_PinMnPoutVsPhi_mode", "ele track inner p - outer p vs phi, mode of GSF components" ,nbinphi2D, phimin,phimax,nbinp2D,0.,100.);
483  h_ele_PinMnPoutVsPt_mode = new TH2F( "h_ele_PinMnPoutVsPt_mode", "ele track inner p - outer p vs pt, mode of GSF components" ,nbinpt2D, 0.,ptmax,nbinp2D,0.,100.);
484  h_ele_PinMnPoutVsE_mode = new TH2F( "h_ele_PinMnPoutVsE_mode", "ele track inner p - outer p vs E, mode of GSF components" ,nbinp2D, 0.,200.,nbinp2D,0.,100.);
485  h_ele_PinMnPoutVsChi2_mode = new TH2F( "h_ele_PinMnPoutVsChi2_mode", "ele track inner p - outer p vs track chi2, mode of GSF components" ,50, 0.,20.,nbinp2D,0.,100.);
486  h_ele_outerP = new TH1F( "h_ele_outerP", "ele track outer p, mean of GSF components", nbinp,0.,pmax);
487  h_ele_outerP->Sumw2();
488  h_ele_outerP_mode = new TH1F( "h_ele_outerP_mode", "ele track outer p, mode of GSF components", nbinp,0.,pmax);
489  h_ele_outerP_mode->Sumw2();
490  h_ele_outerPVsEta_mode = new TH2F( "h_ele_outerPVsEta_mode", "ele track outer p vs eta mode", nbineta2D,etamin,etamax,50,0.,pmax);
491  h_ele_outerPt = new TH1F( "h_ele_outerPt", "ele track outer p_{T}, mean of GSF components", nbinpt,0.,ptmax);
492  h_ele_outerPt->Sumw2();
493  h_ele_outerPt_mode = new TH1F( "h_ele_outerPt_mode", "ele track outer p_{T}, mode of GSF components", nbinpt,0.,ptmax);
494  h_ele_outerPt_mode->Sumw2();
495  h_ele_outerPtVsEta_mode = new TH2F( "h_ele_outerPtVsEta_mode", "ele track outer p_{T} vs eta, mode of GSF components", nbineta2D,etamin,etamax,nbinpt2D,0.,ptmax);
496  h_ele_outerPtVsPhi_mode = new TH2F( "h_ele_outerPtVsPhi_mode", "ele track outer p_{T} vs phi, mode of GSF components", nbinphi2D,phimin,phimax,nbinpt2D,0.,ptmax);
497  h_ele_outerPtVsPt_mode = new TH2F( "h_ele_outerPtVsPt_mode", "ele track outer p_{T} vs pt, mode of GSF components", nbinpt2D,0.,100.,nbinpt2D,0.,ptmax);
498 
499  // matched electrons, matching
500  h_ele_EoP = new TH1F( "h_ele_EoP", "ele E/P_{vertex}", nbineop,0.,eopmax);
501  h_ele_EoP->Sumw2();
502  h_ele_EoP_eg = new TH1F( "h_ele_EoP_eg", "ele E/P_{vertex}, ecal driven", nbineop,0.,eopmax);
503  h_ele_EoP_eg->Sumw2();
504  h_ele_EoP_barrel = new TH1F( "h_ele_EoP_barrel", "ele E/P_{vertex} barrel", nbineop,0.,eopmax);
505  h_ele_EoP_barrel->Sumw2();
506  h_ele_EoP_eg_barrel = new TH1F( "h_ele_EoP_eg_barrel", "ele E/P_{vertex}, ecal driven barrel", nbineop,0.,eopmax);
507  h_ele_EoP_eg_barrel->Sumw2();
508  h_ele_EoP_endcaps = new TH1F( "h_ele_EoP_endcaps", "ele E/P_{vertex} endcaps", nbineop,0.,eopmax);
509  h_ele_EoP_endcaps->Sumw2();
510  h_ele_EoP_eg_endcaps = new TH1F( "h_ele_EoP_eg_endcaps", "ele E/P_{vertex}, ecal driven endcaps", nbineop,0.,eopmax);
511  h_ele_EoP_eg_endcaps->Sumw2();
512  h_ele_EoPVsEta = new TH2F( "h_ele_EoPVsEta", "ele E/P_{vertex} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
513  h_ele_EoPVsPhi = new TH2F( "h_ele_EoPVsPhi", "ele E/P_{vertex} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
514  h_ele_EoPVsE = new TH2F( "h_ele_EoPVsE", "ele E/P_{vertex} vs E", 50,0.,pmax ,50,0.,5.);
515  h_ele_EseedOP = new TH1F( "h_ele_EseedOP", "ele E_{seed}/P_{vertex}", nbineop,0.,eopmax);
516  h_ele_EseedOP->Sumw2();
517  h_ele_EseedOP_eg = new TH1F( "h_ele_EseedOP_eg", "ele E_{seed}/P_{vertex}, ecal driven", nbineop,0.,eopmax);
518  h_ele_EseedOP_eg->Sumw2();
519  h_ele_EseedOP_barrel = new TH1F( "h_ele_EseedOP_barrel", "ele E_{seed}/P_{vertex} barrel", nbineop,0.,eopmax);
520  h_ele_EseedOP_barrel->Sumw2();
521  h_ele_EseedOP_eg_barrel = new TH1F( "h_ele_EseedOP_eg_barrel", "ele E_{seed}/P_{vertex}, ecal driven barrel", nbineop,0.,eopmax);
522  h_ele_EseedOP_eg_barrel->Sumw2();
523  h_ele_EseedOP_endcaps = new TH1F( "h_ele_EseedOP_endcaps", "ele E_{seed}/P_{vertex} endcaps", nbineop,0.,eopmax);
524  h_ele_EseedOP_endcaps->Sumw2();
525  h_ele_EseedOP_eg_endcaps = new TH1F( "h_ele_EseedOP_eg_endcaps", "ele E_{seed}/P_{vertex}, ecal driven, endcaps", nbineop,0.,eopmax);
526  h_ele_EseedOP_eg_endcaps->Sumw2();
527  h_ele_EseedOPVsEta = new TH2F( "h_ele_EseedOPVsEta", "ele E_{seed}/P_{vertex} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
528  h_ele_EseedOPVsPhi = new TH2F( "h_ele_EseedOPVsPhi", "ele E_{seed}/P_{vertex} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
529  h_ele_EseedOPVsE = new TH2F( "h_ele_EseedOPVsE", "ele E_{seed}/P_{vertex} vs E", 50,0.,pmax ,50,0.,5.);
530  h_ele_EoPout = new TH1F( "h_ele_EoPout", "ele E_{seed}/P_{out}", nbineop,0.,eopmax);
531  h_ele_EoPout->Sumw2();
532  h_ele_EoPout_eg = new TH1F( "h_ele_EoPout_eg", "ele E_{seed}/P_{out}, ecal driven", nbineop,0.,eopmax);
533  h_ele_EoPout_eg->Sumw2();
534  h_ele_EoPout_barrel = new TH1F( "h_ele_EoPout_barrel", "ele E_{seed}/P_{out} barrel", nbineop,0.,eopmax);
535  h_ele_EoPout_barrel->Sumw2();
536  h_ele_EoPout_eg_barrel = new TH1F( "h_ele_EoPout_eg_barrel", "ele E_{seed}/P_{out}, ecal driven, barrel", nbineop,0.,eopmax);
537  h_ele_EoPout_eg_barrel->Sumw2();
538  h_ele_EoPout_endcaps = new TH1F( "h_ele_EoPout_endcaps", "ele E_{seed}/P_{out} endcaps", nbineop,0.,eopmax);
539  h_ele_EoPout_endcaps->Sumw2();
540  h_ele_EoPout_eg_endcaps = new TH1F( "h_ele_EoPout_eg_endcaps", "ele E_{seed}/P_{out}, ecal driven, endcaps", nbineop,0.,eopmax);
541  h_ele_EoPout_eg_endcaps->Sumw2();
542  h_ele_EoPoutVsEta = new TH2F( "h_ele_EoPoutVsEta", "ele E_{seed}/P_{out} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
543  h_ele_EoPoutVsPhi = new TH2F( "h_ele_EoPoutVsPhi", "ele E_{seed}/P_{out} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
544  h_ele_EoPoutVsE = new TH2F( "h_ele_EoPoutVsE", "ele E_{seed}/P_{out} vs E", nbinp2D,0.,pmax,nbineop2D,0.,eopmaxsht);
545  h_ele_EeleOPout = new TH1F( "h_ele_EeleOPout", "ele E_{ele}/P_{out}", nbineop,0.,eopmax);
546  h_ele_EeleOPout->Sumw2();
547  h_ele_EeleOPout_eg = new TH1F( "h_ele_EeleOPout_eg", "ele E_{ele}/P_{out}, ecal driven", nbineop,0.,eopmax);
548  h_ele_EeleOPout_eg->Sumw2();
549  h_ele_EeleOPout_barrel = new TH1F( "h_ele_EeleOPout_barrel", "ele E_{ele}/P_{out} barrel", nbineop,0.,eopmax);
550  h_ele_EeleOPout_barrel->Sumw2();
551  h_ele_EeleOPout_eg_barrel = new TH1F( "h_ele_EeleOPout_eg_barrel", "ele E_{ele}/P_{out}, ecal driven, barrel", nbineop,0.,eopmax);
552  h_ele_EeleOPout_eg_barrel->Sumw2();
553  h_ele_EeleOPout_endcaps = new TH1F( "h_ele_EeleOPout_endcaps", "ele E_{ele}/P_{out} endcaps", nbineop,0.,eopmax);
554  h_ele_EeleOPout_endcaps->Sumw2();
555  h_ele_EeleOPout_eg_endcaps = new TH1F( "h_ele_EeleOPout_eg_endcaps", "ele E_{ele}/P_{out}, ecal driven, endcaps", nbineop,0.,eopmax);
556  h_ele_EeleOPout_eg_endcaps->Sumw2();
557  h_ele_EeleOPoutVsEta = new TH2F( "h_ele_EeleOPoutVsEta", "ele E_{ele}/P_{out} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
558  h_ele_EeleOPoutVsPhi = new TH2F( "h_ele_EeleOPoutVsPhi", "ele E_{ele}/P_{out} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
559  h_ele_EeleOPoutVsE = new TH2F( "h_ele_EeleOPoutVsE", "ele E_{ele}/P_{out} vs E", nbinp2D,0.,pmax,nbineop2D,0.,eopmaxsht);
560  h_ele_dEtaSc_propVtx = new TH1F( "h_ele_dEtaSc_propVtx", "ele #eta_{sc} - #eta_{tr}, prop from vertex", nbindetamatch,detamatchmin,detamatchmax);
561  h_ele_dEtaSc_propVtx->Sumw2();
562  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);
563  h_ele_dEtaSc_propVtx_eg->Sumw2();
564  h_ele_dEtaSc_propVtx_barrel = new TH1F( "h_ele_dEtaSc_propVtx_barrel", "ele #eta_{sc} - #eta_{tr}, prop from vertex, barrel", nbindetamatch,detamatchmin,detamatchmax);
565  h_ele_dEtaSc_propVtx_barrel->Sumw2();
566  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);
567  h_ele_dEtaSc_propVtx_eg_barrel->Sumw2();
568  h_ele_dEtaSc_propVtx_endcaps = new TH1F( "h_ele_dEtaSc_propVtx_endcaps", "ele #eta_{sc} - #eta_{tr}, prop from vertex, endcaps", nbindetamatch,detamatchmin,detamatchmax);
569  h_ele_dEtaSc_propVtx_endcaps->Sumw2();
570  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);
571  h_ele_dEtaSc_propVtx_eg_endcaps->Sumw2();
572  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);
573  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);
574  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);
575  h_ele_dPhiSc_propVtx = new TH1F( "h_ele_dPhiSc_propVtx", "ele #phi_{sc} - #phi_{tr}, prop from vertex", nbindphimatch,dphimatchmin,dphimatchmax);
576  h_ele_dPhiSc_propVtx->Sumw2();
577  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);
578  h_ele_dPhiSc_propVtx_eg->Sumw2();
579  h_ele_dPhiSc_propVtx_barrel = new TH1F( "h_ele_dPhiSc_propVtx_barrel", "ele #phi_{sc} - #phi_{tr}, prop from vertex, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
580  h_ele_dPhiSc_propVtx_barrel->Sumw2();
581  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);
582  h_ele_dPhiSc_propVtx_eg_barrel->Sumw2();
583  h_ele_dPhiSc_propVtx_endcaps = new TH1F( "h_ele_dPhiSc_propVtx_endcaps", "ele #phi_{sc} - #phi_{tr}, prop from vertex, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
584  h_ele_dPhiSc_propVtx_endcaps->Sumw2();
585  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);
586  h_ele_dPhiSc_propVtx_eg_endcaps->Sumw2();
587  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);
588  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);
589  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);
590  h_ele_dEtaCl_propOut = new TH1F( "h_ele_dEtaCl_propOut", "ele #eta_{cl} - #eta_{tr}, prop from outermost", nbindetamatch,detamatchmin,detamatchmax);
591  h_ele_dEtaCl_propOut->Sumw2();
592  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);
593  h_ele_dEtaCl_propOut_eg->Sumw2();
594  h_ele_dEtaCl_propOut_barrel = new TH1F( "h_ele_dEtaCl_propOut_barrel", "ele #eta_{cl} - #eta_{tr}, prop from outermost, barrel", nbindetamatch,detamatchmin,detamatchmax);
595  h_ele_dEtaCl_propOut_barrel->Sumw2();
596  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);
597  h_ele_dEtaCl_propOut_eg_barrel->Sumw2();
598  h_ele_dEtaCl_propOut_endcaps = new TH1F( "h_ele_dEtaCl_propOut_endcaps", "ele #eta_{cl} - #eta_{tr}, prop from outermost, endcaps", nbindetamatch,detamatchmin,detamatchmax);
599  h_ele_dEtaCl_propOut_endcaps->Sumw2();
600  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);
601  h_ele_dEtaCl_propOut_eg_endcaps->Sumw2();
602  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);
603  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);
604  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);
605  h_ele_dPhiCl_propOut = new TH1F( "h_ele_dPhiCl_propOut", "ele #phi_{cl} - #phi_{tr}, prop from outermost", nbindphimatch,dphimatchmin,dphimatchmax);
606  h_ele_dPhiCl_propOut->Sumw2();
607  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);
608  h_ele_dPhiCl_propOut_eg->Sumw2();
609  h_ele_dPhiCl_propOut_barrel = new TH1F( "h_ele_dPhiCl_propOut_barrel", "ele #phi_{cl} - #phi_{tr}, prop from outermost, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
610  h_ele_dPhiCl_propOut_barrel->Sumw2();
611  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);
612  h_ele_dPhiCl_propOut_eg_barrel->Sumw2();
613  h_ele_dPhiCl_propOut_endcaps = new TH1F( "h_ele_dPhiCl_propOut_endcaps", "ele #phi_{cl} - #phi_{tr}, prop from outermost, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
614  h_ele_dPhiCl_propOut_endcaps->Sumw2();
615  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);
616  h_ele_dPhiCl_propOut_eg_endcaps->Sumw2();
617  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);
618  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);
619  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);
620  h_ele_dEtaEleCl_propOut = new TH1F( "h_ele_dEtaEleCl_propOut", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost", nbindetamatch,detamatchmin,detamatchmax);
621  h_ele_dEtaEleCl_propOut->Sumw2();
622  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);
623  h_ele_dEtaEleCl_propOut_eg->Sumw2();
624  h_ele_dEtaEleCl_propOut_barrel = new TH1F( "h_ele_dEtaEleCl_propOut_barrel", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, barrel", nbindetamatch,detamatchmin,detamatchmax);
625  h_ele_dEtaEleCl_propOut_barrel->Sumw2();
626  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);
627  h_ele_dEtaEleCl_propOut_eg_barrel->Sumw2();
628  h_ele_dEtaEleCl_propOut_endcaps = new TH1F( "h_ele_dEtaEleCl_propOut_endcaps", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, endcaps", nbindetamatch,detamatchmin,detamatchmax);
629  h_ele_dEtaEleCl_propOut_endcaps->Sumw2();
630  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);
631  h_ele_dEtaEleCl_propOut_eg_endcaps->Sumw2();
632  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);
633  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);
634  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);
635  h_ele_dPhiEleCl_propOut = new TH1F( "h_ele_dPhiEleCl_propOut", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost", nbindphimatch,dphimatchmin,dphimatchmax);
636  h_ele_dPhiEleCl_propOut->Sumw2();
637  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);
638  h_ele_dPhiEleCl_propOut_eg->Sumw2();
639  h_ele_dPhiEleCl_propOut_barrel = new TH1F( "h_ele_dPhiEleCl_propOut_barrel", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
640  h_ele_dPhiEleCl_propOut_barrel->Sumw2();
641  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);
642  h_ele_dPhiEleCl_propOut_eg_barrel->Sumw2();
643  h_ele_dPhiEleCl_propOut_endcaps = new TH1F( "h_ele_dPhiEleCl_propOut_endcaps", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
644  h_ele_dPhiEleCl_propOut_endcaps->Sumw2();
645  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);
646  h_ele_dPhiEleCl_propOut_eg_endcaps->Sumw2();
647  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);
648  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);
649  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);
650 
651  h_ele_HoE = new TH1F("h_ele_HoE", "ele hadronic energy / em energy", nbinhoe, hoemin, hoemax) ;
652  h_ele_HoE->Sumw2();
653  h_ele_HoE_eg = new TH1F("h_ele_HoE_eg", "ele hadronic energy / em energy, ecal driven", nbinhoe, hoemin, hoemax) ;
654  h_ele_HoE_eg->Sumw2();
655  h_ele_HoE_barrel = new TH1F("h_ele_HoE_barrel", "ele hadronic energy / em energy, barrel", nbinhoe, hoemin, hoemax) ;
656  h_ele_HoE_barrel->Sumw2();
657  h_ele_HoE_eg_barrel = new TH1F("h_ele_HoE_eg_barrel", "ele hadronic energy / em energy, ecal driven, barrel", nbinhoe, hoemin, hoemax) ;
658  h_ele_HoE_eg_barrel->Sumw2();
659  h_ele_HoE_endcaps = new TH1F("h_ele_HoE_endcaps", "ele hadronic energy / em energy, endcaps", nbinhoe, hoemin, hoemax) ;
660  h_ele_HoE_endcaps->Sumw2();
661  h_ele_HoE_eg_endcaps = new TH1F("h_ele_HoE_eg_endcaps", "ele hadronic energy / em energy, ecal driven, endcaps", nbinhoe, hoemin, hoemax) ;
662  h_ele_HoE_eg_endcaps->Sumw2();
663  h_ele_HoE_fiducial = new TH1F("h_ele_HoE_fiducial", "ele hadronic energy / em energy, fiducial region", nbinhoe, hoemin, hoemax) ;
664  h_ele_HoE_fiducial->Sumw2();
665  h_ele_HoEVsEta = new TH2F("h_ele_HoEVsEta", "ele hadronic energy / em energy vs eta", nbineta,etamin,etamax,nbinhoe, hoemin, hoemax) ;
666  h_ele_HoEVsPhi = new TH2F("h_ele_HoEVsPhi", "ele hadronic energy / em energy vs phi", nbinphi2D,phimin,phimax,nbinhoe, hoemin, hoemax) ;
667  h_ele_HoEVsE = new TH2F("h_ele_HoEVsE", "ele hadronic energy / em energy vs E", nbinp, 0.,300.,nbinhoe, hoemin, hoemax) ;
668 
669  h_ele_seed_dphi2_ = new TH1F("h_ele_seedDphi2", "ele seed dphi 2nd layer", 50,-0.003,+0.003) ;
670  h_ele_seed_dphi2_->Sumw2();
671  h_ele_seed_dphi2VsEta_ = new TH2F("h_ele_seedDphi2VsEta", "ele seed dphi 2nd layer vs eta", nbineta2D,etamin,etamax,50,-0.003,+0.003) ;
672  h_ele_seed_dphi2VsPt_ = new TH2F("h_ele_seedDphi2VsPt", "ele seed dphi 2nd layer vs pt", nbinpt2D,0.,ptmax,50,-0.003,+0.003) ;
673  h_ele_seed_drz2_ = new TH1F("h_ele_seedDrz2", "ele seed dr (dz) 2nd layer", 50,-0.03,+0.03) ;
674  h_ele_seed_drz2_->Sumw2();
675  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) ;
676  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) ;
677  h_ele_seed_subdet2_ = new TH1F("h_ele_seedSubdet2", "ele seed subdet 2nd layer", 10,0.,10.) ;
678  h_ele_seed_subdet2_->Sumw2();
679 
680  // classes
681  h_ele_classes = new TH1F( "h_ele_classes", "ele classes", 20,0.0,20.);
682  h_ele_classes->Sumw2();
683  h_ele_eta = new TH1F( "h_ele_eta", "ele electron eta", nbineta/2,0.0,etamax);
684  h_ele_eta->Sumw2();
685  h_ele_eta_golden = new TH1F( "h_ele_eta_golden", "ele electron eta golden", nbineta/2,0.0,etamax);
686  h_ele_eta_golden->Sumw2();
687  h_ele_eta_bbrem = new TH1F( "h_ele_eta_bbrem", "ele electron eta bbrem", nbineta/2,0.0,etamax);
688  h_ele_eta_bbrem->Sumw2();
689  h_ele_eta_narrow = new TH1F( "h_ele_eta_narrow", "ele electron eta narrow", nbineta/2,0.0,etamax);
690  h_ele_eta_narrow->Sumw2();
691  h_ele_eta_shower = new TH1F( "h_ele_eta_show", "ele electron eta showering", nbineta/2,0.0,etamax);
692  h_ele_eta_shower->Sumw2();
693  h_ele_PinVsPoutGolden_mode = new TH2F( "h_ele_PinVsPoutGolden_mode", "ele track inner p vs outer p vs eta, golden, mode of GSF components" ,nbinp2D,0.,pmax,50,0.,pmax);
694  h_ele_PinVsPoutShowering_mode = new TH2F( "h_ele_PinVsPoutShowering_mode", "ele track inner p vs outer p vs eta, showering, mode of GSF components" ,nbinp2D,0.,pmax,50,0.,pmax);
695  h_ele_PinVsPoutGolden_mean = new TH2F( "h_ele_PinVsPoutGolden_mean", "ele track inner p vs outer p vs eta, golden, mean of GSF components" ,nbinp2D,0.,pmax,50,0.,pmax);
696  h_ele_PinVsPoutShowering_mean = new TH2F( "h_ele_PinVsPoutShowering_mean", "ele track inner p vs outer p vs eta, showering, mean of GSF components" ,nbinp2D,0.,pmax,50,0.,pmax);
697  h_ele_PtinVsPtoutGolden_mode = new TH2F( "h_ele_PtinVsPtoutGolden_mode", "ele track inner pt vs outer pt vs eta, golden, mode of GSF components" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
698  h_ele_PtinVsPtoutShowering_mode = new TH2F( "h_ele_PtinVsPtoutShowering_mode", "ele track inner pt vs outer pt vs eta, showering, mode of GSF components" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
699  h_ele_PtinVsPtoutGolden_mean = new TH2F( "h_ele_PtinVsPtoutGolden_mean", "ele track inner pt vs outer pt vs eta, golden, mean of GSF components" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
700  h_ele_PtinVsPtoutShowering_mean = new TH2F( "h_ele_PtinVsPtoutShowering_mean", "ele track inner pt vs outer pt vs eta, showering, mean of GSF components" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
701  histSclEoEtrueGolden_barrel = new TH1F("h_scl_EoEtrue_golden_barrel","ele supercluster energy / gen energy, golden, barrel",nbinpoptrue,poptruemin,poptruemax);
702  histSclEoEtrueGolden_barrel->Sumw2();
703  histSclEoEtrueGolden_endcaps = new TH1F("h_scl_EoEtrue_golden_endcaps","ele supercluster energy / gen energy, golden, endcaps",nbinpoptrue,poptruemin,poptruemax);
704  histSclEoEtrueGolden_endcaps->Sumw2();
705  histSclEoEtrueShowering_barrel = new TH1F("h_scl_EoEtrue_showering_barrel","ele supercluster energy / gen energy, showering, barrel",nbinpoptrue,poptruemin,poptruemax);
706  histSclEoEtrueShowering_barrel->Sumw2();
707  histSclEoEtrueShowering_endcaps = new TH1F("h_scl_EoEtrue_showering_endcaps","ele supercluster energy / gen energy, showering, endcaps",nbinpoptrue,poptruemin,poptruemax);
708  histSclEoEtrueShowering_endcaps->Sumw2();
709 
710  // isolation
711  h_ele_tkSumPt_dr03 = new TH1F("h_ele_tkSumPt_dr03","tk isolation sum, dR=0.3",100,0.0,20.);
712  h_ele_tkSumPt_dr03->Sumw2();
713  h_ele_ecalRecHitSumEt_dr03= new TH1F("h_ele_ecalRecHitSumEt_dr03","ecal isolation sum, dR=0.3",100,0.0,20.);
714  h_ele_ecalRecHitSumEt_dr03->Sumw2();
715  h_ele_hcalDepth1TowerSumEt_dr03= new TH1F("h_ele_hcalDepth1TowerSumEt_dr03","hcal depth1 isolation sum, dR=0.3",100,0.0,20.);
716  h_ele_hcalDepth1TowerSumEt_dr03->Sumw2();
717  h_ele_hcalDepth2TowerSumEt_dr03= new TH1F("h_ele_hcalDepth2TowerSumEt_dr03","hcal depth2 isolation sum, dR=0.3",100,0.0,20.);
718  h_ele_hcalDepth2TowerSumEt_dr03->Sumw2();
719  h_ele_tkSumPt_dr04= new TH1F("h_ele_tkSumPt_dr04","tk isolation sum, dR=0.4",100,0.0,20.);
720  h_ele_tkSumPt_dr04->Sumw2();
721  h_ele_ecalRecHitSumEt_dr04= new TH1F("h_ele_ecalRecHitSumEt_dr04","ecal isolation sum, dR=0.4",100,0.0,20.);
722  h_ele_ecalRecHitSumEt_dr04->Sumw2();
723  h_ele_hcalDepth1TowerSumEt_dr04= new TH1F("h_ele_hcalDepth1TowerSumEt_dr04","hcal depth1 isolation sum, dR=0.4",100,0.0,20.);
724  h_ele_hcalDepth1TowerSumEt_dr04->Sumw2();
725  h_ele_hcalDepth2TowerSumEt_dr04= new TH1F("h_ele_hcalDepth2TowerSumEt_dr04","hcal depth2 isolation sum, dR=0.4",100,0.0,20.);
726  h_ele_hcalDepth2TowerSumEt_dr04->Sumw2();
727 
728  // fbrem
729  h_ele_fbrem = new TH1F( "h_ele_fbrem","ele brem fraction, mode of GSF components",100,0.,1.);
730  h_ele_fbrem->Sumw2();
731  h_ele_fbrem_eg = new TH1F( "h_ele_fbrem_eg","ele brem fraction, mode of GSF components, ecal driven",100,0.,1.);
732  h_ele_fbrem_eg->Sumw2();
733  h_ele_fbremVsEta_mode = new TProfile( "h_ele_fbremvsEtamode","mean ele brem fraction vs eta, mode of GSF components",nbineta2D,etamin,etamax,0.,1.);
734  h_ele_fbremVsEta_mean = new TProfile( "h_ele_fbremvsEtamean","mean ele brem fraction vs eta, mean of GSF components",nbineta2D,etamin,etamax,0.,1.);
735 
736  // e/g et pflow electrons
737  h_ele_mva = new TH1F( "h_ele_mva","ele identification mva",100,-1.,1.);
738  h_ele_mva->Sumw2();
739  h_ele_mva_eg = new TH1F( "h_ele_mva_eg","ele identification mva, ecal driven",100,-1.,1.);
740  h_ele_mva_eg->Sumw2();
741  h_ele_provenance = new TH1F( "h_ele_provenance","ele provenance",5,-2.,3.);
742  h_ele_provenance->Sumw2();
743 
744  // histos titles
745  h_mcNum -> GetXaxis()-> SetTitle("N_{gen}");
746  h_mcNum -> GetYaxis()-> SetTitle("Events");
747  h_eleNum -> GetXaxis()-> SetTitle("# gen ele");
748  h_eleNum -> GetYaxis()-> SetTitle("Events");
749  h_gamNum -> GetXaxis()-> SetTitle("N_{gen #gamma}");
750  h_gamNum -> GetYaxis()-> SetTitle("Events");
751  h_simEta -> GetXaxis()-> SetTitle("#eta");
752  h_simEta -> GetYaxis()-> SetTitle("Events");
753  h_simP -> GetXaxis()-> SetTitle("p (GeV/c)");
754  h_simP -> GetYaxis()-> SetTitle("Events");
755  h_ele_foundHits -> GetXaxis()-> SetTitle("N_{hits}");
756  h_ele_foundHits -> GetYaxis()-> SetTitle("Events");
757  h_ele_foundHits_barrel -> GetXaxis()-> SetTitle("N_{hits}");
758  h_ele_foundHits_barrel -> GetYaxis()-> SetTitle("Events");
759  h_ele_foundHits_endcaps -> GetXaxis()-> SetTitle("N_{hits}");
760  h_ele_foundHits_endcaps -> GetYaxis()-> SetTitle("Events");
761  h_ele_ambiguousTracks -> GetXaxis()-> SetTitle("N_{ambiguous tracks}");
762  h_ele_ambiguousTracks -> GetYaxis()-> SetTitle("Events");
763  h_ele_lostHits -> GetXaxis()-> SetTitle("N_{lost hits}");
764  h_ele_lostHits -> GetYaxis()-> SetTitle("Events");
765  h_ele_lostHits_barrel -> GetXaxis()-> SetTitle("N_{lost hits}");
766  h_ele_lostHits_barrel -> GetYaxis()-> SetTitle("Events");
767  h_ele_lostHits_endcaps -> GetXaxis()-> SetTitle("N_{lost hits}");
768  h_ele_lostHits_endcaps -> GetYaxis()-> SetTitle("Events");
769  h_ele_chi2 -> GetXaxis()-> SetTitle("#Chi^{2}");
770  h_ele_chi2 -> GetYaxis()-> SetTitle("Events");
771  h_ele_chi2_barrel -> GetXaxis()-> SetTitle("#Chi^{2}");
772  h_ele_chi2_barrel -> GetYaxis()-> SetTitle("Events");
773  h_ele_chi2_endcaps -> GetXaxis()-> SetTitle("#Chi^{2}");
774  h_ele_chi2_endcaps -> GetYaxis()-> SetTitle("Events");
775  h_ele_charge -> GetXaxis()-> SetTitle("charge");
776  h_ele_charge -> GetYaxis()-> SetTitle("Events");
777  h_ele_vertexP -> GetXaxis()-> SetTitle("p_{vertex} (GeV/c)");
778  h_ele_vertexP -> GetYaxis()-> SetTitle("Events");
779  h_ele_vertexPt -> GetXaxis()-> SetTitle("p_{T vertex} (GeV/c)");
780  h_ele_vertexPt -> GetYaxis()-> SetTitle("Events");
781  h_ele_Et -> GetXaxis()-> SetTitle("E_{T} (GeV)");
782  h_ele_Et -> GetYaxis()-> SetTitle("Events");
783  h_ele_Et_all -> GetXaxis()-> SetTitle("E_{T} (GeV)");
784  h_ele_Et_all -> GetYaxis()-> SetTitle("Events");
785  h_ele_vertexEta -> GetXaxis()-> SetTitle("#eta");
786  h_ele_vertexEta -> GetYaxis()-> SetTitle("Events");
787  h_ele_vertexPhi -> GetXaxis()-> SetTitle("#phi (rad)");
788  h_ele_vertexPhi -> GetYaxis()-> SetTitle("Events");
789  h_ele_PoPtrue -> GetXaxis()-> SetTitle("P/P_{gen}");
790  h_ele_PoPtrue -> GetYaxis()-> SetTitle("Events");
791  h_ele_PoPtrue_barrel -> GetXaxis()-> SetTitle("P/P_{gen}");
792  h_ele_PoPtrue_barrel -> GetYaxis()-> SetTitle("Events");
793  h_ele_PoPtrue_endcaps -> GetXaxis()-> SetTitle("P/P_{gen}");
794  h_ele_PoPtrue_endcaps -> GetYaxis()-> SetTitle("Events");
795  h_ele_PoPtrue_golden_barrel -> GetXaxis()-> SetTitle("P/P_{gen}");
796  h_ele_PoPtrue_golden_barrel -> GetYaxis()-> SetTitle("Events");
797  h_ele_PoPtrue_showering_barrel -> GetXaxis()-> SetTitle("P/P_{gen}");
798  h_ele_PoPtrue_showering_barrel -> GetYaxis()-> SetTitle("Events");
799  h_ele_PoPtrue_golden_endcaps -> GetXaxis()-> SetTitle("P/P_{gen}");
800  h_ele_PoPtrue_golden_endcaps -> GetYaxis()-> SetTitle("Events");
801  h_ele_PoPtrue_showering_endcaps -> GetXaxis()-> SetTitle("P/P_{gen}");
802  h_ele_PoPtrue_showering_endcaps -> GetYaxis()-> SetTitle("Events");
803  h_ele_PtoPttrue -> GetXaxis()-> SetTitle("P_{T}/P_{T}^{gen}");
804  h_ele_PtoPttrue -> GetYaxis()-> SetTitle("Events");
805  h_ele_PtoPttrue_barrel -> GetXaxis()-> SetTitle("P_{T}/P_{T}^{gen}");
806  h_ele_PtoPttrue_barrel -> GetYaxis()-> SetTitle("Events");
807  h_ele_PtoPttrue_endcaps -> GetXaxis()-> SetTitle("P_{T}/P_{T}^{gen}");
808  h_ele_PtoPttrue_endcaps -> GetYaxis()-> SetTitle("Events");
809  histSclEoEtrue_barrel -> GetXaxis()-> SetTitle("E/E_{gen}") ;
810  histSclEoEtrue_barrel -> GetYaxis()-> SetTitle("Events") ;
811  histSclEoEtrue_endcaps -> GetXaxis()-> SetTitle("E/E_{gen}") ;
812  histSclEoEtrue_endcaps -> GetYaxis()-> SetTitle("Events") ;
813  histSclEoEtrueGolden_barrel -> GetXaxis()-> SetTitle("E/E_{gen}") ;
814  histSclEoEtrueGolden_barrel -> GetYaxis()-> SetTitle("Events") ;
815  histSclEoEtrueShowering_barrel -> GetXaxis()-> SetTitle("E/E_{gen}") ;
816  histSclEoEtrueShowering_barrel -> GetYaxis()-> SetTitle("Events") ;
817  histSclEoEtrueGolden_endcaps -> GetXaxis()-> SetTitle("E/E_{gen}") ;
818  histSclEoEtrueGolden_endcaps -> GetYaxis()-> SetTitle("Events") ;
819  histSclEoEtrueShowering_endcaps -> GetXaxis()-> SetTitle("E/E_{gen}") ;
820  histSclEoEtrueShowering_endcaps -> GetYaxis()-> SetTitle("Events") ;
821  histSclEoEtrue_barrel_etagap -> GetXaxis()-> SetTitle("E/E_{gen}");
822  histSclEoEtrue_barrel_etagap -> GetYaxis()-> SetTitle("Events");
823  histSclEoEtrue_barrel_phigap -> GetXaxis()-> SetTitle("E/E_{gen}");
824  histSclEoEtrue_barrel_phigap -> GetYaxis()-> SetTitle("Events");
825  histSclEoEtrue_ebeegap -> GetXaxis()-> SetTitle("E/E_{gen}");
826  histSclEoEtrue_ebeegap -> GetYaxis()-> SetTitle("Events");
827  histSclEoEtrue_endcaps_deegap -> GetXaxis()-> SetTitle("E/E_{gen}");
828  histSclEoEtrue_endcaps_deegap -> GetYaxis()-> SetTitle("Events");
829  histSclEoEtrue_endcaps_ringgap -> GetXaxis()-> SetTitle("E/E_{gen}");
830  histSclEoEtrue_endcaps_ringgap -> GetYaxis()-> SetTitle("Events");
831  histSclEoEtrue_barrel_new -> GetXaxis()-> SetTitle("E/E_{gen}") ;
832  histSclEoEtrue_barrel_new -> GetYaxis()-> SetTitle("Events") ;
833  histSclEoEtrue_endcaps_new -> GetXaxis()-> SetTitle("E/E_{gen}") ;
834  histSclEoEtrue_endcaps_new -> GetYaxis()-> SetTitle("Events") ;
835  histSclEoEtrue_barrel_etagap_new -> GetXaxis()-> SetTitle("E/E_{gen}");
836  histSclEoEtrue_barrel_etagap_new -> GetYaxis()-> SetTitle("Events");
837  histSclEoEtrue_barrel_phigap_new -> GetXaxis()-> SetTitle("E/E_{gen}");
838  histSclEoEtrue_barrel_phigap_new -> GetYaxis()-> SetTitle("Events");
839  histSclEoEtrue_ebeegap_new -> GetXaxis()-> SetTitle("E/E_{gen}");
840  histSclEoEtrue_ebeegap_new -> GetYaxis()-> SetTitle("Events");
841  histSclEoEtrue_endcaps_deegap_new -> GetXaxis()-> SetTitle("E/E_{gen}");
842  histSclEoEtrue_endcaps_deegap_new -> GetYaxis()-> SetTitle("Events");
843  histSclEoEtrue_endcaps_ringgap_new -> GetXaxis()-> SetTitle("E/E_{gen}");
844  histSclEoEtrue_endcaps_ringgap_new -> GetYaxis()-> SetTitle("Events");
845  histSclSigEtaEta_-> GetXaxis()-> SetTitle("#sigma_{#eta #eta}") ;
846  histSclSigEtaEta_-> GetYaxis()-> SetTitle("Events") ;
847  histSclSigEtaEta_barrel_-> GetXaxis()-> SetTitle("#sigma_{#eta #eta}") ;
848  histSclSigEtaEta_barrel_-> GetYaxis()-> SetTitle("Events") ;
849  histSclSigEtaEta_endcaps_-> GetXaxis()-> SetTitle("#sigma_{#eta #eta}") ;
850  histSclSigEtaEta_endcaps_-> GetYaxis()-> SetTitle("Events") ;
851  histSclSigIEtaIEta_-> GetXaxis()-> SetTitle("#sigma_{i#eta i#eta}") ;
852  histSclSigIEtaIEta_-> GetYaxis()-> SetTitle("Events") ;
853  histSclSigIEtaIEta_barrel_-> GetXaxis()-> SetTitle("#sigma_{i#eta i#eta}") ;
854  histSclSigIEtaIEta_barrel_-> GetYaxis()-> SetTitle("Events") ;
855  histSclSigIEtaIEta_endcaps_-> GetXaxis()-> SetTitle("#sigma_{i#eta i#eta}") ;
856  histSclSigIEtaIEta_endcaps_-> GetYaxis()-> SetTitle("Events") ;
857  histSclE1x5_-> GetXaxis()-> SetTitle("E1x5 (GeV)") ;
858  histSclE1x5_-> GetYaxis()-> SetTitle("Events") ;
859  histSclE1x5_barrel_-> GetXaxis()-> SetTitle("E1x5 (GeV)") ;
860  histSclE1x5_barrel_-> GetYaxis()-> SetTitle("Events") ;
861  histSclE1x5_endcaps_-> GetXaxis()-> SetTitle("E1x5 (GeV)") ;
862  histSclE1x5_endcaps_-> GetYaxis()-> SetTitle("Events") ;
863  histSclE2x5max_-> GetXaxis()-> SetTitle("E2x5 (GeV)") ;
864  histSclE2x5max_-> GetYaxis()-> SetTitle("Events") ;
865  histSclE2x5max_barrel_-> GetXaxis()-> SetTitle("E2x5 (GeV)") ;
866  histSclE2x5max_barrel_-> GetYaxis()-> SetTitle("Events") ;
867  histSclE2x5max_endcaps_-> GetXaxis()-> SetTitle("E2x5 (GeV)") ;
868  histSclE2x5max_endcaps_-> GetYaxis()-> SetTitle("Events") ;
869  histSclE5x5_-> GetXaxis()-> SetTitle("E5x5 (GeV)") ;
870  histSclE5x5_-> GetYaxis()-> SetTitle("Events") ;
871  histSclE5x5_barrel_-> GetXaxis()-> SetTitle("E5x5 (GeV)") ;
872  histSclE5x5_barrel_-> GetYaxis()-> SetTitle("Events") ;
873  histSclE5x5_endcaps_-> GetXaxis()-> SetTitle("E5x5 (GeV)") ;
874  histSclE5x5_endcaps_-> GetYaxis()-> SetTitle("Events") ;
875  histSclEoEtruePfVsEg->GetXaxis()->SetTitle("E/E_{gen} (e/g)") ;
876  histSclEoEtruePfVsEg->GetYaxis()->SetTitle("E/E_{gen} (pflow)") ;
877  h_ele_ChargeMnChargeTrue -> GetXaxis()-> SetTitle("q_{rec} - q_{gen}");
878  h_ele_ChargeMnChargeTrue -> GetYaxis()-> SetTitle("Events");
879  h_ele_EtaMnEtaTrue -> GetXaxis()-> SetTitle("#eta_{rec} - #eta_{gen}");
880  h_ele_EtaMnEtaTrue -> GetYaxis()-> SetTitle("Events");
881  h_ele_EtaMnEtaTrue_barrel -> GetXaxis()-> SetTitle("#eta_{rec} - #eta_{gen}");
882  h_ele_EtaMnEtaTrue_barrel -> GetYaxis()-> SetTitle("Events");
883  h_ele_EtaMnEtaTrue_endcaps -> GetXaxis()-> SetTitle("#eta_{rec} - #eta_{gen}");
884  h_ele_EtaMnEtaTrue_endcaps -> GetYaxis()-> SetTitle("Events");
885  h_ele_PhiMnPhiTrue -> GetXaxis()-> SetTitle("#phi_{rec} - #phi_{gen} (rad)");
886  h_ele_PhiMnPhiTrue -> GetYaxis()-> SetTitle("Events");
887  h_ele_PhiMnPhiTrue_barrel -> GetXaxis()-> SetTitle("#phi_{rec} - #phi_{gen} (rad)");
888  h_ele_PhiMnPhiTrue_barrel -> GetYaxis()-> SetTitle("Events");
889  h_ele_PhiMnPhiTrue_endcaps -> GetXaxis()-> SetTitle("#phi_{rec} - #phi_{gen} (rad)");
890  h_ele_PhiMnPhiTrue_endcaps -> GetYaxis()-> SetTitle("Events");
891  h_ele_PinMnPout -> GetXaxis()-> SetTitle("P_{vertex} - P_{out} (GeV/c)");
892  h_ele_PinMnPout -> GetYaxis()-> SetTitle("Events");
893  h_ele_PinMnPout_mode -> GetXaxis()-> SetTitle("P_{vertex} - P_{out}, mode of GSF components (GeV/c)");
894  h_ele_PinMnPout_mode -> GetYaxis()-> SetTitle("Events");
895  h_ele_outerP -> GetXaxis()-> SetTitle("P_{out} (GeV/c)");
896  h_ele_outerP -> GetYaxis()-> SetTitle("Events");
897  h_ele_outerP_mode -> GetXaxis()-> SetTitle("P_{out} (GeV/c)");
898  h_ele_outerP_mode -> GetYaxis()-> SetTitle("Events");
899  h_ele_outerPt -> GetXaxis()-> SetTitle("P_{T out} (GeV/c)");
900  h_ele_outerPt -> GetYaxis()-> SetTitle("Events");
901  h_ele_outerPt_mode -> GetXaxis()-> SetTitle("P_{T out} (GeV/c)");
902  h_ele_outerPt_mode -> GetYaxis()-> SetTitle("Events");
903  h_ele_EoP -> GetXaxis()-> SetTitle("E/P_{vertex}");
904  h_ele_EoP -> GetYaxis()-> SetTitle("Events");
905  h_ele_EseedOP -> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
906  h_ele_EseedOP -> GetYaxis()-> SetTitle("Events");
907  h_ele_EoPout -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
908  h_ele_EoPout -> GetYaxis()-> SetTitle("Events");
909  h_ele_EeleOPout -> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
910  h_ele_EeleOPout -> GetYaxis()-> SetTitle("Events");
911  h_ele_EoP_barrel -> GetXaxis()-> SetTitle("E/P_{vertex}");
912  h_ele_EoP_barrel -> GetYaxis()-> SetTitle("Events");
913  h_ele_EseedOP_barrel -> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
914  h_ele_EseedOP_barrel -> GetYaxis()-> SetTitle("Events");
915  h_ele_EoPout_barrel -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
916  h_ele_EoPout_barrel -> GetYaxis()-> SetTitle("Events");
917  h_ele_EeleOPout_barrel -> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
918  h_ele_EeleOPout_barrel -> GetYaxis()-> SetTitle("Events");
919  h_ele_EoP_endcaps -> GetXaxis()-> SetTitle("E/P_{vertex}");
920  h_ele_EoP_endcaps -> GetYaxis()-> SetTitle("Events");
921  h_ele_EseedOP_endcaps -> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
922  h_ele_EseedOP_endcaps -> GetYaxis()-> SetTitle("Events");
923  h_ele_EoPout_endcaps -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
924  h_ele_EoPout_endcaps -> GetYaxis()-> SetTitle("Events");
925  h_ele_EeleOPout_endcaps -> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
926  h_ele_EeleOPout_endcaps -> GetYaxis()-> SetTitle("Events");
927  h_ele_vertexX-> GetXaxis()-> SetTitle("x (cm)");
928  h_ele_vertexX-> GetYaxis()-> SetTitle("Events");
929  h_ele_vertexY-> GetXaxis()-> SetTitle("y (cm)");
930  h_ele_vertexY-> GetYaxis()-> SetTitle("Events");
931  h_ele_vertexZ-> GetXaxis()-> SetTitle("z (cm)");
932  h_ele_vertexZ-> GetYaxis()-> SetTitle("Events");
933  h_ele_vertexTIP-> GetXaxis()-> SetTitle("TIP (cm)");
934  h_ele_vertexTIP-> GetYaxis()-> SetTitle("Events");
935  h_ele_TIP_all-> GetXaxis()-> SetTitle("r_{T} (cm)");
936  h_ele_TIP_all-> GetYaxis()-> SetTitle("Events");
937  h_ele_vertexTIPVsEta-> GetYaxis()-> SetTitle("TIP (cm)");
938  h_ele_vertexTIPVsEta-> GetXaxis()-> SetTitle("#eta");
939  h_ele_vertexTIPVsPhi-> GetYaxis()-> SetTitle("TIP (cm)");
940  h_ele_vertexTIPVsPhi-> GetXaxis()-> SetTitle("#phi (rad)");
941  h_ele_vertexTIPVsPt-> GetYaxis()-> SetTitle("TIP (cm)");
942  h_ele_vertexTIPVsPt-> GetXaxis()-> SetTitle("p_{T} (GeV/c)");
943  h_ele_dEtaSc_propVtx-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
944  h_ele_dEtaSc_propVtx-> GetYaxis()-> SetTitle("Events");
945  h_ele_dEtaCl_propOut-> GetXaxis()-> SetTitle("#eta_{seedcl} - #eta_{tr}");
946  h_ele_dEtaCl_propOut-> GetYaxis()-> SetTitle("Events");
947  h_ele_dEtaEleCl_propOut-> GetXaxis()-> SetTitle("#eta_{elecl} - #eta_{tr}");
948  h_ele_dEtaEleCl_propOut-> GetYaxis()-> SetTitle("Events");
949  h_ele_dPhiSc_propVtx-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
950  h_ele_dPhiSc_propVtx-> GetYaxis()-> SetTitle("Events");
951  h_ele_dPhiCl_propOut-> GetXaxis()-> SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
952  h_ele_dPhiCl_propOut-> GetYaxis()-> SetTitle("Events");
953  h_ele_dPhiEleCl_propOut-> GetXaxis()-> SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
954  h_ele_dPhiEleCl_propOut-> GetYaxis()-> SetTitle("Events");
955  h_ele_dEtaSc_propVtx_barrel-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
956  h_ele_dEtaSc_propVtx_barrel-> GetYaxis()-> SetTitle("Events");
957  h_ele_dEtaCl_propOut_barrel-> GetXaxis()-> SetTitle("#eta_{seedcl} - #eta_{tr}");
958  h_ele_dEtaCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
959  h_ele_dEtaEleCl_propOut_barrel-> GetXaxis()-> SetTitle("#eta_{elecl} - #eta_{tr}");
960  h_ele_dEtaEleCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
961  h_ele_dPhiSc_propVtx_barrel-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
962  h_ele_dPhiSc_propVtx_barrel-> GetYaxis()-> SetTitle("Events");
963  h_ele_dPhiCl_propOut_barrel-> GetXaxis()-> SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
964  h_ele_dPhiCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
965  h_ele_dPhiEleCl_propOut_barrel-> GetXaxis()-> SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
966  h_ele_dPhiEleCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
967  h_ele_dEtaSc_propVtx_endcaps-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
968  h_ele_dEtaSc_propVtx_endcaps-> GetYaxis()-> SetTitle("Events");
969  h_ele_dEtaCl_propOut_endcaps-> GetXaxis()-> SetTitle("#eta_{seedcl} - #eta_{tr}");
970  h_ele_dEtaCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
971  h_ele_dEtaEleCl_propOut_endcaps-> GetXaxis()-> SetTitle("#eta_{elecl} - #eta_{tr}");
972  h_ele_dEtaEleCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
973  h_ele_dPhiSc_propVtx_endcaps-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
974  h_ele_dPhiSc_propVtx_endcaps-> GetYaxis()-> SetTitle("Events");
975  h_ele_dPhiCl_propOut_endcaps-> GetXaxis()-> SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
976  h_ele_dPhiCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
977  h_ele_dPhiEleCl_propOut_endcaps-> GetXaxis()-> SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
978  h_ele_dPhiEleCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
979  h_ele_HoE-> GetXaxis()-> SetTitle("H/E") ;
980  h_ele_HoE-> GetYaxis()-> SetTitle("Events") ;
981  h_ele_HoE_barrel-> GetXaxis()-> SetTitle("H/E") ;
982  h_ele_HoE_barrel-> GetYaxis()-> SetTitle("Events") ;
983  h_ele_HoE_endcaps-> GetXaxis()-> SetTitle("H/E") ;
984  h_ele_HoE_endcaps-> GetYaxis()-> SetTitle("Events") ;
985  h_ele_HoE_fiducial-> GetXaxis()-> SetTitle("H/E") ;
986  h_ele_HoE_fiducial-> GetYaxis()-> SetTitle("Events") ;
987  h_ele_fbrem-> GetXaxis()-> SetTitle("P_{in} - P_{out} / P_{in}");
988  h_ele_fbrem-> GetYaxis()-> SetTitle("Events");
989  h_ele_seed_dphi2_-> GetXaxis()-> SetTitle("#phi_{hit}-#phi_{pred} (rad)") ;
990  h_ele_seed_dphi2_-> GetYaxis()-> SetTitle("Events") ;
991  h_ele_seed_drz2_-> GetXaxis()-> SetTitle("r(z)_{hit}-r(z)_{pred} (cm)") ;
992  h_ele_seed_drz2_-> GetYaxis()-> SetTitle("Events") ;
993  h_ele_seed_subdet2_-> GetXaxis()-> SetTitle("2nd hit subdet Id") ;
994  h_ele_seed_subdet2_-> GetYaxis()-> SetTitle("Events") ;
995  h_ele_classes-> GetXaxis()-> SetTitle("class Id") ;
996  h_ele_classes-> GetYaxis()-> SetTitle("Events") ;
997  h_ele_EoverP_all-> GetXaxis()-> SetTitle("E/P_{vertex}");
998  h_ele_EoverP_all-> GetYaxis()-> SetTitle("Events");
999  h_ele_EseedOP_all-> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
1000  h_ele_EseedOP_all-> GetYaxis()-> SetTitle("Events");
1001  h_ele_EoPout_all -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
1002  h_ele_EoPout_all-> GetYaxis()-> SetTitle("Events");
1003  h_ele_EeleOPout_all-> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
1004  h_ele_EeleOPout_all-> GetYaxis()-> SetTitle("Events");
1005  h_ele_dEtaSc_propVtx_all-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
1006  h_ele_dEtaSc_propVtx_all-> GetYaxis()-> SetTitle("Events");
1007  h_ele_dPhiSc_propVtx_all-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
1008  h_ele_dPhiSc_propVtx_all-> GetYaxis()-> SetTitle("Events");
1009  h_ele_dEtaCl_propOut_all-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
1010  h_ele_dEtaCl_propOut_all-> GetYaxis()-> SetTitle("Events");
1011  h_ele_dPhiCl_propOut_all-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
1012  h_ele_dPhiCl_propOut_all-> GetYaxis()-> SetTitle("Events");
1013  h_ele_HoE_all-> GetXaxis()-> SetTitle("H/E") ;
1014  h_ele_HoE_all-> GetYaxis()-> SetTitle("Events");
1015  h_ele_mee_all-> GetXaxis()-> SetTitle("m_{ee} (GeV/c^{2})");
1016  h_ele_mee_all-> GetYaxis()-> SetTitle("Events");
1017  h_ele_mee_os-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1018  h_ele_mee_os-> GetYaxis()-> SetTitle("Events");
1019  h_ele_mee_os_ebeb-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1020  h_ele_mee_os_ebeb-> GetYaxis()-> SetTitle("Events");
1021  h_ele_mee_os_ebee-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1022  h_ele_mee_os_ebee-> GetYaxis()-> SetTitle("Events");
1023  h_ele_mee_os_eeee-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1024  h_ele_mee_os_eeee-> GetYaxis()-> SetTitle("Events");
1025  h_ele_mee_os_gg-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1026  h_ele_mee_os_gg-> GetYaxis()-> SetTitle("Events");
1027  h_ele_mee_os_gb-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1028  h_ele_mee_os_gb-> GetYaxis()-> SetTitle("Events");
1029  h_ele_mee_os_bb-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1030  h_ele_mee_os_bb-> GetYaxis()-> SetTitle("Events");
1031  h_ele_E2mnE1vsMee_all-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1032  h_ele_E2mnE1vsMee_all-> GetYaxis()-> SetTitle("E2 - E1 (GeV)");
1033  h_ele_E2mnE1vsMee_egeg_all-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
1034  h_ele_E2mnE1vsMee_egeg_all-> GetYaxis()-> SetTitle("E2 - E1 (GeV)");
1035  histNum_-> GetXaxis()-> SetTitle("N_{ele}");
1036  histNum_-> GetYaxis()-> SetTitle("Events");
1037  h_ele_fbremVsEta_mode-> GetXaxis()-> SetTitle("#eta");
1038  h_ele_fbremVsEta_mean-> GetXaxis()-> SetTitle("#eta");
1039 
1040 }
1041 
1042 void
1043 GsfElectronMCAnalyzer::endJob(){
1044 
1045  histfile_->cd();
1046 
1047  std::cout << "[GsfElectronMCAnalyzer] efficiency calculation " << std::endl;
1048  // efficiency vs eta
1049  TH1F *h_ele_etaEff = (TH1F*)h_ele_simEta_matched->Clone("h_ele_etaEff");
1050  h_ele_etaEff->Reset();
1051  h_ele_etaEff->Divide(h_ele_simEta_matched,h_simEta,1,1,"b");
1052  h_ele_etaEff->Print();
1053  h_ele_etaEff->GetXaxis()->SetTitle("#eta");
1054  h_ele_etaEff->GetYaxis()->SetTitle("Efficiency");
1055 
1056  // efficiency vs z
1057  TH1F *h_ele_zEff = (TH1F*)h_ele_simZ_matched->Clone("h_ele_zEff");
1058  h_ele_zEff->Reset();
1059  h_ele_zEff->Divide(h_ele_simZ_matched,h_simZ,1,1,"b");
1060  h_ele_zEff->Print();
1061  h_ele_zEff->GetXaxis()->SetTitle("z (cm)");
1062  h_ele_zEff->GetYaxis()->SetTitle("Efficiency");
1063 
1064  // efficiency vs |eta|
1065  TH1F *h_ele_absetaEff = (TH1F*)h_ele_simAbsEta_matched->Clone("h_ele_absetaEff");
1066  h_ele_absetaEff->Reset();
1067  h_ele_absetaEff->Divide(h_ele_simAbsEta_matched,h_simAbsEta,1,1,"b");
1068  h_ele_absetaEff->GetXaxis()->SetTitle("|#eta|");
1069  h_ele_absetaEff->GetYaxis()->SetTitle("Efficiency");
1070 
1071  // efficiency vs pt
1072  TH1F *h_ele_ptEff = (TH1F*)h_ele_simPt_matched->Clone("h_ele_ptEff");
1073  h_ele_ptEff->Reset();
1074  h_ele_ptEff->Divide(h_ele_simPt_matched,h_simPt,1,1,"b");
1075  h_ele_ptEff->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1076  h_ele_ptEff->GetYaxis()->SetTitle("Efficiency");
1077 
1078  // efficiency vs phi
1079  TH1F *h_ele_phiEff = (TH1F*)h_ele_simPhi_matched->Clone("h_ele_phiEff");
1080  h_ele_phiEff->Reset();
1081  h_ele_phiEff->Divide(h_ele_simPhi_matched,h_simPhi,1,1,"b");
1082  h_ele_phiEff->GetXaxis()->SetTitle("#phi (rad)");
1083  h_ele_phiEff->GetYaxis()->SetTitle("Efficiency");
1084 
1085  // efficiency vs pt eta
1086  TH2F *h_ele_ptEtaEff = (TH2F*)h_ele_simPtEta_matched->Clone("h_ele_ptEtaEff");
1087  h_ele_ptEtaEff->Reset();
1088  h_ele_ptEtaEff->Divide(h_ele_simPtEta_matched,h_simPtEta,1,1,"b");
1089  h_ele_ptEtaEff->GetYaxis()->SetTitle("p_{T} (GeV/c)");
1090  h_ele_ptEtaEff->GetXaxis()->SetTitle("#eta");
1091 
1092  std::cout << "[GsfElectronMCAnalyzer] q-misid calculation " << std::endl;
1093  // misid vs eta
1094  TH1F *h_ele_etaQmisid = (TH1F*)h_ele_simEta_matched_qmisid->Clone("h_ele_etaQmisid");
1095  h_ele_etaQmisid->Reset();
1096  h_ele_etaQmisid->Divide(h_ele_simEta_matched_qmisid,h_simEta,1,1,"b");
1097  h_ele_etaQmisid->Print();
1098  h_ele_etaQmisid->GetXaxis()->SetTitle("#eta");
1099  h_ele_etaQmisid->GetYaxis()->SetTitle("q misId");
1100 
1101  // misid vs z
1102  TH1F *h_ele_zQmisid = (TH1F*)h_ele_simZ_matched_qmisid->Clone("h_ele_zQmisid");
1103  h_ele_zQmisid->Reset();
1104  h_ele_zQmisid->Divide(h_ele_simZ_matched_qmisid,h_simZ,1,1,"b");
1105  h_ele_zQmisid->Print();
1106  h_ele_zQmisid->GetXaxis()->SetTitle("z (cm)");
1107  h_ele_zQmisid->GetYaxis()->SetTitle("q misId");
1108 
1109  // misid vs |eta|
1110  TH1F *h_ele_absetaQmisid = (TH1F*)h_ele_simAbsEta_matched_qmisid->Clone("h_ele_absetaQmisid");
1111  h_ele_absetaQmisid->Reset();
1112  h_ele_absetaQmisid->Divide(h_ele_simAbsEta_matched_qmisid,h_simAbsEta,1,1,"b");
1113  h_ele_absetaQmisid->GetXaxis()->SetTitle("|#eta|");
1114  h_ele_absetaQmisid->GetYaxis()->SetTitle("q misId");
1115 
1116  // misid vs pt
1117  TH1F *h_ele_ptQmisid = (TH1F*)h_ele_simPt_matched_qmisid->Clone("h_ele_ptQmisid");
1118  h_ele_ptQmisid->Reset();
1119  h_ele_ptQmisid->Divide(h_ele_simPt_matched_qmisid,h_simPt,1,1,"b");
1120  h_ele_ptQmisid->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1121  h_ele_ptQmisid->GetYaxis()->SetTitle("q misId");
1122 
1123  std::cout << "[GsfElectronMCAnalyzer] all reco electrons " << std::endl;
1124  // rec/gen all electrons
1125  TH1F *h_ele_etaEff_all = (TH1F*)h_ele_vertexEta_all->Clone("h_ele_etaEff_all");
1126  h_ele_etaEff_all->Reset();
1127  h_ele_etaEff_all->Divide(h_ele_vertexEta_all,h_simEta,1,1,"b");
1128  h_ele_etaEff_all->Print();
1129  h_ele_etaEff_all->GetXaxis()->SetTitle("#eta");
1130  h_ele_etaEff_all->GetYaxis()->SetTitle("N_{rec}/N_{gen}");
1131  TH1F *h_ele_ptEff_all = (TH1F*)h_ele_vertexPt_all->Clone("h_ele_ptEff_all");
1132  h_ele_ptEff_all->Reset();
1133  h_ele_ptEff_all->Divide(h_ele_vertexPt_all,h_simPt,1,1,"b");
1134  h_ele_ptEff_all->Print();
1135  h_ele_ptEff_all->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1136  h_ele_ptEff_all->GetYaxis()->SetTitle("N_{rec}/N_{gen}");
1137 
1138  // classes
1139  TH1F *h_ele_eta_goldenFrac = (TH1F*)h_ele_eta_golden->Clone("h_ele_eta_goldenFrac");
1140  h_ele_eta_goldenFrac->Reset();
1141  h_ele_eta_goldenFrac->Divide(h_ele_eta_golden,h_ele_eta,1,1,"b");
1142  h_ele_eta_goldenFrac->GetXaxis()->SetTitle("|#eta|");
1143  h_ele_eta_goldenFrac->GetYaxis()->SetTitle("Fraction of electrons");
1144  h_ele_eta_goldenFrac->SetTitle("fraction of golden electrons vs eta");
1145  TH1F *h_ele_eta_bbremFrac = (TH1F*)h_ele_eta_bbrem->Clone("h_ele_eta_bbremFrac");
1146  h_ele_eta_bbremFrac->Reset();
1147  h_ele_eta_bbremFrac->GetXaxis()->SetTitle("|#eta|");
1148  h_ele_eta_bbremFrac->GetYaxis()->SetTitle("Fraction of electrons");
1149  h_ele_eta_bbremFrac->Divide(h_ele_eta_bbrem,h_ele_eta,1,1,"b");
1150  h_ele_eta_bbremFrac->SetTitle("fraction of big brem electrons vs eta");
1151  TH1F *h_ele_eta_narrowFrac = (TH1F*)h_ele_eta_narrow->Clone("h_ele_eta_narrowFrac");
1152  h_ele_eta_narrowFrac->Reset();
1153  h_ele_eta_narrowFrac->Divide(h_ele_eta_narrow,h_ele_eta,1,1,"b");
1154  h_ele_eta_narrowFrac->GetXaxis()->SetTitle("|#eta|");
1155  h_ele_eta_narrowFrac->GetYaxis()->SetTitle("Fraction of electrons");
1156  h_ele_eta_narrowFrac->SetTitle("fraction of narrow electrons vs eta");
1157  TH1F *h_ele_eta_showerFrac = (TH1F*)h_ele_eta_shower->Clone("h_ele_eta_showerFrac");
1158  h_ele_eta_showerFrac->Reset();
1159  h_ele_eta_showerFrac->Divide(h_ele_eta_shower,h_ele_eta,1,1,"b");
1160  h_ele_eta_showerFrac->GetXaxis()->SetTitle("|#eta|");
1161  h_ele_eta_showerFrac->GetYaxis()->SetTitle("Fraction of electrons");
1162  h_ele_eta_showerFrac->SetTitle("fraction of showering electrons vs eta");
1163 
1164  // fbrem
1165  TH1F *h_ele_xOverX0VsEta = new TH1F( "h_ele_xOverx0VsEta","mean X/X_0 vs eta",nbineta/2,0.0,2.5);
1166  for (int ibin=1;ibin<h_ele_fbremVsEta_mean->GetNbinsX()+1;ibin++) {
1167  double xOverX0 = 0.;
1168  if (h_ele_fbremVsEta_mean->GetBinContent(ibin)>0.) xOverX0 = -log(h_ele_fbremVsEta_mean->GetBinContent(ibin));
1169  h_ele_xOverX0VsEta->SetBinContent(ibin,xOverX0);
1170  }
1171 
1172  //profiles from 2D histos
1173  TProfile *p_ele_PoPtrueVsEta = h_ele_PoPtrueVsEta->ProfileX();
1174  p_ele_PoPtrueVsEta->SetTitle("mean ele momentum / gen momentum vs eta");
1175  p_ele_PoPtrueVsEta->GetXaxis()->SetTitle("#eta");
1176  p_ele_PoPtrueVsEta->GetYaxis()->SetTitle("<P/P_{gen}>");
1177  p_ele_PoPtrueVsEta->Write();
1178  TProfile *p_ele_PoPtrueVsPhi = h_ele_PoPtrueVsPhi->ProfileX();
1179  p_ele_PoPtrueVsPhi->SetTitle("mean ele momentum / gen momentum vs phi");
1180  p_ele_PoPtrueVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1181  p_ele_PoPtrueVsPhi->GetYaxis()->SetTitle("<P/P_{gen}>");
1182  p_ele_PoPtrueVsPhi->Write();
1183  TProfile *p_ele_EoEtruePfVsEg_x = histSclEoEtruePfVsEg->ProfileX();
1184  p_ele_EoEtruePfVsEg_x->SetTitle("mean pflow sc energy / true energy vs e/g sc energy");
1185  p_ele_EoEtruePfVsEg_x->GetXaxis()->SetTitle("E/E_{gen} (e/g)") ;
1186  p_ele_EoEtruePfVsEg_x->GetYaxis()->SetTitle("<E/E_{gen}> (pflow)") ;
1187  p_ele_EoEtruePfVsEg_x->Write();
1188  TProfile *p_ele_EoEtruePfVsEg_y = histSclEoEtruePfVsEg->ProfileY();
1189  p_ele_EoEtruePfVsEg_y->SetTitle("mean e/g sc energy / true energy vs pflow sc energy");
1190  p_ele_EoEtruePfVsEg_y->GetXaxis()->SetTitle("E/E_{gen} (pflow)") ;
1191  p_ele_EoEtruePfVsEg_y->GetYaxis()->SetTitle("<E/E_{gen}> (eg)") ;
1192  p_ele_EoEtruePfVsEg_y->Write();
1193  TProfile *p_ele_EtaMnEtaTrueVsEta = h_ele_EtaMnEtaTrueVsEta->ProfileX();
1194  p_ele_EtaMnEtaTrueVsEta->SetTitle("mean ele eta - gen eta vs eta");
1195  p_ele_EtaMnEtaTrueVsEta->GetXaxis()->SetTitle("#eta");
1196  p_ele_EtaMnEtaTrueVsEta->GetYaxis()->SetTitle("<#eta_{rec} - #eta_{gen}>");
1197  p_ele_EtaMnEtaTrueVsEta->Write();
1198  TProfile *p_ele_EtaMnEtaTrueVsPhi = h_ele_EtaMnEtaTrueVsPhi->ProfileX();
1199  p_ele_EtaMnEtaTrueVsPhi->SetTitle("mean ele eta - gen eta vs phi");
1200  p_ele_EtaMnEtaTrueVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1201  p_ele_EtaMnEtaTrueVsPhi->GetYaxis()->SetTitle("<#eta_{rec} - #eta_{gen}>");
1202  p_ele_EtaMnEtaTrueVsPhi->Write();
1203  TProfile *p_ele_PhiMnPhiTrueVsEta = h_ele_PhiMnPhiTrueVsEta->ProfileX();
1204  p_ele_PhiMnPhiTrueVsEta->SetTitle("mean ele phi - gen phi vs eta");
1205  p_ele_PhiMnPhiTrueVsEta->GetXaxis()->SetTitle("#eta");
1206  p_ele_PhiMnPhiTrueVsEta->GetYaxis()->SetTitle("<#phi_{rec} - #phi_{gen}> (rad)");
1207  p_ele_PhiMnPhiTrueVsEta->Write();
1208  TProfile *p_ele_PhiMnPhiTrueVsPhi = h_ele_PhiMnPhiTrueVsPhi->ProfileX();
1209  p_ele_PhiMnPhiTrueVsPhi->SetTitle("mean ele phi - gen phi vs phi");
1210  p_ele_PhiMnPhiTrueVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1211  p_ele_PhiMnPhiTrueVsPhi->Write();
1212  TProfile *p_ele_vertexPtVsEta = h_ele_vertexPtVsEta->ProfileX();
1213  p_ele_vertexPtVsEta->SetTitle("mean ele transverse momentum vs eta");
1214  p_ele_vertexPtVsEta->GetXaxis()->SetTitle("#eta");
1215  p_ele_vertexPtVsEta->GetYaxis()->SetTitle("<p_{T}> (GeV/c)");
1216  p_ele_vertexPtVsEta->Write();
1217  TProfile *p_ele_vertexPtVsPhi = h_ele_vertexPtVsPhi->ProfileX();
1218  p_ele_vertexPtVsPhi->SetTitle("mean ele transverse momentum vs phi");
1219  p_ele_vertexPtVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1220  p_ele_vertexPtVsPhi->GetYaxis()->SetTitle("<p_{T}> (GeV/c)");
1221  p_ele_vertexPtVsPhi->Write();
1222  TProfile *p_ele_EoPVsEta = h_ele_EoPVsEta->ProfileX();
1223  p_ele_EoPVsEta->SetTitle("mean ele E/p vs eta");
1224  p_ele_EoPVsEta->GetXaxis()->SetTitle("#eta");
1225  p_ele_EoPVsEta->GetYaxis()->SetTitle("<E/P_{vertex}>");
1226  p_ele_EoPVsEta->Write();
1227  TProfile *p_ele_EoPVsPhi = h_ele_EoPVsPhi->ProfileX();
1228  p_ele_EoPVsPhi->SetTitle("mean ele E/p vs phi");
1229  p_ele_EoPVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1230  p_ele_EoPVsPhi->GetYaxis()->SetTitle("<E/P_{vertex}>");
1231  p_ele_EoPVsPhi->Write();
1232  TProfile *p_ele_EoPoutVsEta = h_ele_EoPoutVsEta->ProfileX();
1233  p_ele_EoPoutVsEta->SetTitle("mean ele E/pout vs eta");
1234  p_ele_EoPoutVsEta->GetXaxis()->SetTitle("#eta");
1235  p_ele_EoPoutVsEta->GetYaxis()->SetTitle("<E_{seed}/P_{out}>");
1236  p_ele_EoPoutVsEta->Write();
1237  TProfile *p_ele_EoPoutVsPhi = h_ele_EoPoutVsPhi->ProfileX();
1238  p_ele_EoPoutVsPhi->SetTitle("mean ele E/pout vs phi");
1239  p_ele_EoPoutVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1240  p_ele_EoPoutVsPhi->GetYaxis()->SetTitle("<E_{seed}/P_{out}>");
1241  p_ele_EoPoutVsPhi->Write();
1242  TProfile *p_ele_EeleOPoutVsEta = h_ele_EeleOPoutVsEta->ProfileX();
1243  p_ele_EeleOPoutVsEta->SetTitle("mean ele Eele/pout vs eta");
1244  p_ele_EeleOPoutVsEta->GetXaxis()->SetTitle("#eta");
1245  p_ele_EeleOPoutVsEta->GetYaxis()->SetTitle("<E_{ele}/P_{out}>");
1246  p_ele_EeleOPoutVsEta->Write();
1247  TProfile *p_ele_EeleOPoutVsPhi = h_ele_EeleOPoutVsPhi->ProfileX();
1248  p_ele_EeleOPoutVsPhi->SetTitle("mean ele Eele/pout vs phi");
1249  p_ele_EeleOPoutVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1250  p_ele_EeleOPoutVsPhi->GetYaxis()->SetTitle("<E_{ele}/P_{out}>");
1251  p_ele_EeleOPoutVsPhi->Write();
1252  TProfile *p_ele_HoEVsEta = h_ele_HoEVsEta->ProfileX();
1253  p_ele_HoEVsEta->SetTitle("mean ele H/E vs eta");
1254  p_ele_HoEVsEta->GetXaxis()->SetTitle("#eta");
1255  p_ele_HoEVsEta->GetYaxis()->SetTitle("<H/E>");
1256  p_ele_HoEVsEta->Write();
1257  TProfile *p_ele_HoEVsPhi = h_ele_HoEVsPhi->ProfileX();
1258  p_ele_HoEVsPhi->SetTitle("mean ele H/E vs phi");
1259  p_ele_HoEVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1260  p_ele_HoEVsPhi->GetYaxis()->SetTitle("<H/E>");
1261  p_ele_HoEVsPhi->Write();
1262  TProfile *p_ele_chi2VsEta = h_ele_chi2VsEta->ProfileX();
1263  p_ele_chi2VsEta->SetTitle("mean ele track chi2 vs eta");
1264  p_ele_chi2VsEta->GetXaxis()->SetTitle("#eta");
1265  p_ele_chi2VsEta->GetYaxis()->SetTitle("<#Chi^{2}>");
1266  p_ele_chi2VsEta->Write();
1267  TProfile *p_ele_chi2VsPhi = h_ele_chi2VsPhi->ProfileX();
1268  p_ele_chi2VsPhi->SetTitle("mean ele track chi2 vs phi");
1269  p_ele_chi2VsPhi->GetXaxis()->SetTitle("#phi (rad)");
1270  p_ele_chi2VsPhi->GetYaxis()->SetTitle("<#Chi^{2}>");
1271  p_ele_chi2VsPhi->Write();
1272  TProfile *p_ele_foundHitsVsEta = h_ele_foundHitsVsEta->ProfileX();
1273  p_ele_foundHitsVsEta->SetTitle("mean ele track # found hits vs eta");
1274  p_ele_foundHitsVsEta->GetXaxis()->SetTitle("#eta");
1275  p_ele_foundHitsVsEta->GetYaxis()->SetTitle("<N_{hits}>");
1276  p_ele_foundHitsVsEta->Write();
1277  TProfile *p_ele_foundHitsVsPhi = h_ele_foundHitsVsPhi->ProfileX();
1278  p_ele_foundHitsVsPhi->SetTitle("mean ele track # found hits vs phi");
1279  p_ele_foundHitsVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1280  p_ele_foundHitsVsPhi->GetYaxis()->SetTitle("<N_{hits}>");
1281  p_ele_foundHitsVsPhi->Write();
1282  TProfile *p_ele_lostHitsVsEta = h_ele_lostHitsVsEta->ProfileX();
1283  p_ele_lostHitsVsEta->SetTitle("mean ele track # lost hits vs eta");
1284  p_ele_lostHitsVsEta->GetXaxis()->SetTitle("#eta");
1285  p_ele_lostHitsVsEta->GetYaxis()->SetTitle("<N_{hits}>");
1286  p_ele_lostHitsVsEta->Write();
1287  TProfile *p_ele_lostHitsVsPhi = h_ele_lostHitsVsPhi->ProfileX();
1288  p_ele_lostHitsVsPhi->SetTitle("mean ele track # lost hits vs phi");
1289  p_ele_lostHitsVsPhi->GetXaxis()->SetTitle("#phi (rad)");
1290  p_ele_lostHitsVsPhi->GetYaxis()->SetTitle("<N_{hits}>");
1291  p_ele_lostHitsVsPhi->Write();
1292  TProfile *p_ele_vertexTIPVsEta = h_ele_vertexTIPVsEta->ProfileX();
1293  p_ele_vertexTIPVsEta->SetTitle("mean tip (wrt gen vtx) vs eta");
1294  p_ele_vertexTIPVsEta->GetXaxis()->SetTitle("#eta");
1295  p_ele_vertexTIPVsEta->GetYaxis()->SetTitle("<TIP> (cm)");
1296  p_ele_vertexTIPVsEta->Write();
1297  TProfile *p_ele_vertexTIPVsPhi = h_ele_vertexTIPVsPhi->ProfileX();
1298  p_ele_vertexTIPVsPhi->SetTitle("mean tip (wrt gen vtx) vs phi");
1299  p_ele_vertexTIPVsPhi->GetXaxis()->SetTitle("#phi");
1300  p_ele_vertexTIPVsPhi->GetYaxis()->SetTitle("<TIP> (cm)");
1301  p_ele_vertexTIPVsPhi->Write();
1302  TProfile *p_ele_vertexTIPVsPt = h_ele_vertexTIPVsPt->ProfileX();
1303  p_ele_vertexTIPVsPt->SetTitle("mean tip (wrt gen vtx) vs phi");
1304  p_ele_vertexTIPVsPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1305  p_ele_vertexTIPVsPt->GetYaxis()->SetTitle("<TIP> (cm)");
1306  p_ele_vertexTIPVsPt->Write();
1307 
1308  // mc truth
1309  h_mcNum->Write();
1310  h_eleNum->Write();
1311  h_gamNum->Write();
1312 
1313  // rec event
1314  histNum_->Write();
1315 
1316  // mc
1317  h_simEta->Write();
1318  h_simAbsEta->Write();
1319  h_simP->Write();
1320  h_simPt->Write();
1321  h_simZ->Write();
1322  h_simPhi->Write();
1323  h_simPtEta->Write();
1324 
1325  // all electrons
1326  h_ele_EoverP_all->Write();
1327  h_ele_EseedOP_all->Write();
1328  h_ele_EoPout_all->Write();
1329  h_ele_EeleOPout_all ->Write();
1330  h_ele_dEtaSc_propVtx_all->Write();
1331  h_ele_dPhiSc_propVtx_all->Write();
1332  h_ele_dEtaCl_propOut_all ->Write();
1333  h_ele_dPhiCl_propOut_all->Write();
1334  h_ele_HoE_all->Write();
1335  h_ele_TIP_all->Write();
1336  h_ele_vertexPt_all->Write();
1337  h_ele_Et_all->Write();
1338  h_ele_vertexEta_all->Write();
1339  h_ele_mee_all->Write();
1340  h_ele_mee_os->Write();
1341  h_ele_mee_os_ebeb->Write();
1342  h_ele_mee_os_ebee->Write();
1343  h_ele_mee_os_eeee->Write();
1344  h_ele_mee_os_gg->Write();
1345  h_ele_mee_os_gb->Write();
1346  h_ele_mee_os_bb->Write();
1347  h_ele_E2mnE1vsMee_all ->Write();
1348  h_ele_E2mnE1vsMee_egeg_all->Write();
1349 
1350  // charge ID
1351  h_ele_charge->Write();
1352  h_ele_simEta_matched_qmisid->Write();
1353  h_ele_simAbsEta_matched_qmisid->Write();
1354  h_ele_simPt_matched_qmisid->Write();
1355  h_ele_simPhi_matched_qmisid->Write();
1356  h_ele_simZ_matched_qmisid->Write();
1357 
1358  // matched electrons
1359  h_ele_vertexP->Write();
1360  h_ele_vertexPt->Write();
1361  h_ele_Et->Write();
1362  h_ele_vertexPtVsEta->Write();
1363  h_ele_vertexPtVsPhi->Write();
1364  h_ele_simPt_matched->Write();
1365  h_ele_vertexEta->Write();
1366  h_ele_vertexEtaVsPhi->Write();
1367  h_ele_simAbsEta_matched->Write();
1368  h_ele_simEta_matched->Write();
1369  h_ele_simPhi_matched->Write();
1370  h_ele_simPtEta_matched->Write();
1371  h_ele_vertexPhi->Write();
1372  h_ele_vertexX->Write();
1373  h_ele_vertexY ->Write();
1374  h_ele_vertexZ->Write();
1375  h_ele_vertexTIP->Write();
1376  h_ele_simZ_matched->Write();
1377  h_ele_vertexTIPVsEta->Write();
1378  h_ele_vertexTIPVsPhi->Write();
1379  h_ele_vertexTIPVsPt->Write();
1380  h_ele_PoPtrue->Write();
1381  h_ele_PoPtrueVsEta ->Write();
1382  h_ele_PoPtrueVsPhi->Write();
1383  h_ele_PoPtrueVsPt->Write();
1384  h_ele_PoPtrue_barrel ->Write();
1385  h_ele_PoPtrue_endcaps->Write();
1386  h_ele_PoPtrue_golden_barrel ->Write();
1387  h_ele_PoPtrue_golden_endcaps->Write();
1388  h_ele_PoPtrue_showering_barrel ->Write();
1389  h_ele_PoPtrue_showering_endcaps->Write();
1390  h_ele_PtoPttrue->Write();
1391  h_ele_PtoPttrue_barrel ->Write();
1392  h_ele_PtoPttrue_endcaps->Write();
1393  h_ele_ChargeMnChargeTrue->Write();
1394  h_ele_EtaMnEtaTrue->Write();
1395  h_ele_EtaMnEtaTrue_barrel->Write();
1396  h_ele_EtaMnEtaTrue_endcaps->Write();
1397  h_ele_EtaMnEtaTrueVsEta ->Write();
1398  h_ele_EtaMnEtaTrueVsPhi->Write();
1399  h_ele_EtaMnEtaTrueVsPt->Write();
1400  h_ele_PhiMnPhiTrue ->Write();
1401  h_ele_PhiMnPhiTrue_barrel ->Write();
1402  h_ele_PhiMnPhiTrue_endcaps ->Write();
1403  h_ele_PhiMnPhiTrue2 ->Write();
1404  h_ele_PhiMnPhiTrueVsEta->Write();
1405  h_ele_PhiMnPhiTrueVsPhi->Write();
1406  h_ele_PhiMnPhiTrueVsPt->Write();
1407 
1408  // matched electron, superclusters
1409  histSclEn_->Write();
1410  histSclEoEtrue_barrel->Write();
1411  histSclEoEtrue_endcaps->Write();
1412  histSclEoEtrue_barrel_eg->Write();
1413  histSclEoEtrue_endcaps_eg->Write();
1414  histSclEoEtrue_barrel_etagap->Write();
1415  histSclEoEtrue_barrel_phigap->Write();
1416  histSclEoEtrue_ebeegap->Write();
1417  histSclEoEtrue_endcaps->Write();
1418  histSclEoEtrue_endcaps_deegap->Write();
1419  histSclEoEtrue_endcaps_ringgap->Write();
1420  histSclEoEtruePfVsEg->Write();
1421  histSclEoEtrue_barrel_new->Write();
1422  histSclEoEtrue_endcaps_new->Write();
1423  histSclEoEtrue_barrel_eg_new->Write();
1424  histSclEoEtrue_endcaps_eg_new->Write();
1425  histSclEoEtrue_barrel_etagap_new->Write();
1426  histSclEoEtrue_barrel_phigap_new->Write();
1427  histSclEoEtrue_ebeegap_new->Write();
1428  histSclEoEtrue_endcaps_new->Write();
1429  histSclEoEtrue_endcaps_deegap_new->Write();
1430  histSclEoEtrue_endcaps_ringgap_new->Write();
1431  histSclEoEtruePfVsEg->Write();
1432  histSclEt_->Write();
1433  histSclEtVsEta_->Write();
1434  histSclEtVsPhi_->Write();
1435  histSclEtaVsPhi_ ->Write();
1436  histSclEta_->Write();
1437  histSclPhi_->Write();
1438  histSclSigEtaEta_->Write();
1439  histSclSigEtaEta_barrel_->Write();
1440  histSclSigEtaEta_endcaps_->Write();
1441  histSclSigIEtaIEta_->Write();
1442  histSclSigIEtaIEta_barrel_->Write();
1443  histSclSigIEtaIEta_endcaps_->Write();
1444  histSclE1x5_->Write();
1445  histSclE1x5_barrel_->Write();
1446  histSclE1x5_endcaps_->Write();
1447  histSclE2x5max_->Write();
1448  histSclE2x5max_barrel_->Write();
1449  histSclE2x5max_endcaps_->Write();
1450  histSclE5x5_->Write();
1451  histSclE5x5_barrel_->Write();
1452  histSclE5x5_endcaps_->Write();
1453  histSclSigEtaEta_eg_->Write();
1454  histSclSigEtaEta_eg_barrel_->Write();
1455  histSclSigEtaEta_eg_endcaps_->Write();
1456  histSclSigIEtaIEta_eg_->Write();
1457  histSclSigIEtaIEta_eg_barrel_->Write();
1458  histSclSigIEtaIEta_eg_endcaps_->Write();
1459  histSclE1x5_eg_->Write();
1460  histSclE1x5_eg_barrel_->Write();
1461  histSclE1x5_eg_endcaps_->Write();
1462  histSclE2x5max_eg_->Write();
1463  histSclE2x5max_eg_barrel_->Write();
1464  histSclE2x5max_eg_endcaps_->Write();
1465  histSclE5x5_eg_->Write();
1466  histSclE5x5_eg_barrel_->Write();
1467  histSclE5x5_eg_endcaps_->Write();
1468 
1469  // matched electron, gsf tracks
1470  h_ele_ambiguousTracks->Write();
1471  h_ele_ambiguousTracksVsEta->Write();
1472  h_ele_ambiguousTracksVsPhi->Write();
1473  h_ele_ambiguousTracksVsPt->Write();
1474  h_ele_foundHits->Write();
1475  h_ele_foundHits_barrel->Write();
1476  h_ele_foundHits_endcaps->Write();
1477  h_ele_foundHitsVsEta->Write();
1478  h_ele_foundHitsVsPhi->Write();
1479  h_ele_foundHitsVsPt->Write();
1480  h_ele_lostHits->Write();
1481  h_ele_lostHits_barrel->Write();
1482  h_ele_lostHits_endcaps->Write();
1483  h_ele_lostHitsVsEta->Write();
1484  h_ele_lostHitsVsPhi->Write();
1485  h_ele_lostHitsVsPt->Write();
1486  h_ele_chi2 ->Write();
1487  h_ele_chi2_barrel ->Write();
1488  h_ele_chi2_endcaps ->Write();
1489  h_ele_chi2VsEta ->Write();
1490  h_ele_chi2VsPhi ->Write();
1491  h_ele_chi2VsPt->Write();
1492  h_ele_PinMnPout->Write();
1493  h_ele_PinMnPout_mode->Write();
1494  h_ele_PinMnPoutVsEta_mode->Write();
1495  h_ele_PinMnPoutVsPhi_mode->Write();
1496  h_ele_PinMnPoutVsPt_mode->Write();
1497  h_ele_PinMnPoutVsE_mode->Write();
1498  h_ele_PinMnPoutVsChi2_mode->Write();
1499  h_ele_outerP ->Write();
1500  h_ele_outerP_mode->Write();
1501  h_ele_outerPVsEta_mode->Write();
1502  h_ele_outerPt->Write();
1503  h_ele_outerPt_mode ->Write();
1504  h_ele_outerPtVsEta_mode->Write();
1505  h_ele_outerPtVsPhi_mode->Write();
1506  h_ele_outerPtVsPt_mode->Write();
1507 
1508  // matched electrons, matching
1509  h_ele_EoP ->Write();
1510  h_ele_EoP_eg ->Write();
1511  h_ele_EoP_barrel ->Write();
1512  h_ele_EoP_eg_barrel ->Write();
1513  h_ele_EoP_endcaps ->Write();
1514  h_ele_EoP_eg_endcaps ->Write();
1515  h_ele_EoPVsEta ->Write();
1516  h_ele_EoPVsPhi->Write();
1517  h_ele_EoPVsE->Write();
1518  h_ele_EseedOP ->Write();
1519  h_ele_EseedOP_eg ->Write();
1520  h_ele_EseedOP_barrel ->Write();
1521  h_ele_EseedOP_eg_barrel ->Write();
1522  h_ele_EseedOP_endcaps ->Write();
1523  h_ele_EseedOP_eg_endcaps ->Write();
1524  h_ele_EseedOPVsEta ->Write();
1525  h_ele_EseedOPVsPhi->Write();
1526  h_ele_EseedOPVsE->Write();
1527  h_ele_EoPout->Write();
1528  h_ele_EoPout_eg->Write();
1529  h_ele_EoPout_barrel->Write();
1530  h_ele_EoPout_eg_barrel->Write();
1531  h_ele_EoPout_endcaps->Write();
1532  h_ele_EoPout_eg_endcaps->Write();
1533  h_ele_EoPoutVsEta->Write();
1534  h_ele_EoPoutVsPhi->Write();
1535  h_ele_EoPoutVsE ->Write();
1536  h_ele_EeleOPout->Write();
1537  h_ele_EeleOPout_eg->Write();
1538  h_ele_EeleOPout_barrel->Write();
1539  h_ele_EeleOPout_eg_barrel->Write();
1540  h_ele_EeleOPout_endcaps->Write();
1541  h_ele_EeleOPout_eg_endcaps->Write();
1542  h_ele_EeleOPoutVsEta->Write();
1543  h_ele_EeleOPoutVsPhi->Write();
1544  h_ele_EeleOPoutVsE ->Write();
1545  h_ele_dEtaSc_propVtx->Write();
1546  h_ele_dEtaSc_propVtx_eg->Write();
1547  h_ele_dEtaSc_propVtx_barrel->Write();
1548  h_ele_dEtaSc_propVtx_eg_barrel->Write();
1549  h_ele_dEtaSc_propVtx_endcaps->Write();
1550  h_ele_dEtaSc_propVtx_eg_endcaps->Write();
1551  h_ele_dEtaScVsEta_propVtx->Write();
1552  h_ele_dEtaScVsPhi_propVtx->Write();
1553  h_ele_dEtaScVsPt_propVtx ->Write();
1554  h_ele_dPhiSc_propVtx->Write();
1555  h_ele_dPhiSc_propVtx_eg->Write();
1556  h_ele_dPhiSc_propVtx_barrel->Write();
1557  h_ele_dPhiSc_propVtx_eg_barrel->Write();
1558  h_ele_dPhiSc_propVtx_endcaps->Write();
1559  h_ele_dPhiSc_propVtx_eg_endcaps->Write();
1560  h_ele_dPhiScVsEta_propVtx ->Write();
1561  h_ele_dPhiScVsPhi_propVtx->Write();
1562  h_ele_dPhiScVsPt_propVtx->Write();
1563  h_ele_dEtaCl_propOut->Write();
1564  h_ele_dEtaCl_propOut_eg->Write();
1565  h_ele_dEtaCl_propOut_barrel->Write();
1566  h_ele_dEtaCl_propOut_eg_barrel->Write();
1567  h_ele_dEtaCl_propOut_endcaps->Write();
1568  h_ele_dEtaCl_propOut_eg_endcaps->Write();
1569  h_ele_dEtaClVsEta_propOut->Write();
1570  h_ele_dEtaClVsPhi_propOut->Write();
1571  h_ele_dEtaClVsPt_propOut->Write();
1572  h_ele_dPhiCl_propOut->Write();
1573  h_ele_dPhiCl_propOut_eg->Write();
1574  h_ele_dPhiCl_propOut_barrel->Write();
1575  h_ele_dPhiCl_propOut_eg_barrel->Write();
1576  h_ele_dPhiCl_propOut_endcaps->Write();
1577  h_ele_dPhiCl_propOut_eg_endcaps->Write();
1578  h_ele_dPhiClVsEta_propOut->Write();
1579  h_ele_dPhiClVsPhi_propOut->Write();
1580  h_ele_dPhiClVsPt_propOut->Write();
1581  h_ele_dEtaEleCl_propOut->Write();
1582  h_ele_dEtaEleCl_propOut_eg->Write();
1583  h_ele_dEtaEleCl_propOut_barrel->Write();
1584  h_ele_dEtaEleCl_propOut_eg_barrel->Write();
1585  h_ele_dEtaEleCl_propOut_endcaps->Write();
1586  h_ele_dEtaEleCl_propOut_eg_endcaps->Write();
1587  h_ele_dEtaEleClVsEta_propOut->Write();
1588  h_ele_dEtaEleClVsPhi_propOut->Write();
1589  h_ele_dEtaEleClVsPt_propOut->Write();
1590  h_ele_dPhiEleCl_propOut->Write();
1591  h_ele_dPhiEleCl_propOut_eg->Write();
1592  h_ele_dPhiEleCl_propOut_barrel->Write();
1593  h_ele_dPhiEleCl_propOut_eg_barrel->Write();
1594  h_ele_dPhiEleCl_propOut_endcaps->Write();
1595  h_ele_dPhiEleCl_propOut_eg_endcaps->Write();
1596  h_ele_dPhiEleClVsEta_propOut->Write();
1597  h_ele_dPhiEleClVsPhi_propOut->Write();
1598  h_ele_dPhiEleClVsPt_propOut->Write();
1599  h_ele_HoE->Write();
1600  h_ele_HoE_eg->Write();
1601  h_ele_HoE_barrel->Write();
1602  h_ele_HoE_eg_barrel->Write();
1603  h_ele_HoE_endcaps->Write();
1604  h_ele_HoE_eg_endcaps->Write();
1605  h_ele_HoE_fiducial->Write();
1606  h_ele_HoEVsEta->Write();
1607  h_ele_HoEVsPhi->Write();
1608  h_ele_HoEVsE->Write();
1609 
1610  h_ele_seed_dphi2_->Write();
1611  h_ele_seed_subdet2_->Write();
1612  TProfile *p_ele_seed_dphi2VsEta_ = h_ele_seed_dphi2VsEta_->ProfileX();
1613  p_ele_seed_dphi2VsEta_->SetTitle("mean ele seed dphi 2nd layer vs eta");
1614  p_ele_seed_dphi2VsEta_->GetXaxis()->SetTitle("#eta");
1615  p_ele_seed_dphi2VsEta_->GetYaxis()->SetTitle("<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)");
1616  p_ele_seed_dphi2VsEta_->SetMinimum(-0.004);
1617  p_ele_seed_dphi2VsEta_->SetMaximum(0.004);
1618  p_ele_seed_dphi2VsEta_->Write();
1619  TProfile *p_ele_seed_dphi2VsPt_ = h_ele_seed_dphi2VsPt_->ProfileX();
1620  p_ele_seed_dphi2VsPt_->SetTitle("mean ele seed dphi 2nd layer vs pt");
1621  p_ele_seed_dphi2VsPt_->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1622  p_ele_seed_dphi2VsPt_->GetYaxis()->SetTitle("<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)");
1623  p_ele_seed_dphi2VsPt_->SetMinimum(-0.004);
1624  p_ele_seed_dphi2VsPt_->SetMaximum(0.004);
1625  p_ele_seed_dphi2VsPt_->Write();
1626  h_ele_seed_drz2_->Write();
1627  TProfile *p_ele_seed_drz2VsEta_ = h_ele_seed_drz2VsEta_->ProfileX();
1628  p_ele_seed_drz2VsEta_->SetTitle("mean ele seed dr(dz) 2nd layer vs eta");
1629  p_ele_seed_drz2VsEta_->GetXaxis()->SetTitle("#eta");
1630  p_ele_seed_drz2VsEta_->GetYaxis()->SetTitle("<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)");
1631  p_ele_seed_drz2VsEta_->SetMinimum(-0.15);
1632  p_ele_seed_drz2VsEta_->SetMaximum(0.15);
1633  p_ele_seed_drz2VsEta_->Write();
1634  TProfile *p_ele_seed_drz2VsPt_ = h_ele_seed_drz2VsPt_->ProfileX();
1635  p_ele_seed_drz2VsPt_->SetTitle("mean ele seed dr(dz) 2nd layer vs pt");
1636  p_ele_seed_drz2VsPt_->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1637  p_ele_seed_drz2VsPt_->GetYaxis()->SetTitle("<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)");
1638  p_ele_seed_drz2VsPt_->SetMinimum(-0.15);
1639  p_ele_seed_drz2VsPt_->SetMaximum(0.15);
1640  p_ele_seed_drz2VsPt_->Write();
1641 
1642  // classes
1643  h_ele_classes->Write();
1644  h_ele_eta->Write();
1645  h_ele_eta_golden->Write();
1646  h_ele_eta_bbrem->Write();
1647  h_ele_eta_narrow->Write();
1648  h_ele_eta_shower->Write();
1649  h_ele_PinVsPoutGolden_mode->Write();
1650  h_ele_PinVsPoutShowering_mode->Write();
1651  h_ele_PinVsPoutGolden_mean->Write();
1652  h_ele_PinVsPoutShowering_mean->Write();
1653  h_ele_PtinVsPtoutGolden_mode->Write();
1654  h_ele_PtinVsPtoutShowering_mode->Write();
1655  h_ele_PtinVsPtoutGolden_mean->Write();
1656  h_ele_PtinVsPtoutShowering_mean->Write();
1657  histSclEoEtrueGolden_barrel->Write();
1658  histSclEoEtrueGolden_endcaps->Write();
1659  histSclEoEtrueShowering_barrel->Write();
1660  histSclEoEtrueShowering_endcaps->Write();
1661 
1662  // fbrem
1663  h_ele_fbrem->Write();
1664  h_ele_fbrem_eg->Write();
1665  h_ele_fbremVsEta_mode->GetXaxis()->SetTitle("#eta");
1666  h_ele_fbremVsEta_mode->GetYaxis()->SetTitle("<P_{in} - P_{out} / P_{in}>");
1667  h_ele_fbremVsEta_mode->Write();
1668  h_ele_fbremVsEta_mean->GetXaxis()->SetTitle("#eta");
1669  h_ele_fbremVsEta_mean->GetYaxis()->SetTitle("<P_{in} - P_{out} / P_{in}>");
1670  h_ele_fbremVsEta_mean->Write();
1671  h_ele_eta_goldenFrac->Write();
1672  h_ele_eta_bbremFrac->Write();
1673  h_ele_eta_narrowFrac->Write();
1674  h_ele_eta_showerFrac->Write();
1675  h_ele_xOverX0VsEta->Write();
1676 
1677  // efficiencies
1678  h_ele_etaEff->Write();
1679  h_ele_zEff->Write();
1680  h_ele_phiEff->Write();
1681  h_ele_absetaEff->Write();
1682  h_ele_ptEff->Write();
1683  h_ele_ptEtaEff->Write();
1684  h_ele_etaEff_all->Write();
1685  h_ele_ptEff_all->Write();
1686 
1687  // q misid
1688  h_ele_etaQmisid->Write();
1689  h_ele_zQmisid->Write();
1690  h_ele_absetaQmisid->Write();
1691  h_ele_ptQmisid->Write();
1692 
1693  // e/g et pflow electrons
1694  h_ele_mva->Write();
1695  h_ele_mva_eg->Write();
1696  h_ele_provenance->Write();
1697 
1698  // isolation
1699  h_ele_tkSumPt_dr03->GetXaxis()->SetTitle("TkIsoSum, cone 0.3 (GeV/c)");
1700  h_ele_tkSumPt_dr03->GetYaxis()->SetTitle("Events");
1701  h_ele_tkSumPt_dr03->Write();
1702  h_ele_ecalRecHitSumEt_dr03->GetXaxis()->SetTitle("EcalIsoSum, cone 0.3 (GeV)");
1703  h_ele_ecalRecHitSumEt_dr03->GetYaxis()->SetTitle("Events");
1704  h_ele_ecalRecHitSumEt_dr03->Write();
1705  h_ele_hcalDepth1TowerSumEt_dr03->GetXaxis()->SetTitle("Hcal1IsoSum, cone 0.3 (GeV)");
1706  h_ele_hcalDepth1TowerSumEt_dr03->GetYaxis()->SetTitle("Events");
1707  h_ele_hcalDepth1TowerSumEt_dr03->Write();
1708  h_ele_hcalDepth2TowerSumEt_dr03->GetXaxis()->SetTitle("Hcal2IsoSum, cone 0.3 (GeV)");
1709  h_ele_hcalDepth2TowerSumEt_dr03->GetYaxis()->SetTitle("Events");
1710  h_ele_hcalDepth2TowerSumEt_dr03->Write();
1711  h_ele_tkSumPt_dr04->GetXaxis()->SetTitle("TkIsoSum, cone 0.4 (GeV/c)");
1712  h_ele_tkSumPt_dr04->GetYaxis()->SetTitle("Events");
1713  h_ele_tkSumPt_dr04->Write();
1714  h_ele_ecalRecHitSumEt_dr04->GetXaxis()->SetTitle("EcalIsoSum, cone 0.4 (GeV)");
1715  h_ele_ecalRecHitSumEt_dr04->GetYaxis()->SetTitle("Events");
1716  h_ele_ecalRecHitSumEt_dr04->Write();
1717  h_ele_hcalDepth1TowerSumEt_dr04->GetXaxis()->SetTitle("Hcal1IsoSum, cone 0.4 (GeV)");
1718  h_ele_hcalDepth1TowerSumEt_dr04->GetYaxis()->SetTitle("Events");
1719  h_ele_hcalDepth1TowerSumEt_dr04->Write();
1720  h_ele_hcalDepth2TowerSumEt_dr04->GetXaxis()->SetTitle("Hcal2IsoSum, cone 0.4 (GeV)");
1721  h_ele_hcalDepth2TowerSumEt_dr04->GetYaxis()->SetTitle("Events");
1722  h_ele_hcalDepth2TowerSumEt_dr04->Write();
1723 
1724 }
1725 
1726 GsfElectronMCAnalyzer::~GsfElectronMCAnalyzer()
1727 {
1728 
1729  // do anything here that needs to be done at desctruction time
1730  // (e.g. close files, deallocate resources etc.)
1731  histfile_->Write();
1732  histfile_->Close();
1733 
1734 }
1735 
1736 
1737 //=========================================================================
1738 // Main method
1739 //=========================================================================
1740 
1741 void
1742 GsfElectronMCAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
1743 {
1744  std::cout << "analyzing new event " << std::endl;
1745  // get electrons
1746 
1747  edm::Handle<GsfElectronCollection> gsfElectrons;
1748  iEvent.getByLabel(electronCollection_,gsfElectrons);
1749  edm::LogInfo("")<<"\n\n =================> Treating event "<<iEvent.id()<<" Number of electrons "<<gsfElectrons.product()->size();
1750 
1751  edm::Handle<GenParticleCollection> genParticles;
1752  iEvent.getByLabel(mcTruthCollection_, genParticles);
1753 
1754  histNum_->Fill((*gsfElectrons).size());
1755 
1756  // all rec electrons
1757  for (reco::GsfElectronCollection::const_iterator gsfIter=gsfElectrons->begin();
1758  gsfIter!=gsfElectrons->end(); gsfIter++){
1759  // preselect electrons
1760  if (gsfIter->pt()>maxPt_ || std::abs(gsfIter->eta())>maxAbsEta_) continue;
1761  h_ele_EoverP_all -> Fill( gsfIter->eSuperClusterOverP() );
1762  h_ele_EseedOP_all -> Fill( gsfIter->eSeedClusterOverP() );
1763  h_ele_EoPout_all -> Fill( gsfIter->eSeedClusterOverPout() );
1764  h_ele_EeleOPout_all -> Fill( gsfIter->eEleClusterOverPout() );
1765  h_ele_dEtaSc_propVtx_all -> Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
1766  h_ele_dPhiSc_propVtx_all -> Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
1767  h_ele_dEtaCl_propOut_all -> Fill(gsfIter->deltaEtaSeedClusterTrackAtCalo());
1768  h_ele_dPhiCl_propOut_all -> Fill(gsfIter->deltaPhiSeedClusterTrackAtCalo());
1769  h_ele_HoE_all -> Fill( gsfIter->hadronicOverEm() );
1770  double d = gsfIter->vertex().x()*gsfIter->vertex().x()+gsfIter->vertex().y()*gsfIter->vertex().y();
1771  h_ele_TIP_all -> Fill( sqrt(d) );
1772  h_ele_vertexEta_all -> Fill( gsfIter->eta() );
1773  h_ele_vertexPt_all -> Fill( gsfIter->pt() );
1774  h_ele_Et_all -> Fill( gsfIter->superCluster()->energy()/cosh(gsfIter->superCluster()->eta()));
1775  float enrj1=gsfIter->superCluster()->energy();
1776  // mee
1777  for (reco::GsfElectronCollection::const_iterator gsfIter2=gsfIter+1;
1778  gsfIter2!=gsfElectrons->end(); gsfIter2++){
1779  math::XYZTLorentzVector p12 = (*gsfIter).p4()+(*gsfIter2).p4();
1780  float mee2 = p12.Dot(p12);
1781  float enrj2=gsfIter2->superCluster()->energy();
1782  h_ele_mee_all -> Fill(sqrt(mee2));
1783  h_ele_E2mnE1vsMee_all->Fill(sqrt(mee2),enrj2-enrj1);
1784  if (gsfIter->ecalDrivenSeed() && gsfIter2->ecalDrivenSeed()) h_ele_E2mnE1vsMee_egeg_all->Fill(sqrt(mee2),enrj2-enrj1);
1785  if (gsfIter->charge()*gsfIter2->charge()<0.) {
1786  h_ele_mee_os -> Fill(sqrt(mee2));
1787  if (gsfIter->isEB() && gsfIter2->isEB()) h_ele_mee_os_ebeb -> Fill(sqrt(mee2));
1788  if ((gsfIter->isEB() && gsfIter2->isEE()) || (gsfIter->isEE() && gsfIter2->isEB())) h_ele_mee_os_ebee -> Fill(sqrt(mee2));
1789  if (gsfIter->isEE() && gsfIter2->isEE()) h_ele_mee_os_eeee -> Fill(sqrt(mee2));
1790  if ((gsfIter->classification()==GsfElectron::GOLDEN && gsfIter2->classification()==GsfElectron::GOLDEN) ||
1791  (gsfIter->classification()==GsfElectron::GOLDEN && gsfIter2->classification()==GsfElectron::BIGBREM) ||
1792  //(gsfIter->classification()==GsfElectron::GOLDEN && gsfIter2->classification()==GsfElectron::OLDNARROW) ||
1793  (gsfIter->classification()==GsfElectron::BIGBREM && gsfIter2->classification()==GsfElectron::GOLDEN) ||
1794  (gsfIter->classification()==GsfElectron::BIGBREM && gsfIter2->classification()==GsfElectron::BIGBREM)/* ||
1795  (gsfIter->classification()==GsfElectron::BIGBREM && gsfIter2->classification()==GsfElectron::OLDNARROW) ||
1796  (gsfIter->classification()==GsfElectron::OLDNARROW && gsfIter2->classification()==GsfElectron::GOLDEN) ||
1797  (gsfIter->classification()==GsfElectron::OLDNARROW && gsfIter2->classification()==GsfElectron::BIGBREM) ||
1798  (gsfIter->classification()==GsfElectron::OLDNARROW && gsfIter2->classification()==GsfElectron::OLDNARROW)*/ )
1799  { h_ele_mee_os_gg -> Fill(sqrt(mee2));}
1800  else if (
1801  (gsfIter->classification()==GsfElectron::SHOWERING && gsfIter2->classification()==GsfElectron::SHOWERING) ||
1802  (gsfIter->classification()==GsfElectron::SHOWERING && gsfIter2->isGap()) ||
1803  (gsfIter->isGap() && gsfIter2->classification()==GsfElectron::SHOWERING) ||
1804  (gsfIter->isGap() && gsfIter2->isGap()))
1805  { h_ele_mee_os_bb -> Fill(sqrt(mee2));}
1806  else
1807  { h_ele_mee_os_gb -> Fill(sqrt(mee2));}
1808  }
1809  }
1810  }
1811 
1812  int mcNum=0, gamNum=0, eleNum=0;
1813  bool matchingID, matchingMotherID;
1814 
1815  // charge mis-ID
1816  for (reco::GenParticleCollection::const_iterator mcIter=genParticles->begin(); mcIter != genParticles->end(); mcIter++ ) {
1817 
1818  // select requested matching gen particle
1819  matchingID=false;
1820  for (unsigned int i=0; i<matchingIDs_.size(); i++)
1821  if ( mcIter->pdgId() == matchingIDs_[i] ) matchingID=true;
1822 
1823  if (matchingID) {
1824 
1825  // select requested mother matching gen particle
1826  // always include single particle with no mother
1827  const Candidate * mother = mcIter->mother();
1828  matchingMotherID=false;
1829  for (unsigned int i=0; i<matchingMotherIDs_.size(); i++)
1830  if ((mother == 0) || ((mother != 0) && mother->pdgId() == matchingMotherIDs_[i]) ) matchingMotherID=true;
1831 
1832  if (matchingMotherID) {
1833 
1834  if (mcIter->pt()> maxPt_ || std::abs(mcIter->eta())> maxAbsEta_) continue;
1835 
1836  // suppress the endcaps
1837  //if (std::abs(mcIter->eta()) > 1.5) continue;
1838  // select central z
1839  //if ( std::abs(mcIter->production_vertex()->position().z())>50.) continue;
1840 
1841  // looking for the best matching gsf electron
1842  bool okGsfFound = false;
1843  double gsfOkRatio = 999999.;
1844 
1845  // find best matched electron
1846  reco::GsfElectron bestGsfElectron;
1847  for (reco::GsfElectronCollection::const_iterator gsfIter=gsfElectrons->begin();
1848  gsfIter!=gsfElectrons->end(); gsfIter++){
1849 
1850  double dphi = gsfIter->phi()-mcIter->phi();
1851  if (std::abs(dphi)>CLHEP::pi)
1852  dphi = dphi < 0? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
1853  double deltaR = sqrt(std::pow((gsfIter->eta()-mcIter->eta()),2) + std::pow(dphi,2));
1854  if ( deltaR < deltaR_ ){
1855  double mc_charge = mcIter->pdgId() == 11 ? -1. : 1. ;
1856  h_ele_ChargeMnChargeTrue -> Fill( std::abs(gsfIter->charge()-mc_charge));
1857  // require here a charge mismatch
1858  if ( ( (mcIter->pdgId() == 11) && (gsfIter->charge() > 0.) ) ||
1859  ( (mcIter->pdgId() == -11) && (gsfIter->charge() < 0.) ) )
1860  {
1861  double tmpGsfRatio = gsfIter->p()/mcIter->p();
1862  if ( std::abs(tmpGsfRatio-1) < std::abs(gsfOkRatio-1) ) {
1863  gsfOkRatio = tmpGsfRatio;
1864  bestGsfElectron=*gsfIter;
1865  okGsfFound = true;
1866  }
1867  }
1868  }
1869  } // loop over rec ele to look for the best one
1870 
1871  // analysis when the mc track is found
1872  if (okGsfFound)
1873  {
1874 
1875  // generated distributions for matched electrons
1876  h_ele_simPt_matched_qmisid -> Fill( mcIter->pt() );
1877  h_ele_simPhi_matched_qmisid -> Fill( mcIter->phi() );
1878  h_ele_simAbsEta_matched_qmisid -> Fill( std::abs(mcIter->eta()) );
1879  h_ele_simEta_matched_qmisid -> Fill( mcIter->eta() );
1880  h_ele_simZ_matched_qmisid -> Fill( mcIter->vz() );
1881 
1882  }
1883  }
1884  }
1885 
1886  }
1887 
1888  // association mc-reco
1889  for (reco::GenParticleCollection::const_iterator mcIter=genParticles->begin(); mcIter != genParticles->end(); mcIter++ ) {
1890 
1891  // number of mc particles
1892  mcNum++;
1893 
1894  // counts photons
1895  if (mcIter->pdgId() == 22 ){ gamNum++; }
1896 
1897  // select requested matching gen particle
1898  matchingID=false;
1899  for (unsigned int i=0; i<matchingIDs_.size(); i++)
1900  if ( mcIter->pdgId() == matchingIDs_[i] ) matchingID=true;
1901 
1902  if (matchingID) {
1903 
1904  // select requested mother matching gen particle
1905  // always include single particle with no mother
1906  const Candidate * mother = mcIter->mother();
1907  matchingMotherID=false;
1908  for (unsigned int i=0; i<matchingMotherIDs_.size(); i++)
1909  if ((mother == 0) || ((mother != 0) && mother->pdgId() == matchingMotherIDs_[i]) ) matchingMotherID=true;
1910 
1911  if (matchingMotherID) {
1912 
1913  if (mcIter->pt()> maxPt_ || std::abs(mcIter->eta())> maxAbsEta_) continue;
1914 
1915  // suppress the endcaps
1916  //if (std::abs(mcIter->eta()) > 1.5) continue;
1917  // select central z
1918  //if ( std::abs(mcIter->production_vertex()->position().z())>50.) continue;
1919 
1920  eleNum++;
1921  h_simEta -> Fill( mcIter->eta() );
1922  h_simAbsEta -> Fill( std::abs(mcIter->eta()) );
1923  h_simP -> Fill( mcIter->p() );
1924  h_simPt -> Fill( mcIter->pt() );
1925  h_simPhi -> Fill( mcIter->phi() );
1926  h_simZ -> Fill( mcIter->vz() );
1927  h_simPtEta -> Fill( mcIter->eta(),mcIter->pt() );
1928 
1929  // looking for the best matching gsf electron
1930  bool okGsfFound = false;
1931  double gsfOkRatio = 999999.;
1932 
1933  // find best matched electron
1934  reco::GsfElectron bestGsfElectron;
1935  for (reco::GsfElectronCollection::const_iterator gsfIter=gsfElectrons->begin();
1936  gsfIter!=gsfElectrons->end(); gsfIter++){
1937 
1938  double dphi = gsfIter->phi()-mcIter->phi();
1939  if (std::abs(dphi)>CLHEP::pi)
1940  dphi = dphi < 0? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
1941  double deltaR = sqrt(std::pow((gsfIter->eta()-mcIter->eta()),2) + std::pow(dphi,2));
1942  if ( deltaR < deltaR_ ){
1943  if ( ( (mcIter->pdgId() == 11) && (gsfIter->charge() < 0.) ) ||
1944  ( (mcIter->pdgId() == -11) && (gsfIter->charge() > 0.) ) )
1945  {
1946  double tmpGsfRatio = gsfIter->p()/mcIter->p();
1947  if ( std::abs(tmpGsfRatio-1) < std::abs(gsfOkRatio-1) ) {
1948  gsfOkRatio = tmpGsfRatio;
1949  bestGsfElectron=*gsfIter;
1950  okGsfFound = true;
1951  }
1952  }
1953  }
1954  } // loop over rec ele to look for the best one
1955 
1956  // analysis when the mc track is found
1957  if (okGsfFound){
1958 
1959  // electron related distributions
1960  h_ele_charge -> Fill( bestGsfElectron.charge() );
1961  h_ele_chargeVsEta -> Fill( bestGsfElectron.eta(),bestGsfElectron.charge() );
1962  h_ele_chargeVsPhi -> Fill( bestGsfElectron.phi(),bestGsfElectron.charge() );
1963  h_ele_chargeVsPt -> Fill( bestGsfElectron.pt(),bestGsfElectron.charge() );
1964  h_ele_vertexP -> Fill( bestGsfElectron.p() );
1965  h_ele_vertexPt -> Fill( bestGsfElectron.pt() );
1966  h_ele_Et -> Fill( bestGsfElectron.superCluster()->energy()/cosh(bestGsfElectron.superCluster()->eta()));
1967  h_ele_vertexPtVsEta -> Fill( bestGsfElectron.eta(),bestGsfElectron.pt() );
1968  h_ele_vertexPtVsPhi -> Fill( bestGsfElectron.phi(),bestGsfElectron.pt() );
1969  h_ele_vertexEta -> Fill( bestGsfElectron.eta() );
1970  // generated distributions for matched electrons
1971  h_ele_simPt_matched -> Fill( mcIter->pt() );
1972  h_ele_simPhi_matched -> Fill( mcIter->phi() );
1973  h_ele_simAbsEta_matched -> Fill( std::abs(mcIter->eta()) );
1974  h_ele_simEta_matched -> Fill( mcIter->eta() );
1975  h_ele_simPtEta_matched -> Fill( mcIter->eta(),mcIter->pt() );
1976  h_ele_vertexEtaVsPhi -> Fill( bestGsfElectron.phi(),bestGsfElectron.eta() );
1977  h_ele_vertexPhi -> Fill( bestGsfElectron.phi() );
1978  h_ele_vertexX -> Fill( bestGsfElectron.vertex().x() );
1979  h_ele_vertexY -> Fill( bestGsfElectron.vertex().y() );
1980  h_ele_vertexZ -> Fill( bestGsfElectron.vertex().z() );
1981  h_ele_simZ_matched -> Fill( mcIter->vz() );
1982  double d = (bestGsfElectron.vertex().x()-mcIter->vx())
1983  *(bestGsfElectron.vertex().x()-mcIter->vx())+
1984  (bestGsfElectron.vertex().y()-mcIter->vy())
1985  *(bestGsfElectron.vertex().y()-mcIter->vy());
1986  d = sqrt(d);
1987  h_ele_vertexTIP -> Fill( d );
1988  h_ele_vertexTIPVsEta -> Fill( bestGsfElectron.eta(), d );
1989  h_ele_vertexTIPVsPhi -> Fill( bestGsfElectron.phi(), d );
1990  h_ele_vertexTIPVsPt -> Fill( bestGsfElectron.pt(), d );
1991  h_ele_EtaMnEtaTrue -> Fill( bestGsfElectron.eta()-mcIter->eta());
1992  if (bestGsfElectron.isEB()) h_ele_EtaMnEtaTrue_barrel -> Fill( bestGsfElectron.eta()-mcIter->eta());
1993  if (bestGsfElectron.isEE()) h_ele_EtaMnEtaTrue_endcaps -> Fill( bestGsfElectron.eta()-mcIter->eta());
1994  h_ele_EtaMnEtaTrueVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.eta()-mcIter->eta());
1995  h_ele_EtaMnEtaTrueVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.eta()-mcIter->eta());
1996  h_ele_EtaMnEtaTrueVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.eta()-mcIter->eta());
1997  h_ele_PhiMnPhiTrue -> Fill( bestGsfElectron.phi()-mcIter->phi());
1998  if (bestGsfElectron.isEB()) h_ele_PhiMnPhiTrue_barrel -> Fill( bestGsfElectron.phi()-mcIter->phi());
1999  if (bestGsfElectron.isEE()) h_ele_PhiMnPhiTrue_endcaps -> Fill( bestGsfElectron.phi()-mcIter->phi());
2000  h_ele_PhiMnPhiTrue2 -> Fill( bestGsfElectron.phi()-mcIter->phi());
2001  h_ele_PhiMnPhiTrueVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.phi()-mcIter->phi());
2002  h_ele_PhiMnPhiTrueVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.phi()-mcIter->phi());
2003  h_ele_PhiMnPhiTrueVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.phi()-mcIter->phi());
2004  h_ele_PoPtrue -> Fill( bestGsfElectron.p()/mcIter->p());
2005  h_ele_PtoPttrue -> Fill( bestGsfElectron.pt()/mcIter->pt());
2006  h_ele_PoPtrueVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.p()/mcIter->p());
2007  h_ele_PoPtrueVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.p()/mcIter->p());
2008  h_ele_PoPtrueVsPt -> Fill( bestGsfElectron.py(), bestGsfElectron.p()/mcIter->p());
2009  if (bestGsfElectron.isEB()) h_ele_PoPtrue_barrel -> Fill( bestGsfElectron.p()/mcIter->p());
2010  if (bestGsfElectron.isEE()) h_ele_PoPtrue_endcaps -> Fill( bestGsfElectron.p()/mcIter->p());
2011  if (bestGsfElectron.isEB() && bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PoPtrue_golden_barrel -> Fill( bestGsfElectron.p()/mcIter->p());
2012  if (bestGsfElectron.isEE() && bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PoPtrue_golden_endcaps -> Fill( bestGsfElectron.p()/mcIter->p());
2013  if (bestGsfElectron.isEB() && bestGsfElectron.classification() == GsfElectron::SHOWERING) h_ele_PoPtrue_showering_barrel -> Fill( bestGsfElectron.p()/mcIter->p());
2014  if (bestGsfElectron.isEE() && bestGsfElectron.classification() == GsfElectron::SHOWERING) h_ele_PoPtrue_showering_endcaps -> Fill( bestGsfElectron.p()/mcIter->p());
2015  if (bestGsfElectron.isEB()) h_ele_PtoPttrue_barrel -> Fill( bestGsfElectron.pt()/mcIter->pt());
2016  if (bestGsfElectron.isEE()) h_ele_PtoPttrue_endcaps -> Fill( bestGsfElectron.pt()/mcIter->pt());
2017 
2018  // supercluster related distributions
2019  reco::SuperClusterRef sclRef = bestGsfElectron.superCluster();
2020  if (!bestGsfElectron.ecalDrivenSeed()&&bestGsfElectron.trackerDrivenSeed()) sclRef = bestGsfElectron.parentSuperCluster();
2021  histSclEn_->Fill(sclRef->energy());
2022  double R=TMath::Sqrt(sclRef->x()*sclRef->x() + sclRef->y()*sclRef->y() +sclRef->z()*sclRef->z());
2023  double Rt=TMath::Sqrt(sclRef->x()*sclRef->x() + sclRef->y()*sclRef->y());
2024  histSclEt_->Fill(sclRef->energy()*(Rt/R));
2025  histSclEtVsEta_->Fill(sclRef->eta(),sclRef->energy()*(Rt/R));
2026  histSclEtVsPhi_->Fill(sclRef->phi(),sclRef->energy()*(Rt/R));
2027  if (bestGsfElectron.isEB()) histSclEoEtrue_barrel->Fill(sclRef->energy()/mcIter->p());
2028  if (bestGsfElectron.isEE()) histSclEoEtrue_endcaps->Fill(sclRef->energy()/mcIter->p());
2029  if (bestGsfElectron.isEB() && bestGsfElectron.ecalDrivenSeed()) histSclEoEtrue_barrel_eg->Fill(sclRef->energy()/mcIter->p());
2030  if (bestGsfElectron.isEE() && bestGsfElectron.ecalDrivenSeed()) histSclEoEtrue_endcaps_eg->Fill(sclRef->energy()/mcIter->p());
2031  if (bestGsfElectron.isEB() && bestGsfElectron.isEBEtaGap()) histSclEoEtrue_barrel_etagap->Fill(sclRef->energy()/mcIter->p());
2032  if (bestGsfElectron.isEB() && bestGsfElectron.isEBPhiGap()) histSclEoEtrue_barrel_phigap->Fill(sclRef->energy()/mcIter->p());
2033  if (bestGsfElectron.isEBEEGap()) histSclEoEtrue_ebeegap->Fill(sclRef->energy()/mcIter->p());
2034  //if (bestGsfElectron.isEE()) histSclEoEtrue_endcaps->Fill(sclRef->energy()/mcIter->p());
2035  if (bestGsfElectron.isEE() && bestGsfElectron.isEEDeeGap()) histSclEoEtrue_endcaps_deegap->Fill(sclRef->energy()/mcIter->p());
2036  if (bestGsfElectron.isEE() && bestGsfElectron.isEERingGap()) histSclEoEtrue_endcaps_ringgap->Fill(sclRef->energy()/mcIter->p());
2037  if (bestGsfElectron.isEB()) histSclEoEtrue_barrel_new->Fill(sclRef->energy()/mcIter->p());
2038  if (bestGsfElectron.isEE()) histSclEoEtrue_endcaps_new->Fill(sclRef->energy()/mcIter->p());
2039  if (bestGsfElectron.isEB() && bestGsfElectron.ecalDrivenSeed()) histSclEoEtrue_barrel_eg_new->Fill(sclRef->energy()/mcIter->p());
2040  if (bestGsfElectron.isEE() && bestGsfElectron.ecalDrivenSeed()) histSclEoEtrue_endcaps_eg_new->Fill(sclRef->energy()/mcIter->p());
2041  if (bestGsfElectron.isEB() && bestGsfElectron.isEBEtaGap()) histSclEoEtrue_barrel_etagap_new->Fill(sclRef->energy()/mcIter->p());
2042  if (bestGsfElectron.isEB() && bestGsfElectron.isEBPhiGap()) histSclEoEtrue_barrel_phigap_new->Fill(sclRef->energy()/mcIter->p());
2043  if (bestGsfElectron.isEBEEGap()) histSclEoEtrue_ebeegap_new->Fill(sclRef->energy()/mcIter->p());
2044  //if (bestGsfElectron.isEE()) histSclEoEtrue_endcaps_new->Fill(sclRef->energy()/mcIter->p());
2045  if (bestGsfElectron.isEE() && bestGsfElectron.isEEDeeGap()) histSclEoEtrue_endcaps_deegap_new->Fill(sclRef->energy()/mcIter->p());
2046  if (bestGsfElectron.isEE() && bestGsfElectron.isEERingGap()) histSclEoEtrue_endcaps_ringgap_new->Fill(sclRef->energy()/mcIter->p());
2047  histSclEta_->Fill(sclRef->eta());
2048  histSclEtaVsPhi_->Fill(sclRef->phi(),sclRef->eta());
2049  histSclPhi_->Fill(sclRef->phi());
2050  histSclSigEtaEta_->Fill(bestGsfElectron.scSigmaEtaEta());
2051  if (bestGsfElectron.isEB()) histSclSigEtaEta_barrel_->Fill(bestGsfElectron.scSigmaEtaEta());
2052  if (bestGsfElectron.isEE()) histSclSigEtaEta_endcaps_->Fill(bestGsfElectron.scSigmaEtaEta());
2053  histSclSigIEtaIEta_->Fill(bestGsfElectron.scSigmaIEtaIEta());
2054  if (bestGsfElectron.isEB()) histSclSigIEtaIEta_barrel_->Fill(bestGsfElectron.scSigmaIEtaIEta());
2055  if (bestGsfElectron.isEE()) histSclSigIEtaIEta_endcaps_->Fill(bestGsfElectron.scSigmaIEtaIEta());
2056  histSclE1x5_->Fill(bestGsfElectron.scE1x5());
2057  if (bestGsfElectron.isEB()) histSclE1x5_barrel_->Fill(bestGsfElectron.scE1x5());
2058  if (bestGsfElectron.isEE()) histSclE1x5_endcaps_->Fill(bestGsfElectron.scE1x5());
2059  histSclE2x5max_->Fill(bestGsfElectron.scE2x5Max());
2060  if (bestGsfElectron.isEB()) histSclE2x5max_barrel_->Fill(bestGsfElectron.scE2x5Max());
2061  if (bestGsfElectron.isEE()) histSclE2x5max_endcaps_->Fill(bestGsfElectron.scE2x5Max());
2062  histSclE5x5_->Fill(bestGsfElectron.scE5x5());
2063  if (bestGsfElectron.isEB()) histSclE5x5_barrel_->Fill(bestGsfElectron.scE5x5());
2064  if (bestGsfElectron.isEE()) histSclE5x5_endcaps_->Fill(bestGsfElectron.scE5x5());
2065  if (bestGsfElectron.ecalDrivenSeed()) histSclSigIEtaIEta_eg_->Fill(bestGsfElectron.scSigmaIEtaIEta());
2066  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) histSclSigIEtaIEta_eg_barrel_->Fill(bestGsfElectron.scSigmaIEtaIEta());
2067  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) histSclSigIEtaIEta_eg_endcaps_->Fill(bestGsfElectron.scSigmaIEtaIEta());
2068  if (bestGsfElectron.ecalDrivenSeed())histSclE1x5_eg_->Fill(bestGsfElectron.scE1x5());
2069  if (bestGsfElectron.isEB() && bestGsfElectron.ecalDrivenSeed())histSclE1x5_eg_barrel_->Fill(bestGsfElectron.scE1x5());
2070  if (bestGsfElectron.isEE() && bestGsfElectron.ecalDrivenSeed())histSclE1x5_eg_endcaps_->Fill(bestGsfElectron.scE1x5());
2071  if (bestGsfElectron.ecalDrivenSeed())histSclE2x5max_eg_->Fill(bestGsfElectron.scE2x5Max());
2072  if (bestGsfElectron.isEB() && bestGsfElectron.ecalDrivenSeed())histSclE2x5max_eg_barrel_->Fill(bestGsfElectron.scE2x5Max());
2073  if (bestGsfElectron.isEE() && bestGsfElectron.ecalDrivenSeed())histSclE2x5max_eg_endcaps_->Fill(bestGsfElectron.scE2x5Max());
2074  if (bestGsfElectron.ecalDrivenSeed())histSclE5x5_eg_->Fill(bestGsfElectron.scE5x5());
2075  if (bestGsfElectron.isEB() && bestGsfElectron.ecalDrivenSeed())histSclE5x5_eg_barrel_->Fill(bestGsfElectron.scE5x5());
2076  if (bestGsfElectron.isEE() && bestGsfElectron.ecalDrivenSeed())histSclE5x5_eg_endcaps_->Fill(bestGsfElectron.scE5x5());
2077  float pfEnergy=0., egEnergy=0.;
2078  if (!bestGsfElectron.superCluster().isNull()) egEnergy = bestGsfElectron.superCluster()->energy();
2079  if (!bestGsfElectron.parentSuperCluster().isNull()) pfEnergy = bestGsfElectron.parentSuperCluster()->energy();
2080  histSclEoEtruePfVsEg->Fill(egEnergy/mcIter->p(),pfEnergy/mcIter->p());
2081 
2082  // track related distributions
2083  h_ele_ambiguousTracks -> Fill( bestGsfElectron.ambiguousGsfTracksSize() );
2084  h_ele_ambiguousTracksVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.ambiguousGsfTracksSize() );
2085  h_ele_ambiguousTracksVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.ambiguousGsfTracksSize() );
2086  h_ele_ambiguousTracksVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.ambiguousGsfTracksSize() );
2087  if (!readAOD_) { // track extra does not exist in AOD
2088  h_ele_foundHits -> Fill( bestGsfElectron.gsfTrack()->numberOfValidHits() );
2089  if (bestGsfElectron.isEB()) h_ele_foundHits_barrel -> Fill( bestGsfElectron.gsfTrack()->numberOfValidHits() );
2090  if (bestGsfElectron.isEE()) h_ele_foundHits_endcaps -> Fill( bestGsfElectron.gsfTrack()->numberOfValidHits() );
2091  h_ele_foundHitsVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->numberOfValidHits() );
2092  h_ele_foundHitsVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->numberOfValidHits() );
2093  h_ele_foundHitsVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->numberOfValidHits() );
2094  h_ele_lostHits -> Fill( bestGsfElectron.gsfTrack()->numberOfLostHits() );
2095  if (bestGsfElectron.isEB()) h_ele_lostHits_barrel -> Fill( bestGsfElectron.gsfTrack()->numberOfLostHits() );
2096  if (bestGsfElectron.isEE()) h_ele_lostHits_endcaps -> Fill( bestGsfElectron.gsfTrack()->numberOfLostHits() );
2097  h_ele_lostHitsVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->numberOfLostHits() );
2098  h_ele_lostHitsVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->numberOfLostHits() );
2099  h_ele_lostHitsVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->numberOfLostHits() );
2100  h_ele_chi2 -> Fill( bestGsfElectron.gsfTrack()->normalizedChi2() );
2101  if (bestGsfElectron.isEB()) h_ele_chi2_barrel -> Fill( bestGsfElectron.gsfTrack()->normalizedChi2() );
2102  if (bestGsfElectron.isEE()) h_ele_chi2_endcaps -> Fill( bestGsfElectron.gsfTrack()->normalizedChi2() );
2103  h_ele_chi2VsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->normalizedChi2() );
2104  h_ele_chi2VsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->normalizedChi2() );
2105  h_ele_chi2VsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->normalizedChi2() );
2106  }
2107  // from gsf track interface, hence using mean
2108  if (!readAOD_) { // track extra does not exist in AOD
2109  h_ele_PinMnPout -> Fill( bestGsfElectron.gsfTrack()->innerMomentum().R() - bestGsfElectron.gsfTrack()->outerMomentum().R() );
2110  h_ele_outerP -> Fill( bestGsfElectron.gsfTrack()->outerMomentum().R() );
2111  h_ele_outerPt -> Fill( bestGsfElectron.gsfTrack()->outerMomentum().Rho() );
2112  }
2113  // from electron interface, hence using mode
2114  h_ele_PinMnPout_mode -> Fill( bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
2115  h_ele_PinMnPoutVsEta_mode -> Fill( bestGsfElectron.eta(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
2116  h_ele_PinMnPoutVsPhi_mode -> Fill( bestGsfElectron.phi(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
2117  h_ele_PinMnPoutVsPt_mode -> Fill( bestGsfElectron.pt(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
2118  h_ele_PinMnPoutVsE_mode -> Fill( bestGsfElectron.caloEnergy(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
2119  if (!readAOD_) // track extra does not exist in AOD
2120  h_ele_PinMnPoutVsChi2_mode -> Fill( bestGsfElectron.gsfTrack()->normalizedChi2(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
2121  h_ele_outerP_mode -> Fill( bestGsfElectron.trackMomentumOut().R() );
2122  h_ele_outerPVsEta_mode -> Fill(bestGsfElectron.eta(), bestGsfElectron.trackMomentumOut().R() );
2123  h_ele_outerPt_mode -> Fill( bestGsfElectron.trackMomentumOut().Rho() );
2124  h_ele_outerPtVsEta_mode -> Fill(bestGsfElectron.eta(), bestGsfElectron.trackMomentumOut().Rho() );
2125  h_ele_outerPtVsPhi_mode -> Fill(bestGsfElectron.phi(), bestGsfElectron.trackMomentumOut().Rho() );
2126  h_ele_outerPtVsPt_mode -> Fill(bestGsfElectron.pt(), bestGsfElectron.trackMomentumOut().Rho() );
2127 
2128  if (!readAOD_) { // track extra does not exist in AOD
2129  edm::RefToBase<TrajectorySeed> seed = bestGsfElectron.gsfTrack()->extra()->seedRef();
2130  ElectronSeedRef elseed=seed.castTo<ElectronSeedRef>();
2131  h_ele_seed_dphi2_-> Fill(elseed->dPhi2());
2132  h_ele_seed_dphi2VsEta_-> Fill(bestGsfElectron.eta(), elseed->dPhi2());
2133  h_ele_seed_dphi2VsPt_-> Fill(bestGsfElectron.pt(), elseed->dPhi2()) ;
2134  h_ele_seed_drz2_-> Fill(elseed->dRz2());
2135  h_ele_seed_drz2VsEta_-> Fill(bestGsfElectron.eta(), elseed->dRz2());
2136  h_ele_seed_drz2VsPt_-> Fill(bestGsfElectron.pt(), elseed->dRz2());
2137  h_ele_seed_subdet2_-> Fill(elseed->subDet2());
2138  }
2139  // match distributions
2140  h_ele_EoP -> Fill( bestGsfElectron.eSuperClusterOverP() );
2141  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EoP_eg -> Fill( bestGsfElectron.eSuperClusterOverP() );
2142  if (bestGsfElectron.isEB()) h_ele_EoP_barrel -> Fill( bestGsfElectron.eSuperClusterOverP() );
2143  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoP_eg_barrel -> Fill( bestGsfElectron.eSuperClusterOverP() );
2144  if (bestGsfElectron.isEE()) h_ele_EoP_endcaps -> Fill( bestGsfElectron.eSuperClusterOverP() );
2145  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoP_eg_endcaps -> Fill( bestGsfElectron.eSuperClusterOverP() );
2146  h_ele_EoPVsEta -> Fill(bestGsfElectron.eta(), bestGsfElectron.eSuperClusterOverP() );
2147  h_ele_EoPVsPhi -> Fill(bestGsfElectron.phi(), bestGsfElectron.eSuperClusterOverP() );
2148  h_ele_EoPVsE -> Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.eSuperClusterOverP() );
2149  h_ele_EseedOP -> Fill( bestGsfElectron.eSeedClusterOverP() );
2150  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EseedOP_eg -> Fill( bestGsfElectron.eSeedClusterOverP() );
2151  if (bestGsfElectron.isEB()) h_ele_EseedOP_barrel -> Fill( bestGsfElectron.eSeedClusterOverP() );
2152  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EseedOP_eg_barrel -> Fill( bestGsfElectron.eSeedClusterOverP() );
2153  if (bestGsfElectron.isEE()) h_ele_EseedOP_endcaps -> Fill( bestGsfElectron.eSeedClusterOverP() );
2154  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EseedOP_eg_endcaps -> Fill( bestGsfElectron.eSeedClusterOverP() );
2155  h_ele_EseedOPVsEta -> Fill(bestGsfElectron.eta(), bestGsfElectron.eSeedClusterOverP() );
2156  h_ele_EseedOPVsPhi -> Fill(bestGsfElectron.phi(), bestGsfElectron.eSeedClusterOverP() );
2157  h_ele_EseedOPVsE -> Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.eSeedClusterOverP() );
2158  h_ele_EoPout -> Fill( bestGsfElectron.eSeedClusterOverPout() );
2159  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EoPout_eg -> Fill( bestGsfElectron.eSeedClusterOverPout() );
2160  if (bestGsfElectron.isEB()) h_ele_EoPout_barrel -> Fill( bestGsfElectron.eSeedClusterOverPout() );
2161  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoPout_eg_barrel -> Fill( bestGsfElectron.eSeedClusterOverPout() );
2162  if (bestGsfElectron.isEE()) h_ele_EoPout_endcaps -> Fill( bestGsfElectron.eSeedClusterOverPout() );
2163  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoPout_eg_endcaps -> Fill( bestGsfElectron.eSeedClusterOverPout() );
2164  h_ele_EoPoutVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.eSeedClusterOverPout() );
2165  h_ele_EoPoutVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.eSeedClusterOverPout() );
2166  h_ele_EoPoutVsE -> Fill( bestGsfElectron.caloEnergy(), bestGsfElectron.eSeedClusterOverPout() );
2167  h_ele_EeleOPout -> Fill( bestGsfElectron.eEleClusterOverPout() );
2168  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EeleOPout_eg -> Fill( bestGsfElectron.eEleClusterOverPout() );
2169  if (bestGsfElectron.isEB()) h_ele_EeleOPout_barrel -> Fill( bestGsfElectron.eEleClusterOverPout() );
2170  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EeleOPout_eg_barrel -> Fill( bestGsfElectron.eEleClusterOverPout() );
2171  if (bestGsfElectron.isEE()) h_ele_EeleOPout_endcaps -> Fill( bestGsfElectron.eEleClusterOverPout() );
2172  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EeleOPout_eg_endcaps -> Fill( bestGsfElectron.eEleClusterOverPout() );
2173  h_ele_EeleOPoutVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.eEleClusterOverPout() );
2174  h_ele_EeleOPoutVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.eEleClusterOverPout() );
2175  h_ele_EeleOPoutVsE -> Fill( bestGsfElectron.caloEnergy(), bestGsfElectron.eEleClusterOverPout() );
2176  h_ele_dEtaSc_propVtx -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2177  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaSc_propVtx_eg -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2178  if (bestGsfElectron.isEB()) h_ele_dEtaSc_propVtx_barrel -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2179  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaSc_propVtx_eg_barrel -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2180  if (bestGsfElectron.isEE())h_ele_dEtaSc_propVtx_endcaps -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2181  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaSc_propVtx_eg_endcaps -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2182  h_ele_dEtaScVsEta_propVtx -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2183  h_ele_dEtaScVsPhi_propVtx -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2184  h_ele_dEtaScVsPt_propVtx -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
2185  h_ele_dPhiSc_propVtx -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2186  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiSc_propVtx_eg -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2187  if (bestGsfElectron.isEB()) h_ele_dPhiSc_propVtx_barrel -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2188  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiSc_propVtx_eg_barrel -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2189  if (bestGsfElectron.isEE())h_ele_dPhiSc_propVtx_endcaps -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2190  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiSc_propVtx_eg_endcaps -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2191  h_ele_dPhiScVsEta_propVtx -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2192  h_ele_dPhiScVsPhi_propVtx -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2193  h_ele_dPhiScVsPt_propVtx -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
2194  h_ele_dEtaCl_propOut -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2195  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaCl_propOut_eg -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2196  if (bestGsfElectron.isEB()) h_ele_dEtaCl_propOut_barrel -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2197  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2198  if (bestGsfElectron.isEE()) h_ele_dEtaCl_propOut_endcaps -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2199  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2200  h_ele_dEtaClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2201  h_ele_dEtaClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2202  h_ele_dEtaClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
2203  h_ele_dPhiCl_propOut -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2204  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiCl_propOut_eg -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2205  if (bestGsfElectron.isEB()) h_ele_dPhiCl_propOut_barrel -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2206  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2207  if (bestGsfElectron.isEE()) h_ele_dPhiCl_propOut_endcaps -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2208  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2209  h_ele_dPhiClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2210  h_ele_dPhiClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2211  h_ele_dPhiClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
2212  h_ele_dEtaEleCl_propOut -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2213  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaEleCl_propOut_eg -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2214  if (bestGsfElectron.isEB()) h_ele_dEtaEleCl_propOut_barrel -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2215  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaEleCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2216  if (bestGsfElectron.isEE()) h_ele_dEtaEleCl_propOut_endcaps -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2217  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaEleCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2218  h_ele_dEtaEleClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2219  h_ele_dEtaEleClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2220  h_ele_dEtaEleClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
2221  h_ele_dPhiEleCl_propOut -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2222  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiEleCl_propOut_eg -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2223  if (bestGsfElectron.isEB()) h_ele_dPhiEleCl_propOut_barrel -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2224  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiEleCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2225  if (bestGsfElectron.isEE()) h_ele_dPhiEleCl_propOut_endcaps -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2226  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiEleCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2227  h_ele_dPhiEleClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2228  h_ele_dPhiEleClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2229  h_ele_dPhiEleClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
2230  h_ele_HoE -> Fill(bestGsfElectron.hadronicOverEm());
2231  if (bestGsfElectron.ecalDrivenSeed()) h_ele_HoE_eg -> Fill(bestGsfElectron.hadronicOverEm());
2232  if (bestGsfElectron.isEB()) h_ele_HoE_barrel -> Fill(bestGsfElectron.hadronicOverEm());
2233  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_HoE_eg_barrel -> Fill(bestGsfElectron.hadronicOverEm());
2234  if (bestGsfElectron.isEE()) h_ele_HoE_endcaps -> Fill(bestGsfElectron.hadronicOverEm());
2235  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_HoE_eg_endcaps -> Fill(bestGsfElectron.hadronicOverEm());
2236  if (!bestGsfElectron.isEBEtaGap() && !bestGsfElectron.isEBPhiGap()&& !bestGsfElectron.isEBEEGap() &&
2237  !bestGsfElectron.isEERingGap() && !bestGsfElectron.isEEDeeGap()) h_ele_HoE_fiducial -> Fill(bestGsfElectron.hadronicOverEm());
2238  h_ele_HoEVsEta -> Fill( bestGsfElectron.eta(),bestGsfElectron.hadronicOverEm());
2239  h_ele_HoEVsPhi -> Fill(bestGsfElectron.phi(),bestGsfElectron.hadronicOverEm());
2240  h_ele_HoEVsE -> Fill(bestGsfElectron.caloEnergy(),bestGsfElectron.hadronicOverEm());
2241 
2242  //classes
2243  int eleClass = bestGsfElectron.classification();
2244  if (bestGsfElectron.isEE()) eleClass+=10;
2245  h_ele_classes ->Fill(eleClass);
2246 
2247  if (bestGsfElectron.classification() == GsfElectron::GOLDEN && bestGsfElectron.isEB()) histSclEoEtrueGolden_barrel->Fill(sclRef->energy()/mcIter->p());
2248  if (bestGsfElectron.classification() == GsfElectron::GOLDEN && bestGsfElectron.isEE()) histSclEoEtrueGolden_endcaps->Fill(sclRef->energy()/mcIter->p());
2249  if (bestGsfElectron.classification() == GsfElectron::SHOWERING && bestGsfElectron.isEB()) histSclEoEtrueShowering_barrel->Fill(sclRef->energy()/mcIter->p());
2250  if (bestGsfElectron.classification() == GsfElectron::SHOWERING && bestGsfElectron.isEE()) histSclEoEtrueShowering_endcaps->Fill(sclRef->energy()/mcIter->p());
2251 
2252  //eleClass = eleClass%100; // get rid of barrel/endcap distinction
2253  h_ele_eta->Fill(std::abs(bestGsfElectron.eta()));
2254  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_eta_golden ->Fill(std::abs(bestGsfElectron.eta()));
2255  if (bestGsfElectron.classification() == GsfElectron::BIGBREM) h_ele_eta_bbrem ->Fill(std::abs(bestGsfElectron.eta()));
2256  //if (bestGsfElectron.classification() == GsfElectron::NARROW) h_ele_eta_narrow ->Fill(std::abs(bestGsfElectron.eta()));
2257  if (bestGsfElectron.classification() == GsfElectron::SHOWERING) h_ele_eta_shower ->Fill(std::abs(bestGsfElectron.eta()));
2258 
2259  //fbrem
2260  double fbrem_mean=0.;
2261  if (!readAOD_) // track extra does not exist in AOD
2262  fbrem_mean = 1. - bestGsfElectron.gsfTrack()->outerMomentum().R()/bestGsfElectron.gsfTrack()->innerMomentum().R();
2263  double fbrem_mode = bestGsfElectron.fbrem();
2264  h_ele_fbrem->Fill(fbrem_mode);
2265  if (bestGsfElectron.ecalDrivenSeed()) h_ele_fbrem_eg->Fill(fbrem_mode);
2266  h_ele_fbremVsEta_mode->Fill(bestGsfElectron.eta(),fbrem_mode);
2267  if (!readAOD_) // track extra does not exist in AOD
2268  h_ele_fbremVsEta_mean->Fill(bestGsfElectron.eta(),fbrem_mean);
2269 
2270  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PinVsPoutGolden_mode -> Fill(bestGsfElectron.trackMomentumOut().R(), bestGsfElectron.trackMomentumAtVtx().R());
2271  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
2272  h_ele_PinVsPoutShowering_mode -> Fill(bestGsfElectron.trackMomentumOut().R(), bestGsfElectron.trackMomentumAtVtx().R());
2273  if (!readAOD_) { // track extra not available in AOD
2274  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PinVsPoutGolden_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().R(), bestGsfElectron.gsfTrack()->innerMomentum().R());
2275  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
2276  h_ele_PinVsPoutShowering_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().R(), bestGsfElectron.gsfTrack()->innerMomentum().R());
2277  }
2278  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PtinVsPtoutGolden_mode -> Fill(bestGsfElectron.trackMomentumOut().Rho(), bestGsfElectron.trackMomentumAtVtx().Rho());
2279  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
2280  h_ele_PtinVsPtoutShowering_mode -> Fill(bestGsfElectron.trackMomentumOut().Rho(), bestGsfElectron.trackMomentumAtVtx().Rho());
2281  if (!readAOD_) { // track extra not available in AOD
2282  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PtinVsPtoutGolden_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().Rho(), bestGsfElectron.gsfTrack()->innerMomentum().Rho());
2283  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
2284  h_ele_PtinVsPtoutShowering_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().Rho(), bestGsfElectron.gsfTrack()->innerMomentum().Rho());
2285  }
2286 
2287  h_ele_mva->Fill(bestGsfElectron.mva_e_pi());
2288  if (bestGsfElectron.ecalDrivenSeed()) h_ele_mva_eg->Fill(bestGsfElectron.mva_e_pi());
2289  if (bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(1.);
2290  if (bestGsfElectron.trackerDrivenSeed()) h_ele_provenance->Fill(-1.);
2291  if (bestGsfElectron.trackerDrivenSeed()||bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(0.);
2292  if (bestGsfElectron.trackerDrivenSeed()&&!bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(-2.);
2293  if (!bestGsfElectron.trackerDrivenSeed()&&bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(2.);
2294 
2295  h_ele_tkSumPt_dr03->Fill(bestGsfElectron.dr03TkSumPt());
2296  h_ele_ecalRecHitSumEt_dr03->Fill(bestGsfElectron.dr03EcalRecHitSumEt());
2297  h_ele_hcalDepth1TowerSumEt_dr03->Fill(bestGsfElectron.dr03HcalDepth1TowerSumEt());
2298  h_ele_hcalDepth2TowerSumEt_dr03->Fill(bestGsfElectron.dr03HcalDepth2TowerSumEt());
2299  h_ele_tkSumPt_dr04->Fill(bestGsfElectron.dr04TkSumPt());
2300  h_ele_ecalRecHitSumEt_dr04->Fill(bestGsfElectron.dr04EcalRecHitSumEt());
2301  h_ele_hcalDepth1TowerSumEt_dr04->Fill(bestGsfElectron.dr04HcalDepth1TowerSumEt());
2302  h_ele_hcalDepth2TowerSumEt_dr04->Fill(bestGsfElectron.dr04HcalDepth2TowerSumEt());
2303 
2304  } // gsf electron found
2305 
2306  } // mc particle found
2307 
2308  }
2309 
2310  } // loop over mc particle
2311 
2312  h_mcNum->Fill(mcNum);
2313  h_eleNum->Fill(eleNum);
2314 
2315 }
2316 
2317 
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
i
int i
Definition: DBlmapReader.cc:9
reco::LeafCandidate::p
virtual double p() const
magnitude of momentum vector
Definition: LeafCandidate.h:108
reco::GsfElectron::parentSuperCluster
SuperClusterRef parentSuperCluster() const
Definition: GsfElectron.h:188
reco::GsfElectron::scSigmaEtaEta
float scSigmaEtaEta() const
Definition: GsfElectron.h:438
reco::GsfElectron::isEBEtaGap
bool isEBEtaGap() const
Definition: GsfElectron.h:355
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:243
reco::GsfElectron::isEBEEGap
bool isEBEEGap() const
Definition: GsfElectron.h:353
reco::GsfElectron::dr04HcalDepth1TowerSumEt
float dr04HcalDepth1TowerSumEt() const
Definition: GsfElectron.h:490
GsfElectronMCAnalyzer::GsfElectronMCAnalyzer
GsfElectronMCAnalyzer(const edm::ParameterSet &conf)
Definition: GsfElectronMCAnalyzer.cc:51
reco::LeafCandidate::vertex
virtual const Point & vertex() const
vertex position (overwritten by PF...)
Definition: LeafCandidate.h:150
reco::GsfElectron::trackMomentumAtVtx
math::XYZVectorF trackMomentumAtVtx() const
Definition: GsfElectron.h:289
create_public_lumi_plots.log
log
Definition: create_public_lumi_plots.py:1102
MCTruthHelper::mother
const reco::GenParticle * mother(const reco::GenParticle &p, unsigned int imoth=0)
Definition: MCTruthHelper.h:594
deltaR_
double deltaR_
Definition: RecoTauIsolationMasking.cc:46
reco::GsfElectron::dr04TkSumPt
float dr04TkSumPt() const
Definition: GsfElectron.h:488
reco::GsfElectron::isEERingGap
bool isEERingGap() const
Definition: GsfElectron.h:359
reco::GsfElectron::dr03HcalDepth2TowerSumEt
float dr03HcalDepth2TowerSumEt() const
Definition: GsfElectron.h:480
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:684
reco::deltaR
double deltaR(const T1 &t1, const T2 &t2)
Definition: deltaR.h:48
reco::GsfElectron::scE2x5Max
float scE2x5Max() const
Definition: GsfElectron.h:441
genParticleCandidates2GenParticles_cfi.genParticles
tuple genParticles
Definition: genParticleCandidates2GenParticles_cfi.py:3
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107
reco::GsfElectron::isEEDeeGap
bool isEEDeeGap() const
Definition: GsfElectron.h:358
reco::LeafCandidate::eta
virtual double eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:137
reco::LeafCandidate::pt
virtual double pt() const
transverse momentum
Definition: LeafCandidate.h:131
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
reco::GsfElectron::isEE
bool isEE() const
Definition: GsfElectron.h:351
reco::GsfElectron::isEB
bool isEB() const
Definition: GsfElectron.h:350
pi
const Double_t pi
Definition: trackSplitPlot.h:36
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
reco::GsfElectron::dr04HcalDepth2TowerSumEt
float dr04HcalDepth2TowerSumEt() const
Definition: GsfElectron.h:491
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:247
iEvent
int iEvent
Definition: GenABIO.cc:230
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:443
reco::GsfElectron::eSeedClusterOverP
float eSeedClusterOverP() const
Definition: GsfElectron.h:244
reco::GsfElectron::dr04EcalRecHitSumEt
float dr04EcalRecHitSumEt() const
Definition: GsfElectron.h:489
reco::GsfElectron::deltaPhiSeedClusterTrackAtCalo
float deltaPhiSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:251
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
reco::GsfElectron::superCluster
virtual SuperClusterRef superCluster() const
reference to a SuperCluster
Definition: GsfElectron.h:182
reco::GsfElectron::trackMomentumOut
math::XYZVectorF trackMomentumOut() const
Definition: GsfElectron.h:291
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:250
reco::LeafCandidate::charge
virtual int charge() const
electric charge
Definition: LeafCandidate.h:91
reco::GsfElectron::dr03TkSumPt
float dr03TkSumPt() const
Definition: GsfElectron.h:477
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
reco::GsfElectron::dr03HcalDepth1TowerSumEt
float dr03HcalDepth1TowerSumEt() const
Definition: GsfElectron.h:479
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:246
reco::GsfElectron::eSeedClusterOverPout
float eSeedClusterOverPout() const
Definition: GsfElectron.h:245
GsfElectronMCAnalyzer::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: GsfElectronMCAnalyzer.cc:1742
reco::GsfElectron::scSigmaIEtaIEta
float scSigmaIEtaIEta() const
Definition: GsfElectron.h:439
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:405
edm::Ref::isNull
bool isNull() const
Checks for null.
Definition: Ref.h:247
reco::GsfElectron::deltaPhiEleClusterTrackAtCalo
float deltaPhiEleClusterTrackAtCalo() const
Definition: GsfElectron.h:252
dbtoconf.conf
tuple conf
Definition: dbtoconf.py:185
reco::Candidate::pdgId
virtual int pdgId() const =0
PDG identifier.
edm::Handle::product
T const * product() const
Definition: Handle.h:81
reco::GsfElectron::deltaEtaEleClusterTrackAtCalo
float deltaEtaEleClusterTrackAtCalo() const
Definition: GsfElectron.h:249
reco::GsfElectron::classification
Classification classification() const
Definition: GsfElectron.h:680
edm::RefToBase::castTo
REF castTo() const
cast to a concrete type
Definition: RefToBase.h:241
reco::GsfElectron::ambiguousGsfTracksSize
GsfTrackRefVector::size_type ambiguousGsfTracksSize() const
Definition: GsfElectron.h:635
reco::GsfElectron::dr03EcalRecHitSumEt
float dr03EcalRecHitSumEt() const
Definition: GsfElectron.h:478
reco::GsfElectron::mva_e_pi
float mva_e_pi() const
Definition: GsfElectron.h:615
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:60
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:248
reco::GsfElectron::isEBPhiGap
bool isEBPhiGap() const
Definition: GsfElectron.h:356
reco::GsfElectron::trackerDrivenSeed
bool trackerDrivenSeed() const
Definition: GsfElectron.h:187
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
reco::GsfElectron::scE5x5
float scE5x5() const
Definition: GsfElectron.h:442
reco::GsfElectron::caloEnergy
float caloEnergy() const
Definition: GsfElectron.h:774
reco::LeafCandidate::phi
virtual double phi() const
momentum azimuthal angle
Definition: LeafCandidate.h:133
reco::CompositeRefCandidateT::mother
virtual const Candidate * mother(size_type=0) const
return mother at a given position, i = 0, ... numberOfMothers() - 1 (read only mode) ...
reco::LeafCandidate::py
virtual double py() const
y coordinate of momentum vector
Definition: LeafCandidate.h:127
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
reco::GsfElectron::scE1x5
float scE1x5() const
Definition: GsfElectron.h:440
reco::GsfElectron::gsfTrack
virtual GsfTrackRef gsfTrack() const
reference to a GsfTrack
Definition: GsfElectron.h:183
reco::GsfElectron::ecalDrivenSeed
bool ecalDrivenSeed() const
Definition: GsfElectron.h:186
gsfElectrons_cfi.gsfElectrons
tuple gsfElectrons
Definition: gsfElectrons_cfi.py:174
Generated for CMSSW Reference Manual by doxygen 1.8.5