CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
MTVHistoProducerAlgoForTracker.cc
Go to the documentation of this file.
1 #include "Validation/RecoTrack/interface/MTVHistoProducerAlgoForTracker.h"
2 #include "FWCore/MessageLogger/interface/MessageLogger.h"
3 #include "DataFormats/TrackReco/interface/Track.h"
4 #include "FWCore/ParameterSet/interface/ParameterSet.h"
5 
6 #include "DataFormats/TrackReco/interface/DeDxData.h"
7 #include "DataFormats/Common/interface/ValueMap.h"
8 #include "DataFormats/Common/interface/Ref.h"
9 #include <DataFormats/TrackReco/interface/TrackFwd.h>
10 
11 #include "DQMServices/ClientConfig/interface/FitSlicesYTool.h"
12 
13 #include "TMath.h"
14 #include <TF1.h>
15 
16 using namespace std;
17 
18 MTVHistoProducerAlgoForTracker::MTVHistoProducerAlgoForTracker(const edm::ParameterSet& pset): MTVHistoProducerAlgo(pset){
19  //parameters for _vs_eta plots
20  minEta = pset.getParameter<double>("minEta");
21  maxEta = pset.getParameter<double>("maxEta");
22  nintEta = pset.getParameter<int>("nintEta");
23  useFabsEta = pset.getParameter<bool>("useFabsEta");
24 
25  //parameters for _vs_pt plots
26  minPt = pset.getParameter<double>("minPt");
27  maxPt = pset.getParameter<double>("maxPt");
28  nintPt = pset.getParameter<int>("nintPt");
29  useInvPt = pset.getParameter<bool>("useInvPt");
30  useLogPt = pset.getUntrackedParameter<bool>("useLogPt",false);
31 
32  //parameters for _vs_Hit plots
33  minHit = pset.getParameter<double>("minHit");
34  maxHit = pset.getParameter<double>("maxHit");
35  nintHit = pset.getParameter<int>("nintHit");
36 
37  //parameters for _vs_phi plots
38  minPhi = pset.getParameter<double>("minPhi");
39  maxPhi = pset.getParameter<double>("maxPhi");
40  nintPhi = pset.getParameter<int>("nintPhi");
41 
42  //parameters for _vs_Dxy plots
43  minDxy = pset.getParameter<double>("minDxy");
44  maxDxy = pset.getParameter<double>("maxDxy");
45  nintDxy = pset.getParameter<int>("nintDxy");
46 
47  //parameters for _vs_Dz plots
48  minDz = pset.getParameter<double>("minDz");
49  maxDz = pset.getParameter<double>("maxDz");
50  nintDz = pset.getParameter<int>("nintDz");
51 
52  //parameters for _vs_ProductionVertexTransvPosition plots
53  minVertpos = pset.getParameter<double>("minVertpos");
54  maxVertpos = pset.getParameter<double>("maxVertpos");
55  nintVertpos = pset.getParameter<int>("nintVertpos");
56 
57  //parameters for _vs_ProductionVertexZPosition plots
58  minZpos = pset.getParameter<double>("minZpos");
59  maxZpos = pset.getParameter<double>("maxZpos");
60  nintZpos = pset.getParameter<int>("nintZpos");
61 
62  //parameters for dE/dx plots
63  minDeDx = pset.getParameter<double>("minDeDx");
64  maxDeDx = pset.getParameter<double>("maxDeDx");
65  nintDeDx = pset.getParameter<int>("nintDeDx");
66 
67  //parameters for Pileup plots
68  minVertcount = pset.getParameter<double>("minVertcount");
69  maxVertcount = pset.getParameter<double>("maxVertcount");
70  nintVertcount = pset.getParameter<int>("nintVertcount");
71 
72  //parameters for resolution plots
73  ptRes_rangeMin = pset.getParameter<double>("ptRes_rangeMin");
74  ptRes_rangeMax = pset.getParameter<double>("ptRes_rangeMax");
75  ptRes_nbin = pset.getParameter<int>("ptRes_nbin");
76 
77  phiRes_rangeMin = pset.getParameter<double>("phiRes_rangeMin");
78  phiRes_rangeMax = pset.getParameter<double>("phiRes_rangeMax");
79  phiRes_nbin = pset.getParameter<int>("phiRes_nbin");
80 
81  cotThetaRes_rangeMin = pset.getParameter<double>("cotThetaRes_rangeMin");
82  cotThetaRes_rangeMax = pset.getParameter<double>("cotThetaRes_rangeMax");
83  cotThetaRes_nbin = pset.getParameter<int>("cotThetaRes_nbin");
84 
85  dxyRes_rangeMin = pset.getParameter<double>("dxyRes_rangeMin");
86  dxyRes_rangeMax = pset.getParameter<double>("dxyRes_rangeMax");
87  dxyRes_nbin = pset.getParameter<int>("dxyRes_nbin");
88 
89  dzRes_rangeMin = pset.getParameter<double>("dzRes_rangeMin");
90  dzRes_rangeMax = pset.getParameter<double>("dzRes_rangeMax");
91  dzRes_nbin = pset.getParameter<int>("dzRes_nbin");
92 
93 
94  //--- tracking particle selectors for efficiency measurements
95  using namespace edm;
96 
97  ParameterSet generalTpSelectorPSet = pset.getParameter<ParameterSet>("generalTpSelector");
98  ParameterSet TpSelectorForEfficiencyVsEtaPSet = pset.getParameter<ParameterSet>("TpSelectorForEfficiencyVsEta");
99  ParameterSet TpSelectorForEfficiencyVsPhiPSet = pset.getParameter<ParameterSet>("TpSelectorForEfficiencyVsPhi");
100  ParameterSet TpSelectorForEfficiencyVsPtPSet = pset.getParameter<ParameterSet>("TpSelectorForEfficiencyVsPt");
101  ParameterSet TpSelectorForEfficiencyVsVTXRPSet = pset.getParameter<ParameterSet>("TpSelectorForEfficiencyVsVTXR");
102  ParameterSet TpSelectorForEfficiencyVsVTXZPSet = pset.getParameter<ParameterSet>("TpSelectorForEfficiencyVsVTXZ");
103 
104  ParameterSet generalGpSelectorPSet = pset.getParameter<ParameterSet>("generalGpSelector");
105  ParameterSet GpSelectorForEfficiencyVsEtaPSet = pset.getParameter<ParameterSet>("GpSelectorForEfficiencyVsEta");
106  ParameterSet GpSelectorForEfficiencyVsPhiPSet = pset.getParameter<ParameterSet>("GpSelectorForEfficiencyVsPhi");
107  ParameterSet GpSelectorForEfficiencyVsPtPSet = pset.getParameter<ParameterSet>("GpSelectorForEfficiencyVsPt");
108  ParameterSet GpSelectorForEfficiencyVsVTXRPSet = pset.getParameter<ParameterSet>("GpSelectorForEfficiencyVsVTXR");
109  ParameterSet GpSelectorForEfficiencyVsVTXZPSet = pset.getParameter<ParameterSet>("GpSelectorForEfficiencyVsVTXZ");
110 
111  ParameterSet TpSelectorForEfficiencyVsConPSet = TpSelectorForEfficiencyVsEtaPSet;
112  Entry name("signalOnly",false,true);
113  TpSelectorForEfficiencyVsConPSet.insert(true,"signalOnly",name);
114 
115  using namespace reco::modules;
116  generalTpSelector = new TrackingParticleSelector(ParameterAdapter<TrackingParticleSelector>::make(generalTpSelectorPSet));
117  TpSelectorForEfficiencyVsEta = new TrackingParticleSelector(ParameterAdapter<TrackingParticleSelector>::make(TpSelectorForEfficiencyVsEtaPSet));
118  TpSelectorForEfficiencyVsCon = new TrackingParticleSelector(ParameterAdapter<TrackingParticleSelector>::make(TpSelectorForEfficiencyVsConPSet));
119  TpSelectorForEfficiencyVsPhi = new TrackingParticleSelector(ParameterAdapter<TrackingParticleSelector>::make(TpSelectorForEfficiencyVsPhiPSet));
120  TpSelectorForEfficiencyVsPt = new TrackingParticleSelector(ParameterAdapter<TrackingParticleSelector>::make(TpSelectorForEfficiencyVsPtPSet));
121  TpSelectorForEfficiencyVsVTXR = new TrackingParticleSelector(ParameterAdapter<TrackingParticleSelector>::make(TpSelectorForEfficiencyVsVTXRPSet));
122  TpSelectorForEfficiencyVsVTXZ = new TrackingParticleSelector(ParameterAdapter<TrackingParticleSelector>::make(TpSelectorForEfficiencyVsVTXZPSet));
123 
124  generalGpSelector = new GenParticleCustomSelector(ParameterAdapter<GenParticleCustomSelector>::make(generalGpSelectorPSet));
125  GpSelectorForEfficiencyVsEta = new GenParticleCustomSelector(ParameterAdapter<GenParticleCustomSelector>::make(GpSelectorForEfficiencyVsEtaPSet));
126  GpSelectorForEfficiencyVsCon = new GenParticleCustomSelector(ParameterAdapter<GenParticleCustomSelector>::make(GpSelectorForEfficiencyVsEtaPSet));
127  GpSelectorForEfficiencyVsPhi = new GenParticleCustomSelector(ParameterAdapter<GenParticleCustomSelector>::make(GpSelectorForEfficiencyVsPhiPSet));
128  GpSelectorForEfficiencyVsPt = new GenParticleCustomSelector(ParameterAdapter<GenParticleCustomSelector>::make(GpSelectorForEfficiencyVsPtPSet));
129  GpSelectorForEfficiencyVsVTXR = new GenParticleCustomSelector(ParameterAdapter<GenParticleCustomSelector>::make(GpSelectorForEfficiencyVsVTXRPSet));
130  GpSelectorForEfficiencyVsVTXZ = new GenParticleCustomSelector(ParameterAdapter<GenParticleCustomSelector>::make(GpSelectorForEfficiencyVsVTXZPSet));
131 
132  // fix for the LogScale by Ryan
133  if(useLogPt){
134  maxPt=log10(maxPt);
135  if(minPt > 0){
136  minPt=log10(minPt);
137  }
138  else{
139  edm::LogWarning("MultiTrackValidator")
140  << "minPt = "
141  << minPt << " <= 0 out of range while requesting log scale. Using minPt = 0.1.";
142  minPt=log10(0.1);
143  }
144  }
145 
146 }
147 
148 MTVHistoProducerAlgoForTracker::~MTVHistoProducerAlgoForTracker(){
149  delete generalTpSelector;
150  delete TpSelectorForEfficiencyVsEta;
151  delete TpSelectorForEfficiencyVsCon;
152  delete TpSelectorForEfficiencyVsPhi;
153  delete TpSelectorForEfficiencyVsPt;
154  delete TpSelectorForEfficiencyVsVTXR;
155  delete TpSelectorForEfficiencyVsVTXZ;
156 
157  delete generalGpSelector;
158  delete GpSelectorForEfficiencyVsEta;
159  delete GpSelectorForEfficiencyVsCon;
160  delete GpSelectorForEfficiencyVsPhi;
161  delete GpSelectorForEfficiencyVsPt;
162  delete GpSelectorForEfficiencyVsVTXR;
163  delete GpSelectorForEfficiencyVsVTXZ;
164 }
165 
166 
167 void MTVHistoProducerAlgoForTracker::setUpVectors(){
168  std::vector<double> etaintervalsv;
169  std::vector<double> phiintervalsv;
170  std::vector<double> pTintervalsv;
171  std::vector<double> dxyintervalsv;
172  std::vector<double> dzintervalsv;
173  std::vector<double> vertposintervalsv;
174  std::vector<double> zposintervalsv;
175  std::vector<double> vertcountintervalsv;
176 
177  std::vector<int> totSIMveta,totASSveta,totASS2veta,totloopveta,totmisidveta,totASS2vetaSig,totRECveta;
178  std::vector<int> totSIMvpT,totASSvpT,totASS2vpT,totRECvpT,totloopvpT,totmisidvpT;
179  std::vector<int> totSIMv_hit,totASSv_hit,totASS2v_hit,totRECv_hit,totloopv_hit,totmisidv_hit;
180  std::vector<int> totSIMv_phi,totASSv_phi,totASS2v_phi,totRECv_phi,totloopv_phi,totmisidv_phi;
181  std::vector<int> totSIMv_dxy,totASSv_dxy,totASS2v_dxy,totRECv_dxy,totloopv_dxy,totmisidv_dxy;
182  std::vector<int> totSIMv_dz,totASSv_dz,totASS2v_dz,totRECv_dz,totloopv_dz,totmisidv_dz;
183 
184  std::vector<int> totSIMv_vertpos,totASSv_vertpos,totSIMv_zpos,totASSv_zpos;
185  std::vector<int> totSIMv_vertcount,totASSv_vertcount,totRECv_vertcount,totASS2v_vertcount;
186  std::vector<int> totRECv_algo;
187 
188  double step=(maxEta-minEta)/nintEta;
189  //std::ostringstream title,name; ///BM, what is this?
190  etaintervalsv.push_back(minEta);
191  for (int k=1;k<nintEta+1;k++) {
192  double d=minEta+k*step;
193  etaintervalsv.push_back(d);
194  totSIMveta.push_back(0);
195  totASSveta.push_back(0);
196  totASS2veta.push_back(0);
197  totloopveta.push_back(0);
198  totmisidveta.push_back(0);
199  totASS2vetaSig.push_back(0);
200  totRECveta.push_back(0);
201  }
202  etaintervals.push_back(etaintervalsv);
203  totSIMeta.push_back(totSIMveta);
204  totCONeta.push_back(totASSveta);
205  totASSeta.push_back(totASSveta);
206  totASS2eta.push_back(totASS2veta);
207  totloopeta.push_back(totloopveta);
208  totmisideta.push_back(totmisidveta);
209  totASS2etaSig.push_back(totASS2vetaSig);
210  totRECeta.push_back(totRECveta);
211  totFOMT_eta.push_back(totASSveta);
212 
213  totASS2_itpu_eta_entire.push_back(totASS2veta);
214  totASS2_itpu_eta_entire_signal.push_back(totASS2vetaSig);
215  totASS2_ootpu_eta_entire.push_back(totASS2veta);
216 
217  for (size_t i = 0; i < 15; i++) {
218  totRECv_algo.push_back(0);
219  }
220  totREC_algo.push_back(totRECv_algo);
221 
222  double stepPt = (maxPt-minPt)/nintPt;
223  pTintervalsv.push_back(minPt);
224  for (int k=1;k<nintPt+1;k++) {
225  double d=0;
226  if(useLogPt)d=pow(10,minPt+k*stepPt);
227  else d=minPt+k*stepPt;
228  pTintervalsv.push_back(d);
229  totSIMvpT.push_back(0);
230  totASSvpT.push_back(0);
231  totASS2vpT.push_back(0);
232  totRECvpT.push_back(0);
233  totloopvpT.push_back(0);
234  totmisidvpT.push_back(0);
235  }
236  pTintervals.push_back(pTintervalsv);
237  totSIMpT.push_back(totSIMvpT);
238  totASSpT.push_back(totASSvpT);
239  totASS2pT.push_back(totASS2vpT);
240  totRECpT.push_back(totRECvpT);
241  totlooppT.push_back(totloopvpT);
242  totmisidpT.push_back(totmisidvpT);
243 
244  for (int k=1;k<nintHit+1;k++) {
245  totSIMv_hit.push_back(0);
246  totASSv_hit.push_back(0);
247  totASS2v_hit.push_back(0);
248  totRECv_hit.push_back(0);
249  totloopv_hit.push_back(0);
250  totmisidv_hit.push_back(0);
251  }
252  totSIM_hit.push_back(totSIMv_hit);
253  totASS_hit.push_back(totASSv_hit);
254  totASS2_hit.push_back(totASS2v_hit);
255  totREC_hit.push_back(totRECv_hit);
256  totloop_hit.push_back(totloopv_hit);
257  totmisid_hit.push_back(totmisidv_hit);
258 
259  double stepPhi = (maxPhi-minPhi)/nintPhi;
260  phiintervalsv.push_back(minPhi);
261  for (int k=1;k<nintPhi+1;k++) {
262  double d=minPhi+k*stepPhi;
263  phiintervalsv.push_back(d);
264  totSIMv_phi.push_back(0);
265  totASSv_phi.push_back(0);
266  totASS2v_phi.push_back(0);
267  totRECv_phi.push_back(0);
268  totloopv_phi.push_back(0);
269  totmisidv_phi.push_back(0);
270  }
271  phiintervals.push_back(phiintervalsv);
272  totSIM_phi.push_back(totSIMv_phi);
273  totASS_phi.push_back(totASSv_phi);
274  totASS2_phi.push_back(totASS2v_phi);
275  totREC_phi.push_back(totRECv_phi);
276  totloop_phi.push_back(totloopv_phi);
277  totmisid_phi.push_back(totmisidv_phi);
278 
279  double stepDxy = (maxDxy-minDxy)/nintDxy;
280  dxyintervalsv.push_back(minDxy);
281  for (int k=1;k<nintDxy+1;k++) {
282  double d=minDxy+k*stepDxy;
283  dxyintervalsv.push_back(d);
284  totSIMv_dxy.push_back(0);
285  totASSv_dxy.push_back(0);
286  totASS2v_dxy.push_back(0);
287  totRECv_dxy.push_back(0);
288  totloopv_dxy.push_back(0);
289  totmisidv_dxy.push_back(0);
290  }
291  dxyintervals.push_back(dxyintervalsv);
292  totSIM_dxy.push_back(totSIMv_dxy);
293  totASS_dxy.push_back(totASSv_dxy);
294  totASS2_dxy.push_back(totASS2v_dxy);
295  totREC_dxy.push_back(totRECv_dxy);
296  totloop_dxy.push_back(totloopv_dxy);
297  totmisid_dxy.push_back(totmisidv_dxy);
298 
299 
300  double stepDz = (maxDz-minDz)/nintDz;
301  dzintervalsv.push_back(minDz);
302  for (int k=1;k<nintDz+1;k++) {
303  double d=minDz+k*stepDz;
304  dzintervalsv.push_back(d);
305  totSIMv_dz.push_back(0);
306  totASSv_dz.push_back(0);
307  totASS2v_dz.push_back(0);
308  totRECv_dz.push_back(0);
309  totloopv_dz.push_back(0);
310  totmisidv_dz.push_back(0);
311  }
312  dzintervals.push_back(dzintervalsv);
313  totSIM_dz.push_back(totSIMv_dz);
314  totASS_dz.push_back(totASSv_dz);
315  totASS2_dz.push_back(totASS2v_dz);
316  totREC_dz.push_back(totRECv_dz);
317  totloop_dz.push_back(totloopv_dz);
318  totmisid_dz.push_back(totmisidv_dz);
319 
320  double stepVertpos = (maxVertpos-minVertpos)/nintVertpos;
321  vertposintervalsv.push_back(minVertpos);
322  for (int k=1;k<nintVertpos+1;k++) {
323  double d=minVertpos+k*stepVertpos;
324  vertposintervalsv.push_back(d);
325  totSIMv_vertpos.push_back(0);
326  totASSv_vertpos.push_back(0);
327  }
328  vertposintervals.push_back(vertposintervalsv);
329  totSIM_vertpos.push_back(totSIMv_vertpos);
330  totASS_vertpos.push_back(totASSv_vertpos);
331 
332  double stepZpos = (maxZpos-minZpos)/nintZpos;
333  zposintervalsv.push_back(minZpos);
334  for (int k=1;k<nintZpos+1;k++) {
335  double d=minZpos+k*stepZpos;
336  zposintervalsv.push_back(d);
337  totSIMv_zpos.push_back(0);
338  totASSv_zpos.push_back(0);
339  }
340  zposintervals.push_back(zposintervalsv);
341  totSIM_zpos.push_back(totSIMv_zpos);
342  totCONzpos.push_back(totSIMv_zpos);
343  totASS_zpos.push_back(totASSv_zpos);
344  totSIM_vertz_entire.push_back(totSIMv_zpos);
345  totASS_vertz_entire.push_back(totASSv_zpos);
346  totSIM_vertz_barrel.push_back(totSIMv_zpos);
347  totASS_vertz_barrel.push_back(totASSv_zpos);
348  totSIM_vertz_fwdpos.push_back(totSIMv_zpos);
349  totASS_vertz_fwdpos.push_back(totASSv_zpos);
350  totSIM_vertz_fwdneg.push_back(totSIMv_zpos);
351  totASS_vertz_fwdneg.push_back(totASSv_zpos);
352 
353  double stepVertcount=(maxVertcount-minVertcount)/nintVertcount;
354  vertcountintervalsv.push_back(minVertcount);
355  for (int k=1;k<nintVertcount+1;k++) {
356  double d=minVertcount+k*stepVertcount;
357  vertcountintervalsv.push_back(d);
358  totSIMv_vertcount.push_back(0);
359  totASSv_vertcount.push_back(0);
360  totASS2v_vertcount.push_back(0);
361  totRECv_vertcount.push_back(0);
362  }
363  vertcountintervals.push_back(vertcountintervalsv);
364  totSIM_vertcount_entire.push_back(totSIMv_vertcount);
365  totCONvertcount.push_back(totSIMv_vertcount);
366  totASS_vertcount_entire.push_back(totASSv_vertcount);
367  totASS2_vertcount_entire.push_back(totASS2v_vertcount);
368  totASS2_vertcount_entire_signal.push_back(totASS2v_vertcount);
369  totREC_vertcount_entire.push_back(totRECv_vertcount);
370  totSIM_vertcount_barrel.push_back(totSIMv_vertcount);
371  totASS_vertcount_barrel.push_back(totASSv_vertcount);
372  totASS2_vertcount_barrel.push_back(totASS2v_vertcount);
373  totREC_vertcount_barrel.push_back(totRECv_vertcount);
374  totSIM_vertcount_fwdpos.push_back(totSIMv_vertcount);
375  totASS_vertcount_fwdpos.push_back(totASSv_vertcount);
376  totASS2_vertcount_fwdpos.push_back(totASS2v_vertcount);
377  totREC_vertcount_fwdpos.push_back(totRECv_vertcount);
378  totSIM_vertcount_fwdneg.push_back(totSIMv_vertcount);
379  totASS_vertcount_fwdneg.push_back(totASSv_vertcount);
380  totASS2_vertcount_fwdneg.push_back(totASS2v_vertcount);
381  totREC_vertcount_fwdneg.push_back(totRECv_vertcount);
382  totFOMT_vertcount.push_back(totASSv_vertcount);
383 
384  totASS2_itpu_vertcount_entire.push_back(totASS2v_vertcount);
385  totASS2_itpu_vertcount_entire_signal.push_back(totASS2v_vertcount);
386 
387  totASS2_ootpu_entire.push_back(totASS2v_vertcount);
388  totREC_ootpu_entire.push_back(totRECv_vertcount);
389  totASS2_ootpu_barrel.push_back(totASS2v_vertcount);
390  totREC_ootpu_barrel.push_back(totRECv_vertcount);
391  totASS2_ootpu_fwdpos.push_back(totASS2v_vertcount);
392  totREC_ootpu_fwdpos.push_back(totRECv_vertcount);
393  totASS2_ootpu_fwdneg.push_back(totASS2v_vertcount);
394  totREC_ootpu_fwdneg.push_back(totRECv_vertcount);
395 
396 }
397 
398 void MTVHistoProducerAlgoForTracker::bookSimHistos(){
399  h_ptSIM.push_back( dbe_->book1D("ptSIM", "generated p_{t}", 5500, 0, 110 ) );
400  h_etaSIM.push_back( dbe_->book1D("etaSIM", "generated pseudorapidity", 500, -2.5, 2.5 ) );
401  h_tracksSIM.push_back( dbe_->book1D("tracksSIM","number of simulated tracks",200,-0.5,99.5) );
402  h_vertposSIM.push_back( dbe_->book1D("vertposSIM","Transverse position of sim vertices",100,0.,120.) );
403  h_bunchxSIM.push_back( dbe_->book1D("bunchxSIM", "bunch crossing", 22, -5, 5 ) );
404 }
405 
406 
407 void MTVHistoProducerAlgoForTracker::bookRecoHistos(){
408  h_tracks.push_back( dbe_->book1D("tracks","number of reconstructed tracks",401,-0.5,400.5) );
409  h_fakes.push_back( dbe_->book1D("fakes","number of fake reco tracks",101,-0.5,100.5) );
410  h_charge.push_back( dbe_->book1D("charge","charge",3,-1.5,1.5) );
411 
412  h_hits.push_back( dbe_->book1D("hits", "number of hits per track", nintHit,minHit,maxHit ) );
413  h_losthits.push_back( dbe_->book1D("losthits", "number of lost hits per track", nintHit,minHit,maxHit) );
414  h_nchi2.push_back( dbe_->book1D("chi2", "normalized #chi^{2}", 200, 0, 20 ) );
415  h_nchi2_prob.push_back( dbe_->book1D("chi2_prob", "normalized #chi^{2} probability",100,0,1));
416 
417  h_algo.push_back( dbe_->book1D("h_algo","Tracks by algo",15,0.0,15.0) );
418 
420  h_recoeta.push_back( dbe_->book1D("num_reco_eta","N of reco track vs eta",nintEta,minEta,maxEta) );
421  h_assoceta.push_back( dbe_->book1D("num_assoc(simToReco)_eta","N of associated tracks (simToReco) vs eta",nintEta,minEta,maxEta) );
422  h_assoc2eta.push_back( dbe_->book1D("num_assoc(recoToSim)_eta","N of associated (recoToSim) tracks vs eta",nintEta,minEta,maxEta) );
423  h_simuleta.push_back( dbe_->book1D("num_simul_eta","N of simulated tracks vs eta",nintEta,minEta,maxEta) );
424  h_loopereta.push_back( dbe_->book1D("num_duplicate_eta","N of associated (recoToSim) duplicate tracks vs eta",nintEta,minEta,maxEta) );
425  h_misideta.push_back( dbe_->book1D("num_chargemisid_eta","N of associated (recoToSim) charge misIDed tracks vs eta",nintEta,minEta,maxEta) );
426  h_recopT.push_back( dbe_->book1D("num_reco_pT","N of reco track vs pT",nintPt,minPt,maxPt) );
427  h_assocpT.push_back( dbe_->book1D("num_assoc(simToReco)_pT","N of associated tracks (simToReco) vs pT",nintPt,minPt,maxPt) );
428  h_assoc2pT.push_back( dbe_->book1D("num_assoc(recoToSim)_pT","N of associated (recoToSim) tracks vs pT",nintPt,minPt,maxPt) );
429  h_simulpT.push_back( dbe_->book1D("num_simul_pT","N of simulated tracks vs pT",nintPt,minPt,maxPt) );
430  h_looperpT.push_back( dbe_->book1D("num_duplicate_pT","N of associated (recoToSim) duplicate tracks vs pT",nintPt,minPt,maxPt) );
431  h_misidpT.push_back( dbe_->book1D("num_chargemisid_pT","N of associated (recoToSim) charge misIDed tracks vs pT",nintPt,minPt,maxPt) );
432  //
433  h_recohit.push_back( dbe_->book1D("num_reco_hit","N of reco track vs hit",nintHit,minHit,maxHit) );
434  h_assochit.push_back( dbe_->book1D("num_assoc(simToReco)_hit","N of associated tracks (simToReco) vs hit",nintHit,minHit,maxHit) );
435  h_assoc2hit.push_back( dbe_->book1D("num_assoc(recoToSim)_hit","N of associated (recoToSim) tracks vs hit",nintHit,minHit,maxHit) );
436  h_simulhit.push_back( dbe_->book1D("num_simul_hit","N of simulated tracks vs hit",nintHit,minHit,maxHit) );
437  h_looperhit.push_back( dbe_->book1D("num_duplicate_hit","N of associated (recoToSim) duplicate tracks vs hit",nintHit,minHit,maxHit) );
438  h_misidhit.push_back( dbe_->book1D("num_chargemisid_hit","N of associated (recoToSim) charge misIDed tracks vs hit",nintHit,minHit,maxHit) );
439  //
440  h_recophi.push_back( dbe_->book1D("num_reco_phi","N of reco track vs phi",nintPhi,minPhi,maxPhi) );
441  h_assocphi.push_back( dbe_->book1D("num_assoc(simToReco)_phi","N of associated tracks (simToReco) vs phi",nintPhi,minPhi,maxPhi) );
442  h_assoc2phi.push_back( dbe_->book1D("num_assoc(recoToSim)_phi","N of associated (recoToSim) tracks vs phi",nintPhi,minPhi,maxPhi) );
443  h_simulphi.push_back( dbe_->book1D("num_simul_phi","N of simulated tracks vs phi",nintPhi,minPhi,maxPhi) );
444  h_looperphi.push_back( dbe_->book1D("num_duplicate_phi","N of associated (recoToSim) duplicate tracks vs phi",nintPhi,minPhi,maxPhi) );
445  h_misidphi.push_back( dbe_->book1D("num_chargemisid_phi","N of associated (recoToSim) charge misIDed tracks vs phi",nintPhi,minPhi,maxPhi) );
446 
447  h_recodxy.push_back( dbe_->book1D("num_reco_dxy","N of reco track vs dxy",nintDxy,minDxy,maxDxy) );
448  h_assocdxy.push_back( dbe_->book1D("num_assoc(simToReco)_dxy","N of associated tracks (simToReco) vs dxy",nintDxy,minDxy,maxDxy) );
449  h_assoc2dxy.push_back( dbe_->book1D("num_assoc(recoToSim)_dxy","N of associated (recoToSim) tracks vs dxy",nintDxy,minDxy,maxDxy) );
450  h_simuldxy.push_back( dbe_->book1D("num_simul_dxy","N of simulated tracks vs dxy",nintDxy,minDxy,maxDxy) );
451  h_looperdxy.push_back( dbe_->book1D("num_duplicate_dxy","N of associated (recoToSim) looper tracks vs dxy",nintDxy,minDxy,maxDxy) );
452  h_misiddxy.push_back( dbe_->book1D("num_chargemisid_dxy","N of associated (recoToSim) charge misIDed tracks vs dxy",nintDxy,minDxy,maxDxy) );
453 
454  h_recodz.push_back( dbe_->book1D("num_reco_dz","N of reco track vs dz",nintDz,minDz,maxDz) );
455  h_assocdz.push_back( dbe_->book1D("num_assoc(simToReco)_dz","N of associated tracks (simToReco) vs dz",nintDz,minDz,maxDz) );
456  h_assoc2dz.push_back( dbe_->book1D("num_assoc(recoToSim)_dz","N of associated (recoToSim) tracks vs dz",nintDz,minDz,maxDz) );
457  h_simuldz.push_back( dbe_->book1D("num_simul_dz","N of simulated tracks vs dz",nintDz,minDz,maxDz) );
458  h_looperdz.push_back( dbe_->book1D("num_duplicate_dz","N of associated (recoToSim) looper tracks vs dz",nintDz,minDz,maxDz) );
459  h_misiddz.push_back( dbe_->book1D("num_chargemisid_versus_dz","N of associated (recoToSim) charge misIDed tracks vs dz",nintDz,minDz,maxDz) );
460 
461  h_assocvertpos.push_back( dbe_->book1D("num_assoc(simToReco)_vertpos",
462  "N of associated tracks (simToReco) vs transverse vert position",
463  nintVertpos,minVertpos,maxVertpos) );
464  h_simulvertpos.push_back( dbe_->book1D("num_simul_vertpos","N of simulated tracks vs transverse vert position",
465  nintVertpos,minVertpos,maxVertpos) );
466 
467  h_assoczpos.push_back( dbe_->book1D("num_assoc(simToReco)_zpos","N of associated tracks (simToReco) vs z vert position",
468  nintZpos,minZpos,maxZpos) );
469  h_simulzpos.push_back( dbe_->book1D("num_simul_zpos","N of simulated tracks vs z vert position",nintZpos,minZpos,maxZpos) );
470 
471 
472  h_reco_vertcount_entire.push_back( dbe_->book1D("num_reco_vertcount_entire","N of reco tracks vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
473  h_assoc_vertcount_entire.push_back( dbe_->book1D("num_assoc(simToReco)_vertcount_entire","N of associated tracks (simToReco) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
474  h_assoc2_vertcount_entire.push_back( dbe_->book1D("num_assoc(recoToSim)_vertcount_entire","N of associated (recoToSim) tracks vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
475  h_simul_vertcount_entire.push_back( dbe_->book1D("num_simul_vertcount_entire","N of simulated tracks vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
476 
477  h_reco_vertcount_barrel.push_back( dbe_->book1D("num_reco_vertcount_barrel","N of reco tracks in barrel vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
478  h_assoc_vertcount_barrel.push_back( dbe_->book1D("num_assoc(simToReco)_vertcount_barrel","N of associated tracks (simToReco) in barrel vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
479  h_assoc2_vertcount_barrel.push_back( dbe_->book1D("num_assoc(recoToSim)_vertcount_barrel","N of associated (recoToSim) tracks in barrel vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
480  h_simul_vertcount_barrel.push_back( dbe_->book1D("num_simul_vertcount_barrel","N of simulated tracks in barrel vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
481 
482  h_reco_vertcount_fwdpos.push_back( dbe_->book1D("num_reco_vertcount_fwdpos","N of reco tracks in endcap(+) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
483  h_assoc_vertcount_fwdpos.push_back( dbe_->book1D("num_assoc(simToReco)_vertcount_fwdpos","N of associated tracks (simToReco) in endcap(+) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
484  h_assoc2_vertcount_fwdpos.push_back( dbe_->book1D("num_assoc(recoToSim)_vertcount_fwdpos","N of associated (recoToSim) tracks in endcap(+) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
485  h_simul_vertcount_fwdpos.push_back( dbe_->book1D("num_simul_vertcount_fwdpos","N of simulated tracks in endcap(+) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
486 
487  h_reco_vertcount_fwdneg.push_back( dbe_->book1D("num_reco_vertcount_fwdneg","N of reco tracks in endcap(-) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
488  h_assoc_vertcount_fwdneg.push_back( dbe_->book1D("num_assoc(simToReco)_vertcount_fwdneg","N of associated tracks (simToReco) in endcap(-) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
489  h_assoc2_vertcount_fwdneg.push_back( dbe_->book1D("num_assoc(recoToSim)_vertcount_fwdneg","N of associated (recoToSim) tracks in endcap(-) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
490  h_simul_vertcount_fwdneg.push_back( dbe_->book1D("num_simul_vertcount_fwdneg","N of simulated tracks in endcap(-) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
491 
492  h_assoc_vertz_entire.push_back( dbe_->book1D("num_assoc(simToReco)_vertz_entire","N of associated tracks (simToReco) in entire vs z of primary intercation vertex",nintZpos,minZpos,maxZpos) );
493  h_simul_vertz_entire.push_back( dbe_->book1D("num_simul_vertz_entire","N of simulated tracks in entire vs N of pileup vertices",nintZpos,minZpos,maxZpos) );
494 
495  h_assoc_vertz_barrel.push_back( dbe_->book1D("num_assoc(simToReco)_vertz_barrel","N of associated tracks (simToReco) in barrel vs z of primary intercation vertex",nintZpos,minZpos,maxZpos) );
496  h_simul_vertz_barrel.push_back( dbe_->book1D("num_simul_vertz_barrel","N of simulated tracks in barrel vs N of pileup vertices",nintZpos,minZpos,maxZpos) );
497 
498  h_assoc_vertz_fwdpos.push_back( dbe_->book1D("num_assoc(simToReco)_vertz_fwdpos","N of associated tracks (simToReco) in endcap(+) vs z of primary interaction vertex",nintZpos,minZpos,maxZpos) );
499  h_simul_vertz_fwdpos.push_back( dbe_->book1D("num_simul_vertz_fwdpos","N of simulated tracks in endcap(+) vs z of primary interaction vertex",nintZpos,minZpos,maxZpos) );
500 
501  h_assoc_vertz_fwdneg.push_back( dbe_->book1D("num_assoc(simToReco)_vertz_fwdneg","N of associated tracks (simToReco) in endcap(-) vs N of pileup vertices",nintZpos,minZpos,maxZpos) );
502  h_simul_vertz_fwdneg.push_back( dbe_->book1D("num_simul_vertz_fwdneg","N of simulated tracks in endcap(-) vs z of primary interaction vertex",nintZpos,minZpos,maxZpos) );
503 /*
504  h_assoc2_itpu_eta.push_back( dbe_->book1D("num_assoc(recoToSim)_itpu_eta_entire","N of associated tracks (simToReco) from in time pileup vs eta",nintEta,minEta,maxEta) );
505  h_assoc2_itpu_sig_eta.push_back( dbe_->book1D("num_assoc(recoToSim)_itpu_eta_entire_signal","N of associated tracks (simToReco) from in time pileup vs eta",nintEta,minEta,maxEta) );
506  h_assoc2_itpu_vertcount.push_back( dbe_->book1D("num_assoc(recoToSim)_itpu_vertcount_entire","N of associated tracks (simToReco) from in time pileup vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
507  h_assoc2_itpu_sig_vertcount.push_back( dbe_->book1D("num_assoc(recoToSim)_itpu_vertcount_entire_signal","N of associated tracks (simToReco) from in time pileup vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );*/
508 
509  h_reco_ootpu_eta.push_back( dbe_->book1D("num_reco_ootpu_eta_entire","N of reco tracks vs eta",nintEta,minEta,maxEta) );
510  h_reco_ootpu_vertcount.push_back( dbe_->book1D("num_reco_ootpu_vertcount_entire","N of reco tracks vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
511 
512  h_assoc2_ootpu_entire.push_back( dbe_->book1D("num_assoc(recoToSim)_ootpu_eta_entire","N of associated tracks (simToReco) from out of time pileup vs eta",nintEta,minEta,maxEta) );
513  h_assoc2_ootpu_vertcount.push_back( dbe_->book1D("num_assoc(recoToSim)_ootpu_vertcount_entire","N of associated tracks (simToReco) from out of time pileup vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
514 
515  h_reco_ootpu_entire.push_back( dbe_->book1D("num_reco_ootpu_entire","N of reco tracks vs z of primary interaction vertex",nintVertcount,minVertcount,maxVertcount) );
516  h_assoc2_ootpu_barrel.push_back( dbe_->book1D("num_assoc(recoToSim)_ootpu_barrel","N of associated tracks (simToReco) from out of time pileup in barrel vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
517  h_reco_ootpu_barrel.push_back( dbe_->book1D("num_reco_ootpu_barrel","N of reco tracks in barrel vs z of primary interaction vertex",nintVertcount,minVertcount,maxVertcount) );
518  h_assoc2_ootpu_fwdpos.push_back( dbe_->book1D("num_assoc(recoToSim)_ootpu_fwdpos","N of associated tracks (simToReco) from out of time pileup in endcap(+) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
519  h_reco_ootpu_fwdpos.push_back( dbe_->book1D("num_reco_ootpu_fwdpos","N of reco tracks in endcap(+) vs z of primary interaction vertex",nintVertcount,minVertcount,maxVertcount) );
520  h_assoc2_ootpu_fwdneg.push_back( dbe_->book1D("num_assoc(recoToSim)_ootpu_fwdneg","N of associated tracks (simToReco) from out of time pileup in endcap(-) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
521  h_reco_ootpu_fwdneg.push_back( dbe_->book1D("num_reco_ootpu_fwdneg","N of reco tracks in endcap(-) vs z of primary interaction vertex",nintVertcount,minVertcount,maxVertcount) );
522 
523 
525 
526  h_eta.push_back( dbe_->book1D("eta", "pseudorapidity residue", 1000, -0.1, 0.1 ) );
527  h_pt.push_back( dbe_->book1D("pullPt", "pull of p_{t}", 100, -10, 10 ) );
528  h_pullTheta.push_back( dbe_->book1D("pullTheta","pull of #theta parameter",250,-25,25) );
529  h_pullPhi.push_back( dbe_->book1D("pullPhi","pull of #phi parameter",250,-25,25) );
530  h_pullDxy.push_back( dbe_->book1D("pullDxy","pull of dxy parameter",250,-25,25) );
531  h_pullDz.push_back( dbe_->book1D("pullDz","pull of dz parameter",250,-25,25) );
532  h_pullQoverp.push_back( dbe_->book1D("pullQoverp","pull of qoverp parameter",250,-25,25) );
533 
534  /* TO BE FIXED -----------
535  if (associators[ww]=="TrackAssociatorByChi2"){
536  h_assochi2.push_back( dbe_->book1D("assocChi2","track association #chi^{2}",1000000,0,100000) );
537  h_assochi2_prob.push_back(dbe_->book1D("assocChi2_prob","probability of association #chi^{2}",100,0,1));
538  } else if (associators[ww]=="quickTrackAssociatorByHits"){
539  h_assocFraction.push_back( dbe_->book1D("assocFraction","fraction of shared hits",200,0,2) );
540  h_assocSharedHit.push_back(dbe_->book1D("assocSharedHit","number of shared hits",20,0,20));
541  }
542  */
543  h_assocFraction.push_back( dbe_->book1D("assocFraction","fraction of shared hits",200,0,2) );
544  h_assocSharedHit.push_back(dbe_->book1D("assocSharedHit","number of shared hits",41,-0.5,40.5));
545  // ----------------------
546 
547  chi2_vs_nhits.push_back( dbe_->book2D("chi2_vs_nhits","#chi^{2} vs nhits",25,0,25,100,0,10) );
548 
549  etares_vs_eta.push_back( dbe_->book2D("etares_vs_eta","etaresidue vs eta",nintEta,minEta,maxEta,200,-0.1,0.1) );
550  nrec_vs_nsim.push_back( dbe_->book2D("nrec_vs_nsim","nrec vs nsim",401,-0.5,400.5,401,-0.5,400.5) );
551 
552  chi2_vs_eta.push_back( dbe_->book2D("chi2_vs_eta","chi2_vs_eta",nintEta,minEta,maxEta, 200, 0, 20 ));
553  chi2_vs_phi.push_back( dbe_->book2D("chi2_vs_phi","#chi^{2} vs #phi",nintPhi,minPhi,maxPhi, 200, 0, 20 ) );
554 
555  nhits_vs_eta.push_back( dbe_->book2D("nhits_vs_eta","nhits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
556  nPXBhits_vs_eta.push_back( dbe_->book2D("nPXBhits_vs_eta","# PXB its vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
557  nPXFhits_vs_eta.push_back( dbe_->book2D("nPXFhits_vs_eta","# PXF hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
558  nTIBhits_vs_eta.push_back( dbe_->book2D("nTIBhits_vs_eta","# TIB hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
559  nTIDhits_vs_eta.push_back( dbe_->book2D("nTIDhits_vs_eta","# TID hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
560  nTOBhits_vs_eta.push_back( dbe_->book2D("nTOBhits_vs_eta","# TOB hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
561  nTEChits_vs_eta.push_back( dbe_->book2D("nTEChits_vs_eta","# TEC hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
562 
563  nLayersWithMeas_vs_eta.push_back( dbe_->book2D("nLayersWithMeas_vs_eta","# Layers with measurement vs eta",
564  nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
565  nPXLlayersWithMeas_vs_eta.push_back( dbe_->book2D("nPXLlayersWithMeas_vs_eta","# PXL Layers with measurement vs eta",
566  nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
567  nSTRIPlayersWithMeas_vs_eta.push_back( dbe_->book2D("nSTRIPlayersWithMeas_vs_eta","# STRIP Layers with measurement vs eta",
568  nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
569  nSTRIPlayersWith1dMeas_vs_eta.push_back( dbe_->book2D("nSTRIPlayersWith1dMeas_vs_eta","# STRIP Layers with 1D measurement vs eta",
570  nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
571  nSTRIPlayersWith2dMeas_vs_eta.push_back( dbe_->book2D("nSTRIPlayersWith2dMeas_vs_eta","# STRIP Layers with 2D measurement vs eta",
572  nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
573 
574  nhits_vs_phi.push_back( dbe_->book2D("nhits_vs_phi","#hits vs #phi",nintPhi,minPhi,maxPhi,nintHit,minHit,maxHit) );
575 
576  nlosthits_vs_eta.push_back( dbe_->book2D("nlosthits_vs_eta","nlosthits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
577 
578  //resolution of track parameters
579  // dPt/Pt cotTheta Phi TIP LIP
580  // log10(pt)<0.5 100,0.1 240,0.08 100,0.015 100,0.1000 150,0.3000
581  // 0.5<log10(pt)<1.5 100,0.1 120,0.01 100,0.003 100,0.0100 150,0.0500
582  // >1.5 100,0.3 100,0.005 100,0.0008 100,0.0060 120,0.0300
583 
584  ptres_vs_eta.push_back(dbe_->book2D("ptres_vs_eta","ptres_vs_eta",
585  nintEta,minEta,maxEta, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
586 
587  ptres_vs_phi.push_back( dbe_->book2D("ptres_vs_phi","p_{t} res vs #phi",
588  nintPhi,minPhi,maxPhi, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
589 
590  ptres_vs_pt.push_back(dbe_->book2D("ptres_vs_pt","ptres_vs_pt",nintPt,minPt,maxPt, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
591 
592  cotThetares_vs_eta.push_back(dbe_->book2D("cotThetares_vs_eta","cotThetares_vs_eta",
593  nintEta,minEta,maxEta,cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
594 
595 
596  cotThetares_vs_pt.push_back(dbe_->book2D("cotThetares_vs_pt","cotThetares_vs_pt",
597  nintPt,minPt,maxPt, cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
598 
599 
600  phires_vs_eta.push_back(dbe_->book2D("phires_vs_eta","phires_vs_eta",
601  nintEta,minEta,maxEta, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
602 
603  phires_vs_pt.push_back(dbe_->book2D("phires_vs_pt","phires_vs_pt",
604  nintPt,minPt,maxPt, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
605 
606  phires_vs_phi.push_back(dbe_->book2D("phires_vs_phi","#phi res vs #phi",
607  nintPhi,minPhi,maxPhi,phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
608 
609  dxyres_vs_eta.push_back(dbe_->book2D("dxyres_vs_eta","dxyres_vs_eta",
610  nintEta,minEta,maxEta,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
611 
612  dxyres_vs_pt.push_back( dbe_->book2D("dxyres_vs_pt","dxyres_vs_pt",
613  nintPt,minPt,maxPt,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
614 
615  dzres_vs_eta.push_back(dbe_->book2D("dzres_vs_eta","dzres_vs_eta",
616  nintEta,minEta,maxEta,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
617 
618  dzres_vs_pt.push_back(dbe_->book2D("dzres_vs_pt","dzres_vs_pt",nintPt,minPt,maxPt,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
619 
620  ptmean_vs_eta_phi.push_back(dbe_->bookProfile2D("ptmean_vs_eta_phi","mean p_{t} vs #eta and #phi",
621  nintPhi,minPhi,maxPhi,nintEta,minEta,maxEta,1000,0,1000));
622  phimean_vs_eta_phi.push_back(dbe_->bookProfile2D("phimean_vs_eta_phi","mean #phi vs #eta and #phi",
623  nintPhi,minPhi,maxPhi,nintEta,minEta,maxEta,nintPhi,minPhi,maxPhi));
624 
625  //pulls of track params vs eta: to be used with fitslicesytool
626  dxypull_vs_eta.push_back(dbe_->book2D("dxypull_vs_eta","dxypull_vs_eta",nintEta,minEta,maxEta,100,-10,10));
627  ptpull_vs_eta.push_back(dbe_->book2D("ptpull_vs_eta","ptpull_vs_eta",nintEta,minEta,maxEta,100,-10,10));
628  dzpull_vs_eta.push_back(dbe_->book2D("dzpull_vs_eta","dzpull_vs_eta",nintEta,minEta,maxEta,100,-10,10));
629  phipull_vs_eta.push_back(dbe_->book2D("phipull_vs_eta","phipull_vs_eta",nintEta,minEta,maxEta,100,-10,10));
630  thetapull_vs_eta.push_back(dbe_->book2D("thetapull_vs_eta","thetapull_vs_eta",nintEta,minEta,maxEta,100,-10,10));
631 
632  // h_ptshiftetamean.push_back( dbe_->book1D("h_ptshifteta_Mean","<#deltapT/pT>[%] vs #eta",nintEta,minEta,maxEta) );
633 
634 
635  //pulls of track params vs phi
636  ptpull_vs_phi.push_back(dbe_->book2D("ptpull_vs_phi","p_{t} pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
637  phipull_vs_phi.push_back(dbe_->book2D("phipull_vs_phi","#phi pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
638  thetapull_vs_phi.push_back(dbe_->book2D("thetapull_vs_phi","#theta pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));
639 
640 
641  nrecHit_vs_nsimHit_sim2rec.push_back( dbe_->book2D("nrecHit_vs_nsimHit_sim2rec","nrecHit vs nsimHit (Sim2RecAssoc)",
642  nintHit,minHit,maxHit, nintHit,minHit,maxHit ));
643  nrecHit_vs_nsimHit_rec2sim.push_back( dbe_->book2D("nrecHit_vs_nsimHit_rec2sim","nrecHit vs nsimHit (Rec2simAssoc)",
644  nintHit,minHit,maxHit, nintHit,minHit,maxHit ));
645 
646  // dE/dx stuff
647  // FIXME: it would be nice to have an array
648  h_dedx_estim1.push_back( dbe_->book1D("h_dedx_estim1","dE/dx estimator 1",nintDeDx,minDeDx,maxDeDx) );
649  h_dedx_estim2.push_back( dbe_->book1D("h_dedx_estim2","dE/dx estimator 2",nintDeDx,minDeDx,maxDeDx) );
650  h_dedx_nom1.push_back( dbe_->book1D("h_dedx_nom1","dE/dx number of measurements",nintHit,minHit,maxHit) );
651  h_dedx_nom2.push_back( dbe_->book1D("h_dedx_nom2","dE/dx number of measurements",nintHit,minHit,maxHit) );
652  h_dedx_sat1.push_back( dbe_->book1D("h_dedx_sat1","dE/dx number of measurements with saturation",nintHit,minHit,maxHit) );
653  h_dedx_sat2.push_back( dbe_->book1D("h_dedx_sat2","dE/dx number of measurements with saturation",nintHit,minHit,maxHit) );
654 
655  // PU special stuff
656  h_con_eta.push_back( dbe_->book1D("num_con_eta","N of PU tracks vs eta",nintEta,minEta,maxEta) );
657  h_con_vertcount.push_back( dbe_->book1D("num_con_vertcount","N of PU tracks vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
658  h_con_zpos.push_back( dbe_->book1D("num_con_zpos","N of PU tracks vs z of primary interaction vertex",nintZpos,minZpos,maxZpos) );
659 
660 
661  if(useLogPt){
662  BinLogX(dzres_vs_pt.back()->getTH2F());
663  BinLogX(dxyres_vs_pt.back()->getTH2F());
664  BinLogX(phires_vs_pt.back()->getTH2F());
665  BinLogX(cotThetares_vs_pt.back()->getTH2F());
666  BinLogX(ptres_vs_pt.back()->getTH2F());
667  BinLogX(h_looperpT.back()->getTH1F());
668  BinLogX(h_misidpT.back()->getTH1F());
669  BinLogX(h_recopT.back()->getTH1F());
670  BinLogX(h_assocpT.back()->getTH1F());
671  BinLogX(h_assoc2pT.back()->getTH1F());
672  BinLogX(h_simulpT.back()->getTH1F());
673  }
674 }
675 
676 void MTVHistoProducerAlgoForTracker::bookRecoHistosForStandaloneRunning(){
677  h_effic.push_back( dbe_->book1D("effic","efficiency vs #eta",nintEta,minEta,maxEta) );
678  h_efficPt.push_back( dbe_->book1D("efficPt","efficiency vs pT",nintPt,minPt,maxPt) );
679  h_effic_vs_hit.push_back( dbe_->book1D("effic_vs_hit","effic vs hit",nintHit,minHit,maxHit) );
680  h_effic_vs_phi.push_back( dbe_->book1D("effic_vs_phi","effic vs phi",nintPhi,minPhi,maxPhi) );
681  h_effic_vs_dxy.push_back( dbe_->book1D("effic_vs_dxy","effic vs dxy",nintDxy,minDxy,maxDxy) );
682  h_effic_vs_dz.push_back( dbe_->book1D("effic_vs_dz","effic vs dz",nintDz,minDz,maxDz) );
683  h_effic_vs_vertpos.push_back( dbe_->book1D("effic_vs_vertpos","effic vs vertpos",nintVertpos,minVertpos,maxVertpos) );
684  h_effic_vs_zpos.push_back( dbe_->book1D("effic_vs_zpos","effic vs zpos",nintZpos,minZpos,maxZpos) );
685  h_effic_vertcount_entire.push_back( dbe_->book1D("effic_vertcount_entire","efficiency vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
686  h_effic_vertcount_barrel.push_back( dbe_->book1D("effic_vertcount_barrel","efficiency in barrel vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
687  h_effic_vertcount_fwdpos.push_back( dbe_->book1D("effic_vertcount_fwdpos","efficiency in endcap(+) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
688  h_effic_vertcount_fwdneg.push_back( dbe_->book1D("effic_vertcount_fwdneg","efficiency in endcap(-) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
689  h_effic_vertz_entire.push_back( dbe_->book1D("effic_vertz_entire","efficiency vs z of primary interaction vertex",nintZpos,minZpos,maxZpos) );
690  h_effic_vertz_barrel.push_back( dbe_->book1D("effic_vertz_barrel","efficiency in barrel vs z of primary interaction vertex",nintZpos,minZpos,maxZpos) );
691  h_effic_vertz_fwdpos.push_back( dbe_->book1D("effic_vertz_fwdpos","efficiency in endcap(+) vs z of primary interaction vertex",nintZpos,minZpos,maxZpos) );
692  h_effic_vertz_fwdneg.push_back( dbe_->book1D("effic_vertz_fwdneg","efficiency in endcap(-) vs z of primary interaction vertex",nintZpos,minZpos,maxZpos) );
693 
694  h_looprate.push_back( dbe_->book1D("duplicatesRate","loop rate vs #eta",nintEta,minEta,maxEta) );
695  h_misidrate.push_back( dbe_->book1D("chargeMisIdRate","misid rate vs #eta",nintEta,minEta,maxEta) );
696  h_fakerate.push_back( dbe_->book1D("fakerate","fake rate vs #eta",nintEta,minEta,maxEta) );
697  h_loopratepT.push_back( dbe_->book1D("duplicatesRate_Pt","loop rate vs pT",nintPt,minPt,maxPt) );
698  h_misidratepT.push_back( dbe_->book1D("chargeMisIdRate_Pt","misid rate vs pT",nintPt,minPt,maxPt) );
699  h_fakeratePt.push_back( dbe_->book1D("fakeratePt","fake rate vs pT",nintPt,minPt,maxPt) );
700  h_loopratehit.push_back( dbe_->book1D("duplicatesRate_hit","loop rate vs hit",nintHit,minHit,maxHit) );
701  h_misidratehit.push_back( dbe_->book1D("chargeMisIdRate_hit","misid rate vs hit",nintHit,minHit,maxHit) );
702  h_fake_vs_hit.push_back( dbe_->book1D("fakerate_vs_hit","fake rate vs hit",nintHit,minHit,maxHit) );
703  h_loopratephi.push_back( dbe_->book1D("duplicatesRate_phi","loop rate vs #phi",nintPhi,minPhi,maxPhi) );
704  h_misidratephi.push_back( dbe_->book1D("chargeMisIdRate_phi","misid rate vs #phi",nintPhi,minPhi,maxPhi) );
705  h_fake_vs_phi.push_back( dbe_->book1D("fakerate_vs_phi","fake vs #phi",nintPhi,minPhi,maxPhi) );
706  h_loopratedxy.push_back( dbe_->book1D("duplicatesRate_dxy","loop rate vs dxy",nintDxy,minDxy,maxDxy) );
707  h_misidratedxy.push_back( dbe_->book1D("chargeMisIdRate_dxy","misid rate vs dxy",nintDxy,minDxy,maxDxy) );
708  h_fake_vs_dxy.push_back( dbe_->book1D("fakerate_vs_dxy","fake rate vs dxy",nintDxy,minDxy,maxDxy) );
709  h_loopratedz.push_back( dbe_->book1D("duplicatesRate_dz","loop rate vs dz",nintDz,minDz,maxDz) );
710  h_misidratedz.push_back( dbe_->book1D("chargeMisIdRate_dz","misid rate vs dz",nintDz,minDz,maxDz) );
711  h_fake_vs_dz.push_back( dbe_->book1D("fakerate_vs_dz","fake vs dz",nintDz,minDz,maxDz) );
712  h_fakerate_vertcount_entire.push_back( dbe_->book1D("fakerate_vertcount_entire","fake rate vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
713  h_fakerate_vertcount_barrel.push_back( dbe_->book1D("fakerate_vertcount_barrel","fake rate in barrel vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
714  h_fakerate_vertcount_fwdpos.push_back( dbe_->book1D("fakerate_vertcount_fwdpos","fake rate in endcap(+) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
715  h_fakerate_vertcount_fwdneg.push_back( dbe_->book1D("fakerate_vertcount_fwdneg","fake rate in endcap(-) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
716  h_fakerate_ootpu_entire.push_back( dbe_->book1D("fakerate_ootpu_eta_entire","Out of time pileup fake rate vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
717  h_fakerate_ootpu_barrel.push_back( dbe_->book1D("fakerate_ootpu_barrel","Out of time pileup fake rate in barrel vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
718  h_fakerate_ootpu_fwdpos.push_back( dbe_->book1D("fakerate_ootpu_fwdpos","Out of time pileup fake rate in endcap(+) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
719  h_fakerate_ootpu_fwdneg.push_back( dbe_->book1D("fakerate_ootpu_fwdneg","Out of time pileup fake rate in endcap(-) vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
720 
721  h_fomt_eta.push_back( dbe_->book1D("fomt","fraction of misreconstructed tracks vs #eta",nintEta,minEta,maxEta) );
722  h_fomt_sig_eta.push_back( dbe_->book1D("fomt_signal","fraction of misreconstructed tracks vs #eta",nintEta,minEta,maxEta) );
723  h_fomt_vertcount.push_back( dbe_->book1D("fomt_vertcount","fraction of misreconstructed tracks vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
724  h_fomt_sig_vertcount.push_back( dbe_->book1D("fomt_vertcount_signal","fraction of misreconstructed tracks vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
725  h_fomt_ootpu_eta.push_back( dbe_->book1D("fomt_ootpu_eta","Out of time pileup fraction of misreconstructed tracks vs #eta",nintEta,minEta,maxEta) );
726  h_fomt_ootpu_vertcount.push_back( dbe_->book1D("fomt_ootpu_vertcount","Out of time pileup fraction of misreconstructed tracks vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
727  h_fomt_itpu_eta.push_back( dbe_->book1D("fomt_itpu_eta","In time pileup fraction of misreconstructed tracks vs eta",nintEta,minEta,maxEta) );
728  h_fomt_itpu_vertcount.push_back( dbe_->book1D("fomt_itpu_vertcount","In time pileup fraction of misreconstructed tracks vs N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
729 
730  h_effic_PU_eta.push_back( dbe_->book1D("effic_PU_eta","PU efficiency vs #eta",nintEta,minEta,maxEta) );
731  h_effic_PU_vertcount.push_back( dbe_->book1D("effic_PU_vertcount","PU efficiency s N of pileup vertices",nintVertcount,minVertcount,maxVertcount) );
732  h_effic_PU_zpos.push_back( dbe_->book1D("effic_PU_zpos","PU efficiency vs z of primary interaction vertex",nintZpos,minZpos,maxZpos) );
733 
734  h_chi2meanhitsh.push_back( dbe_->bookProfile("chi2mean_vs_nhits","mean #chi^{2} vs nhits",25,0,25,100,0,10) );
735  h_chi2meanh.push_back( dbe_->bookProfile("chi2mean","mean #chi^{2} vs #eta",nintEta,minEta,maxEta, 200, 0, 20) );
736  h_chi2mean_vs_phi.push_back( dbe_->bookProfile("chi2mean_vs_phi","mean of #chi^{2} vs #phi",nintPhi,minPhi,maxPhi, 200, 0, 20) );
737 
738  h_hits_eta.push_back( dbe_->bookProfile("hits_eta","mean #hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
739  h_PXBhits_eta.push_back( dbe_->bookProfile("PXBhits_eta","mean # PXB hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
740  h_PXFhits_eta.push_back( dbe_->bookProfile("PXFhits_eta","mean # PXF hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
741  h_TIBhits_eta.push_back( dbe_->bookProfile("TIBhits_eta","mean # TIB hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
742  h_TIDhits_eta.push_back( dbe_->bookProfile("TIDhits_eta","mean # TID hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
743  h_TOBhits_eta.push_back( dbe_->bookProfile("TOBhits_eta","mean # TOB hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
744  h_TEChits_eta.push_back( dbe_->bookProfile("TEChits_eta","mean # TEC hits vs eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
745 
746  h_LayersWithMeas_eta.push_back(dbe_->bookProfile("LayersWithMeas_eta","mean # LayersWithMeas vs eta",
747  nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
748  h_PXLlayersWithMeas_eta.push_back(dbe_->bookProfile("PXLlayersWith2dMeas_eta","mean # PXLlayersWithMeas vs eta",
749  nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
750  h_STRIPlayersWithMeas_eta.push_back(dbe_->bookProfile("STRIPlayersWithMeas_eta","mean # STRIPlayersWithMeas vs eta",
751  nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
752  h_STRIPlayersWith1dMeas_eta.push_back(dbe_->bookProfile("STRIPlayersWith1dMeas_eta","mean # STRIPlayersWith1dMeas vs eta",
753  nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
754  h_STRIPlayersWith2dMeas_eta.push_back(dbe_->bookProfile("STRIPlayersWith2dMeas_eta","mean # STRIPlayersWith2dMeas vs eta",
755  nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
756  h_hits_phi.push_back( dbe_->bookProfile("hits_phi","mean #hits vs #phi",nintPhi,minPhi,maxPhi, nintHit,minHit,maxHit) );
757  h_losthits_eta.push_back( dbe_->bookProfile("losthits_eta","losthits_eta",nintEta,minEta,maxEta,nintHit,minHit,maxHit) );
758 
759  h_ptrmsh.push_back( dbe_->book1D("ptres_vs_eta_Sigma","#sigma(#deltap_{t}/p_{t}) vs #eta",nintEta,minEta,maxEta) );
760  h_ptmeanhPhi.push_back( dbe_->book1D("ptres_vs_phi_Mean","mean of p_{t} resolution vs #phi",nintPhi,minPhi,maxPhi));
761  h_ptrmshPhi.push_back( dbe_->book1D("ptres_vs_phi_Sigma","#sigma(#deltap_{t}/p_{t}) vs #phi",nintPhi,minPhi,maxPhi) );
762  h_ptmeanhPt.push_back( dbe_->book1D("ptres_vs_pt_Mean","mean of p_{t} resolution vs p_{t}",nintPt,minPt,maxPt));
763  h_ptrmshPt.push_back( dbe_->book1D("ptres_vs_pt_Sigma","#sigma(#deltap_{t}/p_{t}) vs pT",nintPt,minPt,maxPt) );
764  h_cotThetameanh.push_back( dbe_->book1D("cotThetares_vs_eta_Mean","#sigma(cot(#theta)) vs #eta Mean",nintEta,minEta,maxEta) );
765  h_cotThetarmsh.push_back( dbe_->book1D("cotThetares_vs_eta_Sigma","#sigma(cot(#theta)) vs #eta Sigma",nintEta,minEta,maxEta) );
766  h_cotThetameanhPt.push_back( dbe_->book1D("cotThetares_vs_pt_Mean","#sigma(cot(#theta)) vs pT Mean",nintPt,minPt,maxPt) );
767  h_cotThetarmshPt.push_back( dbe_->book1D("cotThetares_vs_pt_Sigma","#sigma(cot(#theta)) vs pT Sigma",nintPt,minPt,maxPt) );
768  h_phimeanh.push_back(dbe_->book1D("phires_vs_eta_Mean","mean of #phi res vs #eta",nintEta,minEta,maxEta));
769  h_phirmsh.push_back( dbe_->book1D("phires_vs_eta_Sigma","#sigma(#delta#phi) vs #eta",nintEta,minEta,maxEta) );
770  h_phimeanhPt.push_back(dbe_->book1D("phires_vs_pt_Mean","mean of #phi res vs pT",nintPt,minPt,maxPt));
771  h_phirmshPt.push_back( dbe_->book1D("phires_vs_pt_Sigma","#sigma(#delta#phi) vs pT",nintPt,minPt,maxPt) );
772  h_phimeanhPhi.push_back(dbe_->book1D("phires_vs_phi_Mean","mean of #phi res vs #phi",nintPhi,minPhi,maxPhi));
773  h_phirmshPhi.push_back( dbe_->book1D("phires_vs_phi_Sigma","#sigma(#delta#phi) vs #phi",nintPhi,minPhi,maxPhi) );
774  h_dxymeanh.push_back( dbe_->book1D("dxyres_vs_eta_Mean","mean of dxyres vs #eta",nintEta,minEta,maxEta) );
775  h_dxyrmsh.push_back( dbe_->book1D("dxyres_vs_eta_Sigma","#sigma(#deltadxy) vs #eta",nintEta,minEta,maxEta) );
776  h_dxymeanhPt.push_back( dbe_->book1D("dxyres_vs_pt_Mean","mean of dxyres vs pT",nintPt,minPt,maxPt) );
777  h_dxyrmshPt.push_back( dbe_->book1D("dxyres_vs_pt_Sigma","#sigmadxy vs pT",nintPt,minPt,maxPt) );
778  h_dzmeanh.push_back( dbe_->book1D("dzres_vs_eta_Mean","mean of dzres vs #eta",nintEta,minEta,maxEta) );
779  h_dzrmsh.push_back( dbe_->book1D("dzres_vs_eta_Sigma","#sigma(#deltadz) vs #eta",nintEta,minEta,maxEta) );
780  h_dzmeanhPt.push_back( dbe_->book1D("dzres_vs_pt_Mean","mean of dzres vs pT",nintPt,minPt,maxPt) );
781  h_dzrmshPt.push_back( dbe_->book1D("dzres_vs_pt_Sigma","#sigma(#deltadz vs pT",nintPt,minPt,maxPt) );
782  h_dxypulletamean.push_back( dbe_->book1D("h_dxypulleta_Mean","mean of dxy pull vs #eta",nintEta,minEta,maxEta) );
783  h_ptpulletamean.push_back( dbe_->book1D("h_ptpulleta_Mean","mean of p_{t} pull vs #eta",nintEta,minEta,maxEta) );
784  h_dzpulletamean.push_back( dbe_->book1D("h_dzpulleta_Mean","mean of dz pull vs #eta",nintEta,minEta,maxEta) );
785  h_phipulletamean.push_back( dbe_->book1D("h_phipulleta_Mean","mean of #phi pull vs #eta",nintEta,minEta,maxEta) );
786  h_thetapulletamean.push_back( dbe_->book1D("h_thetapulleta_Mean","mean of #theta pull vs #eta",nintEta,minEta,maxEta) );
787  h_dxypulleta.push_back( dbe_->book1D("h_dxypulleta_Sigma","#sigma of dxy pull vs #eta",nintEta,minEta,maxEta) );
788  h_ptpulleta.push_back( dbe_->book1D("h_ptpulleta_Sigma","#sigma of p_{t} pull vs #eta",nintEta,minEta,maxEta) );
789  h_dzpulleta.push_back( dbe_->book1D("h_dzpulleta_Sigma","#sigma of dz pull vs #eta",nintEta,minEta,maxEta) );
790  h_phipulleta.push_back( dbe_->book1D("h_phipulleta_Sigma","#sigma of #phi pull vs #eta",nintEta,minEta,maxEta) );
791  h_thetapulleta.push_back( dbe_->book1D("h_thetapulleta_Sigma","#sigma of #theta pull vs #eta",nintEta,minEta,maxEta) );
792  h_ptshifteta.push_back( dbe_->book1D("ptres_vs_eta_Mean","<#deltapT/pT>[%] vs #eta",nintEta,minEta,maxEta) );
793  h_ptpullphimean.push_back( dbe_->book1D("h_ptpullphi_Mean","mean of p_{t} pull vs #phi",nintPhi,minPhi,maxPhi) );
794  h_phipullphimean.push_back( dbe_->book1D("h_phipullphi_Mean","mean of #phi pull vs #phi",nintPhi,minPhi,maxPhi) );
795  h_thetapullphimean.push_back( dbe_->book1D("h_thetapullphi_Mean","mean of #theta pull vs #phi",nintPhi,minPhi,maxPhi) );
796  h_ptpullphi.push_back( dbe_->book1D("h_ptpullphi_Sigma","#sigma of p_{t} pull vs #phi",nintPhi,minPhi,maxPhi) );
797  h_phipullphi.push_back( dbe_->book1D("h_phipullphi_Sigma","#sigma of #phi pull vs #phi",nintPhi,minPhi,maxPhi) );
798  h_thetapullphi.push_back( dbe_->book1D("h_thetapullphi_Sigma","#sigma of #theta pull vs #phi",nintPhi,minPhi,maxPhi) );
799 
800  if(useLogPt){
801  BinLogX(h_dzmeanhPt.back()->getTH1F());
802  BinLogX(h_dzrmshPt.back()->getTH1F());
803  BinLogX(h_dxymeanhPt.back()->getTH1F());
804  BinLogX(h_dxyrmshPt.back()->getTH1F());
805  BinLogX(h_phimeanhPt.back()->getTH1F());
806  BinLogX(h_phirmshPt.back()->getTH1F());
807  BinLogX(h_cotThetameanhPt.back()->getTH1F());
808  BinLogX(h_cotThetarmshPt.back()->getTH1F());
809  BinLogX(h_ptmeanhPt.back()->getTH1F());
810  BinLogX(h_ptrmshPt.back()->getTH1F());
811  BinLogX(h_efficPt.back()->getTH1F());
812  BinLogX(h_fakeratePt.back()->getTH1F());
813  BinLogX(h_loopratepT.back()->getTH1F());
814  BinLogX(h_misidratepT.back()->getTH1F());
815  }
816 }
817 
818 void MTVHistoProducerAlgoForTracker::fill_generic_simTrack_histos(int count,
819  const TrackingParticle::Vector& momentumTP,
820  const TrackingParticle::Point& vertexTP,
821  int bx){
822  h_ptSIM[count]->Fill(sqrt(momentumTP.perp2()));
823  h_etaSIM[count]->Fill(momentumTP.eta());
824  h_vertposSIM[count]->Fill(sqrt(vertexTP.perp2()));
825  h_bunchxSIM[count]->Fill(bx);
826 }
827 
828 
829 
830 // TO BE FIXED USING PLAIN HISTOGRAMS INSTEAD OF RE-IMPLEMENTATION OF HISTOGRAMS (i.d. vectors<int/double>)
831 void MTVHistoProducerAlgoForTracker::fill_recoAssociated_simTrack_histos(int count,
832  const TrackingParticle& tp,
833  const TrackingParticle::Vector& momentumTP,
834  const TrackingParticle::Point& vertexTP,
835  double dxySim, double dzSim, int nSimHits,
836  const reco::Track* track,
837  int numVertices, double vertz){
838  bool isMatched = track;
839 
840  if((*TpSelectorForEfficiencyVsEta)(tp)){
841  //effic vs hits
842  int nSimHitsInBounds = std::min((int)nSimHits,int(maxHit-1));
843  totSIM_hit[count][nSimHitsInBounds]++;
844  if(isMatched) {
845  totASS_hit[count][nSimHitsInBounds]++;
846  nrecHit_vs_nsimHit_sim2rec[count]->Fill( track->numberOfValidHits(),nSimHits);
847  }
848 
849  //effic vs eta
850  for (unsigned int f=0; f<etaintervals[count].size()-1; f++){
851  if (getEta(momentumTP.eta())>etaintervals[count][f]&&
852  getEta(momentumTP.eta())<etaintervals[count][f+1]) {
853  totSIMeta[count][f]++;
854  if (isMatched) {
855  totASSeta[count][f]++;
856  }
857  }
858 
859  } // END for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
860 
861  //effic vs num pileup vertices
862  for (unsigned int f=0; f<vertcountintervals[count].size()-1; f++){
863  if (numVertices == vertcountintervals[count][f]) {
864  totSIM_vertcount_entire[count][f]++;
865  if (isMatched) {
866  totASS_vertcount_entire[count][f]++;
867  }
868  }
869  if (numVertices == vertcountintervals[count][f] && momentumTP.eta() <= 0.9 && momentumTP.eta() >= -0.9) {
870  totSIM_vertcount_barrel[count][f]++;
871  if (isMatched) {
872  totASS_vertcount_barrel[count][f]++;
873  }
874  }
875  if (numVertices == vertcountintervals[count][f] && momentumTP.eta() > 0.9) {
876  totSIM_vertcount_fwdpos[count][f]++;
877  if (isMatched) {
878  totASS_vertcount_fwdpos[count][f]++;
879  }
880  }
881  if (numVertices == vertcountintervals[count][f] && momentumTP.eta() < -0.9) {
882  totSIM_vertcount_fwdneg[count][f]++;
883  if (isMatched) {
884  totASS_vertcount_fwdneg[count][f]++;
885  }
886  }
887  }
888 
889  }
890 
891  if((*TpSelectorForEfficiencyVsPhi)(tp)){
892  for (unsigned int f=0; f<phiintervals[count].size()-1; f++){
893  if (momentumTP.phi() > phiintervals[count][f]&&
894  momentumTP.phi() <phiintervals[count][f+1]) {
895  totSIM_phi[count][f]++;
896  if (isMatched) {
897  totASS_phi[count][f]++;
898  }
899  }
900  } // END for (unsigned int f=0; f<phiintervals[count].size()-1; f++){
901  }
902 
903  if((*TpSelectorForEfficiencyVsPt)(tp)){
904  for (unsigned int f=0; f<pTintervals[count].size()-1; f++){
905  if (getPt(sqrt(momentumTP.perp2()))>pTintervals[count][f]&&
906  getPt(sqrt(momentumTP.perp2()))<pTintervals[count][f+1]) {
907  totSIMpT[count][f]++;
908  if (isMatched) {
909  totASSpT[count][f]++;
910  }
911  }
912  } // END for (unsigned int f=0; f<pTintervals[count].size()-1; f++){
913  }
914 
915  if((*TpSelectorForEfficiencyVsVTXR)(tp)){
916  for (unsigned int f=0; f<dxyintervals[count].size()-1; f++){
917  if (dxySim>dxyintervals[count][f]&&
918  dxySim<dxyintervals[count][f+1]) {
919  totSIM_dxy[count][f]++;
920  if (isMatched) {
921  totASS_dxy[count][f]++;
922  }
923  }
924  } // END for (unsigned int f=0; f<dxyintervals[count].size()-1; f++){
925 
926  for (unsigned int f=0; f<vertposintervals[count].size()-1; f++){
927  if (sqrt(vertexTP.perp2())>vertposintervals[count][f]&&
928  sqrt(vertexTP.perp2())<vertposintervals[count][f+1]) {
929  totSIM_vertpos[count][f]++;
930  if (isMatched) {
931  totASS_vertpos[count][f]++;
932  }
933  }
934  } // END for (unsigned int f=0; f<vertposintervals[count].size()-1; f++){
935  }
936 
937  if((*TpSelectorForEfficiencyVsVTXZ)(tp)){
938  for (unsigned int f=0; f<dzintervals[count].size()-1; f++){
939  if (dzSim>dzintervals[count][f]&&
940  dzSim<dzintervals[count][f+1]) {
941  totSIM_dz[count][f]++;
942  if (isMatched) {
943  totASS_dz[count][f]++;
944  }
945  }
946  } // END for (unsigned int f=0; f<dzintervals[count].size()-1; f++){
947 
948 
949  for (unsigned int f=0; f<zposintervals[count].size()-1; f++){
950  if (vertexTP.z()>zposintervals[count][f]&&vertexTP.z()<zposintervals[count][f+1]) {
951  totSIM_zpos[count][f]++;
952  if (isMatched) totASS_zpos[count][f]++;
953  }
954  if (vertz>zposintervals[count][f]&&vertz<zposintervals[count][f+1]) {
955  totSIM_vertz_entire[count][f]++;
956  if (isMatched) totASS_vertz_entire[count][f]++;
957  }
958  if (vertz>zposintervals[count][f]&&vertz<zposintervals[count][f+1] && fabs(momentumTP.eta())<0.9) {
959  totSIM_vertz_barrel[count][f]++;
960  if (isMatched) totASS_vertz_barrel[count][f]++;
961  }
962  if (vertz>zposintervals[count][f]&&vertz<zposintervals[count][f+1] && momentumTP.eta()>0.9) {
963  totSIM_vertz_fwdpos[count][f]++;
964  if (isMatched) totASS_vertz_fwdpos[count][f]++;
965  }
966  if (vertz>zposintervals[count][f]&&vertz<zposintervals[count][f+1] && momentumTP.eta()<-0.9) {
967  totSIM_vertz_fwdneg[count][f]++;
968  if (isMatched) totASS_vertz_fwdneg[count][f]++;
969  }
970  } // END for (unsigned int f=0; f<zposintervals[count].size()-1; f++){
971  }
972 
973  //Special investigations for PU
974  if(((*TpSelectorForEfficiencyVsCon)(tp)) && (!((*TpSelectorForEfficiencyVsEta)(tp)))){
975 
976  //efficPU vs eta
977  for (unsigned int f=0; f<etaintervals[count].size()-1; f++){
978  if (getEta(momentumTP.eta())>etaintervals[count][f]&&
979  getEta(momentumTP.eta())<etaintervals[count][f+1]) {
980  totCONeta[count][f]++;
981  }
982  } // END for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
983 
984  //efficPU vs num pileup vertices
985  for (unsigned int f=0; f<vertcountintervals[count].size()-1; f++){
986  if (numVertices == vertcountintervals[count][f]) {
987  totCONvertcount[count][f]++;
988  }
989  } // END for (unsigned int f=0; f<vertcountintervals[count].size()-1; f++){
990 
991  for (unsigned int f=0; f<zposintervals[count].size()-1; f++){
992  if (vertexTP.z()>zposintervals[count][f]&&vertexTP.z()<zposintervals[count][f+1]) {
993  totCONzpos[count][f]++;
994  }
995  } // END for (unsigned int f=0; f<zposintervals[count].size()-1; f++){
996 
997  }
998 
999 }
1000 
1001 // dE/dx
1002 void MTVHistoProducerAlgoForTracker::fill_dedx_recoTrack_histos(int count, edm::RefToBase<reco::Track>& trackref, const std::vector< edm::ValueMap<reco::DeDxData> >& v_dEdx) {
1003 //void MTVHistoProducerAlgoForTracker::fill_dedx_recoTrack_histos(reco::TrackRef trackref, std::vector< edm::ValueMap<reco::DeDxData> > v_dEdx) {
1004  double dedx;
1005  int nom;
1006  int sat;
1007  edm::ValueMap<reco::DeDxData> dEdxTrack;
1008  for (unsigned int i=0; i<v_dEdx.size(); i++) {
1009  dEdxTrack = v_dEdx.at(i);
1010  dedx = dEdxTrack[trackref].dEdx();
1011  nom = dEdxTrack[trackref].numberOfMeasurements();
1012  sat = dEdxTrack[trackref].numberOfSaturatedMeasurements();
1013  if (i==0) {
1014  h_dedx_estim1[count]->Fill(dedx);
1015  h_dedx_nom1[count]->Fill(nom);
1016  h_dedx_sat1[count]->Fill(sat);
1017  } else if (i==1) {
1018  h_dedx_estim2[count]->Fill(dedx);
1019  h_dedx_nom2[count]->Fill(nom);
1020  h_dedx_sat2[count]->Fill(sat);
1021  }
1022  }
1023 }
1024 
1025 
1026 // TO BE FIXED USING PLAIN HISTOGRAMS INSTEAD OF RE-IMPLEMENTATION OF HISTOGRAMS (i.d. vectors<int/double>)
1027 void MTVHistoProducerAlgoForTracker::fill_generic_recoTrack_histos(int count,
1028  const reco::Track& track,
1029  const math::XYZPoint& bsPosition,
1030  bool isMatched,
1031  bool isSigMatched,
1032  bool isChargeMatched,
1033  int numAssocRecoTracks,
1034  int numVertices,
1035  int tpbunchcrossing, int nSimHits,
1036  double sharedFraction){
1037 
1038  //Fill track algo histogram
1039 
1040  if (track.algo()>=4 && track.algo()<=14) totREC_algo[count][track.algo()-4]++;
1041  int sharedHits = sharedFraction * track.numberOfValidHits();
1042 
1043  //Compute fake rate vs eta
1044  for (unsigned int f=0; f<etaintervals[count].size()-1; f++){
1045  if (getEta(track.momentum().eta())>etaintervals[count][f]&&
1046  getEta(track.momentum().eta())<etaintervals[count][f+1]) {
1047  totRECeta[count][f]++;
1048  if (isMatched) {
1049  totASS2eta[count][f]++;
1050  if (!isChargeMatched) totmisideta[count][f]++;
1051  if (numAssocRecoTracks>1) totloopeta[count][f]++;
1052  if (tpbunchcrossing==0) totASS2_itpu_eta_entire[count][f]++;
1053  if (tpbunchcrossing!=0) totASS2_ootpu_eta_entire[count][f]++;
1054  nrecHit_vs_nsimHit_rec2sim[count]->Fill( track.numberOfValidHits(),nSimHits);
1055  h_assocFraction[count]->Fill( sharedFraction);
1056  h_assocSharedHit[count]->Fill( sharedHits);
1057  }
1058  if (isSigMatched) {
1059  totASS2etaSig[count][f]++;
1060  if (tpbunchcrossing==0) totASS2_itpu_eta_entire_signal[count][f]++;
1061  }
1062  }
1063  } // End for (unsigned int f=0; f<etaintervals[count].size()-1; f++){
1064 
1065  for (unsigned int f=0; f<phiintervals[count].size()-1; f++){
1066  if (track.momentum().phi()>phiintervals[count][f]&&
1067  track.momentum().phi()<phiintervals[count][f+1]) {
1068  totREC_phi[count][f]++;
1069  if (isMatched) {
1070  totASS2_phi[count][f]++;
1071  if (!isChargeMatched) totmisid_phi[count][f]++;
1072  if (numAssocRecoTracks>1) totloop_phi[count][f]++;
1073  }
1074  }
1075  } // End for (unsigned int f=0; f<phiintervals[count].size()-1; f++){
1076 
1077 
1078  for (unsigned int f=0; f<pTintervals[count].size()-1; f++){
1079  if (getPt(sqrt(track.momentum().perp2()))>pTintervals[count][f]&&
1080  getPt(sqrt(track.momentum().perp2()))<pTintervals[count][f+1]) {
1081  totRECpT[count][f]++;
1082  if (isMatched) {
1083  totASS2pT[count][f]++;
1084  if (!isChargeMatched) totmisidpT[count][f]++;
1085  if (numAssocRecoTracks>1) totlooppT[count][f]++;
1086  }
1087  }
1088  } // End for (unsigned int f=0; f<pTintervals[count].size()-1; f++){
1089 
1090  for (unsigned int f=0; f<dxyintervals[count].size()-1; f++){
1091  if (track.dxy(bsPosition)>dxyintervals[count][f]&&
1092  track.dxy(bsPosition)<dxyintervals[count][f+1]) {
1093  totREC_dxy[count][f]++;
1094  if (isMatched) {
1095  totASS2_dxy[count][f]++;
1096  if (!isChargeMatched) totmisid_dxy[count][f]++;
1097  if (numAssocRecoTracks>1) totloop_dxy[count][f]++;
1098  }
1099  }
1100  } // End for (unsigned int f=0; f<dxyintervals[count].size()-1; f++){
1101 
1102  for (unsigned int f=0; f<dzintervals[count].size()-1; f++){
1103  if (track.dz(bsPosition)>dzintervals[count][f]&&
1104  track.dz(bsPosition)<dzintervals[count][f+1]) {
1105  totREC_dz[count][f]++;
1106  if (isMatched) {
1107  totASS2_dz[count][f]++;
1108  if (!isChargeMatched) totmisid_dz[count][f]++;
1109  if (numAssocRecoTracks>1) totloop_dz[count][f]++;
1110  }
1111  }
1112  } // End for (unsigned int f=0; f<dzintervals[count].size()-1; f++){
1113 
1114  int tmp = std::min((int)track.found(),int(maxHit-1));
1115  totREC_hit[count][tmp]++;
1116  if (isMatched) {
1117  totASS2_hit[count][tmp]++;
1118  if (!isChargeMatched) totmisid_hit[count][tmp]++;
1119  if (numAssocRecoTracks>1) totloop_hit[count][tmp]++;
1120  }
1121 
1122  for (unsigned int f=0; f<vertcountintervals[count].size()-1; f++){
1123  if (numVertices == vertcountintervals[count][f]) {
1124  totREC_vertcount_entire[count][f]++;
1125  totREC_ootpu_entire[count][f]++;
1126  if (isMatched) {
1127  totASS2_vertcount_entire[count][f]++;
1128  if (tpbunchcrossing==0) totASS2_itpu_vertcount_entire[count][f]++;
1129  if (tpbunchcrossing!=0) totASS2_ootpu_entire[count][f]++;
1130  }
1131  if (isSigMatched) {
1132  totASS2_vertcount_entire_signal[count][f]++;
1133  if (tpbunchcrossing==0) totASS2_itpu_vertcount_entire_signal[count][f]++;
1134  }
1135  }
1136  if (numVertices == vertcountintervals[count][f] && track.eta() <= 0.9 && track.eta() >= -0.9) {
1137  totREC_vertcount_barrel[count][f]++;
1138  totREC_ootpu_barrel[count][f]++;
1139  if (isMatched) {
1140  totASS2_vertcount_barrel[count][f]++;
1141  if (isMatched && tpbunchcrossing!=0) totASS2_ootpu_barrel[count][f]++;
1142  }
1143  }
1144  if (numVertices == vertcountintervals[count][f] && track.eta() > 0.9) {
1145  totREC_vertcount_fwdpos[count][f]++;
1146  totREC_ootpu_fwdpos[count][f]++;
1147  if (isMatched) {
1148  totASS2_vertcount_fwdpos[count][f]++;
1149  if (isMatched && tpbunchcrossing!=0) totASS2_ootpu_fwdpos[count][f]++;
1150  }
1151  }
1152  if (numVertices == vertcountintervals[count][f] && track.eta() < -0.9) {
1153  totREC_vertcount_fwdneg[count][f]++;
1154  totREC_ootpu_fwdneg[count][f]++;
1155  if (isMatched) {
1156  totASS2_vertcount_fwdneg[count][f]++;
1157  if (isMatched && tpbunchcrossing!=0) totASS2_ootpu_fwdneg[count][f]++;
1158  }
1159  }
1160 
1161  }
1162 
1163 }
1164 
1165 
1166 void MTVHistoProducerAlgoForTracker::fill_simAssociated_recoTrack_histos(int count,
1167  const reco::Track& track){
1168  //nchi2 and hits global distributions
1169  h_nchi2[count]->Fill(track.normalizedChi2());
1170  h_nchi2_prob[count]->Fill(TMath::Prob(track.chi2(),(int)track.ndof()));
1171  h_hits[count]->Fill(track.numberOfValidHits());
1172  h_losthits[count]->Fill(track.numberOfLostHits());
1173  chi2_vs_nhits[count]->Fill(track.numberOfValidHits(),track.normalizedChi2());
1174  h_charge[count]->Fill( track.charge() );
1175 
1176  //chi2 and #hit vs eta: fill 2D histos
1177  chi2_vs_eta[count]->Fill(getEta(track.eta()),track.normalizedChi2());
1178  nhits_vs_eta[count]->Fill(getEta(track.eta()),track.numberOfValidHits());
1179  nPXBhits_vs_eta[count]->Fill(getEta(track.eta()),track.hitPattern().numberOfValidPixelBarrelHits());
1180  nPXFhits_vs_eta[count]->Fill(getEta(track.eta()),track.hitPattern().numberOfValidPixelEndcapHits());
1181  nTIBhits_vs_eta[count]->Fill(getEta(track.eta()),track.hitPattern().numberOfValidStripTIBHits());
1182  nTIDhits_vs_eta[count]->Fill(getEta(track.eta()),track.hitPattern().numberOfValidStripTIDHits());
1183  nTOBhits_vs_eta[count]->Fill(getEta(track.eta()),track.hitPattern().numberOfValidStripTOBHits());
1184  nTEChits_vs_eta[count]->Fill(getEta(track.eta()),track.hitPattern().numberOfValidStripTECHits());
1185  nLayersWithMeas_vs_eta[count]->Fill(getEta(track.eta()),track.hitPattern().trackerLayersWithMeasurement());
1186  nPXLlayersWithMeas_vs_eta[count]->Fill(getEta(track.eta()),track.hitPattern().pixelLayersWithMeasurement());
1187  int LayersAll = track.hitPattern().stripLayersWithMeasurement();
1188  int Layers2D = track.hitPattern().numberOfValidStripLayersWithMonoAndStereo();
1189  int Layers1D = LayersAll - Layers2D;
1190  nSTRIPlayersWithMeas_vs_eta[count]->Fill(getEta(track.eta()),LayersAll);
1191  nSTRIPlayersWith1dMeas_vs_eta[count]->Fill(getEta(track.eta()),Layers1D);
1192  nSTRIPlayersWith2dMeas_vs_eta[count]->Fill(getEta(track.eta()),Layers2D);
1193 
1194  nlosthits_vs_eta[count]->Fill(getEta(track.eta()),track.numberOfLostHits());
1195 }
1196 
1197 
1198 void MTVHistoProducerAlgoForTracker::fill_trackBased_histos(int count, int assTracks, int numRecoTracks, int numSimTracks){
1199 
1200  h_tracks[count]->Fill(assTracks);
1201  h_fakes[count]->Fill(numRecoTracks-assTracks);
1202  nrec_vs_nsim[count]->Fill(numRecoTracks,numSimTracks);
1203 
1204 }
1205 
1206 
1207 
1208 void MTVHistoProducerAlgoForTracker::fill_ResoAndPull_recoTrack_histos(int count,
1209  const TrackingParticle::Vector& momentumTP,
1210  const TrackingParticle::Point& vertexTP,
1211  int chargeTP,
1212  const reco::Track& track,
1213  const math::XYZPoint& bsPosition){
1214 
1215  // evaluation of TP parameters
1216  double qoverpSim = chargeTP/sqrt(momentumTP.x()*momentumTP.x()+momentumTP.y()*momentumTP.y()+momentumTP.z()*momentumTP.z());
1217  double lambdaSim = M_PI/2-momentumTP.theta();
1218  double phiSim = momentumTP.phi();
1219  double dxySim = (-vertexTP.x()*sin(momentumTP.phi())+vertexTP.y()*cos(momentumTP.phi()));
1220  double dzSim = vertexTP.z() - (vertexTP.x()*momentumTP.x()+vertexTP.y()*momentumTP.y())/sqrt(momentumTP.perp2())
1221  * momentumTP.z()/sqrt(momentumTP.perp2());
1222 
1223 
1224  // reco::Track::ParameterVector rParameters = track.parameters(); // UNUSED
1225 
1226  double qoverpRec(0);
1227  double qoverpErrorRec(0);
1228  double ptRec(0);
1229  double ptErrorRec(0);
1230  double lambdaRec(0);
1231  double lambdaErrorRec(0);
1232  double phiRec(0);
1233  double phiErrorRec(0);
1234 
1235  /* TO BE FIXED LATER -----------
1236  //loop to decide whether to take gsfTrack (utilisation of mode-function) or common track
1237  const GsfTrack* gsfTrack(0);
1238  if(useGsf){
1239  gsfTrack = dynamic_cast<const GsfTrack*>(&(*track));
1240  if (gsfTrack==0) edm::LogInfo("TrackValidator") << "Trying to access mode for a non-GsfTrack";
1241  }
1242 
1243  if (gsfTrack) {
1244  // get values from mode
1245  getRecoMomentum(*gsfTrack, ptRec, ptErrorRec, qoverpRec, qoverpErrorRec,
1246  lambdaRec,lambdaErrorRec, phiRec, phiErrorRec);
1247  }
1248 
1249  else {
1250  // get values from track (without mode)
1251  getRecoMomentum(*track, ptRec, ptErrorRec, qoverpRec, qoverpErrorRec,
1252  lambdaRec,lambdaErrorRec, phiRec, phiErrorRec);
1253  }
1254  */
1255  getRecoMomentum(track, ptRec, ptErrorRec, qoverpRec, qoverpErrorRec,
1256  lambdaRec,lambdaErrorRec, phiRec, phiErrorRec);
1257  // -------------
1258 
1259  double ptError = ptErrorRec;
1260  double ptres=ptRec-sqrt(momentumTP.perp2());
1261  double etares=track.eta()-momentumTP.Eta();
1262 
1263 
1264  double dxyRec = track.dxy(bsPosition);
1265  double dzRec = track.dz(bsPosition);
1266 
1267  // eta residue; pt, k, theta, phi, dxy, dz pulls
1268  double qoverpPull=(qoverpRec-qoverpSim)/qoverpErrorRec;
1269  double thetaPull=(lambdaRec-lambdaSim)/lambdaErrorRec;
1270  double phiPull=(phiRec-phiSim)/phiErrorRec;
1271  double dxyPull=(dxyRec-dxySim)/track.dxyError();
1272  double dzPull=(dzRec-dzSim)/track.dzError();
1273 
1274  double contrib_Qoverp = ((qoverpRec-qoverpSim)/qoverpErrorRec)*
1275  ((qoverpRec-qoverpSim)/qoverpErrorRec)/5;
1276  double contrib_dxy = ((dxyRec-dxySim)/track.dxyError())*((dxyRec-dxySim)/track.dxyError())/5;
1277  double contrib_dz = ((dzRec-dzSim)/track.dzError())*((dzRec-dzSim)/track.dzError())/5;
1278  double contrib_theta = ((lambdaRec-lambdaSim)/lambdaErrorRec)*
1279  ((lambdaRec-lambdaSim)/lambdaErrorRec)/5;
1280  double contrib_phi = ((phiRec-phiSim)/phiErrorRec)*
1281  ((phiRec-phiSim)/phiErrorRec)/5;
1282 
1283  LogTrace("TrackValidatorTEST")
1284  //<< "assocChi2=" << tp.begin()->second << "\n"
1285  << "" << "\n"
1286  << "ptREC=" << ptRec << "\n" << "etaREC=" << track.eta() << "\n" << "qoverpREC=" << qoverpRec << "\n"
1287  << "dxyREC=" << dxyRec << "\n" << "dzREC=" << dzRec << "\n"
1288  << "thetaREC=" << track.theta() << "\n" << "phiREC=" << phiRec << "\n"
1289  << "" << "\n"
1290  << "qoverpError()=" << qoverpErrorRec << "\n" << "dxyError()=" << track.dxyError() << "\n"<< "dzError()="
1291  << track.dzError() << "\n"
1292  << "thetaError()=" << lambdaErrorRec << "\n" << "phiError()=" << phiErrorRec << "\n"
1293  << "" << "\n"
1294  << "ptSIM=" << sqrt(momentumTP.perp2()) << "\n"<< "etaSIM=" << momentumTP.Eta() << "\n"<< "qoverpSIM=" << qoverpSim << "\n"
1295  << "dxySIM=" << dxySim << "\n"<< "dzSIM=" << dzSim << "\n" << "thetaSIM=" << M_PI/2-lambdaSim << "\n"
1296  << "phiSIM=" << phiSim << "\n"
1297  << "" << "\n"
1298  << "contrib_Qoverp=" << contrib_Qoverp << "\n"<< "contrib_dxy=" << contrib_dxy << "\n"<< "contrib_dz=" << contrib_dz << "\n"
1299  << "contrib_theta=" << contrib_theta << "\n"<< "contrib_phi=" << contrib_phi << "\n"
1300  << "" << "\n"
1301  <<"chi2PULL="<<contrib_Qoverp+contrib_dxy+contrib_dz+contrib_theta+contrib_phi<<"\n";
1302 
1303  h_pullQoverp[count]->Fill(qoverpPull);
1304  h_pullTheta[count]->Fill(thetaPull);
1305  h_pullPhi[count]->Fill(phiPull);
1306  h_pullDxy[count]->Fill(dxyPull);
1307  h_pullDz[count]->Fill(dzPull);
1308 
1309 
1310  h_pt[count]->Fill(ptres/ptError);
1311  h_eta[count]->Fill(etares);
1312  etares_vs_eta[count]->Fill(getEta(track.eta()),etares);
1313 
1314 
1315 
1316  //resolution of track params: fill 2D histos
1317  dxyres_vs_eta[count]->Fill(getEta(track.eta()),dxyRec-dxySim);
1318  ptres_vs_eta[count]->Fill(getEta(track.eta()),(ptRec-sqrt(momentumTP.perp2()))/ptRec);
1319  dzres_vs_eta[count]->Fill(getEta(track.eta()),dzRec-dzSim);
1320  phires_vs_eta[count]->Fill(getEta(track.eta()),phiRec-phiSim);
1321  cotThetares_vs_eta[count]->Fill(getEta(track.eta()),1/tan(M_PI*0.5-lambdaRec)-1/tan(M_PI*0.5-lambdaSim));
1322 
1323  //same as before but vs pT
1324  dxyres_vs_pt[count]->Fill(getPt(ptRec),dxyRec-dxySim);
1325  ptres_vs_pt[count]->Fill(getPt(ptRec),(ptRec-sqrt(momentumTP.perp2()))/ptRec);
1326  dzres_vs_pt[count]->Fill(getPt(ptRec),dzRec-dzSim);
1327  phires_vs_pt[count]->Fill(getPt(ptRec),phiRec-phiSim);
1328  cotThetares_vs_pt[count]->Fill(getPt(ptRec),1/tan(M_PI*0.5-lambdaRec)-1/tan(M_PI*0.5-lambdaSim));
1329 
1330  //pulls of track params vs eta: fill 2D histos
1331  dxypull_vs_eta[count]->Fill(getEta(track.eta()),dxyPull);
1332  ptpull_vs_eta[count]->Fill(getEta(track.eta()),ptres/ptError);
1333  dzpull_vs_eta[count]->Fill(getEta(track.eta()),dzPull);
1334  phipull_vs_eta[count]->Fill(getEta(track.eta()),phiPull);
1335  thetapull_vs_eta[count]->Fill(getEta(track.eta()),thetaPull);
1336 
1337  //plots vs phi
1338  nhits_vs_phi[count]->Fill(phiRec,track.numberOfValidHits());
1339  chi2_vs_phi[count]->Fill(phiRec,track.normalizedChi2());
1340  ptmean_vs_eta_phi[count]->Fill(phiRec,getEta(track.eta()),ptRec);
1341  phimean_vs_eta_phi[count]->Fill(phiRec,getEta(track.eta()),phiRec);
1342  ptres_vs_phi[count]->Fill(phiRec,(ptRec-sqrt(momentumTP.perp2()))/ptRec);
1343  phires_vs_phi[count]->Fill(phiRec,phiRec-phiSim);
1344  ptpull_vs_phi[count]->Fill(phiRec,ptres/ptError);
1345  phipull_vs_phi[count]->Fill(phiRec,phiPull);
1346  thetapull_vs_phi[count]->Fill(phiRec,thetaPull);
1347 
1348 
1349 }
1350 
1351 
1352 
1353 void
1354 MTVHistoProducerAlgoForTracker::getRecoMomentum (const reco::Track& track, double& pt, double& ptError,
1355  double& qoverp, double& qoverpError, double& lambda,double& lambdaError,
1356  double& phi, double& phiError ) const {
1357  pt = track.pt();
1358  ptError = track.ptError();
1359  qoverp = track.qoverp();
1360  qoverpError = track.qoverpError();
1361  lambda = track.lambda();
1362  lambdaError = track.lambdaError();
1363  phi = track.phi();
1364  phiError = track.phiError();
1365  // cout <<"test1" << endl;
1366 
1367 
1368 
1369 }
1370 
1371 void
1372 MTVHistoProducerAlgoForTracker::getRecoMomentum (const reco::GsfTrack& gsfTrack, double& pt, double& ptError,
1373  double& qoverp, double& qoverpError, double& lambda,double& lambdaError,
1374  double& phi, double& phiError ) const {
1375 
1376  pt = gsfTrack.ptMode();
1377  ptError = gsfTrack.ptModeError();
1378  qoverp = gsfTrack.qoverpMode();
1379  qoverpError = gsfTrack.qoverpModeError();
1380  lambda = gsfTrack.lambdaMode();
1381  lambdaError = gsfTrack.lambdaModeError();
1382  phi = gsfTrack.phiMode();
1383  phiError = gsfTrack.phiModeError();
1384  // cout <<"test2" << endl;
1385 
1386 }
1387 
1388 double
1389 MTVHistoProducerAlgoForTracker::getEta(double eta) {
1390  if (useFabsEta) return fabs(eta);
1391  else return eta;
1392 }
1393 
1394 double
1395 MTVHistoProducerAlgoForTracker::getPt(double pt) {
1396  if (useInvPt && pt!=0) return 1/pt;
1397  else return pt;
1398 }
1399 
1400 
1401 void MTVHistoProducerAlgoForTracker::finalHistoFits(int counter){
1402  //resolution of track params: get sigma from 2D histos
1403  FitSlicesYTool fsyt_dxy(dxyres_vs_eta[counter]);
1404  fsyt_dxy.getFittedSigmaWithError(h_dxyrmsh[counter]);
1405  fsyt_dxy.getFittedMeanWithError(h_dxymeanh[counter]);
1406  FitSlicesYTool fsyt_dxyPt(dxyres_vs_pt[counter]);
1407  fsyt_dxyPt.getFittedSigmaWithError(h_dxyrmshPt[counter]);
1408  fsyt_dxyPt.getFittedMeanWithError(h_dxymeanhPt[counter]);
1409  FitSlicesYTool fsyt_pt(ptres_vs_eta[counter]);
1410  fsyt_pt.getFittedSigmaWithError(h_ptrmsh[counter]);
1411  fsyt_pt.getFittedMeanWithError(h_ptshifteta[counter]);
1412  FitSlicesYTool fsyt_ptPt(ptres_vs_pt[counter]);
1413  fsyt_ptPt.getFittedSigmaWithError(h_ptrmshPt[counter]);
1414  fsyt_ptPt.getFittedMeanWithError(h_ptmeanhPt[counter]);
1415  FitSlicesYTool fsyt_ptPhi(ptres_vs_phi[counter]);
1416  fsyt_ptPhi.getFittedSigmaWithError(h_ptrmshPhi[counter]);
1417  fsyt_ptPhi.getFittedMeanWithError(h_ptmeanhPhi[counter]);
1418  FitSlicesYTool fsyt_dz(dzres_vs_eta[counter]);
1419  fsyt_dz.getFittedSigmaWithError(h_dzrmsh[counter]);
1420  fsyt_dz.getFittedMeanWithError(h_dzmeanh[counter]);
1421  FitSlicesYTool fsyt_dzPt(dzres_vs_pt[counter]);
1422  fsyt_dzPt.getFittedSigmaWithError(h_dzrmshPt[counter]);
1423  fsyt_dzPt.getFittedMeanWithError(h_dzmeanhPt[counter]);
1424  FitSlicesYTool fsyt_phi(phires_vs_eta[counter]);
1425  fsyt_phi.getFittedSigmaWithError(h_phirmsh[counter]);
1426  fsyt_phi.getFittedMeanWithError(h_phimeanh[counter]);
1427  FitSlicesYTool fsyt_phiPt(phires_vs_pt[counter]);
1428  fsyt_phiPt.getFittedSigmaWithError(h_phirmshPt[counter]);
1429  fsyt_phiPt.getFittedMeanWithError(h_phimeanhPt[counter]);
1430  FitSlicesYTool fsyt_phiPhi(phires_vs_phi[counter]);
1431  fsyt_phiPhi.getFittedSigmaWithError(h_phirmshPhi[counter]);
1432  fsyt_phiPhi.getFittedMeanWithError(h_phimeanhPhi[counter]);
1433  FitSlicesYTool fsyt_cotTheta(cotThetares_vs_eta[counter]);
1434  fsyt_cotTheta.getFittedSigmaWithError(h_cotThetarmsh[counter]);
1435  fsyt_cotTheta.getFittedMeanWithError(h_cotThetameanh[counter]);
1436  FitSlicesYTool fsyt_cotThetaPt(cotThetares_vs_pt[counter]);
1437  fsyt_cotThetaPt.getFittedSigmaWithError(h_cotThetarmshPt[counter]);
1438  fsyt_cotThetaPt.getFittedMeanWithError(h_cotThetameanhPt[counter]);
1439 
1440  //pulls of track params vs eta: get sigma from 2D histos
1441  FitSlicesYTool fsyt_dxyp(dxypull_vs_eta[counter]);
1442  fsyt_dxyp.getFittedSigmaWithError(h_dxypulleta[counter]);
1443  fsyt_dxyp.getFittedMeanWithError(h_dxypulletamean[counter]);
1444  FitSlicesYTool fsyt_ptp(ptpull_vs_eta[counter]);
1445  fsyt_ptp.getFittedSigmaWithError(h_ptpulleta[counter]);
1446  fsyt_ptp.getFittedMeanWithError(h_ptpulletamean[counter]);
1447  FitSlicesYTool fsyt_dzp(dzpull_vs_eta[counter]);
1448  fsyt_dzp.getFittedSigmaWithError(h_dzpulleta[counter]);
1449  fsyt_dzp.getFittedMeanWithError(h_dzpulletamean[counter]);
1450  FitSlicesYTool fsyt_phip(phipull_vs_eta[counter]);
1451  fsyt_phip.getFittedSigmaWithError(h_phipulleta[counter]);
1452  fsyt_phip.getFittedMeanWithError(h_phipulletamean[counter]);
1453  FitSlicesYTool fsyt_thetap(thetapull_vs_eta[counter]);
1454  fsyt_thetap.getFittedSigmaWithError(h_thetapulleta[counter]);
1455  fsyt_thetap.getFittedMeanWithError(h_thetapulletamean[counter]);
1456  //vs phi
1457  FitSlicesYTool fsyt_ptpPhi(ptpull_vs_phi[counter]);
1458  fsyt_ptpPhi.getFittedSigmaWithError(h_ptpullphi[counter]);
1459  fsyt_ptpPhi.getFittedMeanWithError(h_ptpullphimean[counter]);
1460  FitSlicesYTool fsyt_phipPhi(phipull_vs_phi[counter]);
1461  fsyt_phipPhi.getFittedSigmaWithError(h_phipullphi[counter]);
1462  fsyt_phipPhi.getFittedMeanWithError(h_phipullphimean[counter]);
1463  FitSlicesYTool fsyt_thetapPhi(thetapull_vs_phi[counter]);
1464  fsyt_thetapPhi.getFittedSigmaWithError(h_thetapullphi[counter]);
1465  fsyt_thetapPhi.getFittedMeanWithError(h_thetapullphimean[counter]);
1466 
1467  //effic&fake;
1468  for (unsigned int ite = 0;ite<totASSeta[counter].size();ite++) totFOMT_eta[counter][ite]=totASS2eta[counter][ite]-totASS2etaSig[counter][ite];
1469  for (unsigned int ite = 0;ite<totASS2_vertcount_entire[counter].size();ite++)
1470  totFOMT_vertcount[counter][ite]=totASS2_vertcount_entire[counter][ite]-totASS2_vertcount_entire_signal[counter][ite];
1471  for (unsigned int ite = 0;ite<totASS2_itpu_eta_entire[counter].size();ite++) totASS2_itpu_eta_entire[counter][ite]-=totASS2_itpu_eta_entire_signal[counter][ite];
1472  for (unsigned int ite = 0;ite<totASS2_itpu_vertcount_entire[counter].size();ite++) totASS2_itpu_vertcount_entire[counter][ite]-=totASS2_itpu_vertcount_entire_signal[counter][ite];
1473 
1474  fillPlotFromVectors(h_effic[counter],totASSeta[counter],totSIMeta[counter],"effic");
1475  fillPlotFromVectors(h_fakerate[counter],totASS2eta[counter],totRECeta[counter],"fakerate");
1476  fillPlotFromVectors(h_looprate[counter],totloopeta[counter],totRECeta[counter],"effic");
1477  fillPlotFromVectors(h_misidrate[counter],totmisideta[counter],totRECeta[counter],"effic");
1478  fillPlotFromVectors(h_efficPt[counter],totASSpT[counter],totSIMpT[counter],"effic");
1479  fillPlotFromVectors(h_fakeratePt[counter],totASS2pT[counter],totRECpT[counter],"fakerate");
1480  fillPlotFromVectors(h_loopratepT[counter],totlooppT[counter],totRECpT[counter],"effic");
1481  fillPlotFromVectors(h_misidratepT[counter],totmisidpT[counter],totRECpT[counter],"effic");
1482  fillPlotFromVectors(h_effic_vs_hit[counter],totASS_hit[counter],totSIM_hit[counter],"effic");
1483  fillPlotFromVectors(h_fake_vs_hit[counter],totASS2_hit[counter],totREC_hit[counter],"fakerate");
1484  fillPlotFromVectors(h_loopratehit[counter],totloop_hit[counter],totREC_hit[counter],"effic");
1485  fillPlotFromVectors(h_misidratehit[counter],totmisid_hit[counter],totREC_hit[counter],"effic");
1486  fillPlotFromVectors(h_effic_vs_phi[counter],totASS_phi[counter],totSIM_phi[counter],"effic");
1487  fillPlotFromVectors(h_fake_vs_phi[counter],totASS2_phi[counter],totREC_phi[counter],"fakerate");
1488  fillPlotFromVectors(h_loopratephi[counter],totloop_phi[counter],totREC_phi[counter],"effic");
1489  fillPlotFromVectors(h_misidratephi[counter],totmisid_phi[counter],totREC_phi[counter],"effic");
1490  fillPlotFromVectors(h_effic_vs_dxy[counter],totASS_dxy[counter],totSIM_dxy[counter],"effic");
1491  fillPlotFromVectors(h_fake_vs_dxy[counter],totASS2_dxy[counter],totREC_dxy[counter],"fakerate");
1492  fillPlotFromVectors(h_loopratedxy[counter],totloop_dxy[counter],totREC_dxy[counter],"effic");
1493  fillPlotFromVectors(h_misidratedxy[counter],totmisid_dxy[counter],totREC_dxy[counter],"effic");
1494  fillPlotFromVectors(h_effic_vs_dz[counter],totASS_dz[counter],totSIM_dz[counter],"effic");
1495  fillPlotFromVectors(h_fake_vs_dz[counter],totASS2_dz[counter],totREC_dz[counter],"fakerate");
1496  fillPlotFromVectors(h_loopratedz[counter],totloop_dz[counter],totREC_dz[counter],"effic");
1497  fillPlotFromVectors(h_misidratedz[counter],totmisid_dz[counter],totREC_dz[counter],"effic");
1498  fillPlotFromVectors(h_effic_vs_vertpos[counter],totASS_vertpos[counter],totSIM_vertpos[counter],"effic");
1499  fillPlotFromVectors(h_effic_vs_zpos[counter],totASS_zpos[counter],totSIM_zpos[counter],"effic");
1500  fillPlotFromVectors(h_effic_vertcount_entire[counter],totASS_vertcount_entire[counter],totSIM_vertcount_entire[counter],"effic");
1501  fillPlotFromVectors(h_effic_vertcount_barrel[counter],totASS_vertcount_barrel[counter],totSIM_vertcount_barrel[counter],"effic");
1502  fillPlotFromVectors(h_effic_vertcount_fwdpos[counter],totASS_vertcount_fwdpos[counter],totSIM_vertcount_fwdpos[counter],"effic");
1503  fillPlotFromVectors(h_effic_vertcount_fwdneg[counter],totASS_vertcount_fwdneg[counter],totSIM_vertcount_fwdneg[counter],"effic");
1504  fillPlotFromVectors(h_fakerate_vertcount_entire[counter],totASS2_vertcount_entire[counter],totREC_vertcount_entire[counter],"fakerate");
1505  fillPlotFromVectors(h_fakerate_vertcount_barrel[counter],totASS2_vertcount_barrel[counter],totREC_vertcount_barrel[counter],"fakerate");
1506  fillPlotFromVectors(h_fakerate_vertcount_fwdpos[counter],totASS2_vertcount_fwdpos[counter],totREC_vertcount_fwdpos[counter],"fakerate");
1507  fillPlotFromVectors(h_fakerate_vertcount_fwdneg[counter],totASS2_vertcount_fwdneg[counter],totREC_vertcount_fwdneg[counter],"fakerate");
1508  fillPlotFromVectors(h_effic_vertz_entire[counter],totASS_vertz_entire[counter],totSIM_vertz_entire[counter],"effic");
1509  fillPlotFromVectors(h_effic_vertz_barrel[counter],totASS_vertz_barrel[counter],totSIM_vertz_barrel[counter],"effic");
1510  fillPlotFromVectors(h_effic_vertz_fwdpos[counter],totASS_vertz_fwdpos[counter],totSIM_vertz_fwdpos[counter],"effic");
1511  fillPlotFromVectors(h_effic_vertz_fwdneg[counter],totASS_vertz_fwdneg[counter],totSIM_vertz_fwdneg[counter],"effic");
1512  fillPlotFromVectors(h_fakerate_ootpu_entire[counter],totASS2_ootpu_entire[counter],totREC_ootpu_entire[counter],"effic");
1513  fillPlotFromVectors(h_fakerate_ootpu_barrel[counter],totASS2_ootpu_barrel[counter],totREC_ootpu_barrel[counter],"effic");
1514  fillPlotFromVectors(h_fakerate_ootpu_fwdpos[counter],totASS2_ootpu_fwdpos[counter],totREC_ootpu_fwdpos[counter],"effic");
1515  fillPlotFromVectors(h_fakerate_ootpu_fwdneg[counter],totASS2_ootpu_fwdneg[counter],totREC_ootpu_fwdneg[counter],"effic");
1516 
1517  fillPlotFromVectors(h_fomt_eta[counter],totFOMT_eta[counter],totASS2etaSig[counter],"pileup");
1518  fillPlotFromVectors(h_fomt_vertcount[counter],totFOMT_vertcount[counter],totASS2_vertcount_entire_signal[counter],"pileup");
1519  fillPlotFromVectors(h_fomt_sig_eta[counter],totASS2etaSig[counter],totRECeta[counter],"fakerate");
1520  fillPlotFromVectors(h_fomt_sig_vertcount[counter],totASS2_vertcount_entire_signal[counter],totREC_vertcount_entire[counter],"fakerate");
1521  fillPlotFromVectors(h_fomt_itpu_eta[counter],totASS2_itpu_eta_entire[counter],totRECeta[counter],"effic");
1522  fillPlotFromVectors(h_fomt_itpu_vertcount[counter],totASS2_itpu_vertcount_entire[counter],totREC_vertcount_entire[counter],"effic");
1523  fillPlotFromVectors(h_fomt_ootpu_eta[counter],totASS2_ootpu_eta_entire[counter],totRECeta[counter],"effic");
1524  fillPlotFromVectors(h_fomt_ootpu_vertcount[counter],totASS2_ootpu_entire[counter],totREC_ootpu_entire[counter],"effic");
1525 
1526  fillPlotFromVectors(h_effic_PU_eta[counter],totASSeta[counter],totCONeta[counter],"pileup");
1527  fillPlotFromVectors(h_effic_PU_vertcount[counter],totASS_vertcount_entire[counter],totCONvertcount[counter],"pileup");
1528  fillPlotFromVectors(h_effic_PU_zpos[counter],totASS_zpos[counter],totCONzpos[counter],"pileup");
1529 
1530 }
1531 
1532 void MTVHistoProducerAlgoForTracker::fillProfileHistosFromVectors(int counter){
1533  //chi2 and #hit vs eta: get mean from 2D histos
1534  doProfileX(chi2_vs_eta[counter],h_chi2meanh[counter]);
1535  doProfileX(nhits_vs_eta[counter],h_hits_eta[counter]);
1536  doProfileX(nPXBhits_vs_eta[counter],h_PXBhits_eta[counter]);
1537  doProfileX(nPXFhits_vs_eta[counter],h_PXFhits_eta[counter]);
1538  doProfileX(nTIBhits_vs_eta[counter],h_TIBhits_eta[counter]);
1539  doProfileX(nTIDhits_vs_eta[counter],h_TIDhits_eta[counter]);
1540  doProfileX(nTOBhits_vs_eta[counter],h_TOBhits_eta[counter]);
1541  doProfileX(nTEChits_vs_eta[counter],h_TEChits_eta[counter]);
1542 
1543  doProfileX(nLayersWithMeas_vs_eta[counter],h_LayersWithMeas_eta[counter]);
1544  doProfileX(nPXLlayersWithMeas_vs_eta[counter],h_PXLlayersWithMeas_eta[counter]);
1545  doProfileX(nSTRIPlayersWithMeas_vs_eta[counter],h_STRIPlayersWithMeas_eta[counter]);
1546  doProfileX(nSTRIPlayersWith1dMeas_vs_eta[counter],h_STRIPlayersWith1dMeas_eta[counter]);
1547  doProfileX(nSTRIPlayersWith2dMeas_vs_eta[counter],h_STRIPlayersWith2dMeas_eta[counter]);
1548 
1549 
1550 
1551  doProfileX(nlosthits_vs_eta[counter],h_losthits_eta[counter]);
1552  //vs phi
1553  doProfileX(chi2_vs_nhits[counter],h_chi2meanhitsh[counter]);
1554  // doProfileX(ptres_vs_eta[counter],h_ptresmean_vs_eta[counter]);
1555  // doProfileX(phires_vs_eta[counter],h_phiresmean_vs_eta[counter]);
1556  doProfileX(chi2_vs_phi[counter],h_chi2mean_vs_phi[counter]);
1557  doProfileX(nhits_vs_phi[counter],h_hits_phi[counter]);
1558  // doProfileX(ptres_vs_phi[counter],h_ptresmean_vs_phi[counter]);
1559  // doProfileX(phires_vs_phi[counter],h_phiresmean_vs_phi[counter]);
1560 }
1561 
1562 void MTVHistoProducerAlgoForTracker::fillHistosFromVectors(int counter){
1563  fillPlotFromVector(h_algo[counter],totREC_algo[counter]);
1564 
1565  fillPlotFromVector(h_recoeta[counter],totRECeta[counter]);
1566  fillPlotFromVector(h_simuleta[counter],totSIMeta[counter]);
1567  fillPlotFromVector(h_assoceta[counter],totASSeta[counter]);
1568  fillPlotFromVector(h_assoc2eta[counter],totASS2eta[counter]);
1569  fillPlotFromVector(h_loopereta[counter],totloopeta[counter]);
1570  fillPlotFromVector(h_misideta[counter],totmisideta[counter]);
1571 
1572  fillPlotFromVector(h_recopT[counter],totRECpT[counter]);
1573  fillPlotFromVector(h_simulpT[counter],totSIMpT[counter]);
1574  fillPlotFromVector(h_assocpT[counter],totASSpT[counter]);
1575  fillPlotFromVector(h_assoc2pT[counter],totASS2pT[counter]);
1576  fillPlotFromVector(h_looperpT[counter],totlooppT[counter]);
1577  fillPlotFromVector(h_misidpT[counter],totmisidpT[counter]);
1578 
1579  fillPlotFromVector(h_recohit[counter],totREC_hit[counter]);
1580  fillPlotFromVector(h_simulhit[counter],totSIM_hit[counter]);
1581  fillPlotFromVector(h_assochit[counter],totASS_hit[counter]);
1582  fillPlotFromVector(h_assoc2hit[counter],totASS2_hit[counter]);
1583  fillPlotFromVector(h_looperhit[counter],totloop_hit[counter]);
1584  fillPlotFromVector(h_misidhit[counter],totmisid_hit[counter]);
1585 
1586  fillPlotFromVector(h_recophi[counter],totREC_phi[counter]);
1587  fillPlotFromVector(h_simulphi[counter],totSIM_phi[counter]);
1588  fillPlotFromVector(h_assocphi[counter],totASS_phi[counter]);
1589  fillPlotFromVector(h_assoc2phi[counter],totASS2_phi[counter]);
1590  fillPlotFromVector(h_looperphi[counter],totloop_phi[counter]);
1591  fillPlotFromVector(h_misidphi[counter],totmisid_phi[counter]);
1592 
1593  fillPlotFromVector(h_recodxy[counter],totREC_dxy[counter]);
1594  fillPlotFromVector(h_simuldxy[counter],totSIM_dxy[counter]);
1595  fillPlotFromVector(h_assocdxy[counter],totASS_dxy[counter]);
1596  fillPlotFromVector(h_assoc2dxy[counter],totASS2_dxy[counter]);
1597  fillPlotFromVector(h_looperdxy[counter],totloop_dxy[counter]);
1598  fillPlotFromVector(h_misiddxy[counter],totmisid_dxy[counter]);
1599 
1600  fillPlotFromVector(h_recodz[counter],totREC_dz[counter]);
1601  fillPlotFromVector(h_simuldz[counter],totSIM_dz[counter]);
1602  fillPlotFromVector(h_assocdz[counter],totASS_dz[counter]);
1603  fillPlotFromVector(h_assoc2dz[counter],totASS2_dz[counter]);
1604  fillPlotFromVector(h_looperdz[counter],totloop_dz[counter]);
1605  fillPlotFromVector(h_misiddz[counter],totmisid_dz[counter]);
1606 
1607  fillPlotFromVector(h_simulvertpos[counter],totSIM_vertpos[counter]);
1608  fillPlotFromVector(h_assocvertpos[counter],totASS_vertpos[counter]);
1609 
1610  fillPlotFromVector(h_simulzpos[counter],totSIM_zpos[counter]);
1611  fillPlotFromVector(h_assoczpos[counter],totASS_zpos[counter]);
1612 
1613  fillPlotFromVector(h_reco_vertcount_entire[counter],totREC_vertcount_entire[counter]);
1614  fillPlotFromVector(h_simul_vertcount_entire[counter],totSIM_vertcount_entire[counter]);
1615  fillPlotFromVector(h_assoc_vertcount_entire[counter],totASS_vertcount_entire[counter]);
1616  fillPlotFromVector(h_assoc2_vertcount_entire[counter],totASS2_vertcount_entire[counter]);
1617 
1618  fillPlotFromVector(h_reco_vertcount_barrel[counter],totREC_vertcount_barrel[counter]);
1619  fillPlotFromVector(h_simul_vertcount_barrel[counter],totSIM_vertcount_barrel[counter]);
1620  fillPlotFromVector(h_assoc_vertcount_barrel[counter],totASS_vertcount_barrel[counter]);
1621  fillPlotFromVector(h_assoc2_vertcount_barrel[counter],totASS2_vertcount_barrel[counter]);
1622 
1623  fillPlotFromVector(h_reco_vertcount_fwdpos[counter],totREC_vertcount_fwdpos[counter]);
1624  fillPlotFromVector(h_simul_vertcount_fwdpos[counter],totSIM_vertcount_fwdpos[counter]);
1625  fillPlotFromVector(h_assoc_vertcount_fwdpos[counter],totASS_vertcount_fwdpos[counter]);
1626  fillPlotFromVector(h_assoc2_vertcount_fwdpos[counter],totASS2_vertcount_fwdpos[counter]);
1627 
1628  fillPlotFromVector(h_reco_vertcount_fwdneg[counter],totREC_vertcount_fwdneg[counter]);
1629  fillPlotFromVector(h_simul_vertcount_fwdneg[counter],totSIM_vertcount_fwdneg[counter]);
1630  fillPlotFromVector(h_assoc_vertcount_fwdneg[counter],totASS_vertcount_fwdneg[counter]);
1631  fillPlotFromVector(h_assoc2_vertcount_fwdneg[counter],totASS2_vertcount_fwdneg[counter]);
1632 
1633  fillPlotFromVector(h_simul_vertz_entire[counter],totSIM_vertz_entire[counter]);
1634  fillPlotFromVector(h_assoc_vertz_entire[counter],totASS_vertz_entire[counter]);
1635 
1636  fillPlotFromVector(h_simul_vertz_barrel[counter],totSIM_vertz_barrel[counter]);
1637  fillPlotFromVector(h_assoc_vertz_barrel[counter],totASS_vertz_barrel[counter]);
1638 
1639  fillPlotFromVector(h_simul_vertz_fwdpos[counter],totSIM_vertz_fwdpos[counter]);
1640  fillPlotFromVector(h_assoc_vertz_fwdpos[counter],totASS_vertz_fwdpos[counter]);
1641 
1642  fillPlotFromVector(h_simul_vertz_fwdneg[counter],totSIM_vertz_fwdneg[counter]);
1643  fillPlotFromVector(h_assoc_vertz_fwdneg[counter],totASS_vertz_fwdneg[counter]);
1644 
1645  fillPlotFromVector(h_reco_ootpu_entire[counter],totREC_ootpu_entire[counter]);
1646  fillPlotFromVector(h_assoc2_ootpu_entire[counter],totASS2_ootpu_entire[counter]);
1647 
1648  fillPlotFromVector(h_reco_ootpu_barrel[counter],totREC_ootpu_barrel[counter]);
1649  fillPlotFromVector(h_assoc2_ootpu_barrel[counter],totASS2_ootpu_barrel[counter]);
1650 
1651  fillPlotFromVector(h_reco_ootpu_fwdpos[counter],totREC_ootpu_fwdpos[counter]);
1652  fillPlotFromVector(h_assoc2_ootpu_fwdpos[counter],totASS2_ootpu_fwdpos[counter]);
1653 
1654  fillPlotFromVector(h_reco_ootpu_fwdneg[counter],totREC_ootpu_fwdneg[counter]);
1655  fillPlotFromVector(h_assoc2_ootpu_fwdneg[counter],totASS2_ootpu_fwdneg[counter]);
1656 
1657  fillPlotFromVector(h_con_eta[counter],totCONeta[counter]);
1658  fillPlotFromVector(h_con_vertcount[counter],totCONvertcount[counter]);
1659  fillPlotFromVector(h_con_zpos[counter],totCONzpos[counter]);
1660 
1661 }
1662 
1663 
1664 
1665 
1666 
1667 
1668 void MTVHistoProducerAlgoForTracker::fill_recoAssociated_simTrack_histos(int count,
1669  const reco::GenParticle& tp,
1670  const TrackingParticle::Vector& momentumTP,
1671  const TrackingParticle::Point& vertexTP,
1672  double dxySim, double dzSim, int nSimHits,
1673  const reco::Track* track,
1674  int numVertices, double vertz){
1675  bool isMatched = track;
1676 
1677  if((*GpSelectorForEfficiencyVsEta)(tp)){
1678  //effic vs hits
1679  int nSimHitsInBounds = std::min((int)nSimHits,int(maxHit-1));
1680  totSIM_hit[count][nSimHitsInBounds]++;
1681  if(isMatched) {
1682  totASS_hit[count][nSimHitsInBounds]++;
1683  nrecHit_vs_nsimHit_sim2rec[count]->Fill( track->numberOfValidHits(),nSimHits);
1684  }
1685 
1686  //effic vs eta
1687  for (unsigned int f=0; f<etaintervals[count].size()-1; f++){
1688  if (getEta(momentumTP.eta())>etaintervals[count][f]&&
1689  getEta(momentumTP.eta())<etaintervals[count][f+1]) {
1690  totSIMeta[count][f]++;
1691  if (isMatched) {
1692  totASSeta[count][f]++;
1693  }
1694  }
1695 
1696  } // END for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
1697 
1698  //effic vs num pileup vertices
1699  for (unsigned int f=0; f<vertcountintervals[count].size()-1; f++){
1700  if (numVertices == vertcountintervals[count][f]) {
1701  totSIM_vertcount_entire[count][f]++;
1702  if (isMatched) {
1703  totASS_vertcount_entire[count][f]++;
1704  }
1705  }
1706  if (numVertices == vertcountintervals[count][f] && momentumTP.eta() <= 0.9 && momentumTP.eta() >= -0.9) {
1707  totSIM_vertcount_barrel[count][f]++;
1708  if (isMatched) {
1709  totASS_vertcount_barrel[count][f]++;
1710  }
1711  }
1712  if (numVertices == vertcountintervals[count][f] && momentumTP.eta() > 0.9) {
1713  totSIM_vertcount_fwdpos[count][f]++;
1714  if (isMatched) {
1715  totASS_vertcount_fwdpos[count][f]++;
1716  }
1717  }
1718  if (numVertices == vertcountintervals[count][f] && momentumTP.eta() < -0.9) {
1719  totSIM_vertcount_fwdneg[count][f]++;
1720  if (isMatched) {
1721  totASS_vertcount_fwdneg[count][f]++;
1722  }
1723  }
1724  }
1725 
1726  }
1727 
1728  if((*GpSelectorForEfficiencyVsPhi)(tp)){
1729  for (unsigned int f=0; f<phiintervals[count].size()-1; f++){
1730  if (momentumTP.phi() > phiintervals[count][f]&&
1731  momentumTP.phi() <phiintervals[count][f+1]) {
1732  totSIM_phi[count][f]++;
1733  if (isMatched) {
1734  totASS_phi[count][f]++;
1735  }
1736  }
1737  } // END for (unsigned int f=0; f<phiintervals[count].size()-1; f++){
1738  }
1739 
1740  if((*GpSelectorForEfficiencyVsPt)(tp)){
1741  for (unsigned int f=0; f<pTintervals[count].size()-1; f++){
1742  if (getPt(sqrt(momentumTP.perp2()))>pTintervals[count][f]&&
1743  getPt(sqrt(momentumTP.perp2()))<pTintervals[count][f+1]) {
1744  totSIMpT[count][f]++;
1745  if (isMatched) {
1746  totASSpT[count][f]++;
1747  }
1748  }
1749  } // END for (unsigned int f=0; f<pTintervals[count].size()-1; f++){
1750  }
1751 
1752  if((*GpSelectorForEfficiencyVsVTXR)(tp)){
1753  for (unsigned int f=0; f<dxyintervals[count].size()-1; f++){
1754  if (dxySim>dxyintervals[count][f]&&
1755  dxySim<dxyintervals[count][f+1]) {
1756  totSIM_dxy[count][f]++;
1757  if (isMatched) {
1758  totASS_dxy[count][f]++;
1759  }
1760  }
1761  } // END for (unsigned int f=0; f<dxyintervals[count].size()-1; f++){
1762 
1763  for (unsigned int f=0; f<vertposintervals[count].size()-1; f++){
1764  if (sqrt(vertexTP.perp2())>vertposintervals[count][f]&&
1765  sqrt(vertexTP.perp2())<vertposintervals[count][f+1]) {
1766  totSIM_vertpos[count][f]++;
1767  if (isMatched) {
1768  totASS_vertpos[count][f]++;
1769  }
1770  }
1771  } // END for (unsigned int f=0; f<vertposintervals[count].size()-1; f++){
1772  }
1773 
1774  if((*GpSelectorForEfficiencyVsVTXZ)(tp)){
1775  for (unsigned int f=0; f<dzintervals[count].size()-1; f++){
1776  if (dzSim>dzintervals[count][f]&&
1777  dzSim<dzintervals[count][f+1]) {
1778  totSIM_dz[count][f]++;
1779  if (isMatched) {
1780  totASS_dz[count][f]++;
1781  }
1782  }
1783  } // END for (unsigned int f=0; f<dzintervals[count].size()-1; f++){
1784 
1785 
1786  for (unsigned int f=0; f<zposintervals[count].size()-1; f++){
1787  if (vertexTP.z()>zposintervals[count][f]&&vertexTP.z()<zposintervals[count][f+1]) {
1788  totSIM_zpos[count][f]++;
1789  if (isMatched) totASS_zpos[count][f]++;
1790  }
1791  if (vertz>zposintervals[count][f]&&vertz<zposintervals[count][f+1]) {
1792  totSIM_vertz_entire[count][f]++;
1793  if (isMatched) totASS_vertz_entire[count][f]++;
1794  }
1795  if (vertz>zposintervals[count][f]&&vertz<zposintervals[count][f+1] && fabs(momentumTP.eta())<0.9) {
1796  totSIM_vertz_barrel[count][f]++;
1797  if (isMatched) totASS_vertz_barrel[count][f]++;
1798  }
1799  if (vertz>zposintervals[count][f]&&vertz<zposintervals[count][f+1] && momentumTP.eta()>0.9) {
1800  totSIM_vertz_fwdpos[count][f]++;
1801  if (isMatched) totASS_vertz_fwdpos[count][f]++;
1802  }
1803  if (vertz>zposintervals[count][f]&&vertz<zposintervals[count][f+1] && momentumTP.eta()<-0.9) {
1804  totSIM_vertz_fwdneg[count][f]++;
1805  if (isMatched) totASS_vertz_fwdneg[count][f]++;
1806  }
1807  } // END for (unsigned int f=0; f<zposintervals[count].size()-1; f++){
1808  }
1809 
1810  //Special investigations for PU
1811  if(((*GpSelectorForEfficiencyVsCon)(tp)) && (!((*GpSelectorForEfficiencyVsEta)(tp)))){
1812 
1813  //efficPU vs eta
1814  for (unsigned int f=0; f<etaintervals[count].size()-1; f++){
1815  if (getEta(momentumTP.eta())>etaintervals[count][f]&&
1816  getEta(momentumTP.eta())<etaintervals[count][f+1]) {
1817  totCONeta[count][f]++;
1818  }
1819  } // END for (unsigned int f=0; f<etaintervals[w].size()-1; f++){
1820 
1821  //efficPU vs num pileup vertices
1822  for (unsigned int f=0; f<vertcountintervals[count].size()-1; f++){
1823  if (numVertices == vertcountintervals[count][f]) {
1824  totCONvertcount[count][f]++;
1825  }
1826  } // END for (unsigned int f=0; f<vertcountintervals[count].size()-1; f++){
1827 
1828  for (unsigned int f=0; f<zposintervals[count].size()-1; f++){
1829  if (vertexTP.z()>zposintervals[count][f]&&vertexTP.z()<zposintervals[count][f+1]) {
1830  totCONzpos[count][f]++;
1831  }
1832  } // END for (unsigned int f=0; f<zposintervals[count].size()-1; f++){
1833 
1834  }
1835 
1836 }
MTVHistoProducerAlgoForTracker::h_assoc_vertz_fwdpos
std::vector< MonitorElement * > h_assoc_vertz_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:169
MTVHistoProducerAlgoForTracker::totFOMT_vertcount
std::vector< std::vector< int > > totFOMT_vertcount
Definition: MTVHistoProducerAlgoForTracker.h:251
reco::TrackBase::qoverp
double qoverp() const
q/p
Definition: TrackBase.h:115
MTVHistoProducerAlgoForTracker::h_fakerate_ootpu_entire
std::vector< MonitorElement * > h_fakerate_ootpu_entire
Definition: MTVHistoProducerAlgoForTracker.h:172
MTVHistoProducerAlgoForTracker::vertcountintervals
std::vector< std::vector< double > > vertcountintervals
Definition: MTVHistoProducerAlgoForTracker.h:227
reco::GsfTrack::phiModeError
double phiModeError() const
error on phi from mode
Definition: GsfTrack.h:93
MTVHistoProducerAlgoForTracker::h_simul_vertz_entire
std::vector< MonitorElement * > h_simul_vertz_entire
Definition: MTVHistoProducerAlgoForTracker.h:167
MTVHistoProducerAlgoForTracker::h_looperhit
std::vector< MonitorElement * > h_looperhit
Definition: MTVHistoProducerAlgoForTracker.h:156
MTVHistoProducerAlgoForTracker::h_loopratephi
std::vector< MonitorElement * > h_loopratephi
Definition: MTVHistoProducerAlgoForTracker.h:157
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
MTVHistoProducerAlgoForTracker::h_fakerate_vertcount_fwdpos
std::vector< MonitorElement * > h_fakerate_vertcount_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:165
MTVHistoProducerAlgoForTracker::h_assoc2hit
std::vector< MonitorElement * > h_assoc2hit
Definition: MTVHistoProducerAlgoForTracker.h:156
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
MTVHistoProducerAlgoForTracker::h_dedx_estim2
std::vector< MonitorElement * > h_dedx_estim2
Definition: MTVHistoProducerAlgoForTracker.h:190
i
int i
Definition: DBlmapReader.cc:9
MTVHistoProducerAlgoForTracker::h_dzpulletamean
std::vector< MonitorElement * > h_dzpulletamean
Definition: MTVHistoProducerAlgoForTracker.h:295
MTVHistoProducerAlgoForTracker::h_recodxy
std::vector< MonitorElement * > h_recodxy
Definition: MTVHistoProducerAlgoForTracker.h:158
MTVHistoProducerAlgoForTracker::totASS_vertcount_fwdpos
std::vector< std::vector< int > > totASS_vertcount_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:238
MTVHistoProducerAlgoForTracker::h_fomt_itpu_eta
std::vector< MonitorElement * > h_fomt_itpu_eta
Definition: MTVHistoProducerAlgoForTracker.h:177
MTVHistoProducerAlgoForTracker::h_assochit
std::vector< MonitorElement * > h_assochit
Definition: MTVHistoProducerAlgoForTracker.h:156
MTVHistoProducerAlgoForTracker::h_efficPt
std::vector< MonitorElement * > h_efficPt
Definition: MTVHistoProducerAlgoForTracker.h:155
MTVHistoProducerAlgoForTracker::h_fakes
std::vector< MonitorElement * > h_fakes
Definition: MTVHistoProducerAlgoForTracker.h:153
MTVHistoProducerAlgoForTracker::zposintervals
std::vector< std::vector< double > > zposintervals
Definition: MTVHistoProducerAlgoForTracker.h:226
MTVHistoProducerAlgoForTracker::h_loopratedxy
std::vector< MonitorElement * > h_loopratedxy
Definition: MTVHistoProducerAlgoForTracker.h:158
MTVHistoProducerAlgoForTracker::h_PXFhits_eta
std::vector< MonitorElement * > h_PXFhits_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::totASS2_ootpu_fwdpos
std::vector< std::vector< int > > totASS2_ootpu_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:247
MTVHistoProducerAlgoForTracker::h_nchi2_prob
std::vector< MonitorElement * > h_nchi2_prob
Definition: MTVHistoProducerAlgoForTracker.h:258
MTVHistoProducerAlgoForTracker::fill_dedx_recoTrack_histos
void fill_dedx_recoTrack_histos(int count, edm::RefToBase< reco::Track > &trackref, const std::vector< edm::ValueMap< reco::DeDxData > > &v_dEdx)
Definition: MTVHistoProducerAlgoForTracker.cc:1002
MTVHistoProducerAlgoForTracker::h_looperpT
std::vector< MonitorElement * > h_looperpT
Definition: MTVHistoProducerAlgoForTracker.h:155
MTVHistoProducerAlgoForTracker::h_pullDz
std::vector< MonitorElement * > h_pullDz
Definition: MTVHistoProducerAlgoForTracker.h:162
FitSlicesYTool::getFittedSigmaWithError
void getFittedSigmaWithError(MonitorElement *)
Fill the ME with the sigma value (with error) of the gaussian fit in each slice.
Definition: FitSlicesYTool.cc:71
MTVHistoProducerAlgoForTracker::TpSelectorForEfficiencyVsPhi
TrackingParticleSelector * TpSelectorForEfficiencyVsPhi
Definition: MTVHistoProducerAlgoForTracker.h:117
MTVHistoProducerAlgo::fillPlotFromVectors
void fillPlotFromVectors(MonitorElement *h, std::vector< int > &numerator, std::vector< int > &denominator, std::string type)
Definition: MTVHistoProducerAlgo.cc:20
MTVHistoProducerAlgoForTracker::ptpull_vs_phi
std::vector< MonitorElement * > ptpull_vs_phi
Definition: MTVHistoProducerAlgoForTracker.h:292
MTVHistoProducerAlgoForTracker::h_assoczpos
std::vector< MonitorElement * > h_assoczpos
Definition: MTVHistoProducerAlgoForTracker.h:161
MTVHistoProducerAlgoForTracker::totASS_phi
std::vector< std::vector< int > > totASS_phi
Definition: MTVHistoProducerAlgoForTracker.h:231
MTVHistoProducerAlgoForTracker::nPXFhits_vs_eta
std::vector< MonitorElement * > nPXFhits_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::totSIM_vertz_fwdpos
std::vector< std::vector< int > > totSIM_vertz_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:242
reco::TrackBase::momentum
const Vector & momentum() const
track momentum vector
Definition: TrackBase.h:150
MTVHistoProducerAlgoForTracker::h_fake_vs_hit
std::vector< MonitorElement * > h_fake_vs_hit
Definition: MTVHistoProducerAlgoForTracker.h:156
MTVHistoProducerAlgoForTracker::h_dzmeanh
std::vector< MonitorElement * > h_dzmeanh
Definition: MTVHistoProducerAlgoForTracker.h:286
MTVHistoProducerAlgoForTracker::h_assoc2_ootpu_vertcount
std::vector< MonitorElement * > h_assoc2_ootpu_vertcount
Definition: MTVHistoProducerAlgoForTracker.h:181
MTVHistoProducerAlgoForTracker::totRECeta
std::vector< std::vector< int > > totRECeta
Definition: MTVHistoProducerAlgoForTracker.h:228
MTVHistoProducerAlgoForTracker::h_misidratephi
std::vector< MonitorElement * > h_misidratephi
Definition: MTVHistoProducerAlgoForTracker.h:157
MTVHistoProducerAlgo::doProfileX
void doProfileX(TH2 *th2, MonitorElement *me)
Definition: MTVHistoProducerAlgo.cc:3
MTVHistoProducerAlgoForTracker::h_simulzpos
std::vector< MonitorElement * > h_simulzpos
Definition: MTVHistoProducerAlgoForTracker.h:161
MTVHistoProducerAlgoForTracker::totASSeta
std::vector< std::vector< int > > totASSeta
Definition: MTVHistoProducerAlgoForTracker.h:228
MTVHistoProducerAlgoForTracker::nrecHit_vs_nsimHit_rec2sim
std::vector< MonitorElement * > nrecHit_vs_nsimHit_rec2sim
Definition: MTVHistoProducerAlgoForTracker.h:199
MTVHistoProducerAlgoForTracker::vertposintervals
std::vector< std::vector< double > > vertposintervals
Definition: MTVHistoProducerAlgoForTracker.h:225
MTVHistoProducerAlgoForTracker::totASS2_itpu_eta_entire
std::vector< std::vector< int > > totASS2_itpu_eta_entire
Definition: MTVHistoProducerAlgoForTracker.h:250
MTVHistoProducerAlgoForTracker::totCONeta
std::vector< std::vector< int > > totCONeta
Definition: MTVHistoProducerAlgoForTracker.h:252
MTVHistoProducerAlgoForTracker::h_effic_vertz_fwdpos
std::vector< MonitorElement * > h_effic_vertz_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:169
MTVHistoProducerAlgoForTracker::h_assoc2pT
std::vector< MonitorElement * > h_assoc2pT
Definition: MTVHistoProducerAlgoForTracker.h:155
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:722
reco::GsfTrack::lambdaMode
double lambdaMode() const
Lambda angle from mode.
Definition: GsfTrack.h:44
MTVHistoProducerAlgoForTracker::h_simulpT
std::vector< MonitorElement * > h_simulpT
Definition: MTVHistoProducerAlgoForTracker.h:155
MTVHistoProducerAlgoForTracker::h_effic_vertz_barrel
std::vector< MonitorElement * > h_effic_vertz_barrel
Definition: MTVHistoProducerAlgoForTracker.h:168
MTVHistoProducerAlgo::fillPlotFromVector
void fillPlotFromVector(MonitorElement *h, std::vector< int > &vec)
Definition: MTVHistoProducerAlgo.cc:14
MTVHistoProducerAlgoForTracker::h_ptshifteta
std::vector< MonitorElement * > h_ptshifteta
Definition: MTVHistoProducerAlgoForTracker.h:262
reco::TrackBase::normalizedChi2
double normalizedChi2() const
chi-squared divided by n.d.o.f. (or chi-squared * 1e6 if n.d.o.f. is zero)
Definition: TrackBase.h:111
MTVHistoProducerAlgoForTracker::getRecoMomentum
void getRecoMomentum(const reco::Track &track, double &pt, double &ptError, double &qoverp, double &qoverpError, double &lambda, double &lambdaError, double &phi, double &phiError) const
retrieval of reconstructed momentum components from reco::Track (== mean values for GSF) ...
Definition: MTVHistoProducerAlgoForTracker.cc:1354
MTVHistoProducerAlgoForTracker::GpSelectorForEfficiencyVsPhi
GenParticleCustomSelector * GpSelectorForEfficiencyVsPhi
Definition: MTVHistoProducerAlgoForTracker.h:125
MTVHistoProducerAlgoForTracker::nTEChits_vs_eta
std::vector< MonitorElement * > nTEChits_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::h_reco_ootpu_entire
std::vector< MonitorElement * > h_reco_ootpu_entire
Definition: MTVHistoProducerAlgoForTracker.h:172
MTVHistoProducerAlgoForTracker::h_fakerate_vertcount_barrel
std::vector< MonitorElement * > h_fakerate_vertcount_barrel
Definition: MTVHistoProducerAlgoForTracker.h:164
MTVHistoProducerAlgoForTracker::totSIM_vertz_entire
std::vector< std::vector< int > > totSIM_vertz_entire
Definition: MTVHistoProducerAlgoForTracker.h:240
MTVHistoProducerAlgoForTracker::h_fake_vs_dz
std::vector< MonitorElement * > h_fake_vs_dz
Definition: MTVHistoProducerAlgoForTracker.h:159
MTVHistoProducerAlgoForTracker::h_hits
std::vector< MonitorElement * > h_hits
Definition: MTVHistoProducerAlgoForTracker.h:153
MTVHistoProducerAlgoForTracker::h_assoc_vertcount_fwdneg
std::vector< MonitorElement * > h_assoc_vertcount_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:166
MTVHistoProducerAlgoForTracker::h_effic_vs_vertpos
std::vector< MonitorElement * > h_effic_vs_vertpos
Definition: MTVHistoProducerAlgoForTracker.h:161
MTVHistoProducerAlgoForTracker::dxyintervals
std::vector< std::vector< double > > dxyintervals
Definition: MTVHistoProducerAlgoForTracker.h:223
reco::TrackBase::theta
double theta() const
polar angle
Definition: TrackBase.h:117
reco::TrackBase::dxyError
double dxyError() const
error on dxy
Definition: TrackBase.h:209
MTVHistoProducerAlgoForTracker::nLayersWithMeas_vs_eta
std::vector< MonitorElement * > nLayersWithMeas_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::h_pullPhi
std::vector< MonitorElement * > h_pullPhi
Definition: MTVHistoProducerAlgoForTracker.h:162
MTVHistoProducerAlgoForTracker::h_fomt_vertcount
std::vector< MonitorElement * > h_fomt_vertcount
Definition: MTVHistoProducerAlgoForTracker.h:176
MTVHistoProducerAlgoForTracker::totASS_vertcount_barrel
std::vector< std::vector< int > > totASS_vertcount_barrel
Definition: MTVHistoProducerAlgoForTracker.h:237
MTVHistoProducerAlgoForTracker::h_fake_vs_dxy
std::vector< MonitorElement * > h_fake_vs_dxy
Definition: MTVHistoProducerAlgoForTracker.h:158
MTVHistoProducerAlgoForTracker::phipull_vs_phi
std::vector< MonitorElement * > phipull_vs_phi
Definition: MTVHistoProducerAlgoForTracker.h:292
MTVHistoProducerAlgoForTracker::h_dedx_sat1
std::vector< MonitorElement * > h_dedx_sat1
Definition: MTVHistoProducerAlgoForTracker.h:193
MTVHistoProducerAlgoForTracker::h_hits_eta
std::vector< MonitorElement * > h_hits_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::h_effic_vertcount_fwdpos
std::vector< MonitorElement * > h_effic_vertcount_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:165
MTVHistoProducerAlgoForTracker::h_thetapullphi
std::vector< MonitorElement * > h_thetapullphi
Definition: MTVHistoProducerAlgoForTracker.h:294
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
MTVHistoProducerAlgoForTracker::nrec_vs_nsim
std::vector< MonitorElement * > nrec_vs_nsim
Definition: MTVHistoProducerAlgoForTracker.h:197
MTVHistoProducerAlgoForTracker::h_ptrmshPt
std::vector< MonitorElement * > h_ptrmshPt
Definition: MTVHistoProducerAlgoForTracker.h:282
MTVHistoProducerAlgoForTracker::fill_generic_recoTrack_histos
void fill_generic_recoTrack_histos(int count, const reco::Track &track, const math::XYZPoint &bsPosition, bool isMatched, bool isSigMatched, bool isChargeMatched, int numAssocRecoTracks, int numVertices, int tpbunchcrossing, int nSimHits, double sharedFraction)
Definition: MTVHistoProducerAlgoForTracker.cc:1027
MTVHistoProducerAlgoForTracker::h_fakerate_ootpu_barrel
std::vector< MonitorElement * > h_fakerate_ootpu_barrel
Definition: MTVHistoProducerAlgoForTracker.h:173
MTVHistoProducerAlgoForTracker::h_algo
std::vector< MonitorElement * > h_algo
Definition: MTVHistoProducerAlgoForTracker.h:171
MTVHistoProducerAlgoForTracker::h_assocdz
std::vector< MonitorElement * > h_assocdz
Definition: MTVHistoProducerAlgoForTracker.h:159
MTVHistoProducerAlgoForTracker::totREC_dz
std::vector< std::vector< int > > totREC_dz
Definition: MTVHistoProducerAlgoForTracker.h:233
MTVHistoProducerAlgoForTracker::totASS2_phi
std::vector< std::vector< int > > totASS2_phi
Definition: MTVHistoProducerAlgoForTracker.h:231
reco::HitPattern::pixelLayersWithMeasurement
int pixelLayersWithMeasurement() const
Definition: HitPattern.h:726
MTVHistoProducerAlgoForTracker::h_assoc_vertcount_fwdpos
std::vector< MonitorElement * > h_assoc_vertcount_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:165
MTVHistoProducerAlgoForTracker::nTIDhits_vs_eta
std::vector< MonitorElement * > nTIDhits_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::totREC_vertcount_entire
std::vector< std::vector< int > > totREC_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:236
MTVHistoProducerAlgoForTracker::h_ptSIM
std::vector< MonitorElement * > h_ptSIM
Definition: MTVHistoProducerAlgoForTracker.h:150
MTVHistoProducerAlgoForTracker::h_dxypulletamean
std::vector< MonitorElement * > h_dxypulletamean
Definition: MTVHistoProducerAlgoForTracker.h:295
MTVHistoProducerAlgoForTracker::totASS2_vertcount_entire
std::vector< std::vector< int > > totASS2_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:236
MTVHistoProducerAlgoForTracker::totREC_algo
std::vector< std::vector< int > > totREC_algo
Definition: MTVHistoProducerAlgoForTracker.h:244
MTVHistoProducerAlgoForTracker::phires_vs_eta
std::vector< MonitorElement * > phires_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:278
MTVHistoProducerAlgoForTracker::h_pullDxy
std::vector< MonitorElement * > h_pullDxy
Definition: MTVHistoProducerAlgoForTracker.h:162
MTVHistoProducerAlgoForTracker::h_assoc2phi
std::vector< MonitorElement * > h_assoc2phi
Definition: MTVHistoProducerAlgoForTracker.h:157
MTVHistoProducerAlgoForTracker::totASS_dxy
std::vector< std::vector< int > > totASS_dxy
Definition: MTVHistoProducerAlgoForTracker.h:232
MTVHistoProducerAlgoForTracker::h_phipullphi
std::vector< MonitorElement * > h_phipullphi
Definition: MTVHistoProducerAlgoForTracker.h:294
MTVHistoProducerAlgoForTracker::h_nchi2
std::vector< MonitorElement * > h_nchi2
Definition: MTVHistoProducerAlgoForTracker.h:258
MTVHistoProducerAlgoForTracker::h_fakeratePt
std::vector< MonitorElement * > h_fakeratePt
Definition: MTVHistoProducerAlgoForTracker.h:155
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:139
reco::TrackBase::numberOfLostHits
unsigned short numberOfLostHits() const
number of cases where track crossed a layer without getting a hit.
Definition: TrackBase.h:234
MTVHistoProducerAlgoForTracker::totSIM_vertcount_entire
std::vector< std::vector< int > > totSIM_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:236
MTVHistoProducerAlgoForTracker::TpSelectorForEfficiencyVsVTXZ
TrackingParticleSelector * TpSelectorForEfficiencyVsVTXZ
Definition: MTVHistoProducerAlgoForTracker.h:120
MTVHistoProducerAlgoForTracker::h_dxyrmsh
std::vector< MonitorElement * > h_dxyrmsh
Definition: MTVHistoProducerAlgoForTracker.h:281
min
#define min(a, b)
Definition: mlp_lapack.h:161
MTVHistoProducerAlgoForTracker::h_simul_vertcount_fwdpos
std::vector< MonitorElement * > h_simul_vertcount_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:165
MTVHistoProducerAlgoForTracker::totmisid_dxy
std::vector< std::vector< int > > totmisid_dxy
Definition: MTVHistoProducerAlgoForTracker.h:232
MTVHistoProducerAlgoForTracker::h_dzrmshPt
std::vector< MonitorElement * > h_dzrmshPt
Definition: MTVHistoProducerAlgoForTracker.h:282
MTVHistoProducerAlgoForTracker::h_vertposSIM
std::vector< MonitorElement * > h_vertposSIM
Definition: MTVHistoProducerAlgoForTracker.h:150
MTVHistoProducerAlgoForTracker::h_looperdxy
std::vector< MonitorElement * > h_looperdxy
Definition: MTVHistoProducerAlgoForTracker.h:158
MTVHistoProducerAlgoForTracker::h_simuleta
std::vector< MonitorElement * > h_simuleta
Definition: MTVHistoProducerAlgoForTracker.h:154
MTVHistoProducerAlgoForTracker::h_misidratepT
std::vector< MonitorElement * > h_misidratepT
Definition: MTVHistoProducerAlgoForTracker.h:155
MTVHistoProducerAlgoForTracker::totSIM_dz
std::vector< std::vector< int > > totSIM_dz
Definition: MTVHistoProducerAlgoForTracker.h:233
MTVHistoProducerAlgoForTracker::fill_recoAssociated_simTrack_histos
void fill_recoAssociated_simTrack_histos(int count, const TrackingParticle &tp, const TrackingParticle::Vector &momentumTP, const TrackingParticle::Point &vertexTP, double dxy, double dz, int nSimHits, const reco::Track *track, int numVertices, double vertz)
Definition: MTVHistoProducerAlgoForTracker.cc:831
MTVHistoProducerAlgoForTracker::h_ptpulleta
std::vector< MonitorElement * > h_ptpulleta
Definition: MTVHistoProducerAlgoForTracker.h:293
MTVHistoProducerAlgoForTracker::h_effic_PU_eta
std::vector< MonitorElement * > h_effic_PU_eta
Definition: MTVHistoProducerAlgoForTracker.h:183
MTVHistoProducerAlgoForTracker::h_assoc2_ootpu_fwdneg
std::vector< MonitorElement * > h_assoc2_ootpu_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:175
MTVHistoProducerAlgoForTracker::h_assoc2_vertcount_fwdneg
std::vector< MonitorElement * > h_assoc2_vertcount_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:166
MTVHistoProducerAlgoForTracker::h_fomt_sig_vertcount
std::vector< MonitorElement * > h_fomt_sig_vertcount
Definition: MTVHistoProducerAlgoForTracker.h:176
eta
T eta() const
Definition: Basic3DVectorLD.h:179
MTVHistoProducerAlgoForTracker::h_assoc_vertz_fwdneg
std::vector< MonitorElement * > h_assoc_vertz_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:170
reco::HitPattern::numberOfValidStripLayersWithMonoAndStereo
int numberOfValidStripLayersWithMonoAndStereo() const
Definition: HitPattern.cc:224
MTVHistoProducerAlgoForTracker::h_TIBhits_eta
std::vector< MonitorElement * > h_TIBhits_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::h_thetapullphimean
std::vector< MonitorElement * > h_thetapullphimean
Definition: MTVHistoProducerAlgoForTracker.h:296
MTVHistoProducerAlgoForTracker::h_fomt_ootpu_eta
std::vector< MonitorElement * > h_fomt_ootpu_eta
Definition: MTVHistoProducerAlgoForTracker.h:178
MTVHistoProducerAlgoForTracker::h_cotThetameanhPt
std::vector< MonitorElement * > h_cotThetameanhPt
Definition: MTVHistoProducerAlgoForTracker.h:287
MTVHistoProducerAlgoForTracker::totFOMT_eta
std::vector< std::vector< int > > totFOMT_eta
Definition: MTVHistoProducerAlgoForTracker.h:251
MTVHistoProducerAlgoForTracker::h_phimeanhPhi
std::vector< MonitorElement * > h_phimeanhPhi
Definition: MTVHistoProducerAlgoForTracker.h:288
MTVHistoProducerAlgoForTracker::fill_trackBased_histos
void fill_trackBased_histos(int count, int assTracks, int numRecoTracks, int numSimTracks)
Definition: MTVHistoProducerAlgoForTracker.cc:1198
MTVHistoProducerAlgoForTracker::h_misidratedz
std::vector< MonitorElement * > h_misidratedz
Definition: MTVHistoProducerAlgoForTracker.h:159
MTVHistoProducerAlgoForTracker::h_assoc2dz
std::vector< MonitorElement * > h_assoc2dz
Definition: MTVHistoProducerAlgoForTracker.h:159
MTVHistoProducerAlgoForTracker::chi2_vs_phi
std::vector< MonitorElement * > chi2_vs_phi
Definition: MTVHistoProducerAlgoForTracker.h:263
MTVHistoProducerAlgoForTracker::h_fomt_eta
std::vector< MonitorElement * > h_fomt_eta
Definition: MTVHistoProducerAlgoForTracker.h:176
MTVHistoProducerAlgoForTracker::h_LayersWithMeas_eta
std::vector< MonitorElement * > h_LayersWithMeas_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
reco::HitPattern::numberOfValidStripTOBHits
int numberOfValidStripTOBHits() const
Definition: HitPattern.h:606
MTVHistoProducerAlgoForTracker::totASS_zpos
std::vector< std::vector< int > > totASS_zpos
Definition: MTVHistoProducerAlgoForTracker.h:235
MTVHistoProducerAlgoForTracker::totloop_dxy
std::vector< std::vector< int > > totloop_dxy
Definition: MTVHistoProducerAlgoForTracker.h:232
MTVHistoProducerAlgoForTracker::totASS2etaSig
std::vector< std::vector< int > > totASS2etaSig
Definition: MTVHistoProducerAlgoForTracker.h:228
MTVHistoProducerAlgoForTracker::totASS2_vertcount_barrel
std::vector< std::vector< int > > totASS2_vertcount_barrel
Definition: MTVHistoProducerAlgoForTracker.h:237
MTVHistoProducerAlgoForTracker::h_chi2mean_vs_phi
std::vector< MonitorElement * > h_chi2mean_vs_phi
Definition: MTVHistoProducerAlgoForTracker.h:275
MTVHistoProducerAlgoForTracker::h_assocphi
std::vector< MonitorElement * > h_assocphi
Definition: MTVHistoProducerAlgoForTracker.h:157
EgAmbiguityTools::sharedHits
int sharedHits(const reco::GsfTrackRef &, const reco::GsfTrackRef &)
Definition: EgAmbiguityTools.cc:66
MTVHistoProducerAlgoForTracker::h_cotThetarmsh
std::vector< MonitorElement * > h_cotThetarmsh
Definition: MTVHistoProducerAlgoForTracker.h:281
MTVHistoProducerAlgoForTracker::phires_vs_phi
std::vector< MonitorElement * > phires_vs_phi
Definition: MTVHistoProducerAlgoForTracker.h:263
MTVHistoProducerAlgoForTracker::dzres_vs_eta
std::vector< MonitorElement * > dzres_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:278
MTVHistoProducerAlgoForTracker::h_losthits_eta
std::vector< MonitorElement * > h_losthits_eta
Definition: MTVHistoProducerAlgoForTracker.h:273
reco::GsfTrack::qoverpMode
double qoverpMode() const
q/p from mode
Definition: GsfTrack.h:40
MTVHistoProducerAlgoForTracker::h_ptmeanhPhi
std::vector< MonitorElement * > h_ptmeanhPhi
Definition: MTVHistoProducerAlgoForTracker.h:288
MTVHistoProducerAlgoForTracker::totSIM_hit
std::vector< std::vector< int > > totSIM_hit
Definition: MTVHistoProducerAlgoForTracker.h:230
MTVHistoProducerAlgoForTracker::h_misideta
std::vector< MonitorElement * > h_misideta
Definition: MTVHistoProducerAlgoForTracker.h:154
MTVHistoProducerAlgoForTracker::nlosthits_vs_eta
std::vector< MonitorElement * > nlosthits_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:272
MTVHistoProducerAlgoForTracker::GpSelectorForEfficiencyVsVTXZ
GenParticleCustomSelector * GpSelectorForEfficiencyVsVTXZ
Definition: MTVHistoProducerAlgoForTracker.h:128
MTVHistoProducerAlgoForTracker::nSTRIPlayersWith1dMeas_vs_eta
std::vector< MonitorElement * > nSTRIPlayersWith1dMeas_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::h_phimeanh
std::vector< MonitorElement * > h_phimeanh
Definition: MTVHistoProducerAlgoForTracker.h:286
reco::GsfTrack::ptModeError
double ptModeError() const
error on Pt (set to 1000 TeV if charge==0 for safety) from mode
Definition: GsfTrack.h:80
MTVHistoProducerAlgoForTracker::h_recohit
std::vector< MonitorElement * > h_recohit
Definition: MTVHistoProducerAlgoForTracker.h:156
reco::TrackBase::algo
TrackAlgorithm algo() const
Definition: TrackBase.h:332
MTVHistoProducerAlgoForTracker::h_reco_vertcount_entire
std::vector< MonitorElement * > h_reco_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:163
MTVHistoProducerAlgoForTracker::h_thetapulleta
std::vector< MonitorElement * > h_thetapulleta
Definition: MTVHistoProducerAlgoForTracker.h:293
TrackingParticle::Point
math::XYZPointD Point
point in the space
Definition: TrackingParticle.h:36
MTVHistoProducerAlgoForTracker::h_fakerate_vertcount_entire
std::vector< MonitorElement * > h_fakerate_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:163
MTVHistoProducerAlgoForTracker::totSIM_zpos
std::vector< std::vector< int > > totSIM_zpos
Definition: MTVHistoProducerAlgoForTracker.h:235
MTVHistoProducerAlgoForTracker::totASSpT
std::vector< std::vector< int > > totASSpT
Definition: MTVHistoProducerAlgoForTracker.h:229
MTVHistoProducerAlgoForTracker::totSIM_vertz_fwdneg
std::vector< std::vector< int > > totSIM_vertz_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:243
MTVHistoProducerAlgoForTracker::totASS2_vertcount_fwdneg
std::vector< std::vector< int > > totASS2_vertcount_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:239
MTVHistoProducerAlgoForTracker::h_phipullphimean
std::vector< MonitorElement * > h_phipullphimean
Definition: MTVHistoProducerAlgoForTracker.h:296
MTVHistoProducerAlgoForTracker::h_dedx_estim1
std::vector< MonitorElement * > h_dedx_estim1
Definition: MTVHistoProducerAlgoForTracker.h:189
MTVHistoProducerAlgoForTracker::h_hits_phi
std::vector< MonitorElement * > h_hits_phi
Definition: MTVHistoProducerAlgoForTracker.h:274
MTVHistoProducerAlgoForTracker::totASS2_ootpu_fwdneg
std::vector< std::vector< int > > totASS2_ootpu_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:248
MTVHistoProducerAlgoForTracker::phires_vs_pt
std::vector< MonitorElement * > phires_vs_pt
Definition: MTVHistoProducerAlgoForTracker.h:279
MTVHistoProducerAlgoForTracker::h_reco_ootpu_fwdneg
std::vector< MonitorElement * > h_reco_ootpu_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:175
reco::modules::TrackingParticleSelector
SingleObjectSelector< TrackingParticleCollection,::TrackingParticleSelector > TrackingParticleSelector
Definition: TrackingParticleSelector.cc:18
MTVHistoProducerAlgoForTracker::totlooppT
std::vector< std::vector< int > > totlooppT
Definition: MTVHistoProducerAlgoForTracker.h:229
MTVHistoProducerAlgoForTracker::totloop_hit
std::vector< std::vector< int > > totloop_hit
Definition: MTVHistoProducerAlgoForTracker.h:230
MTVHistoProducerAlgoForTracker::fill_ResoAndPull_recoTrack_histos
void fill_ResoAndPull_recoTrack_histos(int count, const TrackingParticle::Vector &momentumTP, const TrackingParticle::Point &vertexTP, int chargeTP, const reco::Track &track, const math::XYZPoint &bsPosition)
Definition: MTVHistoProducerAlgoForTracker.cc:1208
MTVHistoProducerAlgoForTracker::h_STRIPlayersWithMeas_eta
std::vector< MonitorElement * > h_STRIPlayersWithMeas_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::totASS2eta
std::vector< std::vector< int > > totASS2eta
Definition: MTVHistoProducerAlgoForTracker.h:228
MTVHistoProducerAlgoForTracker::h_assoc2_ootpu_fwdpos
std::vector< MonitorElement * > h_assoc2_ootpu_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:174
MTVHistoProducerAlgoForTracker::nPXBhits_vs_eta
std::vector< MonitorElement * > nPXBhits_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:141
MTVHistoProducerAlgoForTracker::totASS_hit
std::vector< std::vector< int > > totASS_hit
Definition: MTVHistoProducerAlgoForTracker.h:230
MTVHistoProducerAlgoForTracker::totASS2_vertcount_entire_signal
std::vector< std::vector< int > > totASS2_vertcount_entire_signal
Definition: MTVHistoProducerAlgoForTracker.h:236
MTVHistoProducerAlgoForTracker::totASS2_itpu_vertcount_entire
std::vector< std::vector< int > > totASS2_itpu_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:250
MTVHistoProducerAlgoForTracker::h_assoc2dxy
std::vector< MonitorElement * > h_assoc2dxy
Definition: MTVHistoProducerAlgoForTracker.h:158
MTVHistoProducerAlgoForTracker::h_dxyrmshPt
std::vector< MonitorElement * > h_dxyrmshPt
Definition: MTVHistoProducerAlgoForTracker.h:282
MTVHistoProducerAlgoForTracker::h_STRIPlayersWith1dMeas_eta
std::vector< MonitorElement * > h_STRIPlayersWith1dMeas_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::h_pullQoverp
std::vector< MonitorElement * > h_pullQoverp
Definition: MTVHistoProducerAlgoForTracker.h:162
reco::HitPattern::numberOfValidPixelBarrelHits
int numberOfValidPixelBarrelHits() const
Definition: HitPattern.h:586
MTVHistoProducerAlgoForTracker::totSIM_phi
std::vector< std::vector< int > > totSIM_phi
Definition: MTVHistoProducerAlgoForTracker.h:231
MTVHistoProducerAlgoForTracker::totASS2pT
std::vector< std::vector< int > > totASS2pT
Definition: MTVHistoProducerAlgoForTracker.h:229
MTVHistoProducerAlgoForTracker::h_fake_vs_phi
std::vector< MonitorElement * > h_fake_vs_phi
Definition: MTVHistoProducerAlgoForTracker.h:157
MTVHistoProducerAlgoForTracker::h_ptpulletamean
std::vector< MonitorElement * > h_ptpulletamean
Definition: MTVHistoProducerAlgoForTracker.h:295
reco::TrackBase::chi2
double chi2() const
chi-squared of the fit
Definition: TrackBase.h:107
MTVHistoProducerAlgoForTracker::h_misidhit
std::vector< MonitorElement * > h_misidhit
Definition: MTVHistoProducerAlgoForTracker.h:156
MTVHistoProducerAlgoForTracker::nSTRIPlayersWithMeas_vs_eta
std::vector< MonitorElement * > nSTRIPlayersWithMeas_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::nTOBhits_vs_eta
std::vector< MonitorElement * > nTOBhits_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::totASS_dz
std::vector< std::vector< int > > totASS_dz
Definition: MTVHistoProducerAlgoForTracker.h:233
reco::TrackBase::ndof
double ndof() const
number of degrees of freedom of the fit
Definition: TrackBase.h:109
reco::HitPattern::trackerLayersWithMeasurement
int trackerLayersWithMeasurement() const
Definition: HitPattern.h:721
MTVHistoProducerAlgoForTracker::totASS_vertz_barrel
std::vector< std::vector< int > > totASS_vertz_barrel
Definition: MTVHistoProducerAlgoForTracker.h:241
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
MTVHistoProducerAlgoForTracker::ptpull_vs_eta
std::vector< MonitorElement * > ptpull_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:291
MTVHistoProducerAlgoForTracker::etares_vs_eta
std::vector< MonitorElement * > etares_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:261
MTVHistoProducerAlgoForTracker::h_fakerate_ootpu_fwdneg
std::vector< MonitorElement * > h_fakerate_ootpu_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:175
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:131
MTVHistoProducerAlgoForTracker::h_recophi
std::vector< MonitorElement * > h_recophi
Definition: MTVHistoProducerAlgoForTracker.h:157
MTVHistoProducerAlgoForTracker::h_misidratehit
std::vector< MonitorElement * > h_misidratehit
Definition: MTVHistoProducerAlgoForTracker.h:156
MTVHistoProducerAlgoForTracker::h_assocpT
std::vector< MonitorElement * > h_assocpT
Definition: MTVHistoProducerAlgoForTracker.h:155
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
MTVHistoProducerAlgoForTracker::h_ptpullphimean
std::vector< MonitorElement * > h_ptpullphimean
Definition: MTVHistoProducerAlgoForTracker.h:296
reco::TrackBase::ptError
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:194
MTVHistoProducerAlgoForTracker::h_simul_vertz_barrel
std::vector< MonitorElement * > h_simul_vertz_barrel
Definition: MTVHistoProducerAlgoForTracker.h:168
MTVHistoProducerAlgoForTracker::totASS_vertcount_entire
std::vector< std::vector< int > > totASS_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:236
MTVHistoProducerAlgoForTracker::nSTRIPlayersWith2dMeas_vs_eta
std::vector< MonitorElement * > nSTRIPlayersWith2dMeas_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::h_tracks
std::vector< MonitorElement * > h_tracks
Definition: MTVHistoProducerAlgoForTracker.h:153
MTVHistoProducerAlgoForTracker::h_chi2meanh
std::vector< MonitorElement * > h_chi2meanh
Definition: MTVHistoProducerAlgoForTracker.h:273
reco::TrackBase::phiError
double phiError() const
error on phi
Definition: TrackBase.h:207
MTVHistoProducerAlgoForTracker::h_TIDhits_eta
std::vector< MonitorElement * > h_TIDhits_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::h_tracksSIM
std::vector< MonitorElement * > h_tracksSIM
Definition: MTVHistoProducerAlgoForTracker.h:150
MTVHistoProducerAlgoForTracker::h_effic_vertcount_entire
std::vector< MonitorElement * > h_effic_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:163
MTVHistoProducerAlgoForTracker::h_effic
std::vector< MonitorElement * > h_effic
Definition: MTVHistoProducerAlgoForTracker.h:154
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
MTVHistoProducerAlgoForTracker::nhits_vs_phi
std::vector< MonitorElement * > nhits_vs_phi
Definition: MTVHistoProducerAlgoForTracker.h:263
MTVHistoProducerAlgoForTracker::totloopeta
std::vector< std::vector< int > > totloopeta
Definition: MTVHistoProducerAlgoForTracker.h:228
MTVHistoProducerAlgoForTracker::totASS2_dz
std::vector< std::vector< int > > totASS2_dz
Definition: MTVHistoProducerAlgoForTracker.h:233
reco::HitPattern::numberOfValidStripTIDHits
int numberOfValidStripTIDHits() const
Definition: HitPattern.h:602
reco::TrackBase::lambda
double lambda() const
Lambda angle.
Definition: TrackBase.h:119
MTVHistoProducerAlgoForTracker::totSIM_vertpos
std::vector< std::vector< int > > totSIM_vertpos
Definition: MTVHistoProducerAlgoForTracker.h:235
MTVHistoProducerAlgoForTracker::h_simul_vertcount_barrel
std::vector< MonitorElement * > h_simul_vertcount_barrel
Definition: MTVHistoProducerAlgoForTracker.h:164
MTVHistoProducerAlgoForTracker::totREC_ootpu_entire
std::vector< std::vector< int > > totREC_ootpu_entire
Definition: MTVHistoProducerAlgoForTracker.h:245
MTVHistoProducerAlgoForTracker::h_ptmeanhPt
std::vector< MonitorElement * > h_ptmeanhPt
Definition: MTVHistoProducerAlgoForTracker.h:287
MTVHistoProducerAlgoForTracker::h_phirmsh
std::vector< MonitorElement * > h_phirmsh
Definition: MTVHistoProducerAlgoForTracker.h:281
f
double f[11][100]
Definition: MuScleFitUtils.cc:79
MTVHistoProducerAlgoForTracker::nrecHit_vs_nsimHit_sim2rec
std::vector< MonitorElement * > nrecHit_vs_nsimHit_sim2rec
Definition: MTVHistoProducerAlgoForTracker.h:198
trackerHitRTTI::isMatched
bool isMatched(TrackingRecHit const &hit)
Definition: BaseTrackerRecHit.h:17
MTVHistoProducerAlgoForTracker::totRECpT
std::vector< std::vector< int > > totRECpT
Definition: MTVHistoProducerAlgoForTracker.h:229
MTVHistoProducerAlgoForTracker::totmisid_phi
std::vector< std::vector< int > > totmisid_phi
Definition: MTVHistoProducerAlgoForTracker.h:231
MTVHistoProducerAlgoForTracker::totREC_vertcount_fwdneg
std::vector< std::vector< int > > totREC_vertcount_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:239
reco::TrackBase::numberOfValidHits
unsigned short numberOfValidHits() const
number of valid hits found
Definition: TrackBase.h:232
MTVHistoProducerAlgoForTracker::fill_simAssociated_recoTrack_histos
void fill_simAssociated_recoTrack_histos(int count, const reco::Track &track)
Definition: MTVHistoProducerAlgoForTracker.cc:1166
reco::HitPattern::numberOfValidStripTECHits
int numberOfValidStripTECHits() const
Definition: HitPattern.h:610
MTVHistoProducerAlgoForTracker::totmisidpT
std::vector< std::vector< int > > totmisidpT
Definition: MTVHistoProducerAlgoForTracker.h:229
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:1036
MTVHistoProducerAlgoForTracker::h_reco_ootpu_barrel
std::vector< MonitorElement * > h_reco_ootpu_barrel
Definition: MTVHistoProducerAlgoForTracker.h:173
MTVHistoProducerAlgoForTracker::h_assoc_vertz_barrel
std::vector< MonitorElement * > h_assoc_vertz_barrel
Definition: MTVHistoProducerAlgoForTracker.h:168
MTVHistoProducerAlgoForTracker::h_effic_PU_zpos
std::vector< MonitorElement * > h_effic_PU_zpos
Definition: MTVHistoProducerAlgoForTracker.h:183
MTVHistoProducerAlgoForTracker::h_phirmshPt
std::vector< MonitorElement * > h_phirmshPt
Definition: MTVHistoProducerAlgoForTracker.h:282
MTVHistoProducerAlgoForTracker::totloop_dz
std::vector< std::vector< int > > totloop_dz
Definition: MTVHistoProducerAlgoForTracker.h:233
MTVHistoProducerAlgoForTracker::phimean_vs_eta_phi
std::vector< MonitorElement * > phimean_vs_eta_phi
Definition: MTVHistoProducerAlgoForTracker.h:266
MTVHistoProducerAlgoForTracker::h_charge
std::vector< MonitorElement * > h_charge
Definition: MTVHistoProducerAlgoForTracker.h:153
reco::TrackBase::hitPattern
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:223
MTVHistoProducerAlgoForTracker::totASS2_itpu_vertcount_entire_signal
std::vector< std::vector< int > > totASS2_itpu_vertcount_entire_signal
Definition: MTVHistoProducerAlgoForTracker.h:250
MTVHistoProducerAlgoForTracker::h_reco_ootpu_eta
std::vector< MonitorElement * > h_reco_ootpu_eta
Definition: MTVHistoProducerAlgoForTracker.h:182
MTVHistoProducerAlgoForTracker::totASS2_hit
std::vector< std::vector< int > > totASS2_hit
Definition: MTVHistoProducerAlgoForTracker.h:230
MTVHistoProducerAlgoForTracker::h_loopratehit
std::vector< MonitorElement * > h_loopratehit
Definition: MTVHistoProducerAlgoForTracker.h:156
MTVHistoProducerAlgoForTracker::h_loopereta
std::vector< MonitorElement * > h_loopereta
Definition: MTVHistoProducerAlgoForTracker.h:154
MTVHistoProducerAlgoForTracker::h_effic_vertcount_fwdneg
std::vector< MonitorElement * > h_effic_vertcount_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:166
MTVHistoProducerAlgoForTracker::h_misiddz
std::vector< MonitorElement * > h_misiddz
Definition: MTVHistoProducerAlgoForTracker.h:159
MTVHistoProducerAlgoForTracker::totREC_dxy
std::vector< std::vector< int > > totREC_dxy
Definition: MTVHistoProducerAlgoForTracker.h:232
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
MTVHistoProducerAlgoForTracker::phiintervals
std::vector< std::vector< double > > phiintervals
Definition: MTVHistoProducerAlgoForTracker.h:222
MTVHistoProducerAlgoForTracker::h_assoc_vertcount_entire
std::vector< MonitorElement * > h_assoc_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:163
MTVHistoProducerAlgoForTracker::chi2_vs_nhits
std::vector< MonitorElement * > chi2_vs_nhits
Definition: MTVHistoProducerAlgoForTracker.h:261
MTVHistoProducerAlgoForTracker::totSIM_vertcount_fwdpos
std::vector< std::vector< int > > totSIM_vertcount_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:238
MTVHistoProducerAlgoForTracker::totSIM_vertz_barrel
std::vector< std::vector< int > > totSIM_vertz_barrel
Definition: MTVHistoProducerAlgoForTracker.h:241
MTVHistoProducerAlgoForTracker::totASS2_itpu_eta_entire_signal
std::vector< std::vector< int > > totASS2_itpu_eta_entire_signal
Definition: MTVHistoProducerAlgoForTracker.h:250
MTVHistoProducerAlgoForTracker::dzpull_vs_eta
std::vector< MonitorElement * > dzpull_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:291
MTVHistoProducerAlgoForTracker::totREC_vertcount_barrel
std::vector< std::vector< int > > totREC_vertcount_barrel
Definition: MTVHistoProducerAlgoForTracker.h:237
MTVHistoProducerAlgoForTracker::h_simuldz
std::vector< MonitorElement * > h_simuldz
Definition: MTVHistoProducerAlgoForTracker.h:159
MTVHistoProducerAlgoForTracker::h_misiddxy
std::vector< MonitorElement * > h_misiddxy
Definition: MTVHistoProducerAlgoForTracker.h:158
MTVHistoProducerAlgoForTracker::totASS_vertz_fwdpos
std::vector< std::vector< int > > totASS_vertz_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:242
MTVHistoProducerAlgoForTracker::h_bunchxSIM
std::vector< MonitorElement * > h_bunchxSIM
Definition: MTVHistoProducerAlgoForTracker.h:150
gen::k
int k[5][pyjets_maxn]
Definition: Cascade2Hadronizer.cc:81
MTVHistoProducerAlgoForTracker::h_fakerate
std::vector< MonitorElement * > h_fakerate
Definition: MTVHistoProducerAlgoForTracker.h:154
reco::TrackBase::qoverpError
double qoverpError() const
error on signed transverse curvature
Definition: TrackBase.h:192
reco::TrackBase::dz
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:127
reco::TrackBase::dzError
double dzError() const
error on dz
Definition: TrackBase.h:215
MTVHistoProducerAlgoForTracker::totmisid_hit
std::vector< std::vector< int > > totmisid_hit
Definition: MTVHistoProducerAlgoForTracker.h:230
MTVHistoProducerAlgoForTracker::h_ptpullphi
std::vector< MonitorElement * > h_ptpullphi
Definition: MTVHistoProducerAlgoForTracker.h:294
MTVHistoProducerAlgoForTracker::totCONzpos
std::vector< std::vector< int > > totCONzpos
Definition: MTVHistoProducerAlgoForTracker.h:252
MTVHistoProducerAlgoForTracker::thetapull_vs_phi
std::vector< MonitorElement * > thetapull_vs_phi
Definition: MTVHistoProducerAlgoForTracker.h:292
MTVHistoProducerAlgoForTracker::h_assoc2_ootpu_barrel
std::vector< MonitorElement * > h_assoc2_ootpu_barrel
Definition: MTVHistoProducerAlgoForTracker.h:173
MTVHistoProducerAlgoForTracker::TpSelectorForEfficiencyVsPt
TrackingParticleSelector * TpSelectorForEfficiencyVsPt
Definition: MTVHistoProducerAlgoForTracker.h:118
MTVHistoProducerAlgoForTracker::totREC_ootpu_barrel
std::vector< std::vector< int > > totREC_ootpu_barrel
Definition: MTVHistoProducerAlgoForTracker.h:246
MTVHistoProducerAlgoForTracker::h_simulvertpos
std::vector< MonitorElement * > h_simulvertpos
Definition: MTVHistoProducerAlgoForTracker.h:161
MTVHistoProducerAlgoForTracker::h_STRIPlayersWith2dMeas_eta
std::vector< MonitorElement * > h_STRIPlayersWith2dMeas_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::dxyres_vs_eta
std::vector< MonitorElement * > dxyres_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:278
MTVHistoProducerAlgoForTracker::h_effic_vs_phi
std::vector< MonitorElement * > h_effic_vs_phi
Definition: MTVHistoProducerAlgoForTracker.h:157
MTVHistoProducerAlgoForTracker::totREC_vertcount_fwdpos
std::vector< std::vector< int > > totREC_vertcount_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:238
MTVHistoProducerAlgoForTracker::h_ptrmsh
std::vector< MonitorElement * > h_ptrmsh
Definition: MTVHistoProducerAlgoForTracker.h:281
MTVHistoProducerAlgoForTracker::totASS_vertpos
std::vector< std::vector< int > > totASS_vertpos
Definition: MTVHistoProducerAlgoForTracker.h:235
MTVHistoProducerAlgoForTracker::h_chi2meanhitsh
std::vector< MonitorElement * > h_chi2meanhitsh
Definition: MTVHistoProducerAlgoForTracker.h:275
MTVHistoProducerAlgoForTracker::h_fomt_ootpu_vertcount
std::vector< MonitorElement * > h_fomt_ootpu_vertcount
Definition: MTVHistoProducerAlgoForTracker.h:178
MTVHistoProducerAlgoForTracker::h_dxymeanhPt
std::vector< MonitorElement * > h_dxymeanhPt
Definition: MTVHistoProducerAlgoForTracker.h:287
MTVHistoProducerAlgoForTracker::h_dzpulleta
std::vector< MonitorElement * > h_dzpulleta
Definition: MTVHistoProducerAlgoForTracker.h:293
M_PI
#define M_PI
Definition: BFit3D.cc:3
MTVHistoProducerAlgoForTracker::h_misidpT
std::vector< MonitorElement * > h_misidpT
Definition: MTVHistoProducerAlgoForTracker.h:155
MTVHistoProducerAlgoForTracker::nPXLlayersWithMeas_vs_eta
std::vector< MonitorElement * > nPXLlayersWithMeas_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::totSIMpT
std::vector< std::vector< int > > totSIMpT
Definition: MTVHistoProducerAlgoForTracker.h:229
MTVHistoProducerAlgoForTracker::totASS_vertz_entire
std::vector< std::vector< int > > totASS_vertz_entire
Definition: MTVHistoProducerAlgoForTracker.h:240
MTVHistoProducerAlgoForTracker::totASS_vertcount_fwdneg
std::vector< std::vector< int > > totASS_vertcount_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:239
MTVHistoProducerAlgoForTracker::nTIBhits_vs_eta
std::vector< MonitorElement * > nTIBhits_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::h_fakerate_ootpu_fwdpos
std::vector< MonitorElement * > h_fakerate_ootpu_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:174
MTVHistoProducerAlgoForTracker::h_dxypulleta
std::vector< MonitorElement * > h_dxypulleta
Definition: MTVHistoProducerAlgoForTracker.h:293
MTVHistoProducerAlgoForTracker::etaintervals
std::vector< std::vector< double > > etaintervals
Definition: MTVHistoProducerAlgoForTracker.h:220
MTVHistoProducerAlgoForTracker::GpSelectorForEfficiencyVsCon
GenParticleCustomSelector * GpSelectorForEfficiencyVsCon
Definition: MTVHistoProducerAlgoForTracker.h:124
MTVHistoProducerAlgoForTracker::h_effic_vs_dz
std::vector< MonitorElement * > h_effic_vs_dz
Definition: MTVHistoProducerAlgoForTracker.h:159
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:13
MTVHistoProducerAlgoForTracker::h_phirmshPhi
std::vector< MonitorElement * > h_phirmshPhi
Definition: MTVHistoProducerAlgoForTracker.h:283
MTVHistoProducerAlgoForTracker::totREC_ootpu_fwdpos
std::vector< std::vector< int > > totREC_ootpu_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:247
reco::GsfTrack::qoverpModeError
double qoverpModeError() const
error on signed transverse curvature from mode
Definition: GsfTrack.h:78
reco::HitPattern::stripLayersWithMeasurement
int stripLayersWithMeasurement() const
Definition: HitPattern.h:731
MTVHistoProducerAlgoForTracker::totASS2_ootpu_entire
std::vector< std::vector< int > > totASS2_ootpu_entire
Definition: MTVHistoProducerAlgoForTracker.h:245
MTVHistoProducerAlgoForTracker::h_eta
std::vector< MonitorElement * > h_eta
Definition: MTVHistoProducerAlgoForTracker.h:162
FitSlicesYTool::getFittedMeanWithError
void getFittedMeanWithError(MonitorElement *)
Fill the ME with the mean value (with error) of the gaussian fit in each slice.
Definition: FitSlicesYTool.cc:59
MTVHistoProducerAlgoForTracker::h_effic_PU_vertcount
std::vector< MonitorElement * > h_effic_PU_vertcount
Definition: MTVHistoProducerAlgoForTracker.h:183
MTVHistoProducerAlgoForTracker::h_pt
std::vector< MonitorElement * > h_pt
Definition: MTVHistoProducerAlgoForTracker.h:162
MTVHistoProducerAlgoForTracker::h_PXBhits_eta
std::vector< MonitorElement * > h_PXBhits_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::totASS2_ootpu_eta_entire
std::vector< std::vector< int > > totASS2_ootpu_eta_entire
Definition: MTVHistoProducerAlgoForTracker.h:249
MTVHistoProducerAlgoForTracker::h_phimeanhPt
std::vector< MonitorElement * > h_phimeanhPt
Definition: MTVHistoProducerAlgoForTracker.h:287
MTVHistoProducerAlgoForTracker::h_dedx_nom2
std::vector< MonitorElement * > h_dedx_nom2
Definition: MTVHistoProducerAlgoForTracker.h:192
MTVHistoProducerAlgoForTracker::cotThetares_vs_pt
std::vector< MonitorElement * > cotThetares_vs_pt
Definition: MTVHistoProducerAlgoForTracker.h:279
MTVHistoProducerAlgoForTracker::h_effic_vs_zpos
std::vector< MonitorElement * > h_effic_vs_zpos
Definition: MTVHistoProducerAlgoForTracker.h:161
MTVHistoProducerAlgoForTracker::h_etaSIM
std::vector< MonitorElement * > h_etaSIM
Definition: MTVHistoProducerAlgoForTracker.h:150
reco::HitPattern::numberOfValidStripTIBHits
int numberOfValidStripTIBHits() const
Definition: HitPattern.h:598
MTVHistoProducerAlgoForTracker::h_simulphi
std::vector< MonitorElement * > h_simulphi
Definition: MTVHistoProducerAlgoForTracker.h:157
MTVHistoProducerAlgoForTracker::h_reco_vertcount_fwdpos
std::vector< MonitorElement * > h_reco_vertcount_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:165
reco::HitPattern::numberOfValidPixelEndcapHits
int numberOfValidPixelEndcapHits() const
Definition: HitPattern.h:590
MTVHistoProducerAlgoForTracker::h_pullTheta
std::vector< MonitorElement * > h_pullTheta
Definition: MTVHistoProducerAlgoForTracker.h:162
relval_parameters_module.step
string step
Definition: relval_parameters_module.py:36
MTVHistoProducerAlgoForTracker::totREC_hit
std::vector< std::vector< int > > totREC_hit
Definition: MTVHistoProducerAlgoForTracker.h:230
MTVHistoProducerAlgoForTracker::ptmean_vs_eta_phi
std::vector< MonitorElement * > ptmean_vs_eta_phi
Definition: MTVHistoProducerAlgoForTracker.h:266
MTVHistoProducerAlgoForTracker::h_looperphi
std::vector< MonitorElement * > h_looperphi
Definition: MTVHistoProducerAlgoForTracker.h:157
tmp
std::vector< std::vector< double > > tmp
Definition: MVATrainer.cc:100
reco::TrackBase::lambdaError
double lambdaError() const
error on lambda
Definition: TrackBase.h:203
MTVHistoProducerAlgoForTracker::totSIMeta
std::vector< std::vector< int > > totSIMeta
Definition: MTVHistoProducerAlgoForTracker.h:228
MTVHistoProducerAlgoForTracker::h_reco_vertcount_fwdneg
std::vector< MonitorElement * > h_reco_vertcount_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:166
MTVHistoProducerAlgoForTracker::totREC_ootpu_fwdneg
std::vector< std::vector< int > > totREC_ootpu_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:248
MTVHistoProducerAlgoForTracker::totSIM_vertcount_barrel
std::vector< std::vector< int > > totSIM_vertcount_barrel
Definition: MTVHistoProducerAlgoForTracker.h:237
MTVHistoProducerAlgoForTracker::dxypull_vs_eta
std::vector< MonitorElement * > dxypull_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:291
MTVHistoProducerAlgoForTracker::h_fomt_itpu_vertcount
std::vector< MonitorElement * > h_fomt_itpu_vertcount
Definition: MTVHistoProducerAlgoForTracker.h:177
MTVHistoProducerAlgoForTracker::ptres_vs_eta
std::vector< MonitorElement * > ptres_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:278
MTVHistoProducerAlgoForTracker::h_cotThetameanh
std::vector< MonitorElement * > h_cotThetameanh
Definition: MTVHistoProducerAlgoForTracker.h:286
MTVHistoProducerAlgoForTracker::h_assoc_vertz_entire
std::vector< MonitorElement * > h_assoc_vertz_entire
Definition: MTVHistoProducerAlgoForTracker.h:167
MTVHistoProducerAlgoForTracker::h_assoc2_vertcount_fwdpos
std::vector< MonitorElement * > h_assoc2_vertcount_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:165
MTVHistoProducerAlgoForTracker::h_dzrmsh
std::vector< MonitorElement * > h_dzrmsh
Definition: MTVHistoProducerAlgoForTracker.h:281
MTVHistoProducerAlgoForTracker::h_con_vertcount
std::vector< MonitorElement * > h_con_vertcount
Definition: MTVHistoProducerAlgoForTracker.h:184
MTVHistoProducerAlgoForTracker::generalGpSelector
GenParticleCustomSelector * generalGpSelector
Definition: MTVHistoProducerAlgoForTracker.h:122
MTVHistoProducerAlgoForTracker::h_effic_vs_dxy
std::vector< MonitorElement * > h_effic_vs_dxy
Definition: MTVHistoProducerAlgoForTracker.h:158
MTVHistoProducerAlgoForTracker::phipull_vs_eta
std::vector< MonitorElement * > phipull_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:291
reco::Track::found
unsigned short found() const
Number of valid hits on track.
Definition: Track.h:100
MTVHistoProducerAlgoForTracker::h_TOBhits_eta
std::vector< MonitorElement * > h_TOBhits_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::GpSelectorForEfficiencyVsPt
GenParticleCustomSelector * GpSelectorForEfficiencyVsPt
Definition: MTVHistoProducerAlgoForTracker.h:126
MTVHistoProducerAlgoForTracker::h_misidrate
std::vector< MonitorElement * > h_misidrate
Definition: MTVHistoProducerAlgoForTracker.h:154
MTVHistoProducerAlgoForTracker::h_TEChits_eta
std::vector< MonitorElement * > h_TEChits_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::h_reco_vertcount_barrel
std::vector< MonitorElement * > h_reco_vertcount_barrel
Definition: MTVHistoProducerAlgoForTracker.h:164
MTVHistoProducerAlgoForTracker::h_simulhit
std::vector< MonitorElement * > h_simulhit
Definition: MTVHistoProducerAlgoForTracker.h:156
MTVHistoProducerAlgoForTracker::h_effic_vertz_entire
std::vector< MonitorElement * > h_effic_vertz_entire
Definition: MTVHistoProducerAlgoForTracker.h:167
MTVHistoProducerAlgoForTracker::GpSelectorForEfficiencyVsVTXR
GenParticleCustomSelector * GpSelectorForEfficiencyVsVTXR
Definition: MTVHistoProducerAlgoForTracker.h:127
MTVHistoProducerAlgoForTracker::h_con_zpos
std::vector< MonitorElement * > h_con_zpos
Definition: MTVHistoProducerAlgoForTracker.h:184
MTVHistoProducerAlgoForTracker::h_dzmeanhPt
std::vector< MonitorElement * > h_dzmeanhPt
Definition: MTVHistoProducerAlgoForTracker.h:287
MTVHistoProducerAlgoForTracker::totASS_vertz_fwdneg
std::vector< std::vector< int > > totASS_vertz_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:243
MTVHistoProducerAlgoForTracker::h_recopT
std::vector< MonitorElement * > h_recopT
Definition: MTVHistoProducerAlgoForTracker.h:155
MTVHistoProducerAlgoForTracker::h_looperdz
std::vector< MonitorElement * > h_looperdz
Definition: MTVHistoProducerAlgoForTracker.h:159
MTVHistoProducerAlgoForTracker::h_simul_vertcount_fwdneg
std::vector< MonitorElement * > h_simul_vertcount_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:166
MTVHistoProducerAlgoForTracker::pTintervals
std::vector< std::vector< double > > pTintervals
Definition: MTVHistoProducerAlgoForTracker.h:221
MTVHistoProducerAlgoForTracker::h_assocdxy
std::vector< MonitorElement * > h_assocdxy
Definition: MTVHistoProducerAlgoForTracker.h:158
reco::GsfTrack::phiMode
double phiMode() const
azimuthal angle of momentum vector from mode
Definition: GsfTrack.h:56
TrackingParticle
Monte Carlo truth information used for tracking validation.
Definition: TrackingParticle.h:29
reco::TrackBase::charge
int charge() const
track electric charge
Definition: TrackBase.h:113
MTVHistoProducerAlgoForTracker::cotThetares_vs_eta
std::vector< MonitorElement * > cotThetares_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:278
MTVHistoProducerAlgoForTracker::h_assocFraction
std::vector< MonitorElement * > h_assocFraction
Definition: MTVHistoProducerAlgoForTracker.h:202
MTVHistoProducerAlgoForTracker::fill_generic_simTrack_histos
void fill_generic_simTrack_histos(int counter, const TrackingParticle::Vector &, const TrackingParticle::Point &vertex, int bx)
Definition: MTVHistoProducerAlgoForTracker.cc:818
MTVHistoProducerAlgoForTracker::h_loopratedz
std::vector< MonitorElement * > h_loopratedz
Definition: MTVHistoProducerAlgoForTracker.h:159
MTVHistoProducerAlgoForTracker::h_fomt_sig_eta
std::vector< MonitorElement * > h_fomt_sig_eta
Definition: MTVHistoProducerAlgoForTracker.h:176
MTVHistoProducerAlgoForTracker::h_cotThetarmshPt
std::vector< MonitorElement * > h_cotThetarmshPt
Definition: MTVHistoProducerAlgoForTracker.h:282
MTVHistoProducerAlgoForTracker::totSIM_vertcount_fwdneg
std::vector< std::vector< int > > totSIM_vertcount_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:239
MTVHistoProducerAlgoForTracker::totREC_phi
std::vector< std::vector< int > > totREC_phi
Definition: MTVHistoProducerAlgoForTracker.h:231
MTVHistoProducerAlgoForTracker::GpSelectorForEfficiencyVsEta
GenParticleCustomSelector * GpSelectorForEfficiencyVsEta
Definition: MTVHistoProducerAlgoForTracker.h:123
reco::GsfTrack::lambdaModeError
double lambdaModeError() const
error on lambda from mode
Definition: GsfTrack.h:89
MTVHistoProducerAlgoForTracker::totASS2_ootpu_barrel
std::vector< std::vector< int > > totASS2_ootpu_barrel
Definition: MTVHistoProducerAlgoForTracker.h:246
MTVHistoProducerAlgoForTracker::h_assoc_vertcount_barrel
std::vector< MonitorElement * > h_assoc_vertcount_barrel
Definition: MTVHistoProducerAlgoForTracker.h:164
MTVHistoProducerAlgoForTracker::h_simul_vertz_fwdpos
std::vector< MonitorElement * > h_simul_vertz_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:169
MTVHistoProducerAlgoForTracker::h_recodz
std::vector< MonitorElement * > h_recodz
Definition: MTVHistoProducerAlgoForTracker.h:159
MTVHistoProducerAlgoForTracker::ptres_vs_pt
std::vector< MonitorElement * > ptres_vs_pt
Definition: MTVHistoProducerAlgoForTracker.h:279
reco::modules::GenParticleCustomSelector
SingleObjectSelector< GenParticleCollection,::GenParticleCustomSelector > GenParticleCustomSelector
Definition: GenParticleCustomSelector.cc:16
MTVHistoProducerAlgoForTracker::h_reco_ootpu_fwdpos
std::vector< MonitorElement * > h_reco_ootpu_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:174
MTVHistoProducerAlgoForTracker::h_dedx_sat2
std::vector< MonitorElement * > h_dedx_sat2
Definition: MTVHistoProducerAlgoForTracker.h:194
TrackingParticle::Vector
math::XYZVectorD Vector
point in the space
Definition: TrackingParticle.h:37
MTVHistoProducerAlgoForTracker::h_reco_ootpu_vertcount
std::vector< MonitorElement * > h_reco_ootpu_vertcount
Definition: MTVHistoProducerAlgoForTracker.h:182
MTVHistoProducerAlgoForTracker::totloop_phi
std::vector< std::vector< int > > totloop_phi
Definition: MTVHistoProducerAlgoForTracker.h:231
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:850
MTVHistoProducerAlgoForTracker::h_assoc2_ootpu_entire
std::vector< MonitorElement * > h_assoc2_ootpu_entire
Definition: MTVHistoProducerAlgoForTracker.h:172
MTVHistoProducerAlgoForTracker::h_phipulleta
std::vector< MonitorElement * > h_phipulleta
Definition: MTVHistoProducerAlgoForTracker.h:293
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:121
MTVHistoProducerAlgoForTracker::h_simul_vertcount_entire
std::vector< MonitorElement * > h_simul_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:163
MTVHistoProducerAlgoForTracker::totCONvertcount
std::vector< std::vector< int > > totCONvertcount
Definition: MTVHistoProducerAlgoForTracker.h:252
MTVHistoProducerAlgoForTracker::h_assoceta
std::vector< MonitorElement * > h_assoceta
Definition: MTVHistoProducerAlgoForTracker.h:154
MTVHistoProducerAlgoForTracker::h_effic_vertcount_barrel
std::vector< MonitorElement * > h_effic_vertcount_barrel
Definition: MTVHistoProducerAlgoForTracker.h:164
MTVHistoProducerAlgoForTracker::h_looprate
std::vector< MonitorElement * > h_looprate
Definition: MTVHistoProducerAlgoForTracker.h:154
MTVHistoProducerAlgoForTracker::ptres_vs_phi
std::vector< MonitorElement * > ptres_vs_phi
Definition: MTVHistoProducerAlgoForTracker.h:263
MTVHistoProducerAlgoForTracker::TpSelectorForEfficiencyVsVTXR
TrackingParticleSelector * TpSelectorForEfficiencyVsVTXR
Definition: MTVHistoProducerAlgoForTracker.h:119
MTVHistoProducerAlgoForTracker::h_con_eta
std::vector< MonitorElement * > h_con_eta
Definition: MTVHistoProducerAlgoForTracker.h:184
MTVHistoProducerAlgoForTracker::h_assocSharedHit
std::vector< MonitorElement * > h_assocSharedHit
Definition: MTVHistoProducerAlgoForTracker.h:202
MTVHistoProducerAlgoForTracker::h_dxymeanh
std::vector< MonitorElement * > h_dxymeanh
Definition: MTVHistoProducerAlgoForTracker.h:286
MTVHistoProducerAlgoForTracker::nhits_vs_eta
std::vector< MonitorElement * > nhits_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:205
MTVHistoProducerAlgoForTracker::h_ptrmshPhi
std::vector< MonitorElement * > h_ptrmshPhi
Definition: MTVHistoProducerAlgoForTracker.h:283
reco::GsfTrack::ptMode
double ptMode() const
track transverse momentum from mode
Definition: GsfTrack.h:48
MTVHistoProducerAlgoForTracker::h_assoc2_vertcount_barrel
std::vector< MonitorElement * > h_assoc2_vertcount_barrel
Definition: MTVHistoProducerAlgoForTracker.h:164
MTVHistoProducerAlgoForTracker::h_effic_vs_hit
std::vector< MonitorElement * > h_effic_vs_hit
Definition: MTVHistoProducerAlgoForTracker.h:156
MTVHistoProducerAlgoForTracker::h_simuldxy
std::vector< MonitorElement * > h_simuldxy
Definition: MTVHistoProducerAlgoForTracker.h:158
MTVHistoProducerAlgoForTracker::totmisideta
std::vector< std::vector< int > > totmisideta
Definition: MTVHistoProducerAlgoForTracker.h:228
MTVHistoProducerAlgoForTracker::h_assocvertpos
std::vector< MonitorElement * > h_assocvertpos
Definition: MTVHistoProducerAlgoForTracker.h:161
MTVHistoProducerAlgoForTracker::h_loopratepT
std::vector< MonitorElement * > h_loopratepT
Definition: MTVHistoProducerAlgoForTracker.h:155
MTVHistoProducerAlgoForTracker::h_effic_vertz_fwdneg
std::vector< MonitorElement * > h_effic_vertz_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:170
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
MTVHistoProducerAlgoForTracker::totASS2_vertcount_fwdpos
std::vector< std::vector< int > > totASS2_vertcount_fwdpos
Definition: MTVHistoProducerAlgoForTracker.h:238
MTVHistoProducerAlgoForTracker::h_assoc2_vertcount_entire
std::vector< MonitorElement * > h_assoc2_vertcount_entire
Definition: MTVHistoProducerAlgoForTracker.h:163
MTVHistoProducerAlgoForTracker::MTVHistoProducerAlgoForTracker
MTVHistoProducerAlgoForTracker(const edm::ParameterSet &pset)
Definition: MTVHistoProducerAlgoForTracker.cc:18
MTVHistoProducerAlgoForTracker::h_simul_vertz_fwdneg
std::vector< MonitorElement * > h_simul_vertz_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:170
MTVHistoProducerAlgoForTracker::h_recoeta
std::vector< MonitorElement * > h_recoeta
Definition: MTVHistoProducerAlgoForTracker.h:154
MTVHistoProducerAlgoForTracker::thetapull_vs_eta
std::vector< MonitorElement * > thetapull_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:291
MTVHistoProducerAlgoForTracker::totSIM_dxy
std::vector< std::vector< int > > totSIM_dxy
Definition: MTVHistoProducerAlgoForTracker.h:232
MTVHistoProducerAlgoForTracker::chi2_vs_eta
std::vector< MonitorElement * > chi2_vs_eta
Definition: MTVHistoProducerAlgoForTracker.h:272
MTVHistoProducerAlgoForTracker::h_phipulletamean
std::vector< MonitorElement * > h_phipulletamean
Definition: MTVHistoProducerAlgoForTracker.h:295
MTVHistoProducerAlgoForTracker::h_PXLlayersWithMeas_eta
std::vector< MonitorElement * > h_PXLlayersWithMeas_eta
Definition: MTVHistoProducerAlgoForTracker.h:213
MTVHistoProducerAlgoForTracker::dzintervals
std::vector< std::vector< double > > dzintervals
Definition: MTVHistoProducerAlgoForTracker.h:224
MTVHistoProducerAlgoForTracker::dxyres_vs_pt
std::vector< MonitorElement * > dxyres_vs_pt
Definition: MTVHistoProducerAlgoForTracker.h:279
MTVHistoProducerAlgoForTracker::TpSelectorForEfficiencyVsCon
TrackingParticleSelector * TpSelectorForEfficiencyVsCon
Definition: MTVHistoProducerAlgoForTracker.h:116
MTVHistoProducerAlgoForTracker::h_misidratedxy
std::vector< MonitorElement * > h_misidratedxy
Definition: MTVHistoProducerAlgoForTracker.h:158
MTVHistoProducerAlgoForTracker::h_dedx_nom1
std::vector< MonitorElement * > h_dedx_nom1
Definition: MTVHistoProducerAlgoForTracker.h:191
MTVHistoProducerAlgoForTracker::totASS2_dxy
std::vector< std::vector< int > > totASS2_dxy
Definition: MTVHistoProducerAlgoForTracker.h:232
MTVHistoProducerAlgoForTracker::h_losthits
std::vector< MonitorElement * > h_losthits
Definition: MTVHistoProducerAlgoForTracker.h:258
MTVHistoProducerAlgoForTracker::TpSelectorForEfficiencyVsEta
TrackingParticleSelector * TpSelectorForEfficiencyVsEta
Definition: MTVHistoProducerAlgoForTracker.h:115
MTVHistoProducerAlgoForTracker::h_assoc2eta
std::vector< MonitorElement * > h_assoc2eta
Definition: MTVHistoProducerAlgoForTracker.h:154
MTVHistoProducerAlgoForTracker::totmisid_dz
std::vector< std::vector< int > > totmisid_dz
Definition: MTVHistoProducerAlgoForTracker.h:233
MTVHistoProducerAlgoForTracker::h_thetapulletamean
std::vector< MonitorElement * > h_thetapulletamean
Definition: MTVHistoProducerAlgoForTracker.h:295
MTVHistoProducerAlgoForTracker::h_misidphi
std::vector< MonitorElement * > h_misidphi
Definition: MTVHistoProducerAlgoForTracker.h:157
MTVHistoProducerAlgoForTracker::h_fakerate_vertcount_fwdneg
std::vector< MonitorElement * > h_fakerate_vertcount_fwdneg
Definition: MTVHistoProducerAlgoForTracker.h:166
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:1180
MTVHistoProducerAlgoForTracker::dzres_vs_pt
std::vector< MonitorElement * > dzres_vs_pt
Definition: MTVHistoProducerAlgoForTracker.h:279
Generated for CMSSW Reference Manual by doxygen 1.8.5