CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
GsfElectronMCFakeAnalyzer.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: RecoEgamma/Examples
4 // Class: GsfElectronMCFakeAnalyzer
5 //
13 //
14 // Original Author: Ursula Berthon
15 // Created: Mon Mar 27 13:22:06 CEST 2006
16 // $Id: GsfElectronMCFakeAnalyzer.cc,v 1.9 2012/09/13 20:08:32 wdd Exp $
17 //
18 //
19 
20 // user include files
21 #include "RecoEgamma/Examples/plugins/GsfElectronMCFakeAnalyzer.h"
22 
23 #include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
24 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
25 #include "DataFormats/EgammaCandidates/interface/GsfElectronFwd.h"
26 #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
27 #include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
28 #include "DataFormats/EgammaReco/interface/ElectronSeed.h"
29 #include "DataFormats/EgammaReco/interface/ElectronSeedFwd.h"
30 #include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
31 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
32 #include "DataFormats/JetReco/interface/GenJetCollection.h"
33 
34 #include "FWCore/ParameterSet/interface/ParameterSet.h"
35 #include "FWCore/Framework/interface/EDAnalyzer.h"
36 #include "FWCore/Framework/interface/Event.h"
37 #include "FWCore/Framework/interface/MakerMacros.h"
38 #include "FWCore/MessageLogger/interface/MessageLogger.h"
39 
40 #include "CLHEP/Units/GlobalPhysicalConstants.h"
41 #include "TMath.h"
42 #include "TFile.h"
43 #include "TH1F.h"
44 #include "TH1I.h"
45 #include "TH2F.h"
46 #include "TProfile.h"
47 #include "TTree.h"
48 #include <iostream>
49 
50 using namespace reco;
51 
52 GsfElectronMCFakeAnalyzer::GsfElectronMCFakeAnalyzer(const edm::ParameterSet& conf) :
53  beamSpot_(conf.getParameter<edm::InputTag>("beamSpot"))
54 {
55 
56  outputFile_ = conf.getParameter<std::string>("outputFile");
57  histfile_ = new TFile(outputFile_.c_str(),"RECREATE");
58  electronCollection_=conf.getParameter<edm::InputTag>("electronCollection");
59  matchingObjectCollection_ = conf.getParameter<edm::InputTag>("matchingObjectCollection");
60  readAOD_ = conf.getParameter<bool>("readAOD");
61  maxPt_ = conf.getParameter<double>("MaxPt");
62  maxAbsEta_ = conf.getParameter<double>("MaxAbsEta");
63  deltaR_ = conf.getParameter<double>("DeltaR");
64 
65  edm::ParameterSet pset =
66  conf.getParameter<edm::ParameterSet>("HistosConfigurationFake") ;
67 
68  etamin=pset.getParameter<double>("Etamin");
69  etamax=pset.getParameter<double>("Etamax");
70  phimin=pset.getParameter<double>("Phimin");
71  phimax=pset.getParameter<double>("Phimax");
72  ptmax=pset.getParameter<double>("Ptmax");
73  pmax=pset.getParameter<double>("Pmax");
74  eopmax=pset.getParameter<double>("Eopmax");
75  eopmaxsht=pset.getParameter<double>("Eopmaxsht");
76  detamin=pset.getParameter<double>("Detamin");
77  detamax=pset.getParameter<double>("Detamax");
78  dphimin=pset.getParameter<double>("Dphimin");
79  dphimax=pset.getParameter<double>("Dphimax");
80  detamatchmin=pset.getParameter<double>("Detamatchmin");
81  detamatchmax=pset.getParameter<double>("Detamatchmax");
82  dphimatchmin=pset.getParameter<double>("Dphimatchmin");
83  dphimatchmax=pset.getParameter<double>("Dphimatchmax");
84  fhitsmax=pset.getParameter<double>("Fhitsmax");
85  lhitsmax=pset.getParameter<double>("Lhitsmax");
86  nbineta=pset.getParameter<int>("Nbineta");
87  nbineta2D=pset.getParameter<int>("Nbineta2D");
88  nbinp=pset.getParameter<int>("Nbinp");
89  nbinpt=pset.getParameter<int>("Nbinpt");
90  nbinp2D=pset.getParameter<int>("Nbinp2D");
91  nbinpt2D=pset.getParameter<int>("Nbinpt2D");
92  nbinpteff=pset.getParameter<int>("Nbinpteff");
93  nbinphi=pset.getParameter<int>("Nbinphi");
94  nbinphi2D=pset.getParameter<int>("Nbinphi2D");
95  nbineop=pset.getParameter<int>("Nbineop");
96  nbineop2D=pset.getParameter<int>("Nbineop2D");
97  nbinfhits=pset.getParameter<int>("Nbinfhits");
98  nbinlhits=pset.getParameter<int>("Nbinlhits");
99  nbinxyz=pset.getParameter<int>("Nbinxyz");
100  nbindeta=pset.getParameter<int>("Nbindeta");
101  nbindphi=pset.getParameter<int>("Nbindphi");
102  nbindetamatch=pset.getParameter<int>("Nbindetamatch");
103  nbindphimatch=pset.getParameter<int>("Nbindphimatch");
104  nbindetamatch2D=pset.getParameter<int>("Nbindetamatch2D");
105  nbindphimatch2D=pset.getParameter<int>("Nbindphimatch2D");
106  nbinmee= pset.getParameter<int>("Nbinmee");
107  meemin=pset.getParameter<double>("Meemin");
108  meemax=pset.getParameter<double>("Meemax");
109  nbinhoe= pset.getParameter<int>("Nbinhoe");
110  hoemin=pset.getParameter<double>("Hoemin");
111  hoemax=pset.getParameter<double>("Hoemax");
112 
113 }
114 
115 GsfElectronMCFakeAnalyzer::~GsfElectronMCFakeAnalyzer()
116 {
117 
118  // do anything here that needs to be done at desctruction time
119  // (e.g. close files, deallocate resources etc.)
120  histfile_->Write();
121  histfile_->Close();
122 }
123 
124 void GsfElectronMCFakeAnalyzer::beginJob(){
125 
126  histfile_->cd();
127 
128 
129  // matching object
130  std::string::size_type locJet = matchingObjectCollection_.label().find("iterativeCone5GenJets",0) ;
131  std::string type_;
132  if ( locJet != std::string::npos ) {
133  std::cout << "Matching objects are GenJets " << std::endl;
134  type_ = "GenJet";
135  } else {
136  std::cout << "Didn't recognize input matching objects!! " << std::endl;
137  }
138 
139  std::string htitle, hlabel;
140  hlabel="h_"+type_+"Num"; htitle="# "+type_+"s";
141  h_matchingObjectNum = new TH1F( hlabel.c_str(), htitle.c_str(), nbinfhits,0.,fhitsmax );
142 
143  // rec event
144  histNum_= new TH1F("h_recEleNum","# rec electrons",20, 0.,20.);
145 
146  // matching object distributions
147  hlabel="h_matchingObject_eta"; htitle=type_+" #eta";
148  h_matchingObjectEta = new TH1F( hlabel.c_str(), htitle.c_str(), nbineta,etamin,etamax);
149  hlabel="h_matchingObject_abseta"; htitle=type_+" |#eta|";
150  h_matchingObjectAbsEta = new TH1F( hlabel.c_str(), htitle.c_str(), nbineta/2,0.,etamax);
151  hlabel="h_matchingObject_P"; htitle=type_+" p";
152  h_matchingObjectP = new TH1F( hlabel.c_str(), htitle.c_str(), nbinp,0.,pmax);
153  hlabel="h_matchingObject_Pt"; htitle=type_+" pt";
154  h_matchingObjectPt = new TH1F( hlabel.c_str(),htitle.c_str(), nbinpteff,5.,ptmax);
155  hlabel="h_matchingObject_phi"; htitle=type_+" phi";
156  h_matchingObjectPhi = new TH1F( hlabel.c_str(), htitle.c_str(), nbinphi,phimin,phimax);
157  hlabel="h_matchingObject_z"; htitle=type_+" z";
158  h_matchingObjectZ = new TH1F( hlabel.c_str(), htitle.c_str(), nbinxyz, -25, 25 );
159 
160  // all electrons
161  h_ele_EoverP_all = new TH1F( "h_ele_EoverP_all", "ele E/P_{vertex}, all reco electrons", nbineop,0.,eopmax);
162  h_ele_EseedOP_all = new TH1F( "h_ele_EseedOP_all", "ele E_{seed}/P_{vertex}, all reco electrons", nbineop,0.,eopmax);
163  h_ele_EoPout_all = new TH1F( "h_ele_EoPout_all", "ele E_{seed}/P_{out}, all reco electrons", nbineop,0.,eopmax);
164  h_ele_EeleOPout_all = new TH1F( "h_ele_EeleOPout_all", "ele E_{ele}/P_{out}, all reco electrons", nbineop,0.,eopmax);
165  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);
166  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);
167  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);
168  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);
169  h_ele_TIP_all = new TH1F( "h_ele_TIP_all", "ele vertex transverse radius, all reco electrons", 100,0.,0.2);
170  h_ele_HoE_all = new TH1F("h_ele_HoE_all", "ele hadronic energy / em energy, all reco electrons", nbinhoe, hoemin, hoemax) ;
171  h_ele_vertexEta_all = new TH1F( "h_ele_vertexEta_all", "ele eta, all reco electrons", nbineta,etamin,etamax);
172  h_ele_vertexEta_all->Sumw2();
173  h_ele_vertexPt_all = new TH1F( "h_ele_vertexPt_all", "ele p_{T}, all reco electrons", nbinpteff,5.,ptmax);
174  h_ele_vertexPt_all->Sumw2();
175  h_ele_mee_all = new TH1F( "h_ele_mee_all", "ele pairs invariant mass, all reco electrons", nbinmee, meemin, meemax );
176  h_ele_mee_os = new TH1F( "h_ele_mee_os", "ele pairs invariant mass, opp. sign", nbinmee, meemin, meemax );
177 
178  // duplicates
179  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. );
180  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. );
181 
182  // matched electrons
183 
184  h_ele_matchingObjectPt_matched = new TH1F( "h_ele_matchingObjectPt_matched", "Efficiency vs matching jet p_{T}", nbinpteff,5.,ptmax);
185  h_ele_matchingObjectPt_matched->Sumw2();
186  h_ele_matchingObjectAbsEta_matched = new TH1F( "h_ele_matchingObjectAbsEta_matched", "Efficiency vs matching jet |#eta|", nbineta/2,0.,2.5);
187  h_ele_matchingObjectAbsEta_matched->Sumw2();
188  h_ele_matchingObjectEta_matched = new TH1F( "h_ele_matchingObjectEta_matched", "Efficiency vs matching jet #eta", nbineta,etamin,etamax);
189  h_ele_matchingObjectEta_matched ->Sumw2();
190  h_ele_matchingObjectPhi_matched = new TH1F( "h_ele_matchingObjectPhi_matched", "Efficiency vs matching jet phi", nbinphi,phimin,phimax);
191  h_ele_matchingObjectPhi_matched->Sumw2();
192  h_ele_matchingObjectZ_matched = new TH1F( "h_ele_matchingObjectZ_matched", "Efficiency vs matching jet z", nbinxyz,-25,25);
193  h_ele_matchingObjectZ_matched->Sumw2();
194 
195  h_ele_charge = new TH1F( "h_ele_charge", "ele charge", 5,-2.,2.);
196  h_ele_chargeVsEta = new TH2F( "h_ele_chargeVsEta", "ele charge vs eta", nbineta2D,etamin,etamax,5,-2.,2.);
197  h_ele_chargeVsPhi = new TH2F( "h_ele_chargeVsPhi", "ele charge vs phi", nbinphi2D,phimin,phimax,5,-2.,2.);
198  h_ele_chargeVsPt = new TH2F( "h_ele_chargeVsPt", "ele charge vs pt", nbinpt,0.,100.,5,-2.,2.);
199  h_ele_vertexP = new TH1F( "h_ele_vertexP", "ele momentum", nbinp,0.,pmax);
200  h_ele_vertexPt = new TH1F( "h_ele_vertexPt", "ele transverse momentum", nbinpt,0.,ptmax);
201  h_ele_vertexPtVsEta = new TH2F( "h_ele_vertexPtVsEta", "ele transverse momentum vs eta",nbineta2D,etamin,etamax,nbinpt2D,0.,ptmax);
202  h_ele_vertexPtVsPhi = new TH2F( "h_ele_vertexPtVsPhi", "ele transverse momentum vs phi",nbinphi2D,phimin,phimax,nbinpt2D,0.,ptmax);
203  h_ele_vertexEta = new TH1F( "h_ele_vertexEta", "ele momentum eta", nbineta,etamin,etamax);
204  h_ele_vertexEtaVsPhi = new TH2F( "h_ele_vertexEtaVsPhi", "ele momentum eta vs phi",nbineta2D,etamin,etamax,nbinphi2D,phimin,phimax );
205  h_ele_vertexPhi = new TH1F( "h_ele_vertexPhi", "ele momentum #phi", nbinphi,phimin,phimax);
206  h_ele_vertexX = new TH1F( "h_ele_vertexX", "ele vertex x", nbinxyz,-0.1,0.1 );
207  h_ele_vertexY = new TH1F( "h_ele_vertexY", "ele vertex y", nbinxyz,-0.1,0.1 );
208  h_ele_vertexZ = new TH1F( "h_ele_vertexZ", "ele vertex z", nbinxyz,-25, 25 );
209  h_ele_vertexTIP = new TH1F( "h_ele_vertexTIP", "ele transverse impact parameter (wrt bs)", 90,0.,0.15);
210  h_ele_vertexTIPVsEta = new TH2F( "h_ele_vertexTIPVsEta", "ele transverse impact parameter (wrt bs) vs eta", nbineta2D,etamin,etamax,45,0.,0.15);
211  h_ele_vertexTIPVsPhi = new TH2F( "h_ele_vertexTIPVsPhi", "ele transverse impact parameter (wrt bs) vs phi", nbinphi2D,phimin,phimax,45,0.,0.15);
212  h_ele_vertexTIPVsPt = new TH2F( "h_ele_vertexTIPVsPt", "ele transverse impact parameter (wrt bs) vs Pt", nbinpt2D,0.,ptmax,45,0.,0.15);
213 
214  h_ele_PoPmatchingObject = new TH1F( "h_ele_PoPmatchingObject", "ele momentum / matching jet momentum", 75,0.,1.5);
215  h_ele_PoPmatchingObjectVsEta = new TH2F( "h_ele_PoPmatchingObjectVsEta", "ele momentum / matching jet momentum vs eta", nbineta2D,etamin,etamax,50,0.,1.5);
216  h_ele_PoPmatchingObjectVsPhi = new TH2F( "h_ele_PoPmatchingObjectVsPhi", "ele momentum / matching jet momentum vs phi", nbinphi2D,phimin,phimax,50,0.,1.5);
217  h_ele_PoPmatchingObjectVsPt = new TH2F( "h_ele_PoPmatchingObjectVsPt", "ele momentum / matching jet momentum vs eta", nbinpt2D,0.,ptmax,50,0.,1.5);
218  h_ele_PoPmatchingObject_barrel = new TH1F( "h_ele_PoPmatchingObject_barrel", "ele momentum / matching jet momentum, barrel",75,0.,1.5);
219  h_ele_PoPmatchingObject_endcaps = new TH1F( "h_ele_PoPmatchingObject_endcaps", "ele momentum / matching jet momentum, endcaps",75,0.,1.5);
220 
221  h_ele_EtaMnEtamatchingObject = new TH1F( "h_ele_EtaMnEtamatchingObject", "ele momentum eta - matching jet eta",nbindeta,detamin,detamax);
222  h_ele_EtaMnEtamatchingObjectVsEta = new TH2F( "h_ele_EtaMnEtamatchingObjectVsEta", "ele momentum eta - matching jet eta vs eta",nbineta2D,etamin,etamax,nbindeta/2,detamin,detamax);
223  h_ele_EtaMnEtamatchingObjectVsPhi = new TH2F( "h_ele_EtaMnEtamatchingObjectVsPhi", "ele momentum eta - matching jet eta vs phi",nbinphi2D,phimin,phimax,nbindeta/2,detamin,detamax);
224  h_ele_EtaMnEtamatchingObjectVsPt = new TH2F( "h_ele_EtaMnEtamatchingObjectVsPt", "ele momentum eta - matching jet eta vs pt",nbinpt,0.,ptmax,nbindeta/2,detamin,detamax);
225  h_ele_PhiMnPhimatchingObject = new TH1F( "h_ele_PhiMnPhimatchingObject", "ele momentum phi - matching jet phi",nbindphi,dphimin,dphimax);
226  h_ele_PhiMnPhimatchingObject2 = new TH1F( "h_ele_PhiMnPhimatchingObject2", "ele momentum phi - matching jet phi",nbindphimatch2D,dphimatchmin,dphimatchmax);
227  h_ele_PhiMnPhimatchingObjectVsEta = new TH2F( "h_ele_PhiMnPhimatchingObjectVsEta", "ele momentum phi - matching SC phi vs eta",nbineta2D,etamin,etamax,nbindphi/2,dphimin,dphimax);
228  h_ele_PhiMnPhimatchingObjectVsPhi = new TH2F( "h_ele_PhiMnPhimatchingObjectVsPhi", "ele momentum phi - matching SC phi vs phi",nbinphi2D,phimin,phimax,nbindphi/2,dphimin,dphimax);
229  h_ele_PhiMnPhimatchingObjectVsPt = new TH2F( "h_ele_PhiMnPhimatchingObjectVsPt", "ele momentum phi - matching SC phi vs pt",nbinpt2D,0.,ptmax,nbindphi/2,dphimin,dphimax);
230 
231  // matched electron, superclusters
232  histSclEn_ = new TH1F("h_scl_energy","ele supercluster energy",nbinp,0.,pmax);
233  histSclEoEmatchingObject_barrel = new TH1F("h_scl_EoEmatchingObject_barrel","ele supercluster energy / matching jet energy, barrel",50,0.2,1.2);
234  histSclEoEmatchingObject_endcaps = new TH1F("h_scl_EoEmatchingObject_endcaps","ele supercluster energy / matching jet energy, endcaps",50,0.2,1.2);
235  histSclEt_ = new TH1F("h_scl_et","ele supercluster transverse energy",nbinpt,0.,ptmax);
236  histSclEtVsEta_ = new TH2F("h_scl_etVsEta","ele supercluster transverse energy vs eta",nbineta2D,etamin,etamax,nbinpt,0.,ptmax);
237  histSclEtVsPhi_ = new TH2F("h_scl_etVsPhi","ele supercluster transverse energy vs phi",nbinphi2D,phimin,phimax,nbinpt,0.,ptmax);
238  histSclEtaVsPhi_ = new TH2F("h_scl_etaVsPhi","ele supercluster eta vs phi",nbinphi2D,phimin,phimax,nbineta2D,etamin,etamax);
239  histSclEta_ = new TH1F("h_scl_eta","ele supercluster eta",nbineta,etamin,etamax);
240  histSclPhi_ = new TH1F("h_scl_phi","ele supercluster phi",nbinphi,phimin,phimax);
241 
242  histSclSigIEtaIEta_ = new TH1F("h_scl_sigietaieta","ele supercluster sigma ieta ieta",100,0.,0.05);
243  histSclSigIEtaIEta_barrel_ = new TH1F("h_scl_sigietaieta_barrel","ele supercluster sigma ieta ieta, barrel",100,0.,0.05);
244  histSclSigIEtaIEta_endcaps_ = new TH1F("h_scl_sigietaieta_endcaps","ele supercluster sigma ieta ieta, endcaps",100,0.,0.05);
245  histSclE1x5_ = new TH1F("h_scl_E1x5","ele supercluster energy in 1x5",nbinp,0.,pmax);
246  histSclE1x5_barrel_ = new TH1F("h_scl_E1x5_barrel","ele supercluster energy in 1x5 barrel",nbinp,0.,pmax);
247  histSclE1x5_endcaps_ = new TH1F("h_scl_E1x5_endcaps","ele supercluster energy in 1x5 endcaps",nbinp,0.,pmax);
248  histSclE2x5max_ = new TH1F("h_scl_E2x5max","ele supercluster energy in 2x5 max",nbinp,0.,pmax);
249  histSclE2x5max_barrel_ = new TH1F("h_scl_E2x5max_barrel","ele supercluster energy in 2x5 max barrel",nbinp,0.,pmax);
250  histSclE2x5max_endcaps_ = new TH1F("h_scl_E2x5max_endcaps","ele supercluster energy in 2x5 max endcaps",nbinp,0.,pmax);
251  histSclE5x5_ = new TH1F("h_scl_E5x5","ele supercluster energy in 5x5",nbinp,0.,pmax);
252  histSclE5x5_barrel_ = new TH1F("h_scl_E5x5_barrel","ele supercluster energy in 5x5 barrel",nbinp,0.,pmax);
253  histSclE5x5_endcaps_ = new TH1F("h_scl_E5x5_endcaps","ele supercluster energy in 5x5 endcaps",nbinp,0.,pmax);
254 
255  // matched electron, gsf tracks
256  h_ele_ambiguousTracks = new TH1F( "h_ele_ambiguousTracks", "ele # ambiguous tracks", 5,0.,5.);
257  h_ele_ambiguousTracksVsEta = new TH2F( "h_ele_ambiguousTracksVsEta","ele # ambiguous tracks vs eta", nbineta2D,etamin,etamax,5,0.,5.);
258  h_ele_ambiguousTracksVsPhi = new TH2F( "h_ele_ambiguousTracksVsPhi", "ele # ambiguous tracks vs phi", nbinphi2D,phimin,phimax,5,0.,5.);
259  h_ele_ambiguousTracksVsPt = new TH2F( "h_ele_ambiguousTracksVsPt", "ele # ambiguous tracks vs pt", nbinpt2D,0.,ptmax,5,0.,5.);
260  h_ele_foundHits = new TH1F( "h_ele_foundHits", "ele track # found hits", nbinfhits,0.,fhitsmax);
261  h_ele_foundHitsVsEta = new TH2F( "h_ele_foundHitsVsEta", "ele track # found hits vs eta", nbineta2D,etamin,etamax,nbinfhits,0.,fhitsmax);
262  h_ele_foundHitsVsPhi = new TH2F( "h_ele_foundHitsVsPhi", "ele track # found hits vs phi", nbinphi2D,phimin,phimax,nbinfhits,0.,fhitsmax);
263  h_ele_foundHitsVsPt = new TH2F( "h_ele_foundHitsVsPt", "ele track # found hits vs pt", nbinpt2D,0.,ptmax,nbinfhits,0.,fhitsmax);
264  h_ele_lostHits = new TH1F( "h_ele_lostHits", "ele track # lost hits", 5,0.,5.);
265  h_ele_lostHitsVsEta = new TH2F( "h_ele_lostHitsVsEta", "ele track # lost hits vs eta", nbineta2D,etamin,etamax,nbinlhits,0.,lhitsmax);
266  h_ele_lostHitsVsPhi = new TH2F( "h_ele_lostHitsVsPhi", "ele track # lost hits vs eta", nbinphi2D,phimin,phimax,nbinlhits,0.,lhitsmax);
267  h_ele_lostHitsVsPt = new TH2F( "h_ele_lostHitsVsPt", "ele track # lost hits vs eta", nbinpt2D,0.,ptmax,nbinlhits,0.,lhitsmax);
268  h_ele_chi2 = new TH1F( "h_ele_chi2", "ele track #chi^{2}", 100,0.,15.);
269  h_ele_chi2VsEta = new TH2F( "h_ele_chi2VsEta", "ele track #chi^{2} vs eta", nbineta2D,etamin,etamax,50,0.,15.);
270  h_ele_chi2VsPhi = new TH2F( "h_ele_chi2VsPhi", "ele track #chi^{2} vs phi", nbinphi2D,phimin,phimax,50,0.,15.);
271  h_ele_chi2VsPt = new TH2F( "h_ele_chi2VsPt", "ele track #chi^{2} vs pt", nbinpt2D,0.,ptmax,50,0.,15.);
272  h_ele_PinMnPout = new TH1F( "h_ele_PinMnPout", "ele track inner p - outer p, mean" ,nbinp,0.,200.);
273  h_ele_PinMnPout_mode = new TH1F( "h_ele_PinMnPout_mode", "ele track inner p - outer p, mode" ,nbinp,0.,100.);
274  h_ele_PinMnPoutVsEta_mode = new TH2F( "h_ele_PinMnPoutVsEta_mode", "ele track inner p - outer p vs eta, mode" ,nbineta2D, etamin,etamax,nbinp2D,0.,100.);
275  h_ele_PinMnPoutVsPhi_mode = new TH2F( "h_ele_PinMnPoutVsPhi_mode", "ele track inner p - outer p vs phi, mode" ,nbinphi2D, phimin,phimax,nbinp2D,0.,100.);
276  h_ele_PinMnPoutVsPt_mode = new TH2F( "h_ele_PinMnPoutVsPt_mode", "ele track inner p - outer p vs pt, mode" ,nbinpt2D, 0.,ptmax,nbinp2D,0.,100.);
277  h_ele_PinMnPoutVsE_mode = new TH2F( "h_ele_PinMnPoutVsE_mode", "ele track inner p - outer p vs E, mode" ,nbinp2D, 0.,200.,nbinp2D,0.,100.);
278  h_ele_PinMnPoutVsChi2_mode = new TH2F( "h_ele_PinMnPoutVsChi2_mode", "ele track inner p - outer p vs track chi2, mode" ,50, 0.,20.,nbinp2D,0.,100.);
279  h_ele_outerP = new TH1F( "h_ele_outerP", "ele track outer p, mean", nbinp,0.,pmax);
280  h_ele_outerP_mode = new TH1F( "h_ele_outerP_mode", "ele track outer p, mode", nbinp,0.,pmax);
281  h_ele_outerPVsEta_mode = new TH2F( "h_ele_outerPVsEta_mode", "ele track outer p vs eta mode", nbineta2D,etamin,etamax,50,0.,pmax);
282  h_ele_outerPt = new TH1F( "h_ele_outerPt", "ele track outer p_{T}, mean", nbinpt,0.,ptmax);
283  h_ele_outerPt_mode = new TH1F( "h_ele_outerPt_mode", "ele track outer p_{T}, mode", nbinpt,0.,ptmax);
284  h_ele_outerPtVsEta_mode = new TH2F( "h_ele_outerPtVsEta_mode", "ele track outer p_{T} vs eta, mode", nbineta2D,etamin,etamax,nbinpt2D,0.,ptmax);
285  h_ele_outerPtVsPhi_mode = new TH2F( "h_ele_outerPtVsPhi_mode", "ele track outer p_{T} vs phi, mode", nbinphi2D,phimin,phimax,nbinpt2D,0.,ptmax);
286  h_ele_outerPtVsPt_mode = new TH2F( "h_ele_outerPtVsPt_mode", "ele track outer p_{T} vs pt, mode", nbinpt2D,0.,100.,nbinpt2D,0.,ptmax);
287 
288  // matched electrons,cluster-track matching
289  h_ele_EoP = new TH1F( "h_ele_EoP", "ele E/P_{vertex}", nbineop,0.,eopmax); h_ele_EoP->Sumw2();
290  h_ele_EoP_eg = new TH1F( "h_ele_EoP_eg", "ele E/P_{vertex}, ecal driven", nbineop,0.,eopmax);
291  h_ele_EoP_eg->Sumw2();
292  h_ele_EoP_barrel = new TH1F( "h_ele_EoP_barrel", "ele E/P_{vertex} barrel", nbineop,0.,eopmax);
293  h_ele_EoP_barrel->Sumw2();
294  h_ele_EoP_eg_barrel = new TH1F( "h_ele_EoP_eg_barrel", "ele E/P_{vertex}, ecal driven barrel", nbineop,0.,eopmax);
295  h_ele_EoP_eg_barrel->Sumw2();
296  h_ele_EoP_endcaps = new TH1F( "h_ele_EoP_endcaps", "ele E/P_{vertex} endcaps", nbineop,0.,eopmax);
297  h_ele_EoP_endcaps->Sumw2();
298  h_ele_EoP_eg_endcaps = new TH1F( "h_ele_EoP_eg_endcaps", "ele E/P_{vertex}, ecal driven endcaps", nbineop,0.,eopmax);
299  h_ele_EoP_eg_endcaps->Sumw2();
300  h_ele_EoPVsEta = new TH2F( "h_ele_EoPVsEta", "ele E/P_{vertex} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
301  h_ele_EoPVsPhi = new TH2F( "h_ele_EoPVsPhi", "ele E/P_{vertex} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
302  h_ele_EoPVsE = new TH2F( "h_ele_EoPVsE", "ele E/P_{vertex} vs E", 50,0.,pmax ,50,0.,5.);
303  h_ele_EseedOP = new TH1F( "h_ele_EseedOP", "ele E_{seed}/P_{vertex}", nbineop,0.,eopmax);
304  h_ele_EseedOP->Sumw2();
305  h_ele_EseedOP_eg = new TH1F( "h_ele_EseedOP_eg", "ele E_{seed}/P_{vertex}, ecal driven", nbineop,0.,eopmax);
306  h_ele_EseedOP_eg->Sumw2();
307  h_ele_EseedOP_barrel = new TH1F( "h_ele_EseedOP_barrel", "ele E_{seed}/P_{vertex} barrel", nbineop,0.,eopmax);
308  h_ele_EseedOP_barrel->Sumw2();
309  h_ele_EseedOP_eg_barrel = new TH1F( "h_ele_EseedOP_eg_barrel", "ele E_{seed}/P_{vertex}, ecal driven barrel", nbineop,0.,eopmax);
310  h_ele_EseedOP_eg_barrel->Sumw2();
311  h_ele_EseedOP_endcaps = new TH1F( "h_ele_EseedOP_endcaps", "ele E_{seed}/P_{vertex} endcaps", nbineop,0.,eopmax);
312  h_ele_EseedOP_endcaps->Sumw2();
313  h_ele_EseedOP_eg_endcaps = new TH1F( "h_ele_EseedOP_eg_endcaps", "ele E_{seed}/P_{vertex}, ecal driven, endcaps", nbineop,0.,eopmax);
314  h_ele_EseedOP_eg_endcaps->Sumw2();
315  h_ele_EseedOPVsEta = new TH2F( "h_ele_EseedOPVsEta", "ele E_{seed}/P_{vertex} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
316  h_ele_EseedOPVsPhi = new TH2F( "h_ele_EseedOPVsPhi", "ele E_{seed}/P_{vertex} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
317  h_ele_EseedOPVsE = new TH2F( "h_ele_EseedOPVsE", "ele E_{seed}/P_{vertex} vs E", 50,0.,pmax ,50,0.,5.);
318  h_ele_EoPout = new TH1F( "h_ele_EoPout", "ele E_{seed}/P_{out}", nbineop,0.,eopmax);
319  h_ele_EoPout->Sumw2();
320  h_ele_EoPout_eg = new TH1F( "h_ele_EoPout_eg", "ele E_{seed}/P_{out}, ecal driven", nbineop,0.,eopmax);
321  h_ele_EoPout_eg->Sumw2();
322  h_ele_EoPout_barrel = new TH1F( "h_ele_EoPout_barrel", "ele E_{seed}/P_{out} barrel", nbineop,0.,eopmax);
323  h_ele_EoPout_barrel->Sumw2();
324  h_ele_EoPout_eg_barrel = new TH1F( "h_ele_EoPout_eg_barrel", "ele E_{seed}/P_{out}, ecal driven, barrel", nbineop,0.,eopmax);
325  h_ele_EoPout_eg_barrel->Sumw2();
326  h_ele_EoPout_endcaps = new TH1F( "h_ele_EoPout_endcaps", "ele E_{seed}/P_{out} endcaps", nbineop,0.,eopmax);
327  h_ele_EoPout_endcaps->Sumw2();
328  h_ele_EoPout_eg_endcaps = new TH1F( "h_ele_EoPout_eg_endcaps", "ele E_{seed}/P_{out}, ecal driven, endcaps", nbineop,0.,eopmax);
329  h_ele_EoPout_eg_endcaps->Sumw2();
330  h_ele_EoPoutVsEta = new TH2F( "h_ele_EoPoutVsEta", "ele E/P_{out} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
331  h_ele_EoPoutVsPhi = new TH2F( "h_ele_EoPoutVsPhi", "ele E/P_{out} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
332  h_ele_EoPoutVsE = new TH2F( "h_ele_EoPoutVsE", "ele E/P_{out} vs E", nbinp2D,0.,pmax,nbineop2D,0.,eopmaxsht);
333  h_ele_EeleOPout = new TH1F( "h_ele_EeleOPout", "ele E_{ele}/P_{out}", nbineop,0.,eopmax);
334  h_ele_EeleOPout->Sumw2();
335  h_ele_EeleOPout_eg = new TH1F( "h_ele_EeleOPout_eg", "ele E_{ele}/P_{out}, ecal driven", nbineop,0.,eopmax);
336  h_ele_EeleOPout_eg->Sumw2();
337  h_ele_EeleOPout_barrel = new TH1F( "h_ele_EeleOPout_barrel", "ele E_{ele}/P_{out} barrel", nbineop,0.,eopmax);
338  h_ele_EeleOPout_barrel->Sumw2();
339  h_ele_EeleOPout_eg_barrel = new TH1F( "h_ele_EeleOPout_eg_barrel", "ele E_{ele}/P_{out}, ecal driven, barrel", nbineop,0.,eopmax);
340  h_ele_EeleOPout_eg_barrel->Sumw2();
341  h_ele_EeleOPout_endcaps = new TH1F( "h_ele_EeleOPout_endcaps", "ele E_{ele}/P_{out} endcaps", nbineop,0.,eopmax);
342  h_ele_EeleOPout_endcaps->Sumw2();
343  h_ele_EeleOPout_eg_endcaps = new TH1F( "h_ele_EeleOPout_eg_endcaps", "ele E_{ele}/P_{out}, ecal driven, endcaps", nbineop,0.,eopmax);
344  h_ele_EeleOPout_eg_endcaps->Sumw2();
345  h_ele_EeleOPoutVsEta = new TH2F( "h_ele_EeleOPoutVsEta", "ele E_{ele}/P_{out} vs eta", nbineta2D,etamin,etamax,nbineop2D,0.,eopmaxsht);
346  h_ele_EeleOPoutVsPhi = new TH2F( "h_ele_EeleOPoutVsPhi", "ele E_{ele}/P_{out} vs phi", nbinphi2D,phimin,phimax,nbineop2D,0.,eopmaxsht);
347  h_ele_EeleOPoutVsE = new TH2F( "h_ele_EeleOPoutVsE", "ele E_{ele}/P_{out} vs E", nbinp2D,0.,pmax,nbineop2D,0.,eopmaxsht);
348  h_ele_dEtaSc_propVtx = new TH1F( "h_ele_dEtaSc_propVtx", "ele #eta_{sc} - #eta_{tr}, prop from vertex", nbindetamatch,detamatchmin,detamatchmax);
349  h_ele_dEtaSc_propVtx->Sumw2();
350  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);
351  h_ele_dEtaSc_propVtx_eg->Sumw2();
352  h_ele_dEtaSc_propVtx_barrel = new TH1F( "h_ele_dEtaSc_propVtx_barrel", "ele #eta_{sc} - #eta_{tr}, prop from vertex, barrel", nbindetamatch,detamatchmin,detamatchmax);
353  h_ele_dEtaSc_propVtx_barrel->Sumw2();
354  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);
355  h_ele_dEtaSc_propVtx_eg_barrel->Sumw2();
356  h_ele_dEtaSc_propVtx_endcaps = new TH1F( "h_ele_dEtaSc_propVtx_endcaps", "ele #eta_{sc} - #eta_{tr}, prop from vertex, endcaps", nbindetamatch,detamatchmin,detamatchmax);
357  h_ele_dEtaSc_propVtx_endcaps->Sumw2();
358  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);
359  h_ele_dEtaSc_propVtx_eg_endcaps->Sumw2();
360  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);
361  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);
362  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);
363  h_ele_dPhiSc_propVtx = new TH1F( "h_ele_dPhiSc_propVtx", "ele #phi_{sc} - #phi_{tr}, prop from vertex", nbindphimatch,dphimatchmin,dphimatchmax);
364  h_ele_dPhiSc_propVtx->Sumw2();
365  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);
366  h_ele_dPhiSc_propVtx_eg->Sumw2();
367  h_ele_dPhiSc_propVtx_barrel = new TH1F( "h_ele_dPhiSc_propVtx_barrel", "ele #phi_{sc} - #phi_{tr}, prop from vertex, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
368  h_ele_dPhiSc_propVtx_barrel->Sumw2();
369  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);
370  h_ele_dPhiSc_propVtx_eg_barrel->Sumw2();
371  h_ele_dPhiSc_propVtx_endcaps = new TH1F( "h_ele_dPhiSc_propVtx_endcaps", "ele #phi_{sc} - #phi_{tr}, prop from vertex, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
372  h_ele_dPhiSc_propVtx_endcaps->Sumw2();
373  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);
374  h_ele_dPhiSc_propVtx_eg_endcaps->Sumw2();
375  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);
376  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);
377  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);
378  h_ele_dEtaCl_propOut = new TH1F( "h_ele_dEtaCl_propOut", "ele #eta_{cl} - #eta_{tr}, prop from outermost", nbindetamatch,detamatchmin,detamatchmax);
379  h_ele_dEtaCl_propOut->Sumw2();
380  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);
381  h_ele_dEtaCl_propOut_eg->Sumw2();
382  h_ele_dEtaCl_propOut_barrel = new TH1F( "h_ele_dEtaCl_propOut_barrel", "ele #eta_{cl} - #eta_{tr}, prop from outermost, barrel", nbindetamatch,detamatchmin,detamatchmax);
383  h_ele_dEtaCl_propOut_barrel->Sumw2();
384  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);
385  h_ele_dEtaCl_propOut_eg_barrel->Sumw2();
386  h_ele_dEtaCl_propOut_endcaps = new TH1F( "h_ele_dEtaCl_propOut_endcaps", "ele #eta_{cl} - #eta_{tr}, prop from outermost, endcaps", nbindetamatch,detamatchmin,detamatchmax);
387  h_ele_dEtaCl_propOut_endcaps->Sumw2();
388  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);
389  h_ele_dEtaCl_propOut_eg_endcaps->Sumw2();
390  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);
391  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);
392  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);
393  h_ele_dPhiCl_propOut = new TH1F( "h_ele_dPhiCl_propOut", "ele #phi_{cl} - #phi_{tr}, prop from outermost", nbindphimatch,dphimatchmin,dphimatchmax);
394  h_ele_dPhiCl_propOut->Sumw2();
395  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);
396  h_ele_dPhiCl_propOut_eg->Sumw2();
397  h_ele_dPhiCl_propOut_barrel = new TH1F( "h_ele_dPhiCl_propOut_barrel", "ele #phi_{cl} - #phi_{tr}, prop from outermost, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
398  h_ele_dPhiCl_propOut_barrel->Sumw2();
399  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);
400  h_ele_dPhiCl_propOut_eg_barrel->Sumw2();
401  h_ele_dPhiCl_propOut_endcaps = new TH1F( "h_ele_dPhiCl_propOut_endcaps", "ele #phi_{cl} - #phi_{tr}, prop from outermost, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
402  h_ele_dPhiCl_propOut_endcaps->Sumw2();
403  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);
404  h_ele_dPhiCl_propOut_eg_endcaps->Sumw2();
405  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);
406  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);
407  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);
408  h_ele_dEtaEleCl_propOut = new TH1F( "h_ele_dEtaEleCl_propOut", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost", nbindetamatch,detamatchmin,detamatchmax);
409  h_ele_dEtaEleCl_propOut->Sumw2();
410  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);
411  h_ele_dEtaEleCl_propOut_eg->Sumw2();
412  h_ele_dEtaEleCl_propOut_barrel = new TH1F( "h_ele_dEtaEleCl_propOut_barrel", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, barrel", nbindetamatch,detamatchmin,detamatchmax);
413  h_ele_dEtaEleCl_propOut_barrel->Sumw2();
414  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);
415  h_ele_dEtaEleCl_propOut_eg_barrel->Sumw2();
416  h_ele_dEtaEleCl_propOut_endcaps = new TH1F( "h_ele_dEtaEleCl_propOut_endcaps", "ele #eta_{EleCl} - #eta_{tr}, prop from outermost, endcaps", nbindetamatch,detamatchmin,detamatchmax);
417  h_ele_dEtaEleCl_propOut_endcaps->Sumw2();
418  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);
419  h_ele_dEtaEleCl_propOut_eg_endcaps->Sumw2();
420  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);
421  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);
422  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);
423  h_ele_dPhiEleCl_propOut = new TH1F( "h_ele_dPhiEleCl_propOut", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost", nbindphimatch,dphimatchmin,dphimatchmax);
424  h_ele_dPhiEleCl_propOut->Sumw2();
425  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);
426  h_ele_dPhiEleCl_propOut_eg->Sumw2();
427  h_ele_dPhiEleCl_propOut_barrel = new TH1F( "h_ele_dPhiEleCl_propOut_barrel", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, barrel", nbindphimatch,dphimatchmin,dphimatchmax);
428  h_ele_dPhiEleCl_propOut_barrel->Sumw2();
429  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);
430  h_ele_dPhiEleCl_propOut_eg_barrel->Sumw2();
431  h_ele_dPhiEleCl_propOut_endcaps = new TH1F( "h_ele_dPhiEleCl_propOut_endcaps", "ele #phi_{EleCl} - #phi_{tr}, prop from outermost, endcaps", nbindphimatch,dphimatchmin,dphimatchmax);
432  h_ele_dPhiEleCl_propOut_endcaps->Sumw2();
433  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);
434  h_ele_dPhiEleCl_propOut_eg_endcaps->Sumw2();
435  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);
436  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);
437  h_ele_dPhiEleClVsPt_propOut = new TH2F( "h_ele_dPhiSEleClsPt_propOut", "ele #phi_{EleCl} - #phi_{tr} vs pt, prop from out", nbinpt2D,0.,ptmax,nbindphimatch2D,dphimatchmin,dphimatchmax);
438 
439  h_ele_HoE = new TH1F("h_ele_HoE", "ele hadronic energy / em energy", nbinhoe, hoemin, hoemax) ;
440  h_ele_HoE->Sumw2();
441  h_ele_HoE_eg = new TH1F("h_ele_HoE_eg", "ele hadronic energy / em energy, ecal driven", nbinhoe, hoemin, hoemax) ;
442  h_ele_HoE_eg->Sumw2();
443  h_ele_HoE_barrel = new TH1F("h_ele_HoE_barrel", "ele hadronic energy / em energy, barrel", nbinhoe, hoemin, hoemax) ;
444  h_ele_HoE_barrel->Sumw2();
445  h_ele_HoE_eg_barrel = new TH1F("h_ele_HoE_eg_barrel", "ele hadronic energy / em energy, ecal driven, barrel", nbinhoe, hoemin, hoemax) ;
446  h_ele_HoE_eg_barrel->Sumw2();
447  h_ele_HoE_endcaps = new TH1F("h_ele_HoE_endcaps", "ele hadronic energy / em energy, endcaps", nbinhoe, hoemin, hoemax) ;
448  h_ele_HoE_endcaps->Sumw2();
449  h_ele_HoE_eg_endcaps = new TH1F("h_ele_HoE_eg_endcaps", "ele hadronic energy / em energy, ecal driven, endcaps", nbinhoe, hoemin, hoemax) ;
450  h_ele_HoE_eg_endcaps->Sumw2();
451  h_ele_HoE_fiducial = new TH1F("h_ele_HoE_fiducial", "ele hadronic energy / em energy",nbinhoe, hoemin, hoemax) ;
452  h_ele_HoEVsEta = new TH2F("h_ele_HoEVsEta", "ele hadronic energy / em energy vs eta", nbineta,etamin,etamax,nbinhoe, hoemin, hoemax) ;
453  h_ele_HoEVsPhi = new TH2F("h_ele_HoEVsPhi", "ele hadronic energy / em energy vs phi", nbinphi2D,phimin,phimax,nbinhoe, hoemin, hoemax) ;
454  h_ele_HoEVsE = new TH2F("h_ele_HoEVsE", "ele hadronic energy / em energy vs E", nbinp, 0.,300.,nbinhoe, hoemin, hoemax) ;
455 
456  h_ele_seed_dphi2_ = new TH1F("h_ele_seedDphi2", "ele seed dphi 2nd layer", 50,-0.003,+0.003) ;
457  h_ele_seed_dphi2VsEta_ = new TH2F("h_ele_seedDphi2VsEta", "ele seed dphi 2nd layer vs eta", nbineta2D,etamin,etamax,50,-0.003,+0.003) ;
458  h_ele_seed_dphi2VsPt_ = new TH2F("h_ele_seedDphi2VsPt", "ele seed dphi 2nd layer vs pt", nbinpt2D,0.,ptmax,50,-0.003,+0.003) ;
459  h_ele_seed_drz2_ = new TH1F("h_ele_seedDrz2", "ele seed dr/dz 2nd layer", 50,-0.03,+0.03) ;
460  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) ;
461  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) ;
462  h_ele_seed_subdet2_ = new TH1F("h_ele_seedSubdet2", "ele seed subdet 2nd layer", 10,0.,10.) ;
463 
464  // classes
465  h_ele_classes = new TH1F( "h_ele_classes", "electron classes", 20,0.0,20.);
466  h_ele_eta = new TH1F( "h_ele_eta", "ele electron eta", nbineta/2,0.0,etamax);
467  h_ele_eta_golden = new TH1F( "h_ele_eta_golden", "ele electron eta golden", nbineta/2,0.0,etamax);
468  h_ele_eta_bbrem = new TH1F( "h_ele_eta_bbrem", "ele electron eta bbrem", nbineta/2,0.0,etamax);
469  h_ele_eta_narrow = new TH1F( "h_ele_eta_narrow", "ele electron eta narrow", nbineta/2,0.0,etamax);
470  h_ele_eta_shower = new TH1F( "h_ele_eta_show", "el, electron eta showering", nbineta/2,0.0,etamax);
471  h_ele_PinVsPoutGolden_mode = new TH2F( "h_ele_PinVsPoutGolden_mode", "ele track inner p vs outer p, golden, mode" ,nbinp2D,0.,pmax,50,0.,pmax);
472  h_ele_PinVsPoutShowering_mode = new TH2F( "h_ele_PinVsPoutShowering_mode", "ele track inner p vs outer p vs eta, showering, mode" ,nbinp2D,0.,pmax,50,0.,pmax);
473  h_ele_PinVsPoutGolden_mean = new TH2F( "h_ele_PinVsPoutGolden_mean", "ele track inner p vs outer p, golden, mean" ,nbinp2D,0.,pmax,50,0.,pmax);
474  h_ele_PinVsPoutShowering_mean = new TH2F( "h_ele_PinVsPoutShowering_mean", "ele track inner p vs outer p, showering, mean" ,nbinp2D,0.,pmax,50,0.,pmax);
475  h_ele_PtinVsPtoutGolden_mode = new TH2F( "h_ele_PtinVsPtoutGolden_mode", "ele track inner pt vs outer pt, golden, mode" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
476  h_ele_PtinVsPtoutShowering_mode = new TH2F( "h_ele_PtinVsPtoutShowering_mode", "ele track inner pt vs outer pt, showering, mode" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
477  h_ele_PtinVsPtoutGolden_mean = new TH2F( "h_ele_PtinVsPtoutGolden_mean", "ele track inner pt vs outer pt, golden, mean" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
478  h_ele_PtinVsPtoutShowering_mean = new TH2F( "h_ele_PtinVsPtoutShowering_mean", "ele track inner pt vs outer pt, showering, mean" ,nbinpt2D,0.,ptmax,50,0.,ptmax);
479  histSclEoEmatchingObjectGolden_barrel = new TH1F("h_scl_EoEmatchingObject golden, barrel","ele SC energy over matching jet energy, golden, barrel",100,0.2,1.2);
480  histSclEoEmatchingObjectGolden_endcaps = new TH1F("h_scl_EoEmatchingObject golden, endcaps","ele SC energy over matching jet energy, golden, endcaps",100,0.2,1.2);
481  histSclEoEmatchingObjectShowering_barrel = new TH1F("h_scl_EoEmatchingObject showering, barrel","ele SC energy over matching jet energy, showering, barrel",100,0.2,1.2);
482  histSclEoEmatchingObjectShowering_endcaps = new TH1F("h_scl_EoEmatchingObject showering, endcaps","ele SC energy over matching jet energy, showering, endcaps",100,0.2,1.2);
483 
484  // isolation
485  h_ele_tkSumPt_dr03 = new TH1F("h_ele_tkSumPt_dr03","tk isolation sum, dR=0.3",100,0.0,20.);
486  h_ele_ecalRecHitSumEt_dr03= new TH1F("h_ele_ecalRecHitSumEt_dr03","ecal isolation sum, dR=0.3",100,0.0,20.);
487  h_ele_hcalDepth1TowerSumEt_dr03= new TH1F("h_ele_hcalDepth1TowerSumEt_dr03","hcal depth1 isolation sum, dR=0.3",100,0.0,20.);
488  h_ele_hcalDepth2TowerSumEt_dr03= new TH1F("h_ele_hcalDepth2TowerSumEt_dr03","hcal depth2 isolation sum, dR=0.3",100,0.0,20.);
489  h_ele_tkSumPt_dr04= new TH1F("h_ele_tkSumPt_dr04","hcal isolation sum",100,0.0,20.);
490  h_ele_ecalRecHitSumEt_dr04= new TH1F("h_ele_ecalRecHitSumEt_dr04","ecal isolation sum, dR=0.4",100,0.0,20.);
491  h_ele_hcalDepth1TowerSumEt_dr04= new TH1F("h_ele_hcalDepth1TowerSumEt_dr04","hcal depth1 isolation sum, dR=0.4",100,0.0,20.);
492  h_ele_hcalDepth2TowerSumEt_dr04= new TH1F("h_ele_hcalDepth2TowerSumEt_dr04","hcal depth2 isolation sum, dR=0.4",100,0.0,20.);
493 
494  // fbrem
495  h_ele_fbrem = new TH1F( "h_ele_fbrem","ele brem fraction, mode",100,0.,1.);
496  h_ele_fbremVsEta_mode = new TProfile( "h_ele_fbremvsEtamode","mean ele brem fraction vs eta, mode",nbineta2D,etamin,etamax,0.,1.);
497  h_ele_fbremVsEta_mean = new TProfile( "h_ele_fbremvsEtamean","mean ele brem fraction vs eta, mean",nbineta2D,etamin,etamax,0.,1.);
498 
499  // e/g et pflow electrons
500  h_ele_mva = new TH1F( "h_ele_mva","ele mva",100,-1.,1.);
501  h_ele_provenance = new TH1F( "h_ele_provenance","ele provenance",5,-2.,3.);
502 
503  // histos titles
504  h_matchingObjectNum -> GetXaxis()-> SetTitle("N_{GenJet}");
505  h_matchingObjectNum -> GetYaxis()-> SetTitle("Events");
506  h_matchingObjectEta -> GetXaxis()-> SetTitle("#eta_{GenJet}");
507  h_matchingObjectEta -> GetYaxis()-> SetTitle("Events");
508  h_matchingObjectP -> GetXaxis()-> SetTitle("E_{GenJet} (GeV)");
509  h_matchingObjectP -> GetYaxis()-> SetTitle("Events");
510  h_ele_foundHits -> GetXaxis()-> SetTitle("N_{hits}");
511  h_ele_foundHits -> GetYaxis()-> SetTitle("Events");
512  h_ele_ambiguousTracks -> GetXaxis()-> SetTitle("N_{ambiguous tracks}");
513  h_ele_ambiguousTracks -> GetYaxis()-> SetTitle("Events");
514  h_ele_lostHits -> GetXaxis()-> SetTitle("N_{lost hits}");
515  h_ele_lostHits -> GetYaxis()-> SetTitle("Events");
516  h_ele_chi2 -> GetXaxis()-> SetTitle("#Chi^{2}");
517  h_ele_chi2 -> GetYaxis()-> SetTitle("Events");
518  h_ele_charge -> GetXaxis()-> SetTitle("charge");
519  h_ele_charge -> GetYaxis()-> SetTitle("Events");
520  h_ele_vertexP -> GetXaxis()-> SetTitle("p_{vertex} (GeV/c)");
521  h_ele_vertexP -> GetYaxis()-> SetTitle("Events");
522  h_ele_vertexPt -> GetXaxis()-> SetTitle("p_{T vertex} (GeV/c)");
523  h_ele_vertexPt -> GetYaxis()-> SetTitle("Events");
524  h_ele_vertexEta -> GetXaxis()-> SetTitle("#eta");
525  h_ele_vertexEta -> GetYaxis()-> SetTitle("Events");
526  h_ele_vertexPhi -> GetXaxis()-> SetTitle("#phi (rad)");
527  h_ele_vertexPhi -> GetYaxis()-> SetTitle("Events");
528  h_ele_PoPmatchingObject -> GetXaxis()-> SetTitle("P/P_{gen jet}");
529  h_ele_PoPmatchingObject -> GetYaxis()-> SetTitle("Events");
530  h_ele_PoPmatchingObject_barrel -> GetXaxis()-> SetTitle("P/P_{gen jet}");
531  h_ele_PoPmatchingObject_barrel -> GetYaxis()-> SetTitle("Events");
532  h_ele_PoPmatchingObject_endcaps -> GetXaxis()-> SetTitle("P/P_{gen jet}");
533  h_ele_PoPmatchingObject_endcaps -> GetYaxis()-> SetTitle("Events");
534  histSclSigIEtaIEta_-> GetXaxis()-> SetTitle("#sigma_{i#eta i#eta}") ;
535  histSclSigIEtaIEta_-> GetYaxis()-> SetTitle("Events") ;
536  histSclSigIEtaIEta_barrel_-> GetXaxis()-> SetTitle("#sigma_{i#eta i#eta}") ;
537  histSclSigIEtaIEta_barrel_-> GetYaxis()-> SetTitle("Events") ;
538  histSclSigIEtaIEta_endcaps_-> GetXaxis()-> SetTitle("#sigma_{i#eta i#eta}") ;
539  histSclSigIEtaIEta_endcaps_-> GetYaxis()-> SetTitle("Events") ;
540  histSclE1x5_-> GetXaxis()-> SetTitle("E1x5 (GeV)") ;
541  histSclE1x5_-> GetYaxis()-> SetTitle("Events") ;
542  histSclE1x5_barrel_-> GetXaxis()-> SetTitle("E1x5 (GeV)") ;
543  histSclE1x5_barrel_-> GetYaxis()-> SetTitle("Events") ;
544  histSclE1x5_endcaps_-> GetXaxis()-> SetTitle("E1x5 (GeV)") ;
545  histSclE1x5_endcaps_-> GetYaxis()-> SetTitle("Events") ;
546  histSclE2x5max_-> GetXaxis()-> SetTitle("E2x5 (GeV)") ;
547  histSclE2x5max_-> GetYaxis()-> SetTitle("Events") ;
548  histSclE2x5max_barrel_-> GetXaxis()-> SetTitle("E2x5 (GeV)") ;
549  histSclE2x5max_barrel_-> GetYaxis()-> SetTitle("Events") ;
550  histSclE2x5max_endcaps_-> GetXaxis()-> SetTitle("E2x5 (GeV)") ;
551  histSclE2x5max_endcaps_-> GetYaxis()-> SetTitle("Events") ;
552  histSclE5x5_-> GetXaxis()-> SetTitle("E5x5 (GeV)") ;
553  histSclE5x5_-> GetYaxis()-> SetTitle("Events") ;
554  histSclE5x5_barrel_-> GetXaxis()-> SetTitle("E5x5 (GeV)") ;
555  histSclE5x5_barrel_-> GetYaxis()-> SetTitle("Events") ;
556  histSclE5x5_endcaps_-> GetXaxis()-> SetTitle("E5x5 (GeV)") ;
557  histSclE5x5_endcaps_-> GetYaxis()-> SetTitle("Events") ;
558  h_ele_EtaMnEtamatchingObject -> GetXaxis()-> SetTitle("#eta_{rec} - #eta_{jet}");
559  h_ele_EtaMnEtamatchingObject -> GetYaxis()-> SetTitle("Events");
560  h_ele_PhiMnPhimatchingObject -> GetXaxis()-> SetTitle("#phi_{rec} - #phi_{jet} (rad)");
561  h_ele_PhiMnPhimatchingObject -> GetYaxis()-> SetTitle("Events");
562  h_ele_PinMnPout -> GetXaxis()-> SetTitle("P_{vertex} - P_{out} (GeV/c)");
563  h_ele_PinMnPout -> GetYaxis()-> SetTitle("Events");
564  h_ele_PinMnPout_mode -> GetXaxis()-> SetTitle("P_{vertex} - P_{out}, mode (GeV/c)");
565  h_ele_PinMnPout_mode -> GetYaxis()-> SetTitle("Events");
566  h_ele_outerP -> GetXaxis()-> SetTitle("P_{out} (GeV/c)");
567  h_ele_outerP -> GetYaxis()-> SetTitle("Events");
568  h_ele_outerP_mode -> GetXaxis()-> SetTitle("P_{out} (GeV/c)");
569  h_ele_outerP_mode -> GetYaxis()-> SetTitle("Events");
570  h_ele_outerPt -> GetXaxis()-> SetTitle("P_{T out} (GeV/c)");
571  h_ele_outerPt -> GetYaxis()-> SetTitle("Events");
572  h_ele_outerPt_mode -> GetXaxis()-> SetTitle("P_{T out} (GeV/c)");
573  h_ele_outerPt_mode -> GetYaxis()-> SetTitle("Events");
574  h_ele_EoP -> GetXaxis()-> SetTitle("E/P_{vertex}");
575  h_ele_EoP -> GetYaxis()-> SetTitle("Events");
576  h_ele_EseedOP -> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
577  h_ele_EseedOP -> GetYaxis()-> SetTitle("Events");
578  h_ele_EoPout -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
579  h_ele_EoPout -> GetYaxis()-> SetTitle("Events");
580  h_ele_EeleOPout -> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
581  h_ele_EeleOPout -> GetYaxis()-> SetTitle("Events");
582  h_ele_EoP_barrel -> GetXaxis()-> SetTitle("E/P_{vertex}");
583  h_ele_EoP_barrel -> GetYaxis()-> SetTitle("Events");
584  h_ele_EseedOP_barrel -> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
585  h_ele_EseedOP_barrel -> GetYaxis()-> SetTitle("Events");
586  h_ele_EoPout_barrel -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
587  h_ele_EoPout_barrel -> GetYaxis()-> SetTitle("Events");
588  h_ele_EeleOPout_barrel -> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
589  h_ele_EeleOPout_barrel -> GetYaxis()-> SetTitle("Events");
590  h_ele_EoP_endcaps -> GetXaxis()-> SetTitle("E/P_{vertex}");
591  h_ele_EoP_endcaps -> GetYaxis()-> SetTitle("Events");
592  h_ele_EseedOP_endcaps -> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
593  h_ele_EseedOP_endcaps -> GetYaxis()-> SetTitle("Events");
594  h_ele_EoPout_endcaps -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
595  h_ele_EoPout_endcaps -> GetYaxis()-> SetTitle("Events");
596  h_ele_EeleOPout_endcaps -> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
597  h_ele_EeleOPout_endcaps -> GetYaxis()-> SetTitle("Events");
598  h_ele_vertexX-> GetXaxis()-> SetTitle("x (cm)");
599  h_ele_vertexX-> GetYaxis()-> SetTitle("Events");
600  h_ele_vertexY-> GetXaxis()-> SetTitle("y (cm)");
601  h_ele_vertexY-> GetYaxis()-> SetTitle("Events");
602  h_ele_vertexZ-> GetXaxis()-> SetTitle("z (cm)");
603  h_ele_vertexZ-> GetYaxis()-> SetTitle("Events");
604  h_ele_vertexTIP-> GetXaxis()-> SetTitle("TIP (cm)");
605  h_ele_vertexTIP-> GetYaxis()-> SetTitle("Events");
606  h_ele_TIP_all-> GetXaxis()-> SetTitle("r_{T} (cm)");
607  h_ele_TIP_all-> GetYaxis()-> SetTitle("Events");
608  h_ele_vertexTIPVsEta-> GetYaxis()-> SetTitle("TIP (cm)");
609  h_ele_vertexTIPVsEta-> GetXaxis()-> SetTitle("#eta");
610  h_ele_vertexTIPVsPhi-> GetYaxis()-> SetTitle("TIP (cm)");
611  h_ele_vertexTIPVsPhi-> GetXaxis()-> SetTitle("#phi (rad)");
612  h_ele_vertexTIPVsPt-> GetYaxis()-> SetTitle("TIP (cm)");
613  h_ele_vertexTIPVsPt-> GetXaxis()-> SetTitle("p_{T} (GeV/c)");
614  h_ele_dEtaSc_propVtx-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
615  h_ele_dEtaSc_propVtx-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
616  h_ele_dEtaSc_propVtx-> GetYaxis()-> SetTitle("Events");
617  h_ele_dEtaCl_propOut-> GetXaxis()-> SetTitle("#eta_{seedcl} - #eta_{tr}");
618  h_ele_dEtaCl_propOut-> GetYaxis()-> SetTitle("Events");
619  h_ele_dEtaEleCl_propOut-> GetXaxis()-> SetTitle("#eta_{elecl} - #eta_{tr}");
620  h_ele_dEtaEleCl_propOut-> GetYaxis()-> SetTitle("Events");
621  h_ele_dPhiSc_propVtx-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
622  h_ele_dPhiSc_propVtx-> GetYaxis()-> SetTitle("Events");
623  h_ele_dPhiCl_propOut-> GetXaxis()-> SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
624  h_ele_dPhiCl_propOut-> GetYaxis()-> SetTitle("Events");
625  h_ele_dPhiEleCl_propOut-> GetXaxis()-> SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
626  h_ele_dPhiEleCl_propOut-> GetYaxis()-> SetTitle("Events");
627  h_ele_dEtaSc_propVtx_barrel-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
628  h_ele_dEtaSc_propVtx_barrel-> GetYaxis()-> SetTitle("Events");
629  h_ele_dEtaCl_propOut_barrel-> GetXaxis()-> SetTitle("#eta_{seedcl} - #eta_{tr}");
630  h_ele_dEtaCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
631  h_ele_dEtaEleCl_propOut_barrel-> GetXaxis()-> SetTitle("#eta_{elecl} - #eta_{tr}");
632  h_ele_dEtaEleCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
633  h_ele_dPhiSc_propVtx_barrel-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
634  h_ele_dPhiSc_propVtx_barrel-> GetYaxis()-> SetTitle("Events");
635  h_ele_dPhiCl_propOut_barrel-> GetXaxis()-> SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
636  h_ele_dPhiCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
637  h_ele_dPhiEleCl_propOut_barrel-> GetXaxis()-> SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
638  h_ele_dPhiEleCl_propOut_barrel-> GetYaxis()-> SetTitle("Events");
639  h_ele_dEtaSc_propVtx_endcaps-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
640  h_ele_dEtaSc_propVtx_endcaps-> GetYaxis()-> SetTitle("Events");
641  h_ele_dEtaCl_propOut_endcaps-> GetXaxis()-> SetTitle("#eta_{seedcl} - #eta_{tr}");
642  h_ele_dEtaCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
643  h_ele_dEtaEleCl_propOut_endcaps-> GetXaxis()-> SetTitle("#eta_{elecl} - #eta_{tr}");
644  h_ele_dEtaEleCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
645  h_ele_dPhiSc_propVtx_endcaps-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
646  h_ele_dPhiSc_propVtx_endcaps-> GetYaxis()-> SetTitle("Events");
647  h_ele_dPhiCl_propOut_endcaps-> GetXaxis()-> SetTitle("#phi_{seedcl} - #phi_{tr} (rad)");
648  h_ele_dPhiCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
649  h_ele_dPhiEleCl_propOut_endcaps-> GetXaxis()-> SetTitle("#phi_{elecl} - #phi_{tr} (rad)");
650  h_ele_dPhiEleCl_propOut_endcaps-> GetYaxis()-> SetTitle("Events");
651  h_ele_HoE-> GetXaxis()-> SetTitle("H/E") ;
652  h_ele_HoE-> GetYaxis()-> SetTitle("Events") ;
653  h_ele_HoE_barrel-> GetXaxis()-> SetTitle("H/E") ;
654  h_ele_HoE_barrel-> GetYaxis()-> SetTitle("Events") ;
655  h_ele_HoE_endcaps-> GetXaxis()-> SetTitle("H/E") ;
656  h_ele_HoE_endcaps-> GetYaxis()-> SetTitle("Events") ;
657  h_ele_HoE_fiducial-> GetXaxis()-> SetTitle("H/E") ;
658  h_ele_HoE_fiducial-> GetYaxis()-> SetTitle("Events") ;
659  h_ele_fbrem-> GetXaxis()-> SetTitle("P_{in} - P_{out} / P_{in}");
660  h_ele_fbrem-> GetYaxis()-> SetTitle("Events");
661  h_ele_seed_dphi2_-> GetXaxis()-> SetTitle("#phi_{hit}-#phi_{pred} (rad)") ;
662  h_ele_seed_dphi2_-> GetYaxis()-> SetTitle("Events") ;
663  h_ele_seed_drz2_-> GetXaxis()-> SetTitle("r(z)_{hit}-r(z)_{pred} (cm)") ;
664  h_ele_seed_drz2_-> GetYaxis()-> SetTitle("Events") ;
665  h_ele_seed_subdet2_-> GetXaxis()-> SetTitle("2nd hit subdet Id") ;
666  h_ele_seed_subdet2_-> GetYaxis()-> SetTitle("Events") ;
667  h_ele_classes-> GetXaxis()-> SetTitle("class Id") ;
668  h_ele_classes-> GetYaxis()-> SetTitle("Events") ;
669  h_ele_EoverP_all-> GetXaxis()-> SetTitle("E/P_{vertex}");
670  h_ele_EoverP_all-> GetYaxis()-> SetTitle("Events");
671  h_ele_EseedOP_all-> GetXaxis()-> SetTitle("E_{seed}/P_{vertex}");
672  h_ele_EseedOP_all-> GetYaxis()-> SetTitle("Events");
673  h_ele_EoPout_all -> GetXaxis()-> SetTitle("E_{seed}/P_{out}");
674  h_ele_EoPout_all-> GetYaxis()-> SetTitle("Events");
675  h_ele_EeleOPout_all-> GetXaxis()-> SetTitle("E_{ele}/P_{out}");
676  h_ele_EeleOPout_all-> GetYaxis()-> SetTitle("Events");
677  h_ele_dEtaSc_propVtx_all-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
678  h_ele_dEtaSc_propVtx_all-> GetYaxis()-> SetTitle("Events");
679  h_ele_dPhiSc_propVtx_all-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
680  h_ele_dPhiSc_propVtx_all-> GetYaxis()-> SetTitle("Events");
681  h_ele_dEtaCl_propOut_all-> GetXaxis()-> SetTitle("#eta_{sc} - #eta_{tr}");
682  h_ele_dEtaCl_propOut_all-> GetYaxis()-> SetTitle("Events");
683  h_ele_dPhiCl_propOut_all-> GetXaxis()-> SetTitle("#phi_{sc} - #phi_{tr} (rad)");
684  h_ele_dPhiCl_propOut_all-> GetYaxis()-> SetTitle("Events");
685  h_ele_HoE_all-> GetXaxis()-> SetTitle("H/E") ;
686  h_ele_HoE_all-> GetYaxis()-> SetTitle("Events");
687  h_ele_mee_all-> GetXaxis()-> SetTitle("m_{ee} (GeV/c^{2})");
688  h_ele_mee_all-> GetYaxis()-> SetTitle("Events");
689  h_ele_mee_os-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
690  h_ele_mee_os-> GetYaxis()-> SetTitle("Events");
691  h_ele_E2mnE1vsMee_all-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
692  h_ele_E2mnE1vsMee_all-> GetYaxis()-> SetTitle("E2 - E1 (GeV)");
693  h_ele_E2mnE1vsMee_egeg_all-> GetXaxis()-> SetTitle("m_{e^{+}e^{-}} (GeV/c^{2})");
694  h_ele_E2mnE1vsMee_egeg_all-> GetYaxis()-> SetTitle("E2 - E1 (GeV)");
695  histNum_-> GetXaxis()-> SetTitle("N_{ele}");
696  histNum_-> GetYaxis()-> SetTitle("Events");
697  h_ele_fbremVsEta_mode-> GetXaxis()-> SetTitle("#eta");
698  h_ele_fbremVsEta_mean-> GetXaxis()-> SetTitle("#eta");
699 
700 }
701 
702 void
703 GsfElectronMCFakeAnalyzer::endJob(){
704 
705  histfile_->cd();
706  std::cout << "efficiency calculation " << std::endl;
707  // efficiency vs eta
708  TH1F *h_ele_etaEff = (TH1F*)h_ele_matchingObjectEta_matched->Clone("h_ele_etaEff");
709  h_ele_etaEff->Reset();
710  h_ele_etaEff->Divide(h_ele_matchingObjectEta_matched,h_matchingObjectEta,1,1,"b");
711  h_ele_etaEff->Print();
712  h_ele_etaEff->GetXaxis()->SetTitle("#eta");
713  h_ele_etaEff->GetYaxis()->SetTitle("Efficiency");
714 
715  // efficiency vs z
716  TH1F *h_ele_zEff = (TH1F*)h_ele_matchingObjectZ_matched->Clone("h_ele_zEff");
717  h_ele_zEff->Reset();
718  h_ele_zEff->Divide(h_ele_matchingObjectZ_matched,h_matchingObjectZ,1,1,"b");
719  h_ele_zEff->Print();
720  h_ele_zEff->GetXaxis()->SetTitle("z (cm)");
721  h_ele_zEff->GetYaxis()->SetTitle("Efficiency");
722 
723  // efficiency vs |eta|
724  TH1F *h_ele_absetaEff = (TH1F*)h_ele_matchingObjectAbsEta_matched->Clone("h_ele_absetaEff");
725  h_ele_absetaEff->Reset();
726  h_ele_absetaEff->Divide(h_ele_matchingObjectAbsEta_matched,h_matchingObjectAbsEta,1,1,"b");
727  h_ele_absetaEff->GetXaxis()->SetTitle("|#eta|");
728  h_ele_absetaEff->GetYaxis()->SetTitle("Efficiency");
729 
730  // efficiency vs pt
731  TH1F *h_ele_ptEff = (TH1F*)h_ele_matchingObjectPt_matched->Clone("h_ele_ptEff");
732  h_ele_ptEff->Reset();
733  h_ele_ptEff->Divide(h_ele_matchingObjectPt_matched,h_matchingObjectPt,1,1,"b");
734  h_ele_ptEff->GetXaxis()->SetTitle("p_{T} (GeV/c)");
735  h_ele_ptEff->GetYaxis()->SetTitle("Efficiency");
736 
737  // efficiency vs phi
738  TH1F *h_ele_phiEff = (TH1F*)h_ele_matchingObjectPhi_matched->Clone("h_ele_phiEff");
739  h_ele_phiEff->Reset();
740  h_ele_phiEff->Divide(h_ele_matchingObjectPhi_matched,h_matchingObjectPhi,1,1,"b");
741  h_ele_phiEff->GetXaxis()->SetTitle("#phi (rad)");
742  h_ele_phiEff->GetYaxis()->SetTitle("Efficiency");
743 
744  // rec/matching objects all electrons
745  TH1F *h_ele_etaEff_all = (TH1F*)h_ele_vertexEta_all->Clone("h_ele_etaEff_all");
746  h_ele_etaEff_all->Reset();
747  h_ele_etaEff_all->Divide(h_ele_vertexEta_all,h_matchingObjectEta,1,1,"b");
748  h_ele_etaEff_all->Print();
749  h_ele_etaEff_all->GetXaxis()->SetTitle("#eta");
750  h_ele_etaEff_all->GetYaxis()->SetTitle("N_{rec}/N_{matching jet}");
751 
752  TH1F *h_ele_ptEff_all = (TH1F*)h_ele_vertexPt_all->Clone("h_ele_ptEff_all");
753  h_ele_ptEff_all->Reset();
754  h_ele_ptEff_all->Divide(h_ele_vertexPt_all,h_matchingObjectPt,1,1,"b");
755  h_ele_ptEff_all->Print();
756  h_ele_ptEff_all->GetXaxis()->SetTitle("p_{T} (GeV/c)");
757  h_ele_ptEff_all->GetYaxis()->SetTitle("N_{rec}/N_{matching gen jet}");
758 
759  // classes
760  TH1F *h_ele_eta_goldenFrac = (TH1F*)h_ele_eta_golden->Clone("h_ele_eta_goldenFrac");
761  h_ele_eta_goldenFrac->Reset();
762  h_ele_eta_goldenFrac->GetXaxis()->SetTitle("|#eta|");
763  h_ele_eta_goldenFrac->GetYaxis()->SetTitle("Fraction of electrons");
764  h_ele_eta_goldenFrac->Divide(h_ele_eta_golden,h_ele_eta,1,1);
765  h_ele_eta_goldenFrac->SetTitle("fraction of golden electrons vs eta");
766  TH1F *h_ele_eta_bbremFrac = (TH1F*)h_ele_eta_bbrem->Clone("h_ele_eta_bbremFrac");
767  h_ele_eta_bbremFrac->Reset();
768  h_ele_eta_bbremFrac->Divide(h_ele_eta_bbrem,h_ele_eta,1,1);
769  h_ele_eta_bbremFrac->GetXaxis()->SetTitle("|#eta|");
770  h_ele_eta_bbremFrac->GetYaxis()->SetTitle("Fraction of electrons");
771  h_ele_eta_bbremFrac->SetTitle("fraction of big brem electrons vs eta");
772  TH1F *h_ele_eta_narrowFrac = (TH1F*)h_ele_eta_narrow->Clone("h_ele_eta_narrowFrac");
773  h_ele_eta_narrowFrac->Reset();
774  h_ele_eta_narrowFrac->Divide(h_ele_eta_narrow,h_ele_eta,1,1);
775  h_ele_eta_narrowFrac->GetXaxis()->SetTitle("|#eta|");
776  h_ele_eta_narrowFrac->GetYaxis()->SetTitle("Fraction of electrons");
777  h_ele_eta_narrowFrac->SetTitle("fraction of narrow electrons vs eta");
778  TH1F *h_ele_eta_showerFrac = (TH1F*)h_ele_eta_shower->Clone("h_ele_eta_showerFrac");
779  h_ele_eta_showerFrac->Reset();
780  h_ele_eta_showerFrac->Divide(h_ele_eta_shower,h_ele_eta,1,1);
781  h_ele_eta_showerFrac->GetXaxis()->SetTitle("|#eta|");
782  h_ele_eta_showerFrac->GetYaxis()->SetTitle("Fraction of electrons");
783  h_ele_eta_showerFrac->SetTitle("fraction of showering electrons vs eta");
784 
785  // fbrem
786  TH1F *h_ele_xOverX0VsEta = new TH1F( "h_ele_xOverx0VsEta","mean X/X_0 vs eta",nbineta/2,0.0,2.5);
787  for (int ibin=1;ibin<h_ele_fbremVsEta_mean->GetNbinsX()+1;ibin++) {
788  double xOverX0 = 0.;
789  if (h_ele_fbremVsEta_mean->GetBinContent(ibin)>0.) xOverX0 = -log(h_ele_fbremVsEta_mean->GetBinContent(ibin));
790  h_ele_xOverX0VsEta->SetBinContent(ibin,xOverX0);
791  }
792 
793  //profiles from 2D histos
794  TProfile *p_ele_PoPmatchingObjectVsEta = h_ele_PoPmatchingObjectVsEta->ProfileX();
795  p_ele_PoPmatchingObjectVsEta->SetTitle("mean ele momentum / matching SC energy vs eta");
796  p_ele_PoPmatchingObjectVsEta->GetXaxis()->SetTitle("#eta");
797  p_ele_PoPmatchingObjectVsEta->GetYaxis()->SetTitle("<P/E_{matching gen jet}>");
798  p_ele_PoPmatchingObjectVsEta->Write();
799  TProfile *p_ele_PoPmatchingObjectVsPhi = h_ele_PoPmatchingObjectVsPhi->ProfileX();
800  p_ele_PoPmatchingObjectVsPhi->SetTitle("mean ele momentum / gen momentum vs phi");
801  p_ele_PoPmatchingObjectVsPhi->GetXaxis()->SetTitle("#phi (rad)");
802  p_ele_PoPmatchingObjectVsPhi->GetYaxis()->SetTitle("<P/E_{matching gen jet}>");
803  p_ele_PoPmatchingObjectVsPhi->Write();
804  TProfile *p_ele_EtaMnEtamatchingObjectVsEta = h_ele_EtaMnEtamatchingObjectVsEta->ProfileX();
805  p_ele_EtaMnEtamatchingObjectVsEta->GetXaxis()->SetTitle("#eta");
806  p_ele_EtaMnEtamatchingObjectVsEta->GetYaxis()->SetTitle("<#eta_{rec} - #eta_{matching gen jet}>");
807  p_ele_EtaMnEtamatchingObjectVsEta->Write();
808  TProfile *p_ele_EtaMnEtamatchingObjectVsPhi = h_ele_EtaMnEtamatchingObjectVsPhi->ProfileX();
809  p_ele_EtaMnEtamatchingObjectVsPhi-> GetXaxis()-> SetTitle("#phi");
810  p_ele_EtaMnEtamatchingObjectVsPhi-> GetYaxis()-> SetTitle("<#eta_{rec} - #eta_{matching gen jet}>");
811  p_ele_EtaMnEtamatchingObjectVsPhi->Write();
812  TProfile *p_ele_PhiMnPhimatchingObjectVsEta = h_ele_PhiMnPhimatchingObjectVsEta->ProfileX();
813  p_ele_PhiMnPhimatchingObjectVsEta-> GetXaxis()-> SetTitle("#eta");
814  p_ele_PhiMnPhimatchingObjectVsEta-> GetYaxis()-> SetTitle("<#phi_{rec} - #phi_{matching gen jet}> (rad)");
815  p_ele_PhiMnPhimatchingObjectVsEta->Write();
816  TProfile *p_ele_PhiMnPhimatchingObjectVsPhi = h_ele_PhiMnPhimatchingObjectVsPhi->ProfileX();
817  p_ele_PhiMnPhimatchingObjectVsPhi-> GetXaxis()-> SetTitle("#phi");
818  p_ele_PhiMnPhimatchingObjectVsPhi-> GetYaxis()-> SetTitle("<#phi_{rec} - #phi_{matching gen jet}> (rad)");
819  p_ele_PhiMnPhimatchingObjectVsPhi->Write();
820  TProfile *p_ele_vertexPtVsEta = h_ele_vertexPtVsEta->ProfileX();
821  p_ele_vertexPtVsEta->GetXaxis()->SetTitle("#eta");
822  p_ele_vertexPtVsEta->GetYaxis()->SetTitle("<p_{T}> (GeV/c)");
823  p_ele_vertexPtVsEta->Write();
824  TProfile *p_ele_vertexPtVsPhi = h_ele_vertexPtVsPhi->ProfileX();
825  p_ele_vertexPtVsPhi->GetXaxis()->SetTitle("#phi (rad)");
826  p_ele_vertexPtVsPhi->GetYaxis()->SetTitle("<p_{T}> (GeV/c)");
827  p_ele_vertexPtVsPhi->Write();
828  TProfile *p_ele_EoPVsEta = h_ele_EoPVsEta->ProfileX();
829  p_ele_EoPVsEta->GetXaxis()->SetTitle("#eta");
830  p_ele_EoPVsEta->GetYaxis()->SetTitle("<E/P_{vertex}>");
831  p_ele_EoPVsEta->Write();
832  TProfile *p_ele_EoPVsPhi = h_ele_EoPVsPhi->ProfileX();
833  p_ele_EoPVsPhi->GetXaxis()->SetTitle("#phi (rad)");
834  p_ele_EoPVsPhi->GetYaxis()->SetTitle("<E/P_{vertex}>");
835  p_ele_EoPVsPhi->Write();
836  TProfile *p_ele_EoPoutVsEta = h_ele_EoPoutVsEta->ProfileX();
837  p_ele_EoPoutVsEta->GetXaxis()->SetTitle("#eta");
838  p_ele_EoPoutVsEta->GetYaxis()->SetTitle("<E_{seed}/P_{out}>");
839  p_ele_EoPoutVsEta->Write();
840  TProfile *p_ele_EoPoutVsPhi = h_ele_EoPoutVsPhi->ProfileX();
841  p_ele_EoPoutVsPhi->GetXaxis()->SetTitle("#phi (rad)");
842  p_ele_EoPoutVsPhi->GetYaxis()->SetTitle("<E_{seed}/P_{out}>");
843  p_ele_EoPoutVsPhi->Write();
844  TProfile *p_ele_EeleOPoutVsEta = h_ele_EeleOPoutVsEta->ProfileX();
845  p_ele_EeleOPoutVsEta->SetTitle("mean ele Eele/pout vs eta");
846  p_ele_EeleOPoutVsEta->GetXaxis()->SetTitle("#eta");
847  p_ele_EeleOPoutVsEta->GetYaxis()->SetTitle("<E_{ele}/P_{out}>");
848  p_ele_EeleOPoutVsEta->Write();
849  TProfile *p_ele_EeleOPoutVsPhi = h_ele_EeleOPoutVsPhi->ProfileX();
850  p_ele_EeleOPoutVsPhi->SetTitle("mean ele Eele/pout vs phi");
851  p_ele_EeleOPoutVsPhi->GetXaxis()->SetTitle("#phi (rad)");
852  p_ele_EeleOPoutVsPhi->GetYaxis()->SetTitle("<E_{ele}/P_{out}>");
853  p_ele_EeleOPoutVsPhi->Write();
854  TProfile *p_ele_HoEVsEta = h_ele_HoEVsEta->ProfileX();
855  p_ele_HoEVsEta->GetXaxis()->SetTitle("#eta");
856  p_ele_HoEVsEta->GetYaxis()->SetTitle("<H/E>");
857  p_ele_HoEVsEta->Write();
858  TProfile *p_ele_HoEVsPhi = h_ele_HoEVsPhi->ProfileX();
859  p_ele_HoEVsPhi->GetXaxis()->SetTitle("#phi (rad)");
860  p_ele_HoEVsPhi->GetYaxis()->SetTitle("<H/E>");
861  p_ele_HoEVsPhi->Write();
862  TProfile *p_ele_chi2VsEta = h_ele_chi2VsEta->ProfileX();
863  p_ele_chi2VsEta->GetXaxis()->SetTitle("#eta");
864  p_ele_chi2VsEta->GetYaxis()->SetTitle("<#Chi^{2}>");
865  p_ele_chi2VsEta->Write();
866  TProfile *p_ele_chi2VsPhi = h_ele_chi2VsPhi->ProfileX();
867  p_ele_chi2VsPhi->GetXaxis()->SetTitle("#phi (rad)");
868  p_ele_chi2VsPhi->GetYaxis()->SetTitle("<#Chi^{2}>");
869  p_ele_chi2VsPhi->Write();
870  TProfile *p_ele_foundHitsVsEta = h_ele_foundHitsVsEta->ProfileX();
871  p_ele_foundHitsVsEta->GetXaxis()->SetTitle("#eta");
872  p_ele_foundHitsVsEta->GetYaxis()->SetTitle("<N_{hits}>");
873  p_ele_foundHitsVsEta->Write();
874  TProfile *p_ele_foundHitsVsPhi = h_ele_foundHitsVsPhi->ProfileX();
875  p_ele_foundHitsVsPhi->GetXaxis()->SetTitle("#phi (rad)");
876  p_ele_foundHitsVsPhi->GetYaxis()->SetTitle("<N_{hits}>");
877  p_ele_foundHitsVsPhi->Write();
878  TProfile *p_ele_lostHitsVsEta = h_ele_lostHitsVsEta->ProfileX();
879  p_ele_lostHitsVsEta->GetXaxis()->SetTitle("#eta");
880  p_ele_lostHitsVsEta->GetYaxis()->SetTitle("<N_{hits}>");
881  p_ele_lostHitsVsEta->Write();
882  TProfile *p_ele_lostHitsVsPhi = h_ele_lostHitsVsPhi->ProfileX();
883  p_ele_lostHitsVsPhi->GetXaxis()->SetTitle("#phi (rad)");
884  p_ele_lostHitsVsPhi->GetYaxis()->SetTitle("<N_{hits}>");
885  p_ele_lostHitsVsPhi->Write();
886  TProfile *p_ele_vertexTIPVsEta = h_ele_vertexTIPVsEta->ProfileX();
887  p_ele_vertexTIPVsEta->SetTitle("mean tip (wrt gen vtx) vs eta");
888  p_ele_vertexTIPVsEta->GetXaxis()->SetTitle("#eta");
889  p_ele_vertexTIPVsEta->GetYaxis()->SetTitle("<TIP> (cm)");
890  p_ele_vertexTIPVsEta->Write();
891  TProfile *p_ele_vertexTIPVsPhi = h_ele_vertexTIPVsPhi->ProfileX();
892  p_ele_vertexTIPVsPhi->SetTitle("mean tip (wrt gen vtx) vs phi");
893  p_ele_vertexTIPVsPhi->GetXaxis()->SetTitle("#phi");
894  p_ele_vertexTIPVsPhi->GetYaxis()->SetTitle("<TIP> (cm)");
895  p_ele_vertexTIPVsPhi->Write();
896  TProfile *p_ele_vertexTIPVsPt = h_ele_vertexTIPVsPt->ProfileX();
897  p_ele_vertexTIPVsPt->SetTitle("mean tip (wrt gen vtx) vs phi");
898  p_ele_vertexTIPVsPt->GetXaxis()->SetTitle("p_{T} (GeV/c)");
899  p_ele_vertexTIPVsPt->GetYaxis()->SetTitle("<TIP> (cm)");
900  p_ele_vertexTIPVsPt->Write();
901 
902  // mc truth
903 
904  h_matchingObjectNum->Write();
905 
906  // rec event
907 
908  histNum_->Write();
909 
910  // mc
911  h_matchingObjectEta->Write();
912  h_matchingObjectAbsEta->Write();
913  h_matchingObjectP->Write();
914  h_matchingObjectPt->Write();
915  h_matchingObjectPhi->Write();
916  h_matchingObjectZ->Write();
917 
918  // all electrons
919  h_ele_EoverP_all->Write();
920  h_ele_EseedOP_all->Write();
921  h_ele_EoPout_all->Write();
922  h_ele_EeleOPout_all ->Write();
923  h_ele_dEtaSc_propVtx_all->Write();
924  h_ele_dPhiSc_propVtx_all->Write();
925  h_ele_dEtaCl_propOut_all ->Write();
926  h_ele_dPhiCl_propOut_all->Write();
927  h_ele_HoE_all->Write();
928  h_ele_TIP_all->Write();
929  h_ele_vertexPt_all->Write();
930  h_ele_vertexEta_all->Write();
931  h_ele_mee_all->Write();
932  h_ele_mee_os->Write();
933  h_ele_E2mnE1vsMee_all ->Write();
934  h_ele_E2mnE1vsMee_egeg_all->Write();
935 
936  // matched electrons
937  h_ele_charge->Write();
938  h_ele_chargeVsEta->Write();
939  h_ele_chargeVsPhi->Write();
940  h_ele_chargeVsPt->Write();
941  h_ele_vertexP->Write();
942  h_ele_vertexPt->Write();
943  h_ele_vertexPtVsEta->Write();
944  h_ele_vertexPtVsPhi->Write();
945  h_ele_matchingObjectPt_matched->Write();
946  h_ele_vertexEta->Write();
947  h_ele_vertexEtaVsPhi->Write();
948  h_ele_matchingObjectAbsEta_matched->Write();
949  h_ele_matchingObjectEta_matched->Write();
950  h_ele_vertexPhi->Write();
951  h_ele_matchingObjectPhi_matched->Write();
952  h_ele_vertexX->Write();
953  h_ele_vertexY ->Write();
954  h_ele_vertexZ->Write();
955  h_ele_vertexTIP->Write();
956  h_ele_matchingObjectZ_matched->Write();
957  h_ele_vertexTIPVsEta->Write();
958  h_ele_vertexTIPVsPhi->Write();
959  h_ele_vertexTIPVsPt->Write();
960  h_ele_PoPmatchingObject->Write();
961  h_ele_PoPmatchingObjectVsEta ->Write();
962  h_ele_PoPmatchingObjectVsPhi->Write();
963  h_ele_PoPmatchingObjectVsPt->Write();
964  h_ele_PoPmatchingObject_barrel ->Write();
965  h_ele_PoPmatchingObject_endcaps->Write();
966  h_ele_EtaMnEtamatchingObject->Write();
967  h_ele_EtaMnEtamatchingObjectVsEta ->Write();
968  h_ele_EtaMnEtamatchingObjectVsPhi->Write();
969  h_ele_EtaMnEtamatchingObjectVsPt->Write();
970  h_ele_PhiMnPhimatchingObject ->Write();
971  h_ele_PhiMnPhimatchingObject2 ->Write();
972  h_ele_PhiMnPhimatchingObjectVsEta->Write();
973  h_ele_PhiMnPhimatchingObjectVsPhi->Write();
974  h_ele_PhiMnPhimatchingObjectVsPt->Write();
975 
976  // matched electron, superclusters
977  histSclEn_->Write();
978  histSclEoEmatchingObject_barrel->Write();
979  histSclEoEmatchingObject_endcaps->Write();
980  histSclEt_->Write();
981  histSclEtVsEta_->Write();
982  histSclEtVsPhi_->Write();
983  histSclEtaVsPhi_ ->Write();
984  histSclEta_->Write();
985  histSclPhi_->Write();
986  histSclSigIEtaIEta_->Write();
987  histSclSigIEtaIEta_barrel_->Write();
988  histSclSigIEtaIEta_endcaps_->Write();
989  histSclE1x5_->Write();
990  histSclE1x5_barrel_->Write();
991  histSclE1x5_endcaps_->Write();
992  histSclE2x5max_->Write();
993  histSclE2x5max_barrel_->Write();
994  histSclE2x5max_endcaps_->Write();
995  histSclE5x5_->Write();
996  histSclE5x5_barrel_->Write();
997  histSclE5x5_endcaps_->Write();
998 
999  // matched electron, gsf tracks
1000  h_ele_ambiguousTracks->Write();
1001  h_ele_ambiguousTracksVsEta->Write();
1002  h_ele_ambiguousTracksVsPhi->Write();
1003  h_ele_ambiguousTracksVsPt->Write();
1004  h_ele_foundHits->Write();
1005  h_ele_foundHitsVsEta->Write();
1006  h_ele_foundHitsVsPhi->Write();
1007  h_ele_foundHitsVsPt->Write();
1008  h_ele_lostHits->Write();
1009  h_ele_lostHitsVsEta->Write();
1010  h_ele_lostHitsVsPhi->Write();
1011  h_ele_lostHitsVsPt->Write();
1012  h_ele_chi2 ->Write();
1013  h_ele_chi2VsEta ->Write();
1014  h_ele_chi2VsPhi ->Write();
1015  h_ele_chi2VsPt->Write();
1016  h_ele_PinMnPout->Write();
1017  h_ele_PinMnPout_mode->Write();
1018  h_ele_PinMnPoutVsEta_mode->Write();
1019  h_ele_PinMnPoutVsPhi_mode->Write();
1020  h_ele_PinMnPoutVsPt_mode->Write();
1021  h_ele_PinMnPoutVsE_mode->Write();
1022  h_ele_PinMnPoutVsChi2_mode->Write();
1023  h_ele_outerP ->Write();
1024  h_ele_outerP_mode->Write();
1025  h_ele_outerPVsEta_mode->Write();
1026  h_ele_outerPt->Write();
1027  h_ele_outerPt_mode ->Write();
1028  h_ele_outerPtVsEta_mode->Write();
1029  h_ele_outerPtVsPhi_mode->Write();
1030  h_ele_outerPtVsPt_mode->Write();
1031 
1032  // matched electrons, matching
1033  h_ele_EoP ->Write();
1034  h_ele_EoP_eg ->Write();
1035  h_ele_EoP_barrel ->Write();
1036  h_ele_EoP_eg_barrel ->Write();
1037  h_ele_EoP_endcaps ->Write();
1038  h_ele_EoP_eg_endcaps ->Write();
1039  h_ele_EoPVsEta ->Write();
1040  h_ele_EoPVsPhi->Write();
1041  h_ele_EoPVsE->Write();
1042  h_ele_EseedOP ->Write();
1043  h_ele_EseedOP_eg ->Write();
1044  h_ele_EseedOP_barrel ->Write();
1045  h_ele_EseedOP_eg_barrel ->Write();
1046  h_ele_EseedOP_endcaps ->Write();
1047  h_ele_EseedOP_eg_endcaps ->Write();
1048  h_ele_EseedOPVsEta ->Write();
1049  h_ele_EseedOPVsPhi->Write();
1050  h_ele_EseedOPVsE->Write();
1051  h_ele_EoPout->Write();
1052  h_ele_EoPout_eg->Write();
1053  h_ele_EoPout_barrel->Write();
1054  h_ele_EoPout_eg_barrel->Write();
1055  h_ele_EoPout_endcaps->Write();
1056  h_ele_EoPout_eg_endcaps->Write();
1057  h_ele_EoPoutVsEta->Write();
1058  h_ele_EoPoutVsPhi->Write();
1059  h_ele_EoPoutVsE ->Write();
1060  h_ele_EeleOPout->Write();
1061  h_ele_EeleOPout_eg->Write();
1062  h_ele_EeleOPout_barrel->Write();
1063  h_ele_EeleOPout_eg_barrel->Write();
1064  h_ele_EeleOPout_endcaps->Write();
1065  h_ele_EeleOPout_eg_endcaps->Write();
1066  h_ele_EeleOPoutVsEta->Write();
1067  h_ele_EeleOPoutVsPhi->Write();
1068  h_ele_EeleOPoutVsE ->Write();
1069  h_ele_dEtaSc_propVtx->Write();
1070  h_ele_dEtaSc_propVtx_eg->Write();
1071  h_ele_dEtaSc_propVtx_barrel->Write();
1072  h_ele_dEtaSc_propVtx_eg_barrel->Write();
1073  h_ele_dEtaSc_propVtx_endcaps->Write();
1074  h_ele_dEtaSc_propVtx_eg_endcaps->Write();
1075  h_ele_dEtaScVsEta_propVtx->Write();
1076  h_ele_dEtaScVsPhi_propVtx->Write();
1077  h_ele_dEtaScVsPt_propVtx ->Write();
1078  h_ele_dPhiSc_propVtx->Write();
1079  h_ele_dPhiSc_propVtx_eg->Write();
1080  h_ele_dPhiSc_propVtx_barrel->Write();
1081  h_ele_dPhiSc_propVtx_eg_barrel->Write();
1082  h_ele_dPhiSc_propVtx_endcaps->Write();
1083  h_ele_dPhiSc_propVtx_eg_endcaps->Write();
1084  h_ele_dPhiScVsEta_propVtx ->Write();
1085  h_ele_dPhiScVsPhi_propVtx->Write();
1086  h_ele_dPhiScVsPt_propVtx->Write();
1087  h_ele_dEtaCl_propOut->Write();
1088  h_ele_dEtaCl_propOut_eg->Write();
1089  h_ele_dEtaCl_propOut_barrel->Write();
1090  h_ele_dEtaCl_propOut_eg_barrel->Write();
1091  h_ele_dEtaCl_propOut_endcaps->Write();
1092  h_ele_dEtaCl_propOut_eg_endcaps->Write();
1093  h_ele_dEtaClVsEta_propOut->Write();
1094  h_ele_dEtaClVsPhi_propOut->Write();
1095  h_ele_dEtaClVsPt_propOut->Write();
1096  h_ele_dPhiCl_propOut->Write();
1097  h_ele_dPhiCl_propOut_eg->Write();
1098  h_ele_dPhiCl_propOut_barrel->Write();
1099  h_ele_dPhiCl_propOut_eg_barrel->Write();
1100  h_ele_dPhiCl_propOut_endcaps->Write();
1101  h_ele_dPhiCl_propOut_eg_endcaps->Write();
1102  h_ele_dPhiClVsEta_propOut->Write();
1103  h_ele_dPhiClVsPhi_propOut->Write();
1104  h_ele_dPhiClVsPt_propOut->Write();
1105  h_ele_dEtaEleCl_propOut->Write();
1106  h_ele_dEtaEleCl_propOut_eg->Write();
1107  h_ele_dEtaEleCl_propOut_barrel->Write();
1108  h_ele_dEtaEleCl_propOut_eg_barrel->Write();
1109  h_ele_dEtaEleCl_propOut_endcaps->Write();
1110  h_ele_dEtaEleCl_propOut_eg_endcaps->Write();
1111  h_ele_dEtaEleClVsEta_propOut->Write();
1112  h_ele_dEtaEleClVsPhi_propOut->Write();
1113  h_ele_dEtaEleClVsPt_propOut->Write();
1114  h_ele_dPhiEleCl_propOut->Write();
1115  h_ele_dPhiEleCl_propOut_eg->Write();
1116  h_ele_dPhiEleCl_propOut_barrel->Write();
1117  h_ele_dPhiEleCl_propOut_eg_barrel->Write();
1118  h_ele_dPhiEleCl_propOut_endcaps->Write();
1119  h_ele_dPhiEleCl_propOut_eg_endcaps->Write();
1120  h_ele_dPhiEleClVsEta_propOut->Write();
1121  h_ele_dPhiEleClVsPhi_propOut->Write();
1122  h_ele_dPhiEleClVsPt_propOut->Write();
1123  h_ele_HoE->Write();
1124  h_ele_HoE_eg->Write();
1125  h_ele_HoE_barrel->Write();
1126  h_ele_HoE_eg_barrel->Write();
1127  h_ele_HoE_endcaps->Write();
1128  h_ele_HoE_eg_endcaps->Write();
1129  h_ele_HoE_fiducial->Write();
1130  h_ele_HoEVsEta->Write();
1131  h_ele_HoEVsPhi->Write();
1132  h_ele_HoEVsE->Write();
1133 
1134  h_ele_seed_dphi2_->Write();
1135  h_ele_seed_subdet2_->Write();
1136  TProfile *p_ele_seed_dphi2VsEta_ = h_ele_seed_dphi2VsEta_->ProfileX();
1137  p_ele_seed_dphi2VsEta_->SetTitle("mean ele seed dphi 2nd layer vs eta");
1138  p_ele_seed_dphi2VsEta_->GetXaxis()->SetTitle("#eta");
1139  p_ele_seed_dphi2VsEta_->GetYaxis()->SetTitle("<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)");
1140  p_ele_seed_dphi2VsEta_->SetMinimum(-0.004);
1141  p_ele_seed_dphi2VsEta_->SetMaximum(0.004);
1142  p_ele_seed_dphi2VsEta_->Write();
1143  TProfile *p_ele_seed_dphi2VsPt_ = h_ele_seed_dphi2VsPt_->ProfileX();
1144  p_ele_seed_dphi2VsPt_->SetTitle("mean ele seed dphi 2nd layer vs pt");
1145  p_ele_seed_dphi2VsPt_->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1146  p_ele_seed_dphi2VsPt_->GetYaxis()->SetTitle("<#phi_{pred} - #phi_{hit}, 2nd layer> (rad)");
1147  p_ele_seed_dphi2VsPt_->SetMinimum(-0.004);
1148  p_ele_seed_dphi2VsPt_->SetMaximum(0.004);
1149  p_ele_seed_dphi2VsPt_->Write();
1150  h_ele_seed_drz2_->Write();
1151  TProfile *p_ele_seed_drz2VsEta_ = h_ele_seed_drz2VsEta_->ProfileX();
1152  p_ele_seed_drz2VsEta_->SetTitle("mean ele seed dr(dz) 2nd layer vs eta");
1153  p_ele_seed_drz2VsEta_->GetXaxis()->SetTitle("#eta");
1154  p_ele_seed_drz2VsEta_->GetYaxis()->SetTitle("<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)");
1155  p_ele_seed_drz2VsEta_->SetMinimum(-0.15);
1156  p_ele_seed_drz2VsEta_->SetMaximum(0.15);
1157  p_ele_seed_drz2VsEta_->Write();
1158  TProfile *p_ele_seed_drz2VsPt_ = h_ele_seed_drz2VsPt_->ProfileX();
1159  p_ele_seed_drz2VsPt_->SetTitle("mean ele seed dr(dz) 2nd layer vs pt");
1160  p_ele_seed_drz2VsPt_->GetXaxis()->SetTitle("p_{T} (GeV/c)");
1161  p_ele_seed_drz2VsPt_->GetYaxis()->SetTitle("<r(z)_{pred} - r(z)_{hit}, 2nd layer> (cm)");
1162  p_ele_seed_drz2VsPt_->SetMinimum(-0.15);
1163  p_ele_seed_drz2VsPt_->SetMaximum(0.15);
1164  p_ele_seed_drz2VsPt_->Write();
1165 
1166  // classes
1167  h_ele_classes->Write();
1168  h_ele_eta->Write();
1169  h_ele_eta_golden->Write();
1170  h_ele_eta_bbrem->Write();
1171  h_ele_eta_narrow->Write();
1172  h_ele_eta_shower->Write();
1173  h_ele_PinVsPoutGolden_mode->Write();
1174  h_ele_PinVsPoutShowering_mode->Write();
1175  h_ele_PinVsPoutGolden_mean->Write();
1176  h_ele_PinVsPoutShowering_mean->Write();
1177  h_ele_PtinVsPtoutGolden_mode->Write();
1178  h_ele_PtinVsPtoutShowering_mode->Write();
1179  h_ele_PtinVsPtoutGolden_mean->Write();
1180  h_ele_PtinVsPtoutShowering_mean->Write();
1181  histSclEoEmatchingObjectGolden_barrel->Write();
1182  histSclEoEmatchingObjectGolden_endcaps->Write();
1183  histSclEoEmatchingObjectShowering_barrel->Write();
1184  histSclEoEmatchingObjectShowering_endcaps->Write();
1185 
1186  // fbrem
1187  h_ele_fbrem->Write();
1188  h_ele_fbremVsEta_mode->Write();
1189  h_ele_fbremVsEta_mean->Write();
1190  h_ele_etaEff->Write();
1191  h_ele_zEff->Write();
1192  h_ele_phiEff->Write();
1193  h_ele_absetaEff->Write();
1194  h_ele_ptEff->Write();
1195  h_ele_etaEff_all->Write();
1196  h_ele_ptEff_all->Write();
1197  h_ele_eta_goldenFrac->Write();
1198  h_ele_eta_bbremFrac->Write();
1199  h_ele_eta_narrowFrac->Write();
1200  h_ele_eta_showerFrac->Write();
1201  h_ele_xOverX0VsEta->Write();
1202 
1203  // e/g et pflow electrons
1204  h_ele_mva->Write();
1205  h_ele_provenance->Write();
1206 
1207  // isolation
1208  h_ele_tkSumPt_dr03->GetXaxis()->SetTitle("TkIsoSum, cone 0.3 (GeV/c)");
1209  h_ele_tkSumPt_dr03->GetYaxis()->SetTitle("Events");
1210  h_ele_tkSumPt_dr03->Write();
1211  h_ele_ecalRecHitSumEt_dr03->GetXaxis()->SetTitle("EcalIsoSum, cone 0.3 (GeV)");
1212  h_ele_ecalRecHitSumEt_dr03->GetYaxis()->SetTitle("Events");
1213  h_ele_ecalRecHitSumEt_dr03->Write();
1214  h_ele_hcalDepth1TowerSumEt_dr03->GetXaxis()->SetTitle("Hcal1IsoSum, cone 0.3 (GeV)");
1215  h_ele_hcalDepth1TowerSumEt_dr03->GetYaxis()->SetTitle("Events");
1216  h_ele_hcalDepth1TowerSumEt_dr03->Write();
1217  h_ele_hcalDepth2TowerSumEt_dr03->GetXaxis()->SetTitle("Hcal2IsoSum, cone 0.3 (GeV)");
1218  h_ele_hcalDepth2TowerSumEt_dr03->GetYaxis()->SetTitle("Events");
1219  h_ele_hcalDepth2TowerSumEt_dr03->Write();
1220  h_ele_tkSumPt_dr04->GetXaxis()->SetTitle("TkIsoSum, cone 0.4 (GeV/c)");
1221  h_ele_tkSumPt_dr04->GetYaxis()->SetTitle("Events");
1222  h_ele_tkSumPt_dr04->Write();
1223  h_ele_ecalRecHitSumEt_dr04->GetXaxis()->SetTitle("EcalIsoSum, cone 0.4 (GeV)");
1224  h_ele_ecalRecHitSumEt_dr04->GetYaxis()->SetTitle("Events");
1225  h_ele_ecalRecHitSumEt_dr04->Write();
1226  h_ele_hcalDepth1TowerSumEt_dr04->GetXaxis()->SetTitle("Hcal1IsoSum, cone 0.4 (GeV)");
1227  h_ele_hcalDepth1TowerSumEt_dr04->GetYaxis()->SetTitle("Events");
1228  h_ele_hcalDepth1TowerSumEt_dr04->Write();
1229  h_ele_hcalDepth2TowerSumEt_dr04->GetXaxis()->SetTitle("Hcal2IsoSum, cone 0.4 (GeV)");
1230  h_ele_hcalDepth2TowerSumEt_dr04->GetYaxis()->SetTitle("Events");
1231  h_ele_hcalDepth2TowerSumEt_dr04->Write();
1232 
1233 }
1234 
1235 void
1236 GsfElectronMCFakeAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
1237 {
1238  std::cout << "analyzing new event " << std::endl;
1239 
1240  // get reco electrons
1241  edm::Handle<reco::GsfElectronCollection> gsfElectrons;
1242  iEvent.getByLabel(electronCollection_,gsfElectrons);
1243  edm::LogInfo("")<<"\n\n =================> Treating event "<<iEvent.id()<<" Number of electrons "<<gsfElectrons.product()->size();
1244 
1245  // get gen jets
1246  edm::Handle<reco::GenJetCollection> genJets ;
1247  iEvent.getByLabel(matchingObjectCollection_,genJets);
1248 
1249  // get the beamspot from the Event:
1250  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
1251  iEvent.getByLabel(beamSpot_, recoBeamSpotHandle);
1252  const BeamSpot bs = *recoBeamSpotHandle;
1253 
1254  histNum_->Fill((*gsfElectrons).size());
1255 
1256  // all rec electrons
1257  for (reco::GsfElectronCollection::const_iterator gsfIter=gsfElectrons->begin();
1258  gsfIter!=gsfElectrons->end(); gsfIter++){
1259  // preselect electrons
1260  if (gsfIter->pt()>maxPt_ || std::abs(gsfIter->eta())>maxAbsEta_) continue;
1261  h_ele_EoverP_all -> Fill( gsfIter->eSuperClusterOverP() );
1262  h_ele_EseedOP_all -> Fill( gsfIter->eSeedClusterOverP() );
1263  h_ele_EoPout_all -> Fill( gsfIter->eSeedClusterOverPout() );
1264  h_ele_EeleOPout_all -> Fill( gsfIter->eEleClusterOverPout() );
1265  h_ele_dEtaSc_propVtx_all -> Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
1266  h_ele_dPhiSc_propVtx_all -> Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
1267  h_ele_dEtaCl_propOut_all -> Fill(gsfIter->deltaEtaSeedClusterTrackAtCalo());
1268  h_ele_dPhiCl_propOut_all -> Fill(gsfIter->deltaPhiSeedClusterTrackAtCalo());
1269  h_ele_HoE_all -> Fill( gsfIter->hadronicOverEm() );
1270  double d = gsfIter->vertex().x()*gsfIter->vertex().x()
1271  +gsfIter->vertex().y()*gsfIter->vertex().y();
1272  h_ele_TIP_all -> Fill( sqrt(d) );
1273  h_ele_vertexEta_all -> Fill( gsfIter->eta() );
1274  h_ele_vertexPt_all -> Fill( gsfIter->pt() );
1275  float enrj1=gsfIter->superCluster()->energy();
1276  // mee
1277  for (reco::GsfElectronCollection::const_iterator gsfIter2=gsfIter+1;
1278  gsfIter2!=gsfElectrons->end(); gsfIter2++){
1279  math::XYZTLorentzVector p12 = (*gsfIter).p4()+(*gsfIter2).p4();
1280  float mee2 = p12.Dot(p12);
1281  h_ele_mee_all -> Fill(sqrt(mee2));
1282  float enrj2=gsfIter2->superCluster()->energy();
1283  h_ele_E2mnE1vsMee_all->Fill(sqrt(mee2),enrj2-enrj1);
1284  if (gsfIter->ecalDrivenSeed() && gsfIter2->ecalDrivenSeed()) h_ele_E2mnE1vsMee_egeg_all->Fill(sqrt(mee2),enrj2-enrj1);
1285  if (gsfIter->charge()*gsfIter2->charge()<0.) h_ele_mee_os -> Fill(sqrt(mee2));
1286  }
1287  }
1288 
1289  // association matching object-reco electrons
1290  int matchingObjectNum=0;
1291 
1292  for (reco::GenJetCollection::const_iterator moIter = genJets->begin();
1293  moIter != genJets->end(); ++moIter) {
1294 
1295  // number of matching objects
1296  matchingObjectNum++;
1297 
1298  if (moIter->energy()/cosh(moIter->eta())> maxPt_ || std::abs(moIter->eta())> maxAbsEta_) continue;
1299 
1300  // suppress the endcaps
1301  //if (std::abs(moIter->eta()) > 1.5) continue;
1302  // select central z
1303  //if ( std::abs((*mcIter)->production_vertex()->position().z())>50.) continue;
1304 
1305  h_matchingObjectEta -> Fill( moIter->eta() );
1306  h_matchingObjectAbsEta -> Fill( std::abs(moIter->eta()) );
1307  h_matchingObjectP -> Fill( moIter->energy() );
1308  h_matchingObjectPt -> Fill( moIter->energy()/cosh(moIter->eta()) );
1309  h_matchingObjectPhi -> Fill( moIter->phi() );
1310  h_matchingObjectZ -> Fill( moIter->vz() );
1311 
1312  // looking for the best matching gsf electron
1313  bool okGsfFound = false;
1314  double gsfOkRatio = 999999.;
1315 
1316  // find best matched electron
1317  reco::GsfElectron bestGsfElectron;
1318  for (reco::GsfElectronCollection::const_iterator gsfIter=gsfElectrons->begin();
1319  gsfIter!=gsfElectrons->end(); gsfIter++){
1320 
1321  double dphi = gsfIter->phi()-moIter->phi();
1322  if (std::abs(dphi)>CLHEP::pi)
1323  dphi = dphi < 0? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi;
1324  double deltaR = sqrt(std::pow((gsfIter->eta()-moIter->eta()),2) + std::pow(dphi,2));
1325  if ( deltaR < deltaR_ ){
1326  //if ( (genPc->pdg_id() == 11) && (gsfIter->charge() < 0.) || (genPc->pdg_id() == -11) &&
1327  //(gsfIter->charge() > 0.) ){
1328  double tmpGsfRatio = gsfIter->p()/moIter->energy();
1329  if ( std::abs(tmpGsfRatio-1) < std::abs(gsfOkRatio-1) ) {
1330  gsfOkRatio = tmpGsfRatio;
1331  bestGsfElectron=*gsfIter;
1332  okGsfFound = true;
1333  }
1334  //}
1335  }
1336  } // loop over rec ele to look for the best one
1337 
1338  // analysis when the matching object is matched by a rec electron
1339  if (okGsfFound){
1340 
1341  // electron related distributions
1342  h_ele_charge -> Fill( bestGsfElectron.charge() );
1343  h_ele_chargeVsEta -> Fill( bestGsfElectron.eta(),bestGsfElectron.charge() );
1344  h_ele_chargeVsPhi -> Fill( bestGsfElectron.phi(),bestGsfElectron.charge() );
1345  h_ele_chargeVsPt -> Fill( bestGsfElectron.pt(),bestGsfElectron.charge() );
1346  h_ele_vertexP -> Fill( bestGsfElectron.p() );
1347  h_ele_vertexPt -> Fill( bestGsfElectron.pt() );
1348  h_ele_vertexPtVsEta -> Fill( bestGsfElectron.eta(),bestGsfElectron.pt() );
1349  h_ele_vertexPtVsPhi -> Fill( bestGsfElectron.phi(),bestGsfElectron.pt() );
1350  h_ele_vertexEta -> Fill( bestGsfElectron.eta() );
1351  // generated distributions for matched electrons
1352  h_ele_matchingObjectPt_matched -> Fill( moIter->energy()/cosh(moIter->eta()) );
1353  h_ele_matchingObjectPhi_matched -> Fill( moIter->phi() );
1354  h_ele_matchingObjectAbsEta_matched -> Fill( std::abs(moIter->eta()) );
1355  h_ele_matchingObjectEta_matched -> Fill( moIter->eta() );
1356  h_ele_vertexEtaVsPhi -> Fill( bestGsfElectron.phi(),bestGsfElectron.eta() );
1357  h_ele_vertexPhi -> Fill( bestGsfElectron.phi() );
1358  h_ele_vertexX -> Fill( bestGsfElectron.vertex().x() );
1359  h_ele_vertexY -> Fill( bestGsfElectron.vertex().y() );
1360  h_ele_vertexZ -> Fill( bestGsfElectron.vertex().z() );
1361  h_ele_matchingObjectZ_matched -> Fill( moIter->vz() );
1362  double d = (bestGsfElectron.vertex().x()-bs.position().x())
1363  *(bestGsfElectron.vertex().x()-bs.position().x())+
1364  (bestGsfElectron.vertex().y()-bs.position().y())
1365  *(bestGsfElectron.vertex().y()-bs.position().y());
1366  d = sqrt(d);
1367  h_ele_vertexTIP -> Fill( d );
1368  h_ele_vertexTIPVsEta -> Fill( bestGsfElectron.eta(), d );
1369  h_ele_vertexTIPVsPhi -> Fill( bestGsfElectron.phi(), d );
1370  h_ele_vertexTIPVsPt -> Fill( bestGsfElectron.pt(), d );
1371  h_ele_EtaMnEtamatchingObject -> Fill( bestGsfElectron.eta()-moIter->eta());
1372  h_ele_EtaMnEtamatchingObjectVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.eta()-moIter->eta());
1373  h_ele_EtaMnEtamatchingObjectVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.eta()-moIter->eta());
1374  h_ele_EtaMnEtamatchingObjectVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.eta()-moIter->eta());
1375  h_ele_PhiMnPhimatchingObject -> Fill( bestGsfElectron.phi()-moIter->phi());
1376  h_ele_PhiMnPhimatchingObject2 -> Fill( bestGsfElectron.phi()-moIter->phi());
1377  h_ele_PhiMnPhimatchingObjectVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.phi()-moIter->phi());
1378  h_ele_PhiMnPhimatchingObjectVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.phi()-moIter->phi());
1379  h_ele_PhiMnPhimatchingObjectVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.phi()-moIter->phi());
1380  h_ele_PoPmatchingObject -> Fill( bestGsfElectron.p()/moIter->energy());
1381  h_ele_PoPmatchingObjectVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.p()/moIter->energy());
1382  h_ele_PoPmatchingObjectVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.p()/moIter->energy());
1383  h_ele_PoPmatchingObjectVsPt -> Fill( bestGsfElectron.py(), bestGsfElectron.p()/moIter->energy());
1384  if (bestGsfElectron.isEB()) h_ele_PoPmatchingObject_barrel -> Fill( bestGsfElectron.p()/moIter->energy());
1385  if (bestGsfElectron.isEE()) h_ele_PoPmatchingObject_endcaps -> Fill( bestGsfElectron.p()/moIter->energy());
1386 
1387  // supercluster related distributions
1388  reco::SuperClusterRef sclRef = bestGsfElectron.superCluster();
1389  if (!bestGsfElectron.ecalDrivenSeed()&&bestGsfElectron.trackerDrivenSeed()) sclRef = bestGsfElectron.pflowSuperCluster();
1390  histSclEn_->Fill(sclRef->energy());
1391  double R=TMath::Sqrt(sclRef->x()*sclRef->x() + sclRef->y()*sclRef->y() +sclRef->z()*sclRef->z());
1392  double Rt=TMath::Sqrt(sclRef->x()*sclRef->x() + sclRef->y()*sclRef->y());
1393  histSclEt_->Fill(sclRef->energy()*(Rt/R));
1394  histSclEtVsEta_->Fill(sclRef->eta(),sclRef->energy()*(Rt/R));
1395  histSclEtVsPhi_->Fill(sclRef->phi(),sclRef->energy()*(Rt/R));
1396  if (bestGsfElectron.classification() < 100) histSclEoEmatchingObject_barrel->Fill(sclRef->energy()/moIter->energy());
1397  if (bestGsfElectron.classification() >= 100) histSclEoEmatchingObject_endcaps->Fill(sclRef->energy()/moIter->energy());
1398  histSclEta_->Fill(sclRef->eta());
1399  histSclEtaVsPhi_->Fill(sclRef->phi(),sclRef->eta());
1400  histSclPhi_->Fill(sclRef->phi());
1401  histSclSigIEtaIEta_->Fill(bestGsfElectron.scSigmaIEtaIEta());
1402  if (bestGsfElectron.isEB()) histSclSigIEtaIEta_barrel_->Fill(bestGsfElectron.scSigmaIEtaIEta());
1403  if (bestGsfElectron.isEE()) histSclSigIEtaIEta_endcaps_->Fill(bestGsfElectron.scSigmaIEtaIEta());
1404  histSclE1x5_->Fill(bestGsfElectron.scE1x5());
1405  if (bestGsfElectron.isEB()) histSclE1x5_barrel_->Fill(bestGsfElectron.scE1x5());
1406  if (bestGsfElectron.isEE()) histSclE1x5_endcaps_->Fill(bestGsfElectron.scE1x5());
1407  histSclE2x5max_->Fill(bestGsfElectron.scE2x5Max());
1408  if (bestGsfElectron.isEB()) histSclE2x5max_barrel_->Fill(bestGsfElectron.scE2x5Max());
1409  if (bestGsfElectron.isEE()) histSclE2x5max_endcaps_->Fill(bestGsfElectron.scE2x5Max());
1410  histSclE5x5_->Fill(bestGsfElectron.scE5x5());
1411  if (bestGsfElectron.isEB()) histSclE5x5_barrel_->Fill(bestGsfElectron.scE5x5());
1412  if (bestGsfElectron.isEE()) histSclE5x5_endcaps_->Fill(bestGsfElectron.scE5x5());
1413 
1414  // track related distributions
1415  h_ele_ambiguousTracks -> Fill( bestGsfElectron.ambiguousGsfTracksSize() );
1416  h_ele_ambiguousTracksVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.ambiguousGsfTracksSize() );
1417  h_ele_ambiguousTracksVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.ambiguousGsfTracksSize() );
1418  h_ele_ambiguousTracksVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.ambiguousGsfTracksSize() );
1419  if (!readAOD_) { // track extra does not exist in AOD
1420  h_ele_foundHits -> Fill( bestGsfElectron.gsfTrack()->numberOfValidHits() );
1421  h_ele_foundHitsVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->numberOfValidHits() );
1422  h_ele_foundHitsVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->numberOfValidHits() );
1423  h_ele_foundHitsVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->numberOfValidHits() );
1424  h_ele_lostHits -> Fill( bestGsfElectron.gsfTrack()->numberOfLostHits() );
1425  h_ele_lostHitsVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->numberOfLostHits() );
1426  h_ele_lostHitsVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->numberOfLostHits() );
1427  h_ele_lostHitsVsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->numberOfLostHits() );
1428  h_ele_chi2 -> Fill( bestGsfElectron.gsfTrack()->normalizedChi2() );
1429  h_ele_chi2VsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.gsfTrack()->normalizedChi2() );
1430  h_ele_chi2VsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.gsfTrack()->normalizedChi2() );
1431  h_ele_chi2VsPt -> Fill( bestGsfElectron.pt(), bestGsfElectron.gsfTrack()->normalizedChi2() );
1432  }
1433  // from gsf track interface, hence using mean
1434  if (!readAOD_) { // track extra does not exist in AOD
1435  h_ele_PinMnPout -> Fill( bestGsfElectron.gsfTrack()->innerMomentum().R() - bestGsfElectron.gsfTrack()->outerMomentum().R() );
1436  h_ele_outerP -> Fill( bestGsfElectron.gsfTrack()->outerMomentum().R() );
1437  h_ele_outerPt -> Fill( bestGsfElectron.gsfTrack()->outerMomentum().Rho() );
1438  }
1439  // from electron interface, hence using mode
1440  h_ele_PinMnPout_mode -> Fill( bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1441  h_ele_PinMnPoutVsEta_mode -> Fill( bestGsfElectron.eta(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1442  h_ele_PinMnPoutVsPhi_mode -> Fill( bestGsfElectron.phi(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1443  h_ele_PinMnPoutVsPt_mode -> Fill( bestGsfElectron.pt(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1444  h_ele_PinMnPoutVsE_mode -> Fill( bestGsfElectron.caloEnergy(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1445  if (!readAOD_) // track extra does not exist in AOD
1446  h_ele_PinMnPoutVsChi2_mode -> Fill( bestGsfElectron.gsfTrack()->normalizedChi2(), bestGsfElectron.trackMomentumAtVtx().R() - bestGsfElectron.trackMomentumOut().R() );
1447  h_ele_outerP_mode -> Fill( bestGsfElectron.trackMomentumOut().R() );
1448  h_ele_outerPVsEta_mode -> Fill(bestGsfElectron.eta(), bestGsfElectron.trackMomentumOut().R() );
1449  h_ele_outerPt_mode -> Fill( bestGsfElectron.trackMomentumOut().Rho() );
1450  h_ele_outerPtVsEta_mode -> Fill(bestGsfElectron.eta(), bestGsfElectron.trackMomentumOut().Rho() );
1451  h_ele_outerPtVsPhi_mode -> Fill(bestGsfElectron.phi(), bestGsfElectron.trackMomentumOut().Rho() );
1452  h_ele_outerPtVsPt_mode -> Fill(bestGsfElectron.pt(), bestGsfElectron.trackMomentumOut().Rho() );
1453 
1454  if (!readAOD_) { // track extra does not exist in AOD
1455  edm::RefToBase<TrajectorySeed> seed = bestGsfElectron.gsfTrack()->extra()->seedRef();
1456  ElectronSeedRef elseed=seed.castTo<ElectronSeedRef>();
1457  h_ele_seed_dphi2_-> Fill(elseed->dPhi2());
1458  h_ele_seed_dphi2VsEta_-> Fill(bestGsfElectron.eta(), elseed->dPhi2());
1459  h_ele_seed_dphi2VsPt_-> Fill(bestGsfElectron.pt(), elseed->dPhi2()) ;
1460  h_ele_seed_drz2_-> Fill(elseed->dRz2());
1461  h_ele_seed_drz2VsEta_-> Fill(bestGsfElectron.eta(), elseed->dRz2());
1462  h_ele_seed_drz2VsPt_-> Fill(bestGsfElectron.pt(), elseed->dRz2());
1463  h_ele_seed_subdet2_-> Fill(elseed->subDet2());
1464  }
1465  // match distributions
1466  h_ele_EoP -> Fill( bestGsfElectron.eSuperClusterOverP() );
1467  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EoP_eg -> Fill( bestGsfElectron.eSuperClusterOverP() );
1468  if (bestGsfElectron.isEB()) h_ele_EoP_barrel -> Fill( bestGsfElectron.eSuperClusterOverP() );
1469  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoP_eg_barrel -> Fill( bestGsfElectron.eSuperClusterOverP() );
1470  if (bestGsfElectron.isEE()) h_ele_EoP_endcaps -> Fill( bestGsfElectron.eSuperClusterOverP() );
1471  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoP_eg_endcaps -> Fill( bestGsfElectron.eSuperClusterOverP() );
1472  h_ele_EoPVsEta -> Fill(bestGsfElectron.eta(), bestGsfElectron.eSuperClusterOverP() );
1473  h_ele_EoPVsPhi -> Fill(bestGsfElectron.phi(), bestGsfElectron.eSuperClusterOverP() );
1474  h_ele_EoPVsE -> Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.eSuperClusterOverP() );
1475  h_ele_EseedOP -> Fill( bestGsfElectron.eSeedClusterOverP() );
1476  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EseedOP_eg -> Fill( bestGsfElectron.eSeedClusterOverP() );
1477  if (bestGsfElectron.isEB()) h_ele_EseedOP_barrel -> Fill( bestGsfElectron.eSeedClusterOverP() );
1478  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EseedOP_eg_barrel -> Fill( bestGsfElectron.eSeedClusterOverP() );
1479  if (bestGsfElectron.isEE()) h_ele_EseedOP_endcaps -> Fill( bestGsfElectron.eSeedClusterOverP() );
1480  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EseedOP_eg_endcaps -> Fill( bestGsfElectron.eSeedClusterOverP() );
1481  h_ele_EseedOPVsEta -> Fill(bestGsfElectron.eta(), bestGsfElectron.eSeedClusterOverP() );
1482  h_ele_EseedOPVsPhi -> Fill(bestGsfElectron.phi(), bestGsfElectron.eSeedClusterOverP() );
1483  h_ele_EseedOPVsE -> Fill(bestGsfElectron.caloEnergy(), bestGsfElectron.eSeedClusterOverP() );
1484  h_ele_EoPout -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1485  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EoPout_eg -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1486  if (bestGsfElectron.isEB()) h_ele_EoPout_barrel -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1487  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoPout_eg_barrel -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1488  if (bestGsfElectron.isEE()) h_ele_EoPout_endcaps -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1489  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EoPout_eg_endcaps -> Fill( bestGsfElectron.eSeedClusterOverPout() );
1490  h_ele_EoPoutVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.eSeedClusterOverPout() );
1491  h_ele_EoPoutVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.eSeedClusterOverPout() );
1492  h_ele_EoPoutVsE -> Fill( bestGsfElectron.caloEnergy(), bestGsfElectron.eSeedClusterOverPout() );
1493  h_ele_EeleOPout -> Fill( bestGsfElectron.eEleClusterOverPout() );
1494  if (bestGsfElectron.ecalDrivenSeed()) h_ele_EeleOPout_eg -> Fill( bestGsfElectron.eEleClusterOverPout() );
1495  if (bestGsfElectron.isEB()) h_ele_EeleOPout_barrel -> Fill( bestGsfElectron.eEleClusterOverPout() );
1496  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EeleOPout_eg_barrel -> Fill( bestGsfElectron.eEleClusterOverPout() );
1497  if (bestGsfElectron.isEE()) h_ele_EeleOPout_endcaps -> Fill( bestGsfElectron.eEleClusterOverPout() );
1498  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_EeleOPout_eg_endcaps -> Fill( bestGsfElectron.eEleClusterOverPout() );
1499  h_ele_EeleOPoutVsEta -> Fill( bestGsfElectron.eta(), bestGsfElectron.eEleClusterOverPout() );
1500  h_ele_EeleOPoutVsPhi -> Fill( bestGsfElectron.phi(), bestGsfElectron.eEleClusterOverPout() );
1501  h_ele_EeleOPoutVsE -> Fill( bestGsfElectron.caloEnergy(), bestGsfElectron.eEleClusterOverPout() );
1502  h_ele_dEtaSc_propVtx -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1503  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaSc_propVtx_eg -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1504  if (bestGsfElectron.isEB()) h_ele_dEtaSc_propVtx_barrel -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1505  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaSc_propVtx_eg_barrel -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1506  if (bestGsfElectron.isEE())h_ele_dEtaSc_propVtx_endcaps -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1507  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaSc_propVtx_eg_endcaps -> Fill(bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1508  h_ele_dEtaScVsEta_propVtx -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1509  h_ele_dEtaScVsPhi_propVtx -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1510  h_ele_dEtaScVsPt_propVtx -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaEtaSuperClusterTrackAtVtx());
1511  h_ele_dPhiSc_propVtx -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1512  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiSc_propVtx_eg -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1513  if (bestGsfElectron.isEB()) h_ele_dPhiSc_propVtx_barrel -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1514  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiSc_propVtx_eg_barrel -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1515  if (bestGsfElectron.isEE())h_ele_dPhiSc_propVtx_endcaps -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1516  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiSc_propVtx_eg_endcaps -> Fill(bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1517  h_ele_dPhiScVsEta_propVtx -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1518  h_ele_dPhiScVsPhi_propVtx -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1519  h_ele_dPhiScVsPt_propVtx -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaPhiSuperClusterTrackAtVtx());
1520  h_ele_dEtaCl_propOut -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1521  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaCl_propOut_eg -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1522  if (bestGsfElectron.isEB()) h_ele_dEtaCl_propOut_barrel -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1523  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1524  if (bestGsfElectron.isEE()) h_ele_dEtaCl_propOut_endcaps -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1525  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1526  h_ele_dEtaClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1527  h_ele_dEtaClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1528  h_ele_dEtaClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaEtaSeedClusterTrackAtCalo());
1529  h_ele_dPhiCl_propOut -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1530  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiCl_propOut_eg -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1531  if (bestGsfElectron.isEB()) h_ele_dPhiCl_propOut_barrel -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1532  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1533  if (bestGsfElectron.isEE()) h_ele_dPhiCl_propOut_endcaps -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1534  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1535  h_ele_dPhiClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1536  h_ele_dPhiClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1537  h_ele_dPhiClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaPhiSeedClusterTrackAtCalo());
1538  h_ele_dEtaEleCl_propOut -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1539  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaEleCl_propOut_eg -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1540  if (bestGsfElectron.isEB()) h_ele_dEtaEleCl_propOut_barrel -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1541  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaEleCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1542  if (bestGsfElectron.isEE()) h_ele_dEtaEleCl_propOut_endcaps -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1543  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dEtaEleCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1544  h_ele_dEtaEleClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1545  h_ele_dEtaEleClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1546  h_ele_dEtaEleClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaEtaEleClusterTrackAtCalo());
1547  h_ele_dPhiEleCl_propOut -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1548  if (bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiEleCl_propOut_eg -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1549  if (bestGsfElectron.isEB()) h_ele_dPhiEleCl_propOut_barrel -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1550  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiEleCl_propOut_eg_barrel -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1551  if (bestGsfElectron.isEE()) h_ele_dPhiEleCl_propOut_endcaps -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1552  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_dPhiEleCl_propOut_eg_endcaps -> Fill(bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1553  h_ele_dPhiEleClVsEta_propOut -> Fill( bestGsfElectron.eta(),bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1554  h_ele_dPhiEleClVsPhi_propOut -> Fill(bestGsfElectron.phi(),bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1555  h_ele_dPhiEleClVsPt_propOut -> Fill(bestGsfElectron.pt(),bestGsfElectron.deltaPhiEleClusterTrackAtCalo());
1556  h_ele_HoE -> Fill(bestGsfElectron.hadronicOverEm());
1557  if (bestGsfElectron.ecalDrivenSeed()) h_ele_HoE_eg -> Fill(bestGsfElectron.hadronicOverEm());
1558  if (bestGsfElectron.isEB()) h_ele_HoE_barrel -> Fill(bestGsfElectron.hadronicOverEm());
1559  if (bestGsfElectron.isEB()&&bestGsfElectron.ecalDrivenSeed()) h_ele_HoE_eg_barrel -> Fill(bestGsfElectron.hadronicOverEm());
1560  if (bestGsfElectron.isEE()) h_ele_HoE_endcaps -> Fill(bestGsfElectron.hadronicOverEm());
1561  if (bestGsfElectron.isEE()&&bestGsfElectron.ecalDrivenSeed()) h_ele_HoE_eg_endcaps -> Fill(bestGsfElectron.hadronicOverEm());
1562  if (!bestGsfElectron.isEBEtaGap() && !bestGsfElectron.isEBPhiGap()&& !bestGsfElectron.isEBEEGap() &&
1563  !bestGsfElectron.isEERingGap() && !bestGsfElectron.isEEDeeGap()) h_ele_HoE_fiducial -> Fill(bestGsfElectron.hadronicOverEm());
1564  h_ele_HoEVsEta -> Fill( bestGsfElectron.eta(),bestGsfElectron.hadronicOverEm());
1565  h_ele_HoEVsPhi -> Fill(bestGsfElectron.phi(),bestGsfElectron.hadronicOverEm());
1566  h_ele_HoEVsE -> Fill(bestGsfElectron.caloEnergy(),bestGsfElectron.hadronicOverEm());
1567 
1568  //classes
1569  int eleClass = bestGsfElectron.classification();
1570  if (bestGsfElectron.isEE()) eleClass+=10;
1571  h_ele_classes ->Fill(eleClass);
1572 
1573  h_ele_eta->Fill(std::abs(bestGsfElectron.eta()));
1574  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_eta_golden ->Fill(std::abs(bestGsfElectron.eta()));
1575  if (bestGsfElectron.classification() == GsfElectron::BIGBREM) h_ele_eta_bbrem ->Fill(std::abs(bestGsfElectron.eta()));
1576  //if (bestGsfElectron.classification() == GsfElectron::OLDNARROW) h_ele_eta_narrow ->Fill(std::abs(bestGsfElectron.eta()));
1577  if (bestGsfElectron.classification() == GsfElectron::SHOWERING) h_ele_eta_shower ->Fill(std::abs(bestGsfElectron.eta()));
1578 
1579  //fbrem
1580  double fbrem_mean=0.;
1581  if (!readAOD_) // track extra does not exist in AOD
1582  fbrem_mean = 1. - bestGsfElectron.gsfTrack()->outerMomentum().R()/bestGsfElectron.gsfTrack()->innerMomentum().R();
1583  double fbrem_mode = bestGsfElectron.fbrem();
1584  h_ele_fbrem->Fill(fbrem_mode);
1585  h_ele_fbremVsEta_mode->Fill(bestGsfElectron.eta(),fbrem_mode);
1586  if (!readAOD_) // track extra does not exist in AOD
1587  h_ele_fbremVsEta_mean->Fill(bestGsfElectron.eta(),fbrem_mean);
1588 
1589  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PinVsPoutGolden_mode -> Fill(bestGsfElectron.trackMomentumOut().R(), bestGsfElectron.trackMomentumAtVtx().R());
1590  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
1591  h_ele_PinVsPoutShowering_mode -> Fill(bestGsfElectron.trackMomentumOut().R(), bestGsfElectron.trackMomentumAtVtx().R());
1592  if (!readAOD_) // track extra does not exist in AOD
1593  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PinVsPoutGolden_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().R(), bestGsfElectron.gsfTrack()->innerMomentum().R());
1594  if (!readAOD_) // track extra does not exist in AOD
1595  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
1596  h_ele_PinVsPoutShowering_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().R(), bestGsfElectron.gsfTrack()->innerMomentum().R());
1597  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PtinVsPtoutGolden_mode -> Fill(bestGsfElectron.trackMomentumOut().Rho(), bestGsfElectron.trackMomentumAtVtx().Rho());
1598  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
1599  h_ele_PtinVsPtoutShowering_mode -> Fill(bestGsfElectron.trackMomentumOut().Rho(), bestGsfElectron.trackMomentumAtVtx().Rho());
1600  if (!readAOD_) // track extra does not exist in AOD
1601  if (bestGsfElectron.classification() == GsfElectron::GOLDEN) h_ele_PtinVsPtoutGolden_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().Rho(), bestGsfElectron.gsfTrack()->innerMomentum().Rho());
1602  if (!readAOD_) // track extra does not exist in AOD
1603  if (bestGsfElectron.classification() == GsfElectron::SHOWERING)
1604  h_ele_PtinVsPtoutShowering_mean -> Fill(bestGsfElectron.gsfTrack()->outerMomentum().Rho(), bestGsfElectron.gsfTrack()->innerMomentum().Rho());
1605 
1606  h_ele_mva->Fill(bestGsfElectron.mva());
1607  if (bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(1.);
1608  if (bestGsfElectron.trackerDrivenSeed()) h_ele_provenance->Fill(-1.);
1609  if (bestGsfElectron.trackerDrivenSeed()||bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(0.);
1610  if (bestGsfElectron.trackerDrivenSeed()&&!bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(-2.);
1611  if (!bestGsfElectron.trackerDrivenSeed()&&bestGsfElectron.ecalDrivenSeed()) h_ele_provenance->Fill(2.);
1612 
1613  h_ele_tkSumPt_dr03->Fill(bestGsfElectron.dr03TkSumPt());
1614  h_ele_ecalRecHitSumEt_dr03->Fill(bestGsfElectron.dr03EcalRecHitSumEt());
1615  h_ele_hcalDepth1TowerSumEt_dr03->Fill(bestGsfElectron.dr03HcalDepth1TowerSumEt());
1616  h_ele_hcalDepth2TowerSumEt_dr03->Fill(bestGsfElectron.dr03HcalDepth2TowerSumEt());
1617  h_ele_tkSumPt_dr04->Fill(bestGsfElectron.dr04TkSumPt());
1618  h_ele_ecalRecHitSumEt_dr04->Fill(bestGsfElectron.dr04EcalRecHitSumEt());
1619  h_ele_hcalDepth1TowerSumEt_dr04->Fill(bestGsfElectron.dr04HcalDepth1TowerSumEt());
1620  h_ele_hcalDepth2TowerSumEt_dr04->Fill(bestGsfElectron.dr04HcalDepth2TowerSumEt());
1621 
1622  } // gsf electron found
1623 
1624 // } // matching object found
1625 
1626 // }
1627 
1628  } // loop overmatching object
1629 
1630  h_matchingObjectNum->Fill(matchingObjectNum);
1631 
1632 }
1633 
1634 
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
reco::GsfElectron::isEBEtaGap
bool isEBEtaGap() const
Definition: GsfElectron.h:335
reco::GsfElectron::eSuperClusterOverP
float eSuperClusterOverP() const
Definition: GsfElectron.h:229
reco::LeafCandidate::p
virtual double p() const GCC11_FINAL
magnitude of momentum vector
Definition: LeafCandidate.h:132
reco::GsfElectron::isEBEEGap
bool isEBEEGap() const
Definition: GsfElectron.h:333
reco::GsfElectron::dr04HcalDepth1TowerSumEt
float dr04HcalDepth1TowerSumEt() const
Definition: GsfElectron.h:449
reco::LeafCandidate::vertex
virtual const Point & vertex() const
vertex position (overwritten by PF...)
Definition: LeafCandidate.h:200
reco::GsfElectron::trackMomentumAtVtx
math::XYZVectorF trackMomentumAtVtx() const
Definition: GsfElectron.h:272
create_public_lumi_plots.log
log
Definition: create_public_lumi_plots.py:1119
abs
#define abs(x)
Definition: mlp_lapack.h:159
reco::GsfElectron::dr04TkSumPt
float dr04TkSumPt() const
Definition: GsfElectron.h:447
reco::GsfElectron::isEERingGap
bool isEERingGap() const
Definition: GsfElectron.h:339
reco::GsfElectron::dr03HcalDepth2TowerSumEt
float dr03HcalDepth2TowerSumEt() const
Definition: GsfElectron.h:439
reco::GsfElectron::fbrem
float fbrem() const
Definition: GsfElectron.h:639
reco::GsfElectron::scE2x5Max
float scE2x5Max() const
Definition: GsfElectron.h:401
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:21
reco::GsfElectron::isEEDeeGap
bool isEEDeeGap() const
Definition: GsfElectron.h:338
reco::LeafCandidate::py
virtual double py() const GCC11_FINAL
y coordinate of momentum vector
Definition: LeafCandidate.h:149
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
reco::GsfElectron::isEE
bool isEE() const
Definition: GsfElectron.h:331
reco::GsfElectron::isEB
bool isEB() const
Definition: GsfElectron.h:330
GsfElectronMCFakeAnalyzer::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: GsfElectronMCFakeAnalyzer.cc:1236
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:30
reco::GsfElectron::dr04HcalDepth2TowerSumEt
float dr04HcalDepth2TowerSumEt() const
Definition: GsfElectron.h:450
reco::GsfElectron::mva
float mva() const
Definition: GsfElectron.h:567
reco::LeafCandidate::phi
virtual float phi() const GCC11_FINAL
momentum azimuthal angle
Definition: LeafCandidate.h:155
reco::GsfElectron::deltaEtaSuperClusterTrackAtVtx
float deltaEtaSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:233
iEvent
int iEvent
Definition: GenABIO.cc:243
reco::GsfElectron::hadronicOverEm
float hadronicOverEm() const
Definition: GsfElectron.h:403
reco::GsfElectron::eSeedClusterOverP
float eSeedClusterOverP() const
Definition: GsfElectron.h:230
reco::GsfElectron::dr04EcalRecHitSumEt
float dr04EcalRecHitSumEt() const
Definition: GsfElectron.h:448
reco::GsfElectron::deltaPhiSeedClusterTrackAtCalo
float deltaPhiSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:237
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
L1Trigger_dataformats.reco
dictionary reco
Definition: L1Trigger_dataformats.py:9
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
reco::GsfElectron::superCluster
virtual SuperClusterRef superCluster() const
reference to a SuperCluster
Definition: GsfElectron.h:168
reco::GsfElectron::trackMomentumOut
math::XYZVectorF trackMomentumOut() const
Definition: GsfElectron.h:274
reco::GsfElectron::deltaPhiSuperClusterTrackAtVtx
float deltaPhiSuperClusterTrackAtVtx() const
Definition: GsfElectron.h:236
reco::GsfElectron::dr03TkSumPt
float dr03TkSumPt() const
Definition: GsfElectron.h:436
reco::GsfElectron::dr03HcalDepth1TowerSumEt
float dr03HcalDepth1TowerSumEt() const
Definition: GsfElectron.h:438
reco::GsfElectron::eEleClusterOverPout
float eEleClusterOverPout() const
Definition: GsfElectron.h:232
reco::GsfElectron::eSeedClusterOverPout
float eSeedClusterOverPout() const
Definition: GsfElectron.h:231
reco::GsfElectron::scSigmaIEtaIEta
float scSigmaIEtaIEta() const
Definition: GsfElectron.h:399
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:361
reco::GsfElectron::deltaPhiEleClusterTrackAtCalo
float deltaPhiEleClusterTrackAtCalo() const
Definition: GsfElectron.h:238
dbtoconf.conf
tuple conf
Definition: dbtoconf.py:185
reco::LeafCandidate::eta
virtual float eta() const GCC11_FINAL
momentum pseudorapidity
Definition: LeafCandidate.h:159
deltaR
double deltaR(double eta1, double eta2, double phi1, double phi2)
Definition: TreeUtility.cc:17
reco::GsfElectron::deltaEtaEleClusterTrackAtCalo
float deltaEtaEleClusterTrackAtCalo() const
Definition: GsfElectron.h:235
reco::LeafCandidate::charge
virtual int charge() const GCC11_FINAL
electric charge
Definition: LeafCandidate.h:115
reco::GsfElectron::classification
Classification classification() const
Definition: GsfElectron.h:635
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:588
edm::Handle::product
T const * product() const
Definition: Handle.h:74
reco::GsfElectron::dr03EcalRecHitSumEt
float dr03EcalRecHitSumEt() const
Definition: GsfElectron.h:437
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:42
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
reco::GsfElectron::deltaEtaSeedClusterTrackAtCalo
float deltaEtaSeedClusterTrackAtCalo() const
Definition: GsfElectron.h:234
reco::GsfElectron::isEBPhiGap
bool isEBPhiGap() const
Definition: GsfElectron.h:336
reco::GsfElectron::trackerDrivenSeed
bool trackerDrivenSeed() const
Definition: GsfElectron.h:173
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:63
reco::GsfElectron::scE5x5
float scE5x5() const
Definition: GsfElectron.h:402
pi
double pi
Definition: Vx3DHLTAnalyzer.h:63
reco::LeafCandidate::pt
virtual float pt() const GCC11_FINAL
transverse momentum
Definition: LeafCandidate.h:153
reco::GsfElectron::caloEnergy
float caloEnergy() const
Definition: GsfElectron.h:727
reco::GsfElectron::pflowSuperCluster
SuperClusterRef pflowSuperCluster() const
Definition: GsfElectron.h:174
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
GsfElectronMCFakeAnalyzer::GsfElectronMCFakeAnalyzer
GsfElectronMCFakeAnalyzer(const edm::ParameterSet &conf)
Definition: GsfElectronMCFakeAnalyzer.cc:52
reco::GsfElectron::scE1x5
float scE1x5() const
Definition: GsfElectron.h:400
reco::GsfElectron::gsfTrack
virtual GsfTrackRef gsfTrack() const
reference to a GsfTrack
Definition: GsfElectron.h:169
reco::GsfElectron::ecalDrivenSeed
bool ecalDrivenSeed() const
Definition: GsfElectron.h:172
gsfElectrons_cfi.gsfElectrons
tuple gsfElectrons
Definition: gsfElectrons_cfi.py:148
Generated for CMSSW Reference Manual by doxygen 1.8.5