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