CMS 3D CMS Logo
PrimaryVertexValidation.h
Go to the documentation of this file.
1 #ifndef PrimaryVertexValidation_h
2 #define PrimaryVertexValidation_h
3 
4 // system include files
5 #include <string>
6 #include <sstream>
7 #include <vector>
8 #include <map>
9 
10 // ROOT Included
11 #include "TFile.h"
12 #include "TH1D.h"
13 #include "TH1I.h"
14 #include "TH2D.h"
15 #include "TTree.h"
16 
17 // CMSSW includes
18 #include "CommonTools/UtilAlgos/interface/TFileService.h"
19 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
20 #include "DataFormats/GeometryVector/interface/LocalPoint.h"
21 #include "DataFormats/TrackCandidate/interface/TrackCandidate.h"
22 #include "DataFormats/TrackCandidate/interface/TrackCandidateCollection.h"
23 #include "DataFormats/TrackReco/interface/Track.h"
24 #include "DataFormats/TrackReco/interface/TrackFwd.h"
25 #include "DataFormats/TrackerRecHit2D/interface/ProjectedSiStripRecHit2D.h"
26 #include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2DCollection.h"
27 #include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2D.h"
28 #include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2DCollection.h"
29 #include "DataFormats/TrackingRecHit/interface/TrackingRecHit.h"
30 #include "DataFormats/TrajectorySeed/interface/TrajectorySeed.h"
31 #include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"
32 #include "DataFormats/VertexReco/interface/Vertex.h"
33 #include "DataFormats/VertexReco/interface/VertexFwd.h"
34 #include "FWCore/Framework/interface/ConsumesCollector.h"
35 #include "FWCore/Framework/interface/one/EDAnalyzer.h"
36 #include "FWCore/Framework/interface/ESHandle.h"
37 #include "FWCore/Framework/interface/Event.h"
38 #include "FWCore/Framework/interface/EventSetup.h"
39 #include "FWCore/Framework/interface/Frameworkfwd.h"
40 #include "FWCore/Framework/interface/MakerMacros.h"
41 #include "FWCore/ParameterSet/interface/ParameterSet.h"
42 #include "FWCore/Utilities/interface/EDGetToken.h"
43 #include "Geometry/CommonDetUnit/interface/GeomDet.h"
44 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
45 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
46 #include "RecoVertex/AdaptiveVertexFit/interface/AdaptiveVertexFitter.h"
47 #include "RecoVertex/PrimaryVertexProducer/interface/PrimaryVertexProducerAlgorithm.h"
48 #include "RecoVertex/PrimaryVertexProducer/interface/TrackClusterizerInZ.h"
49 #include "RecoVertex/PrimaryVertexProducer/interface/TrackFilterForPVFinding.h"
50 #include "RecoVertex/PrimaryVertexProducer/interface/TrackFilterForPVFindingBase.h"
51 #include "RecoVertex/VertexPrimitives/interface/VertexFitter.h"
52 #include "TrackingTools/IPTools/interface/IPTools.h"
53 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateClosestToPoint.h"
54 #include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
55 
56 //
57 // ancyllary enum for
58 // residuals moments estimation
59 //
60 
61 namespace statmode{
62  enum estimator
63  { MEAN = 1,
64  WIDTH = 2,
65  MEDIAN = 3,
66  MAD = 4,
67  UNKWN = -1
68  };
69 }
70 
71 //
72 // class decleration
73 //
74 
75 class PrimaryVertexValidation : public edm::one::EDAnalyzer<edm::one::SharedResources> {
76 
77  public:
78  explicit PrimaryVertexValidation(const edm::ParameterSet&);
79  ~PrimaryVertexValidation();
80 
81  private:
82  virtual void beginJob();
83  virtual void analyze(const edm::Event&, const edm::EventSetup&);
84  virtual void endJob();
85  bool isBFieldConsistentWithMode(const edm::EventSetup& iSetup) const;
86  bool isHit2D(const TrackingRecHit &hit) const;
87  bool hasFirstLayerPixelHits(const reco::TransientTrack track);
88  std::pair<bool,bool> pixelHitsCheck(const reco::TransientTrack track);
89  std::pair<Double_t,Double_t> getMedian(TH1F *histo);
90  std::pair<Double_t,Double_t> getMAD(TH1F *histo);
91  std::pair<std::pair<Double_t,Double_t>, std::pair<Double_t,Double_t> > fitResiduals(TH1 *hist);
92  void fillTrendPlot(TH1F* trendPlot, TH1F *residualsPlot[100], statmode::estimator fitPar_, TString var_);
93  void fillTrendPlotByIndex(TH1F* trendPlot,std::vector<TH1F*>& h, statmode::estimator fitPar_);
94 
95  static bool vtxSort( const reco::Vertex & a, const reco::Vertex & b );
96  bool passesTrackCuts(const reco::Track & track, const reco::Vertex & vertex,std::string qualityString_, double dxyErrMax_,double dzErrMax_, double ptErrMax_);
97 
98  std::vector<TH1F*> bookResidualsHistogram(TFileDirectory dir,unsigned int theNOfBins,TString resType,TString varType);
99  std::map<std::string, TH1*> bookVertexHistograms(TFileDirectory dir);
100  void fillTrackHistos(std::map<std::string, TH1*> & h, const std::string & ttype, const reco::TransientTrack *tt, const reco::Vertex & v,const reco::BeamSpot & beamSpot, double fBfield);
101  void add(std::map<std::string, TH1*>& h, TH1* hist);
102  void fill(std::map<std::string, TH1*>& h, std::string s, double x);
103  void fill(std::map<std::string, TH1*>& h, std::string s, double x, double y);
104  void fillByIndex(std::vector<TH1F*>& h, unsigned int index, double x);
105  void fillMap(TH2F* trendMap, TH1F* residualsMapPlot[100][100], statmode::estimator fitPar_);
106 
107  inline double square(double x){
108  return x*x;
109  }
110 
111  // ----------member data ---------------------------
112  edm::ParameterSet theConfig;
113  int Nevt_;
114 
115  TrackFilterForPVFindingBase* theTrackFilter_;
116  TrackClusterizerInZ* theTrackClusterizer_;
117 
118  // setting of the number of plots
119  static const int nMaxBins_ = 100; // maximum number of bookable histograms
120 
121  // Output
122  bool storeNtuple_;
123  bool lightNtupleSwitch_; // switch to keep only info for daily validation
124  bool useTracksFromRecoVtx_;
125 
126  // requirements on the vertex
127  double vertexZMax_;
128 
129  // integrated lumi (if info available)
130  double intLumi_;
131 
132  // requirements on the probe
133  bool askFirstLayerHit_; // ask hit in the first layer of pixels
134  bool doBPix_;
135  bool doFPix_;
136  double ptOfProbe_;
137  double etaOfProbe_;
138  bool isPhase1_;
139  int nBins_; // actual number of histograms
140  std::vector<unsigned int> runControlNumbers_;
141 
142  bool debug_;
143  bool runControl_;
144 
145  edm::EDGetTokenT<reco::TrackCollection> theTrackCollectionToken;
146  edm::EDGetTokenT<reco::VertexCollection> theVertexCollectionToken;
147  edm::EDGetTokenT<reco::BeamSpot> theBeamspotToken;
148 
149  TTree* rootTree_;
150 
151  // Root-Tuple variables :
152  //=======================
153  void SetVarToZero();
154 
155  static const int nMaxtracks_ = 1000;
156  static const int cmToum = 10000;
157  static const int nPtBins_ = 48;
158 
159  float phiSect_;
160  float etaSect_;
161 
162  // pT binning as in paragraph 3.2 of CMS-PAS-TRK-10-005 (https://cds.cern.ch/record/1279383/files/TRK-10-005-pas.pdf)
163 
164  // 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
165  const float mypT_bins_[nPtBins_+1] = {0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3.0,3.1,3.2,3.3,3.4,3.5,3.6,3.7,3.8,3.9,4.0,4.1,4.25,4.5,4.75,5.0,5.5,6.0,7.0,8.0,9.0,11.0,14.0,20.};
166 
167  // event-related quantities
168  int nTracks_;
169  int nTracksPerClus_;
170  int nClus_;
171  int nOfflineVertices_;
172  unsigned int RunNumber_;
173  unsigned int EventNumber_;
174  unsigned int LuminosityBlockNumber_;
175  double xOfflineVertex_;
176  double yOfflineVertex_;
177  double zOfflineVertex_;
178 
179  double xErrOfflineVertex_;
180  double yErrOfflineVertex_;
181  double zErrOfflineVertex_;
182 
183  double BSx0_;
184  double BSy0_;
185  double BSz0_;
186  double Beamsigmaz_;
187  double Beamdxdz_;
188  double BeamWidthX_;
189  double BeamWidthY_;
190  double wxy2_;
191 
192  // track-related quantities
193  double pt_[nMaxtracks_];
194  double p_[nMaxtracks_];
195  int nhits_[nMaxtracks_];
196  int nhits1D_[nMaxtracks_];
197  int nhits2D_[nMaxtracks_];
198  int nhitsBPIX_[nMaxtracks_];
199  int nhitsFPIX_[nMaxtracks_];
200  int nhitsTIB_[nMaxtracks_];
201  int nhitsTID_[nMaxtracks_];
202  int nhitsTOB_[nMaxtracks_];
203  int nhitsTEC_[nMaxtracks_];
204  int isHighPurity_[nMaxtracks_];
205  double eta_[nMaxtracks_];
206  double theta_[nMaxtracks_];
207  double phi_[nMaxtracks_];
208  double chi2_[nMaxtracks_];
209  double chi2ndof_[nMaxtracks_];
210  int charge_[nMaxtracks_];
211  double qoverp_[nMaxtracks_];
212  double dz_[nMaxtracks_];
213  double dxy_[nMaxtracks_];
214  double dxyBs_[nMaxtracks_];
215  double dzBs_[nMaxtracks_];
216  double xPCA_[nMaxtracks_];
217  double yPCA_[nMaxtracks_];
218  double zPCA_[nMaxtracks_];
219  double xUnbiasedVertex_[nMaxtracks_];
220  double yUnbiasedVertex_[nMaxtracks_];
221  double zUnbiasedVertex_[nMaxtracks_];
222  float chi2normUnbiasedVertex_[nMaxtracks_];
223  float chi2UnbiasedVertex_[nMaxtracks_];
224  float chi2ProbUnbiasedVertex_[nMaxtracks_];
225  float DOFUnbiasedVertex_[nMaxtracks_];
226  float sumOfWeightsUnbiasedVertex_[nMaxtracks_];
227  int tracksUsedForVertexing_[nMaxtracks_];
228 
229  double dxyFromMyVertex_[nMaxtracks_];
230  double dzFromMyVertex_[nMaxtracks_];
231  double d3DFromMyVertex_[nMaxtracks_];
232 
233  double dxyErrorFromMyVertex_[nMaxtracks_];
234  double dzErrorFromMyVertex_[nMaxtracks_];
235  double d3DErrorFromMyVertex_[nMaxtracks_];
236 
237  double IPTsigFromMyVertex_[nMaxtracks_];
238  double IPLsigFromMyVertex_[nMaxtracks_];
239  double IP3DsigFromMyVertex_[nMaxtracks_];
240 
241  int hasRecVertex_[nMaxtracks_];
242  int isGoodTrack_[nMaxtracks_];
243 
244  // histogram for max(eta)
245  TH1F* h_etaMax;
246  TH1F* h_nbins;
247 
248  // ---- directly histograms // ===> unbiased residuals
249 
250  // absolute residuals
251 
252  TH1F* a_dxyPhiResiduals[nMaxBins_];
253  TH1F* a_dxyEtaResiduals[nMaxBins_];
254 
255  TH1F* a_dxPhiResiduals[nMaxBins_];
256  TH1F* a_dxEtaResiduals[nMaxBins_];
257 
258  TH1F* a_dyPhiResiduals[nMaxBins_];
259  TH1F* a_dyEtaResiduals[nMaxBins_];
260 
261  TH1F* a_dzPhiResiduals[nMaxBins_];
262  TH1F* a_dzEtaResiduals[nMaxBins_];
263 
264  TH1F* a_IP2DPhiResiduals[nMaxBins_];
265  TH1F* a_IP2DEtaResiduals[nMaxBins_];
266 
267  TH1F* a_IP3DPhiResiduals[nMaxBins_];
268  TH1F* a_IP3DEtaResiduals[nMaxBins_];
269 
270  TH1F* a_reszPhiResiduals[nMaxBins_];
271  TH1F* a_reszEtaResiduals[nMaxBins_];
272 
273  TH1F* a_d3DPhiResiduals[nMaxBins_];
274  TH1F* a_d3DEtaResiduals[nMaxBins_];
275 
276  // normalized residuals
277 
278  TH1F* n_dxyPhiResiduals[nMaxBins_];
279  TH1F* n_dxyEtaResiduals[nMaxBins_];
280 
281  TH1F* n_dzPhiResiduals[nMaxBins_];
282  TH1F* n_dzEtaResiduals[nMaxBins_];
283 
284  TH1F* n_IP2DPhiResiduals[nMaxBins_];
285  TH1F* n_IP2DEtaResiduals[nMaxBins_];
286 
287  TH1F* n_IP3DPhiResiduals[nMaxBins_];
288  TH1F* n_IP3DEtaResiduals[nMaxBins_];
289 
290  TH1F* n_reszPhiResiduals[nMaxBins_];
291  TH1F* n_reszEtaResiduals[nMaxBins_];
292 
293  TH1F* n_d3DPhiResiduals[nMaxBins_];
294  TH1F* n_d3DEtaResiduals[nMaxBins_];
295 
296  // for the maps
297 
298  TH1F* a_dxyResidualsMap[nMaxBins_][nMaxBins_];
299  TH1F* a_dzResidualsMap[nMaxBins_][nMaxBins_];
300  TH1F* a_d3DResidualsMap[nMaxBins_][nMaxBins_];
301 
302  TH1F* n_dxyResidualsMap[nMaxBins_][nMaxBins_];
303  TH1F* n_dzResidualsMap[nMaxBins_][nMaxBins_];
304  TH1F* n_d3DResidualsMap[nMaxBins_][nMaxBins_];
305 
306  // ---- trends as function of phi and eta
307 
308  TH1F* a_dxyPhiMeanTrend;
309  TH1F* a_dxyPhiWidthTrend;
310  TH1F* a_dzPhiMeanTrend;
311  TH1F* a_dzPhiWidthTrend;
312 
313  TH1F* a_dxyEtaMeanTrend;
314  TH1F* a_dxyEtaWidthTrend;
315  TH1F* a_dzEtaMeanTrend;
316  TH1F* a_dzEtaWidthTrend;
317 
318  TH1F* n_dxyPhiMeanTrend;
319  TH1F* n_dxyPhiWidthTrend;
320  TH1F* n_dzPhiMeanTrend;
321  TH1F* n_dzPhiWidthTrend;
322 
323  TH1F* n_dxyEtaMeanTrend;
324  TH1F* n_dxyEtaWidthTrend;
325  TH1F* n_dzEtaMeanTrend;
326  TH1F* n_dzEtaWidthTrend;
327 
328  // ---- trends as a function of pT
329 
330  TH1F* a_dxypTMeanTrend;
331  TH1F* a_dxypTWidthTrend;
332  TH1F* a_dzpTMeanTrend;
333  TH1F* a_dzpTWidthTrend;
334 
335  TH1F* a_dxypTCentralMeanTrend;
336  TH1F* a_dxypTCentralWidthTrend;
337  TH1F* a_dzpTCentralMeanTrend;
338  TH1F* a_dzpTCentralWidthTrend;
339 
340  TH1F* n_dxypTMeanTrend;
341  TH1F* n_dxypTWidthTrend;
342  TH1F* n_dzpTMeanTrend;
343  TH1F* n_dzpTWidthTrend;
344 
345  TH1F* n_dxypTCentralMeanTrend;
346  TH1F* n_dxypTCentralWidthTrend;
347  TH1F* n_dzpTCentralMeanTrend;
348  TH1F* n_dzpTCentralWidthTrend;
349 
350  // ---- medians and MAD
351 
352  TH1F* a_dxyPhiMedianTrend;
353  TH1F* a_dxyPhiMADTrend;
354  TH1F* a_dzPhiMedianTrend;
355  TH1F* a_dzPhiMADTrend;
356 
357  TH1F* a_dxyEtaMedianTrend;
358  TH1F* a_dxyEtaMADTrend;
359  TH1F* a_dzEtaMedianTrend;
360  TH1F* a_dzEtaMADTrend;
361 
362  TH1F* n_dxyPhiMedianTrend;
363  TH1F* n_dxyPhiMADTrend;
364  TH1F* n_dzPhiMedianTrend;
365  TH1F* n_dzPhiMADTrend;
366 
367  TH1F* n_dxyEtaMedianTrend;
368  TH1F* n_dxyEtaMADTrend;
369  TH1F* n_dzEtaMedianTrend;
370  TH1F* n_dzEtaMADTrend;
371 
372  // 2D residuals
373 
374  TH2F* a_dxyVsPhi;
375  TH2F* a_dzVsPhi;
376 
377  TH2F* n_dxyVsPhi;
378  TH2F* n_dzVsPhi;
379 
380  TH2F* a_dxyVsEta;
381  TH2F* a_dzVsEta;
382 
383  TH2F* n_dxyVsEta;
384  TH2F* n_dzVsEta;
385 
386  // 2D maps
387 
388  TH2F* a_dxyMeanMap;
389  TH2F* a_dzMeanMap;
390 
391  TH2F* n_dxyMeanMap;
392  TH2F* n_dzMeanMap;
393 
394  TH2F* a_dxyWidthMap;
395  TH2F* a_dzWidthMap;
396 
397  TH2F* n_dxyWidthMap;
398  TH2F* n_dzWidthMap;
399 
400  // ---- directly histograms =================> biased residuals
401 
402  // absolute residuals
403 
404  TH1F* a_dxyPhiBiasResiduals[nMaxBins_];
405  TH1F* a_dxyEtaBiasResiduals[nMaxBins_];
406 
407  TH1F* a_dzPhiBiasResiduals[nMaxBins_];
408  TH1F* a_dzEtaBiasResiduals[nMaxBins_];
409 
410  // normalized BiasResiduals
411 
412  TH1F* n_dxyPhiBiasResiduals[nMaxBins_];
413  TH1F* n_dxyEtaBiasResiduals[nMaxBins_];
414 
415  TH1F* n_dzPhiBiasResiduals[nMaxBins_];
416  TH1F* n_dzEtaBiasResiduals[nMaxBins_];
417 
418  // for the maps
419 
420  TH1F* a_dxyBiasResidualsMap[nMaxBins_][nMaxBins_];
421  TH1F* a_dzBiasResidualsMap[nMaxBins_][nMaxBins_];
422 
423  TH1F* n_dxyBiasResidualsMap[nMaxBins_][nMaxBins_];
424  TH1F* n_dzBiasResidualsMap[nMaxBins_][nMaxBins_];
425 
426  // ---- trends as function of phi
427 
428  TH1F* a_dxyPhiMeanBiasTrend;
429  TH1F* a_dxyPhiWidthBiasTrend;
430  TH1F* a_dzPhiMeanBiasTrend;
431  TH1F* a_dzPhiWidthBiasTrend;
432 
433  TH1F* a_dxyEtaMeanBiasTrend;
434  TH1F* a_dxyEtaWidthBiasTrend;
435  TH1F* a_dzEtaMeanBiasTrend;
436  TH1F* a_dzEtaWidthBiasTrend;
437 
438  TH1F* n_dxyPhiMeanBiasTrend;
439  TH1F* n_dxyPhiWidthBiasTrend;
440  TH1F* n_dzPhiMeanBiasTrend;
441  TH1F* n_dzPhiWidthBiasTrend;
442 
443  TH1F* n_dxyEtaMeanBiasTrend;
444  TH1F* n_dxyEtaWidthBiasTrend;
445  TH1F* n_dzEtaMeanBiasTrend;
446  TH1F* n_dzEtaWidthBiasTrend;
447 
448  // ---- medians and MAD
449 
450  TH1F* a_dxyPhiMedianBiasTrend;
451  TH1F* a_dxyPhiMADBiasTrend;
452  TH1F* a_dzPhiMedianBiasTrend;
453  TH1F* a_dzPhiMADBiasTrend;
454 
455  TH1F* a_dxyEtaMedianBiasTrend;
456  TH1F* a_dxyEtaMADBiasTrend;
457  TH1F* a_dzEtaMedianBiasTrend;
458  TH1F* a_dzEtaMADBiasTrend;
459 
460  TH1F* n_dxyPhiMedianBiasTrend;
461  TH1F* n_dxyPhiMADBiasTrend;
462  TH1F* n_dzPhiMedianBiasTrend;
463  TH1F* n_dzPhiMADBiasTrend;
464 
465  TH1F* n_dxyEtaMedianBiasTrend;
466  TH1F* n_dxyEtaMADBiasTrend;
467  TH1F* n_dzEtaMedianBiasTrend;
468  TH1F* n_dzEtaMADBiasTrend;
469 
470  // 2D maps
471 
472  TH2F* a_dxyMeanBiasMap;
473  TH2F* a_dzMeanBiasMap;
474 
475  TH2F* n_dxyMeanBiasMap;
476  TH2F* n_dzMeanBiasMap;
477 
478  TH2F* a_dxyWidthBiasMap;
479  TH2F* a_dzWidthBiasMap;
480 
481  TH2F* n_dxyWidthBiasMap;
482  TH2F* n_dzWidthBiasMap;
483 
484  // check event
485  TH1F* h_nTracks;
486  TH1F* h_nClus;
487  TH1F* h_nOfflineVertices;
488  TH1F* h_runNumber;
489  TH1F* h_lumiFromConfig;
490  TH1I* h_runFromConfig;
491  TH1I* h_runFromEvent;
492  TH1F* h_xOfflineVertex;
493  TH1F* h_yOfflineVertex;
494  TH1F* h_zOfflineVertex;
495  TH1F* h_xErrOfflineVertex;
496  TH1F* h_yErrOfflineVertex;
497  TH1F* h_zErrOfflineVertex;
498  TH1F* h_BSx0;
499  TH1F* h_BSy0;
500  TH1F* h_BSz0;
501  TH1F* h_Beamsigmaz;
502  TH1F* h_BeamWidthX;
503  TH1F* h_BeamWidthY;
504 
505  // check probe
506 
507  TH2F* h2_probeEtaPhi_;
508  TH2F* h2_probeEtaPt_;
509 
510  TH1F* h_probeP_;
511  TH1F* h_probePt_;
512  TH1F* h_probeEta_;
513  TH1F* h_probePhi_;
514  TH1F* h_probeChi2_;
515  TH1F* h_probeNormChi2_;
516  TH1F* h_probeCharge_;
517  TH1F* h_probeQoverP_;
518 
519  TH1F* h_probedzRecoV_;
520  TH1F* h_probedxyRecoV_;
521 
522  TH1F* h_probedzRefitV_;
523  TH1F* h_probedxyRefitV_;
524 
525  TH1F* h_probed0RefitV_;
526  TH1F* h_probez0RefitV_;
527 
528  TH1F* h_probesignIP2DRefitV_;
529  TH1F* h_probed3DRefitV_;
530  TH1F* h_probereszRefitV_;
531 
532  TH1F* h_probeRecoVSigZ_;
533  TH1F* h_probeRecoVSigXY_;
534  TH1F* h_probeRefitVSigZ_;
535  TH1F* h_probeRefitVSigXY_;
536  TH1F* h_probeRefitVSig3D_;
537  TH1F* h_probeRefitVLogSig3D_;
538  TH1F* h_probeRefitVSigResZ_;
539 
540  TH1F* h_probeHits_;
541  TH1F* h_probeHits1D_;
542  TH1F* h_probeHits2D_;
543  TH1F* h_probeHitsInTIB_;
544  TH1F* h_probeHitsInTOB_;
545  TH1F* h_probeHitsInTID_;
546  TH1F* h_probeHitsInTEC_;
547  TH1F* h_probeHitsInBPIX_;
548  TH1F* h_probeHitsInFPIX_;
549 
550  // check vertex
551 
552  TH1F* h_fitVtxNdof_;
553  TH1F* h_fitVtxChi2_;
554  TH1F* h_fitVtxNtracks_;
555  TH1F* h_fitVtxChi2ndf_;
556  TH1F* h_fitVtxChi2Prob_;
557  TH1F* h_fitVtxTrackWeights_;
558  TH1F* h_fitVtxTrackAverageWeight_;
559 
560  TH1F* h_recoVtxNtracks_;
561  TH1F* h_recoVtxChi2ndf_;
562  TH1F* h_recoVtxChi2Prob_;
563  TH1F* h_recoVtxSumPt_;
564 
565  std::map<std::string, TH1*> hDA;
566 
567  // histograms for the plots as function of pT
568 
569  std::vector<TH1F*> h_dxy_pT_;
570  std::vector<TH1F*> h_dz_pT_;
571  std::vector<TH1F*> h_norm_dxy_pT_;
572  std::vector<TH1F*> h_norm_dz_pT_;
573 
574  std::vector<TH1F*> h_dxy_Central_pT_;
575  std::vector<TH1F*> h_dz_Central_pT_;
576  std::vector<TH1F*> h_norm_dxy_Central_pT_;
577  std::vector<TH1F*> h_norm_dz_Central_pT_;
578 
579 };
580 
581 #endif
PrimaryVertexValidation::runControlNumbers_
std::vector< unsigned int > runControlNumbers_
Definition: PrimaryVertexValidation.h:140
Vispa.Share.Profiling.analyze
def analyze(function, filename, filter=None)
Definition: Profiling.py:11
PrimaryVertexValidation::h_dxy_pT_
std::vector< TH1F * > h_dxy_pT_
Definition: PrimaryVertexValidation.h:569
PrimaryVertexValidation::theTrackFilter_
TrackFilterForPVFindingBase * theTrackFilter_
Definition: PrimaryVertexValidation.h:115
PrimaryVertexValidation::hDA
std::map< std::string, TH1 * > hDA
Definition: PrimaryVertexValidation.h:565
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
PrimaryVertexValidation::h_norm_dz_Central_pT_
std::vector< TH1F * > h_norm_dz_Central_pT_
Definition: PrimaryVertexValidation.h:577
Clusterizer1DCommons::add
void add(const std::vector< const T * > &source, std::vector< const T * > &dest)
Definition: Clusterizer1DCommons.h:26
ecalDrivenElectronSeedsParameters_cff.beamSpot
beamSpot
Definition: ecalDrivenElectronSeedsParameters_cff.py:11
bk::beginJob
void beginJob()
Definition: Breakpoints.cc:15
x
T x() const
Cartesian x coordinate.
Definition: Basic3DVectorLD.h:127
PrimaryVertexValidation::h_norm_dxy_pT_
std::vector< TH1F * > h_norm_dxy_pT_
Definition: PrimaryVertexValidation.h:571
PrimaryVertexValidation::theVertexCollectionToken
edm::EDGetTokenT< reco::VertexCollection > theVertexCollectionToken
Definition: PrimaryVertexValidation.h:146
PrimaryVertexValidation::h_norm_dz_pT_
std::vector< TH1F * > h_norm_dz_pT_
Definition: PrimaryVertexValidation.h:572
PrimaryVertexValidation::h_dz_pT_
std::vector< TH1F * > h_dz_pT_
Definition: PrimaryVertexValidation.h:570
PrimaryVertexValidation::h_dz_Central_pT_
std::vector< TH1F * > h_dz_Central_pT_
Definition: PrimaryVertexValidation.h:575
HiIsolationCommonParameters_cff.track
track
Definition: HiIsolationCommonParameters_cff.py:9
PrimaryVertexValidation::theTrackClusterizer_
TrackClusterizerInZ * theTrackClusterizer_
Definition: PrimaryVertexValidation.h:116
b
double b
Definition: hdecay.h:120
PrimaryVertexValidation::h_dxy_Central_pT_
std::vector< TH1F * > h_dxy_Central_pT_
Definition: PrimaryVertexValidation.h:574
PrimaryVertexValidation::theBeamspotToken
edm::EDGetTokenT< reco::BeamSpot > theBeamspotToken
Definition: PrimaryVertexValidation.h:147
PrimaryVertexValidation::theTrackCollectionToken
edm::EDGetTokenT< reco::TrackCollection > theTrackCollectionToken
Definition: PrimaryVertexValidation.h:145
a
double a
Definition: hdecay.h:121
dir
dbl *** dir
Definition: mlp_gen.cc:35
create_public_lumi_plots.hist
hist
Definition: create_public_lumi_plots.py:34
PrimaryVertexValidation::h_norm_dxy_Central_pT_
std::vector< TH1F * > h_norm_dxy_Central_pT_
Definition: PrimaryVertexValidation.h:576
Generated for CMSSW Reference Manual by doxygen 1.8.11