CMS 3D CMS Logo
TkConvValidator.cc
Go to the documentation of this file.
1 #include <iostream>
2 //
3 #include "FWCore/ServiceRegistry/interface/Service.h"
4 #include "CommonTools/UtilAlgos/interface/TFileService.h"
5 //
6 #include "Validation/RecoEgamma/plugins/TkConvValidator.h"
7 
8 //
9 #include "FWCore/Framework/interface/MakerMacros.h"
10 #include "FWCore/MessageLogger/interface/MessageLogger.h"
11 #include "FWCore/Utilities/interface/Exception.h"
12 //
13 #include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
14 #include "SimDataFormats/Track/interface/SimTrack.h"
15 #include "SimDataFormats/Track/interface/SimTrackContainer.h"
16 #include "SimDataFormats/Vertex/interface/SimVertex.h"
17 #include "SimDataFormats/Vertex/interface/SimVertexContainer.h"
18 //
19 #include "SimDataFormats/TrackingAnalysis/interface/TrackingParticleFwd.h"
20 #include "SimDataFormats/TrackingAnalysis/interface/TrackingParticle.h"
21 #include "SimDataFormats/TrackingAnalysis/interface/TrackingVertexContainer.h"
22 //
23 #include "SimTracker/TrackerHitAssociation/interface/TrackerHitAssociator.h"
24 #include "SimTracker/Records/interface/TrackAssociatorRecord.h"
25 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
26 #include "Geometry/CommonDetUnit/interface/GeomDet.h"
27 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
28 #include "Geometry/CommonDetUnit/interface/GeomDetUnit.h"
29 //
30 #include "TrackingTools/TrajectoryState/interface/FreeTrajectoryState.h"
31 #include "TrackingTools/PatternTools/interface/TwoTrackMinimumDistance.h"
32 #include "TrackingTools/TransientTrack/interface/TrackTransientTrack.h"
33 //
34 #include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
35 #include "TrackingTools/Records/interface/TransientTrackRecord.h"
36 #include "CLHEP/Units/GlobalPhysicalConstants.h"
37 #include "CommonTools/Statistics/interface/ChiSquaredProbability.h"
38 //
39 #include "MagneticField/Engine/interface/MagneticField.h"
40 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
41 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
42 #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
43 #include "TrackingTools/Records/interface/TrackingComponentsRecord.h"
44 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
45 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"
46 #include "Geometry/CommonDetUnit/interface/GeomDet.h"
47 #include "TrackingTools/MaterialEffects/interface/PropagatorWithMaterial.h"
48 
49 //
50 #include "DataFormats/Common/interface/Handle.h"
51 #include "DataFormats/TrackReco/interface/Track.h"
52 #include "DataFormats/TrackReco/interface/TrackExtra.h"
53 #include "DataFormats/TrackReco/interface/TrackFwd.h"
54 #include "DataFormats/EgammaCandidates/interface/Conversion.h"
55 #include "DataFormats/EgammaCandidates/interface/Photon.h"
56 #include "DataFormats/EgammaReco/interface/SuperCluster.h"
57 #include "DataFormats/TrackCandidate/interface/TrackCandidateCollection.h"
58 #include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
59 #include "DataFormats/EgammaReco/interface/BasicCluster.h"
60 #include "DataFormats/EgammaReco/interface/BasicCluster.h"
61 #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
62 #include "DataFormats/CaloRecHit/interface/CaloCluster.h"
63 #include "DataFormats/CaloRecHit/interface/CaloClusterFwd.h"
64 #include "DataFormats/DetId/interface/DetId.h"
65 #include "DataFormats/JetReco/interface/GenJetCollection.h"
66 #include "DataFormats/HepMCCandidate/interface/GenParticle.h"
67 #include "DataFormats/GeometrySurface/interface/BoundCylinder.h"
68 #include "DataFormats/GeometrySurface/interface/BoundDisk.h"
69 #include "DataFormats/GeometrySurface/interface/SimpleCylinderBounds.h"
70 #include "DataFormats/GeometrySurface/interface/SimpleDiskBounds.h"
71 #include "DataFormats/Math/interface/deltaPhi.h"
72 
73 //
74 #include "RecoEgamma/EgammaMCTools/interface/PhotonMCTruthFinder.h"
75 #include "RecoEgamma/EgammaMCTools/interface/PhotonMCTruth.h"
76 #include "RecoEgamma/EgammaMCTools/interface/ElectronMCTruth.h"
77 #include "RecoEcal/EgammaCoreTools/interface/EcalClusterTools.h"
78 #include "Geometry/Records/interface/CaloGeometryRecord.h"
79 #include "Geometry/CaloTopology/interface/CaloTopology.h"
80 #include "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"
81 
82 #include "RecoEgamma/EgammaPhotonAlgos/interface/ConversionHitChecker.h"
83 
84 //
85 //
86 #include "TFile.h"
87 #include "TH1.h"
88 #include "TH2.h"
89 #include "TTree.h"
90 #include "TVector3.h"
91 #include "TProfile.h"
92 //
101 using namespace std;
102 
103 
104 TkConvValidator::TkConvValidator( const edm::ParameterSet& pset )
105  {
106 
107  fName_ = pset.getUntrackedParameter<std::string>("Name");
108  verbosity_ = pset.getUntrackedParameter<int>("Verbosity");
109  parameters_ = pset;
110 
111  photonCollectionProducer_ = pset.getParameter<std::string>("phoProducer");
112  photonCollection_ = pset.getParameter<std::string>("photonCollection");
113  photonCollectionPr_Token_ = consumes<reco::PhotonCollection> (
114  edm::InputTag(photonCollectionProducer_,
115  photonCollection_));
116 
117  conversionCollectionProducer_ = pset.getParameter<std::string>("convProducer");
118  conversionCollection_ = pset.getParameter<std::string>("conversionCollection");
119  conversionCollectionPr_Token_ = consumes<reco::ConversionCollection> (
120  edm::InputTag(conversionCollectionProducer_,
121  conversionCollection_));
122 
123  // conversionTrackProducer_ = pset.getParameter<std::string>("trackProducer");
124  dqmpath_ = pset.getParameter<std::string>("dqmpath");
125  minPhoEtCut_ = pset.getParameter<double>("minPhoEtCut");
126  generalTracksOnly_ = pset.getParameter<bool>("generalTracksOnly");
127  arbitratedMerged_ = pset.getParameter<bool>("arbitratedMerged");
128  arbitratedEcalSeeded_ = pset.getParameter<bool>("arbitratedEcalSeeded");
129  ecalalgotracks_ = pset.getParameter<bool>("ecalalgotracks");
130  highPurity_ = pset.getParameter<bool>("highPurity");
131  minProb_ = pset.getParameter<double>("minProb");
132  maxHitsBeforeVtx_ = pset.getParameter<uint>("maxHitsBeforeVtx");
133  minLxy_ = pset.getParameter<double>("minLxy");
134  isRunCentrally_= pset.getParameter<bool>("isRunCentrally");
135 
136  offline_pvToken_ = consumes<reco::VertexCollection>(
137  pset.getUntrackedParameter<edm::InputTag> ("offlinePV",
138  edm::InputTag("offlinePrimaryVertices")));
139  beamspotToken_ = consumes<reco::BeamSpot>(
140  pset.getUntrackedParameter<edm::InputTag> ("beamspot",
141  edm::InputTag("offlineBeamSpot")));
142  g4_simTk_Token_ = consumes<edm::SimTrackContainer>(pset.getParameter<edm::InputTag>("simTracks"));
143  g4_simVtx_Token_ = consumes<edm::SimVertexContainer>(pset.getParameter<edm::InputTag>("simTracks"));
144 
145  tpSelForEff_Token_ = consumes<TrackingParticleRefVector> (
146  edm::InputTag("tpSelecForEfficiency"));
147  tpSelForFake_Token_ = consumes<TrackingParticleRefVector> (
148  edm::InputTag("tpSelecForFakeRate"));
149  //hepMC_Token_ = consumes<edm::HepMCProduct>(edm::InputTag("generatorSmeared"));
150  //genjets_Token_ = consumes<reco::GenJetCollection>(
151  // edm::InputTag("ak4GenJets"));
152 
153  trackAssociator_Token_ = consumes<reco::TrackToTrackingParticleAssociator>(edm::InputTag("trackAssociatorByHitsForConversionValidation"));
154  }
155 
156 
157 
158 TkConvValidator::~TkConvValidator() {}
159 
160 
161 
162 
163 void TkConvValidator::bookHistograms( DQMStore::IBooker & iBooker, edm::Run const & run, edm::EventSetup const & es) {
164 
165  nEvt_=0;
166  nEntry_=0;
167  nRecConv_=0;
168  nRecConvAss_=0;
169  nRecConvAssWithEcal_=0;
170 
171  nInvalidPCA_=0;
172 
173  dbe_ = 0;
174  dbe_ = edm::Service<DQMStore>().operator->();
175 
176 
177  double etMin = parameters_.getParameter<double>("etMin");
178  double etMax = parameters_.getParameter<double>("etMax");
179  int etBin = parameters_.getParameter<int>("etBin");
180 
181 
182  double resMin = parameters_.getParameter<double>("resMin");
183  double resMax = parameters_.getParameter<double>("resMax");
184  int resBin = parameters_.getParameter<int>("resBin");
185 
186  double etaMin = parameters_.getParameter<double>("etaMin");
187  double etaMax = parameters_.getParameter<double>("etaMax");
188  int etaBin = parameters_.getParameter<int>("etaBin");
189  int etaBin2 = parameters_.getParameter<int>("etaBin2");
190 
191 
192  double phiMin = parameters_.getParameter<double>("phiMin");
193  double phiMax = parameters_.getParameter<double>("phiMax");
194  int phiBin = parameters_.getParameter<int>("phiBin");
195 
196 
197  double rMin = parameters_.getParameter<double>("rMin");
198  double rMax = parameters_.getParameter<double>("rMax");
199  int rBin = parameters_.getParameter<int>("rBin");
200 
201  double zMin = parameters_.getParameter<double>("zMin");
202  double zMax = parameters_.getParameter<double>("zMax");
203  int zBin = parameters_.getParameter<int>("zBin");
204 
205  double dPhiTracksMin = parameters_.getParameter<double>("dPhiTracksMin");
206  double dPhiTracksMax = parameters_.getParameter<double>("dPhiTracksMax");
207  int dPhiTracksBin = parameters_.getParameter<int>("dPhiTracksBin");
208 
209  double eoverpMin = parameters_.getParameter<double>("eoverpMin");
210  double eoverpMax = parameters_.getParameter<double>("eoverpMax");
211  int eoverpBin = parameters_.getParameter<int>("eoverpBin");
212 
213 
214  // double dEtaTracksMin = parameters_.getParameter<double>("dEtaTracksMin"); // unused
215  // double dEtaTracksMax = parameters_.getParameter<double>("dEtaTracksMax"); // unused
216  // int dEtaTracksBin = parameters_.getParameter<int>("dEtaTracksBin"); // unused
217 
218  double dCotTracksMin = parameters_.getParameter<double>("dCotTracksMin");
219  double dCotTracksMax = parameters_.getParameter<double>("dCotTracksMax");
220  int dCotTracksBin = parameters_.getParameter<int>("dCotTracksBin");
221 
222 
223  double chi2Min = parameters_.getParameter<double>("chi2Min");
224  double chi2Max = parameters_.getParameter<double>("chi2Max");
225 
226 
227  double rMinForXray = parameters_.getParameter<double>("rMinForXray");
228  double rMaxForXray = parameters_.getParameter<double>("rMaxForXray");
229  int rBinForXray = parameters_.getParameter<int>("rBinForXray");
230  double zMinForXray = parameters_.getParameter<double>("zMinForXray");
231  double zMaxForXray = parameters_.getParameter<double>("zMaxForXray");
232  int zBinForXray = parameters_.getParameter<int>("zBinForXray");
233  int zBin2ForXray = parameters_.getParameter<int>("zBin2ForXray");
234 
235  minPhoPtForEffic = parameters_.getParameter<double>("minPhoPtForEffic");
236  maxPhoEtaForEffic = parameters_.getParameter<double>("maxPhoEtaForEffic");
237  maxPhoZForEffic = parameters_.getParameter<double>("maxPhoZForEffic");
238  maxPhoRForEffic = parameters_.getParameter<double>("maxPhoRForEffic");
239  minPhoPtForPurity = parameters_.getParameter<double>("minPhoPtForPurity");
240  maxPhoEtaForPurity = parameters_.getParameter<double>("maxPhoEtaForPurity");
241  maxPhoZForPurity = parameters_.getParameter<double>("maxPhoZForPurity");
242  maxPhoRForPurity = parameters_.getParameter<double>("maxPhoRForPurity");
243 
244 
245  if (dbe_) {
246 
248  // SC from reco photons
249 
250  //TString simfolder = TString(
251  std::string simpath = dqmpath_ + "SimulationInfo";
252  iBooker.setCurrentFolder(simpath);
253  //
254  // simulation information about conversions
256  std::string histname = "nOfSimConversions";
257  h_nSimConv_[0] = iBooker.book1D(histname,"# of Sim conversions per event ",20,-0.5,19.5);
259  histname = "h_AllSimConvEta";
260  h_AllSimConv_[0] = iBooker.book1D(histname," All conversions: simulated #eta",etaBin2,etaMin,etaMax);
261  histname = "h_AllSimConvPhi";
262  h_AllSimConv_[1] = iBooker.book1D(histname," All conversions: simulated #phi",phiBin,phiMin,phiMax);
263  histname = "h_AllSimConvR";
264  h_AllSimConv_[2] = iBooker.book1D(histname," All conversions: simulated R",rBin,rMin,rMax);
265  histname = "h_AllSimConvZ";
266  h_AllSimConv_[3] = iBooker.book1D(histname," All conversions: simulated Z",zBin,zMin,zMax);
267  histname = "h_AllSimConvEt";
268  h_AllSimConv_[4] = iBooker.book1D(histname," All conversions: simulated Et",etBin,etMin,etMax);
269  //
270  histname = "nOfVisSimConversions";
271  h_nSimConv_[1] = iBooker.book1D(histname,"# of Sim conversions per event ",20,-0.5,19.5);
272  histname = "h_VisSimConvEta";
273  h_VisSimConv_[0] = iBooker.book1D(histname," All vis conversions: simulated #eta",etaBin2,etaMin, etaMax);
274  histname = "h_VisSimConvPhi";
275  h_VisSimConv_[1] = iBooker.book1D(histname," All vis conversions: simulated #phi",phiBin,phiMin, phiMax);
276  histname = "h_VisSimConvR";
277  h_VisSimConv_[2] = iBooker.book1D(histname," All vis conversions: simulated R",rBin,rMin,rMax);
278  histname = "h_VisSimConvZ";
279  h_VisSimConv_[3] = iBooker.book1D(histname," All vis conversions: simulated Z",zBin,zMin, zMax);
280  histname = "h_VisSimConvEt";
281  h_VisSimConv_[4] = iBooker.book1D(histname," All vis conversions: simulated Et",etBin,etMin, etMax);
282 
283  //
284  histname = "h_SimConvTwoMTracksEta";
285  h_SimConvTwoMTracks_[0] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks: simulated #eta",etaBin2,etaMin, etaMax);
286  histname = "h_SimConvTwoMTracksPhi";
287  h_SimConvTwoMTracks_[1] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks: simulated #phi",phiBin,phiMin, phiMax);
288  histname = "h_SimConvTwoMTracksR";
289  h_SimConvTwoMTracks_[2] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks: simulated R",rBin,rMin, rMax);
290  histname = "h_SimConvTwoMTracksZ";
291  h_SimConvTwoMTracks_[3] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks: simulated Z",zBin,zMin, zMax);
292  histname = "h_SimConvTwoMTracksEt";
293  h_SimConvTwoMTracks_[4] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks: simulated Et",etBin,etMin, etMax);
294  //
295  histname = "h_SimConvTwoTracksEta";
296  h_SimConvTwoTracks_[0] = iBooker.book1D(histname," All vis conversions with 2 reco tracks: simulated #eta",etaBin2,etaMin, etaMax);
297  histname = "h_SimConvTwoTracksPhi";
298  h_SimConvTwoTracks_[1] = iBooker.book1D(histname," All vis conversions with 2 reco tracks: simulated #phi",phiBin,phiMin, phiMax);
299  histname = "h_SimConvTwoTracksR";
300  h_SimConvTwoTracks_[2] = iBooker.book1D(histname," All vis conversions with 2 reco tracks: simulated R",rBin,rMin, rMax);
301  histname = "h_SimConvTwoTracksZ";
302  h_SimConvTwoTracks_[3] = iBooker.book1D(histname," All vis conversions with 2 reco tracks: simulated Z",zBin,zMin, zMax);
303  histname = "h_SimConvTwoTracksEt";
304  h_SimConvTwoTracks_[4] = iBooker.book1D(histname," All vis conversions with 2 reco tracks: simulated Et",etBin,etMin, etMax);
305  //
306  histname = "h_SimConvTwoMTracksEtaAndVtxPGT0";
307  h_SimConvTwoMTracksAndVtxPGT0_[0] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks + vertex: simulated #eta",etaBin2,etaMin, etaMax);
308  histname = "h_SimConvTwoMTracksPhiAndVtxPGT0";
309  h_SimConvTwoMTracksAndVtxPGT0_[1] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks + vertex: simulated #phi",phiBin,phiMin, phiMax);
310  histname = "h_SimConvTwoMTracksRAndVtxPGT0";
311  h_SimConvTwoMTracksAndVtxPGT0_[2] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks + vertex: simulated R",rBin,rMin, rMax);
312  histname = "h_SimConvTwoMTracksZAndVtxPGT0";
313  h_SimConvTwoMTracksAndVtxPGT0_[3] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks + vertex: simulated Z",zBin,zMin, zMax);
314  histname = "h_SimConvTwoMTracksEtAndVtxPGT0";
315  h_SimConvTwoMTracksAndVtxPGT0_[4] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks + vertex: simulated Et",etBin,etMin, etMax);
316 
317  //
318  histname = "h_SimConvTwoMTracksEtaAndVtxPGT0005";
319  h_SimConvTwoMTracksAndVtxPGT0005_[0] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks + vertex: simulated #eta",etaBin2,etaMin, etaMax);
320  histname = "h_SimConvTwoMTracksPhiAndVtxPGT0005";
321  h_SimConvTwoMTracksAndVtxPGT0005_[1] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks + vertex: simulated #phi",phiBin,phiMin, phiMax);
322  histname = "h_SimConvTwoMTracksRAndVtxPGT0005";
323  h_SimConvTwoMTracksAndVtxPGT0005_[2] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks + vertex: simulated R",rBin,rMin, rMax);
324  histname = "h_SimConvTwoMTracksZAndVtxPGT0005";
325  h_SimConvTwoMTracksAndVtxPGT0005_[3] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks + vertex: simulated Z",zBin,zMin, zMax);
326  histname = "h_SimConvTwoMTracksEtAndVtxPGT0005";
327  h_SimConvTwoMTracksAndVtxPGT0005_[4] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks + vertex: simulated Et",etBin,etMin, etMax);
328 
329  histname = "h_SimRecConvTwoMTracksEta";
330  h_SimRecConvTwoMTracks_[0] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks: simulated #eta",etaBin2,etaMin, etaMax);
331  histname = "h_SimRecConvTwoMTracksPhi";
332  h_SimRecConvTwoMTracks_[1] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks: simulated #phi",phiBin,phiMin, phiMax);
333  histname = "h_SimRecConvTwoMTracksR";
334  h_SimRecConvTwoMTracks_[2] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks: simulated R",rBin,rMin, rMax);
335  histname = "h_SimRecConvTwoMTracksZ";
336  h_SimRecConvTwoMTracks_[3] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks: simulated Z",zBin,zMin, zMax);
337  histname = "h_SimRecConvTwoMTracksEt";
338  h_SimRecConvTwoMTracks_[4] = iBooker.book1D(histname," All vis conversions with 2 reco-matching tracks: simulated Et",etBin,etMin, etMax);
339  //
340 
341 
342  h_SimConvEtaPix_[0] = iBooker.book1D("simConvEtaPix"," sim converted Photon Eta: Pix ",etaBin,etaMin, etaMax) ;
343  h_simTkPt_ = iBooker.book1D("simTkPt","Sim conversion tracks pt ",etBin*3,0.,etMax);
344  h_simTkEta_ = iBooker.book1D("simTkEta","Sim conversion tracks eta ",etaBin,etaMin,etaMax);
345 
346  h_simConvVtxRvsZ_[0] = iBooker.book2D("simConvVtxRvsZAll"," Photon Sim conversion vtx position",zBinForXray, zMinForXray, zMaxForXray, rBinForXray, rMinForXray, rMaxForXray);
347  h_simConvVtxRvsZ_[1] = iBooker.book2D("simConvVtxRvsZBarrel"," Photon Sim conversion vtx position",zBinForXray, zMinForXray, zMaxForXray, rBinForXray, rMinForXray, rMaxForXray);
348  h_simConvVtxRvsZ_[2] = iBooker.book2D("simConvVtxRvsZEndcap"," Photon Sim conversion vtx position",zBin2ForXray, zMinForXray, zMaxForXray, rBinForXray, rMinForXray, rMaxForXray);
349  h_simConvVtxRvsZ_[3] = iBooker.book2D("simConvVtxRvsZBarrel2"," Photon Sim conversion vtx position when reco R<4cm",zBinForXray, zMinForXray, zMaxForXray, rBinForXray, rMinForXray, rMaxForXray);
350  h_simConvVtxYvsX_ = iBooker.book2D("simConvVtxYvsXTrkBarrel"," Photon Sim conversion vtx position, (x,y) eta<1 ",100, -80., 80., 100, -80., 80.);
351 
352  std::string convpath = dqmpath_ + "ConversionInfo";
353  iBooker.setCurrentFolder(convpath);
354 
355  histname="nConv";
356  h_nConv_[0][0] = iBooker.book1D(histname+"All","Number Of Conversions per isolated candidates per events: All Ecal ",10,-0.5, 9.5);
357  h_nConv_[0][1] = iBooker.book1D(histname+"Barrel","Number Of Conversions per isolated candidates per events: Ecal Barrel ",10,-0.5, 9.5);
358  h_nConv_[0][2] = iBooker.book1D(histname+"Endcap","Number Of Conversions per isolated candidates per events: Ecal Endcap ",10,-0.5, 9.5);
359  h_nConv_[1][0] = iBooker.book1D(histname+"All_Ass","Number Of associated Conversions per isolated candidates per events: All Ecal ",10,-0.5, 9.5);
360 
361  h_convEta_[0][0] = iBooker.book1D("convEta"," converted Photon Eta ",etaBin,etaMin, etaMax) ;
362  h_convEtaMatchSC_[0][0] = iBooker.book1D("convEtaMatchSC"," converted Photon Eta when SC is matched ",etaBin,etaMin, etaMax) ;
363  h_convEta2_[0][0] = iBooker.book1D("convEta2"," converted Photon Eta ",etaBin2,etaMin, etaMax) ;
364  h_convPhi_[0][0] = iBooker.book1D("convPhi"," converted Photon Phi ",phiBin,phiMin,phiMax) ;
365  h_convR_[0][0] = iBooker.book1D("convR"," converted photon R",rBin,rMin, rMax);
366  h_convZ_[0][0] = iBooker.book1D("convZ"," converted photon Z",zBin,zMin, zMax);
367  h_convPt_[0][0] = iBooker.book1D("convPt"," conversions Transverse Energy: all eta ", etBin,etMin, etMax);
368 
369  h_convEta_[1][0] = iBooker.book1D("convEtaAss2"," Matched converted Photon Eta ",etaBin2,etaMin, etaMax) ;
370  h_convEta_[1][1] = iBooker.book1D("convEtaAss"," Matched converted Photon Eta ",etaBin,etaMin, etaMax) ;
371  h_convEtaMatchSC_[1][0] = iBooker.book1D("convEtaMatchSCAss"," converted Photon Eta when SC is matched ",etaBin,etaMin, etaMax) ;
372  h_convPhi_[1][0] = iBooker.book1D("convPhiAss"," Matched converted Photon Phi ",phiBin,phiMin,phiMax) ;
373  h_convR_[1][0] = iBooker.book1D("convRAss"," Matched converted photon R",rBin,rMin, rMax);
374  h_convZ_[1][0] = iBooker.book1D("convZAss"," Matched converted photon Z",zBin,zMin, zMax);
375  h_convPt_[1][0] = iBooker.book1D("convPtAss","Matched conversions Transverse Energy: all eta ", etBin,etMin, etMax);
376 
377  h_convEta_[2][0] = iBooker.book1D("convEtaFake2"," Fake converted Photon Eta ",etaBin2,etaMin, etaMax) ;
378  h_convEta_[2][1] = iBooker.book1D("convEtaFake"," Fake converted Photon Eta ",etaBin,etaMin, etaMax) ;
379  h_convEtaMatchSC_[2][0] = iBooker.book1D("convEtaMatchSCFake"," converted Photon Eta when SC is matched ",etaBin,etaMin, etaMax) ;
380  h_convPhi_[2][0] = iBooker.book1D("convPhiFake"," Fake converted Photon Phi ",phiBin,phiMin,phiMax) ;
381  h_convR_[2][0] = iBooker.book1D("convRFake"," Fake converted photon R",rBin,rMin, rMax);
382  h_convZ_[2][0] = iBooker.book1D("convZFake"," Fake converted photon Z",zBin,zMin, zMax);
383  h_convPt_[2][0] = iBooker.book1D("convPtFake","Fake conversions Transverse Energy: all eta ", etBin,etMin, etMax);
384 
385  h_convRplot_ = iBooker.book1D("convRplot"," converted photon R",600, 0.,120.);
386  h_convZplot_ = iBooker.book1D("convZplot"," converted photon Z",320,-160.,160.);
387 
388  histname = "convSCdPhi";
389  h_convSCdPhi_[0][0] = iBooker.book1D(histname+"All","dPhi between SC and conversion",100, -0.1,0.1);
390  h_convSCdPhi_[0][1] = iBooker.book1D(histname+"Barrel"," dPhi between SC and conversion: Barrel",100, -0.1,0.1);
391  h_convSCdPhi_[0][2] = iBooker.book1D(histname+"Endcap"," dPhi between SC and conversion: Endcap",100, -0.1,0.1);
392  h_convSCdPhi_[1][0] = iBooker.book1D(histname+"All_Ass","dPhi between SC and conversion",100, -0.1,0.1);
393  h_convSCdPhi_[1][1] = iBooker.book1D(histname+"Barrel_Ass"," dPhi between SC and conversion: Barrel",100, -0.1,0.1);
394  h_convSCdPhi_[1][2] = iBooker.book1D(histname+"Endcap_Ass"," dPhi between SC and conversion: Endcap",100, -0.1,0.1);
395  h_convSCdPhi_[2][0] = iBooker.book1D(histname+"All_Fakes","dPhi between SC and conversion",100, -0.1,0.1);
396  h_convSCdPhi_[2][1] = iBooker.book1D(histname+"Barrel_Fakes"," dPhi between SC and conversion: Barrel",100, -0.1,0.1);
397  h_convSCdPhi_[2][2] = iBooker.book1D(histname+"Endcap_Fakes"," dPhi between SC and conversion: Endcap",100, -0.1,0.1);
398  histname = "convSCdEta";
399  h_convSCdEta_[0][0] = iBooker.book1D(histname+"All"," dEta between SC and conversion",100, -0.1,0.1);
400  h_convSCdEta_[0][1] = iBooker.book1D(histname+"Barrel"," dEta between SC and conversion: Barrel",100, -0.1,0.1);
401  h_convSCdEta_[0][2] = iBooker.book1D(histname+"Endcap"," dEta between SC and conversion: Endcap",100, -0.1,0.1);
402  h_convSCdEta_[1][0] = iBooker.book1D(histname+"All_Ass"," dEta between SC and conversion",100, -0.1,0.1);
403  h_convSCdEta_[1][1] = iBooker.book1D(histname+"Barrel_Ass"," dEta between SC and conversion: Barrel",100, -0.1,0.1);
404  h_convSCdEta_[1][2] = iBooker.book1D(histname+"Endcap_Ass"," dEta between SC and conversion: Endcap",100, -0.1,0.1);
405  h_convSCdEta_[2][0] = iBooker.book1D(histname+"All_Fakes"," dEta between SC and conversion",100, -0.1,0.1);
406  h_convSCdEta_[2][1] = iBooker.book1D(histname+"Barrel_Fakes"," dEta between SC and conversion: Barrel",100, -0.1,0.1);
407  h_convSCdEta_[2][2] = iBooker.book1D(histname+"Endcap_Fakes"," dEta between SC and conversion: Endcap",100, -0.1,0.1);
408 
409  histname = "convPtRes";
410  h_convPtRes_[0] = iBooker.book1D(histname+"All"," Conversion Pt rec/true : All ecal ", resBin,resMin, resMax);
411  h_convPtRes_[1] = iBooker.book1D(histname+"Barrel"," Conversion Pt rec/true : Barrel ",resBin,resMin, resMax);
412  h_convPtRes_[2] = iBooker.book1D(histname+"Endcap"," Conversion Pt rec/true : Endcap ",resBin,resMin, resMax);
413 
414 
415  histname="hInvMass";
416  h_invMass_[0][0]= iBooker.book1D(histname+"All_AllTracks"," Photons:Tracks from conversion: Pair invariant mass: all Ecal ",100, 0., 1.5);
417  h_invMass_[0][1]= iBooker.book1D(histname+"Barrel_AllTracks"," Photons:Tracks from conversion: Pair invariant mass: Barrel Ecal ",100, 0., 1.5);
418  h_invMass_[0][2]= iBooker.book1D(histname+"Endcap_AllTracks"," Photons:Tracks from conversion: Pair invariant mass: Endcap Ecal ",100, 0., 1.5);
419  //
420  h_invMass_[1][0]= iBooker.book1D(histname+"All_AssTracks"," Photons:Tracks from conversion: Pair invariant mass: all Ecal ",100, 0., 1.5);
421  h_invMass_[1][1]= iBooker.book1D(histname+"Barrel_AssTracks"," Photons:Tracks from conversion: Pair invariant mass: Barrel Ecal ",100, 0., 1.5);
422  h_invMass_[1][2]= iBooker.book1D(histname+"Endcap_AssTracks"," Photons:Tracks from conversion: Pair invariant mass: Endcap Ecal ",100, 0., 1.5);
423  //
424  h_invMass_[2][0]= iBooker.book1D(histname+"All_FakeTracks"," Photons:Tracks from conversion: Pair invariant mass: all Ecal ",100, 0., 1.5);
425  h_invMass_[2][1]= iBooker.book1D(histname+"Barrel_FakeTracks"," Photons:Tracks from conversion: Pair invariant mass: Barrel Ecal ",100, 0., 1.5);
426  h_invMass_[2][2]= iBooker.book1D(histname+"Endcap_FaleTracks"," Photons:Tracks from conversion: Pair invariant mass: Endcap Ecal ",100, 0., 1.5);
427 
428 
429 
430  histname="hDPhiTracksAtVtx";
431  h_DPhiTracksAtVtx_[0][0] =iBooker.book1D(histname+"All", " Photons:Tracks from conversions: #delta#phi Tracks at vertex: all Ecal",dPhiTracksBin,dPhiTracksMin,dPhiTracksMax);
432  h_DPhiTracksAtVtx_[0][1] =iBooker.book1D(histname+"Barrel", " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Barrel Ecal",dPhiTracksBin,dPhiTracksMin,dPhiTracksMax);
433  h_DPhiTracksAtVtx_[0][2] =iBooker.book1D(histname+"Endcap", " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Endcap Ecal",dPhiTracksBin,dPhiTracksMin,dPhiTracksMax);
434  h_DPhiTracksAtVtx_[1][0] =iBooker.book1D(histname+"All_Ass", " Photons:Tracks from conversions: #delta#phi Tracks at vertex: all Ecal",dPhiTracksBin,dPhiTracksMin,dPhiTracksMax);
435  h_DPhiTracksAtVtx_[1][1] =iBooker.book1D(histname+"Barrel_Ass", " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Barrel Ecal",dPhiTracksBin,dPhiTracksMin,dPhiTracksMax);
436  h_DPhiTracksAtVtx_[1][2] =iBooker.book1D(histname+"Endcap_Ass", " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Endcap Ecal",dPhiTracksBin,dPhiTracksMin,dPhiTracksMax);
437  h_DPhiTracksAtVtx_[2][0] =iBooker.book1D(histname+"All_Fakes", " Photons:Tracks from conversions: #delta#phi Tracks at vertex: all Ecal",dPhiTracksBin,dPhiTracksMin,dPhiTracksMax);
438  h_DPhiTracksAtVtx_[2][1] =iBooker.book1D(histname+"Barrel_Fakes", " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Barrel Ecal",dPhiTracksBin,dPhiTracksMin,dPhiTracksMax);
439  h_DPhiTracksAtVtx_[2][2] =iBooker.book1D(histname+"Endcap_Fakes", " Photons:Tracks from conversions: #delta#phi Tracks at vertex: Endcap Ecal",dPhiTracksBin,dPhiTracksMin,dPhiTracksMax);
440 
441 
442 
443  histname="hDPhiTracksAtVtxVsEta";
444  h2_DPhiTracksAtVtxVsEta_ = iBooker.book2D(histname+"All"," Photons:Tracks from conversions: #delta#phi Tracks at vertex vs #eta",etaBin2,etaMin, etaMax,100, -0.5, 0.5);
445  histname="pDPhiTracksAtVtxVsEta";
446  p_DPhiTracksAtVtxVsEta_ = iBooker.bookProfile(histname+"All"," Photons:Tracks from conversions: #delta#phi Tracks at vertex vs #eta ",etaBin2,etaMin, etaMax, 100, -0.5, 0.5,"");
447 
448  histname="hDPhiTracksAtVtxVsR";
449  h2_DPhiTracksAtVtxVsR_ = iBooker.book2D(histname+"All"," Photons:Tracks from conversions: #delta#phi Tracks at vertex vs R",rBin,rMin, rMax,100, -0.5, 0.5);
450  histname="pDPhiTracksAtVtxVsR";
451  p_DPhiTracksAtVtxVsR_ = iBooker.bookProfile(histname+"All"," Photons:Tracks from conversions: #delta#phi Tracks at vertex vs R ",rBin,rMin, rMax,100, -0.5, 0.5,"");
452 
453 
454  histname="hDCotTracks";
455  h_DCotTracks_[0][0]= iBooker.book1D(histname+"All"," Photons:Tracks from conversions #delta cotg(#Theta) Tracks: all Ecal ",dCotTracksBin,dCotTracksMin,dCotTracksMax);
456  h_DCotTracks_[0][1]= iBooker.book1D(histname+"Barrel"," Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Barrel Ecal ",dCotTracksBin,dCotTracksMin,dCotTracksMax);
457  h_DCotTracks_[0][2]= iBooker.book1D(histname+"Endcap"," Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Endcap Ecal ",dCotTracksBin,dCotTracksMin,dCotTracksMax);
458  h_DCotTracks_[1][0]= iBooker.book1D(histname+"All_Ass"," Photons:Tracks from conversions #delta cotg(#Theta) Tracks: all Ecal ",dCotTracksBin,dCotTracksMin,dCotTracksMax);
459  h_DCotTracks_[1][1]= iBooker.book1D(histname+"Barrel_Ass"," Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Barrel Ecal ",dCotTracksBin,dCotTracksMin,dCotTracksMax);
460  h_DCotTracks_[1][2]= iBooker.book1D(histname+"Endcap_Ass"," Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Endcap Ecal ",dCotTracksBin,dCotTracksMin,dCotTracksMax);
461  h_DCotTracks_[2][0]= iBooker.book1D(histname+"All_Fakes"," Photons:Tracks from conversions #delta cotg(#Theta) Tracks: all Ecal ",dCotTracksBin,dCotTracksMin,dCotTracksMax);
462  h_DCotTracks_[2][1]= iBooker.book1D(histname+"Barrel_Fakes"," Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Barrel Ecal ",dCotTracksBin,dCotTracksMin,dCotTracksMax);
463  h_DCotTracks_[2][2]= iBooker.book1D(histname+"Endcap_Fakes"," Photons:Tracks from conversions #delta cotg(#Theta) Tracks: Endcap Ecal ",dCotTracksBin,dCotTracksMin,dCotTracksMax);
464 
465 
466  histname="hDCotTracksVsEta";
467  h2_DCotTracksVsEta_ = iBooker.book2D(histname+"All"," Photons:Tracks from conversions: #delta cotg(#Theta) Tracks vs #eta",etaBin2,etaMin, etaMax,100, -0.2, 0.2);
468  histname="pDCotTracksVsEta";
469  p_DCotTracksVsEta_ = iBooker.bookProfile(histname+"All"," Photons:Tracks from conversions: #delta cotg(#Theta) Tracks vs #eta ",etaBin2,etaMin, etaMax, 100, -0.2, 0.2,"");
470 
471  histname="hDCotTracksVsR";
472  h2_DCotTracksVsR_ = iBooker.book2D(histname+"All"," Photons:Tracks from conversions: #delta cotg(#Theta) Tracks at vertex vs R",rBin,rMin, rMax,100, -0.2, 0.2);
473  histname="pDCotTracksVsR";
474  p_DCotTracksVsR_ = iBooker.bookProfile(histname+"All"," Photons:Tracks from conversions: #delta cotg(#Theta) Tracks at vertex vs R ",rBin,rMin, rMax,100, -0.2, 0.2,"");
475 
476 
477  histname="hDistMinAppTracks";
478  h_distMinAppTracks_[0][0]= iBooker.book1D(histname+"All"," Photons:Tracks from conversions Min Approach Dist Tracks: all Ecal ",120, -0.5, 1.0);
479  h_distMinAppTracks_[0][1]= iBooker.book1D(histname+"Barrel"," Photons:Tracks from conversions Min Approach Dist Tracks: Barrel Ecal ",120, -0.5, 1.0);
480  h_distMinAppTracks_[0][2]= iBooker.book1D(histname+"Endcap"," Photons:Tracks from conversions Min Approach Dist Tracks: Endcap Ecal ",120, -0.5, 1.0);
481  h_distMinAppTracks_[1][0]= iBooker.book1D(histname+"All_Ass"," Photons:Tracks from conversions Min Approach Dist Tracks: all Ecal ",120, -0.5, 1.0);
482  h_distMinAppTracks_[1][1]= iBooker.book1D(histname+"Barrel_Ass"," Photons:Tracks from conversions Min Approach Dist Tracks: Barrel Ecal ",120, -0.5, 1.0);
483  h_distMinAppTracks_[1][2]= iBooker.book1D(histname+"Endcap_Ass"," Photons:Tracks from conversions Min Approach Dist Tracks: Endcap Ecal ",120, -0.5, 1.0);
484  h_distMinAppTracks_[2][0]= iBooker.book1D(histname+"All_Fakes"," Photons:Tracks from conversions Min Approach Dist Tracks: all Ecal ",120, -0.5, 1.0);
485  h_distMinAppTracks_[2][1]= iBooker.book1D(histname+"Barrel_Fakes"," Photons:Tracks from conversions Min Approach Dist Tracks: Barrel Ecal ",120, -0.5, 1.0);
486  h_distMinAppTracks_[2][2]= iBooker.book1D(histname+"Endcap_Fakes"," Photons:Tracks from conversions Min Approach Dist Tracks: Endcap Ecal ",120, -0.5, 1.0);
487 
488 
489  h_convVtxRvsZ_[0] = iBooker.book2D("convVtxRvsZAll"," Photon Reco conversion vtx position",zBinForXray, zMinForXray, zMaxForXray, rBinForXray, rMinForXray, rMaxForXray);
490  h_convVtxRvsZ_[1] = iBooker.book2D("convVtxRvsZBarrel"," Photon Reco conversion vtx position",zBinForXray, zMinForXray, zMaxForXray, rBinForXray, rMinForXray, rMaxForXray);
491  h_convVtxRvsZ_[2] = iBooker.book2D("convVtxRvsZEndcap"," Photon Reco conversion vtx position",zBin2ForXray, zMinForXray, zMaxForXray, rBinForXray, rMinForXray, rMaxForXray);
492  h_convVtxYvsX_ = iBooker.book2D("convVtxYvsXTrkBarrel"," Photon Reco conversion vtx position, (x,y) eta<1 ", 1000, -60., 60., 1000, -60., 60.);
494  h_convVtxRvsZ_zoom_[0] = iBooker.book2D("convVtxRvsZBarrelZoom1"," Photon Reco conversion vtx position",zBinForXray, zMinForXray, zMaxForXray, rBinForXray, -10., 40.);
495  h_convVtxRvsZ_zoom_[1] = iBooker.book2D("convVtxRvsZBarrelZoom2"," Photon Reco conversion vtx position",zBinForXray, zMinForXray, zMaxForXray, rBinForXray, -10., 20.);
496  h_convVtxYvsX_zoom_[0] = iBooker.book2D("convVtxYvsXTrkBarrelZoom1"," Photon Reco conversion vtx position, (x,y) eta<1 ",100, -40., 40., 100, -40., 40.);
497  h_convVtxYvsX_zoom_[1] = iBooker.book2D("convVtxYvsXTrkBarrelZoom2"," Photon Reco conversion vtx position, (x,y) eta<1 ",100, -20., 20., 100, -20., 20.);
498 
499  h_convVtxdR_ = iBooker.book1D("convVtxdR"," Photon Reco conversion vtx dR",100, -10.,10.);
500  h_convVtxdX_ = iBooker.book1D("convVtxdX"," Photon Reco conversion vtx dX",100, -10.,10.);
501  h_convVtxdY_ = iBooker.book1D("convVtxdY"," Photon Reco conversion vtx dY",100, -10.,10.);
502  h_convVtxdZ_ = iBooker.book1D("convVtxdZ"," Photon Reco conversion vtx dZ",100, -20.,20.);
503 
504  h_convVtxdPhi_ = iBooker.book1D("convVtxdPhi"," Photon Reco conversion vtx dPhi",100, -0.01,0.01);
505  h_convVtxdEta_ = iBooker.book1D("convVtxdEta"," Photon Reco conversion vtx dEta",100, -0.5,0.5);
506 
507  h_convVtxdR_barrel_ = iBooker.book1D("convVtxdR_barrel"," Photon Reco conversion vtx dR, |eta|<=1.2",100, -10.,10.);
508  h_convVtxdX_barrel_ = iBooker.book1D("convVtxdX_barrel"," Photon Reco conversion vtx dX, |eta|<=1.2",100, -10.,10.);
509  h_convVtxdY_barrel_ = iBooker.book1D("convVtxdY_barrel"," Photon Reco conversion vtx dY, |eta|<=1.2 ",100, -10.,10.);
510  h_convVtxdZ_barrel_ = iBooker.book1D("convVtxdZ_barrel"," Photon Reco conversion vtx dZ, |eta|<=1.2,",100, -20.,20.);
511 
512  h_convVtxdR_endcap_ = iBooker.book1D("convVtxdR_endcap"," Photon Reco conversion vtx dR, |eta|>1.2 ",100, -10.,10.);
513  h_convVtxdX_endcap_ = iBooker.book1D("convVtxdX_endcap"," Photon Reco conversion vtx dX, |eta|>1.2",100, -10.,10.);
514  h_convVtxdY_endcap_ = iBooker.book1D("convVtxdY_endcap"," Photon Reco conversion vtx dY, |eta|>1.2",100, -10.,10.);
515  h_convVtxdZ_endcap_ = iBooker.book1D("convVtxdZ_endcap"," Photon Reco conversion vtx dZ, |eta|>1.2",100, -20.,20.);
516 
517 
518 
519  h2_convVtxdRVsR_ = iBooker.book2D("h2ConvVtxdRVsR"," Conversion vtx dR vsR" ,rBin,rMin, rMax,100, -20.,20.);
520  h2_convVtxdRVsEta_ = iBooker.book2D("h2ConvVtxdRVsEta","Conversion vtx dR vs Eta" ,etaBin2,etaMin, etaMax,100, -20.,20.);
521 
522  p_convVtxdRVsR_ = iBooker.bookProfile("pConvVtxdRVsR"," Conversion vtx dR vsR" ,rBin,rMin, rMax ,100, -20.,20., "");
523  p_convVtxdRVsEta_ = iBooker.bookProfile("pConvVtxdRVsEta","Conversion vtx dR vs Eta" ,etaBin2,etaMin, etaMax, 100, -20.,20., "");
524  p_convVtxdXVsX_ = iBooker.bookProfile("pConvVtxdXVsX","Conversion vtx dX vs X" ,120,-60, 60 ,100, -20.,20., "");
525  p_convVtxdYVsY_ = iBooker.bookProfile("pConvVtxdYVsY","Conversion vtx dY vs Y" ,120,-60, 60 ,100, -20.,20., "");
526  p_convVtxdZVsZ_ = iBooker.bookProfile("pConvVtxdZVsZ","Conversion vtx dZ vs Z" ,zBin,zMin,zMax ,100, -20.,20., "");
527 
528  p_convVtxdZVsR_ = iBooker.bookProfile("pConvVtxdZVsR","Conversion vtx dZ vs R" ,rBin,rMin,rMax ,100, -20.,20., "");
529  p2_convVtxdRVsRZ_ = iBooker.bookProfile2D("p2ConvVtxdRVsRZ","Conversion vtx dR vs RZ" ,zBin,zMin, zMax,rBin,rMin,rMax,100, 0.,20.,"s");
530  p2_convVtxdZVsRZ_ = iBooker.bookProfile2D("p2ConvVtxdZVsRZ","Conversion vtx dZ vs RZ" ,zBin,zMin, zMax,rBin,rMin,rMax,100, 0.,20.,"s");
531 
532 
533  histname="EoverPtracks";
534  h_EoverPTracks_[0][0] = iBooker.book1D(histname+"All"," photons conversion E/p: all Ecal ", eoverpBin, eoverpMin, eoverpMax );
535  h_EoverPTracks_[0][1] = iBooker.book1D(histname+"Barrel"," photons conversion E/p: Barrel Ecal", eoverpBin, eoverpMin, eoverpMax);
536  h_EoverPTracks_[0][2] = iBooker.book1D(histname+"Endcap"," photons conversion E/p: Endcap Ecal ", eoverpBin, eoverpMin, eoverpMax);
537  h_EoverPTracks_[1][0] = iBooker.book1D(histname+"All_Ass"," photons conversion E/p: all Ecal ", eoverpBin, eoverpMin, eoverpMax);
538  h_EoverPTracks_[1][1] = iBooker.book1D(histname+"Barrel_Ass"," photons conversion E/p: Barrel Ecal", eoverpBin, eoverpMin, eoverpMax);
539  h_EoverPTracks_[1][2] = iBooker.book1D(histname+"Endcap_Ass"," photons conversion E/p: Endcap Ecal ", eoverpBin, eoverpMin, eoverpMax);
540  h_EoverPTracks_[2][0] = iBooker.book1D(histname+"All_Fakes"," photons conversion E/p: all Ecal ", eoverpBin, eoverpMin, eoverpMax);
541  h_EoverPTracks_[2][1] = iBooker.book1D(histname+"Barrel_Fakes"," photons conversion E/p: Barrel Ecal", eoverpBin, eoverpMin, eoverpMax);
542  h_EoverPTracks_[2][2] = iBooker.book1D(histname+"Endcap_Fakes"," photons conversion E/p: Endcap Ecal ", eoverpBin, eoverpMin, eoverpMax);
543 
544 
545  h2_convVtxRrecVsTrue_ = iBooker.book2D("h2ConvVtxRrecVsTrue","Photon Reco conversion vtx R rec vs true" ,rBin,rMin, rMax,rBin,rMin, rMax);
546 
547  histname="vtxChi2Prob";
548  h_vtxChi2Prob_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 100, 0., 1.);
549  h_vtxChi2Prob_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 100, 0., 1.);
550  h_vtxChi2Prob_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 100, 0., 1.);
551  h_vtxChi2Prob_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 100, 0., 1.);
552  h_vtxChi2Prob_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 100, 0., 1.);
553  h_vtxChi2Prob_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 100, 0., 1.);
554  h_vtxChi2Prob_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 100, 0., 1.);
555  h_vtxChi2Prob_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 100, 0., 1.);
556  h_vtxChi2Prob_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 100, 0., 1.);
557 
558 
559  h_zPVFromTracks_[1] = iBooker.book1D("zPVFromTracks"," Photons: PV z from conversion tracks",100, -25., 25.);
560  h_dzPVFromTracks_[1] = iBooker.book1D("dzPVFromTracks"," Photons: PV Z_rec - Z_true from conversion tracks",100, -5., 5.);
561  h2_dzPVVsR_ = iBooker.book2D("h2dzPVVsR","Photon Reco conversions: dz(PV) vs R" ,rBin,rMin, rMax,100, -3.,3.);
562  p_dzPVVsR_ = iBooker.bookProfile("pdzPVVsR","Photon Reco conversions: dz(PV) vs R" ,rBin,rMin, rMax, 100, -3.,3.,"");
563 
564 
565  histname="lxybs";
566  h_lxybs_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 200, -100., 100.);
567  h_lxybs_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 200, -100., 100.);
568  h_lxybs_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 200, -100., 100.);
569  h_lxybs_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 200, -100., 100.);
570  h_lxybs_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 200, -100., 100.);
571  h_lxybs_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 200, -100., 100.);
572  h_lxybs_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 200, -100., 100.);
573  h_lxybs_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 200, -100., 100.);
574  h_lxybs_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 200, -100., 100.);
575 
576  histname="maxNHitsBeforeVtx";
577  h_maxNHitsBeforeVtx_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 16, -0.5, 15.5);
578  h_maxNHitsBeforeVtx_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 16, -0.5, 15.5);
579  h_maxNHitsBeforeVtx_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 16, -0.5, 15.5);
580  h_maxNHitsBeforeVtx_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 16, -0.5, 15.5);
581  h_maxNHitsBeforeVtx_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 16, -0.5, 15.5);
582  h_maxNHitsBeforeVtx_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 16, -0.5, 15.5);
583  h_maxNHitsBeforeVtx_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 16, -0.5, 15.5);
584  h_maxNHitsBeforeVtx_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 16, -0.5, 15.5);
585  h_maxNHitsBeforeVtx_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 16, -0.5, 15.5);
586 
587  histname="leadNHitsBeforeVtx";
588  h_leadNHitsBeforeVtx_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 16, -0.5, 15.5);
589  h_leadNHitsBeforeVtx_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 16, -0.5, 15.5);
590  h_leadNHitsBeforeVtx_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 16, -0.5, 15.5);
591  h_leadNHitsBeforeVtx_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 16, -0.5, 15.5);
592  h_leadNHitsBeforeVtx_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 16, -0.5, 15.5);
593  h_leadNHitsBeforeVtx_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 16, -0.5, 15.5);
594  h_leadNHitsBeforeVtx_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 16, -0.5, 15.5);
595  h_leadNHitsBeforeVtx_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 16, -0.5, 15.5);
596  h_leadNHitsBeforeVtx_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 16, -0.5, 15.5);
597 
598  histname="trailNHitsBeforeVtx";
599  h_trailNHitsBeforeVtx_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 16, -0.5, 15.5);
600  h_trailNHitsBeforeVtx_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 16, -0.5, 15.5);
601  h_trailNHitsBeforeVtx_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 16, -0.5, 15.5);
602  h_trailNHitsBeforeVtx_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 16, -0.5, 15.5);
603  h_trailNHitsBeforeVtx_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 16, -0.5, 15.5);
604  h_trailNHitsBeforeVtx_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 16, -0.5, 15.5);
605  h_trailNHitsBeforeVtx_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 16, -0.5, 15.5);
606  h_trailNHitsBeforeVtx_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 16, -0.5, 15.5);
607  h_trailNHitsBeforeVtx_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 16, -0.5, 15.5);
608 
609  histname="sumNHitsBeforeVtx";
610  h_sumNHitsBeforeVtx_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 16, -0.5, 15.5);
611  h_sumNHitsBeforeVtx_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 16, -0.5, 15.5);
612  h_sumNHitsBeforeVtx_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 16, -0.5, 15.5);
613  h_sumNHitsBeforeVtx_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 16, -0.5, 15.5);
614  h_sumNHitsBeforeVtx_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 16, -0.5, 15.5);
615  h_sumNHitsBeforeVtx_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 16, -0.5, 15.5);
616  h_sumNHitsBeforeVtx_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 16, -0.5, 15.5);
617  h_sumNHitsBeforeVtx_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 16, -0.5, 15.5);
618  h_sumNHitsBeforeVtx_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 16, -0.5, 15.5);
619 
620  histname="maxDlClosestHitToVtx";
621  h_maxDlClosestHitToVtx_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 100, -10., 10.);
622  h_maxDlClosestHitToVtx_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 100, -10., 10.);
623  h_maxDlClosestHitToVtx_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 100, -10., 10.);
624  h_maxDlClosestHitToVtx_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 100, -10., 10.);
625  h_maxDlClosestHitToVtx_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 100, -10., 10.);
626  h_maxDlClosestHitToVtx_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 100, -10., 10.);
627  h_maxDlClosestHitToVtx_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 100, -10., 10.);
628  h_maxDlClosestHitToVtx_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 100, -10., 10.);
629  h_maxDlClosestHitToVtx_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 100, -10., 10.);
630 
631  histname="maxDlClosestHitToVtxSig";
632  h_maxDlClosestHitToVtxSig_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 100, -8., 8.);
633  h_maxDlClosestHitToVtxSig_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 100, -8., 8.);
634  h_maxDlClosestHitToVtxSig_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 100, -8., 8.);
635  h_maxDlClosestHitToVtxSig_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 100, -8., 8.);
636  h_maxDlClosestHitToVtxSig_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 100, -8., 8.);
637  h_maxDlClosestHitToVtxSig_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 100, -8., 8.);
638  h_maxDlClosestHitToVtxSig_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 100, -8., 8.);
639  h_maxDlClosestHitToVtxSig_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 100, -8., 8.);
640  h_maxDlClosestHitToVtxSig_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 100, -8., 8.);
641 
642  histname="deltaExpectedHitsInner";
643  h_deltaExpectedHitsInner_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 31, -15.5, 15.5);
644  h_deltaExpectedHitsInner_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 31, -15.5, 15.5);
645  h_deltaExpectedHitsInner_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 31, -15.5, 15.5);
646  h_deltaExpectedHitsInner_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 31, -15.5, 15.5);
647  h_deltaExpectedHitsInner_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 31, -15.5, 15.5);
648  h_deltaExpectedHitsInner_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 31, -15.5, 15.5);
649  h_deltaExpectedHitsInner_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 31, -15.5, 15.5);
650  h_deltaExpectedHitsInner_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 31, -15.5, 15.5);
651  h_deltaExpectedHitsInner_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 31, -15.5, 15.5);
652 
653  histname="leadExpectedHitsInner";
654  h_leadExpectedHitsInner_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 16, -0.5, 15.5);
655  h_leadExpectedHitsInner_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 16, -0.5, 15.5);
656  h_leadExpectedHitsInner_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 16, -0.5, 15.5);
657  h_leadExpectedHitsInner_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 16, -0.5, 15.5);
658  h_leadExpectedHitsInner_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 16, -0.5, 15.5);
659  h_leadExpectedHitsInner_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 16, -0.5, 15.5);
660  h_leadExpectedHitsInner_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 16, -0.5, 15.5);
661  h_leadExpectedHitsInner_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 16, -0.5, 15.5);
662  h_leadExpectedHitsInner_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 16, -0.5, 15.5);
663 
664  histname="nSharedHits";
665  h_nSharedHits_[0][0] = iBooker.book1D(histname+"All","vertex #chi^{2} all", 16, -0.5, 15.5);
666  h_nSharedHits_[0][1] = iBooker.book1D(histname+"Barrel","vertex #chi^{2} barrel", 16, -0.5, 15.5);
667  h_nSharedHits_[0][2] = iBooker.book1D(histname+"Endcap","vertex #chi^{2} endcap", 16, -0.5, 15.5);
668  h_nSharedHits_[1][0] = iBooker.book1D(histname+"All_Ass","vertex #chi^{2} all", 16, -0.5, 15.5);
669  h_nSharedHits_[1][1] = iBooker.book1D(histname+"Barrel_Ass","vertex #chi^{2} barrel", 16, -0.5, 15.5);
670  h_nSharedHits_[1][2] = iBooker.book1D(histname+"Endcap_Ass","vertex #chi^{2} endcap", 16, -0.5, 15.5);
671  h_nSharedHits_[2][0] = iBooker.book1D(histname+"All_Fakes","vertex #chi^{2} all", 16, -0.5, 15.5);
672  h_nSharedHits_[2][1] = iBooker.book1D(histname+"Barrel_Fakes","vertex #chi^{2} barrel", 16, -0.5, 15.5);
673  h_nSharedHits_[2][2] = iBooker.book1D(histname+"Endcap_Fakes","vertex #chi^{2} endcap", 16, -0.5, 15.5);
674 
676  histname="nHits";
677  nHits_[0] = iBooker.book2D(histname+"AllTracks","Photons:Tracks from conversions: # of hits all tracks",etaBin,etaMin, etaMax,30,0., 30.);
678  nHits_[1] = iBooker.book2D(histname+"AllTracks_Ass","Photons:Tracks from conversions: # of hits all tracks ass",etaBin,etaMin, etaMax,30,0., 30.);
679  nHits_[2] = iBooker.book2D(histname+"AllTracks_Fakes","Photons:Tracks from conversions: # of hits all tracks fakes",etaBin,etaMin, etaMax,30,0., 30.);
680 
681 
682  histname="nHitsVsEta";
683  nHitsVsEta_[0] = iBooker.book2D(histname+"AllTracks","Photons:Tracks from conversions: # of hits vs #eta all tracks",etaBin,etaMin, etaMax,30,0., 30.);
684  nHitsVsEta_[1] = iBooker.book2D(histname+"AllTracks_Ass","Photons:Tracks from conversions: # of hits vs #eta all tracks",etaBin,etaMin, etaMax,30,0., 30.);
685  nHitsVsEta_[2] = iBooker.book2D(histname+"AllTracks_Fakes","Photons:Tracks from conversions: # of hits vs #eta all tracks",etaBin,etaMin, etaMax,30,0., 30.);
686  histname="h_nHitsVsEta";
687  p_nHitsVsEta_[0] = iBooker.bookProfile(histname+"AllTracks","Photons:Tracks from conversions: # of hits vs #eta all tracks",etaBin,etaMin, etaMax, 30,-0.5, 29.5,"");
688  p_nHitsVsEta_[1] = iBooker.bookProfile(histname+"AllTracks_Ass","Photons:Tracks from conversions: # of hits vs #eta all tracks",etaBin,etaMin, etaMax, 30,-0.5, 29.5,"");
689  p_nHitsVsEta_[2] = iBooker.bookProfile(histname+"AllTracks_Fakes","Photons:Tracks from conversions: # of hits vs #eta all tracks",etaBin,etaMin, etaMax, 30,-0.5, 29.5,"");
690 
691 
692  histname="nHitsVsR";
693  nHitsVsR_[0] = iBooker.book2D(histname+"AllTracks","Photons:Tracks from conversions: # of hits vs radius all tracks" ,rBin,rMin, rMax,30,0., 30.);
694  nHitsVsR_[1] = iBooker.book2D(histname+"AllTracks_Ass","Photons:Tracks from conversions: # of hits vs radius all tracks" ,rBin,rMin, rMax,30,0., 30.);
695  nHitsVsR_[2] = iBooker.book2D(histname+"AllTracks_Fakes","Photons:Tracks from conversions: # of hits vs radius all tracks" ,rBin,rMin, rMax,30,0., 30.);
696 
697  histname="h_nHitsVsR";
698  p_nHitsVsR_[0] = iBooker.bookProfile(histname+"AllTracks","Photons:Tracks from conversions: # of hits vs radius all tracks",rBin,rMin, rMax, 30,-0.5, 29.5,"");
699  p_nHitsVsR_[1] = iBooker.bookProfile(histname+"AllTracks_Ass","Photons:Tracks from conversions: # of hits vs radius all tracks",rBin,rMin, rMax, 30,-0.5, 29.5,"");
700  p_nHitsVsR_[2] = iBooker.bookProfile(histname+"AllTracks_Fakes","Photons:Tracks from conversions: # of hits vs radius all tracks",rBin,rMin, rMax, 30,-0.5, 29.5,"");
701 
702  histname="tkChi2";
703  h_tkChi2_[0] = iBooker.book1D(histname+"AllTracks","Photons:Tracks from conversions: #chi^{2} of all tracks", 100, chi2Min, chi2Max);
704  h_tkChi2_[1] = iBooker.book1D(histname+"AllTracks_Ass","Photons:Tracks from conversions: #chi^{2} of all tracks", 100, chi2Min, chi2Max);
705  h_tkChi2_[2] = iBooker.book1D(histname+"AllTracks_Fakes","Photons:Tracks from conversions: #chi^{2} of all tracks", 100, chi2Min, chi2Max);
706 
707  histname="tkChi2Large";
708  h_tkChi2Large_[0] = iBooker.book1D(histname+"AllTracks","Photons:Tracks from conversions: #chi^{2} of all tracks", 1000, 0., 5000.0);
709  h_tkChi2Large_[1] = iBooker.book1D(histname+"AllTracks_Ass","Photons:Tracks from conversions: #chi^{2} of all tracks", 1000, 0., 5000.0);
710  h_tkChi2Large_[2] = iBooker.book1D(histname+"AllTracks_Fakes","Photons:Tracks from conversions: #chi^{2} of all tracks", 1000, 0., 5000.0);
711 
712 
713  histname="h2Chi2VsEta";
714  h2_Chi2VsEta_[0]=iBooker.book2D(histname+"All"," Reco Track #chi^{2} vs #eta: All ",etaBin2,etaMin, etaMax,100, chi2Min, chi2Max);
715  h2_Chi2VsEta_[1]=iBooker.book2D(histname+"All_Ass"," Reco Track #chi^{2} vs #eta: All ",etaBin2,etaMin, etaMax,100, chi2Min, chi2Max);
716  h2_Chi2VsEta_[2]=iBooker.book2D(histname+"All_Fakes"," Reco Track #chi^{2} vs #eta: All ",etaBin2,etaMin, etaMax,100, chi2Min, chi2Max);
717  histname="pChi2VsEta";
718  p_Chi2VsEta_[0]=iBooker.bookProfile(histname+"All"," Reco Track #chi^{2} vs #eta : All ",etaBin2,etaMin, etaMax, 100, chi2Min, chi2Max,"");
719  p_Chi2VsEta_[1]=iBooker.bookProfile(histname+"All_Ass"," Reco Track #chi^{2} vs #eta : All ",etaBin2,etaMin, etaMax, 100, chi2Min, chi2Max,"");
720  p_Chi2VsEta_[2]=iBooker.bookProfile(histname+"All_Fakes"," Reco Track #chi^{2} vs #eta : All ",etaBin2,etaMin, etaMax, 100, chi2Min, chi2Max,"");
721 
722  histname="h2Chi2VsR";
723  h2_Chi2VsR_[0]=iBooker.book2D(histname+"All"," Reco Track #chi^{2} vs R: All ",rBin,rMin, rMax,100,chi2Min, chi2Max);
724  h2_Chi2VsR_[1]=iBooker.book2D(histname+"All_Ass"," Reco Track #chi^{2} vs R: All ",rBin,rMin, rMax,100,chi2Min, chi2Max);
725  h2_Chi2VsR_[2]=iBooker.book2D(histname+"All_Fakes"," Reco Track #chi^{2} vs R: All ",rBin,rMin, rMax,100,chi2Min, chi2Max);
726  histname="pChi2VsR";
727  p_Chi2VsR_[0]=iBooker.bookProfile(histname+"All"," Reco Track #chi^{2} vas R : All ",rBin,rMin,rMax, 100,chi2Min, chi2Max,"");
728  p_Chi2VsR_[1]=iBooker.bookProfile(histname+"All_Ass"," Reco Track #chi^{2} vas R : All ",rBin,rMin,rMax, 100,chi2Min, chi2Max,"");
729  p_Chi2VsR_[2]=iBooker.bookProfile(histname+"All_Fakes"," Reco Track #chi^{2} vas R : All ",rBin,rMin,rMax, 100,chi2Min, chi2Max,"");
730 
731  histname="hTkD0";
732  h_TkD0_[0]=iBooker.book1D(histname+"All"," Reco Track D0*q: All ",200,-0.1,60);
733  h_TkD0_[1]=iBooker.book1D(histname+"All_Ass"," Reco Track D0*q: Barrel ",200,-0.1,60);
734  h_TkD0_[2]=iBooker.book1D(histname+"All_Fakes"," Reco Track D0*q: Endcap ",200,-0.1,60);
735 
736 
737 
738  histname="hTkPtPull";
739  h_TkPtPull_[0]=iBooker.book1D(histname+"All"," Reco Track Pt pull: All ",100, -20., 10.);
740  histname="hTkPtPull";
741  h_TkPtPull_[1]=iBooker.book1D(histname+"Barrel"," Reco Track Pt pull: Barrel ",100, -20., 10.);
742  histname="hTkPtPull";
743  h_TkPtPull_[2]=iBooker.book1D(histname+"Endcap"," Reco Track Pt pull: Endcap ",100, -20., 10.);
744 
745  histname="h2TkPtPullEta";
746  h2_TkPtPull_[0]=iBooker.book2D(histname+"All"," Reco Track Pt pull: All ",etaBin2,etaMin, etaMax,100, -20., 10.);
747  histname="pTkPtPullEta";
748  p_TkPtPull_[0]=iBooker.bookProfile(histname+"All"," Reco Track Pt pull: All ",etaBin2,etaMin, etaMax, 100, -20., 10., " ");
749 
750 
751  histname="PtRecVsPtSim";
752  h2_PtRecVsPtSim_[0]=iBooker.book2D(histname+"All", "Pt Rec vs Pt sim: All ", etBin,etMin,etMax,etBin,etMin, etMax);
753  h2_PtRecVsPtSim_[1]=iBooker.book2D(histname+"Barrel", "Pt Rec vs Pt sim: Barrel ", etBin,etMin,etMax,etBin,etMin, etMax);
754  h2_PtRecVsPtSim_[2]=iBooker.book2D(histname+"Endcap", "Pt Rec vs Pt sim: Endcap ", etBin,etMin,etMax,etBin,etMin, etMax);
755 
756  histname="photonPtRecVsPtSim";
757  h2_photonPtRecVsPtSim_=iBooker.book2D(histname+"All", "Pt Rec vs Pt sim: All ", etBin,etMin,etMax,etBin,etMin, etMax);
758 
759  histname="nHitsBeforeVtx";
760  h_nHitsBeforeVtx_[0]=iBooker.book1D(histname+"All", "Pt Rec vs Pt sim: All ", 16, -0.5, 15.5);
761  h_nHitsBeforeVtx_[1]=iBooker.book1D(histname+"Barrel", "Pt Rec vs Pt sim: Barrel ", 16, -0.5, 15.5);
762  h_nHitsBeforeVtx_[2]=iBooker.book1D(histname+"Endcap", "Pt Rec vs Pt sim: Endcap ", 16, -0.5, 15.5);
763 
764  histname="dlClosestHitToVtx";
765  h_dlClosestHitToVtx_[0]=iBooker.book1D(histname+"All", "Pt Rec vs Pt sim: All ", 100, -10., 10.);
766  h_dlClosestHitToVtx_[1]=iBooker.book1D(histname+"Barrel", "Pt Rec vs Pt sim: Barrel ", 100, -10., 10.);
767  h_dlClosestHitToVtx_[2]=iBooker.book1D(histname+"Endcap", "Pt Rec vs Pt sim: Endcap ", 100, -10., 10.);
768 
769  histname="dlClosestHitToVtxSig";
770  h_dlClosestHitToVtxSig_[0]=iBooker.book1D(histname+"All", "Pt Rec vs Pt sim: All ", 100, -8., 8.);
771  h_dlClosestHitToVtxSig_[1]=iBooker.book1D(histname+"Barrel", "Pt Rec vs Pt sim: Barrel ", 100, -8., 8.);
772  h_dlClosestHitToVtxSig_[2]=iBooker.book1D(histname+"Endcap", "Pt Rec vs Pt sim: Endcap ", 100, -8., 8.);
773 
774  h_match_= iBooker.book1D("h_match"," ", 3, -0.5,2.5);
775 
776 
777  } // if DQM
778 
779 
780 
781 }
782 
783 
784 
785  void TkConvValidator::dqmBeginRun (edm::Run const & r, edm::EventSetup const & theEventSetup) {
786 
787  //get magnetic field
788  edm::LogInfo("ConvertedPhotonProducer") << " get magnetic field" << "\n";
789  theEventSetup.get<IdealMagneticFieldRecord>().get(theMF_);
790 
791  thePhotonMCTruthFinder_ = new PhotonMCTruthFinder();
792 
793 }
794 
795 void TkConvValidator::endRun (edm::Run const& r, edm::EventSetup const & theEventSetup) {
796 
797  delete thePhotonMCTruthFinder_;
798 
799 }
800 
801 
802 
803 void TkConvValidator::analyze( const edm::Event& e, const edm::EventSetup& esup ) {
804  thePhotonMCTruthFinder_->clear();
805  using namespace edm;
806  // const float etaPhiDistance=0.01;
807  // Fiducial region
808  // const float TRK_BARL =0.9;
809  const float BARL = 1.4442; // DAQ TDR p.290
810  // const float END_LO = 1.566; // unused
811  const float END_HI = 2.5;
812  // Electron mass
813  // const Float_t mElec= 0.000511; // unused
814 
815  Handle<reco::TrackToTrackingParticleAssociator> theTrackAssociator;
816  e.getByToken(trackAssociator_Token_,theTrackAssociator);
817 
818 
819  nEvt_++;
820  LogInfo("TkConvValidator") << "TkConvValidator Analyzing event number: " << e.id() << " Global Counter " << nEvt_ <<"\n";
821  // std::cout << "TkConvValidator Analyzing event number: " << e.id() << " Global Counter " << nEvt_ <<"\n";
822 
823 
824  // get the geometry from the event setup:
825  esup.get<CaloGeometryRecord>().get(theCaloGeom_);
826 
827 
828  // Transform Track into TransientTrack (needed by the Vertex fitter)
829  edm::ESHandle<TransientTrackBuilder> theTTB;
830  esup.get<TransientTrackRecord>().get("TransientTrackBuilder",theTTB);
831 
832 
834  Handle<reco::ConversionCollection> convHandle;
835  e.getByToken(conversionCollectionPr_Token_, convHandle);
836  const reco::ConversionCollection convCollection = *(convHandle.product());
837  if (!convHandle.isValid()) {
838  edm::LogError("ConversionsProducer") << "Error! Can't get the collection "<< std::endl;
839  return;
840  }
841 
843  Handle<reco::PhotonCollection> photonHandle;
844  e.getByToken(photonCollectionPr_Token_, photonHandle);
845  const reco::PhotonCollection photonCollection = *(photonHandle.product());
846  if (!photonHandle.isValid()) {
847  edm::LogError("PhotonProducer") << "Error! Can't get the Photon collection "<< std::endl;
848  return;
849  }
850 
851 
852  // offline Primary vertex
853  edm::Handle<reco::VertexCollection> vertexHandle;
854  reco::VertexCollection vertexCollection;
855  e.getByToken(offline_pvToken_, vertexHandle);
856  if (!vertexHandle.isValid()) {
857  edm::LogError("TrackerOnlyConversionProducer") << "Error! Can't get the product primary Vertex Collection "<< "\n";
858  } else {
859  vertexCollection = *(vertexHandle.product());
860  }
861  reco::Vertex the_pvtx;
862  bool valid_pvtx = false;
863  if (!vertexCollection.empty()){
864  the_pvtx = *(vertexCollection.begin());
865  //asking for one good vertex
866  if (the_pvtx.isValid() && fabs(the_pvtx.position().z())<=15 && the_pvtx.position().Rho()<=2){
867  valid_pvtx = true;
868  }
869  }
870 
871  edm::Handle<reco::BeamSpot> bsHandle;
872  e.getByToken(beamspotToken_, bsHandle);
873  if (!bsHandle.isValid()) {
874  edm::LogError("TrackerOnlyConversionProducer")
875  << "Error! Can't get the product primary Vertex Collection "<< "\n";
876  return;
877  }
878  const reco::BeamSpot &thebs = *bsHandle.product();
879 
880  //get tracker geometry for hits positions
881  edm::ESHandle<TrackerGeometry> tracker;
882  esup.get<TrackerDigiGeometryRecord>().get(tracker);
883  const TrackerGeometry* trackerGeom = tracker.product();
884 
885 
886 
888  //get simtrack info
889  //std::vector<SimTrack> theSimTracks;
890  //std::vector<SimVertex> theSimVertices;
891 
892  edm::Handle<SimTrackContainer> SimTk;
893  edm::Handle<SimVertexContainer> SimVtx;
894  e.getByToken(g4_simTk_Token_, SimTk);
895  e.getByToken(g4_simVtx_Token_, SimVtx);
896 
897  bool useTP = parameters_.getParameter<bool>("useTP");
898  TrackingParticleRefVector dummy;
899  edm::Handle<TrackingParticleRefVector> TPHandleForEff;
900  edm::Handle<TrackingParticleRefVector> TPHandleForFakeRate;
901  if (useTP) {
902  e.getByToken(tpSelForEff_Token_, TPHandleForEff);
903  e.getByToken(tpSelForFake_Token_, TPHandleForFakeRate);
904  }
905 
906  const TrackingParticleRefVector &tpForEfficiency= useTP? *(TPHandleForEff.product()) : dummy;
907  const TrackingParticleRefVector &tpForFakeRate= useTP? *(TPHandleForFakeRate.product()):dummy;
908 
909  const std::vector<SimTrack> &theSimTracks= *SimTk;
910  const std::vector<SimVertex> &theSimVertices= *SimVtx;
911 
912  //theSimTracks.insert(theSimTracks.end(),SimTk->begin(),SimTk->end());
913  //theSimVertices.insert(theSimVertices.end(),SimVtx->begin(),SimVtx->end());
914  std::vector<PhotonMCTruth> mcPhotons=thePhotonMCTruthFinder_->find (theSimTracks, theSimVertices);
915 
916  //edm::Handle<edm::HepMCProduct> hepMC;
917  //e.getByToken(hepMC_Token_, hepMC);
918  // const HepMC::GenEvent *myGenEvent = hepMC->GetEvent(); // unused
919 
920 
921  // get generated jets
922  //edm::Handle<reco::GenJetCollection> GenJetsHandle;
923  //e.getByToken(genjets_Token_, GenJetsHandle);
924  //const reco::GenJetCollection &genJetCollection = *(GenJetsHandle.product());
925 
926  ConversionHitChecker hitChecker;
927 
928  // ################ SIM to RECO ######################### //
929  std::map<const reco::Track*,TrackingParticleRef> myAss;
930  std::map<const reco::Track*,TrackingParticleRef>::const_iterator itAss;
931 
932  for ( std::vector<PhotonMCTruth>::const_iterator mcPho=mcPhotons.begin(); mcPho !=mcPhotons.end(); mcPho++) {
933 
934  mcConvPt_= (*mcPho).fourMomentum().et();
935  float mcPhi= (*mcPho).fourMomentum().phi();
936  mcPhi_= phiNormalization(mcPhi);
937  mcEta_= (*mcPho).fourMomentum().pseudoRapidity();
938  mcEta_ = etaTransformation(mcEta_, (*mcPho).primaryVertex().z() );
939  mcConvR_= (*mcPho).vertex().perp();
940  mcConvX_= (*mcPho).vertex().x();
941  mcConvY_= (*mcPho).vertex().y();
942  mcConvZ_= (*mcPho).vertex().z();
943  mcConvEta_= (*mcPho).vertex().eta();
944  mcConvPhi_= (*mcPho).vertex().phi();
945 
946  if ( fabs(mcEta_) > END_HI ) continue;
947 
948  if (mcConvPt_<minPhoPtForEffic) continue;
949  if (fabs(mcEta_)>maxPhoEtaForEffic) continue;
950  if (fabs(mcConvZ_)>maxPhoZForEffic) continue;
951  if (mcConvR_>maxPhoRForEffic) continue;
953 
954  bool goodSimConversion=false;
955  bool visibleConversion=false;
956  bool visibleConversionsWithTwoSimTracks=false;
957  if ( (*mcPho).isAConversion() == 1 ) {
958  nSimConv_[0]++;
959  h_AllSimConv_[0]->Fill( mcEta_ ) ;
960  h_AllSimConv_[1]->Fill( mcPhi_ );
961  h_AllSimConv_[2]->Fill( mcConvR_ );
962  h_AllSimConv_[3]->Fill( mcConvZ_ );
963  h_AllSimConv_[4]->Fill( (*mcPho).fourMomentum().et());
964 
965  if ( mcConvR_ <15) h_SimConvEtaPix_[0]->Fill( mcEta_ ) ;
966 
967  if ( ( fabs(mcEta_) <= BARL && mcConvR_ <85 ) ||
968  ( fabs(mcEta_) > BARL && fabs(mcEta_) <=END_HI && fabs( (*mcPho).vertex().z() ) < 210 ) ) visibleConversion=true;
969 
970  theConvTP_.clear();
971  // std::cout << " TkConvValidator TrackingParticles TrackingParticleCollection size "<< trackingParticles.size() << "\n";
972  //duplicated TP collections for two associations
973  for(const TrackingParticleRef tp: tpForEfficiency) {
974  if ( fabs( tp->vx() - (*mcPho).vertex().x() ) < 0.0001 &&
975  fabs( tp->vy() - (*mcPho).vertex().y() ) < 0.0001 &&
976  fabs( tp->vz() - (*mcPho).vertex().z() ) < 0.0001) {
977  theConvTP_.push_back( tp );
978  }
979  }
980  //std::cout << " TkConvValidator theConvTP_ size " << theConvTP_.size() << std::endl;
981 
982  if ( theConvTP_.size() == 2 ) visibleConversionsWithTwoSimTracks=true;
983  goodSimConversion=false;
984 
985  if ( visibleConversion && visibleConversionsWithTwoSimTracks ) goodSimConversion=true;
986  if ( goodSimConversion ) {
987  nSimConv_[1]++;
988  h_VisSimConv_[0]->Fill( mcEta_ ) ;
989  h_VisSimConv_[1]->Fill( mcPhi_ );
990  h_VisSimConv_[2]->Fill( mcConvR_ );
991  h_VisSimConv_[3]->Fill( mcConvZ_ );
992  h_VisSimConv_[4]->Fill( (*mcPho).fourMomentum().et());
993 
994  }
995 
996  for ( edm::RefVector<TrackingParticleCollection>::iterator iTrk=theConvTP_.begin(); iTrk!=theConvTP_.end(); ++iTrk) {
997  h_simTkPt_ -> Fill ( (*iTrk)->pt() );
998  h_simTkEta_ -> Fill ( (*iTrk)->eta() );
999  }
1000 
1001 
1002  }
1003 
1004  if ( ! (visibleConversion && visibleConversionsWithTwoSimTracks ) ) continue;
1005 
1006  h_simConvVtxRvsZ_[0] ->Fill ( fabs (mcConvZ_), mcConvR_ ) ;
1007  if ( fabs(mcEta_) <=1.) {
1008  h_simConvVtxRvsZ_[1] ->Fill ( fabs (mcConvZ_), mcConvR_ ) ;
1009  h_simConvVtxYvsX_ ->Fill ( mcConvX_, mcConvY_ ) ;
1010  }
1011  else
1012  h_simConvVtxRvsZ_[2] ->Fill ( fabs (mcConvZ_), mcConvR_ ) ;
1013 
1014  //std::cout << " TkConvValidator theConvTP_ size " << theConvTP_.size() << std::endl;
1015  for ( edm::RefVector<TrackingParticleCollection>::iterator iTP= theConvTP_.begin(); iTP!=theConvTP_.end(); iTP++)
1016  {
1017  // std::cout << " SIM to RECO TP vertex " << (*iTP)->vx() << " " << (*iTP)->vy() << " " << (*iTP)->vz() << " pt " << (*iTP)->pt() << std::endl;
1018  }
1019 
1020  bool recomatch = false;
1021  float chi2Prob = 0.;
1023  // cout << " size of conversions " << convHandle->size() << endl;
1024  for (reco::ConversionCollection::const_iterator conv = convHandle->begin();conv!=convHandle->end();++conv) {
1025 
1026  const reco::Conversion aConv = (*conv);
1027  if ( arbitratedMerged_ && !aConv.quality(reco::Conversion::arbitratedMerged) ) continue;
1028  if ( generalTracksOnly_ && !aConv.quality(reco::Conversion::generalTracksOnly) ) continue;
1029  if ( arbitratedEcalSeeded_ && !aConv.quality(reco::Conversion::arbitratedEcalSeeded) ) continue;
1030 
1031 
1032  if ( highPurity_ && !aConv.quality(reco::Conversion::highPurity) ) continue;
1033 
1034  //problematic?
1035  std::vector<edm::RefToBase<reco::Track> > tracks = aConv.tracks();
1036 
1037 
1038  const reco::Vertex& vtx = aConv.conversionVertex();
1039  //requires two tracks and a valid vertex
1040  if (tracks.size() !=2 || !(vtx.isValid())) continue;
1041 
1042 
1043  if (ChiSquaredProbability( aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof() ) <= minProb_) continue;
1044  if (aConv.nHitsBeforeVtx().size()>1 && max(aConv.nHitsBeforeVtx().at(0),aConv.nHitsBeforeVtx().at(1)) > maxHitsBeforeVtx_ ) continue;
1045 
1046 
1047  //compute transverse decay length with respect to beamspot
1048  math::XYZVectorF themom = aConv.refittedPairMomentum();
1049  double dbsx = aConv.conversionVertex().x() - thebs.x0();
1050  double dbsy = aConv.conversionVertex().y() - thebs.y0();
1051  double lxy = (themom.x()*dbsx + themom.y()*dbsy)/themom.rho();
1052 
1053  if (lxy<minLxy_) continue;
1054 
1055  // bool phoIsInBarrel=false; // unused
1056  // bool phoIsInEndcap=false; // unused
1057  RefToBase<reco::Track> tfrb1 = aConv.tracks().front();
1058  RefToBase<reco::Track> tfrb2 = aConv.tracks().back();
1059 
1060  if ( ecalalgotracks_ && ( !(tfrb1->algo()==reco::TrackBase::outInEcalSeededConv || tfrb1->algo()==reco::TrackBase::inOutEcalSeededConv) || !(tfrb2->algo()==reco::TrackBase::outInEcalSeededConv || tfrb2->algo()==reco::TrackBase::inOutEcalSeededConv) ) ) continue;
1061 
1062 
1063  //reco::TrackRef tk1 = aConv.tracks().front();
1064  //reco::TrackRef tk2 = aConv.tracks().back();
1065  //std::cout << "SIM to RECO conversion track pt " << tk1->pt() << " " << tk2->pt() << endl;
1066  //
1067  //Use two RefToBaseVector and do two association actions to avoid that if two tracks from different collection
1068  RefToBaseVector<reco::Track> tc1, tc2;
1069  tc1.push_back(tfrb1);
1070  tc2.push_back(tfrb2);
1071  bool isAssociated = false;
1072  reco::SimToRecoCollection q1 = theTrackAssociator->associateSimToReco(tc1,theConvTP_);
1073  reco::SimToRecoCollection q2 = theTrackAssociator->associateSimToReco(tc2,theConvTP_);
1074  //try {
1075  std::vector<std::pair<RefToBase<reco::Track>, double> > trackV1, trackV2;
1076 
1077  int tp_1 = 0, tp_2 = 1;//the index of associated tp in theConvTP_ for two tracks
1078  if (q1.find(theConvTP_[0])!=q1.end()){
1079  trackV1 = (std::vector<std::pair<RefToBase<reco::Track>, double> >) q1[theConvTP_[0]];
1080  } else if (q1.find(theConvTP_[1])!=q1.end()){
1081  trackV1 = (std::vector<std::pair<RefToBase<reco::Track>, double> >) q1[theConvTP_[1]];
1082  tp_1 = 1;
1083  }
1084  if (q2.find(theConvTP_[1])!=q2.end()){
1085  trackV2 = (std::vector<std::pair<RefToBase<reco::Track>, double> >) q2[theConvTP_[1]];
1086  } else if (q2.find(theConvTP_[0])!=q2.end()){
1087  trackV2 = (std::vector<std::pair<RefToBase<reco::Track>, double> >) q2[theConvTP_[0]];
1088  tp_2 = 0;
1089  }
1090  if (!(trackV1.size()&&trackV2.size()))
1091  continue;
1092  if (tp_1 == tp_2) continue;
1093 
1094  edm::RefToBase<reco::Track> tr1 = trackV1.front().first;
1095  edm::RefToBase<reco::Track> tr2 = trackV2.front().first;
1096  //std::cout << "associated tp1 " <<theConvTP_[0]->pt() << " to track with pT=" << tr1->pt() << " " << (tr1.get())->pt() << endl;
1097  //std::cout << "associated tp2 " <<theConvTP_[1]->pt() << " to track with pT=" << tr2->pt() << " " << (tr2.get())->pt() << endl;
1098  myAss.insert( std::make_pair (tr1.get(),theConvTP_[tp_1] ) );
1099  myAss.insert( std::make_pair (tr2.get(),theConvTP_[tp_2]) );
1100 
1101  //} catch (Exception event) {
1102  //cout << "continue: " << event.what() << endl;
1103  // continue;
1104  //}
1105 
1106 
1107  isAssociated = true;
1108  recomatch = true;
1109  chi2Prob = ChiSquaredProbability( aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof() );
1110 
1111  if (isAssociated) {
1112  h_SimRecConvTwoMTracks_[0]->Fill( mcEta_ ) ;
1113  h_SimRecConvTwoMTracks_[1]->Fill( mcPhi_ );
1114  h_SimRecConvTwoMTracks_[2]->Fill( mcConvR_ );
1115  h_SimRecConvTwoMTracks_[3]->Fill( mcConvZ_ );
1116  h_SimRecConvTwoMTracks_[4]->Fill( (*mcPho).fourMomentum().et());
1117  }
1118 
1119  // break;
1120  } // loop over reco conversions
1121  if (recomatch) {
1123  h_SimConvTwoMTracks_[0]->Fill( mcEta_ ) ;
1124  h_SimConvTwoMTracks_[1]->Fill( mcPhi_ );
1125  h_SimConvTwoMTracks_[2]->Fill( mcConvR_ );
1126  h_SimConvTwoMTracks_[3]->Fill( mcConvZ_ );
1127  h_SimConvTwoMTracks_[4]->Fill( (*mcPho).fourMomentum().et());
1128 
1129 
1130  if ( chi2Prob > 0) {
1131  h_SimConvTwoMTracksAndVtxPGT0_[0]->Fill( mcEta_ ) ;
1132  h_SimConvTwoMTracksAndVtxPGT0_[1]->Fill( mcPhi_ );
1133  h_SimConvTwoMTracksAndVtxPGT0_[2]->Fill( mcConvR_ );
1134  h_SimConvTwoMTracksAndVtxPGT0_[3]->Fill( mcConvZ_ );
1135  h_SimConvTwoMTracksAndVtxPGT0_[4]->Fill( (*mcPho).fourMomentum().et());
1136  }
1137  if ( chi2Prob > 0.0005) {
1138  h_SimConvTwoMTracksAndVtxPGT0005_[0]->Fill( mcEta_ ) ;
1139  h_SimConvTwoMTracksAndVtxPGT0005_[1]->Fill( mcPhi_ );
1140  h_SimConvTwoMTracksAndVtxPGT0005_[2]->Fill( mcConvR_ );
1141  h_SimConvTwoMTracksAndVtxPGT0005_[3]->Fill( mcConvZ_ );
1142  h_SimConvTwoMTracksAndVtxPGT0005_[4]->Fill( (*mcPho).fourMomentum().et());
1143 
1144  }
1145  }
1146 
1147  } //End loop over simulated conversions
1148 
1149 
1150  // ########################### RECO to SIM ############################## //
1151 
1152  for (reco::ConversionCollection::const_iterator conv = convHandle->begin();conv!=convHandle->end();++conv) {
1153  const reco::Conversion aConv = (*conv);
1154  if ( arbitratedMerged_ && !aConv.quality(reco::Conversion::arbitratedMerged) ) continue;
1155  if ( generalTracksOnly_ && !aConv.quality(reco::Conversion::generalTracksOnly) ) continue;
1156  if ( arbitratedEcalSeeded_ && !aConv.quality(reco::Conversion::arbitratedEcalSeeded) ) continue;
1157 
1158 
1159  if ( highPurity_ && !aConv.quality(reco::Conversion::highPurity) ) continue;
1160 
1161  //problematic?
1162  std::vector<edm::RefToBase<reco::Track> > tracks = aConv.tracks();
1163 
1164  const reco::Vertex& vtx = aConv.conversionVertex();
1165  //requires two tracks and a valid vertex
1166  if (tracks.size() !=2 || !(vtx.isValid())) continue;
1167  //if (tracks.size() !=2) continue;
1168 
1169 
1170  if (ChiSquaredProbability( aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof() ) <= minProb_) continue;
1171  if (aConv.nHitsBeforeVtx().size()>1 && max(aConv.nHitsBeforeVtx().at(0),aConv.nHitsBeforeVtx().at(1)) > maxHitsBeforeVtx_ ) continue;
1172 
1173  //compute transverse decay length with respect to beamspot
1174  math::XYZVectorF themom = aConv.refittedPairMomentum();
1175  double dbsx = aConv.conversionVertex().x() - thebs.x0();
1176  double dbsy = aConv.conversionVertex().y() - thebs.y0();
1177  double lxy = (themom.x()*dbsx + themom.y()*dbsy)/themom.rho();
1178 
1179  if (lxy<minLxy_) continue;
1180 
1181  bool phoIsInBarrel=false;
1182  bool phoIsInEndcap=false;
1183  RefToBase<reco::Track> tk1 = aConv.tracks().front();
1184  RefToBase<reco::Track> tk2 = aConv.tracks().back();
1185  RefToBaseVector<reco::Track> tc1, tc2;
1186  tc1.push_back(tk1);
1187  tc2.push_back(tk2);
1188 
1189  if ( ecalalgotracks_ && ( !(tk1->algo()==reco::TrackBase::outInEcalSeededConv || tk1->algo()==reco::TrackBase::inOutEcalSeededConv) || !(tk2->algo()==reco::TrackBase::outInEcalSeededConv || tk2->algo()==reco::TrackBase::inOutEcalSeededConv) ) ) continue;
1190 
1191 
1192  //std::cout << " RECO to SIM conversion track pt " << tk1->pt() << " " << tk2->pt() << endl;
1193  const reco::Track refTk1 = aConv.conversionVertex().refittedTracks().front();
1194  const reco::Track refTk2 = aConv.conversionVertex().refittedTracks().back();
1195 
1196  //TODO replace it with phi at vertex
1197  float dPhiTracksAtVtx = aConv.dPhiTracksAtVtx();
1198  // override with the phi calculated at the vertex
1199  math::XYZVector p1AtVtx= recalculateMomentumAtFittedVertex ( (*theMF_), *trackerGeom, tk1, aConv.conversionVertex() );
1200  math::XYZVector p2AtVtx= recalculateMomentumAtFittedVertex ( (*theMF_), *trackerGeom, tk2, aConv.conversionVertex() );
1201  if ( sqrt(p1AtVtx.perp2()) > sqrt(p2AtVtx.perp2()) )
1202  dPhiTracksAtVtx = p1AtVtx.phi() - p2AtVtx.phi();
1203  else
1204  dPhiTracksAtVtx = p2AtVtx.phi() - p1AtVtx.phi();
1205 
1206 
1207  math::XYZVectorF refittedMom = aConv.refittedPairMomentum();
1208 
1209 
1210  if (fabs(refittedMom.eta())< 1.479 ) {
1211  phoIsInBarrel=true;
1212  } else {
1213  phoIsInEndcap=true;
1214  }
1215 
1216  nRecConv_++;
1217 
1218  // check matching with reco photon
1219  double Mindeltaeta = 999999;
1220  double Mindeltaphi = 999999;
1221  bool matchConvSC=false;
1222  reco::PhotonCollection::const_iterator iMatchingSC;
1223  for( reco::PhotonCollection::const_iterator iPho = photonCollection.begin(); iPho != photonCollection.end(); iPho++) {
1224  reco::Photon aPho = reco::Photon(*iPho);
1225  const double deltaphi= reco::deltaPhi( aConv.refittedPairMomentum().phi(), aPho.superCluster()->position().phi());
1226  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(),aConv.zOfPrimaryVertexFromTracks());
1227  double deltaeta = abs( aPho.superCluster()->position().eta() -ConvEta);
1228  if (abs(deltaeta)<abs(Mindeltaeta) && abs(deltaphi)<abs(Mindeltaphi)) {
1229  Mindeltaphi=abs(deltaphi);
1230  Mindeltaeta=abs(deltaeta);
1231  iMatchingSC = iPho ;
1232  }
1233  }
1234  if (abs(Mindeltaeta)<0.1 && abs(Mindeltaphi)<0.1) {
1235  matchConvSC=true;
1236  }
1237 
1238 
1240  int match =0;
1241  float invM=aConv.pairInvariantMass();
1242  float chi2Prob = ChiSquaredProbability( aConv.conversionVertex().chi2(), aConv.conversionVertex().ndof() );
1243  uint maxNHitsBeforeVtx = aConv.nHitsBeforeVtx().size()>1 ? max(aConv.nHitsBeforeVtx().at(0),aConv.nHitsBeforeVtx().at(1)) : 0;
1244  uint sumNHitsBeforeVtx = aConv.nHitsBeforeVtx().size()>1 ? aConv.nHitsBeforeVtx().at(0) + aConv.nHitsBeforeVtx().at(1) : 0;
1245  float maxDlClosestHitToVtx = aConv.dlClosestHitToVtx().size()>1 ? max(aConv.dlClosestHitToVtx().at(0).value(),aConv.dlClosestHitToVtx().at(1).value()) : 0;
1246  float maxDlClosestHitToVtxSig = aConv.dlClosestHitToVtx().size()>1 ? max(aConv.dlClosestHitToVtx().at(0).value()/aConv.dlClosestHitToVtx().at(0).error(),aConv.dlClosestHitToVtx().at(1).value()/aConv.dlClosestHitToVtx().at(1).error()) : 0;
1247 
1248  int ilead = 0, itrail = 1;
1249  if (tk2->pt() > tk1->pt()) {
1250  ilead = 1;
1251  itrail = 0;
1252  }
1253  RefToBase<reco::Track> tklead = aConv.tracks().at(ilead);
1254  RefToBase<reco::Track> tktrail = aConv.tracks().at(itrail);
1255 
1256  int deltaExpectedHitsInner = tklead->hitPattern().numberOfHits(reco::HitPattern::MISSING_INNER_HITS)
1257  - tktrail->hitPattern().numberOfHits(reco::HitPattern::MISSING_INNER_HITS);
1258  int leadExpectedHitsInner = tklead->hitPattern().numberOfHits(reco::HitPattern::MISSING_INNER_HITS);
1259  uint leadNHitsBeforeVtx = aConv.nHitsBeforeVtx().size()>1 ? aConv.nHitsBeforeVtx().at(ilead) : 0;
1260  uint trailNHitsBeforeVtx = aConv.nHitsBeforeVtx().size()>1 ? aConv.nHitsBeforeVtx().at(itrail) : 0;
1261 
1262 
1263  h_convEta_[match][0]->Fill( refittedMom.eta() );
1264  h_convEta2_[match][0]->Fill( refittedMom.eta() );
1265 
1266  h_convPhi_[match][0]->Fill( refittedMom.phi() );
1267  h_convR_[match][0]->Fill( sqrt(aConv.conversionVertex().position().perp2()) );
1268  h_convRplot_->Fill( sqrt(aConv.conversionVertex().position().perp2()) );
1269  h_convZ_[match][0]->Fill( aConv.conversionVertex().position().z() );
1270  h_convZplot_->Fill( aConv.conversionVertex().position().z() );
1271  h_convPt_[match][0]->Fill( sqrt(refittedMom.perp2()) );
1272  h_invMass_[match][0] ->Fill( invM);
1273  h_vtxChi2Prob_[match][0] ->Fill (chi2Prob);
1274  h_lxybs_[match][0] ->Fill (lxy);
1275  h_maxNHitsBeforeVtx_[match][0] ->Fill (maxNHitsBeforeVtx);
1276  h_leadNHitsBeforeVtx_[match][0] ->Fill (leadNHitsBeforeVtx);
1277  h_trailNHitsBeforeVtx_[match][0] ->Fill (trailNHitsBeforeVtx);
1278  h_sumNHitsBeforeVtx_[match][0] ->Fill (sumNHitsBeforeVtx);
1279  h_deltaExpectedHitsInner_[match][0] ->Fill (deltaExpectedHitsInner);
1280  h_leadExpectedHitsInner_[match][0] ->Fill (leadExpectedHitsInner);
1281  h_maxDlClosestHitToVtx_[match][0] ->Fill (maxDlClosestHitToVtx);
1282  h_maxDlClosestHitToVtxSig_[match][0] ->Fill (maxDlClosestHitToVtxSig);
1283  h_nSharedHits_[match][0] ->Fill (aConv.nSharedHits());
1284 
1285 
1286  if ( matchConvSC ) {
1287  h_convEtaMatchSC_[match][0]->Fill( refittedMom.eta() );
1288  h_EoverPTracks_[match][0] ->Fill (iMatchingSC->superCluster()->energy()/sqrt(refittedMom.mag2()));
1289  h_convSCdPhi_[match][0]->Fill( iMatchingSC->superCluster()->position().phi() - refittedMom.phi() );
1290  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(),aConv.zOfPrimaryVertexFromTracks());
1291  h_convSCdEta_[match][0]->Fill( iMatchingSC->superCluster()->position().eta() - ConvEta );
1292  }
1293 
1294  h_distMinAppTracks_[match][0] ->Fill (aConv.distOfMinimumApproach());
1295  h_DPhiTracksAtVtx_[match][0]->Fill( dPhiTracksAtVtx);
1296  h2_DPhiTracksAtVtxVsEta_->Fill( mcEta_, dPhiTracksAtVtx);
1297  h2_DPhiTracksAtVtxVsR_->Fill( mcConvR_, dPhiTracksAtVtx);
1298  p_DPhiTracksAtVtxVsEta_->Fill( mcEta_, dPhiTracksAtVtx);
1299  p_DPhiTracksAtVtxVsR_->Fill( mcConvR_, dPhiTracksAtVtx);
1300 
1301  h_DCotTracks_[match][0] ->Fill ( aConv.pairCotThetaSeparation() );
1302  h2_DCotTracksVsEta_->Fill( mcEta_, aConv.pairCotThetaSeparation() );
1303  h2_DCotTracksVsR_->Fill( mcConvR_, aConv.pairCotThetaSeparation() );
1304  p_DCotTracksVsEta_->Fill( mcEta_, aConv.pairCotThetaSeparation() );
1305  p_DCotTracksVsR_->Fill( mcConvR_, aConv.pairCotThetaSeparation() );
1306 
1307  if ( phoIsInBarrel ) {
1308  h_invMass_[match][1] ->Fill(invM);
1309  h_vtxChi2Prob_[match][1] ->Fill (chi2Prob);
1310  h_distMinAppTracks_[match][1] ->Fill (aConv.distOfMinimumApproach());
1311  h_DPhiTracksAtVtx_[match][1]->Fill( dPhiTracksAtVtx);
1312  h_DCotTracks_[match][1] ->Fill ( aConv.pairCotThetaSeparation() );
1313  h_lxybs_[match][1] ->Fill (lxy);
1314  h_maxNHitsBeforeVtx_[match][1] ->Fill (maxNHitsBeforeVtx);
1315  h_leadNHitsBeforeVtx_[match][1] ->Fill (leadNHitsBeforeVtx);
1316  h_trailNHitsBeforeVtx_[match][1] ->Fill (trailNHitsBeforeVtx);
1317  h_sumNHitsBeforeVtx_[match][1] ->Fill (sumNHitsBeforeVtx);
1318  h_deltaExpectedHitsInner_[match][1] ->Fill (deltaExpectedHitsInner);
1319  h_leadExpectedHitsInner_[match][1] ->Fill (leadExpectedHitsInner);
1320  h_maxDlClosestHitToVtx_[match][1] ->Fill (maxDlClosestHitToVtx);
1321  h_maxDlClosestHitToVtxSig_[match][1] ->Fill (maxDlClosestHitToVtxSig);
1322  h_nSharedHits_[match][1] ->Fill (aConv.nSharedHits());
1323 
1324  /*
1325  if ( aConv.caloCluster().size() ) {
1326  h_convSCdPhi_[match][1]->Fill( aConv.caloCluster()[0]->phi() - refittedMom.phi() );
1327  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(),aConv.zOfPrimaryVertexFromTracks());
1328  h_convSCdEta_[match][1]->Fill( aConv.caloCluster()[0]->eta() - ConvEta );
1329  }
1330  */
1331 
1332  if ( matchConvSC ) {
1333  h_EoverPTracks_[match][1] -> Fill(iMatchingSC->superCluster()->energy()/sqrt(refittedMom.mag2()));
1334  h_convSCdPhi_[match][1]->Fill( iMatchingSC->superCluster()->position().phi() - refittedMom.phi() );
1335  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(),aConv.zOfPrimaryVertexFromTracks());
1336  h_convSCdEta_[match][1]->Fill( iMatchingSC->superCluster()->position().eta() - ConvEta );
1337  }
1338  }
1339 
1340 
1341  if ( phoIsInEndcap ) {
1342  h_invMass_[match][2] ->Fill(invM);
1343  h_vtxChi2Prob_[match][2] ->Fill (chi2Prob);
1344  h_distMinAppTracks_[match][2] ->Fill (aConv.distOfMinimumApproach());
1345  h_DPhiTracksAtVtx_[match][2]->Fill( dPhiTracksAtVtx);
1346  h_DCotTracks_[match][2] ->Fill ( aConv.pairCotThetaSeparation() );
1347  h_lxybs_[match][2] ->Fill (lxy);
1348  h_maxNHitsBeforeVtx_[match][2] ->Fill (maxNHitsBeforeVtx);
1349  h_leadNHitsBeforeVtx_[match][2] ->Fill (leadNHitsBeforeVtx);
1350  h_trailNHitsBeforeVtx_[match][2] ->Fill (trailNHitsBeforeVtx);
1351  h_sumNHitsBeforeVtx_[match][2] ->Fill (sumNHitsBeforeVtx);
1352  h_deltaExpectedHitsInner_[match][2] ->Fill (deltaExpectedHitsInner);
1353  h_leadExpectedHitsInner_[match][2] ->Fill (leadExpectedHitsInner);
1354  h_maxDlClosestHitToVtx_[match][2] ->Fill (maxDlClosestHitToVtx);
1355  h_maxDlClosestHitToVtxSig_[match][2] ->Fill (maxDlClosestHitToVtxSig);
1356  h_nSharedHits_[match][2] ->Fill (aConv.nSharedHits());
1357  if ( matchConvSC ) {
1358  h_EoverPTracks_[match][2] ->Fill (iMatchingSC->superCluster()->energy()/sqrt(refittedMom.mag2()));
1359  h_convSCdPhi_[match][2]->Fill( iMatchingSC->superCluster()->position().phi() - refittedMom.phi() );
1360  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(),aConv.zOfPrimaryVertexFromTracks());
1361  h_convSCdEta_[match][2]->Fill( iMatchingSC->superCluster()->position().eta() - ConvEta );
1362 
1363  }
1364  }
1365 
1366  h_convVtxRvsZ_[0] ->Fill ( fabs (aConv.conversionVertex().position().z() ), sqrt(aConv.conversionVertex().position().perp2()) ) ;
1367  h_convVtxYvsX_ ->Fill ( aConv.conversionVertex().position().x(), aConv.conversionVertex().position().y() );
1368  h_convVtxYvsX_zoom_[0] ->Fill ( aConv.conversionVertex().position().x(), aConv.conversionVertex().position().y() );
1369  h_convVtxYvsX_zoom_[1] ->Fill ( aConv.conversionVertex().position().x(), aConv.conversionVertex().position().y() );
1370 
1371 
1372  // quantities per track: all conversions
1373  for (unsigned int i=0; i<tracks.size(); i++) {
1374  double d0;
1375  if (valid_pvtx){
1376  d0 = - tracks[i]->dxy(the_pvtx.position());
1377  } else {
1378  d0 = tracks[i]->d0();
1379  }
1380  h_TkD0_[match]->Fill ( d0* tracks[i]->charge() );
1381  h_nHitsBeforeVtx_[match]->Fill ( aConv.nHitsBeforeVtx().size()>1 ? aConv.nHitsBeforeVtx().at(i) : 0 );
1382  h_dlClosestHitToVtx_[match]->Fill ( aConv.dlClosestHitToVtx().size()>1 ? aConv.dlClosestHitToVtx().at(i).value() : 0 );
1383  h_dlClosestHitToVtxSig_[match]->Fill ( aConv.dlClosestHitToVtx().size()>1 ? aConv.dlClosestHitToVtx().at(i).value()/aConv.dlClosestHitToVtx().at(i).error() : 0 );
1384 
1385  nHitsVsEta_[match] ->Fill (mcEta_, float(tracks[i]->numberOfValidHits()) );
1386  nHitsVsR_[match] ->Fill (mcConvR_, float(tracks[i]->numberOfValidHits()) );
1387  p_nHitsVsEta_[match] ->Fill (mcEta_, float(tracks[i]->numberOfValidHits()) -0.0001);
1388  p_nHitsVsR_[match] ->Fill (mcConvR_, float(tracks[i]->numberOfValidHits()) -0.0001);
1389  h_tkChi2_[match] ->Fill (tracks[i]->normalizedChi2() );
1390  h_tkChi2Large_[match] ->Fill (tracks[i]->normalizedChi2() );
1391  h2_Chi2VsEta_[match] ->Fill( mcEta_, tracks[i]->normalizedChi2() );
1392  h2_Chi2VsR_[match] ->Fill( mcConvR_, tracks[i]->normalizedChi2() );
1393  p_Chi2VsEta_[match] ->Fill( mcEta_, tracks[i]->normalizedChi2() );
1394  p_Chi2VsR_[match] ->Fill( mcConvR_, tracks[i]->normalizedChi2() );
1395 
1396  }
1397 
1398  bool associated = false;
1399  float mcConvPt_= -99999999;
1400  // float mcPhi= 0; // unused
1401  float simPV_Z=0;
1402  for ( std::vector<PhotonMCTruth>::const_iterator mcPho=mcPhotons.begin(); mcPho !=mcPhotons.end(); mcPho++) {
1403  mcConvPt_= (*mcPho).fourMomentum().et();
1404  float mcPhi= (*mcPho).fourMomentum().phi();
1405  simPV_Z = (*mcPho).primaryVertex().z();
1406  mcPhi_= phiNormalization(mcPhi);
1407  mcEta_= (*mcPho).fourMomentum().pseudoRapidity();
1408  mcEta_ = etaTransformation(mcEta_, (*mcPho).primaryVertex().z() );
1409  mcConvR_= (*mcPho).vertex().perp();
1410  mcConvX_= (*mcPho).vertex().x();
1411  mcConvY_= (*mcPho).vertex().y();
1412  mcConvZ_= (*mcPho).vertex().z();
1413  mcConvEta_= (*mcPho).vertex().eta();
1414  mcConvPhi_= (*mcPho).vertex().phi();
1415  if ( fabs(mcEta_) > END_HI ) continue;
1416  if (mcConvPt_<minPhoPtForPurity) continue;
1417  if (fabs(mcEta_)>maxPhoEtaForPurity) continue;
1418  if (fabs(mcConvZ_)>maxPhoZForPurity) continue;
1419  if (mcConvR_>maxPhoRForEffic) continue;
1420 
1421  if ( (*mcPho).isAConversion() != 1 ) continue;
1422  if (!( ( fabs(mcEta_) <= BARL && mcConvR_ <85 ) ||
1423  ( fabs(mcEta_) > BARL && fabs(mcEta_) <=END_HI && fabs( (*mcPho).vertex().z() ) < 210 ) ) )
1424  continue;
1425 
1426 
1427  theConvTP_.clear();
1428  for(const TrackingParticleRef tp: tpForFakeRate) {
1429  if ( fabs( tp->vx() - (*mcPho).vertex().x() ) < 0.0001 &&
1430  fabs( tp->vy() - (*mcPho).vertex().y() ) < 0.0001 &&
1431  fabs( tp->vz() - (*mcPho).vertex().z() ) < 0.0001) {
1432  theConvTP_.push_back( tp );
1433  }
1434  }
1435 
1436  if ( theConvTP_.size() < 2 ) continue;
1437 
1438  //associated = false;
1439  reco::RecoToSimCollection p1 = theTrackAssociator->associateRecoToSim(tc1,theConvTP_);
1440  reco::RecoToSimCollection p2 = theTrackAssociator->associateRecoToSim(tc2,theConvTP_);
1441  try{
1442  std::vector<std::pair<TrackingParticleRef, double> > tp1 = p1[tk1];
1443  std::vector<std::pair<TrackingParticleRef, double> > tp2 = p2[tk2];
1444  if (!(tp1.size()&&tp2.size())){
1445  tp1 = p1[tk2];
1446  tp2 = p2[tk1];
1447  }
1448  if (tp1.size()&&tp2.size()) {
1449  TrackingParticleRef tpr1 = tp1.front().first;
1450  TrackingParticleRef tpr2 = tp2.front().first;
1451  if (abs(tpr1->pdgId())==11&&abs(tpr2->pdgId())==11&& tpr1->pdgId()*tpr2->pdgId()<0) {
1452  if ( (tpr1->parentVertex()->sourceTracks_end()-tpr1->parentVertex()->sourceTracks_begin()==1) &&
1453  (tpr2->parentVertex()->sourceTracks_end()-tpr2->parentVertex()->sourceTracks_begin()==1)) {
1454  if (tpr1->parentVertex().key()==tpr2->parentVertex().key() && ((*tpr1->parentVertex()->sourceTracks_begin())->pdgId()==22)) {
1455  mcConvR_ = sqrt(tpr1->parentVertex()->position().Perp2());
1456  mcConvZ_ = tpr1->parentVertex()->position().z();
1457  mcConvX_ = tpr1->parentVertex()->position().x();
1458  mcConvY_ = tpr1->parentVertex()->position().y();
1459  mcConvEta_ = tpr1->parentVertex()->position().eta();
1460  mcConvPhi_ = tpr1->parentVertex()->position().phi();
1461  mcConvPt_ = sqrt((*tpr1->parentVertex()->sourceTracks_begin())->momentum().Perp2());
1462  //std::cout << " Reco to Sim mcconvpt " << mcConvPt_ << std::endl;
1463  //cout << "associated track1 to " << tpr1->pdgId() << " with p=" << tpr1->p4() << " with pT=" << tpr1->pt() << endl;
1464  //cout << "associated track2 to " << tpr2->pdgId() << " with p=" << tpr2->p4() << " with pT=" << tpr2->pt() << endl;
1465  associated = true;
1466  break;
1467  }
1468  }
1469  }
1470  }
1471  } catch (Exception event) {
1472  //cout << "do not continue: " << event.what() << endl;
1473  //continue;
1474  }
1475 
1476  }// end loop on sim photons
1477 
1478 
1479  if (0) {
1480  theConvTP_.clear();
1481  theConvTP_ = tpForFakeRate;
1482  reco::RecoToSimCollection p1incl = theTrackAssociator->associateRecoToSim(tc1,theConvTP_);
1483  reco::RecoToSimCollection p2incl = theTrackAssociator->associateRecoToSim(tc2,theConvTP_);
1484 
1485 
1486  for ( std::vector<PhotonMCTruth>::const_iterator mcPho=mcPhotons.begin(); mcPho !=mcPhotons.end(); mcPho++) {
1487  mcConvPt_= (*mcPho).fourMomentum().et();
1488  float mcPhi= (*mcPho).fourMomentum().phi();
1489  simPV_Z = (*mcPho).primaryVertex().z();
1490  mcPhi_= phiNormalization(mcPhi);
1491  mcEta_= (*mcPho).fourMomentum().pseudoRapidity();
1492  mcEta_ = etaTransformation(mcEta_, (*mcPho).primaryVertex().z() );
1493  mcConvR_= (*mcPho).vertex().perp();
1494  mcConvX_= (*mcPho).vertex().x();
1495  mcConvY_= (*mcPho).vertex().y();
1496  mcConvZ_= (*mcPho).vertex().z();
1497  mcConvEta_= (*mcPho).vertex().eta();
1498  mcConvPhi_= (*mcPho).vertex().phi();
1499  if ( fabs(mcEta_) > END_HI ) continue;
1500  if (mcConvPt_<minPhoPtForPurity) continue;
1501  if (fabs(mcEta_)>maxPhoEtaForPurity) continue;
1502  if (fabs(mcConvZ_)>maxPhoZForPurity) continue;
1503  if (mcConvR_>maxPhoRForEffic) continue;
1504 
1505  if ( (*mcPho).isAConversion() != 1 ) continue;
1506  if (!( ( fabs(mcEta_) <= BARL && mcConvR_ <85 ) ||
1507  ( fabs(mcEta_) > BARL && fabs(mcEta_) <=END_HI && fabs( (*mcPho).vertex().z() ) < 210 ) ) )
1508  continue;
1509 
1510 
1511  theConvTP_.clear();
1512  for(TrackingParticleRef tp: tpForFakeRate) {
1513  if ( fabs( tp->vx() - (*mcPho).vertex().x() ) < 0.0001 &&
1514  fabs( tp->vy() - (*mcPho).vertex().y() ) < 0.0001 &&
1515  fabs( tp->vz() - (*mcPho).vertex().z() ) < 0.0001) {
1516  theConvTP_.push_back( tp );
1517 
1518 
1519  }
1520  }
1521 
1522  if ( theConvTP_.size() < 2 ) continue;
1523 
1524  //associated = false;
1525  reco::RecoToSimCollection p1 = theTrackAssociator->associateRecoToSim(tc1,theConvTP_);
1526  reco::RecoToSimCollection p2 = theTrackAssociator->associateRecoToSim(tc2,theConvTP_);
1527 
1528 
1529 
1530 
1531 
1532  if ( (p1incl.size() && p2incl.size()) && (p1.size() || p2.size()) ) { // associated = true;
1533  try{
1534  std::vector<std::pair<TrackingParticleRef, double> > tp1 = p1incl[tk1];
1535  std::vector<std::pair<TrackingParticleRef, double> > tp2 = p2incl[tk2];
1536  if (!(tp1.size()&&tp2.size())){
1537  tp1 = p1[tk2];
1538  tp2 = p2[tk1];
1539  }
1540  if (tp1.size()&&tp2.size()) {
1541  TrackingParticleRef tpr1 = tp1.front().first;
1542  TrackingParticleRef tpr2 = tp2.front().first;
1543  if (abs(tpr1->pdgId())==11&&abs(tpr2->pdgId())==11 && tpr1->pdgId()*tpr2->pdgId()<0) {
1544  if ( ((tpr1->parentVertex()->sourceTracks_end()-tpr1->parentVertex()->sourceTracks_begin()>=1) && (*tpr1->parentVertex()->sourceTracks_begin())->pdgId()==22) &&
1545  ((tpr2->parentVertex()->sourceTracks_end()-tpr2->parentVertex()->sourceTracks_begin()>=1) && (*tpr2->parentVertex()->sourceTracks_begin())->pdgId()==22) ) {
1546 
1547  // if ( fabs(tpr1->vx() - tpr2->vx()) < 0.1 && fabs(tpr1->vy() - tpr2->vy()) < 0.1 && fabs(tpr1->vz() - tpr2->vz()) < 0.1) {
1548  //if (((*tpr1->parentVertex()->sourceTracks_begin())->pdgId()==22) || ((*tpr2->parentVertex()->sourceTracks_begin())->pdgId()==22)) {
1549 // mcConvR_ = sqrt(tpr1->parentVertex()->position().Perp2());
1550 // mcConvZ_ = tpr1->parentVertex()->position().z();
1551 // mcConvX_ = tpr1->parentVertex()->position().x();
1552 // mcConvY_ = tpr1->parentVertex()->position().y();
1553 // mcConvEta_ = tpr1->parentVertex()->position().eta();
1554 // mcConvPhi_ = tpr1->parentVertex()->position().phi();
1555 // mcConvPt_ = sqrt((*tpr1->parentVertex()->sourceTracks_begin())->momentum().Perp2());
1556  //std::cout << " Reco to Sim mcconvpt " << mcConvPt_ << std::endl;
1557  //cout << "associated track1 to " << tpr1->pdgId() << " with p=" << tpr1->p4() << " with pT=" << tpr1->pt() << endl;
1558  //cout << "associated track2 to " << tpr2->pdgId() << " with p=" << tpr2->p4() << " with pT=" << tpr2->pt() << endl;
1559  associated = true;
1560  break;
1561  //}
1562  //}
1563  }
1564  }
1565  }
1566  } catch (Exception event) {
1567  //cout << "do not continue: " << event.what() << endl;
1568  //continue;
1569  }
1570 
1571  }
1572 
1573  }
1574  }
1575 
1576  if ( associated ) match=1;
1577  else
1578  match=2;
1579 
1580  h_match_->Fill(float(match));
1582  if ( match == 1) nRecConvAss_++;
1583  h_convEta_[match][0]->Fill( refittedMom.eta() );
1584  h_convEta_[match][1]->Fill( refittedMom.eta() );
1585  if (matchConvSC) h_convEtaMatchSC_[match][0]->Fill( refittedMom.eta() );
1586  h_convPhi_[match][0]->Fill( refittedMom.phi() );
1587  h_convR_[match][0]->Fill( sqrt(aConv.conversionVertex().position().perp2()) );
1588  h_convZ_[match][0]->Fill( aConv.conversionVertex().position().z() );
1589  h_convPt_[match][0]->Fill( sqrt(refittedMom.perp2()) );
1590  h_invMass_[match][0] ->Fill( invM);
1591  h_vtxChi2Prob_[match][0] ->Fill (chi2Prob);
1592  h_DPhiTracksAtVtx_[match][0]->Fill( dPhiTracksAtVtx);
1593  h_DCotTracks_[match][0] ->Fill ( aConv.pairCotThetaSeparation() );
1594  h_distMinAppTracks_[match][0] ->Fill (aConv.distOfMinimumApproach());
1595  h_lxybs_[match][0] ->Fill (lxy);
1596  h_maxNHitsBeforeVtx_[match][0] ->Fill (maxNHitsBeforeVtx);
1597  h_leadNHitsBeforeVtx_[match][0] ->Fill (leadNHitsBeforeVtx);
1598  h_trailNHitsBeforeVtx_[match][0] ->Fill (trailNHitsBeforeVtx);
1599  h_sumNHitsBeforeVtx_[match][0] ->Fill (sumNHitsBeforeVtx);
1600  h_deltaExpectedHitsInner_[match][0] ->Fill (deltaExpectedHitsInner);
1601  h_leadExpectedHitsInner_[match][0] ->Fill (leadExpectedHitsInner);
1602  h_maxDlClosestHitToVtx_[match][0] ->Fill (maxDlClosestHitToVtx);
1603  h_maxDlClosestHitToVtxSig_[match][0] ->Fill (maxDlClosestHitToVtxSig);
1604  h_nSharedHits_[match][0] ->Fill (aConv.nSharedHits());
1605  if ( matchConvSC ) {
1606  //h_EoverPTracks_[match][0] ->Fill (aConv.EoverPrefittedTracks());
1607  h_EoverPTracks_[match][0] ->Fill (iMatchingSC->superCluster()->energy()/sqrt(refittedMom.mag2()));
1608  h_convSCdPhi_[match][0]->Fill( iMatchingSC->superCluster()->position().phi() - refittedMom.phi() );
1609  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(),aConv.zOfPrimaryVertexFromTracks());
1610  h_convSCdEta_[match][0]->Fill( iMatchingSC->superCluster()->position().eta() - ConvEta );
1611 
1612  }
1613  if ( match==1) {
1614  h2_photonPtRecVsPtSim_->Fill ( mcConvPt_, sqrt(refittedMom.perp2()) );
1615  h_convPtRes_[0]->Fill ( sqrt(refittedMom.perp2())/mcConvPt_);
1616  }
1617 
1618  if ( phoIsInBarrel ) {
1619  h_invMass_[match][1] ->Fill(invM);
1620  h_vtxChi2Prob_[match][1] ->Fill (chi2Prob);
1621  h_DPhiTracksAtVtx_[match][1]->Fill( dPhiTracksAtVtx);
1622  h_DCotTracks_[match][1] ->Fill ( aConv.pairCotThetaSeparation() );
1623  h_distMinAppTracks_[match][1] ->Fill (aConv.distOfMinimumApproach());
1624  h_lxybs_[match][1] ->Fill (lxy);
1625  h_maxNHitsBeforeVtx_[match][1] ->Fill (maxNHitsBeforeVtx);
1626  h_leadNHitsBeforeVtx_[match][1] ->Fill (leadNHitsBeforeVtx);
1627  h_trailNHitsBeforeVtx_[match][1] ->Fill (trailNHitsBeforeVtx);
1628  h_sumNHitsBeforeVtx_[match][1] ->Fill (sumNHitsBeforeVtx);
1629  h_deltaExpectedHitsInner_[match][1] ->Fill (deltaExpectedHitsInner);
1630  h_leadExpectedHitsInner_[match][1] ->Fill (leadExpectedHitsInner);
1631  h_maxDlClosestHitToVtx_[match][1] ->Fill (maxDlClosestHitToVtx);
1632  h_maxDlClosestHitToVtxSig_[match][1] ->Fill (maxDlClosestHitToVtxSig);
1633  h_nSharedHits_[match][1] ->Fill (aConv.nSharedHits());
1634  if ( matchConvSC ) {
1635  // h_EoverPTracks_[match][1] ->Fill (aConv.EoverPrefittedTracks());
1636  h_EoverPTracks_[match][1] ->Fill (iMatchingSC->superCluster()->energy()/sqrt(refittedMom.mag2()));
1637  h_convSCdPhi_[match][1]->Fill( iMatchingSC->superCluster()->position().phi() - refittedMom.phi() );
1638  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(),aConv.zOfPrimaryVertexFromTracks());
1639  h_convSCdEta_[match][1]->Fill( iMatchingSC->superCluster()->position().eta() - ConvEta );
1640 
1641  }
1642  if ( match==1) h_convPtRes_[1]->Fill ( sqrt(refittedMom.perp2())/mcConvPt_);
1643  }
1644 
1645 
1646  if ( phoIsInEndcap ) {
1647  h_invMass_[match][2] ->Fill(invM);
1648  h_vtxChi2Prob_[match][2] ->Fill (chi2Prob);
1649  h_DPhiTracksAtVtx_[match][2]->Fill( dPhiTracksAtVtx);
1650  h_DCotTracks_[match][2] ->Fill ( aConv.pairCotThetaSeparation() );
1651  h_distMinAppTracks_[match][2] ->Fill (aConv.distOfMinimumApproach());
1652  h_lxybs_[match][2] ->Fill (lxy);
1653  h_maxNHitsBeforeVtx_[match][2] ->Fill (maxNHitsBeforeVtx);
1654  h_leadNHitsBeforeVtx_[match][2] ->Fill (leadNHitsBeforeVtx);
1655  h_trailNHitsBeforeVtx_[match][2] ->Fill (trailNHitsBeforeVtx);
1656  h_sumNHitsBeforeVtx_[match][2] ->Fill (sumNHitsBeforeVtx);
1657  h_deltaExpectedHitsInner_[match][2] ->Fill (deltaExpectedHitsInner);
1658  h_leadExpectedHitsInner_[match][2] ->Fill (leadExpectedHitsInner);
1659  h_maxDlClosestHitToVtx_[match][2] ->Fill (maxDlClosestHitToVtx);
1660  h_maxDlClosestHitToVtxSig_[match][2] ->Fill (maxDlClosestHitToVtxSig);
1661  h_nSharedHits_[match][2] ->Fill (aConv.nSharedHits());
1662  if ( matchConvSC ) {
1663  // h_EoverPTracks_[match][2] ->Fill (aConv.EoverPrefittedTracks());
1664  h_EoverPTracks_[match][2] ->Fill (iMatchingSC->superCluster()->energy()/sqrt(refittedMom.mag2()));
1665  h_convSCdPhi_[match][2]->Fill( iMatchingSC->superCluster()->position().phi() - refittedMom.phi() );
1666  double ConvEta = etaTransformation(aConv.refittedPairMomentum().eta(),aConv.zOfPrimaryVertexFromTracks());
1667  h_convSCdEta_[match][2]->Fill( iMatchingSC->superCluster()->position().eta() - ConvEta );
1668  }
1669  if ( match==1) h_convPtRes_[2]->Fill ( sqrt(refittedMom.perp2())/mcConvPt_);
1670  }
1671 
1672 
1673  if ( match == 1 ) {
1674  h_convVtxdX_ ->Fill ( aConv.conversionVertex().position().x() - mcConvX_);
1675  h_convVtxdY_ ->Fill ( aConv.conversionVertex().position().y() - mcConvY_);
1676  h_convVtxdZ_ ->Fill ( aConv.conversionVertex().position().z() - mcConvZ_);
1677  h_convVtxdR_ ->Fill ( sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
1678  h_convVtxdPhi_ ->Fill ( aConv.conversionVertex().position().phi() - mcConvPhi_);
1679  h_convVtxdEta_ ->Fill ( aConv.conversionVertex().position().eta() - mcConvEta_);
1680  h2_convVtxdRVsR_ ->Fill (mcConvR_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_ );
1681  h2_convVtxdRVsEta_ ->Fill (mcEta_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_ );
1682  p_convVtxdRVsR_ ->Fill (mcConvR_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_ );
1683  p_convVtxdRVsEta_ ->Fill (mcEta_, sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_ );
1684  p_convVtxdXVsX_ ->Fill (mcConvX_, aConv.conversionVertex().position().x() - mcConvX_ );
1685  p_convVtxdYVsY_ ->Fill (mcConvY_, aConv.conversionVertex().position().y() - mcConvY_ );
1686  p_convVtxdZVsZ_ ->Fill (mcConvZ_, aConv.conversionVertex().position().z() - mcConvZ_ );
1687  p_convVtxdZVsR_ ->Fill (mcConvR_, aConv.conversionVertex().position().z() - mcConvZ_ );
1688 
1689  float dR=sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_;
1690  float dZ=aConv.conversionVertex().position().z() - mcConvZ_;
1691  p2_convVtxdRVsRZ_ ->Fill (mcConvZ_,mcConvR_, dR );
1692  p2_convVtxdZVsRZ_ ->Fill (mcConvZ_,mcConvR_, dZ );
1693 
1694 
1695 
1696 
1697  h2_convVtxRrecVsTrue_ -> Fill (mcConvR_, sqrt(aConv.conversionVertex().position().perp2()) );
1698 
1699 
1700  h_zPVFromTracks_[match]->Fill ( aConv.zOfPrimaryVertexFromTracks() );
1701  h_dzPVFromTracks_[match]->Fill ( aConv.zOfPrimaryVertexFromTracks() - simPV_Z );
1702  h2_dzPVVsR_ ->Fill(mcConvR_, aConv.zOfPrimaryVertexFromTracks() - simPV_Z );
1703  p_dzPVVsR_ ->Fill(mcConvR_, aConv.zOfPrimaryVertexFromTracks() - simPV_Z );
1704 
1705  if ( phoIsInBarrel ) {
1706  h_convVtxdX_barrel_ ->Fill ( aConv.conversionVertex().position().x() - mcConvX_);
1707  h_convVtxdY_barrel_ ->Fill ( aConv.conversionVertex().position().y() - mcConvY_);
1708  h_convVtxdZ_barrel_ ->Fill ( aConv.conversionVertex().position().z() - mcConvZ_);
1709  h_convVtxdR_barrel_ ->Fill ( sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
1710 
1711  }
1712  if ( phoIsInEndcap ) {
1713  h_convVtxdX_endcap_ ->Fill ( aConv.conversionVertex().position().x() - mcConvX_);
1714  h_convVtxdY_endcap_ ->Fill ( aConv.conversionVertex().position().y() - mcConvY_);
1715  h_convVtxdZ_endcap_ ->Fill ( aConv.conversionVertex().position().z() - mcConvZ_);
1716  h_convVtxdR_endcap_ ->Fill ( sqrt(aConv.conversionVertex().position().perp2()) - mcConvR_);
1717 
1718  }
1719 
1720 
1721  }
1722 
1724  for (unsigned int i=0; i<tracks.size(); i++) {
1725  //std::cout << " Loop over tracks pt " << tracks[i]->pt() << std::endl;
1726  RefToBase<reco::Track> tfrb(aConv.tracks()[i] );
1727  itAss= myAss.find( tfrb.get() );
1728 
1729  nHitsVsEta_[match] ->Fill (mcEta_, float(tracks[i]->numberOfValidHits()) );
1730  nHitsVsR_[match] ->Fill (mcConvR_, float(tracks[i]->numberOfValidHits()) );
1731  p_nHitsVsEta_[match] ->Fill (mcEta_, float(tracks[i]->numberOfValidHits()) -0.0001);
1732  p_nHitsVsR_[match] ->Fill (mcConvR_, float(tracks[i]->numberOfValidHits()) -0.0001);
1733  h_tkChi2_[match] ->Fill (tracks[i]->normalizedChi2() );
1734  h_tkChi2Large_[match] ->Fill (tracks[i]->normalizedChi2() );
1735  h2_Chi2VsEta_[match] ->Fill( mcEta_, tracks[i]->normalizedChi2() );
1736  h2_Chi2VsR_[match] ->Fill( mcConvR_, tracks[i]->normalizedChi2() );
1737  p_Chi2VsEta_[match] ->Fill( mcEta_, tracks[i]->normalizedChi2() );
1738  p_Chi2VsR_[match] ->Fill( mcConvR_, tracks[i]->normalizedChi2() );
1739  double d0;
1740  if (valid_pvtx){
1741  d0 = - tracks[i]->dxy(the_pvtx.position());
1742  } else {
1743  d0 = tracks[i]->d0();
1744  }
1745  h_TkD0_[match]->Fill (d0* tracks[i]->charge() );
1746  h_nHitsBeforeVtx_[match]->Fill ( aConv.nHitsBeforeVtx().size()>1 ? aConv.nHitsBeforeVtx().at(i) : 0 );
1747  h_dlClosestHitToVtx_[match]->Fill ( aConv.dlClosestHitToVtx().size()>1 ? aConv.dlClosestHitToVtx().at(i).value() : 0 );
1748  h_dlClosestHitToVtxSig_[match]->Fill ( aConv.dlClosestHitToVtx().size()>1 ? aConv.dlClosestHitToVtx().at(i).value()/aConv.dlClosestHitToVtx().at(i).error() : 0 );
1749 
1750 
1751  if ( itAss == myAss.end() ) continue;
1752  reco::Track refTrack= aConv.conversionVertex().refittedTracks()[i];
1753 
1754  float simPt = sqrt( ((*itAss).second)->momentum().perp2() );
1755  float recPt = refTrack.pt();
1756  float ptres= recPt - simPt ;
1757  //float pterror = aConv.tracks()[i]->ptError();
1758  float pterror = aConv.conversionVertex().refittedTracks()[i].ptError();
1759  h2_PtRecVsPtSim_[0]->Fill ( simPt, recPt);
1760  h_TkPtPull_[0] ->Fill(ptres/pterror);
1761  h2_TkPtPull_[0] ->Fill(mcEta_, ptres/pterror);
1762 
1763  if ( phoIsInBarrel ) {
1764  h_TkPtPull_[1] ->Fill(ptres/pterror);
1765  h2_PtRecVsPtSim_[1]->Fill ( simPt, recPt);
1766  }
1767  if ( phoIsInEndcap ) {
1768  h_TkPtPull_[2] ->Fill(ptres/pterror);
1769  h2_PtRecVsPtSim_[2]->Fill ( simPt, recPt);
1770  }
1771  } // end loop over track
1772 
1773 
1774 
1775  } // loop over reco conversions
1776 
1777 
1778  h_nConv_[0][0]->Fill (float(nRecConv_));
1779  h_nConv_[1][0]->Fill (float(nRecConvAss_));
1780 
1781 
1782 
1783 }
1784 
1785 
1786 
1787 
1788 
1789 void TkConvValidator::endJob() {
1790 
1791 
1792  std::string outputFileName = parameters_.getParameter<std::string>("OutputFileName");
1793  if ( ! isRunCentrally_ ) {
1794  dbe_->save(outputFileName);
1795  }
1796 
1797  edm::LogInfo("TkConvValidator") << "Analyzed " << nEvt_ << "\n";
1798  // std::cout << "::endJob Analyzed " << nEvt_ << " events " << " with total " << nPho_ << " Photons " << "\n";
1799  // std::cout << "TkConvValidator::endJob Analyzed " << nEvt_ << " events " << "\n";
1800 
1801  return ;
1802 }
1803 
1804 
1805 math::XYZVector TkConvValidator::recalculateMomentumAtFittedVertex ( const MagneticField& mf, const TrackerGeometry& trackerGeom, const edm::RefToBase<reco::Track>& tk, const reco::Vertex& vtx) {
1806 
1807  math::XYZVector result;
1808  Surface::RotationType rot;
1809  auto scp = new SimpleCylinderBounds( sqrt(vtx.position().perp2())-0.001f,
1810  sqrt(vtx.position().perp2())+0.001f,
1811  -fabs(vtx.position().z()),
1812  fabs(vtx.position().z())
1813  );
1814  ReferenceCountingPointer<Cylinder> theBarrel_(new Cylinder(Cylinder::computeRadius(*scp), Surface::PositionType(0,0,0), rot,scp));
1815 
1816  ReferenceCountingPointer<Disk> theDisk_(new Disk( Surface::PositionType( 0, 0, vtx.position().z()), rot,
1817  new SimpleDiskBounds( 0, sqrt(vtx.position().perp2()), -0.001, 0.001) )
1818  );
1819 
1820 
1821  const TrajectoryStateOnSurface myTSOS = trajectoryStateTransform::innerStateOnSurface(*tk, trackerGeom, &mf);
1822  PropagatorWithMaterial propag( anyDirection, 0.000511, &mf );
1823  TrajectoryStateOnSurface stateAtVtx;
1824  stateAtVtx = propag.propagate(myTSOS, *theBarrel_);
1825  if (!stateAtVtx.isValid() ) {
1826  stateAtVtx = propag.propagate(myTSOS, *theDisk_);
1827  }
1828  if (stateAtVtx.isValid()){
1829  return math::XYZVector ( double(stateAtVtx.globalMomentum().x()), double(stateAtVtx.globalMomentum().y()), double(stateAtVtx.globalMomentum().z()));
1830  } else {
1831  return math::XYZVector(0.,0.,0.);
1832  }
1833 
1834 
1835 
1836 }
1837 
1838 
1839 float TkConvValidator::phiNormalization(float & phi)
1840 {
1841  //---Definitions
1842  const float PI = 3.1415927;
1843  const float TWOPI = 2.0*PI;
1844 
1845 
1846  if(phi > PI) {phi = phi - TWOPI;}
1847  if(phi < -PI) {phi = phi + TWOPI;}
1848 
1849  // cout << " Float_t PHInormalization out " << PHI << endl;
1850  return phi;
1851 
1852 }
1853 
1854 
1855 float TkConvValidator::etaTransformation( float EtaParticle , float Zvertex) {
1856 
1857  //---Definitions
1858  const float PI = 3.1415927;
1859 
1860  //---Definitions for ECAL
1861  const float R_ECAL = 136.5;
1862  const float Z_Endcap = 328.0;
1863  const float etaBarrelEndcap = 1.479;
1864 
1865  //---ETA correction
1866 
1867  float Theta = 0.0 ;
1868  float ZEcal = R_ECAL*sinh(EtaParticle)+Zvertex;
1869 
1870  if(ZEcal != 0.0) Theta = atan(R_ECAL/ZEcal);
1871  if(Theta<0.0) Theta = Theta+PI ;
1872  float ETA = - log(tan(0.5*Theta));
1873 
1874  if( fabs(ETA) > etaBarrelEndcap )
1875  {
1876  float Zend = Z_Endcap ;
1877  if(EtaParticle<0.0 ) Zend = -Zend ;
1878  float Zlen = Zend - Zvertex ;
1879  float RR = Zlen/sinh(EtaParticle);
1880  Theta = atan(RR/Zend);
1881  if(Theta<0.0) Theta = Theta+PI ;
1882  ETA = - log(tan(0.5*Theta));
1883  }
1884  //---Return the result
1885  return ETA;
1886  //---end
1887 }
1888 
1889 
reco::Conversion::conversionVertex
const reco::Vertex & conversionVertex() const
returns the reco conversion vertex
Definition: Conversion.h:97
edm::RefToBase::get
value_type const * get() const
Definition: RefToBase.h:234
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
photonValidator_cfi.dPhiTracksBin
dPhiTracksBin
Definition: photonValidator_cfi.py:69
photonValidator_cfi.dPhiTracksMin
dPhiTracksMin
Definition: photonValidator_cfi.py:70
btagDijet_cfi.etMin
etMin
Definition: btagDijet_cfi.py:5
mps_fire.i
i
Definition: mps_fire.py:156
conversionPostprocessing_cfi.rBin
rBin
Definition: conversionPostprocessing_cfi.py:27
tkConvValidator_cfi.maxPhoZForEffic
maxPhoZForEffic
Definition: tkConvValidator_cfi.py:60
photonValidationSequence_cff.eoverpMax
eoverpMax
Definition: photonValidationSequence_cff.py:17
tkConvValidator_cfi.maxPhoEtaForEffic
maxPhoEtaForEffic
Definition: tkConvValidator_cfi.py:59
reco::Photon::superCluster
reco::SuperClusterRef superCluster() const
Ref to SuperCluster.
photonValidator_cfi.rMaxForXray
rMaxForXray
Definition: photonValidator_cfi.py:118
mps_fire.result
result
Definition: mps_fire.py:113
edm::AssociationMap::end
const_iterator end() const
last iterator over the map (read only)
Definition: AssociationMap.h:160
DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
CSCSegmentAlgorithmDF_cfi.chi2Max
chi2Max
Definition: CSCSegmentAlgorithmDF_cfi.py:11
conv
static HepMC::IO_HEPEVT conv
Definition: BeamHaloProducer.cc:50
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
photonValidator_cfi.resMax
resMax
Definition: photonValidator_cfi.py:91
hpstanc_transforms.max
max
Definition: hpstanc_transforms.py:5
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
reco::Conversion::quality
bool quality(ConversionQuality q) const
Definition: Conversion.h:181
reco::Vertex::y
double y() const
y coordinate
Definition: Vertex.h:113
TkConvValidator::recalculateMomentumAtFittedVertex
math::XYZVector recalculateMomentumAtFittedVertex(const MagneticField &mf, const TrackerGeometry &trackerGeom, const edm::RefToBase< reco::Track > &tk, const reco::Vertex &vtx)
Definition: TkConvValidator.cc:1805
reco::Vertex::isValid
bool isValid() const
Tells whether the vertex is valid.
Definition: Vertex.h:68
photonValidator_cfi.zMaxForXray
zMaxForXray
Definition: photonValidator_cfi.py:122
reco::Conversion::dlClosestHitToVtx
const std::vector< Measurement1DFloat > & dlClosestHitToVtx() const
Vector of signed decay length with uncertainty from nearest hit on track to the conversion vtx positi...
Definition: Conversion.h:163
edm::AssociationMap::find
const_iterator find(const key_type &k) const
find element with specified reference key
Definition: AssociationMap.h:162
tkConvValidator_cfi.maxPhoRForEffic
maxPhoRForEffic
Definition: tkConvValidator_cfi.py:61
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
reco::Conversion::zOfPrimaryVertexFromTracks
double zOfPrimaryVertexFromTracks(const math::XYZPoint &myBeamSpot=math::XYZPoint()) const
Definition: Conversion.h:145
reco::Conversion::distOfMinimumApproach
double distOfMinimumApproach() const
Definition: Conversion.h:125
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
reco::Conversion::pairCotThetaSeparation
double pairCotThetaSeparation() const
Delta cot(Theta) where Theta is the angle in the (y,z) plane between the two tracks. Original tracks are used.
Definition: Conversion.cc:209
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
photonValidator_cfi.resMin
resMin
Definition: photonValidator_cfi.py:90
edm::Ref::key
key_type key() const
Accessor for product key.
Definition: Ref.h:264
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
reco::Vertex::refittedTracks
const std::vector< Track > & refittedTracks() const
Returns the container of refitted tracks.
Definition: Vertex.h:164
correctedPhotons_cfi.photonCollection
photonCollection
Definition: correctedPhotons_cfi.py:12
ALCARECOTkAlBeamHalo_cff.etaMin
etaMin
GeV.
Definition: ALCARECOTkAlBeamHalo_cff.py:31
reco::Conversion::pairInvariantMass
double pairInvariantMass() const
if nTracks=2 returns the pair invariant mass. Original tracks are used here
Definition: Conversion.cc:190
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:109
reco::Conversion::refittedPairMomentum
math::XYZVectorF refittedPairMomentum() const
Conversion tracks momentum from the tracks refitted with vertex constraint.
Definition: Conversion.cc:248
q2
double q2[4]
Definition: TauolaWrapper.h:88
reco::ConversionCollection
std::vector< Conversion > ConversionCollection
collectin of Conversion objects
Definition: ConversionFwd.h:9
photonValidator_cfi.rBinForXray
rBinForXray
Definition: photonValidator_cfi.py:116
TkConvValidator::~TkConvValidator
virtual ~TkConvValidator()
Definition: TkConvValidator.cc:158
conversionPostprocessing_cfi.zMin
zMin
Definition: conversionPostprocessing_cfi.py:32
reco::TrackBase::algo
TrackAlgorithm algo() const
Definition: TrackBase.h:492
Disk
Definition: BoundDisk.h:19
TkConvValidator::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: TkConvValidator.cc:163
photonValidator_cfi.eoverpBin
eoverpBin
Definition: photonValidator_cfi.py:105
TkConvValidator::endRun
virtual void endRun(edm::Run const &r, edm::EventSetup const &es) override
Definition: TkConvValidator.cc:795
photonValidator_cfi.zMinForXray
zMinForXray
Definition: photonValidator_cfi.py:121
ETA
#define ETA
Definition: GenericBenchmark.cc:32
AlignmentTrackSelector_cfi.phiMin
phiMin
Definition: AlignmentTrackSelector_cfi.py:18
photonValidator_cfi.rMinForXray
rMinForXray
Definition: photonValidator_cfi.py:117
TkConvValidator::etaTransformation
float etaTransformation(float a, float b)
Definition: TkConvValidator.cc:1855
TkConvValidator::phiNormalization
float phiNormalization(float &a)
Definition: TkConvValidator.cc:1839
tkConvValidator_cfi.maxPhoEtaForPurity
maxPhoEtaForPurity
Definition: tkConvValidator_cfi.py:63
DQMStore::IBooker::bookProfile2D
MonitorElement * bookProfile2D(Args &&...args)
Definition: DQMStore.h:163
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
photonValidator_cfi.zBinForXray
zBinForXray
Definition: photonValidator_cfi.py:119
conversionPostprocessing_cfi.zMax
zMax
Definition: conversionPostprocessing_cfi.py:33
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
conversionPostprocessing_cfi.etaBin
etaBin
Definition: conversionPostprocessing_cfi.py:18
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:616
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
photonValidator_cfi.dCotTracksMax
dCotTracksMax
Definition: photonValidator_cfi.py:99
TkConvValidator::TkConvValidator
TkConvValidator(const edm::ParameterSet &)
Definition: TkConvValidator.cc:104
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
photonValidator_cfi.dCotTracksMin
dCotTracksMin
Definition: photonValidator_cfi.py:98
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
postProcessorL1Gen_cff.outputFileName
outputFileName
Definition: postProcessorL1Gen_cff.py:25
reco::Vertex::chi2
double chi2() const
chi-squares
Definition: Vertex.h:98
TWOPI
#define TWOPI
EgammaCoreTools.
Definition: DQMSourcePi0.cc:40
conversionPostprocessing_cfi.zBin
zBin
Definition: conversionPostprocessing_cfi.py:31
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
ChiSquaredProbability
float ChiSquaredProbability(double chiSquared, double nrDOF)
Definition: ChiSquaredProbability.cc:13
conversionPostprocessing_cfi.rMax
rMax
Definition: conversionPostprocessing_cfi.py:29
math::XYZVectorF
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< float > > XYZVectorF
spatial vector with cartesian internal representation
Definition: Vector3D.h:17
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:74
ALCARECOTkAlBeamHalo_cff.etaMax
etaMax
Definition: ALCARECOTkAlBeamHalo_cff.py:32
p2
double p2[4]
Definition: TauolaWrapper.h:90
PI
#define PI
Definition: QcdUeDQM.h:36
dbe_
DQMStore * dbe_
Definition: PFJetBenchmarkAnalyzer.cc:79
reco::Vertex::ndof
double ndof() const
Definition: Vertex.h:105
PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi.dR
dR
Definition: PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi.py:19
reco::deltaPhi
double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:22
reco::Vertex::x
double x() const
x coordinate
Definition: Vertex.h:111
photonValidator_cfi.dCotTracksBin
dCotTracksBin
Definition: photonValidator_cfi.py:97
tkConvValidator_cfi.maxPhoRForPurity
maxPhoRForPurity
Definition: tkConvValidator_cfi.py:65
cosmictrackingParticleSelector_cfi.pdgId
pdgId
Definition: cosmictrackingParticleSelector_cfi.py:6
reco::Conversion::nSharedHits
uint8_t nSharedHits() const
number of shared hits btw the two track
Definition: Conversion.h:165
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
edm::AssociationMap::size
size_type size() const
map size
Definition: AssociationMap.h:147
edm::Handle::product
T const * product() const
Definition: Handle.h:81
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
reco::TrackBase::hitPattern
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:445
reco::TrackToTrackingParticleAssociator::associateRecoToSim
reco::RecoToSimCollection associateRecoToSim(const edm::Handle< edm::View< reco::Track > > &tCH, const edm::Handle< TrackingParticleCollection > &tPCH) const
compare reco to sim the handle of reco::Track and TrackingParticle collections
Definition: TrackToTrackingParticleAssociator.h:63
q1
double q1[4]
Definition: TauolaWrapper.h:87
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:30
reco::PhotonCollection
std::vector< Photon > PhotonCollection
collectin of Photon objects
Definition: PhotonFwd.h:9
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
mixOne_simraw_on_sim_cfi.tracker
tracker
Definition: mixOne_simraw_on_sim_cfi.py:264
photonValidationSequence_cff.eoverpMin
eoverpMin
Definition: photonValidationSequence_cff.py:16
photonValidator_cfi.zBin2ForXray
zBin2ForXray
Definition: photonValidator_cfi.py:120
TkConvValidator::endJob
virtual void endJob()
Definition: TkConvValidator.cc:1789
conversionPostprocessing_cfi.etBin
etBin
Definition: conversionPostprocessing_cfi.py:14
Propagator::propagate
TrajectoryStateOnSurface propagate(STA const &state, SUR const &surface) const
Definition: Propagator.h:53
etaBarrelEndcap
static const float etaBarrelEndcap
Definition: ECALPositionCalculator.cc:12
photonValidator_cfi.resBin
resBin
Definition: photonValidator_cfi.py:89
Z_Endcap
static const float Z_Endcap
Definition: ECALPositionCalculator.cc:11
reco::return
return(e1-e2)*(e1-e2)+dp *dp
particleFlowSuperClusterECAL_cfi.vertexCollection
vertexCollection
Definition: particleFlowSuperClusterECAL_cfi.py:41
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:58
edm
HLT enums.
Definition: AlignableModifier.h:17
p1
double p1[4]
Definition: TauolaWrapper.h:89
TrajectoryStateOnSurface::globalMomentum
GlobalVector globalMomentum() const
Definition: TrajectoryStateOnSurface.h:86
Cylinder::computeRadius
static float computeRadius(Bounds const &bounds)
Definition: Cylinder.h:30
conversionPostprocessing_cfi.etaBin2
etaBin2
Definition: conversionPostprocessing_cfi.py:19
allConversions_cfi.d0
d0
Definition: allConversions_cfi.py:30
tkConvValidator_cfi.minPhoPtForEffic
minPhoPtForEffic
Definition: tkConvValidator_cfi.py:58
edm::RefToBaseVector::push_back
void push_back(const RefToBase< T > &)
Definition: RefToBaseVector.h:264
reco::Conversion::nHitsBeforeVtx
const std::vector< uint8_t > & nHitsBeforeVtx() const
Vector of the number of hits before the vertex along each track trajector.
Definition: Conversion.h:161
reco::BeamSpot::y0
double y0() const
y coordinate
Definition: BeamSpot.h:66
R_ECAL
static const float R_ECAL
Definition: ECALPositionCalculator.cc:10
photonValidator_cfi.dPhiTracksMax
dPhiTracksMax
Definition: photonValidator_cfi.py:71
conversionPostprocessing_cfi.rMin
rMin
Definition: conversionPostprocessing_cfi.py:28
TkConvValidator::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: TkConvValidator.cc:803
recPt
double recPt
Definition: PFJetBenchmarkAnalyzer.cc:77
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10
photonValidationSequence_cff.useTP
useTP
Definition: photonValidationSequence_cff.py:14
badGlobalMuonTaggersAOD_cff.vtx
vtx
Definition: badGlobalMuonTaggersAOD_cff.py:5
reco::TrackToTrackingParticleAssociator::associateSimToReco
reco::SimToRecoCollection associateSimToReco(const edm::Handle< edm::View< reco::Track > > &tCH, const edm::Handle< TrackingParticleCollection > &tPCH) const
compare reco to sim the handle of reco::Track and TrackingParticle collections
Definition: TrackToTrackingParticleAssociator.h:69
tkConvValidator_cfi.maxPhoZForPurity
maxPhoZForPurity
Definition: tkConvValidator_cfi.py:64
AlignmentTrackSelector_cfi.phiMax
phiMax
Definition: AlignmentTrackSelector_cfi.py:17
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
conversionPostprocessing_cfi.phiBin
phiBin
Definition: conversionPostprocessing_cfi.py:23
conversionPostprocessing_cfi.etMax
etMax
Definition: conversionPostprocessing_cfi.py:16
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
tkConvValidator_cfi.minPhoPtForPurity
minPhoPtForPurity
Definition: tkConvValidator_cfi.py:62
muonBadTrackFilter_cfi.chi2Min
chi2Min
Definition: muonBadTrackFilter_cfi.py:8
reco::Conversion::dPhiTracksAtVtx
double dPhiTracksAtVtx() const
Definition: Conversion.cc:313
reco::Conversion::tracks
std::vector< edm::RefToBase< reco::Track > > const & tracks() const
vector of track to base references
Definition: Conversion.cc:176
event
Definition: event.py:1
edm::Run
Definition: Run.h:42
reco::HitPattern::numberOfHits
int numberOfHits(HitCategory category) const
Definition: HitPattern.h:807
trajectoryStateTransform::innerStateOnSurface
TrajectoryStateOnSurface innerStateOnSurface(const reco::Track &tk, const TrackingGeometry &geom, const MagneticField *field, bool withErr=true)
Definition: TrajectoryStateTransform.cc:110
reco::BeamSpot::x0
double x0() const
x coordinate
Definition: BeamSpot.h:64
TkConvValidator::dqmBeginRun
virtual void dqmBeginRun(edm::Run const &r, edm::EventSetup const &theEventSetup) override
Definition: TkConvValidator.cc:785
Generated for CMSSW Reference Manual by doxygen 1.8.11