CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
V0Validator.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: V0Validator
4 // Class: V0Validator
5 //
13 //
14 // Original Author: Brian Drell
15 // Created: Wed Feb 18 17:21:04 MST 2009
16 //
17 //
18 
19 
20 #include "Validation/RecoVertex/interface/V0Validator.h"
21 #include "DataFormats/VertexReco/interface/VertexFwd.h"
22 
23 typedef std::vector<TrackingVertex> TrackingVertexCollection;
24 typedef edm::Ref<TrackingVertexCollection> TrackingVertexRef;
25 typedef edm::RefVector<edm::HepMCProduct, HepMC::GenVertex > GenVertexRefVector;
26 typedef edm::RefVector<edm::HepMCProduct, HepMC::GenParticle > GenParticleRefVector;
27 
28 V0Validator::V0Validator(const edm::ParameterSet& iConfig)
29  : theDQMRootFileName(iConfig.getUntrackedParameter<std::string>("DQMRootFileName"))
30  , dirName(iConfig.getUntrackedParameter<std::string>("dirName"))
31  , recoRecoToSimCollectionToken_( consumes<reco::RecoToSimCollection>( iConfig.getUntrackedParameter<edm::InputTag>("trackAssociatorMap") ) )
32  , recoSimToRecoCollectionToken_( consumes<reco::SimToRecoCollection>( iConfig.getUntrackedParameter<edm::InputTag>("trackAssociatorMap") ) )
33  , trackingParticleCollection_Eff_Token_( consumes<TrackingParticleCollection>( iConfig.getUntrackedParameter<edm::InputTag>("trackingParticleCollectionEff") ) )
34  , vec_recoVertex_Token_( consumes< std::vector<reco::Vertex> >( iConfig.getUntrackedParameter<edm::InputTag>("vertexCollection") ) )
35  , recoVertexCompositeCandidateCollection_k0s_Token_( consumes<reco::VertexCompositeCandidateCollection>( iConfig.getUntrackedParameter<edm::InputTag>( "kShortCollection" ) ) )
36  , recoVertexCompositeCandidateCollection_lambda_Token_( consumes<reco::VertexCompositeCandidateCollection>( iConfig.getUntrackedParameter<edm::InputTag>( "lambdaCollection" ) ) )
37 {
38 
39  genLam = genK0s = realLamFoundEff = realK0sFoundEff = lamCandFound = k0sCandFound = noTPforK0sCand = noTPforLamCand = realK0sFound = realLamFound = 0;
40 
41 }
42 
43 V0Validator::~V0Validator() {
44 
45 }
46 
47 void V0Validator::bookHistograms(DQMStore::IBooker & ibooker, edm::Run const &, edm::EventSetup const &) {
48  ibooker.cd();
49  std::string subDirName = V0Validator::dirName + "/EffFakes";
50  ibooker.setCurrentFolder(subDirName.c_str());
51 
52  ksEffVsR = ibooker.book1D("K0sEffVsR",
53  "K^{0}_{S} Efficiency vs #rho", 40, 0., 40.);
54  ksEffVsEta = ibooker.book1D("K0sEffVsEta",
55  "K^{0}_{S} Efficiency vs #eta", 40, -2.5, 2.5);
56  ksEffVsPt = ibooker.book1D("K0sEffVsPt",
57  "K^{0}_{S} Efficiency vs p_{T}", 70, 0., 20.);;
58 
59  ksTkEffVsR = ibooker.book1D("K0sTkEffVsR",
60  "K^{0}_{S} Tracking Efficiency vs #rho", 40, 0., 40.);
61  ksTkEffVsEta = ibooker.book1D("K0sTkEffVsEta",
62  "K^{0}_{S} Tracking Efficiency vs #eta", 40, -2.5, 2.5);
63  ksTkEffVsPt = ibooker.book1D("K0sTkEffVsPt",
64  "K^{0}_{S} Tracking Efficiency vs p_{T}", 70, 0., 20.);
65 
66  ksEffVsR_num = ibooker.book1D("K0sEffVsR_num",
67  "K^{0}_{S} Efficiency vs #rho", 40, 0., 40.);
68  ksEffVsEta_num = ibooker.book1D("K0sEffVsEta_num",
69  "K^{0}_{S} Efficiency vs #eta", 40, -2.5, 2.5);
70  ksEffVsPt_num = ibooker.book1D("K0sEffVsPt_num",
71  "K^{0}_{S} Efficiency vs p_{T}", 70, 0., 20.);;
72 
73  ksTkEffVsR_num = ibooker.book1D("K0sTkEffVsR_num",
74  "K^{0}_{S} Tracking Efficiency vs #rho", 40, 0., 40.);
75  ksTkEffVsEta_num = ibooker.book1D("K0sTkEffVsEta_num",
76  "K^{0}_{S} Tracking Efficiency vs #eta", 40, -2.5, 2.5);
77  ksTkEffVsPt_num = ibooker.book1D("K0sTkEffVsPt_num",
78  "K^{0}_{S} Tracking Efficiency vs p_{T}", 70, 0., 20.);;
79 
80  ksEffVsR_denom = ibooker.book1D("K0sEffVsR_denom",
81  "K^{0}_{S} Efficiency vs #rho", 40, 0., 40.);
82  ksEffVsEta_denom = ibooker.book1D("K0sEffVsEta_denom",
83  "K^{0}_{S} Efficiency vs #eta", 40, -2.5, 2.5);
84  ksEffVsPt_denom = ibooker.book1D("K0sEffVsPt_denom",
85  "K^{0}_{S} Efficiency vs p_{T}", 70, 0., 20.);;
86 
87  lamEffVsR = ibooker.book1D("LamEffVsR",
88  "#Lambda^{0} Efficiency vs #rho", 40, 0., 40.);
89  lamEffVsEta = ibooker.book1D("LamEffVsEta",
90  "#Lambda^{0} Efficiency vs #eta", 40, -2.5, 2.5);
91  lamEffVsPt = ibooker.book1D("LamEffVsPt",
92  "#Lambda^{0} Efficiency vs p_{T}", 70, 0., 20.);
93 
94  lamTkEffVsR = ibooker.book1D("LamTkEffVsR",
95  "#Lambda^{0} TrackingEfficiency vs #rho", 40, 0., 40.);
96  lamTkEffVsEta = ibooker.book1D("LamTkEffVsEta",
97  "#Lambda^{0} Tracking Efficiency vs #eta", 40, -2.5, 2.5);
98  lamTkEffVsPt = ibooker.book1D("LamTkEffVsPt",
99  "#Lambda^{0} Tracking Efficiency vs p_{T}", 70, 0., 20.);
100 
101  lamEffVsR_num = ibooker.book1D("LamEffVsR_num",
102  "#Lambda^{0} Efficiency vs #rho", 40, 0., 40.);
103  lamEffVsEta_num = ibooker.book1D("LamEffVsEta_num",
104  "#Lambda^{0} Efficiency vs #eta", 40, -2.5, 2.5);
105  lamEffVsPt_num = ibooker.book1D("LamEffVsPt_num",
106  "#Lambda^{0} Efficiency vs p_{T}", 70, 0., 20.);
107 
108  lamTkEffVsR_num = ibooker.book1D("LamTkEffVsR_num",
109  "#Lambda^{0} TrackingEfficiency vs #rho", 40, 0., 40.);
110  lamTkEffVsEta_num = ibooker.book1D("LamTkEffVsEta_num",
111  "#Lambda^{0} Tracking Efficiency vs #eta", 40, -2.5, 2.5);
112  lamTkEffVsPt_num = ibooker.book1D("LamTkEffVsPt_num",
113  "#Lambda^{0} Tracking Efficiency vs p_{T}", 70, 0., 20.);
114 
115  lamEffVsR_denom = ibooker.book1D("LamEffVsR_denom",
116  "#Lambda^{0} Efficiency vs #rho", 40, 0., 40.);
117  lamEffVsEta_denom = ibooker.book1D("LamEffVsEta_denom",
118  "#Lambda^{0} Efficiency vs #eta", 40, -2.5, 2.5);
119  lamEffVsPt_denom = ibooker.book1D("LamEffVsPt_denom",
120  "#Lambda^{0} Efficiency vs p_{T}", 70, 0., 20.);
121 
122  ksFakeVsR = ibooker.book1D("K0sFakeVsR",
123  "K^{0}_{S} Fake Rate vs #rho", 40, 0., 40.);
124  ksFakeVsEta = ibooker.book1D("K0sFakeVsEta",
125  "K^{0}_{S} Fake Rate vs #eta", 40, -2.5, 2.5);
126  ksFakeVsPt = ibooker.book1D("K0sFakeVsPt",
127  "K^{0}_{S} Fake Rate vs p_{T}", 70, 0., 20.);
128  ksTkFakeVsR = ibooker.book1D("K0sTkFakeVsR",
129  "K^{0}_{S} Tracking Fake Rate vs #rho", 40, 0., 40.);
130  ksTkFakeVsEta = ibooker.book1D("K0sTkFakeVsEta",
131  "K^{0}_{S} Tracking Fake Rate vs #eta", 40, -2.5, 2.5);
132  ksTkFakeVsPt = ibooker.book1D("K0sTkFakeVsPt",
133  "K^{0}_{S} Tracking Fake Rate vs p_{T}", 70, 0., 20.);
134 
135  ksFakeVsR_num = ibooker.book1D("K0sFakeVsR_num",
136  "K^{0}_{S} Fake Rate vs #rho", 40, 0., 40.);
137  ksFakeVsEta_num = ibooker.book1D("K0sFakeVsEta_num",
138  "K^{0}_{S} Fake Rate vs #eta", 40, -2.5, 2.5);
139  ksFakeVsPt_num = ibooker.book1D("K0sFakeVsPt_num",
140  "K^{0}_{S} Fake Rate vs p_{T}", 70, 0., 20.);
141  ksTkFakeVsR_num = ibooker.book1D("K0sTkFakeVsR_num",
142  "K^{0}_{S} Tracking Fake Rate vs #rho", 40, 0., 40.);
143  ksTkFakeVsEta_num = ibooker.book1D("K0sTkFakeVsEta_num",
144  "K^{0}_{S} Tracking Fake Rate vs #eta", 40, -2.5, 2.5);
145  ksTkFakeVsPt_num = ibooker.book1D("K0sTkFakeVsPt_num",
146  "K^{0}_{S} Tracking Fake Rate vs p_{T}", 70, 0., 20.);
147 
148  ksFakeVsR_denom = ibooker.book1D("K0sFakeVsR_denom",
149  "K^{0}_{S} Fake Rate vs #rho", 40, 0., 40.);
150  ksFakeVsEta_denom = ibooker.book1D("K0sFakeVsEta_denom",
151  "K^{0}_{S} Fake Rate vs #eta", 40, -2.5, 2.5);
152  ksFakeVsPt_denom = ibooker.book1D("K0sFakeVsPt_denom",
153  "K^{0}_{S} Fake Rate vs p_{T}", 70, 0., 20.);
154 
155  lamFakeVsR = ibooker.book1D("LamFakeVsR",
156  "#Lambda^{0} Fake Rate vs #rho", 40, 0., 40.);
157  lamFakeVsEta = ibooker.book1D("LamFakeVsEta",
158  "#Lambda^{0} Fake Rate vs #eta", 40, -2.5, 2.5);
159  lamFakeVsPt = ibooker.book1D("LamFakeVsPt",
160  "#Lambda^{0} Fake Rate vs p_{T}", 70, 0., 20.);
161  lamTkFakeVsR = ibooker.book1D("LamTkFakeVsR",
162  "#Lambda^{0} Tracking Fake Rate vs #rho", 40, 0., 40.);
163  lamTkFakeVsEta = ibooker.book1D("LamTkFakeVsEta",
164  "#Lambda^{0} Tracking Fake Rate vs #eta", 40, -2.5, 2.5);
165  lamTkFakeVsPt = ibooker.book1D("LamTkFakeVsPt",
166  "#Lambda^{0} Tracking Fake Rate vs p_{T}", 70, 0., 20.);
167 
168  lamFakeVsR_num = ibooker.book1D("LamFakeVsR_num",
169  "#Lambda^{0} Fake Rate vs #rho", 40, 0., 40.);
170  lamFakeVsEta_num = ibooker.book1D("LamFakeVsEta_num",
171  "#Lambda^{0} Fake Rate vs #eta", 40, -2.5, 2.5);
172  lamFakeVsPt_num = ibooker.book1D("LamFakeVsPt_num",
173  "#Lambda^{0} Fake Rate vs p_{T}", 70, 0., 20.);
174  lamTkFakeVsR_num = ibooker.book1D("LamTkFakeVsR_num",
175  "#Lambda^{0} Tracking Fake Rate vs #rho", 40, 0., 40.);
176  lamTkFakeVsEta_num = ibooker.book1D("LamTkFakeVsEta_num",
177  "#Lambda^{0} Tracking Fake Rate vs #eta", 40, -2.5, 2.5);
178  lamTkFakeVsPt_num = ibooker.book1D("LamTkFakeVsPt_num",
179  "#Lambda^{0} Tracking Fake Rate vs p_{T}", 70, 0., 20.);
180 
181  lamFakeVsR_denom = ibooker.book1D("LamFakeVsR_denom",
182  "#Lambda^{0} Fake Rate vs #rho", 40, 0., 40.);
183  lamFakeVsEta_denom = ibooker.book1D("LamFakeVsEta_denom",
184  "#Lambda^{0} Fake Rate vs #eta", 40, -2.5, 2.5);
185  lamFakeVsPt_denom = ibooker.book1D("LamFakeVsPt_denom",
186  "#Lambda^{0} Fake Rate vs p_{T}", 70, 0., 20.);
187 
188  ibooker.cd();
189  subDirName = dirName + "/Other";
190  ibooker.setCurrentFolder(subDirName.c_str());
191 
192  nKs = ibooker.book1D("nK0s",
193  "Number of K^{0}_{S} found per event", 60, 0., 60.);
194  nLam = ibooker.book1D("nLam",
195  "Number of #Lambda^{0} found per event", 60, 0., 60.);
196 
197  ksXResolution = ibooker.book1D("ksXResolution",
198  "Resolution of V0 decay vertex X coordinate", 50, 0., 50.);
199  ksYResolution = ibooker.book1D("ksYResolution",
200  "Resolution of V0 decay vertex Y coordinate", 50, 0., 50.);
201  ksZResolution = ibooker.book1D("ksZResolution",
202  "Resolution of V0 decay vertex Z coordinate", 50, 0., 50.);
203  lamXResolution = ibooker.book1D("lamXResolution",
204  "Resolution of V0 decay vertex X coordinate", 50, 0., 50.);
205  lamYResolution = ibooker.book1D("lamYResolution",
206  "Resolution of V0 decay vertex Y coordinate", 50, 0., 50.);
207  lamZResolution = ibooker.book1D("lamZResolution",
208  "Resolution of V0 decay vertex Z coordinate", 50, 0., 50.);
209  ksAbsoluteDistResolution = ibooker.book1D("ksRResolution",
210  "Resolution of absolute distance from primary vertex to V0 vertex",
211  100, 0., 50.);
212  lamAbsoluteDistResolution = ibooker.book1D("lamRResolution",
213  "Resolution of absolute distance from primary vertex to V0 vertex",
214  100, 0., 50.);
215 
216  ksCandStatus = ibooker.book1D("ksCandStatus",
217  "Fake type by cand status",
218  10, 0., 10.);
219  lamCandStatus = ibooker.book1D("ksCandStatus",
220  "Fake type by cand status",
221  10, 0., 10.);
222 
223  double minKsMass = 0.49767 - 0.07;
224  double maxKsMass = 0.49767 + 0.07;
225  double minLamMass = 1.1156 - 0.05;
226  double maxLamMass = 1.1156 + 0.05;
227  int ksMassNbins = 100;
228  double ksMassXmin = minKsMass;
229  double ksMassXmax = maxKsMass;
230  int lamMassNbins = 100;
231  double lamMassXmin = minLamMass;
232  double lamMassXmax = maxLamMass;
233 
234  fakeKsMass = ibooker.book1D("ksMassFake",
235  "Mass of fake K0S",
236  ksMassNbins, minKsMass, maxKsMass);
237  goodKsMass = ibooker.book1D("ksMassGood",
238  "Mass of good reco K0S",
239  ksMassNbins, minKsMass, maxKsMass);
240  fakeLamMass = ibooker.book1D("lamMassFake",
241  "Mass of fake Lambda",
242  lamMassNbins, minLamMass, maxLamMass);
243  goodLamMass = ibooker.book1D("lamMassGood",
244  "Mass of good Lambda",
245  lamMassNbins, minLamMass, maxLamMass);
246 
247  ksMassAll = ibooker.book1D("ksMassAll",
248  "Invariant mass of all K0S",
249  ksMassNbins, ksMassXmin, ksMassXmax);
250  lamMassAll = ibooker.book1D("lamMassAll",
251  "Invariant mass of all #Lambda^{0}",
252  lamMassNbins, lamMassXmin, lamMassXmax);
253 
254  ksFakeDauRadDist = ibooker.book1D("radDistFakeKs",
255  "Production radius of daughter particle of Ks fake",
256  100, 0., 15.);
257  lamFakeDauRadDist = ibooker.book1D("radDistFakeLam",
258  "Production radius of daughter particle of Lam fake",
259  100, 0., 15.);
260 }
261 
262 void V0Validator::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
263 
264  using std::cout;
265  using std::endl;
266  using namespace edm;
267  using namespace std;
268 
269  // Get event setup info, B-field and tracker geometry
270  ESHandle<MagneticField> bFieldHandle;
271  iSetup.get<IdealMagneticFieldRecord>().get(bFieldHandle);
272  ESHandle<GlobalTrackingGeometry> globTkGeomHandle;
273  iSetup.get<GlobalTrackingGeometryRecord>().get(globTkGeomHandle);
274 
275  // Make matching collections
276  Handle<reco::RecoToSimCollection > recotosimCollectionH;
277  iEvent.getByToken( recoRecoToSimCollectionToken_, recotosimCollectionH );
278 
279  Handle<reco::SimToRecoCollection> simtorecoCollectionH;
280  iEvent.getByToken( recoSimToRecoCollectionToken_, simtorecoCollectionH );
281 
282  edm::Handle<TrackingParticleCollection> TPCollectionEff ;
283  iEvent.getByToken( trackingParticleCollection_Eff_Token_, TPCollectionEff );
284  const TrackingParticleCollection& tPCeff = *( TPCollectionEff.product() );
285 
286  // Select the primary vertex, create a new reco::Vertex to hold it
287  edm::Handle< std::vector<reco::Vertex> > primaryVtxCollectionH;
288  iEvent.getByToken( vec_recoVertex_Token_, primaryVtxCollectionH );
289 
290  std::vector<reco::Vertex>::const_iterator iVtxPH = primaryVtxCollectionH->begin();
291  for(std::vector<reco::Vertex>::const_iterator iVtx = primaryVtxCollectionH->begin();
292  iVtx < primaryVtxCollectionH->end();
293  iVtx++) {
294  if(primaryVtxCollectionH->size() > 1) {
295  if(iVtx->tracksSize() > iVtxPH->tracksSize()) {
296  iVtxPH = iVtx;
297  }
298  }
299  else iVtxPH = iVtx;
300  }
301 
302  //get the V0s;
303  edm::Handle<reco::VertexCompositeCandidateCollection> k0sCollection;
304  edm::Handle<reco::VertexCompositeCandidateCollection> lambdaCollection;
305  iEvent.getByToken( recoVertexCompositeCandidateCollection_k0s_Token_, k0sCollection );
306  iEvent.getByToken( recoVertexCompositeCandidateCollection_lambda_Token_, lambdaCollection );
307 
308  //make vector of pair of trackingParticles to hold good V0 candidates
309  std::vector< pair<TrackingParticleRef, TrackingParticleRef> > trueK0s;
310  std::vector< pair<TrackingParticleRef, TrackingParticleRef> > trueLams;
311  std::vector<double> trueKsMasses;
312  std::vector<double> trueLamMasses;
313 
315  // Do fake rate calculation //
317 
318  // Kshorts
319  double numK0sFound = 0.;
320  double mass = 0.;
321  std::vector<double> radDist;
322  if ( k0sCollection->size() > 0 ) {
323  vector<reco::TrackRef> theDaughterTracks;
324  for( reco::VertexCompositeCandidateCollection::const_iterator iK0s = k0sCollection->begin();
325  iK0s != k0sCollection->end();
326  iK0s++) {
327  // Fill mass of all K0S
328  ksMassAll->Fill( iK0s->mass() );
329  // Fill values to be histogrammed
330  K0sCandpT = (sqrt( iK0s->momentum().perp2() ));
331  K0sCandEta = iK0s->momentum().eta();
332  K0sCandR = (sqrt( iK0s->vertex().perp2() ));
333  K0sCandStatus = 0;
334  mass = iK0s->mass();
335 
336  theDaughterTracks.push_back( (*(dynamic_cast<const reco::RecoChargedCandidate *> (iK0s->daughter(0)) )).track() );
337  theDaughterTracks.push_back( (*(dynamic_cast<const reco::RecoChargedCandidate *> (iK0s->daughter(1)) )).track() );
338 
339  for (int itrack = 0; itrack < 2; itrack++) {
340  K0sPiCandStatus[itrack] = 0;
341  }
342 
343  std::vector< std::pair<TrackingParticleRef, double> > tp;
344  TrackingParticleRef tpref;
345  TrackingParticleRef firstDauTP;
346  TrackingVertexRef k0sVtx;
347 
348  // Loop through K0s candidate daugher tracks
349  for(View<reco::Track>::size_type i=0; i<theDaughterTracks.size(); ++i){
350  // Found track from theDaughterTracks
351  RefToBase<reco::Track> track( theDaughterTracks.at(i) );
352 
353  if(recotosimCollectionH->find(track) != recotosimCollectionH->end()) {
354  tp = (*recotosimCollectionH)[track];
355  if (tp.size() != 0) {
356  K0sPiCandStatus[i] = 1;
357  tpref = tp.begin()->first;
358 
359  if( simtorecoCollectionH->find(tpref) == simtorecoCollectionH->end() ) {
360  K0sPiCandStatus[i] = 3;
361  }
362  TrackingVertexRef parentVertex = tpref->parentVertex();
363  if(parentVertex.isNonnull()) radDist.push_back(parentVertex->position().R());
364 
365  if( parentVertex.isNonnull() ) {
366  if( k0sVtx.isNonnull() ) {
367  if( k0sVtx->position() == parentVertex->position() ) {
368  if( parentVertex->nDaughterTracks() == 2 ) {
369  if( parentVertex->nSourceTracks() == 0 ) {
370  // No source tracks found for K0s vertex; shouldn't happen, but does for evtGen events
371  K0sCandStatus = 6;
372  }
373 
374  for( TrackingVertex::tp_iterator iTP = parentVertex->sourceTracks_begin();
375  iTP != parentVertex->sourceTracks_end(); iTP++) {
376  if( (*iTP)->pdgId() == 310 ) {
377  K0sCandStatus = 1;
378  realK0sFound++;
379  numK0sFound += 1.;
380  std::pair<TrackingParticleRef, TrackingParticleRef> pair(firstDauTP, tpref);
381  // Pushing back a good V0
382  trueK0s.push_back(pair);
383  trueKsMasses.push_back(mass);
384  }
385  else {
386  K0sCandStatus = 2;
387  if( (*iTP)->pdgId() == 3122 ) {
388  K0sCandStatus = 7;
389  }
390  }
391  }
392  }
393  else {
394  // Found a bad match because the mother has too many daughters
395  K0sCandStatus = 3;
396  }
397  }
398  else {
399  // Found a bad match because the parent vertices from the two tracks are different
400  K0sCandStatus = 4;
401  }
402  }
403  else {
404  // if k0sVtx is null, fill it with parentVertex to compare to the parentVertex from the second track
405  k0sVtx = parentVertex;
406  firstDauTP = tpref;
407  }
408  }//parent vertex is null
409  }//tp size zero
410  }
411  else {
412  K0sPiCandStatus[i] = 2;
413  noTPforK0sCand++;
414  K0sCandStatus = 5;
415  theDaughterTracks.clear();
416  }
417  }
418  theDaughterTracks.clear();
419  // fill the fake rate histograms
420  if( K0sCandStatus > 1 ) {
421  ksFakeVsR_num->Fill(K0sCandR);
422  ksFakeVsEta_num->Fill(K0sCandEta);
423  ksFakeVsPt_num->Fill(K0sCandpT);
424  ksCandStatus->Fill((float) K0sCandStatus);
425  fakeKsMass->Fill(mass);
426  for( unsigned int ndx = 0; ndx < radDist.size(); ndx++ ) {
427  ksFakeDauRadDist->Fill(radDist[ndx]);
428  }
429  }
430  if( K0sCandStatus == 5 ) {
431  ksTkFakeVsR_num->Fill(K0sCandR);
432  ksTkFakeVsEta_num->Fill(K0sCandEta);
433  ksTkFakeVsPt_num->Fill(K0sCandpT);
434  }
435  ksFakeVsR_denom->Fill(K0sCandR);
436  ksFakeVsEta_denom->Fill(K0sCandEta);
437  ksFakeVsPt_denom->Fill(K0sCandpT);
438  }
439  }
440  nKs->Fill( (float) numK0sFound );
441  numK0sFound = 0.;
442 
443  double numLamFound = 0.;
444  mass = 0.;
445  radDist.clear();
446  // Lambdas
447  if ( lambdaCollection->size() > 0 ) {
448  vector<reco::TrackRef> theDaughterTracks;
449  for( reco::VertexCompositeCandidateCollection::const_iterator iLam = lambdaCollection->begin();
450  iLam != lambdaCollection->end();
451  iLam++) {
452  // Fill mass plot with ALL lambdas
453  lamMassAll->Fill( iLam->mass() );
454  // Fill values to be histogrammed
455  LamCandpT = (sqrt( iLam->momentum().perp2() ));
456  LamCandEta = iLam->momentum().eta();
457  LamCandR = (sqrt( iLam->vertex().perp2() ));
458  LamCandStatus = 0;
459  mass = iLam->mass();
460 
461  //cout << "Lam daughter tracks" << endl;
462  theDaughterTracks.push_back( (*(dynamic_cast<const reco::RecoChargedCandidate *> (iLam->daughter(0)) )).track() );
463  theDaughterTracks.push_back( (*(dynamic_cast<const reco::RecoChargedCandidate *> (iLam->daughter(1)) )).track() );
464 
465  for (int itrack = 0; itrack < 2; itrack++) {
466  LamPiCandStatus[itrack] = 0;
467  }
468 
469  std::vector< std::pair<TrackingParticleRef, double> > tp;
470  TrackingParticleRef tpref;
471  TrackingParticleRef firstDauTP;
472  TrackingVertexRef LamVtx;
473  // Loop through Lambda candidate daughter tracks
474  for(View<reco::Track>::size_type i=0; i<theDaughterTracks.size(); ++i){
475  // Found track from theDaughterTracks
476  RefToBase<reco::Track> track( theDaughterTracks.at(i) );
477 
478  if(recotosimCollectionH->find(track) != recotosimCollectionH->end()) {
479  tp = (*recotosimCollectionH)[track];
480  if (tp.size() != 0) {
481  LamPiCandStatus[i] = 1;
482  tpref = tp.begin()->first;
483  if( simtorecoCollectionH->find(tpref) == simtorecoCollectionH->end() ) {
484  LamPiCandStatus[i] = 3;
485  }
486  TrackingVertexRef parentVertex = tpref->parentVertex();
487  if( parentVertex.isNonnull() ) radDist.push_back(parentVertex->position().R());
488 
489  if( parentVertex.isNonnull() ) {
490  if( LamVtx.isNonnull() ) {
491  if( LamVtx->position() == parentVertex->position() ) {
492  if( parentVertex->nDaughterTracks() == 2 ) {
493  if( parentVertex->nSourceTracks() == 0 ) {
494  // No source tracks found for K0s vertex; shouldn't happen, but does for evtGen events
495  LamCandStatus = 6;
496  }
497 
498  for( TrackingVertex::tp_iterator iTP = parentVertex->sourceTracks_begin();
499  iTP != parentVertex->sourceTracks_end(); ++iTP) {
500  if( abs((*iTP)->pdgId()) == 3122 ) {
501  LamCandStatus = 1;
502  realLamFound++;
503  numLamFound += 1.;
504  std::pair<TrackingParticleRef, TrackingParticleRef> pair(firstDauTP, tpref);
505  // Pushing back a good V0
506  trueLams.push_back(pair);
507  trueLamMasses.push_back(mass);
508  }
509  else {
510  LamCandStatus = 2;
511  if( abs((*iTP)->pdgId() ) == 310 ) {
512  LamCandStatus = 7;
513  }
514  }
515  }
516  }
517  else {
518  // Found a bad match because the mother has too many daughters
519  LamCandStatus = 3;
520  }
521  }
522  else {
523  // Found a bad match because the parent vertices from the two tracks are different
524  LamCandStatus = 4;
525  }
526  }
527  else {
528  // if lamVtx is null, fill it with parentVertex to compare to the parentVertex from the second track
529  LamVtx = parentVertex;
530  firstDauTP = tpref;
531  }
532  }//parent vertex is null
533  }//tp size zero
534  }
535  else {
536  LamPiCandStatus[i] = 2;
537  noTPforLamCand++;
538  LamCandStatus = 5;
539  theDaughterTracks.clear();
540  }
541  }
542  theDaughterTracks.clear();
543  // fill the fake rate histograms
544  if( LamCandStatus > 1 ) {
545  lamFakeVsR_num->Fill(LamCandR);
546  lamFakeVsEta_num->Fill(LamCandEta);
547  lamFakeVsPt_num->Fill(LamCandpT);
548  lamCandStatus->Fill((float) LamCandStatus);
549  fakeLamMass->Fill(mass);
550  for( unsigned int ndx = 0; ndx < radDist.size(); ndx++ ) {
551  lamFakeDauRadDist->Fill(radDist[ndx]);
552  }
553  }
554  if( K0sCandStatus == 5 ) {
555  lamTkFakeVsR_num->Fill(LamCandR);
556  lamTkFakeVsEta_num->Fill(LamCandEta);
557  lamTkFakeVsPt_num->Fill(LamCandpT);
558  }
559  lamFakeVsR_denom->Fill(LamCandR);
560  lamFakeVsEta_denom->Fill(LamCandEta);
561  lamFakeVsPt_denom->Fill(LamCandpT);
562  }
563  }
564  nLam->Fill( (double) numLamFound );
565  numLamFound = 0.;
566 
567 
569  // Do efficiency calculation //
571  // Lambdas
572  for(TrackingParticleCollection::size_type i = 0; i < tPCeff.size(); i++) {
573  TrackingParticleRef tpr1(TPCollectionEff, i);
574  const TrackingParticle* itp1 = tpr1.get();
575  if( (itp1->pdgId() == 211 || itp1->pdgId() == 2212)
576  && itp1->parentVertex().isNonnull()
577  && abs(itp1->momentum().eta()) < 2.4
578  && sqrt( itp1->momentum().perp2() ) > 0.9) {
579  bool isLambda = false;
580  if( itp1->pdgId() == 2212 ) isLambda = true;
581  if( itp1->parentVertex()->nDaughterTracks() == 2 ) {
582 
583  TrackingVertexRef piCand1Vertex = itp1->parentVertex();
584  for(TrackingVertex::tp_iterator iTP1 = piCand1Vertex->sourceTracks_begin();
585  iTP1 != piCand1Vertex->sourceTracks_end(); iTP1++) {
586  if( abs((*iTP1)->pdgId()) == 3122 ) {
587  // ----->>>>>>Keep going here
588  for(TrackingParticleCollection::size_type j=0;
589  j < tPCeff.size();
590  j++) {
591  TrackingParticleRef tpr2(TPCollectionEff, j);
592  const TrackingParticle* itp2 = tpr2.get();
593  int particle2pdgId;
594  if (isLambda) particle2pdgId = -211;
595  else particle2pdgId = -2212;
596  if( itp2->pdgId() == particle2pdgId
597  && itp2->parentVertex().isNonnull()
598  && abs(itp2->momentum().eta()) < 2.4
599  && sqrt(itp2->momentum().perp2()) > 0.9) {
600  if(itp2->parentVertex() == itp1->parentVertex()) {
601  // Found a good pair of Lambda daughters
602  TrackingVertexRef piCand2Vertex = itp2->parentVertex();
603  for (TrackingVertex::tp_iterator iTP2 = piCand2Vertex->sourceTracks_begin();
604  iTP2 != piCand2Vertex->sourceTracks_end();
605  ++iTP2) {
606  LamGenEta = LamGenpT = LamGenR = 0.;
607  LamGenStatus = 0;
608  for(int ifill = 0;
609  ifill < 2;
610  ifill++) {
611  // do nothing?
612  }
613  if( abs((*iTP2)->pdgId()) == 3122 ) {
614  // found generated Lambda
615  LamGenpT = sqrt((*iTP2)->momentum().perp2());
616  LamGenEta = (*iTP2)->momentum().eta();
617  LamGenR = sqrt(itp2->vertex().perp2());
618  genLam++;
619  if(trueLams.size() > 0) {
620  int loop_1 = 0;
621  for(std::vector< pair<TrackingParticleRef, TrackingParticleRef> >::const_iterator iEffCheck = trueLams.begin();
622  iEffCheck != trueLams.end();
623  iEffCheck++) {
624  if( itp1->parentVertex() == iEffCheck->first->parentVertex()
625  && itp2->parentVertex() == iEffCheck->second->parentVertex() ) {
626  realLamFoundEff++;
627  //V0Producer found the generated Lambda
628  LamGenStatus = 1;
629  goodLamMass->Fill(trueLamMasses[loop_1]);
630  break;
631  }
632  else {
633  //V0Producer didn't find the generated Lambda
634  LamGenStatus = 2;
635  }
636  loop_1++;
637  }
638  }
639  else {
640  //No V0 cand found, so V0Producer didn't find the generated Lambda
641  LamGenStatus = 2;
642  }
643  std::vector< std::pair<RefToBase<reco::Track>, double> > rt1;
644  std::vector< std::pair<RefToBase<reco::Track>, double> > rt2;
645 
646  if( simtorecoCollectionH->find(tpr1) != simtorecoCollectionH->end() ) {
647  rt1 = (std::vector<std::pair<RefToBase<reco::Track>, double> >) (*simtorecoCollectionH)[tpr1];
648  if(rt1.size() != 0) {
649  LamPiEff[0] = 1; //Found the first daughter track
650  edm::RefToBase<reco::Track> t1 = rt1.begin()->first;
651  }
652  }
653  else {
654  LamPiEff[0] = 2;//First daughter not found
655  }
656  if( (simtorecoCollectionH->find(tpr2) != simtorecoCollectionH->end()) ) {
657  rt2 = (std::vector<std::pair<RefToBase<reco::Track>, double> >) (*simtorecoCollectionH)[tpr2];
658  if(rt2.size() != 0) {
659  LamPiEff[1] = 1;//Found the second daughter track
660  edm::RefToBase<reco::Track> t2 = rt2.begin()->first;
661  }
662  }
663  else {
664  LamPiEff[1] = 2;//Second daughter not found
665  }
666 
667  if( LamGenStatus == 1
668  && (LamPiEff[0] == 2 || LamPiEff[1] == 2) ) {
669  // Good Lambda found, but recoTrack->trackingParticle->recoTrack didn't work
670  LamGenStatus = 4;
671  realLamFoundEff--;
672  }
673  if( LamGenStatus == 2
674  && (LamPiEff[0] == 2 || LamPiEff[1] == 2) ) {
675  // Lambda not found because we didn't find a daughter track
676  LamGenStatus = 3;
677  }
678  // Fill histograms
679  if(LamGenR > 0.) {
680  if(LamGenStatus == 1) {
681  lamEffVsR_num->Fill(LamGenR);
682  }
683  if((double) LamGenStatus < 2.5) {
684  lamTkEffVsR_num->Fill(LamGenR);
685  }
686  lamEffVsR_denom->Fill(LamGenR);
687  }
688  if(abs(LamGenEta) > 0.) {
689  if(LamGenStatus == 1) {
690  lamEffVsEta_num->Fill(LamGenEta);
691  }
692  if((double) LamGenStatus < 2.5) {
693  lamTkEffVsEta_num->Fill(LamGenEta);
694  }
695  lamEffVsEta_denom->Fill(LamGenEta);
696  }
697  if(LamGenpT > 0.) {
698  if(LamGenStatus == 1) {
699  lamEffVsPt_num->Fill(LamGenpT);
700  }
701  if((double) LamGenStatus < 2.5) {
702  lamTkEffVsPt_num->Fill(LamGenpT);
703  }
704  lamEffVsPt_denom->Fill(LamGenpT);
705  }
706  }
707  }
708  }
709  }
710  }
711  }
712  }
713  }
714  }
715  }
716 
717  //Kshorts
718 
719  for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
720  TrackingParticleRef tpr1(TPCollectionEff, i);
721  const TrackingParticle* itp1 = tpr1.get();
722  // only count the efficiency for pions with |eta|<2.4 and pT>0.9 GeV. First search for a suitable pi+
723  if ( itp1->pdgId() == 211
724  && itp1->parentVertex().isNonnull()
725  && abs(itp1->momentum().eta()) < 2.4
726  && sqrt(itp1->momentum().perp2()) > 0.9) {
727  if ( itp1->parentVertex()->nDaughterTracks() == 2 ) {
728  TrackingVertexRef piCand1Vertex = itp1->parentVertex();
729  //check trackingParticle pion for a Ks mother
730  for (TrackingVertex::tp_iterator iTP1 = piCand1Vertex->sourceTracks_begin();
731  iTP1 != piCand1Vertex->sourceTracks_end(); ++iTP1) {
732  //iTP1 is a TrackingParticle
733  if ( (*iTP1)->pdgId()==310 ) {
734  //with a Ks mother found for the pi+, loop through trackingParticles again to find a pi-
735  for (TrackingParticleCollection::size_type j=0; j<tPCeff.size(); j++){
736  TrackingParticleRef tpr2(TPCollectionEff, j);
737  const TrackingParticle* itp2 = tpr2.get();
738 
739  if ( itp2->pdgId() == -211 && itp2->parentVertex().isNonnull()
740  && abs(itp2->momentum().eta()) < 2.4
741  && sqrt(itp2->momentum().perp2()) > 0.9) {
742  //check the pi+ and pi- have the same vertex
743  if ( itp2->parentVertex() == itp1->parentVertex() ) {
744  TrackingVertexRef piCand2Vertex = itp2->parentVertex();
745  for (TrackingVertex::tp_iterator iTP2 = piCand2Vertex->sourceTracks_begin();
746  iTP2 != piCand2Vertex->sourceTracks_end(); ++iTP2) {
747  //iTP2 is a TrackingParticle
748  K0sGenEta = K0sGenpT = K0sGenR = 0.;
749  K0sGenStatus = 0;
750  if( (*iTP2)->pdgId() == 310 ) {
751  K0sGenpT = sqrt( (*iTP2)->momentum().perp2() );
752  K0sGenEta = (*iTP2)->momentum().eta();
753  K0sGenR = sqrt(itp2->vertex().perp2());
754  genK0s++;
755  int loop_2 = 0;
756  if( trueK0s.size() > 0 ) {
757  for( std::vector< pair<TrackingParticleRef, TrackingParticleRef> >::const_iterator iEffCheck = trueK0s.begin();
758  iEffCheck != trueK0s.end();
759  iEffCheck++) {
760  //if the parent vertices for the tracks are the same, then the generated Ks was found
761  if (itp1->parentVertex()==iEffCheck->first->parentVertex() &&
762  itp2->parentVertex()==iEffCheck->second->parentVertex()) {
763  realK0sFoundEff++;
764  K0sGenStatus = 1;
765  goodKsMass->Fill(trueKsMasses[loop_2]);
766  break;
767  }
768  else {
769  K0sGenStatus = 2;
770  }
771  }
772  }
773  else {
774  K0sGenStatus = 2;
775  }
776 
777  // Check if the generated Ks tracks were found or not
778  // by searching the recoTracks list for a match to the trackingParticles
779 
780  std::vector<std::pair<RefToBase<reco::Track>, double> > rt1;
781  std::vector<std::pair<RefToBase<reco::Track>, double> > rt2;
782 
783  if( simtorecoCollectionH->find(tpr1) != simtorecoCollectionH->end() ) {
784  rt1 = (std::vector< std::pair<RefToBase<reco::Track>, double> >) (*simtorecoCollectionH)[tpr1];
785  if(rt1.size() != 0) {
786  //First pion found
787  K0sPiEff[0] = 1;
788  edm::RefToBase<reco::Track> t1 = rt1.begin()->first;
789  }
790  }
791  else {
792  //First pion not found
793  K0sPiEff[0] = 2;
794  }
795 
796  if( simtorecoCollectionH->find(tpr2) != simtorecoCollectionH->end() ) {
797  rt2 = (std::vector< std::pair<RefToBase<reco::Track>, double> >) (*simtorecoCollectionH)[tpr2];
798  if(rt2.size() != 0) {
799  //Second pion found
800  K0sPiEff[1] = 1;
801  edm::RefToBase<reco::Track> t2 = rt2.begin()->first;
802  }
803  }
804  else {
805  K0sPiEff[1] = 2;
806  }
807  if(K0sGenStatus == 1
808  && (K0sPiEff[0] == 2 || K0sPiEff[1] == 2)) {
809  K0sGenStatus = 4;
810  realK0sFoundEff--;
811  }
812  if(K0sGenStatus == 2
813  && (K0sPiEff[0] == 2 || K0sPiEff[1] == 2)) {
814  K0sGenStatus = 3;
815  }
816  if(K0sPiEff[0] == 1 && K0sPiEff[1] == 1) {
817  k0sTracksFound++;
818  }
819  //Fill Histograms
820  if(K0sGenR > 0.) {
821  if(K0sGenStatus == 1) {
822  ksEffVsR_num->Fill(K0sGenR);
823  }
824  if((double) K0sGenStatus < 2.5) {
825  ksTkEffVsR_num->Fill(K0sGenR);
826  }
827  ksEffVsR_denom->Fill(K0sGenR);
828  }
829  if(abs(K0sGenEta) > 0.) {
830  if(K0sGenStatus == 1) {
831  ksEffVsEta_num->Fill(K0sGenEta);
832  }
833  if((double) K0sGenStatus < 2.5) {
834  ksTkEffVsEta_num->Fill(K0sGenEta);
835  }
836  ksEffVsEta_denom->Fill(K0sGenEta);
837  }
838  if(K0sGenpT > 0.) {
839  if(K0sGenStatus == 1) {
840  ksEffVsPt_num->Fill(K0sGenpT);
841  }
842  if((double) K0sGenStatus < 2.5) {
843  ksTkEffVsPt_num->Fill(K0sGenpT);
844  }
845  ksEffVsPt_denom->Fill(K0sGenpT);
846  }
847  }
848  }
849  }
850  }
851  }
852  }
853  }
854  }
855  }
856  }
857 }
858 
859 //define this as a plug-in
860 //DEFINE_FWK_MODULE(V0Validator);
V0Validator::ksTkFakeVsPt
MonitorElement * ksTkFakeVsPt
Definition: V0Validator.h:115
edm::View::size_type
unsigned int size_type
Definition: View.h:89
i
int i
Definition: DBlmapReader.cc:9
V0Validator::nLam
MonitorElement * nLam
Definition: V0Validator.h:180
V0Validator::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: V0Validator.cc:47
reco::VertexCompositeCandidateCollection
std::vector< VertexCompositeCandidate > VertexCompositeCandidateCollection
collection of Candidate objects
Definition: VertexCompositeCandidateFwd.h:19
edm::Ref::isNonnull
bool isNonnull() const
Checks for non-null.
Definition: Ref.h:252
V0Validator::lamTkFakeVsR_num
MonitorElement * lamTkFakeVsR_num
Definition: V0Validator.h:166
V0Validator::ksEffVsR
MonitorElement * ksEffVsR
Definition: V0Validator.h:104
V0Validator::K0sCandStatus
unsigned int K0sCandStatus
Definition: V0Validator.h:96
V0Validator::recoSimToRecoCollectionToken_
edm::EDGetTokenT< reco::SimToRecoCollection > recoSimToRecoCollectionToken_
Definition: V0Validator.h:201
V0Validator::V0Validator
V0Validator(const edm::ParameterSet &)
Definition: V0Validator.cc:28
V0Validator::recoVertexCompositeCandidateCollection_k0s_Token_
edm::EDGetTokenT< reco::VertexCompositeCandidateCollection > recoVertexCompositeCandidateCollection_k0s_Token_
Definition: V0Validator.h:204
V0Validator::LamGenR
float LamGenR
Definition: V0Validator.h:93
TrackingParticleCollection
std::vector< TrackingParticle > TrackingParticleCollection
Definition: HiEvtPlaneFlatProducer.cc:111
TrackingParticle::momentum
Vector momentum() const
spatial momentum vector
Definition: TrackingParticle.h:109
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
DQMStore::IBooker::cd
void cd(void)
Definition: DQMStore.cc:266
V0Validator::ksFakeVsPt_denom
MonitorElement * ksFakeVsPt_denom
Definition: V0Validator.h:132
V0Validator::LamCandR
float LamCandR
Definition: V0Validator.h:95
V0Validator::lamTkEffVsR
MonitorElement * lamTkEffVsR
Definition: V0Validator.h:147
TrackingParticle::pdgId
int pdgId() const
PDG ID.
Definition: TrackingParticle.h:61
V0Validator::lamFakeVsPt
MonitorElement * lamFakeVsPt
Definition: V0Validator.h:152
V0Validator::ksTkFakeVsEta
MonitorElement * ksTkFakeVsEta
Definition: V0Validator.h:114
V0Validator::ksFakeVsR_denom
MonitorElement * ksFakeVsR_denom
Definition: V0Validator.h:130
V0Validator::ksFakeVsEta_num
MonitorElement * ksFakeVsEta_num
Definition: V0Validator.h:124
V0Validator::LamPiCandStatus
unsigned int LamPiCandStatus[2]
Definition: V0Validator.h:97
V0Validator::lamEffVsEta_num
MonitorElement * lamEffVsEta_num
Definition: V0Validator.h:158
V0Validator::lamFakeVsPt_denom
MonitorElement * lamFakeVsPt_denom
Definition: V0Validator.h:139
V0Validator::ksEffVsR_num
MonitorElement * ksEffVsR_num
Definition: V0Validator.h:117
V0Validator::ksFakeVsPt_num
MonitorElement * ksFakeVsPt_num
Definition: V0Validator.h:125
V0Validator::lamFakeVsR_denom
MonitorElement * lamFakeVsR_denom
Definition: V0Validator.h:137
V0Validator::lamEffVsPt_denom
MonitorElement * lamEffVsPt_denom
Definition: V0Validator.h:142
V0Validator::K0sPiCandStatus
unsigned int K0sPiCandStatus[2]
Definition: V0Validator.h:97
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:106
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:19
V0Validator::lamEffVsPt_num
MonitorElement * lamEffVsPt_num
Definition: V0Validator.h:159
HLT_25ns14e33_v1_cff.InputTag
tuple InputTag
Definition: HLT_25ns14e33_v1_cff.py:45333
V0Validator::lamTkFakeVsR
MonitorElement * lamTkFakeVsR
Definition: V0Validator.h:153
V0Validator::lamFakeVsR
MonitorElement * lamFakeVsR
Definition: V0Validator.h:150
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
V0Validator::ksFakeVsEta
MonitorElement * ksFakeVsEta
Definition: V0Validator.h:111
V0Validator::lamZResolution
MonitorElement * lamZResolution
Definition: V0Validator.h:176
V0Validator::LamPiEff
unsigned int LamPiEff[2]
Definition: V0Validator.h:97
V0Validator::lamFakeDauRadDist
MonitorElement * lamFakeDauRadDist
Definition: V0Validator.h:195
V0Validator::lamTkEffVsEta
MonitorElement * lamTkEffVsEta
Definition: V0Validator.h:148
V0Validator::lamEffVsR
MonitorElement * lamEffVsR
Definition: V0Validator.h:144
V0Validator::ksFakeVsEta_denom
MonitorElement * ksFakeVsEta_denom
Definition: V0Validator.h:131
iEvent
int iEvent
Definition: GenABIO.cc:230
V0Validator::fakeLamMass
MonitorElement * fakeLamMass
Definition: V0Validator.h:187
V0Validator::ksFakeDauRadDist
MonitorElement * ksFakeDauRadDist
Definition: V0Validator.h:194
V0Validator::lamTkFakeVsEta_num
MonitorElement * lamTkFakeVsEta_num
Definition: V0Validator.h:167
GenVertexRefVector
edm::RefVector< edm::HepMCProduct, HepMC::GenVertex > GenVertexRefVector
Definition: TrackingVertex.cc:5
V0Validator::lamTkFakeVsPt
MonitorElement * lamTkFakeVsPt
Definition: V0Validator.h:155
V0Validator::ksEffVsPt_denom
MonitorElement * ksEffVsPt_denom
Definition: V0Validator.h:135
V0Validator::lamCandFound
int lamCandFound
Definition: V0Validator.h:101
V0Validator::ksTkFakeVsEta_num
MonitorElement * ksTkFakeVsEta_num
Definition: V0Validator.h:127
V0Validator::k0sCandFound
int k0sCandFound
Definition: V0Validator.h:101
V0Validator::noTPforLamCand
int noTPforLamCand
Definition: V0Validator.h:101
TrackerOfflineValidation_Dqm_cff.dirName
tuple dirName
Definition: TrackerOfflineValidation_Dqm_cff.py:56
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
V0Validator::ksTkEffVsR_num
MonitorElement * ksTkEffVsR_num
Definition: V0Validator.h:120
V0Validator::lamFakeVsEta_denom
MonitorElement * lamFakeVsEta_denom
Definition: V0Validator.h:138
V0Validator::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: V0Validator.cc:262
V0Validator::lamXResolution
MonitorElement * lamXResolution
Definition: V0Validator.h:174
V0Validator::ksEffVsEta_denom
MonitorElement * ksEffVsEta_denom
Definition: V0Validator.h:134
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
V0Validator::lamTkEffVsPt
MonitorElement * lamTkEffVsPt
Definition: V0Validator.h:149
V0Validator::K0sCandEta
float K0sCandEta
Definition: V0Validator.h:95
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
j
int j
Definition: DBlmapReader.cc:9
V0Validator::goodKsMass
MonitorElement * goodKsMass
Definition: V0Validator.h:186
edm::Ref::get
T const * get() const
Returns C++ pointer to the item.
Definition: Ref.h:244
V0Validator::LamGenpT
float LamGenpT
Definition: V0Validator.h:93
V0Validator::realLamFound
int realLamFound
Definition: V0Validator.h:100
TrackingParticle::parentVertex
const TrackingVertexRef & parentVertex() const
Definition: TrackingParticle.h:90
V0Validator::fakeKsMass
MonitorElement * fakeKsMass
Definition: V0Validator.h:185
V0Validator::vec_recoVertex_Token_
edm::EDGetTokenT< std::vector< reco::Vertex > > vec_recoVertex_Token_
Definition: V0Validator.h:203
V0Validator::ksTkEffVsPt
MonitorElement * ksTkEffVsPt
Definition: V0Validator.h:109
V0Validator::goodLamMass
MonitorElement * goodLamMass
Definition: V0Validator.h:188
V0Validator::ksFakeVsPt
MonitorElement * ksFakeVsPt
Definition: V0Validator.h:112
TrackingVertexRef
edm::Ref< TrackingVertexCollection > TrackingVertexRef
Definition: TrackingVertexContainer.h:9
V0Validator::LamGenStatus
unsigned int LamGenStatus
Definition: V0Validator.h:96
V0Validator::lamTkFakeVsEta
MonitorElement * lamTkFakeVsEta
Definition: V0Validator.h:154
V0Validator::K0sGenStatus
unsigned int K0sGenStatus
Definition: V0Validator.h:96
V0Validator::ksMassAll
MonitorElement * ksMassAll
Definition: V0Validator.h:190
V0Validator::lamTkEffVsR_num
MonitorElement * lamTkEffVsR_num
Definition: V0Validator.h:160
V0Validator::ksTkEffVsR
MonitorElement * ksTkEffVsR
Definition: V0Validator.h:107
V0Validator::nKs
MonitorElement * nKs
Definition: V0Validator.h:179
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
V0Validator::lamTkEffVsEta_num
MonitorElement * lamTkEffVsEta_num
Definition: V0Validator.h:161
V0Validator::noTPforK0sCand
int noTPforK0sCand
Definition: V0Validator.h:101
V0Validator::K0sGenEta
float K0sGenEta
Definition: V0Validator.h:93
V0Validator::ksXResolution
MonitorElement * ksXResolution
Definition: V0Validator.h:170
V0Validator::K0sCandpT
float K0sCandpT
Definition: V0Validator.h:95
V0Validator::K0sGenpT
float K0sGenpT
Definition: V0Validator.h:93
V0Validator::lamEffVsEta_denom
MonitorElement * lamEffVsEta_denom
Definition: V0Validator.h:141
TrackingVertexCollection
std::vector< TrackingVertex > TrackingVertexCollection
Definition: TrackingVertexContainer.h:8
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
V0Validator::lamFakeVsEta
MonitorElement * lamFakeVsEta
Definition: V0Validator.h:151
V0Validator::ksAbsoluteDistResolution
MonitorElement * ksAbsoluteDistResolution
Definition: V0Validator.h:173
reco::RecoToSimCollection
edm::AssociationMap< edm::OneToManyWithQualityGeneric< edm::View< reco::Track >, TrackingParticleCollection, double > > RecoToSimCollection
Definition: TrackAssociation.h:19
V0Validator::trackingParticleCollection_Eff_Token_
edm::EDGetTokenT< TrackingParticleCollection > trackingParticleCollection_Eff_Token_
Definition: V0Validator.h:202
V0Validator::ksCandStatus
MonitorElement * ksCandStatus
Definition: V0Validator.h:182
V0Validator::lamCandStatus
MonitorElement * lamCandStatus
Definition: V0Validator.h:183
V0Validator::lamEffVsEta
MonitorElement * lamEffVsEta
Definition: V0Validator.h:145
V0Validator::lamTkFakeVsPt_num
MonitorElement * lamTkFakeVsPt_num
Definition: V0Validator.h:168
V0Validator::lamFakeVsEta_num
MonitorElement * lamFakeVsEta_num
Definition: V0Validator.h:164
TrackingParticle::vertex
Point vertex() const
Parent vertex position.
Definition: TrackingParticle.h:199
V0Validator::ksTkEffVsPt_num
MonitorElement * ksTkEffVsPt_num
Definition: V0Validator.h:122
V0Validator::lamFakeVsPt_num
MonitorElement * lamFakeVsPt_num
Definition: V0Validator.h:165
V0Validator::lamMassAll
MonitorElement * lamMassAll
Definition: V0Validator.h:191
V0Validator::K0sCandR
float K0sCandR
Definition: V0Validator.h:95
V0Validator::lamFakeVsR_num
MonitorElement * lamFakeVsR_num
Definition: V0Validator.h:163
V0Validator::LamGenEta
float LamGenEta
Definition: V0Validator.h:93
V0Validator::LamCandEta
float LamCandEta
Definition: V0Validator.h:95
V0Validator::lamEffVsR_num
MonitorElement * lamEffVsR_num
Definition: V0Validator.h:157
V0Validator::ksFakeVsR
MonitorElement * ksFakeVsR
Definition: V0Validator.h:110
V0Validator::ksEffVsEta
MonitorElement * ksEffVsEta
Definition: V0Validator.h:105
V0Validator::lamTkEffVsPt_num
MonitorElement * lamTkEffVsPt_num
Definition: V0Validator.h:162
V0Validator::recoRecoToSimCollectionToken_
edm::EDGetTokenT< reco::RecoToSimCollection > recoRecoToSimCollectionToken_
Definition: V0Validator.h:200
V0Validator::recoVertexCompositeCandidateCollection_lambda_Token_
edm::EDGetTokenT< reco::VertexCompositeCandidateCollection > recoVertexCompositeCandidateCollection_lambda_Token_
Definition: V0Validator.h:204
reco::SimToRecoCollection
edm::AssociationMap< edm::OneToManyWithQualityGeneric< TrackingParticleCollection, edm::View< reco::Track >, double > > SimToRecoCollection
Definition: TrackAssociation.h:16
V0Validator::lamAbsoluteDistResolution
MonitorElement * lamAbsoluteDistResolution
Definition: V0Validator.h:177
V0Validator::ksZResolution
MonitorElement * ksZResolution
Definition: V0Validator.h:172
V0Validator::LamCandStatus
unsigned int LamCandStatus
Definition: V0Validator.h:96
TrackingParticle
Monte Carlo truth information used for tracking validation.
Definition: TrackingParticle.h:29
V0Validator::ksEffVsPt_num
MonitorElement * ksEffVsPt_num
Definition: V0Validator.h:119
V0Validator::ksTkEffVsEta_num
MonitorElement * ksTkEffVsEta_num
Definition: V0Validator.h:121
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
V0Validator::lamEffVsPt
MonitorElement * lamEffVsPt
Definition: V0Validator.h:146
V0Validator::ksEffVsR_denom
MonitorElement * ksEffVsR_denom
Definition: V0Validator.h:133
V0Validator::ksTkFakeVsPt_num
MonitorElement * ksTkFakeVsPt_num
Definition: V0Validator.h:128
V0Validator::realK0sFoundEff
int realK0sFoundEff
Definition: V0Validator.h:100
V0Validator::ksTkEffVsEta
MonitorElement * ksTkEffVsEta
Definition: V0Validator.h:108
V0Validator::ksFakeVsR_num
MonitorElement * ksFakeVsR_num
Definition: V0Validator.h:123
V0Validator::ksEffVsEta_num
MonitorElement * ksEffVsEta_num
Definition: V0Validator.h:118
V0Validator::lamYResolution
MonitorElement * lamYResolution
Definition: V0Validator.h:175
V0Validator::dirName
std::string dirName
Definition: V0Validator.h:199
V0Validator::ksEffVsPt
MonitorElement * ksEffVsPt
Definition: V0Validator.h:106
V0Validator::K0sPiEff
unsigned int K0sPiEff[2]
Definition: V0Validator.h:97
V0Validator::k0sTracksFound
int k0sTracksFound
Definition: V0Validator.h:101
GenParticleRefVector
edm::RefVector< edm::HepMCProduct, HepMC::GenParticle > GenParticleRefVector
Definition: V0Validator.cc:26
V0Validator::realLamFoundEff
int realLamFoundEff
Definition: V0Validator.h:100
V0Validator::K0sGenR
float K0sGenR
Definition: V0Validator.h:93
V0Validator::ksTkFakeVsR
MonitorElement * ksTkFakeVsR
Definition: V0Validator.h:113
V0Validator::ksYResolution
MonitorElement * ksYResolution
Definition: V0Validator.h:171
V0Validator::lamEffVsR_denom
MonitorElement * lamEffVsR_denom
Definition: V0Validator.h:140
edm::Run
Definition: Run.h:41
V0Validator::realK0sFound
int realK0sFound
Definition: V0Validator.h:100
V0Validator::LamCandpT
float LamCandpT
Definition: V0Validator.h:95
V0Validator::ksTkFakeVsR_num
MonitorElement * ksTkFakeVsR_num
Definition: V0Validator.h:126
Generated for CMSSW Reference Manual by doxygen 1.8.5