CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
ElectronAnalyzer.cc
Go to the documentation of this file.
1 #include "DQMOffline/EGamma/plugins/ElectronAnalyzer.h"
2 
3 #include "DQMServices/Core/interface/MonitorElement.h"
4 
5 #include "DataFormats/GsfTrackReco/interface/GsfTrack.h"
6 #include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
7 #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
8 #include "DataFormats/Common/interface/Handle.h"
9 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
10 #include "DataFormats/Common/interface/TriggerResults.h"
11 
12 #include "FWCore/Common/interface/TriggerNames.h"
13 #include "FWCore/ParameterSet/interface/ParameterSet.h"
14 #include "FWCore/Framework/interface/EDAnalyzer.h"
15 #include "FWCore/Framework/interface/Event.h"
16 #include "FWCore/Framework/interface/MakerMacros.h"
17 #include "FWCore/MessageLogger/interface/MessageLogger.h"
18 #include "FWCore/ServiceRegistry/interface/Service.h"
19 
20 #include "CLHEP/Units/GlobalPhysicalConstants.h"
21 #include "TMath.h"
22 
23 #include <iostream>
24 
25 using namespace reco ;
26 
27 ElectronAnalyzer::ElectronAnalyzer( const edm::ParameterSet & conf )
28  : ElectronDqmAnalyzerBase(conf)
29  {
30  // general, collections
31  Selection_ = conf.getParameter<int>("Selection");
32  electronCollection_ = consumes<GsfElectronCollection>(conf.getParameter<edm::InputTag>("ElectronCollection"));
33  matchingObjectCollection_ = consumes<SuperClusterCollection>(conf.getParameter<edm::InputTag>("MatchingObjectCollection"));
34  trackCollection_ = consumes<TrackCollection>(conf.getParameter<edm::InputTag>("TrackCollection"));
35  vertexCollection_ = consumes<VertexCollection>(conf.getParameter<edm::InputTag>("VertexCollection"));
36  gsftrackCollection_ = consumes<GsfTrackCollection>(conf.getParameter<edm::InputTag>("GsfTrackCollection"));
37  beamSpotTag_ = consumes<BeamSpot>(conf.getParameter<edm::InputTag>("BeamSpot"));
38  readAOD_ = conf.getParameter<bool>("ReadAOD");
39 
40  // matching
41  matchingCondition_ = conf.getParameter<std::string>("MatchingCondition");
42  assert (matchingCondition_=="Cone") ;
43  maxPtMatchingObject_ = conf.getParameter<double>("MaxPtMatchingObject");
44  maxAbsEtaMatchingObject_ = conf.getParameter<double>("MaxAbsEtaMatchingObject");
45  deltaR_ = conf.getParameter<double>("DeltaR");
46 
47  // electron selection
48  minEt_ = conf.getParameter<double>("MinEt");
49  minPt_ = conf.getParameter<double>("MinPt");
50  maxAbsEta_ = conf.getParameter<double>("MaxAbsEta");
51  isEB_ = conf.getParameter<bool>("SelectEb");
52  isEE_ = conf.getParameter<bool>("SelectEe");
53  isNotEBEEGap_ = conf.getParameter<bool>("SelectNotEbEeGap");
54  isEcalDriven_ = conf.getParameter<bool>("SelectEcalDriven");
55  isTrackerDriven_ = conf.getParameter<bool>("SelectTrackerDriven");
56  eOverPMinBarrel_ = conf.getParameter<double>("MinEopBarrel");
57  eOverPMaxBarrel_ = conf.getParameter<double>("MaxEopBarrel");
58  eOverPMinEndcaps_ = conf.getParameter<double>("MinEopEndcaps");
59  eOverPMaxEndcaps_ = conf.getParameter<double>("MaxEopEndcaps");
60  dEtaMinBarrel_ = conf.getParameter<double>("MinDetaBarrel");
61  dEtaMaxBarrel_ = conf.getParameter<double>("MaxDetaBarrel");
62  dEtaMinEndcaps_ = conf.getParameter<double>("MinDetaEndcaps");
63  dEtaMaxEndcaps_ = conf.getParameter<double>("MaxDetaEndcaps");
64  dPhiMinBarrel_ = conf.getParameter<double>("MinDphiBarrel");
65  dPhiMaxBarrel_ = conf.getParameter<double>("MaxDphiBarrel");
66  dPhiMinEndcaps_ = conf.getParameter<double>("MinDphiEndcaps");
67  dPhiMaxEndcaps_ = conf.getParameter<double>("MaxDphiEndcaps");
68  sigIetaIetaMinBarrel_ = conf.getParameter<double>("MinSigIetaIetaBarrel");
69  sigIetaIetaMaxBarrel_ = conf.getParameter<double>("MaxSigIetaIetaBarrel");
70  sigIetaIetaMinEndcaps_ = conf.getParameter<double>("MinSigIetaIetaEndcaps");
71  sigIetaIetaMaxEndcaps_ = conf.getParameter<double>("MaxSigIetaIetaEndcaps");
72  hadronicOverEmMaxBarrel_ = conf.getParameter<double>("MaxHoeBarrel");
73  hadronicOverEmMaxEndcaps_ = conf.getParameter<double>("MaxHoeEndcaps");
74  mvaMin_ = conf.getParameter<double>("MinMva");
75  tipMaxBarrel_ = conf.getParameter<double>("MaxTipBarrel");
76  tipMaxEndcaps_ = conf.getParameter<double>("MaxTipEndcaps");
77  tkIso03Max_ = conf.getParameter<double>("MaxTkIso03");
78  hcalIso03Depth1MaxBarrel_ = conf.getParameter<double>("MaxHcalIso03Depth1Barrel");
79  hcalIso03Depth1MaxEndcaps_ = conf.getParameter<double>("MaxHcalIso03Depth1Endcaps");
80  hcalIso03Depth2MaxEndcaps_ = conf.getParameter<double>("MaxHcalIso03Depth2Endcaps");
81  ecalIso03MaxBarrel_ = conf.getParameter<double>("MaxEcalIso03Barrel");
82  ecalIso03MaxEndcaps_ = conf.getParameter<double>("MaxEcalIso03Endcaps");
83 
84  // for trigger
85  triggerResults_ = conf.getParameter<edm::InputTag>("TriggerResults");
86 
87  // histos limits and binning
88 // edm::ParameterSet histosSet = conf.getParameter<edm::ParameterSet>("histosCfg") ;
89 
90  nbineta=conf.getParameter<int>("NbinEta");
91  nbineta2D=conf.getParameter<int>("NbinEta2D");
92  etamin=conf.getParameter<double>("EtaMin");
93  etamax=conf.getParameter<double>("EtaMax");
94  //
95  nbinphi=conf.getParameter<int>("NbinPhi");
96  nbinphi2D=conf.getParameter<int>("NbinPhi2D");
97  phimin=conf.getParameter<double>("PhiMin");
98  phimax=conf.getParameter<double>("PhiMax");
99  //
100  nbinpt=conf.getParameter<int>("NbinPt");
101  nbinpteff=conf.getParameter<int>("NbinPtEff");
102  nbinpt2D=conf.getParameter<int>("NbinPt2D");
103  ptmax=conf.getParameter<double>("PtMax");
104  //
105  nbinp=conf.getParameter<int>("NbinP");
106  nbinp2D=conf.getParameter<int>("NbinP2D");
107  pmax=conf.getParameter<double>("PMax");
108  //
109  nbineop=conf.getParameter<int>("NbinEop");
110  nbineop2D=conf.getParameter<int>("NbinEop2D");
111  eopmax=conf.getParameter<double>("EopMax");
112  eopmaxsht=conf.getParameter<double>("EopMaxSht");
113  //
114  nbindeta=conf.getParameter<int>("NbinDeta");
115  detamin=conf.getParameter<double>("DetaMin");
116  detamax=conf.getParameter<double>("DetaMax");
117  //
118  nbindphi=conf.getParameter<int>("NbinDphi");
119  dphimin=conf.getParameter<double>("DphiMin");
120  dphimax=conf.getParameter<double>("DphiMax");
121  //
122  nbindetamatch=conf.getParameter<int>("NbinDetaMatch");
123  nbindetamatch2D=conf.getParameter<int>("NbinDetaMatch2D");
124  detamatchmin=conf.getParameter<double>("DetaMatchMin");
125  detamatchmax=conf.getParameter<double>("DetaMatchMax");
126  //
127  nbindphimatch=conf.getParameter<int>("NbinDphiMatch");
128  nbindphimatch2D=conf.getParameter<int>("NbinDphiMatch2D");
129  dphimatchmin=conf.getParameter<double>("DphiMatchMin");
130  dphimatchmax=conf.getParameter<double>("DphiMatchMax");
131  //
132  nbinfhits=conf.getParameter<int>("NbinFhits");
133  fhitsmax=conf.getParameter<double>("FhitsMax");
134  //
135  nbinlhits=conf.getParameter<int>("NbinLhits");
136  lhitsmax=conf.getParameter<double>("LhitsMax");
137  //
138  nbinxyz=conf.getParameter<int>("NbinXyz");
139  nbinxyz2D=conf.getParameter<int>("NbinXyz2D");
140  //
141  nbinpoptrue= conf.getParameter<int>("NbinPopTrue");
142  poptruemin=conf.getParameter<double>("PopTrueMin");
143  poptruemax=conf.getParameter<double>("PopTrueMax");
144  //
145  nbinmee= conf.getParameter<int>("NbinMee");
146  meemin=conf.getParameter<double>("MeeMin");
147  meemax=conf.getParameter<double>("MeeMax");
148  //
149  nbinhoe= conf.getParameter<int>("NbinHoe");
150  hoemin=conf.getParameter<double>("HoeMin");
151  hoemax=conf.getParameter<double>("HoeMax");
152 
153 // set_EfficiencyFlag=histosSet.getParameter<bool>("EfficiencyFlag");
154 // set_StatOverflowFlag=histosSet.getParameter<bool>("StatOverflowFlag");
155  }
156 
157 ElectronAnalyzer::~ElectronAnalyzer()
158  {}
159 
160 void ElectronAnalyzer::bookHistograms( DQMStore::IBooker & iBooker, edm::Run const &, edm::EventSetup const & )
161  {
162  iBooker.setCurrentFolder(outputInternalPath_) ;
163 
164  nEvents_ = 0 ;
165 
166  // basic quantities
167  h1_vertexPt_barrel = bookH1(iBooker, "vertexPt_barrel","ele transverse momentum in barrel",nbinpt,0.,ptmax,"p_{T vertex} (GeV/c)");
168  h1_vertexPt_endcaps = bookH1(iBooker, "vertexPt_endcaps","ele transverse momentum in endcaps",nbinpt,0.,ptmax,"p_{T vertex} (GeV/c)");
169  h1_vertexEta = bookH1(iBooker, "vertexEta","ele momentum #eta",nbineta,etamin,etamax,"#eta");
170  h2_vertexEtaVsPhi = bookH2(iBooker, "vertexEtaVsPhi","ele momentum #eta vs #phi",nbineta2D,etamin,etamax,nbinphi2D,phimin,phimax,"#eta","#phi (rad)");
171  h2_vertexXvsY = bookH2(iBooker, "vertexXvsY","ele vertex x vs y",nbinxyz2D,-0.1,0.1,nbinxyz2D,-0.1,0.1,"x (cm)","y (cm)");
172  h1_vertexZ = bookH1(iBooker, "vertexZ","ele vertex z",nbinxyz,-25, 25,"z (cm)");
173 
174  // super-clusters
175  h1_sclEt = bookH1(iBooker, "sclEt","ele supercluster transverse energy",nbinpt,0.,ptmax,"E_{T} (GeV)");
176 
177  // electron track
178  h1_chi2 = bookH1(iBooker, "chi2","ele track #chi^{2}",100,0.,15.,"#Chi^{2}");
179  py_chi2VsEta = bookP1(iBooker, "chi2VsEta","ele track #chi^{2} vs #eta",nbineta2D,etamin,etamax,0.,15.,"#eta","<#chi^{2}>");
180  py_chi2VsPhi = bookP1(iBooker, "chi2VsPhi","ele track #chi^{2} vs #phi",nbinphi2D,phimin,phimax,0.,15.,"#phi (rad)","<#chi^{2}>");
181  h1_foundHits = bookH1(iBooker, "foundHits","ele track # found hits",nbinfhits,0.,fhitsmax,"N_{hits}");
182  py_foundHitsVsEta = bookP1(iBooker, "foundHitsVsEta","ele track # found hits vs #eta",nbineta2D,etamin,etamax,0.,fhitsmax,"#eta","<# hits>");
183  py_foundHitsVsPhi = bookP1(iBooker, "foundHitsVsPhi","ele track # found hits vs #phi",nbinphi2D,phimin,phimax,0.,fhitsmax,"#phi (rad)","<# hits>");
184  h1_lostHits = bookH1(iBooker, "lostHits","ele track # lost hits",5,0.,5.,"N_{lost hits}");
185  py_lostHitsVsEta = bookP1(iBooker, "lostHitsVsEta","ele track # lost hits vs #eta",nbineta2D,etamin,etamax,0.,lhitsmax,"#eta","<# hits>");
186  py_lostHitsVsPhi = bookP1(iBooker, "lostHitsVsPhi","ele track # lost hits vs #eta",nbinphi2D,phimin,phimax,0.,lhitsmax,"#phi (rad)","<# hits>");
187 
188  // electron matching and ID
189  h1_Eop_barrel = bookH1(iBooker, "Eop_barrel","ele E/P_{vertex} in barrel",nbineop,0.,eopmax,"E/P_{vertex}");
190  h1_Eop_endcaps = bookH1(iBooker, "Eop_endcaps","ele E/P_{vertex} in endcaps",nbineop,0.,eopmax,"E/P_{vertex}");
191  py_EopVsPhi = bookP1(iBooker, "EopVsPhi","ele E/P_{vertex} vs #phi",nbinphi2D,phimin,phimax,0.,eopmax,"#phi (rad)","<E/P_{vertex}>");
192  h1_EeleOPout_barrel = bookH1(iBooker, "EeleOPout_barrel","ele E_{ele}/P_{out} in barrel",nbineop,0.,eopmax,"E_{ele}/P_{out}");
193  h1_EeleOPout_endcaps = bookH1(iBooker, "EeleOPout_endcaps","ele E_{ele}/P_{out} in endcaps",nbineop,0.,eopmax,"E_{ele}/P_{out}");
194  h1_dEtaSc_propVtx_barrel = bookH1(iBooker, "dEtaSc_propVtx_barrel","ele #eta_{sc} - #eta_{tr}, prop from vertex, in barrel",nbindetamatch,detamatchmin,detamatchmax,"#eta_{sc} - #eta_{tr}");
195  h1_dEtaSc_propVtx_endcapsPos = bookH1(iBooker, "dEtaSc_propVtx_endcapsPos","ele #eta_{sc} - #eta_{tr}, prop from vertex, in positive endcap",nbindetamatch,detamatchmin,detamatchmax,"#eta_{sc} - #eta_{tr}");
196  h1_dEtaSc_propVtx_endcapsNeg = bookH1(iBooker, "dEtaSc_propVtx_endcapsNeg","ele #eta_{sc} - #eta_{tr}, prop from vertex, in negative endcap",nbindetamatch,detamatchmin,detamatchmax,"#eta_{sc} - #eta_{tr}");
197  py_dEtaSc_propVtxVsPhi = bookP1(iBooker, "dEtaSc_propVtxVsPhi","ele #eta_{sc} - #eta_{tr}, prop from vertex vs #phi",nbinphi2D,phimin,phimax,detamatchmin,detamatchmax,"#phi (rad)","<#eta_{sc} - #eta_{tr}>");
198  h1_dEtaEleCl_propOut_barrel = bookH1(iBooker, "dEtaEleCl_propOut_barrel","ele #eta_{EleCl} - #eta_{tr}, prop from outermost, in barrel",nbindetamatch,detamatchmin,detamatchmax,"#eta_{elecl} - #eta_{tr}");
199  h1_dEtaEleCl_propOut_endcapsPos = bookH1(iBooker, "dEtaEleCl_propOut_endcapsPos","ele #eta_{EleCl} - #eta_{tr}, prop from outermost, in positive endcap",nbindetamatch,detamatchmin,detamatchmax,"#eta_{elecl} - #eta_{tr}");
200  h1_dEtaEleCl_propOut_endcapsNeg = bookH1(iBooker, "dEtaEleCl_propOut_endcapsNeg","ele #eta_{EleCl} - #eta_{tr}, prop from outermost, in negative endcap",nbindetamatch,detamatchmin,detamatchmax,"#eta_{elecl} - #eta_{tr}");
201  h1_dPhiSc_propVtx_barrel = bookH1(iBooker, "dPhiSc_propVtx_barrel","ele #phi_{sc} - #phi_{tr}, prop from vertex, in barrel",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{sc} - #phi_{tr} (rad)");
202  h1_dPhiSc_propVtx_endcapsPos = bookH1(iBooker, "dPhiSc_propVtx_endcapsPos","ele #phi_{sc} - #phi_{tr}, prop from vertex, in positive endcap",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{sc} - #phi_{tr} (rad)");
203  h1_dPhiSc_propVtx_endcapsNeg = bookH1(iBooker, "dPhiSc_propVtx_endcapsNeg","ele #phi_{sc} - #phi_{tr}, prop from vertex, in negative endcap",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{sc} - #phi_{tr} (rad)");
204  py_dPhiSc_propVtxVsPhi = bookP1(iBooker, "dPhiSc_propVtxVsPhi","ele #phi_{sc} - #phi_{tr}, prop from vertex vs #phi",nbinphi2D,phimin,phimax,dphimatchmin,dphimatchmax,"#phi (rad)","<#phi_{sc} - #phi_{tr}> (rad)");
205  h1_dPhiEleCl_propOut_barrel = bookH1(iBooker, "dPhiEleCl_propOut_barrel","ele #phi_{EleCl} - #phi_{tr}, prop from outermost, in barrel",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{elecl} - #phi_{tr} (rad)");
206  h1_dPhiEleCl_propOut_endcapsPos = bookH1(iBooker, "dPhiEleCl_propOut_endcapsPos","ele #phi_{EleCl} - #phi_{tr}, prop from outermost, in positive endcap",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{elecl} - #phi_{tr} (rad)");
207  h1_dPhiEleCl_propOut_endcapsNeg = bookH1(iBooker, "dPhiEleCl_propOut_endcapsNeg","ele #phi_{EleCl} - #phi_{tr}, prop from outermost, in negative endcap",nbindphimatch,dphimatchmin,dphimatchmax,"#phi_{elecl} - #phi_{tr} (rad)");
208  h1_Hoe_barrel = bookH1(iBooker, "Hoe_barrel","ele hadronic energy / em energy, in barrel", nbinhoe, hoemin, hoemax,"H/E","Events","ELE_LOGY E1 P") ;
209  h1_Hoe_endcaps = bookH1(iBooker, "Hoe_endcaps","ele hadronic energy / em energy, in endcaps", nbinhoe, hoemin, hoemax,"H/E","Events","ELE_LOGY E1 P") ;
210  py_HoeVsPhi = bookP1(iBooker, "HoeVsPhi","ele hadronic energy / em energy vs #phi",nbinphi2D,phimin,phimax,hoemin,hoemax,"#phi (rad)","<H/E>","E1 P") ;
211  h1_sclSigEtaEta_barrel = bookH1(iBooker, "sclSigEtaEta_barrel","ele sigma eta eta in barrel",100,0.,0.05,"sietaieta");
212  h1_sclSigEtaEta_endcaps = bookH1(iBooker, "sclSigEtaEta_endcaps","ele sigma eta eta in endcaps",100,0.,0.05,"sietaieta");
213  h1_sigIEtaIEta5x5_barrel = bookH1(iBooker, "sigIEtaIEta5x5_barrel","ele sigma ieta ieta 5x5 in barrel",100,0.,0.05,"sietaieta5x5");
214  h1_sigIEtaIEta5x5_endcaps = bookH1(iBooker, "sigIEtaIEta5x5_endcaps","ele sigma ieta ieta 5x5 in endcaps",100,0.,0.05,"sietaieta5x5");
215 
216  // fbrem
217  h1_fbrem = bookH1(iBooker, "fbrem","ele brem fraction",100,0.,1.,"P_{in} - P_{out} / P_{in}") ;
218  py_fbremVsEta = bookP1(iBooker, "fbremVsEta","ele brem fraction vs #eta",nbineta2D,etamin,etamax,0.,1.,"#eta","<P_{in} - P_{out} / P_{in}>") ;
219  py_fbremVsPhi = bookP1(iBooker, "fbremVsPhi","ele brem fraction vs #phi",nbinphi2D,phimin,phimax,0.,1.,"#phi (rad)","<P_{in} - P_{out} / P_{in}>") ;
220  h1_classes = bookH1(iBooker, "classes","ele electron classes",10,0.0,10.);
221 
222  // pflow
223  h1_mva = bookH1(iBooker, "mva","ele identification mva",100,-1.,1.,"mva");
224  h1_provenance = bookH1(iBooker, "provenance","ele provenance",5,-2.,3.,"provenance");
225 
226  // isolation
227  h1_tkSumPt_dr03 = bookH1(iBooker, "tkSumPt_dr03","tk isolation sum, dR=0.3",100,0.0,20.,"TkIsoSum (GeV/c)","Events","ELE_LOGY E1 P");
228  h1_ecalRecHitSumEt_dr03 = bookH1(iBooker, "ecalRecHitSumEt_dr03","ecal isolation sum, dR=0.3",100,0.0,20.,"EcalIsoSum (GeV)","Events","ELE_LOGY E1 P");
229  h1_hcalTowerSumEt_dr03 = bookH1(iBooker, "hcalTowerSumEt_dr03","hcal isolation sum, dR=0.3",100,0.0,20.,"HcalIsoSum (GeV)","Events","ELE_LOGY E1 P");
230 
231  // pf isolation
232  h1_PFch_dr03 = bookH1(iBooker, "PFch_dr03","Charged PF candidate sum, dR=0.3",100,0.0,20.,"PF charged (GeV/c)","Events","ELE_LOGY E1 P");
233  h1_PFem_dr03 = bookH1(iBooker, "PFem_dr03","Neutral EM PF candidate sum, dR=0.3",100,0.0,20.,"PF neutral EM (GeV)","Events","ELE_LOGY E1 P");
234  h1_PFnh_dr03 = bookH1(iBooker, "PFnh_dr03","Neutral Had PF candidate sum, dR=0.3",100,0.0,20.,"PF neutral Had (GeV)","Events","ELE_LOGY E1 P");
235 
236  // di-electron mass
237  setBookIndex(200) ;
238  h1_mee = bookH1(iBooker, "mee","ele pairs invariant mass", nbinmee, meemin, meemax,"m_{ee} (GeV/c^{2})");
239  h1_mee_os = bookH1(iBooker, "mee_os","ele pairs invariant mass, opposite sign", nbinmee, meemin, meemax,"m_{e^{+}e^{-}} (GeV/c^{2})");
240 
241  //===========================
242  // histos for matching and matched matched objects
243  //===========================
244 
245  // matching object
246  std::string matchingObjectType ;
247  Labels l;
248  labelsForToken(matchingObjectCollection_,l);
249  if (std::string::npos != std::string(l.module).find("SuperCluster",0))
250  { matchingObjectType = "SC" ; }
251  if (matchingObjectType=="")
252  { edm::LogError("ElectronMcFakeValidator::beginJob")<<"Unknown matching object type !" ; }
253  else
254  { edm::LogInfo("ElectronMcFakeValidator::beginJob")<<"Matching object type: "<<matchingObjectType ; }
255 
256  // matching object distributions
257  h1_matchingObject_Eta = bookH1withSumw2(iBooker, "matchingObject_Eta",matchingObjectType+" #eta",nbineta,etamin,etamax,"#eta_{SC}");
258 
259 
260  h1_matchingObject_Pt = bookH1withSumw2(iBooker, "matchingObject_Pt",matchingObjectType+" pt",nbinpteff,5.,ptmax,"pt_{SC} (GeV/c)");
261  h1_matchingObject_Phi = bookH1withSumw2(iBooker, "matchingObject_Phi",matchingObjectType+" #phi",nbinphi,phimin,phimax,"#phi (rad)");
262 
263  h1_matchedObject_Eta = bookH1withSumw2(iBooker, "matchedObject_Eta","Efficiency vs matching SC #eta",nbineta,etamin,etamax,"#eta_{SC}");
264  h1_matchedObject_Pt = bookH1withSumw2(iBooker, "matchedObject_Pt","Efficiency vs matching SC E_{T}",nbinpteff,5.,ptmax,"pt_{SC} (GeV/c)");
265  h1_matchedObject_Phi = bookH1withSumw2(iBooker, "matchedObject_Phi","Efficiency vs matching SC #phi",nbinphi,phimin,phimax,"#phi (rad)");
266 
267 
268  }
269 
270 void ElectronAnalyzer::analyze( const edm::Event& iEvent, const edm::EventSetup & iSetup )
271 {
272  nEvents_++ ;
273 
274  edm::Handle<GsfElectronCollection> gsfElectrons ;
275  iEvent.getByToken(electronCollection_,gsfElectrons) ;
276  edm::Handle<reco::SuperClusterCollection> recoClusters ;
277  iEvent.getByToken(matchingObjectCollection_,recoClusters) ;
278  edm::Handle<reco::TrackCollection> tracks;
279  iEvent.getByToken(trackCollection_,tracks);
280  edm::Handle<reco::GsfTrackCollection> gsfTracks;
281  iEvent.getByToken(gsftrackCollection_,gsfTracks);
282  edm::Handle<reco::VertexCollection> vertices;
283  iEvent.getByToken(vertexCollection_,vertices);
284  edm::Handle<reco::BeamSpot> recoBeamSpotHandle ;
285  iEvent.getByToken(beamSpotTag_,recoBeamSpotHandle) ;
286  const BeamSpot bs = *recoBeamSpotHandle ;
287 
288  edm::EventNumber_t ievt = iEvent.id().event();
289  edm::RunNumber_t irun = iEvent.id().run();
290  edm::LuminosityBlockNumber_t ils = iEvent.luminosityBlock();
291 
292  edm::LogInfo("ElectronAnalyzer::analyze")
293  <<"Treating "<<gsfElectrons.product()->size()<<" electrons"
294  <<" from event "<<ievt<<" in run "<<irun<<" and lumiblock "<<ils ;
295  //h1_num_->Fill((*gsfElectrons).size()) ;
296 
297  // selected rec electrons
298  reco::GsfElectronCollection::const_iterator gsfIter ;
299  for
300  ( gsfIter=gsfElectrons->begin() ;
301  gsfIter!=gsfElectrons->end();
302  gsfIter++ )
303  {
304  // vertex TIP
305  double vertexTIP =
306  (gsfIter->vertex().x()-bs.position().x()) * (gsfIter->vertex().x()-bs.position().x()) +
307  (gsfIter->vertex().y()-bs.position().y()) * (gsfIter->vertex().y()-bs.position().y()) ;
308  vertexTIP = sqrt(vertexTIP) ;
309 
310  // select electrons
311  if (!selected(gsfIter,vertexTIP)) continue ;
312 
313  // invariant mass of 2 selected electrons
314  reco::GsfElectronCollection::const_iterator gsfIter2 ;
315  for
316  ( gsfIter2=gsfIter+1;
317  gsfIter2!=gsfElectrons->end() ;
318  gsfIter2++ )
319  {
320  if (!selected(gsfIter2,vertexTIP)) continue ;
321  float invMass = computeInvMass(*gsfIter,*gsfIter2) ;
322  h1_mee->Fill(invMass) ;
323  if ( ( (gsfIter->charge())*(gsfIter2->charge()) )<0. )
324  { h1_mee_os->Fill(invMass) ; }
325  }
326 
327  // basic quantities
328  if (gsfIter->isEB()) h1_vertexPt_barrel->Fill( gsfIter->pt() );
329  if (gsfIter->isEE()) h1_vertexPt_endcaps->Fill( gsfIter->pt() );
330  h1_vertexEta->Fill( gsfIter->eta() );
331  h2_vertexEtaVsPhi->Fill( gsfIter->eta(), gsfIter->phi() );
332  h2_vertexXvsY->Fill( gsfIter->vertex().x(), gsfIter->vertex().y() );
333  h1_vertexZ->Fill( gsfIter->vertex().z() );
334 
335  // supercluster related distributions
336  reco::SuperClusterRef sclRef = gsfIter->superCluster() ;
337  // ALREADY DONE IN GSF ELECTRON CORE
338  double R=TMath::Sqrt(sclRef->x()*sclRef->x() + sclRef->y()*sclRef->y() +sclRef->z()*sclRef->z());
339  double Rt=TMath::Sqrt(sclRef->x()*sclRef->x() + sclRef->y()*sclRef->y());
340  h1_sclEt->Fill(sclRef->energy()*(Rt/R));
341 
342  // track related distributions
343  if (!readAOD_)
344  { // track extra does not exist in AOD
345  h1_foundHits->Fill( gsfIter->gsfTrack()->numberOfValidHits() );
346  py_foundHitsVsEta->Fill( gsfIter->eta(), gsfIter->gsfTrack()->numberOfValidHits() );
347  py_foundHitsVsPhi->Fill( gsfIter->phi(), gsfIter->gsfTrack()->numberOfValidHits() );
348  h1_lostHits->Fill( gsfIter->gsfTrack()->numberOfLostHits() );
349  py_lostHitsVsEta->Fill( gsfIter->eta(), gsfIter->gsfTrack()->numberOfLostHits() );
350  py_lostHitsVsPhi->Fill( gsfIter->phi(), gsfIter->gsfTrack()->numberOfLostHits() );
351  h1_chi2->Fill( gsfIter->gsfTrack()->normalizedChi2() );
352  py_chi2VsEta->Fill( gsfIter->eta(), gsfIter->gsfTrack()->normalizedChi2() );
353  py_chi2VsPhi->Fill( gsfIter->phi(), gsfIter->gsfTrack()->normalizedChi2() );
354  }
355 
356  // match distributions
357  if (gsfIter->isEB())
358  {
359  h1_Eop_barrel->Fill( gsfIter->eSuperClusterOverP() );
360  h1_EeleOPout_barrel->Fill( gsfIter->eEleClusterOverPout() );
361  h1_dEtaSc_propVtx_barrel->Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
362  h1_dEtaEleCl_propOut_barrel->Fill(gsfIter->deltaEtaEleClusterTrackAtCalo());
363  h1_dPhiSc_propVtx_barrel->Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
364  h1_dPhiEleCl_propOut_barrel->Fill(gsfIter->deltaPhiEleClusterTrackAtCalo());
365  h1_Hoe_barrel->Fill(gsfIter->hadronicOverEm());
366  h1_sclSigEtaEta_barrel->Fill( gsfIter->scSigmaEtaEta() );
367  h1_sigIEtaIEta5x5_barrel->Fill( gsfIter->full5x5_sigmaIetaIeta() );
368  }
369  if (gsfIter->isEE())
370  {
371  h1_Eop_endcaps->Fill( gsfIter->eSuperClusterOverP() );
372  h1_EeleOPout_endcaps->Fill( gsfIter->eEleClusterOverPout() );
373  if (gsfIter->eta() > 0) {
374  h1_dEtaSc_propVtx_endcapsPos->Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
375  h1_dEtaEleCl_propOut_endcapsPos->Fill(gsfIter->deltaEtaEleClusterTrackAtCalo());
376  h1_dPhiSc_propVtx_endcapsPos->Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
377  h1_dPhiEleCl_propOut_endcapsPos->Fill(gsfIter->deltaPhiEleClusterTrackAtCalo());
378  }
379  else {
380  h1_dEtaSc_propVtx_endcapsNeg->Fill(gsfIter->deltaEtaSuperClusterTrackAtVtx());
381  h1_dEtaEleCl_propOut_endcapsNeg->Fill(gsfIter->deltaEtaEleClusterTrackAtCalo());
382  h1_dPhiSc_propVtx_endcapsNeg->Fill(gsfIter->deltaPhiSuperClusterTrackAtVtx());
383  h1_dPhiEleCl_propOut_endcapsNeg->Fill(gsfIter->deltaPhiEleClusterTrackAtCalo());
384  }
385  h1_Hoe_endcaps->Fill(gsfIter->hadronicOverEm());
386  h1_sclSigEtaEta_endcaps->Fill( gsfIter->scSigmaEtaEta() );
387  h1_sigIEtaIEta5x5_endcaps->Fill( gsfIter->full5x5_sigmaIetaIeta() );
388  }
389  py_EopVsPhi->Fill( gsfIter->phi(), gsfIter->eSuperClusterOverP() );
390  py_dEtaSc_propVtxVsPhi->Fill(gsfIter->phi(), gsfIter->deltaEtaSuperClusterTrackAtVtx());
391  py_dPhiSc_propVtxVsPhi->Fill(gsfIter->phi(), gsfIter->deltaPhiSuperClusterTrackAtVtx());
392  py_HoeVsPhi->Fill(gsfIter->phi(), gsfIter->hadronicOverEm());
393 
394  // fbrem, classes
395  h1_fbrem->Fill(gsfIter->fbrem()) ;
396  py_fbremVsEta->Fill(gsfIter->eta(),gsfIter->fbrem()) ;
397  py_fbremVsPhi->Fill(gsfIter->phi(),gsfIter->fbrem()) ;
398  int eleClass = gsfIter->classification() ;
399  if (gsfIter->isEE()) eleClass+=5;
400  h1_classes->Fill(eleClass) ;
401 
402 
403  // pflow
404  h1_mva->Fill(gsfIter->mva_e_pi()) ;
405  if (gsfIter->ecalDrivenSeed()) h1_provenance->Fill(1.) ;
406  if (gsfIter->trackerDrivenSeed()) h1_provenance->Fill(-1.) ;
407  if (gsfIter->trackerDrivenSeed()||gsfIter->ecalDrivenSeed()) h1_provenance->Fill(0.);
408  if (gsfIter->trackerDrivenSeed()&&!gsfIter->ecalDrivenSeed()) h1_provenance->Fill(-2.);
409  if (!gsfIter->trackerDrivenSeed()&&gsfIter->ecalDrivenSeed()) h1_provenance->Fill(2.);
410 
411  // isolation
412  h1_tkSumPt_dr03->Fill(gsfIter->dr03TkSumPt());
413  h1_ecalRecHitSumEt_dr03->Fill(gsfIter->dr03EcalRecHitSumEt());
414  h1_hcalTowerSumEt_dr03->Fill(gsfIter->dr03HcalTowerSumEt());
415 
416  // PF isolation
417  GsfElectron::PflowIsolationVariables pfIso = gsfIter->pfIsolationVariables();
418  h1_PFch_dr03->Fill( pfIso.sumChargedHadronPt );
419  h1_PFem_dr03->Fill( pfIso.sumPhotonEt );
420  h1_PFnh_dr03->Fill( pfIso.sumNeutralHadronEt );
421 
422 
423  }
424 
425  // association matching object-reco electrons
426  int matchingObjectNum=0;
427  reco::SuperClusterCollection::const_iterator moIter ;
428  for
429  ( moIter=recoClusters->begin() ;
430  moIter!=recoClusters->end() ;
431  moIter++ )
432  {
433 // // number of matching objects
434  matchingObjectNum++;
435 
436  if
437  ( moIter->energy()/cosh(moIter->eta())>maxPtMatchingObject_ ||
438  std::abs(moIter->eta())> maxAbsEtaMatchingObject_ )
439  { continue ; }
440 
441 // // suppress the endcaps
442  h1_matchingObject_Eta->Fill( moIter->eta() );
443  h1_matchingObject_Pt->Fill( moIter->energy()/cosh(moIter->eta()) );
444  h1_matchingObject_Phi->Fill( moIter->phi() );
445 
446  bool okGsfFound = false ;
447  double gsfOkRatio = 999999999. ;
448  reco::GsfElectron bestGsfElectron ;
449  reco::GsfElectronCollection::const_iterator gsfIter ;
450  for
451  ( gsfIter=gsfElectrons->begin() ;
452  gsfIter!=gsfElectrons->end() ;
453  gsfIter++ )
454  {
455 
456  double vertexTIP =
457  (gsfIter->vertex().x()-bs.position().x()) * (gsfIter->vertex().x()-bs.position().x()) +
458  (gsfIter->vertex().y()-bs.position().y()) * (gsfIter->vertex().y()-bs.position().y()) ;
459  vertexTIP = sqrt(vertexTIP) ;
460 
461  // select electrons
462  if (!selected(gsfIter,vertexTIP)) continue ;
463 
464  // matching with a cone in eta phi
465  if ( matchingCondition_ == "Cone" )
466  {
467  double dphi = gsfIter->phi()-moIter->phi() ;
468  if (std::abs(dphi)>CLHEP::pi)
469  { dphi = dphi < 0? (CLHEP::twopi) + dphi : dphi - CLHEP::twopi ; }
470  double deltaR = sqrt(pow((moIter->eta()-gsfIter->eta()),2) + pow(dphi,2)) ;
471  if ( deltaR < deltaR_ )
472  {
473  double tmpGsfRatio = gsfIter->p()/moIter->energy() ;
474  if ( std::abs(tmpGsfRatio-1) < std::abs(gsfOkRatio-1) )
475  {
476  gsfOkRatio = tmpGsfRatio;
477  bestGsfElectron=*gsfIter;
478  okGsfFound = true;
479  }
480  }
481  }
482  } // loop over rec ele to look for the best one
483  if (okGsfFound)
484  {
485  // generated distributions for matched electrons
486  h1_matchedObject_Eta->Fill( moIter->eta() );
487  h1_matchedObject_Pt->Fill( moIter->energy()/cosh(moIter->eta()) );
488  h1_matchedObject_Phi->Fill( moIter->phi() );
489 
490  //classes
491  }
492 
493  } // loop overmatching object
494 
495 
496  }
497 
498 float ElectronAnalyzer::computeInvMass
499  ( const reco::GsfElectron & e1,
500  const reco::GsfElectron & e2 )
501  {
502  math::XYZTLorentzVector p12 = e1.p4()+e2.p4() ;
503  float mee2 = p12.Dot(p12) ;
504  float invMass = mee2 > 0. ? sqrt(mee2) : 0;
505  return invMass ;
506  }
507 
508 bool ElectronAnalyzer::selected( const reco::GsfElectronCollection::const_iterator & gsfIter , double vertexTIP )
509  {
510  if ((Selection_>0)&&generalCut(gsfIter)) return false ;
511  if ((Selection_>=1)&&etCut(gsfIter)) return false ;
512  if ((Selection_>=2)&&isolationCut(gsfIter,vertexTIP)) return false ;
513  if ((Selection_>=3)&&idCut(gsfIter)) return false ;
514  return true ;
515  }
516 
517 bool ElectronAnalyzer::generalCut( const reco::GsfElectronCollection::const_iterator & gsfIter)
518  {
519  if (std::abs(gsfIter->eta())>maxAbsEta_) return true ;
520  if (gsfIter->pt()<minPt_) return true ;
521 
522  if (gsfIter->isEB() && isEE_) return true ;
523  if (gsfIter->isEE() && isEB_) return true ;
524  if (gsfIter->isEBEEGap() && isNotEBEEGap_) return true ;
525 
526  if (gsfIter->ecalDrivenSeed() && isTrackerDriven_) return true ;
527  if (gsfIter->trackerDrivenSeed() && isEcalDriven_) return true ;
528 
529  return false ;
530  }
531 
532 bool ElectronAnalyzer::etCut( const reco::GsfElectronCollection::const_iterator & gsfIter )
533  {
534  if (gsfIter->superCluster()->energy()/cosh(gsfIter->superCluster()->eta())<minEt_) return true ;
535 
536  return false ;
537  }
538 
539 bool ElectronAnalyzer::isolationCut( const reco::GsfElectronCollection::const_iterator & gsfIter, double vertexTIP )
540  {
541  if (gsfIter->isEB() && vertexTIP > tipMaxBarrel_) return true ;
542  if (gsfIter->isEE() && vertexTIP > tipMaxEndcaps_) return true ;
543 
544  if (gsfIter->dr03TkSumPt() > tkIso03Max_) return true ;
545  if (gsfIter->isEB() && gsfIter->dr03HcalDepth1TowerSumEt() > hcalIso03Depth1MaxBarrel_) return true ;
546  if (gsfIter->isEE() && gsfIter->dr03HcalDepth1TowerSumEt() > hcalIso03Depth1MaxEndcaps_) return true ;
547  if (gsfIter->isEE() && gsfIter->dr03HcalDepth2TowerSumEt() > hcalIso03Depth2MaxEndcaps_) return true ;
548  if (gsfIter->isEB() && gsfIter->dr03EcalRecHitSumEt() > ecalIso03MaxBarrel_) return true ;
549  if (gsfIter->isEE() && gsfIter->dr03EcalRecHitSumEt() > ecalIso03MaxEndcaps_) return true ;
550 
551  return false ;
552  }
553 
554 bool ElectronAnalyzer::idCut( const reco::GsfElectronCollection::const_iterator & gsfIter )
555  {
556  if (gsfIter->isEB() && gsfIter->eSuperClusterOverP() < eOverPMinBarrel_) return true ;
557  if (gsfIter->isEB() && gsfIter->eSuperClusterOverP() > eOverPMaxBarrel_) return true ;
558  if (gsfIter->isEE() && gsfIter->eSuperClusterOverP() < eOverPMinEndcaps_) return true ;
559  if (gsfIter->isEE() && gsfIter->eSuperClusterOverP() > eOverPMaxEndcaps_) return true ;
560  if (gsfIter->isEB() && std::abs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) < dEtaMinBarrel_) return true ;
561  if (gsfIter->isEB() && std::abs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) > dEtaMaxBarrel_) return true ;
562  if (gsfIter->isEE() && std::abs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) < dEtaMinEndcaps_) return true ;
563  if (gsfIter->isEE() && std::abs(gsfIter->deltaEtaSuperClusterTrackAtVtx()) > dEtaMaxEndcaps_) return true ;
564  if (gsfIter->isEB() && std::abs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) < dPhiMinBarrel_) return true ;
565  if (gsfIter->isEB() && std::abs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) > dPhiMaxBarrel_) return true ;
566  if (gsfIter->isEE() && std::abs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) < dPhiMinEndcaps_) return true ;
567  if (gsfIter->isEE() && std::abs(gsfIter->deltaPhiSuperClusterTrackAtVtx()) > dPhiMaxEndcaps_) return true ;
568  if (gsfIter->isEB() && gsfIter->scSigmaIEtaIEta() < sigIetaIetaMinBarrel_) return true ;
569  if (gsfIter->isEB() && gsfIter->scSigmaIEtaIEta() > sigIetaIetaMaxBarrel_) return true ;
570  if (gsfIter->isEE() && gsfIter->scSigmaIEtaIEta() < sigIetaIetaMinEndcaps_) return true ;
571  if (gsfIter->isEE() && gsfIter->scSigmaIEtaIEta() > sigIetaIetaMaxEndcaps_) return true ;
572  if (gsfIter->isEB() && gsfIter->hadronicOverEm() > hadronicOverEmMaxBarrel_) return true ;
573  if (gsfIter->isEE() && gsfIter->hadronicOverEm() > hadronicOverEmMaxEndcaps_) return true ;
574 
575  return false ;
576  }
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:39
ElectronAnalyzer::h1_foundHits
MonitorElement * h1_foundHits
Definition: ElectronAnalyzer.h:204
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:41
ElectronAnalyzer::h1_dEtaSc_propVtx_endcapsPos
MonitorElement * h1_dEtaSc_propVtx_endcapsPos
Definition: ElectronAnalyzer.h:232
ElectronAnalyzer::hadronicOverEmMaxEndcaps_
double hadronicOverEmMaxEndcaps_
Definition: ElectronAnalyzer.h:100
ElectronAnalyzer::h1_mee_os
MonitorElement * h1_mee_os
Definition: ElectronAnalyzer.h:303
reco::GsfElectron::p4
const LorentzVector & p4(P4Kind kind) const
Definition: GsfElectron.cc:223
ElectronAnalyzer::selected
bool selected(const reco::GsfElectronCollection::const_iterator &gsfIter, double vertexTIP)
Definition: ElectronAnalyzer.cc:508
ElectronAnalyzer::h1_dPhiSc_propVtx_barrel
MonitorElement * h1_dPhiSc_propVtx_barrel
Definition: ElectronAnalyzer.h:245
ElectronAnalyzer::nEvents_
unsigned int nEvents_
Definition: ElectronAnalyzer.h:137
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
ElectronAnalyzer::h1_vertexPt_barrel
MonitorElement * h1_vertexPt_barrel
Definition: ElectronAnalyzer.h:182
ElectronAnalyzer::h1_PFch_dr03
MonitorElement * h1_PFch_dr03
Definition: ElectronAnalyzer.h:290
ElectronAnalyzer::py_lostHitsVsEta
MonitorElement * py_lostHitsVsEta
Definition: ElectronAnalyzer.h:209
ElectronAnalyzer::h1_sclSigEtaEta_barrel
MonitorElement * h1_sclSigEtaEta_barrel
Definition: ElectronAnalyzer.h:263
assert
assert(m_qm.get())
ElectronAnalyzer::py_chi2VsEta
MonitorElement * py_chi2VsEta
Definition: ElectronAnalyzer.h:201
ElectronAnalyzer::matchingObjectCollection_
edm::EDGetTokenT< reco::SuperClusterCollection > matchingObjectCollection_
Definition: ElectronAnalyzer.h:56
ElectronAnalyzer::h1_vertexEta
MonitorElement * h1_vertexEta
Definition: ElectronAnalyzer.h:184
ElectronAnalyzer::py_chi2VsPhi
MonitorElement * py_chi2VsPhi
Definition: ElectronAnalyzer.h:202
ElectronAnalyzer::h1_Eop_endcaps
MonitorElement * h1_Eop_endcaps
Definition: ElectronAnalyzer.h:220
ElectronDqmAnalyzerBase::bookH1withSumw2
MonitorElement * bookH1withSumw2(DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")
Definition: ElectronDqmAnalyzerBase.cc:110
edm::EventNumber_t
unsigned long long EventNumber_t
Definition: RunLumiEventNumber.h:12
ElectronAnalyzer::h1_ecalRecHitSumEt_dr03
MonitorElement * h1_ecalRecHitSumEt_dr03
Definition: ElectronAnalyzer.h:288
edm::EventBase::luminosityBlock
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:63
ElectronAnalyzer::h1_hcalTowerSumEt_dr03
MonitorElement * h1_hcalTowerSumEt_dr03
Definition: ElectronAnalyzer.h:289
ElectronAnalyzer::h1_dPhiEleCl_propOut_endcapsNeg
MonitorElement * h1_dPhiEleCl_propOut_endcapsNeg
Definition: ElectronAnalyzer.h:254
ElectronAnalyzer::h1_matchingObject_Eta
MonitorElement * h1_matchingObject_Eta
Definition: ElectronAnalyzer.h:321
ElectronAnalyzer::h1_Hoe_barrel
MonitorElement * h1_Hoe_barrel
Definition: ElectronAnalyzer.h:259
ElectronDqmAnalyzerBase::bookH2
MonitorElement * bookH2(DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="COLZ")
Definition: ElectronDqmAnalyzerBase.cc:126
ElectronAnalyzer::triggerResults_
edm::InputTag triggerResults_
Definition: ElectronAnalyzer.h:112
ElectronAnalyzer::py_HoeVsPhi
MonitorElement * py_HoeVsPhi
Definition: ElectronAnalyzer.h:261
ElectronAnalyzer::h1_matchedObject_Pt
MonitorElement * h1_matchedObject_Pt
Definition: ElectronAnalyzer.h:316
ElectronAnalyzer::sigIetaIetaMinBarrel_
double sigIetaIetaMinBarrel_
Definition: ElectronAnalyzer.h:95
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:106
ElectronAnalyzer::bookHistograms
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
Definition: ElectronAnalyzer.cc:160
edm::LuminosityBlockNumber_t
unsigned int LuminosityBlockNumber_t
Definition: RunLumiEventNumber.h:13
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
ElectronAnalyzer::h1_matchedObject_Phi
MonitorElement * h1_matchedObject_Phi
Definition: ElectronAnalyzer.h:317
ElectronAnalyzer::h1_PFem_dr03
MonitorElement * h1_PFem_dr03
Definition: ElectronAnalyzer.h:291
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
pi
const Double_t pi
Definition: trackSplitPlot.h:36
ElectronAnalyzer::h1_tkSumPt_dr03
MonitorElement * h1_tkSumPt_dr03
Definition: ElectronAnalyzer.h:287
ElectronAnalyzer::h1_dEtaEleCl_propOut_endcapsPos
MonitorElement * h1_dEtaEleCl_propOut_endcapsPos
Definition: ElectronAnalyzer.h:239
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
ElectronAnalyzer::sigIetaIetaMaxBarrel_
double sigIetaIetaMaxBarrel_
Definition: ElectronAnalyzer.h:96
iEvent
int iEvent
Definition: GenABIO.cc:230
ElectronAnalyzer::h2_vertexXvsY
MonitorElement * h2_vertexXvsY
Definition: ElectronAnalyzer.h:186
ElectronAnalyzer::h1_dPhiSc_propVtx_endcapsPos
MonitorElement * h1_dPhiSc_propVtx_endcapsPos
Definition: ElectronAnalyzer.h:246
ElectronAnalyzer::py_fbremVsEta
MonitorElement * py_fbremVsEta
Definition: ElectronAnalyzer.h:277
ElectronAnalyzer::h1_provenance
MonitorElement * h1_provenance
Definition: ElectronAnalyzer.h:284
ElectronAnalyzer::h1_dEtaEleCl_propOut_barrel
MonitorElement * h1_dEtaEleCl_propOut_barrel
Definition: ElectronAnalyzer.h:238
ElectronAnalyzer::py_EopVsPhi
MonitorElement * py_EopVsPhi
Definition: ElectronAnalyzer.h:221
ElectronAnalyzer::h1_PFnh_dr03
MonitorElement * h1_PFnh_dr03
Definition: ElectronAnalyzer.h:292
ElectronAnalyzer::h1_Eop_barrel
MonitorElement * h1_Eop_barrel
Definition: ElectronAnalyzer.h:219
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
reco::GsfElectron::PflowIsolationVariables::sumPhotonEt
float sumPhotonEt
sum pt of PF photons // old float photonIso ;
Definition: GsfElectron.h:561
ElectronAnalyzer::sigIetaIetaMinEndcaps_
double sigIetaIetaMinEndcaps_
Definition: ElectronAnalyzer.h:97
ElectronAnalyzer::py_dEtaSc_propVtxVsPhi
MonitorElement * py_dEtaSc_propVtxVsPhi
Definition: ElectronAnalyzer.h:234
ElectronAnalyzer::sigIetaIetaMaxEndcaps_
double sigIetaIetaMaxEndcaps_
Definition: ElectronAnalyzer.h:98
HLT_25ns14e33_v1_cff.vertices
tuple vertices
Definition: HLT_25ns14e33_v1_cff.py:4396
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
ElectronAnalyzer::matchingCondition_
std::string matchingCondition_
Definition: ElectronAnalyzer.h:69
reco::GsfElectron::PflowIsolationVariables::sumNeutralHadronEt
float sumNeutralHadronEt
sum pt of neutral hadrons // old float neutralHadronIso ;
Definition: GsfElectron.h:560
ElectronAnalyzer::h1_dPhiEleCl_propOut_endcapsPos
MonitorElement * h1_dPhiEleCl_propOut_endcapsPos
Definition: ElectronAnalyzer.h:253
ElectronAnalyzer::hcalIso03Depth1MaxBarrel_
double hcalIso03Depth1MaxBarrel_
Definition: ElectronAnalyzer.h:105
ElectronAnalyzer::h1_Hoe_endcaps
MonitorElement * h1_Hoe_endcaps
Definition: ElectronAnalyzer.h:260
ElectronDqmAnalyzerBase::bookP1
MonitorElement * bookP1(DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, double lowY, double highY, const std::string &titleX="", const std::string &titleY="", Option_t *option="E1 P")
Definition: ElectronDqmAnalyzerBase.cc:159
ElectronAnalyzer::h1_dEtaEleCl_propOut_endcapsNeg
MonitorElement * h1_dEtaEleCl_propOut_endcapsNeg
Definition: ElectronAnalyzer.h:240
dbtoconf.conf
tuple conf
Definition: dbtoconf.py:185
ElectronAnalyzer::h1_dPhiSc_propVtx_endcapsNeg
MonitorElement * h1_dPhiSc_propVtx_endcapsNeg
Definition: ElectronAnalyzer.h:247
ElectronAnalyzer::ElectronAnalyzer
ElectronAnalyzer(const edm::ParameterSet &conf)
Definition: ElectronAnalyzer.cc:27
ElectronAnalyzer::gsftrackCollection_
edm::EDGetTokenT< reco::GsfTrackCollection > gsftrackCollection_
Definition: ElectronAnalyzer.h:57
ElectronAnalyzer::h1_vertexPt_endcaps
MonitorElement * h1_vertexPt_endcaps
Definition: ElectronAnalyzer.h:183
ElectronAnalyzer::py_lostHitsVsPhi
MonitorElement * py_lostHitsVsPhi
Definition: ElectronAnalyzer.h:210
deltaR
double deltaR(double eta1, double eta2, double phi1, double phi2)
Definition: TreeUtility.cc:17
ElectronAnalyzer::h1_vertexZ
MonitorElement * h1_vertexZ
Definition: ElectronAnalyzer.h:187
ElectronAnalyzer::py_dPhiSc_propVtxVsPhi
MonitorElement * py_dPhiSc_propVtxVsPhi
Definition: ElectronAnalyzer.h:248
ElectronAnalyzer::hcalIso03Depth1MaxEndcaps_
double hcalIso03Depth1MaxEndcaps_
Definition: ElectronAnalyzer.h:106
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
ElectronAnalyzer::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: ElectronAnalyzer.cc:270
ElectronAnalyzer::py_fbremVsPhi
MonitorElement * py_fbremVsPhi
Definition: ElectronAnalyzer.h:278
edm::Handle::product
T const * product() const
Definition: Handle.h:81
ElectronAnalyzer::h1_dEtaSc_propVtx_endcapsNeg
MonitorElement * h1_dEtaSc_propVtx_endcapsNeg
Definition: ElectronAnalyzer.h:233
ElectronAnalyzer::py_foundHitsVsPhi
MonitorElement * py_foundHitsVsPhi
Definition: ElectronAnalyzer.h:206
ElectronAnalyzer::etCut
bool etCut(const reco::GsfElectronCollection::const_iterator &gsfIter)
Definition: ElectronAnalyzer.cc:532
testEve_cfg.tracks
tuple tracks
Definition: testEve_cfg.py:39
ElectronAnalyzer::h1_fbrem
MonitorElement * h1_fbrem
Definition: ElectronAnalyzer.h:276
ElectronAnalyzer::trackCollection_
edm::EDGetTokenT< reco::TrackCollection > trackCollection_
Definition: ElectronAnalyzer.h:58
ElectronAnalyzer::h1_chi2
MonitorElement * h1_chi2
Definition: ElectronAnalyzer.h:200
ElectronAnalyzer::maxAbsEtaMatchingObject_
double maxAbsEtaMatchingObject_
Definition: ElectronAnalyzer.h:71
ElectronAnalyzer::h1_lostHits
MonitorElement * h1_lostHits
Definition: ElectronAnalyzer.h:208
ElectronAnalyzer::h2_vertexEtaVsPhi
MonitorElement * h2_vertexEtaVsPhi
Definition: ElectronAnalyzer.h:185
ElectronAnalyzer::h1_EeleOPout_endcaps
MonitorElement * h1_EeleOPout_endcaps
Definition: ElectronAnalyzer.h:226
ElectronAnalyzer::h1_mee
MonitorElement * h1_mee
Definition: ElectronAnalyzer.h:302
ElectronAnalyzer::electronCollection_
edm::EDGetTokenT< reco::GsfElectronCollection > electronCollection_
Definition: ElectronAnalyzer.h:55
ElectronAnalyzer::h1_classes
MonitorElement * h1_classes
Definition: ElectronAnalyzer.h:280
ElectronAnalyzer::h1_sclEt
MonitorElement * h1_sclEt
Definition: ElectronAnalyzer.h:193
edm::EDConsumerBase::labelsForToken
void labelsForToken(EDGetToken iToken, Labels &oLabels) const
Definition: EDConsumerBase.cc:288
ElectronAnalyzer::h1_sigIEtaIEta5x5_barrel
MonitorElement * h1_sigIEtaIEta5x5_barrel
Definition: ElectronAnalyzer.h:265
ElectronAnalyzer::h1_EeleOPout_barrel
MonitorElement * h1_EeleOPout_barrel
Definition: ElectronAnalyzer.h:225
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:60
ElectronAnalyzer::outputInternalPath_
std::string outputInternalPath_
Definition: ElectronAnalyzer.h:66
ElectronAnalyzer::generalCut
bool generalCut(const reco::GsfElectronCollection::const_iterator &gsfIter)
Definition: ElectronAnalyzer.cc:517
ElectronAnalyzer::h1_dEtaSc_propVtx_barrel
MonitorElement * h1_dEtaSc_propVtx_barrel
Definition: ElectronAnalyzer.h:231
ElectronAnalyzer::computeInvMass
float computeInvMass(const reco::GsfElectron &e1, const reco::GsfElectron &e2)
Definition: ElectronAnalyzer.cc:499
ElectronAnalyzer::h1_matchingObject_Phi
MonitorElement * h1_matchingObject_Phi
Definition: ElectronAnalyzer.h:325
ElectronAnalyzer::beamSpotTag_
edm::EDGetTokenT< reco::BeamSpot > beamSpotTag_
Definition: ElectronAnalyzer.h:60
ElectronAnalyzer::~ElectronAnalyzer
virtual ~ElectronAnalyzer()
Definition: ElectronAnalyzer.cc:157
ElectronAnalyzer::h1_sigIEtaIEta5x5_endcaps
MonitorElement * h1_sigIEtaIEta5x5_endcaps
Definition: ElectronAnalyzer.h:266
ElectronAnalyzer::h1_sclSigEtaEta_endcaps
MonitorElement * h1_sclSigEtaEta_endcaps
Definition: ElectronAnalyzer.h:264
reco::BeamSpot::position
const Point & position() const
position
Definition: BeamSpot.h:62
ElectronDqmAnalyzerBase::bookH1
MonitorElement * bookH1(DQMStore::IBooker &, const std::string &name, const std::string &title, int nchX, double lowX, double highX, const std::string &titleX="", const std::string &titleY="Events", Option_t *option="E1 P")
Definition: ElectronDqmAnalyzerBase.cc:95
edm::RunNumber_t
unsigned int RunNumber_t
Definition: RunLumiEventNumber.h:14
edm::false
volatile std::atomic< bool > shutdown_flag false
Definition: UnixSignalHandlers.cc:22
ElectronAnalyzer::py_foundHitsVsEta
MonitorElement * py_foundHitsVsEta
Definition: ElectronAnalyzer.h:205
ElectronAnalyzer::vertexCollection_
edm::EDGetTokenT< reco::VertexCollection > vertexCollection_
Definition: ElectronAnalyzer.h:59
ElectronAnalyzer::idCut
bool idCut(const reco::GsfElectronCollection::const_iterator &gsfIter)
Definition: ElectronAnalyzer.cc:554
ElectronAnalyzer::hcalIso03Depth2MaxEndcaps_
double hcalIso03Depth2MaxEndcaps_
Definition: ElectronAnalyzer.h:107
reco::GsfElectron::PflowIsolationVariables::sumChargedHadronPt
float sumChargedHadronPt
sum-pt of charged Hadron // old float chargedHadronIso ;
Definition: GsfElectron.h:559
ElectronAnalyzer::hadronicOverEmMaxBarrel_
double hadronicOverEmMaxBarrel_
Definition: ElectronAnalyzer.h:99
ElectronAnalyzer::h1_matchedObject_Eta
MonitorElement * h1_matchedObject_Eta
Definition: ElectronAnalyzer.h:314
ElectronAnalyzer::h1_matchingObject_Pt
MonitorElement * h1_matchingObject_Pt
Definition: ElectronAnalyzer.h:324
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
ElectronAnalyzer::h1_dPhiEleCl_propOut_barrel
MonitorElement * h1_dPhiEleCl_propOut_barrel
Definition: ElectronAnalyzer.h:252
ElectronAnalyzer::h1_mva
MonitorElement * h1_mva
Definition: ElectronAnalyzer.h:283
edm::Run
Definition: Run.h:41
ElectronAnalyzer::isolationCut
bool isolationCut(const reco::GsfElectronCollection::const_iterator &gsfIter, double vertexTIP)
Definition: ElectronAnalyzer.cc:539
gsfElectrons_cfi.gsfElectrons
tuple gsfElectrons
Definition: gsfElectrons_cfi.py:174
Generated for CMSSW Reference Manual by doxygen 1.8.5