CMS 3D CMS Logo
StandaloneTrackMonitor.cc
Go to the documentation of this file.
1 // system includes
2 #include <map>
3 #include <set>
4 #include <string>
5 #include <vector>
6 
7 // user includes
8 #include "DQMServices/Core/interface/DQMEDAnalyzer.h"
9 #include "DQMServices/Core/interface/DQMStore.h"
10 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
11 #include "DataFormats/Common/interface/DetSetVectorNew.h"
12 #include "DataFormats/Common/interface/Handle.h"
13 #include "DataFormats/TrackCandidate/interface/TrajectoryStopReasons.h"
14 #include "DataFormats/DetId/interface/DetId.h"
15 #include "DataFormats/JetReco/interface/PFJet.h"
16 #include "DataFormats/TrackReco/interface/Track.h"
17 #include "DataFormats/TrackReco/interface/TrackFwd.h"
18 #include "DataFormats/TrackerRecHit2D/interface/BaseTrackerRecHit.h"
19 #include "DataFormats/TrackerRecHit2D/interface/OmniClusterRef.h"
20 #include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2D.h"
21 #include "DataFormats/VertexReco/interface/Vertex.h"
22 #include "DataFormats/VertexReco/interface/VertexFwd.h"
23 #include "DataFormats/Math/interface/deltaR.h"
24 #include "FWCore/Framework/interface/ESHandle.h"
25 #include "FWCore/Framework/interface/Event.h"
26 #include "FWCore/Framework/interface/Frameworkfwd.h"
27 #include "FWCore/Framework/interface/MakerMacros.h"
28 #include "FWCore/ParameterSet/interface/ParameterSet.h"
29 #include "FWCore/Utilities/interface/InputTag.h"
30 #include "Geometry/CommonDetUnit/interface/GluedGeomDet.h"
31 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
32 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
33 #include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
34 #include "RecoLocalTracker/SiStripClusterizer/interface/SiStripClusterInfo.h"
35 #include "SimDataFormats/PileupSummaryInfo/interface/PileupSummaryInfo.h"
36 #include "TrackingTools/IPTools/interface/IPTools.h"
37 #include "TrackingTools/Records/interface/TransientTrackRecord.h"
38 #include "TrackingTools/TransientTrack/interface/TransientTrack.h"
39 #include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
40 #include "RecoVertex/VertexPrimitives/interface/TransientVertex.h"
41 #include "RecoVertex/KalmanVertexFit/interface/KalmanVertexFitter.h"
42 #include "RecoVertex/VertexTools/interface/VertexDistance3D.h"
43 
44 // ROOT includes
45 #include "TFile.h"
46 #include "TH1.h"
47 #include "TLorentzVector.h"
48 #include "TMath.h"
49 #include "TPRegexp.h"
50 
51 using MVACollection = std::vector<float>;
52 using QualityMaskCollection = std::vector<unsigned char>;
53 
54 class StandaloneTrackMonitor : public DQMEDAnalyzer {
55 public:
56  StandaloneTrackMonitor(const edm::ParameterSet&);
57  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
58 
59 protected:
60  void analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) override;
61  void addClusterToMap(uint32_t detid, const SiStripCluster* cluster);
62  void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
63  void dqmBeginRun(const edm::Run& iRun, const edm::EventSetup& iSetup) override;
64  void processClusters(edm::Event const& iEvent,
65  edm::EventSetup const& iSetup,
66  const TrackerGeometry& tkGeom,
67  double wfac = 1);
68  void processHit(const TrackingRecHit& recHit,
69  edm::EventSetup const& iSetup,
70  const TrackerGeometry& tkGeom,
71  double wfac = 1);
72 
73 private:
74  const reco::Vertex* findClosestVertex(const TransientVertex aTransVtx, const reco::VertexCollection* vertices)
75  const; //same as in /DQMOffline/Alignment/src/DiMuonVertexSelector.cc
76  const std::string moduleName_;
77  const std::string folderName_;
78 
79  const bool isRECO_;
80 
81  SiStripClusterInfo siStripClusterInfo_;
82 
83  const edm::InputTag trackTag_;
84  const edm::InputTag bsTag_;
85  const edm::InputTag vertexTag_;
86  const edm::InputTag puSummaryTag_;
87  const edm::InputTag clusterTag_;
88  const edm::InputTag jetsTag_;
89  const edm::EDGetTokenT<reco::TrackCollection> trackToken_;
90  const edm::EDGetTokenT<reco::BeamSpot> bsToken_;
91  const edm::EDGetTokenT<reco::VertexCollection> vertexToken_;
92  const edm::EDGetTokenT<std::vector<PileupSummaryInfo> > puSummaryToken_;
93  const edm::EDGetTokenT<edmNew::DetSetVector<SiStripCluster> > clusterToken_;
94  const edm::EDGetTokenT<std::vector<reco::PFJet> > jetsToken_;
95  // track MVA
96  const std::string trackQuality_;
97  const bool doPUCorrection_;
98  const bool doTrackCorrection_;
99  const bool isMC_;
100  const bool haveAllHistograms_;
101  const std::string puScaleFactorFile_;
102  const std::string trackScaleFactorFile_;
103  const std::vector<std::string> mvaProducers_;
104  const edm::InputTag mvaTrackTag_;
105  const edm::EDGetTokenT<edm::View<reco::Track> > mvaTrackToken_;
106  const edm::InputTag tcProducer_;
107  const std::string algoName_;
108 
109  int nevt = 0;
110  int chi2it = 0, chi2itGt = 0, chi2itLt = 0;
111  const bool verbose_;
112  std::vector<std::tuple<edm::EDGetTokenT<MVACollection>, edm::EDGetTokenT<QualityMaskCollection> > > mvaQualityTokens_;
113  const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> geomToken_;
114  const edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> transTrackToken_;
115  const edm::ParameterSet TrackEtaHistoPar_;
116  const edm::ParameterSet TrackPhiHistoPar_;
117  const edm::ParameterSet TrackPtHistoPar_;
118  const TrackerGeometry* tkGeom_ = nullptr;
119 
120  MonitorElement* trackEtaH_;
121  MonitorElement* trackEtaerrH_;
122  MonitorElement* trackCosThetaH_;
123  MonitorElement* trackThetaerrH_;
124  MonitorElement* trackPhiH_;
125  MonitorElement* trackPhierrH_;
126  MonitorElement* trackPH_;
127  MonitorElement* trackPtH_;
128  MonitorElement* trackPt_ZoomH_;
129  MonitorElement* trackPtUpto2GeVH_;
130  MonitorElement* trackPtOver10GeVH_;
131  MonitorElement* trackPterrH_;
132  MonitorElement* trackqOverpH_;
133  MonitorElement* trackqOverperrH_;
134  MonitorElement* trackChargeH_;
135  MonitorElement* trackChi2H_;
136  MonitorElement* tracknDOFH_;
137  MonitorElement* trackChi2ProbH_;
138  MonitorElement* trackChi2oNDFH_;
139  MonitorElement* trackd0H_;
140  MonitorElement* trackChi2bynDOFH_;
141  MonitorElement* trackalgoH_;
142  MonitorElement* trackorigalgoH_;
143  MonitorElement* trackStoppingSourceH_;
144 
145  MonitorElement* DistanceOfClosestApproachToPVH_;
146  MonitorElement* DistanceOfClosestApproachToPVZoomedH_;
147  MonitorElement* DistanceOfClosestApproachToPVVsPhiH_;
148  MonitorElement* xPointOfClosestApproachVsZ0wrtPVH_;
149  MonitorElement* yPointOfClosestApproachVsZ0wrtPVH_;
150  MonitorElement* trackDeltaRwrtClosestTrack_;
151 
152  MonitorElement* ip2dToPVH_;
153  MonitorElement* iperr2dToPVH_;
154 
155  MonitorElement* ip3dToBSH_;
156  MonitorElement* iperr3dToBSH_;
157  MonitorElement* iperr3dToBSWtH_;
158  MonitorElement* iperr2dToBSH_;
159  MonitorElement* iperr2dToPVWtH_;
160 
161  MonitorElement* ip2dToBSH_;
162  MonitorElement* sip2dToBSH_;
163 
164  MonitorElement* ip3dToPVH_;
165  MonitorElement* iperr3dToPVH_;
166  MonitorElement* iperr3dToPVWtH_;
167  MonitorElement* sip3dToPVH_;
168  MonitorElement* sip3dToBSH_;
169  MonitorElement* sip2dToPVH_;
170  MonitorElement* sip2dToPVWtH_;
171  MonitorElement* sipDxyToPVH_;
172  MonitorElement* sipDzToPVH_;
173 
174  MonitorElement* ip3dToPV2validpixelhitsH_;
175  MonitorElement* ip3dToBS2validpixelhitsH_;
176  MonitorElement* iperr3dToPV2validpixelhitsH_;
177  MonitorElement* iperr3dToBS2validpixelhitsH_;
178  MonitorElement* sip3dToPV2validpixelhitsH_;
179  MonitorElement* sip3dToBS2validpixelhitsH_;
180 
181  MonitorElement* nallHitsH_;
182  MonitorElement* ntrackerHitsH_;
183 
184  MonitorElement* nvalidTrackerHitsH_;
185  MonitorElement* nvalidPixelHitsH_;
186  MonitorElement* nvalidPixelBHitsH_;
187  MonitorElement* nvalidPixelEHitsH_;
188  MonitorElement* nvalidStripHitsH_;
189  MonitorElement* nvalidTIBHitsH_;
190  MonitorElement* nvalidTOBHitsH_;
191  MonitorElement* nvalidTIDHitsH_;
192  MonitorElement* nvalidTECHitsH_;
193 
194  MonitorElement* nlostTrackerHitsH_;
195  MonitorElement* nlostPixelHitsH_;
196  MonitorElement* nlostPixelBHitsH_;
197  MonitorElement* nlostPixelEHitsH_;
198  MonitorElement* nlostStripHitsH_;
199  MonitorElement* nlostTIBHitsH_;
200  MonitorElement* nlostTOBHitsH_;
201  MonitorElement* nlostTIDHitsH_;
202  MonitorElement* nlostTECHitsH_;
203 
204  MonitorElement* nMissingInnerHitBH_;
205  MonitorElement* nMissingInnerHitEH_;
206  MonitorElement* nMissingOuterHitBH_;
207  MonitorElement* nMissingOuterHitEH_;
208 
209  MonitorElement* trkLayerwithMeasurementH_;
210  MonitorElement* pixelLayerwithMeasurementH_;
211  MonitorElement* pixelBLayerwithMeasurementH_;
212  MonitorElement* pixelELayerwithMeasurementH_;
213  MonitorElement* stripLayerwithMeasurementH_;
214  MonitorElement* stripTIBLayerwithMeasurementH_;
215  MonitorElement* stripTOBLayerwithMeasurementH_;
216  MonitorElement* stripTIDLayerwithMeasurementH_;
217  MonitorElement* stripTECLayerwithMeasurementH_;
218 
219  MonitorElement* nlostHitsH_;
220  MonitorElement* nMissingExpectedInnerHitsH_;
221  MonitorElement* nMissingExpectedOuterHitsH_;
222 
223  MonitorElement* beamSpotXYposH_;
224  MonitorElement* beamSpotXYposerrH_;
225  MonitorElement* beamSpotZposH_;
226  MonitorElement* beamSpotZposerrH_;
227 
228  MonitorElement* vertexXposH_;
229  MonitorElement* vertexYposH_;
230  MonitorElement* vertexZposH_;
231  MonitorElement* nVertexH_;
232  MonitorElement* nVtxH_;
233 
234  MonitorElement* nPixBarrelH_;
235  MonitorElement* nPixEndcapH_;
236  MonitorElement* nStripTIBH_;
237  MonitorElement* nStripTOBH_;
238  MonitorElement* nStripTECH_;
239  MonitorElement* nStripTIDH_;
240  MonitorElement* nTracksH_;
241 
242  // MC only
243  MonitorElement* bunchCrossingH_;
244  MonitorElement* nPUH_;
245  MonitorElement* trueNIntH_;
246 
247  // Exclusive Quantities
248  MonitorElement* nLostHitByLayerH_;
249  MonitorElement* nLostHitByLayerPixH_;
250  MonitorElement* nLostHitByLayerStripH_;
251  MonitorElement* nLostHitsVspTH_;
252  MonitorElement* nLostHitsVsEtaH_;
253  MonitorElement* nLostHitsVsCosThetaH_;
254  MonitorElement* nLostHitsVsPhiH_;
255  MonitorElement* nLostHitsVsIterationH_;
256 
257  MonitorElement* nHitsTIBSVsEtaH_;
258  MonitorElement* nHitsTOBSVsEtaH_;
259  MonitorElement* nHitsTECSVsEtaH_;
260  MonitorElement* nHitsTIDSVsEtaH_;
261  MonitorElement* nHitsStripSVsEtaH_;
262 
263  MonitorElement* nHitsTIBDVsEtaH_;
264  MonitorElement* nHitsTOBDVsEtaH_;
265  MonitorElement* nHitsTECDVsEtaH_;
266  MonitorElement* nHitsTIDDVsEtaH_;
267  MonitorElement* nHitsStripDVsEtaH_;
268 
269  MonitorElement* nValidHitsVspTH_;
270  MonitorElement* nValidHitsVsnVtxH_;
271  MonitorElement* nValidHitsVsEtaH_;
272  MonitorElement* nValidHitsVsCosThetaH_;
273  MonitorElement* nValidHitsVsPhiH_;
274 
275  MonitorElement* nValidHitsPixVsEtaH_;
276  MonitorElement* nValidHitsPixBVsEtaH_;
277  MonitorElement* nValidHitsPixEVsEtaH_;
278  MonitorElement* nValidHitsStripVsEtaH_;
279  MonitorElement* nValidHitsTIBVsEtaH_;
280  MonitorElement* nValidHitsTOBVsEtaH_;
281  MonitorElement* nValidHitsTECVsEtaH_;
282  MonitorElement* nValidHitsTIDVsEtaH_;
283 
284  MonitorElement* nValidHitsPixVsPhiH_;
285  MonitorElement* nValidHitsPixBVsPhiH_;
286  MonitorElement* nValidHitsPixEVsPhiH_;
287  MonitorElement* nValidHitsStripVsPhiH_;
288  MonitorElement* nValidHitsTIBVsPhiH_;
289  MonitorElement* nValidHitsTOBVsPhiH_;
290  MonitorElement* nValidHitsTECVsPhiH_;
291  MonitorElement* nValidHitsTIDVsPhiH_;
292 
293  MonitorElement* nLostHitsPixVsEtaH_;
294  MonitorElement* nLostHitsPixBVsEtaH_;
295  MonitorElement* nLostHitsPixEVsEtaH_;
296  MonitorElement* nLostHitsStripVsEtaH_;
297  MonitorElement* nLostHitsTIBVsEtaH_;
298  MonitorElement* nLostHitsTOBVsEtaH_;
299  MonitorElement* nLostHitsTECVsEtaH_;
300  MonitorElement* nLostHitsTIDVsEtaH_;
301 
302  MonitorElement* nLostHitsPixVsIterationH_;
303  MonitorElement* nLostHitsPixBVsIterationH_;
304  MonitorElement* nLostHitsPixEVsIterationH_;
305  MonitorElement* nLostHitsStripVsIterationH_;
306  MonitorElement* nLostHitsTIBVsIterationH_;
307  MonitorElement* nLostHitsTOBVsIterationH_;
308  MonitorElement* nLostHitsTECVsIterationH_;
309  MonitorElement* nLostHitsTIDVsIterationH_;
310 
311  MonitorElement* nLostHitsPixVsPhiH_;
312  MonitorElement* nLostHitsPixBVsPhiH_;
313  MonitorElement* nLostHitsPixEVsPhiH_;
314  MonitorElement* nLostHitsStripVsPhiH_;
315  MonitorElement* nLostHitsTIBVsPhiH_;
316  MonitorElement* nLostHitsTOBVsPhiH_;
317  MonitorElement* nLostHitsTECVsPhiH_;
318  MonitorElement* nLostHitsTIDVsPhiH_;
319 
320  MonitorElement* trackChi2oNDFVsEtaH_;
321  MonitorElement* trackChi2oNDFVsPhiH_;
322  MonitorElement* trackChi2probVsEtaH_;
323  MonitorElement* trackChi2probVsPhiH_;
324 
325  MonitorElement* trackIperr3dVsEtaH_;
326  MonitorElement* trackIperr3dVsChi2probH_;
327 
328  MonitorElement* trackSip2dVsEtaH_;
329 
330  MonitorElement* trackIperr3dVsEta2DH_;
331  MonitorElement* trackIperr3dVsChi2prob2DH_;
332 
333  MonitorElement* trackSip2dVsEta2DH_;
334  MonitorElement* trackSip2dVsChi2prob2DH_;
335 
336  MonitorElement* hOnTrkClusChargeThinH_;
337  MonitorElement* hOnTrkClusWidthThinH_;
338  MonitorElement* hOnTrkClusChargeThickH_;
339  MonitorElement* hOnTrkClusWidthThickH_;
340 
341  MonitorElement* hOffTrkClusChargeThinH_;
342  MonitorElement* hOffTrkClusWidthThinH_;
343  MonitorElement* hOffTrkClusChargeThickH_;
344  MonitorElement* hOffTrkClusWidthThickH_;
345 
346  std::vector<MonitorElement*> trackMVAs;
347  std::vector<MonitorElement*> trackMVAsHP;
348  std::vector<MonitorElement*> trackMVAsVsPtProfile;
349  std::vector<MonitorElement*> trackMVAsHPVsPtProfile;
350  std::vector<MonitorElement*> trackMVAsVsEtaProfile;
351  std::vector<MonitorElement*> trackMVAsHPVsEtaProfile;
352 
353  MonitorElement* nJet_;
354  MonitorElement* Jet_pt_;
355  MonitorElement* Jet_eta_;
356  MonitorElement* Jet_energy_;
357  MonitorElement* Jet_chargedMultiplicity_;
358 
359  MonitorElement* Zpt_;
360  MonitorElement* ZInvMass_;
361 
362  MonitorElement* cosPhi3DdileptonH_;
363 
364  std::vector<int> lumivec1;
365  std::vector<int> lumivec2;
366  std::vector<float> vpu_;
367  std::vector<float> vtrack_;
368  std::map<uint32_t, std::set<const SiStripCluster*> > clusterMap_;
369 };
370 
371 void StandaloneTrackMonitor::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
372  edm::ParameterSetDescription desc;
373  desc.addUntracked<std::string>("moduleName", "StandaloneTrackMonitor");
374  desc.addUntracked<std::string>("folderName", "highPurityTracks");
375  desc.addUntracked<bool>("isRECO", false);
376  desc.addUntracked<edm::InputTag>("trackInputTag", edm::InputTag("generalTracks"));
377  desc.addUntracked<edm::InputTag>("offlineBeamSpot", edm::InputTag("offlineBeamSpot"));
378  desc.addUntracked<edm::InputTag>("vertexTag", edm::InputTag("offlinePrimaryVertices"));
379  desc.addUntracked<edm::InputTag>("puTag", edm::InputTag("addPileupInfo"));
380  desc.addUntracked<edm::InputTag>("clusterTag", edm::InputTag("siStripClusters"));
381  desc.addUntracked<edm::InputTag>("PFJetsCollection", edm::InputTag("ak4PFJetsCHS"));
382 
383  desc.addUntracked<std::string>("trackQuality", "highPurity");
384  desc.addUntracked<bool>("doPUCorrection", false);
385  desc.addUntracked<bool>("doTrackCorrection", false);
386  desc.addUntracked<bool>("isMC", false);
387  desc.addUntracked<bool>("haveAllHistograms", false);
388  desc.addUntracked<std::string>("puScaleFactorFile", "PileupScaleFactor.root");
389  desc.addUntracked<std::string>("trackScaleFactorFile", "PileupScaleFactor.root");
390  desc.addUntracked<std::vector<std::string> >("MVAProducers", {"initialStepClassifier1", "initialStepClassifier2"});
391  desc.addUntracked<edm::InputTag>("TrackProducerForMVA", edm::InputTag("initialStepTracks"));
392 
393  desc.addUntracked<edm::InputTag>("TCProducer", edm::InputTag("initialStepTrackCandidates"));
394  desc.addUntracked<std::string>("AlgoName", "GenTk");
395  desc.addUntracked<bool>("verbose", false);
396 
397  {
398  edm::ParameterSetDescription TrackEtaHistoPar;
399  TrackEtaHistoPar.add<int>("Xbins", 60);
400  TrackEtaHistoPar.add<double>("Xmin", -3.0);
401  TrackEtaHistoPar.add<double>("Xmax", 3.0);
402  desc.add<edm::ParameterSetDescription>("trackEtaH", TrackEtaHistoPar);
403  }
404 
405  {
406  edm::ParameterSetDescription TrackPtHistoPar;
407  TrackPtHistoPar.add<int>("Xbins", 60);
408  TrackPtHistoPar.add<double>("Xmin", 0.0);
409  TrackPtHistoPar.add<double>("Xmax", 100.0);
410  desc.add<edm::ParameterSetDescription>("trackPtH", TrackPtHistoPar);
411  }
412 
413  {
414  edm::ParameterSetDescription TrackPhiHistoPar;
415  TrackPhiHistoPar.add<int>("Xbins", 100);
416  TrackPhiHistoPar.add<double>("Xmin", -M_PI);
417  TrackPhiHistoPar.add<double>("Xmax", M_PI);
418  desc.add<edm::ParameterSetDescription>("trackPhiH", TrackPhiHistoPar);
419  }
420 
421  descriptions.add("standaloneTrackMonitorDefault", desc);
422 }
423 
424 // -----------------------------
425 // constructors and destructor
426 // -----------------------------
427 StandaloneTrackMonitor::StandaloneTrackMonitor(const edm::ParameterSet& ps)
428  : moduleName_(ps.getUntrackedParameter<std::string>("moduleName", "StandaloneTrackMonitor")),
429  folderName_(ps.getUntrackedParameter<std::string>("folderName", "highPurityTracks")),
430  isRECO_(ps.getUntrackedParameter<bool>("isRECO", false)),
431  siStripClusterInfo_(consumesCollector()),
432  trackTag_(ps.getUntrackedParameter<edm::InputTag>("trackInputTag", edm::InputTag("generalTracks"))),
433  bsTag_(ps.getUntrackedParameter<edm::InputTag>("offlineBeamSpot", edm::InputTag("offlineBeamSpot"))),
434  vertexTag_(ps.getUntrackedParameter<edm::InputTag>("vertexTag", edm::InputTag("offlinePrimaryVertices"))),
435  puSummaryTag_(ps.getUntrackedParameter<edm::InputTag>("puTag", edm::InputTag("addPileupInfo"))),
436  clusterTag_(ps.getUntrackedParameter<edm::InputTag>("clusterTag", edm::InputTag("siStripClusters"))),
437  jetsTag_(ps.getUntrackedParameter<edm::InputTag>("PFJetsCollection", edm::InputTag("ak4PFJetsCHS"))),
438  trackToken_(consumes<reco::TrackCollection>(trackTag_)),
439  bsToken_(consumes<reco::BeamSpot>(bsTag_)),
440  vertexToken_(consumes<reco::VertexCollection>(vertexTag_)),
441  puSummaryToken_(consumes<std::vector<PileupSummaryInfo> >(puSummaryTag_)),
442  clusterToken_(consumes<edmNew::DetSetVector<SiStripCluster> >(clusterTag_)),
443  jetsToken_(consumes<std::vector<reco::PFJet> >(jetsTag_)),
444  trackQuality_(ps.getUntrackedParameter<std::string>("trackQuality", "highPurity")),
445  doPUCorrection_(ps.getUntrackedParameter<bool>("doPUCorrection", false)),
446  doTrackCorrection_(ps.getUntrackedParameter<bool>("doTrackCorrection", false)),
447  isMC_(ps.getUntrackedParameter<bool>("isMC", false)),
448  haveAllHistograms_(ps.getUntrackedParameter<bool>("haveAllHistograms", false)),
449  puScaleFactorFile_(ps.getUntrackedParameter<std::string>("puScaleFactorFile", "PileupScaleFactor.root")),
450  trackScaleFactorFile_(ps.getUntrackedParameter<std::string>("trackScaleFactorFile", "PileupScaleFactor.root")),
451  mvaProducers_(ps.getUntrackedParameter<std::vector<std::string> >("MVAProducers")),
452  mvaTrackTag_(ps.getUntrackedParameter<edm::InputTag>("TrackProducerForMVA")),
453  mvaTrackToken_(consumes<edm::View<reco::Track> >(mvaTrackTag_)),
454  tcProducer_(ps.getUntrackedParameter<edm::InputTag>("TCProducer")),
455  algoName_(ps.getUntrackedParameter<std::string>("AlgoName")),
456  verbose_(ps.getUntrackedParameter<bool>("verbose", false)),
457  geomToken_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord, edm::Transition::BeginRun>()),
458  transTrackToken_(esConsumes<TransientTrackBuilder, TransientTrackRecord, edm::Transition::Event>(
459  edm::ESInputTag{"", "TransientTrackBuilder"})),
460  TrackEtaHistoPar_(ps.getParameter<edm::ParameterSet>("trackEtaH")),
461  TrackPhiHistoPar_(ps.getParameter<edm::ParameterSet>("trackPhiH")),
462  TrackPtHistoPar_(ps.getParameter<edm::ParameterSet>("trackPtH")) {
463  for (const auto& v : mvaProducers_) {
464  mvaQualityTokens_.push_back(std::make_tuple(consumes<MVACollection>(edm::InputTag(v, "MVAValues")),
465  consumes<QualityMaskCollection>(edm::InputTag(v, "QualityMasks"))));
466  }
467 
468  trackEtaH_ = nullptr;
469  trackEtaerrH_ = nullptr;
470  trackCosThetaH_ = nullptr;
471  trackThetaerrH_ = nullptr;
472  trackPhiH_ = nullptr;
473  trackPhierrH_ = nullptr;
474  trackPH_ = nullptr;
475  trackPtH_ = nullptr;
476  trackPtUpto2GeVH_ = nullptr;
477  trackPtOver10GeVH_ = nullptr;
478  trackPterrH_ = nullptr;
479  trackqOverpH_ = nullptr;
480  trackqOverperrH_ = nullptr;
481  trackChargeH_ = nullptr;
482  nlostHitsH_ = nullptr;
483  nvalidTrackerHitsH_ = nullptr;
484  nvalidPixelHitsH_ = nullptr;
485  nvalidStripHitsH_ = nullptr;
486  trkLayerwithMeasurementH_ = nullptr;
487  pixelLayerwithMeasurementH_ = nullptr;
488  stripLayerwithMeasurementH_ = nullptr;
489  beamSpotXYposH_ = nullptr;
490  beamSpotXYposerrH_ = nullptr;
491  beamSpotZposH_ = nullptr;
492  beamSpotZposerrH_ = nullptr;
493  trackChi2H_ = nullptr;
494  tracknDOFH_ = nullptr;
495  trackd0H_ = nullptr;
496  trackChi2bynDOFH_ = nullptr;
497  vertexXposH_ = nullptr;
498  vertexYposH_ = nullptr;
499  vertexZposH_ = nullptr;
500 
501  nPixBarrelH_ = nullptr;
502  nPixEndcapH_ = nullptr;
503  nStripTIBH_ = nullptr;
504  nStripTOBH_ = nullptr;
505  nStripTECH_ = nullptr;
506  nStripTIDH_ = nullptr;
507 
508  // for MC only
509  nVtxH_ = nullptr;
510  nVertexH_ = nullptr;
511  bunchCrossingH_ = nullptr;
512  nPUH_ = nullptr;
513  trueNIntH_ = nullptr;
514 
515  nLostHitsVspTH_ = nullptr;
516  nLostHitsVsEtaH_ = nullptr;
517  nLostHitsVsCosThetaH_ = nullptr;
518  nLostHitsVsPhiH_ = nullptr;
519  nLostHitsVsIterationH_ = nullptr;
520 
521  nHitsTIBSVsEtaH_ = nullptr;
522  nHitsTOBSVsEtaH_ = nullptr;
523  nHitsTECSVsEtaH_ = nullptr;
524  nHitsTIDSVsEtaH_ = nullptr;
525  nHitsStripSVsEtaH_ = nullptr;
526 
527  nHitsTIBDVsEtaH_ = nullptr;
528  nHitsTOBDVsEtaH_ = nullptr;
529  nHitsTECDVsEtaH_ = nullptr;
530  nHitsTIDDVsEtaH_ = nullptr;
531  nHitsStripDVsEtaH_ = nullptr;
532 
533  nValidHitsVspTH_ = nullptr;
534  nValidHitsVsEtaH_ = nullptr;
535  nValidHitsVsCosThetaH_ = nullptr;
536  nValidHitsVsPhiH_ = nullptr;
537  nValidHitsVsnVtxH_ = nullptr;
538 
539  nValidHitsPixVsEtaH_ = nullptr;
540  nValidHitsPixBVsEtaH_ = nullptr;
541  nValidHitsPixEVsEtaH_ = nullptr;
542  nValidHitsStripVsEtaH_ = nullptr;
543  nValidHitsTIBVsEtaH_ = nullptr;
544  nValidHitsTOBVsEtaH_ = nullptr;
545  nValidHitsTECVsEtaH_ = nullptr;
546  nValidHitsTIDVsEtaH_ = nullptr;
547 
548  nValidHitsPixVsPhiH_ = nullptr;
549  nValidHitsPixBVsPhiH_ = nullptr;
550  nValidHitsPixEVsPhiH_ = nullptr;
551  nValidHitsStripVsPhiH_ = nullptr;
552  nValidHitsTIBVsPhiH_ = nullptr;
553  nValidHitsTOBVsPhiH_ = nullptr;
554  nValidHitsTECVsPhiH_ = nullptr;
555  nValidHitsTIDVsPhiH_ = nullptr;
556 
557  nLostHitsPixVsEtaH_ = nullptr;
558  nLostHitsPixBVsEtaH_ = nullptr;
559  nLostHitsPixEVsEtaH_ = nullptr;
560  nLostHitsStripVsEtaH_ = nullptr;
561  nLostHitsTIBVsEtaH_ = nullptr;
562  nLostHitsTOBVsEtaH_ = nullptr;
563  nLostHitsTECVsEtaH_ = nullptr;
564  nLostHitsTIDVsEtaH_ = nullptr;
565 
566  nLostHitsPixVsPhiH_ = nullptr;
567  nLostHitsPixBVsPhiH_ = nullptr;
568  nLostHitsPixEVsPhiH_ = nullptr;
569  nLostHitsStripVsPhiH_ = nullptr;
570  nLostHitsTIBVsPhiH_ = nullptr;
571  nLostHitsTOBVsPhiH_ = nullptr;
572  nLostHitsTECVsPhiH_ = nullptr;
573  nLostHitsTIDVsPhiH_ = nullptr;
574 
575  nLostHitsPixVsIterationH_ = nullptr;
576  nLostHitsPixBVsIterationH_ = nullptr;
577  nLostHitsPixEVsIterationH_ = nullptr;
578  nLostHitsStripVsIterationH_ = nullptr;
579  nLostHitsTIBVsIterationH_ = nullptr;
580  nLostHitsTOBVsIterationH_ = nullptr;
581  nLostHitsTECVsIterationH_ = nullptr;
582  nLostHitsTIDVsIterationH_ = nullptr;
583 
584  hOnTrkClusChargeThinH_ = nullptr;
585  hOnTrkClusWidthThinH_ = nullptr;
586  hOnTrkClusChargeThickH_ = nullptr;
587  hOnTrkClusWidthThickH_ = nullptr;
588 
589  hOffTrkClusChargeThinH_ = nullptr;
590  hOffTrkClusWidthThinH_ = nullptr;
591  hOffTrkClusChargeThickH_ = nullptr;
592  hOffTrkClusWidthThickH_ = nullptr;
593 
594  cosPhi3DdileptonH_ = nullptr;
595  ip3dToPV2validpixelhitsH_ = nullptr;
596  ip3dToBS2validpixelhitsH_ = nullptr;
597  iperr3dToPV2validpixelhitsH_ = nullptr;
598  iperr3dToBS2validpixelhitsH_ = nullptr;
599  sip3dToPV2validpixelhitsH_ = nullptr;
600  sip3dToBS2validpixelhitsH_ = nullptr;
601 
602  // Read pileup weight factors
603 
604  if (isMC_ && doPUCorrection_ && doTrackCorrection_) {
605  throw std::runtime_error("if isMC is true, only one of doPUCorrection and doTrackCorrection can be true");
606  }
607 
608  if (isMC_ && doPUCorrection_) {
609  vpu_.clear();
610  TFile* f1 = TFile::Open(puScaleFactorFile_.c_str());
611  TH1F* h1 = dynamic_cast<TH1F*>(f1->Get("pileupweight"));
612  for (int i = 1; i <= h1->GetNbinsX(); ++i)
613  vpu_.push_back(h1->GetBinContent(i));
614  f1->Close();
615  }
616 
617  if (isMC_ && doTrackCorrection_) {
618  vtrack_.clear();
619  TFile* f1 = TFile::Open(trackScaleFactorFile_.c_str());
620  TH1F* h1 = dynamic_cast<TH1F*>(f1->Get("trackweight"));
621  for (int i = 1; i <= h1->GetNbinsX(); ++i)
622  vtrack_.push_back(h1->GetBinContent(i));
623  f1->Close();
624  }
625 }
626 
627 void StandaloneTrackMonitor::dqmBeginRun(const edm::Run& iRun, const edm::EventSetup& iSetup) {
628  tkGeom_ = &(iSetup.getData(geomToken_));
629 }
630 
631 void StandaloneTrackMonitor::bookHistograms(DQMStore::IBooker& ibook,
632  edm::Run const& iRun,
633  edm::EventSetup const& iSetup) {
634  std::string currentFolder = moduleName_ + "/" + folderName_;
635  ibook.setCurrentFolder(currentFolder);
636 
637  // The following are common with the official tool
638  if (haveAllHistograms_) {
639  trackEtaH_ = ibook.book1DD("trackEta",
640  "Track Eta",
641  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
642  TrackEtaHistoPar_.getParameter<double>("Xmin"),
643  TrackEtaHistoPar_.getParameter<double>("Xmax"));
644 
645  trackEtaerrH_ = ibook.book1DD("trackEtaerr", "Track Eta Error", 50, 0.0, 1.0);
646  trackCosThetaH_ = ibook.book1DD("trackCosTheta", "Track Cos(Theta)", 50, -1.0, 1.0);
647  trackThetaerrH_ = ibook.book1DD("trackThetaerr", "Track Theta Error", 50, 0.0, 1.0);
648  trackPhiH_ = ibook.book1DD("trackPhi", "Track Phi", 70, -3.5, 3.5);
649  trackPhierrH_ = ibook.book1DD("trackPhierr", "Track Phi Error", 50, 0.0, 1.0);
650 
651  trackPH_ = ibook.book1DD("trackP", "Track 4-momentum", 50, 0.0, 10.0);
652  trackPtH_ = ibook.book1DD("trackPt",
653  "Track Pt",
654  TrackPtHistoPar_.getParameter<int32_t>("Xbins"),
655  TrackPtHistoPar_.getParameter<double>("Xmin"),
656  TrackPtHistoPar_.getParameter<double>("Xmax"));
657  trackPt_ZoomH_ = ibook.book1DD("trackPt_Zoom", "Track Pt", 100, 60, 70);
658  trackPterrH_ = ibook.book1DD("trackPterr", "Track Pt Error", 100, 0.0, 100.0);
659  trackqOverpH_ = ibook.book1DD("trackqOverp", "q Over p", 100, -0.5, 0.5);
660  trackqOverperrH_ = ibook.book1DD("trackqOverperr", "q Over p Error", 50, 0.0, 0.5);
661  trackChargeH_ = ibook.book1DD("trackCharge", "Track Charge", 3, -1.5, 1.5);
662  trackChi2H_ = ibook.book1DD("trackChi2", "Chi2", 100, 0.0, 100.0);
663  tracknDOFH_ = ibook.book1DD("tracknDOF", "nDOF", 100, 0.0, 100.0);
664  trackChi2ProbH_ = ibook.book1DD("trackChi2Prob", "Chi2prob", 50, 0.0, 1.0);
665  trackChi2oNDFH_ = ibook.book1DD("trackChi2oNDF", "Chi2oNDF", 100, 0.0, 100.0);
666  trackd0H_ = ibook.book1DD("trackd0", "Track d0", 100, -1, 1);
667  trackChi2bynDOFH_ = ibook.book1DD("trackChi2bynDOF", "Chi2 Over nDOF", 100, 0.0, 10.0);
668  trackalgoH_ = ibook.book1DD("trackalgo", "Track Algo", 46, 0.0, 46.0);
669  trackorigalgoH_ = ibook.book1DD("trackorigalgo", "Track Original Algo", 46, 0.0, 46.0);
670 
671  for (size_t ibin = 0; ibin < reco::TrackBase::algoSize - 1; ibin++) {
672  trackalgoH_->setBinLabel(ibin + 1, reco::TrackBase::algoNames[ibin], 1);
673  trackorigalgoH_->setBinLabel(ibin + 1, reco::TrackBase::algoNames[ibin], 1);
674  }
675 
676  trackStoppingSourceH_ = ibook.book1DD("trackstoppingsource", "Track Stopping Source", 12, 0.0, 12.0);
677 
678  // DataFormats/TrackCandidate/interface/TrajectoryStopReasons.h
679  size_t StopReasonNameSize = static_cast<size_t>(StopReason::SIZE);
680  for (unsigned int i = 0; i < StopReasonNameSize; ++i) {
681  trackStoppingSourceH_->setBinLabel(i + 1, StopReasonName::StopReasonName[i], 1);
682  }
683 
684  DistanceOfClosestApproachToPVH_ =
685  ibook.book1DD("DistanceOfClosestApproachToPV", "DistanceOfClosestApproachToPV", 1000, -1.0, 1.0);
686  DistanceOfClosestApproachToPVZoomedH_ =
687  ibook.book1DD("DistanceOfClosestApproachToPVZoomed", "DistanceOfClosestApproachToPV", 1000, -0.1, 0.1);
688  DistanceOfClosestApproachToPVVsPhiH_ = ibook.bookProfile(
689  "DistanceOfClosestApproachToPVVsPhi", "DistanceOfClosestApproachToPVVsPhi", 100, -3.5, 3.5, 0.0, 0.0, "g");
690  xPointOfClosestApproachVsZ0wrtPVH_ = ibook.bookProfile(
691  "xPointOfClosestApproachVsZ0wrtPV", "xPointOfClosestApproachVsZ0wrtPV", 120, -60, 60, 0.0, 0.0, "g");
692  yPointOfClosestApproachVsZ0wrtPVH_ = ibook.bookProfile(
693  "yPointOfClosestApproachVsZ0wrtPV", "yPointOfClosestApproachVsZ0wrtPV", 120, -60, 60, 0.0, 0.0, "g");
694  trackDeltaRwrtClosestTrack_ =
695  ibook.book1DD("trackDeltaRwrtClosestTrack", "min#DeltaR(considered track,other tracks)", 500, 0, 10);
696 
697  ip2dToPVH_ = ibook.book1DD("ip2dToPV", "IP in 2d To PV", 1000, -1.0, 1.0);
698  unsigned int niperrbins = 100;
699  float iperrbinning[niperrbins + 1];
700  float miniperr = 0.0001, maxiperr = 5;
701  iperrbinning[0] = miniperr;
702  for (unsigned int i = 1; i != niperrbins + 1; i++) {
703  iperrbinning[i] = iperrbinning[i - 1] * pow(maxiperr / miniperr, 1. / niperrbins);
704  }
705  iperr2dToPVH_ = ibook.book1DD("iperr2dToPV", "IP error in 2d To PV", niperrbins, iperrbinning);
706  iperr2dToPVWtH_ = ibook.book1DD("iperr2dToPVWt", "IP error in 2d To PV", niperrbins, iperrbinning);
707 
708  ip3dToPVH_ = ibook.book1DD("ip3dToPV", "IP in 3d To PV", 200, -20, 20);
709  ip3dToBSH_ = ibook.book1DD("ip3dToBS", "IP in 3d To BS", 200, -20, 20);
710  iperr3dToPVH_ = ibook.book1DD("iperr3dToPV", "IP error in 3d To PV", niperrbins, iperrbinning);
711  iperr3dToBSH_ = ibook.book1DD("iperr3dToBS", "IP error in 3d To BS", niperrbins, iperrbinning);
712  sip3dToPVH_ = ibook.book1DD("sip3dToPV", "IP significance in 3d To PV", 200, -10, 10);
713  sip3dToBSH_ = ibook.book1DD("sip3dToBS", "IP significance in 3d To BS", 200, -10, 10);
714 
715  ip3dToPV2validpixelhitsH_ =
716  ibook.book1DD("ip3dToPV2validpixelhits", "IP in 3d To PV (nValidPixelHits>2)", 200, -0.20, 0.20);
717  ip3dToBS2validpixelhitsH_ =
718  ibook.book1DD("ip3dToBS2validpixelhits", "IP in 3d To BS (nValidPixelHits>2)", 200, -0.20, 0.20);
719  iperr3dToPV2validpixelhitsH_ = ibook.book1DD(
720  "iperr3dToPV2validpixelhits", "IP error in 3d To PV (nValidPixelHits>2)", niperrbins, iperrbinning);
721  iperr3dToBS2validpixelhitsH_ = ibook.book1DD(
722  "iperr3dToBS2validpixelhits", "IP error in 3d To BS (nValidPixelHits>2)", niperrbins, iperrbinning);
723  sip3dToPV2validpixelhitsH_ =
724  ibook.book1DD("sip3dToPV2validpixelhits", "IP significance in 3d To PV (nValidPixelHits>2)", 200, -10, 10);
725  sip3dToBS2validpixelhitsH_ =
726  ibook.book1DD("sip3dToBS2validpixelhits", "IP significance in 3d To BS (nValidPixelHits>2)", 200, -10, 10);
727 
728  ip2dToBSH_ = ibook.book1DD("ip2dToBS", "IP in 2d To BS", 1000, -1., 1.); //Beamspot
729  iperr2dToBSH_ = ibook.book1DD("iperr2dToBS", "IP error in 2d To BS", niperrbins, iperrbinning);
730  sip2dToBSH_ = ibook.book1DD("sip2dToBS", "IP significance in 2d To BS", 200, -10, 10);
731 
732  iperr3dToPVWtH_ = ibook.book1DD("iperr3dToPVWt", "IP error in 3d To PV", niperrbins, iperrbinning);
733  sip2dToPVH_ = ibook.book1DD("sip2dToPV", "IP significance in 2d To PV", 200, -10, 10);
734 
735  sip2dToPVWtH_ = ibook.book1DD("sip2dToPVWt", "IP significance in 2d To PV", 200, -10, 10);
736  sipDxyToPVH_ = ibook.book1DD("sipDxyToPV", "IP significance in dxy To PV", 100, -10, 10);
737  sipDzToPVH_ = ibook.book1DD("sipDzToPV", "IP significance in dz To PV", 100, -10, 10);
738 
739  nallHitsH_ = ibook.book1DD("nallHits", "No. of All Hits", 60, -0.5, 59.5);
740  ntrackerHitsH_ = ibook.book1DD("ntrackerHits", "No. of Tracker Hits", 60, -0.5, 59.5);
741 
742  nvalidTrackerHitsH_ = ibook.book1DD("nvalidTrackerhits", "No. of Valid Tracker Hits", 47, -0.5, 46.5);
743  nvalidPixelHitsH_ = ibook.book1DD("nvalidPixelHits", "No. of Valid Hits in Pixel", 8, -0.5, 7.5);
744  nvalidPixelBHitsH_ = ibook.book1DD("nvalidPixelBarrelHits", "No. of Valid Hits in Pixel Barrel", 6, -0.5, 5.5);
745  nvalidPixelEHitsH_ = ibook.book1DD("nvalidPixelEndcapHits", "No. of Valid Hits in Pixel Endcap", 6, -0.5, 6.5);
746  nvalidStripHitsH_ = ibook.book1DD("nvalidStripHits", "No. of Valid Hits in Strip", 36, -0.5, 35.5);
747  nvalidTIBHitsH_ = ibook.book1DD("nvalidTIBHits", "No. of Valid Hits in Strip TIB", 6, -0.5, 5.5);
748  nvalidTOBHitsH_ = ibook.book1DD("nvalidTOBHits", "No. of Valid Hits in Strip TOB", 11, -0.5, 10.5);
749  nvalidTIDHitsH_ = ibook.book1DD("nvalidTIDHits", "No. of Valid Hits in Strip TID", 6, -0.5, 5.5);
750  nvalidTECHitsH_ = ibook.book1DD("nvalidTECHits", "No. of Valid Hits in Strip TEC", 11, -0.5, 10.5);
751 
752  nlostTrackerHitsH_ = ibook.book1DD("nlostTrackerhits", "No. of Lost Tracker Hits", 15, -0.5, 14.5);
753  nlostPixelHitsH_ = ibook.book1DD("nlostPixelHits", "No. of Lost Hits in Pixel", 8, -0.5, 7.5);
754  nlostPixelBHitsH_ = ibook.book1DD("nlostPixelBarrelHits", "No. of Lost Hits in Pixel Barrel", 5, -0.5, 4.5);
755  nlostPixelEHitsH_ = ibook.book1DD("nlostPixelEndcapHits", "No. of Lost Hits in Pixel Endcap", 4, -0.5, 3.5);
756  nlostStripHitsH_ = ibook.book1DD("nlostStripHits", "No. of Lost Hits in Strip", 10, -0.5, 9.5);
757  nlostTIBHitsH_ = ibook.book1DD("nlostTIBHits", "No. of Lost Hits in Strip TIB", 5, -0.5, 4.5);
758  nlostTOBHitsH_ = ibook.book1DD("nlostTOBHits", "No. of Lost Hits in Strip TOB", 10, -0.5, 9.5);
759  nlostTIDHitsH_ = ibook.book1DD("nlostTIDHits", "No. of Lost Hits in Strip TID", 5, -0.5, 4.5);
760  nlostTECHitsH_ = ibook.book1DD("nlostTECHits", "No. of Lost Hits in Strip TEC", 10, -0.5, 9.5);
761 
762  trkLayerwithMeasurementH_ = ibook.book1DD("trkLayerwithMeasurement", "No. of Layers per Track", 20, 0.0, 20.0);
763  pixelLayerwithMeasurementH_ =
764  ibook.book1DD("pixelLayerwithMeasurement", "No. of Pixel Layers per Track", 10, 0.0, 10.0);
765  pixelBLayerwithMeasurementH_ =
766  ibook.book1DD("pixelBLayerwithMeasurement", "No. of Pixel Barrel Layers per Track", 5, 0.0, 5.0);
767  pixelELayerwithMeasurementH_ =
768  ibook.book1DD("pixelELayerwithMeasurement", "No. of Pixel Endcap Layers per Track", 5, 0.0, 5.0);
769  stripLayerwithMeasurementH_ =
770  ibook.book1DD("stripLayerwithMeasurement", "No. of Strip Layers per Track", 20, 0.0, 20.0);
771  stripTIBLayerwithMeasurementH_ =
772  ibook.book1DD("stripTIBLayerwithMeasurement", "No. of Strip TIB Layers per Track", 10, 0.0, 10.0);
773  stripTOBLayerwithMeasurementH_ =
774  ibook.book1DD("stripTOBLayerwithMeasurement", "No. of Strip TOB Layers per Track", 10, 0.0, 10.0);
775  stripTIDLayerwithMeasurementH_ =
776  ibook.book1DD("stripTIDLayerwithMeasurement", "No. of Strip TID Layers per Track", 5, 0.0, 5.0);
777  stripTECLayerwithMeasurementH_ =
778  ibook.book1DD("stripTECLayerwithMeasurement", "No. of Strip TEC Layers per Track", 15, 0.0, 15.0);
779 
780  nlostHitsH_ = ibook.book1DD("nlostHits", "No. of Lost Hits", 10, -0.5, 9.5);
781  nMissingExpectedInnerHitsH_ =
782  ibook.book1DD("nMissingExpectedInnerHits", "No. of Missing Expected Inner Hits", 10, -0.5, 9.5);
783  nMissingExpectedOuterHitsH_ =
784  ibook.book1DD("nMissingExpectedOuterHits", "No. of Missing Expected Outer Hits", 10, -0.5, 9.5);
785 
786  beamSpotXYposH_ = ibook.book1DD("beamSpotXYpos", "d_{xy} w.r.t beam spot", 100, -1.0, 1.0);
787  beamSpotXYposerrH_ = ibook.book1DD("beamSpotXYposerr", "error on d_{xy} w.r.t beam spot", 50, 0.0, 0.25);
788  beamSpotZposH_ = ibook.book1DD("beamSpotZpos", "d_{z} w.r.t beam spot", 100, -20.0, 20.0);
789  beamSpotZposerrH_ = ibook.book1DD("beamSpotZposerr", "error on d_{z} w.r.t. beam spot", 50, 0.0, 0.25);
790 
791  vertexXposH_ = ibook.book1DD("vertexXpos", "Vertex X position", 100, -0.1, 0.1);
792  vertexYposH_ = ibook.book1DD("vertexYpos", "Vertex Y position", 200, -0.1, 0.1);
793  vertexZposH_ = ibook.book1DD("vertexZpos", "Vertex Z position", 100, -20.0, 20.0);
794  nVertexH_ = ibook.book1DD("nVertex", "# of vertices", 120, -0.5, 119.5);
795  nVtxH_ = ibook.book1DD("nVtx", "# of vtxs", 120, -0.5, 119.5);
796 
797  nMissingInnerHitBH_ = ibook.book1DD("nMissingInnerHitB", "No. missing inner hit per Track in Barrel", 6, -0.5, 5.5);
798  nMissingInnerHitEH_ = ibook.book1DD("nMissingInnerHitE", "No. missing inner hit per Track in Endcap", 6, -0.5, 5.5);
799  nMissingOuterHitBH_ =
800  ibook.book1DD("nMissingOuterHitB", "No. missing outer hit per Track in Barrel", 11, -0.5, 10.5);
801  nMissingOuterHitEH_ =
802  ibook.book1DD("nMissingOuterHitE", "No. missing outer hit per Track in Endcap", 11, -0.5, 10.5);
803  nPixBarrelH_ = ibook.book1DD("nHitPixelBarrel", "No. of hits in Pixel Barrel per Track", 20, 0, 20.0);
804  nPixEndcapH_ = ibook.book1DD("nHitPixelEndcap", "No. of hits in Pixel Endcap per Track", 20, 0, 20.0);
805  nStripTIBH_ = ibook.book1DD("nHitStripTIB", "No. of hits in Strip TIB per Track", 30, 0, 30.0);
806  nStripTOBH_ = ibook.book1DD("nHitStripTOB", "No. of hits in Strip TOB per Track", 30, 0, 30.0);
807  nStripTECH_ = ibook.book1DD("nHitStripTEC", "No. of hits in Strip TEC per Track", 30, 0, 30.0);
808  nStripTIDH_ = ibook.book1DD("nHitStripTID", "No. of hits in Strip TID per Tracks", 30, 0, 30.0);
809  nTracksH_ = ibook.book1DD("nTracks", "No. of Tracks", 1200, -0.5, 1199.5);
810  nJet_ = ibook.book1DD("nJet", "Number of Jets", 101, -0.5, 100.5);
811  Jet_pt_ = ibook.book1DD("Jet_pt", "Jet p_{T}", 200, 0., 200.);
812  Jet_eta_ = ibook.book1DD("Jet_eta", "Jet #eta", 100, -5.2, 5.2);
813  Jet_energy_ = ibook.book1DD("Jet_energy", "Jet Energy", 200, 0., 200.);
814  Jet_chargedMultiplicity_ =
815  ibook.book1DD("Jet_chargedMultiplicity", "Jet charged Hadron Multiplicity", 201, -0.5, 200.5);
816  Zpt_ = ibook.book1DD("Zpt", "Z-boson transverse momentum", 100, 0, 100);
817  ZInvMass_ = ibook.book1DD("ZInvMass", "m_{ll}", 120, 75, 105);
818  cosPhi3DdileptonH_ = ibook.book1DD("cosPhi3Ddilepton", "cos#Phi_{3D,ll}", 202, -1.01, 1.01);
819  }
820  if (isMC_) {
821  bunchCrossingH_ = ibook.book1DD("bunchCrossing", "Bunch Crossing", 60, 0, 60.0);
822  nPUH_ = ibook.book1DD("nPU", "No of Pileup", 100, 0, 100.0);
823  trueNIntH_ = ibook.book1DD("trueNInt", "True no of Interactions", 100, 0, 100.0);
824  }
825  // Exclusive histograms
826 
827  nLostHitByLayerH_ = ibook.book1DD("nLostHitByLayer", "No. of Lost Hit per Layer", 29, 0.5, 29.5);
828 
829  nLostHitByLayerPixH_ =
830  ibook.book1DD("nLostHitByLayerPix", "No. of Lost Hit per Layer for Pixel detector", 7, 0.5, 7.5);
831 
832  nLostHitByLayerStripH_ =
833  ibook.book1DD("nLostHitByLayerStrip", "No. of Lost Hit per Layer for SiStrip detector", 22, 0.5, 22.5);
834 
835  nLostHitsVspTH_ = ibook.bookProfile("nLostHitsVspT",
836  "Number of Lost Hits Vs pT",
837  TrackPtHistoPar_.getParameter<int32_t>("Xbins"),
838  TrackPtHistoPar_.getParameter<double>("Xmin"),
839  TrackPtHistoPar_.getParameter<double>("Xmax"),
840  0.0,
841  0.0,
842  "g");
843  nLostHitsVsEtaH_ = ibook.bookProfile("nLostHitsVsEta",
844  "Number of Lost Hits Vs Eta",
845  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
846  TrackEtaHistoPar_.getParameter<double>("Xmin"),
847  TrackEtaHistoPar_.getParameter<double>("Xmax"),
848  0.0,
849  0.0,
850  "g");
851  nLostHitsVsCosThetaH_ =
852  ibook.bookProfile("nLostHitsVsCosTheta", "Number of Lost Hits Vs Cos(Theta)", 50, -1.0, 1.0, 0.0, 0.0, "g");
853  nLostHitsVsPhiH_ = ibook.bookProfile("nLostHitsVsPhi", "Number of Lost Hits Vs Phi", 100, -3.5, 3.5, 0.0, 0.0, "g");
854  nLostHitsVsIterationH_ =
855  ibook.bookProfile("nLostHitsVsIteration", "Number of Lost Hits Vs Iteration", 47, -0.5, 46.5, 0.0, 0.0, "g");
856 
857  nHitsTIBSVsEtaH_ = ibook.bookProfile("nHitsTIBSVsEta",
858  "Number of Hits in TIB Vs Eta (Single-sided)",
859  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
860  TrackEtaHistoPar_.getParameter<double>("Xmin"),
861  TrackEtaHistoPar_.getParameter<double>("Xmax"),
862  0.0,
863  0.0,
864  "g");
865  nHitsTOBSVsEtaH_ = ibook.bookProfile("nHitsTOBSVsEta",
866  "Number of Hits in TOB Vs Eta (Single-sided)",
867  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
868  TrackEtaHistoPar_.getParameter<double>("Xmin"),
869  TrackEtaHistoPar_.getParameter<double>("Xmax"),
870  0.0,
871  0.0,
872  "g");
873  nHitsTECSVsEtaH_ = ibook.bookProfile("nHitsTECSVsEta",
874  "Number of Hits in TEC Vs Eta (Single-sided)",
875  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
876  TrackEtaHistoPar_.getParameter<double>("Xmin"),
877  TrackEtaHistoPar_.getParameter<double>("Xmax"),
878  0.0,
879  0.0,
880  "g");
881  nHitsTIDSVsEtaH_ = ibook.bookProfile("nHitsTIDSVsEta",
882  "Number of Hits in TID Vs Eta (Single-sided)",
883  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
884  TrackEtaHistoPar_.getParameter<double>("Xmin"),
885  TrackEtaHistoPar_.getParameter<double>("Xmax"),
886  0.0,
887  0.0,
888  "g");
889 
890  nHitsStripSVsEtaH_ = ibook.bookProfile("nHitsStripSVsEta",
891  "Number of Strip Hits Vs Eta (Single-sided)",
892  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
893  TrackEtaHistoPar_.getParameter<double>("Xmin"),
894  TrackEtaHistoPar_.getParameter<double>("Xmax"),
895  0.0,
896  0.0,
897  "g");
898 
899  nHitsTIBDVsEtaH_ = ibook.bookProfile("nHitsTIBDVsEta",
900  "Number of Hits in TIB Vs Eta (Double-sided)",
901  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
902  TrackEtaHistoPar_.getParameter<double>("Xmin"),
903  TrackEtaHistoPar_.getParameter<double>("Xmax"),
904  0.0,
905  0.0,
906  "g");
907  nHitsTOBDVsEtaH_ = ibook.bookProfile("nHitsTOBDVsEta",
908  "Number of Hits in TOB Vs Eta (Double-sided)",
909  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
910  TrackEtaHistoPar_.getParameter<double>("Xmin"),
911  TrackEtaHistoPar_.getParameter<double>("Xmax"),
912  0.0,
913  0.0,
914  "g");
915  nHitsTECDVsEtaH_ = ibook.bookProfile("nHitsTECDVsEta",
916  "Number of Hits in TEC Vs Eta (Double-sided)",
917  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
918  TrackEtaHistoPar_.getParameter<double>("Xmin"),
919  TrackEtaHistoPar_.getParameter<double>("Xmax"),
920  0.0,
921  0.0,
922  "g");
923  nHitsTIDDVsEtaH_ = ibook.bookProfile("nHitsTIDDVsEta",
924  "Number of Hits in TID Vs Eta (Double-sided)",
925  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
926  TrackEtaHistoPar_.getParameter<double>("Xmin"),
927  TrackEtaHistoPar_.getParameter<double>("Xmax"),
928  0.0,
929  0.0,
930  "g");
931  nHitsStripDVsEtaH_ = ibook.bookProfile("nHitsStripDVsEta",
932  "Number of Strip Hits Vs Eta (Double-sided)",
933  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
934  TrackEtaHistoPar_.getParameter<double>("Xmin"),
935  TrackEtaHistoPar_.getParameter<double>("Xmax"),
936  0.0,
937  0.0,
938  "g");
939 
940  nValidHitsVspTH_ = ibook.bookProfile("nValidHitsVspT",
941  "Number of Valid Hits Vs pT",
942  TrackPtHistoPar_.getParameter<int32_t>("Xbins"),
943  TrackPtHistoPar_.getParameter<double>("Xmin"),
944  TrackPtHistoPar_.getParameter<double>("Xmax"),
945  0.0,
946  0.0,
947  "g");
948  nValidHitsVsnVtxH_ =
949  ibook.bookProfile("nValidHitsVsnVtx", "Number of Valid Hits Vs Number of Vertex", 100, 0., 100., 0.0, 0.0, "g");
950  nValidHitsVsEtaH_ = ibook.bookProfile("nValidHitsVsEta",
951  "Number of Hits Vs Eta",
952  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
953  TrackEtaHistoPar_.getParameter<double>("Xmin"),
954  TrackEtaHistoPar_.getParameter<double>("Xmax"),
955  0.0,
956  0.0,
957  "g");
958 
959  nValidHitsVsCosThetaH_ =
960  ibook.bookProfile("nValidHitsVsCosTheta", "Number of Valid Hits Vs Cos(Theta)", 50, -1.0, 1.0, 0.0, 0.0, "g");
961  nValidHitsVsPhiH_ =
962  ibook.bookProfile("nValidHitsVsPhi", "Number of Valid Hits Vs Phi", 100, -3.5, 3.5, 0.0, 0.0, "g");
963 
964  nValidHitsPixVsEtaH_ = ibook.bookProfile("nValidHitsPixVsEta",
965  "Number of Valid Hits in Pixel Vs Eta",
966  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
967  TrackEtaHistoPar_.getParameter<double>("Xmin"),
968  TrackEtaHistoPar_.getParameter<double>("Xmax"),
969  0.0,
970  0.0,
971  "g");
972  nValidHitsPixBVsEtaH_ = ibook.bookProfile("nValidHitsPixBVsEta",
973  "Number of Valid Hits in Pixel Barrel Vs Eta",
974  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
975  TrackEtaHistoPar_.getParameter<double>("Xmin"),
976  TrackEtaHistoPar_.getParameter<double>("Xmax"),
977  0.0,
978  0.0,
979  "g");
980  nValidHitsPixEVsEtaH_ = ibook.bookProfile("nValidHitsPixEVsEta",
981  "Number of Valid Hits in Pixel Endcap Vs Eta",
982  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
983  TrackEtaHistoPar_.getParameter<double>("Xmin"),
984  TrackEtaHistoPar_.getParameter<double>("Xmax"),
985  0.0,
986  0.0,
987  "g");
988  nValidHitsStripVsEtaH_ = ibook.bookProfile("nValidHitsStripVsEta",
989  "Number of Valid Hits in SiStrip Vs Eta",
990  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
991  TrackEtaHistoPar_.getParameter<double>("Xmin"),
992  TrackEtaHistoPar_.getParameter<double>("Xmax"),
993  0.0,
994  0.0,
995  "g");
996  nValidHitsTIBVsEtaH_ = ibook.bookProfile("nValidHitsTIBVsEta",
997  "Number of Valid Hits in TIB Vs Eta",
998  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
999  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1000  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1001  0.0,
1002  0.0,
1003  "g");
1004  nValidHitsTOBVsEtaH_ = ibook.bookProfile("nValidHitsTOBVsEta",
1005  "Number of Valid Hits in TOB Vs Eta",
1006  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1007  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1008  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1009  0.0,
1010  0.0,
1011  "g");
1012  nValidHitsTECVsEtaH_ = ibook.bookProfile("nValidHitsTECVsEta",
1013  "Number of Valid Hits in TEC Vs Eta",
1014  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1015  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1016  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1017  0.0,
1018  0.0,
1019  "g");
1020  nValidHitsTIDVsEtaH_ = ibook.bookProfile("nValidHitsTIDVsEta",
1021  "Number of Valid Hits in TID Vs Eta",
1022  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1023  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1024  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1025  0.0,
1026  0.0,
1027  "g");
1028 
1029  nValidHitsPixVsPhiH_ = ibook.bookProfile("nValidHitsPixVsPhi",
1030  "Number of Valid Hits in Pixel Vs Phi",
1031  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1032  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1033  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1034  0.0,
1035  0.0,
1036  "g");
1037  nValidHitsPixBVsPhiH_ = ibook.bookProfile("nValidHitsPixBVsPhi",
1038  "Number of Valid Hits in Pixel Barrel Vs Phi",
1039  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1040  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1041  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1042  0.0,
1043  0.0,
1044  "g");
1045  nValidHitsPixEVsPhiH_ = ibook.bookProfile("nValidHitsPixEVsPhi",
1046  "Number of Valid Hits in Pixel Endcap Vs Phi",
1047  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1048  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1049  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1050  0.0,
1051  0.0,
1052  "g");
1053  nValidHitsStripVsPhiH_ = ibook.bookProfile("nValidHitsStripVsPhi",
1054  "Number of Valid Hits in SiStrip Vs Phi",
1055  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1056  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1057  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1058  0.0,
1059  0.0,
1060  "g");
1061  nValidHitsTIBVsPhiH_ = ibook.bookProfile("nValidHitsTIBVsPhi",
1062  "Number of Valid Hits in TIB Vs Phi",
1063  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1064  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1065  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1066  0.0,
1067  0.0,
1068  "g");
1069  nValidHitsTOBVsPhiH_ = ibook.bookProfile("nValidHitsTOBVsPhi",
1070  "Number of Valid Hits in TOB Vs Phi",
1071  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1072  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1073  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1074  0.0,
1075  0.0,
1076  "g");
1077  nValidHitsTECVsPhiH_ = ibook.bookProfile("nValidHitsTECVsPhi",
1078  "Number of Valid Hits in TEC Vs Phi",
1079  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1080  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1081  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1082  0.0,
1083  0.0,
1084  "g");
1085  nValidHitsTIDVsPhiH_ = ibook.bookProfile("nValidHitsTIDVsPhi",
1086  "Number of Valid Hits in TID Vs Phi",
1087  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1088  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1089  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1090  0.0,
1091  0.0,
1092  "g");
1093 
1094  nLostHitsPixVsEtaH_ = ibook.bookProfile("nLostHitsPixVsEta",
1095  "Number of Lost Hits in Pixel Vs Eta",
1096  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1097  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1098  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1099  0.0,
1100  0.0,
1101  "g");
1102  nLostHitsPixBVsEtaH_ = ibook.bookProfile("nLostHitsPixBVsEta",
1103  "Number of Lost Hits in Pixel Barrel Vs Eta",
1104  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1105  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1106  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1107  0.0,
1108  0.0,
1109  "g");
1110  nLostHitsPixEVsEtaH_ = ibook.bookProfile("nLostHitsPixEVsEta",
1111  "Number of Lost Hits in Pixel Endcap Vs Eta",
1112  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1113  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1114  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1115  0.0,
1116  0.0,
1117  "g");
1118  nLostHitsStripVsEtaH_ = ibook.bookProfile("nLostHitsStripVsEta",
1119  "Number of Lost Hits in SiStrip Vs Eta",
1120  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1121  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1122  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1123  0.0,
1124  0.0,
1125  "g");
1126  nLostHitsTIBVsEtaH_ = ibook.bookProfile("nLostHitsTIBVsEta",
1127  "Number of Lost Hits in TIB Vs Eta",
1128  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1129  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1130  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1131  0.0,
1132  0.0,
1133  "g");
1134  nLostHitsTOBVsEtaH_ = ibook.bookProfile("nLostHitsTOBVsEta",
1135  "Number of Lost Hits in TOB Vs Eta",
1136  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1137  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1138  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1139  0.0,
1140  0.0,
1141  "g");
1142  nLostHitsTECVsEtaH_ = ibook.bookProfile("nLostHitsTECVsEta",
1143  "Number of Lost Hits in TEC Vs Eta",
1144  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1145  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1146  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1147  0.0,
1148  0.0,
1149  "g");
1150  nLostHitsTIDVsEtaH_ = ibook.bookProfile("nLostHitsTIDVsEta",
1151  "Number of Lost Hits in TID Vs Eta",
1152  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1153  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1154  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1155  0.0,
1156  0.0,
1157  "g");
1158 
1159  nLostHitsPixVsPhiH_ = ibook.bookProfile("nLostHitsPixVsPhi",
1160  "Number of Lost Hits in Pixel Vs Phi",
1161  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1162  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1163  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1164  0.0,
1165  0.0,
1166  "g");
1167  nLostHitsPixBVsPhiH_ = ibook.bookProfile("nLostHitsPixBVsPhi",
1168  "Number of Lost Hits in Pixel Barrel Vs Phi",
1169  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1170  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1171  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1172  0.0,
1173  0.0,
1174  "g");
1175  nLostHitsPixEVsPhiH_ = ibook.bookProfile("nLostHitsPixEVsPhi",
1176  "Number of Lost Hits in Pixel Endcap Vs Phi",
1177  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1178  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1179  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1180  0.0,
1181  0.0,
1182  "g");
1183  nLostHitsStripVsPhiH_ = ibook.bookProfile("nLostHitsStripVsPhi",
1184  "Number of Lost Hits in SiStrip Vs Phi",
1185  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1186  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1187  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1188  0.0,
1189  0.0,
1190  "g");
1191  nLostHitsTIBVsPhiH_ = ibook.bookProfile("nLostHitsTIBVsPhi",
1192  "Number of Lost Hits in TIB Vs Phi",
1193  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1194  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1195  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1196  0.0,
1197  0.0,
1198  "g");
1199  nLostHitsTOBVsPhiH_ = ibook.bookProfile("nLostHitsTOBVsPhi",
1200  "Number of Lost Hits in TOB Vs Phi",
1201  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1202  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1203  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1204  0.0,
1205  0.0,
1206  "g");
1207  nLostHitsTECVsPhiH_ = ibook.bookProfile("nLostHitsTECVsPhi",
1208  "Number of Lost Hits in TEC Vs Phi",
1209  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1210  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1211  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1212  0.0,
1213  0.0,
1214  "g");
1215  nLostHitsTIDVsPhiH_ = ibook.bookProfile("nLostHitsTIDVsPhi",
1216  "Number of Lost Hits in TID Vs Phi",
1217  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1218  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1219  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1220  0.0,
1221  0.0,
1222  "g");
1223 
1224  nLostHitsPixVsIterationH_ = ibook.bookProfile(
1225  "nLostHitsPixVsIteration", "Number of Lost Hits in Pixel Vs Iteration", 47, -0.5, 46.5, 0.0, 0.0, "g");
1226  nLostHitsPixBVsIterationH_ = ibook.bookProfile(
1227  "nLostHitsPixBVsIteration", "Number of Lost Hits in Pixel Barrel Vs Iteration", 47, -0.5, 46.5, 0.0, 0.0, "g");
1228  nLostHitsPixEVsIterationH_ = ibook.bookProfile(
1229  "nLostHitsPixEVsIteration", "Number of Lost Hits in Pixel Endcap Vs Iteration", 47, -0.5, 46.5, 0.0, 0.0, "g");
1230  nLostHitsStripVsIterationH_ = ibook.bookProfile(
1231  "nLostHitsStripVsIteration", "Number of Lost Hits in SiStrip Vs Iteration", 47, -0.5, 46.5, 0.0, 0.0, "g");
1232  nLostHitsTIBVsIterationH_ = ibook.bookProfile(
1233  "nLostHitsTIBVsIteration", "Number of Lost Hits in TIB Vs Iteration", 47, -0.5, 46.5, 0.0, 0.0, "g");
1234  nLostHitsTOBVsIterationH_ = ibook.bookProfile(
1235  "nLostHitsTOBVsIteration", "Number of Lost Hits in TOB Vs Iteration", 47, -0.5, 46.5, 0.0, 0.0, "g");
1236  nLostHitsTECVsIterationH_ = ibook.bookProfile(
1237  "nLostHitsTECVsIteration", "Number of Lost Hits in TEC Vs Iteration", 47, -0.5, 46.5, 0.0, 0.0, "g");
1238  nLostHitsTIDVsIterationH_ = ibook.bookProfile(
1239  "nLostHitsTIDVsIteration", "Number of Lost Hits in TID Vs Iteration", 47, -0.5, 46.5, 0.0, 0.0, "g");
1240 
1241  trackChi2oNDFVsEtaH_ = ibook.bookProfile("trackChi2oNDFVsEta",
1242  "chi2/ndof of Tracks Vs Eta",
1243  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1244  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1245  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1246  0.0,
1247  10.0,
1248  "g");
1249  trackChi2oNDFVsPhiH_ = ibook.bookProfile("trackChi2oNDFVsPhi",
1250  "chi2/ndof of Tracks Vs Phi",
1251  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1252  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1253  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1254  0.0,
1255  10.0,
1256  "g");
1257 
1258  trackChi2probVsEtaH_ = ibook.bookProfile("trackChi2probVsEta",
1259  "chi2 probability of Tracks Vs Eta",
1260  TrackEtaHistoPar_.getParameter<int32_t>("Xbins"),
1261  TrackEtaHistoPar_.getParameter<double>("Xmin"),
1262  TrackEtaHistoPar_.getParameter<double>("Xmax"),
1263  0.0,
1264  1.0,
1265  "g");
1266  trackChi2probVsPhiH_ = ibook.bookProfile("trackChi2probVsPhi",
1267  "chi2 probability of Tracks Vs Phi",
1268  TrackPhiHistoPar_.getParameter<int32_t>("Xbins"),
1269  TrackPhiHistoPar_.getParameter<double>("Xmin"),
1270  TrackPhiHistoPar_.getParameter<double>("Xmax"),
1271  0.0,
1272  1.0,
1273  "g");
1274 
1275  trackIperr3dVsEtaH_ =
1276  ibook.bookProfile("trackIperr3dVsEta", "ip3d error of Tracks Vs Eta", 80, -4., 4., 0.0, 0.0, "g");
1277 
1278  trackSip2dVsEtaH_ = ibook.bookProfile("trackSip2dVsEta", "sip2d of Tracks Vs Eta", 80, -4., 4., 0.0, 0.0, "g");
1279 
1280  trackIperr3dVsEta2DH_ =
1281  ibook.book2D("trackIperr3dVsEta2D", "ip3d error of Tracks Vs Eta 2d", 80, -4., 4., 100, 0., 5.);
1282  trackIperr3dVsChi2prob2DH_ =
1283  ibook.book2D("trackIperr3dVsChi2prob2D", "ip3d error of Tracks Vs chi2prob 2d", 50, 0., 1., 100, 0., 5.);
1284  trackSip2dVsEta2DH_ = ibook.book2D("trackSip2dVsEta2D", "sip2d of Tracks Vs Eta 2d", 80, -4., 4., 200, -10., 10.);
1285  trackSip2dVsChi2prob2DH_ =
1286  ibook.book2D("trackSip2dVsChi2prob2D", "sip2d of Tracks Vs chi2prob 2d", 50, 0., 1., 200, -10., 10.);
1287 
1288  // On and off-track cluster properties
1289  hOnTrkClusChargeThinH_ = ibook.book1DD("hOnTrkClusChargeThin", "On-track Cluster Charge (Thin Sensor)", 100, 0, 1000);
1290  hOnTrkClusWidthThinH_ = ibook.book1DD("hOnTrkClusWidthThin", "On-track Cluster Width (Thin Sensor)", 20, -0.5, 19.5);
1291  hOnTrkClusChargeThickH_ =
1292  ibook.book1DD("hOnTrkClusChargeThick", "On-track Cluster Charge (Thick Sensor)", 100, 0, 1000);
1293  hOnTrkClusWidthThickH_ =
1294  ibook.book1DD("hOnTrkClusWidthThick", "On-track Cluster Width (Thick Sensor)", 20, -0.5, 19.5);
1295 
1296  hOffTrkClusChargeThinH_ =
1297  ibook.book1DD("hOffTrkClusChargeThin", "Off-track Cluster Charge (Thin Sensor)", 100, 0, 1000);
1298  hOffTrkClusWidthThinH_ =
1299  ibook.book1DD("hOffTrkClusWidthThin", "Off-track Cluster Width (Thin Sensor)", 20, -0.5, 19.5);
1300  hOffTrkClusChargeThickH_ =
1301  ibook.book1DD("hOffTrkClusChargeThick", "Off-track Cluster Charge (Thick Sensor)", 100, 0, 1000);
1302  hOffTrkClusWidthThickH_ =
1303  ibook.book1DD("hOffTrkClusWidthThick", "Off-track Cluster Width (Thick Sensor)", 20, -0.5, 19.5);
1304 }
1305 void StandaloneTrackMonitor::analyze(edm::Event const& iEvent, edm::EventSetup const& iSetup) {
1306  if (verbose_)
1307  edm::LogInfo("StandaloneTrackMonitor") << "Begin StandaloneTrackMonitor" << std::endl;
1308 
1309  nevt++;
1310 
1311  // Get event setup (to get global transformation)
1312  const TrackerGeometry& tkGeom = (*tkGeom_);
1313 
1314  // Primary vertex collection
1315  edm::Handle<reco::VertexCollection> vertexColl;
1316  iEvent.getByToken(vertexToken_, vertexColl);
1317  if (!vertexColl.isValid()) {
1318  edm::LogError("StandaloneTrackMonitor") << "Error! Failed to get reco::Vertex Collection, " << vertexTag_;
1319  return;
1320  }
1321  if (vertexColl->empty()) {
1322  edm::LogError("StandalonaTrackMonitor") << "No good vertex in the event!!" << std::endl;
1323  return;
1324  }
1325  const reco::Vertex& pv = (*vertexColl)[0];
1326 
1327  // Beam spot
1328  edm::Handle<reco::BeamSpot> beamSpot;
1329  iEvent.getByToken(bsToken_, beamSpot);
1330  if (!beamSpot.isValid()) {
1331  edm::LogError("StandalonaTrackMonitor") << "Beamspot for input tag: " << bsTag_ << " not found!!";
1332  return;
1333  }
1334 
1335  // Track collection
1336  edm::Handle<reco::TrackCollection> tracks;
1337  iEvent.getByToken(trackToken_, tracks);
1338  if (!tracks.isValid()) {
1339  edm::LogError("StandalonaTrackMonitor") << "TrackCollection for input tag: " << trackTag_ << " not found!!";
1340  return;
1341  }
1342 
1343  // Access PU information
1344  double wfac = 1.0; // for data
1345  if (!iEvent.isRealData()) {
1346  edm::Handle<std::vector<PileupSummaryInfo> > PupInfo;
1347  iEvent.getByToken(puSummaryToken_, PupInfo);
1348 
1349  if (verbose_)
1350  edm::LogInfo("StandaloneTrackMonitor") << "nPUColl = " << PupInfo->size();
1351  if (PupInfo.isValid()) {
1352  for (auto const& v : *PupInfo) {
1353  int bx = v.getBunchCrossing();
1354  if (bunchCrossingH_)
1355  bunchCrossingH_->Fill(bx);
1356  if (bx == 0) {
1357  if (nPUH_)
1358  nPUH_->Fill(v.getPU_NumInteractions());
1359  int ntrueInt = v.getTrueNumInteractions();
1360  int nVertex = (vertexColl.isValid() ? vertexColl->size() : 0);
1361  if (trueNIntH_)
1362  trueNIntH_->Fill(ntrueInt);
1363  if (doPUCorrection_) {
1364  if (nVertex > -1 && nVertex < int(vpu_.size()))
1365  wfac = vpu_.at(nVertex);
1366  else
1367  wfac = 0.0;
1368  }
1369  }
1370  }
1371  } else
1372  edm::LogError("StandalonaTrackMonitor") << "PUSummary for input tag: " << puSummaryTag_ << " not found!!";
1373  if (doTrackCorrection_) {
1374  int ntrack = 0;
1375  for (auto const& track : *tracks) {
1376  if (!track.quality(reco::Track::qualityByName(trackQuality_)))
1377  continue;
1378  ++ntrack;
1379  }
1380  if (ntrack > -1 && ntrack < int(vtrack_.size()))
1381  wfac = vtrack_.at(ntrack);
1382  else
1383  wfac = 0.0;
1384  }
1385  }
1386  if (verbose_)
1387  edm::LogInfo("StandaloneTrackMonitor") << "PU reweight factor = " << wfac;
1388 
1389  if (haveAllHistograms_) {
1390  int nvtx = (vertexColl.isValid() ? vertexColl->size() : 0);
1391  nVertexH_->Fill(nvtx, wfac);
1392  nVtxH_->Fill(nvtx);
1393  }
1394 
1395  // Get MVA and quality mask collections
1396  int ntracks = 0;
1397  std::vector<TLorentzVector> list;
1398 
1399  const TransientTrackBuilder* theB = &iSetup.getData(transTrackToken_);
1400  TransientVertex mumuTransientVtx;
1401  std::vector<reco::TransientTrack> tks;
1402 
1403  if (tracks.isValid()) {
1404  if (verbose_)
1405  edm::LogInfo("StandaloneTrackMonitor") << "Total # of Tracks: " << tracks->size();
1406  reco::Track::TrackQuality quality = reco::Track::qualityByName(trackQuality_);
1407 
1408  for (auto const& track : *tracks) {
1409  if (!track.quality(quality))
1410  continue;
1411  ++ntracks;
1412 
1413  reco::TransientTrack trajectory = theB->build(track);
1414  tks.push_back(trajectory);
1415 
1416  double eta = track.eta();
1417  double theta = track.theta();
1418  double phi = track.phi();
1419  double pt = track.pt();
1420 
1421  const reco::HitPattern& hitp = track.hitPattern();
1422  double nAllHits = hitp.numberOfAllHits(reco::HitPattern::TRACK_HITS);
1423  double nAllTrackerHits = hitp.numberOfAllTrackerHits(reco::HitPattern::TRACK_HITS);
1424 
1425  double trackdeltaR = std::numeric_limits<double>::max();
1426  double muonmass = 0.1056583745;
1427  TLorentzVector trackpmu;
1428  trackpmu.SetPtEtaPhiM(track.pt(), track.eta(), track.phi(), muonmass);
1429  for (auto const& TRACK : *tracks) {
1430  if (&track == &TRACK)
1431  continue;
1432  double tmpdr = reco::deltaR(track.eta(), track.phi(), TRACK.eta(), TRACK.phi());
1433  if (tmpdr < trackdeltaR)
1434  trackdeltaR = tmpdr;
1435  }
1436 
1437  list.push_back(trackpmu);
1438 
1439  double nValidTrackerHits = hitp.numberOfValidTrackerHits();
1440  double nValidPixelHits = hitp.numberOfValidPixelHits();
1441  double nValidPixelBHits = hitp.numberOfValidPixelBarrelHits();
1442  double nValidPixelEHits = hitp.numberOfValidPixelEndcapHits();
1443  double nValidStripHits = hitp.numberOfValidStripHits();
1444  double nValidTIBHits = hitp.numberOfValidStripTIBHits();
1445  double nValidTOBHits = hitp.numberOfValidStripTOBHits();
1446  double nValidTIDHits = hitp.numberOfValidStripTIDHits();
1447  double nValidTECHits = hitp.numberOfValidStripTECHits();
1448 
1449  int missingInnerHit = hitp.numberOfAllHits(reco::HitPattern::MISSING_INNER_HITS);
1450  int missingOuterHit = hitp.numberOfAllHits(reco::HitPattern::MISSING_OUTER_HITS);
1451 
1452  nValidHitsVspTH_->Fill(pt, nValidTrackerHits);
1453  nValidHitsVsEtaH_->Fill(eta, nValidTrackerHits);
1454  nValidHitsVsCosThetaH_->Fill(std::cos(theta), nValidTrackerHits);
1455  nValidHitsVsPhiH_->Fill(phi, nValidTrackerHits);
1456  nValidHitsVsnVtxH_->Fill(vertexColl->size(), nValidTrackerHits);
1457 
1458  nValidHitsPixVsEtaH_->Fill(eta, nValidPixelHits);
1459  nValidHitsPixBVsEtaH_->Fill(eta, nValidPixelBHits);
1460  nValidHitsPixEVsEtaH_->Fill(eta, nValidPixelEHits);
1461  nValidHitsStripVsEtaH_->Fill(eta, nValidStripHits);
1462  nValidHitsTIBVsEtaH_->Fill(eta, nValidTIBHits);
1463  nValidHitsTOBVsEtaH_->Fill(eta, nValidTOBHits);
1464  nValidHitsTECVsEtaH_->Fill(eta, nValidTECHits);
1465  nValidHitsTIDVsEtaH_->Fill(eta, nValidTIDHits);
1466 
1467  nValidHitsPixVsPhiH_->Fill(phi, nValidPixelHits);
1468  nValidHitsPixBVsPhiH_->Fill(phi, nValidPixelBHits);
1469  nValidHitsPixEVsPhiH_->Fill(phi, nValidPixelEHits);
1470  nValidHitsStripVsPhiH_->Fill(phi, nValidStripHits);
1471  nValidHitsTIBVsPhiH_->Fill(phi, nValidTIBHits);
1472  nValidHitsTOBVsPhiH_->Fill(phi, nValidTOBHits);
1473  nValidHitsTECVsPhiH_->Fill(phi, nValidTECHits);
1474  nValidHitsTIDVsPhiH_->Fill(phi, nValidTIDHits);
1475 
1476  int nLostHits = track.numberOfLostHits();
1477  int nMissingExpectedInnerHits = hitp.numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS);
1478  int nMissingExpectedOuterHits = hitp.numberOfLostHits(reco::HitPattern::MISSING_OUTER_HITS);
1479  int nLostTrackerHits = hitp.numberOfLostTrackerHits(reco::HitPattern::TRACK_HITS);
1480  int nLostPixHits = hitp.numberOfLostPixelHits(reco::HitPattern::TRACK_HITS);
1481  int nLostPixBHits = hitp.numberOfLostPixelBarrelHits(reco::HitPattern::TRACK_HITS);
1482  int nLostPixEHits = hitp.numberOfLostPixelEndcapHits(reco::HitPattern::TRACK_HITS);
1483  int nLostStripHits = hitp.numberOfLostStripHits(reco::HitPattern::TRACK_HITS);
1484  int nLostStripTIBHits = hitp.numberOfLostStripTIBHits(reco::HitPattern::TRACK_HITS);
1485  int nLostStripTIDHits = hitp.numberOfLostStripTIDHits(reco::HitPattern::TRACK_HITS);
1486  int nLostStripTOBHits = hitp.numberOfLostStripTOBHits(reco::HitPattern::TRACK_HITS);
1487  int nLostStripTECHits = hitp.numberOfLostStripTECHits(reco::HitPattern::TRACK_HITS);
1488  int nIteration = track.originalAlgo();
1489 
1490  nLostHitsVspTH_->Fill(pt, nLostTrackerHits);
1491  nLostHitsVsEtaH_->Fill(eta, nLostTrackerHits);
1492  nLostHitsVsCosThetaH_->Fill(std::cos(theta), nLostTrackerHits);
1493  nLostHitsVsPhiH_->Fill(phi, nLostTrackerHits);
1494  nLostHitsVsIterationH_->Fill(nIteration, nLostTrackerHits);
1495 
1496  nLostHitsPixVsEtaH_->Fill(eta, nLostPixHits);
1497  nLostHitsPixBVsEtaH_->Fill(eta, nLostPixBHits);
1498  nLostHitsPixEVsEtaH_->Fill(eta, nLostPixEHits);
1499  nLostHitsStripVsEtaH_->Fill(eta, nLostStripHits);
1500  nLostHitsTIBVsEtaH_->Fill(eta, nLostStripTIBHits);
1501  nLostHitsTOBVsEtaH_->Fill(eta, nLostStripTOBHits);
1502  nLostHitsTECVsEtaH_->Fill(eta, nLostStripTECHits);
1503  nLostHitsTIDVsEtaH_->Fill(eta, nLostStripTIDHits);
1504 
1505  nLostHitsPixVsPhiH_->Fill(phi, nLostPixHits);
1506  nLostHitsPixBVsPhiH_->Fill(phi, nLostPixBHits);
1507  nLostHitsPixEVsPhiH_->Fill(phi, nLostPixEHits);
1508  nLostHitsStripVsPhiH_->Fill(phi, nLostStripHits);
1509  nLostHitsTIBVsPhiH_->Fill(phi, nLostStripTIBHits);
1510  nLostHitsTOBVsPhiH_->Fill(phi, nLostStripTOBHits);
1511  nLostHitsTECVsPhiH_->Fill(phi, nLostStripTECHits);
1512  nLostHitsTIDVsPhiH_->Fill(phi, nLostStripTIDHits);
1513 
1514  nLostHitsPixVsIterationH_->Fill(nIteration, nLostPixHits);
1515  nLostHitsPixBVsIterationH_->Fill(nIteration, nLostPixBHits);
1516  nLostHitsPixEVsIterationH_->Fill(nIteration, nLostPixEHits);
1517  nLostHitsStripVsIterationH_->Fill(nIteration, nLostStripHits);
1518  nLostHitsTIBVsIterationH_->Fill(nIteration, nLostStripTIBHits);
1519  nLostHitsTOBVsIterationH_->Fill(nIteration, nLostStripTOBHits);
1520  nLostHitsTECVsIterationH_->Fill(nIteration, nLostStripTECHits);
1521  nLostHitsTIDVsIterationH_->Fill(nIteration, nLostStripTIDHits);
1522 
1523  if (abs(eta) <= 1.4) {
1524  nMissingInnerHitBH_->Fill(missingInnerHit, wfac);
1525  nMissingOuterHitBH_->Fill(missingOuterHit, wfac);
1526  } else {
1527  nMissingInnerHitEH_->Fill(missingInnerHit, wfac);
1528  nMissingOuterHitEH_->Fill(missingOuterHit, wfac);
1529  }
1530 
1531  for (int i = 0; i < hitp.numberOfAllHits(reco::HitPattern::TRACK_HITS); i++) {
1532  uint32_t hit = hitp.getHitPattern(reco::HitPattern::TRACK_HITS, i);
1533  if (hitp.missingHitFilter(hit)) {
1534  double losthitBylayer = -1.0;
1535  double losthitBylayerPix = -1.0;
1536  double losthitBylayerStrip = -1.0;
1537  int layer = hitp.getLayer(hit);
1538  if (hitp.pixelBarrelHitFilter(hit)) {
1539  losthitBylayer = layer;
1540  losthitBylayerPix = layer;
1541  } else if (hitp.pixelEndcapHitFilter(hit)) {
1542  losthitBylayer = layer + 4;
1543  losthitBylayerPix = layer + 4;
1544  } else if (hitp.stripTIBHitFilter(hit)) {
1545  losthitBylayer = layer + 7;
1546  losthitBylayerStrip = layer;
1547  } else if (hitp.stripTIDHitFilter(hit)) {
1548  losthitBylayer = layer + 11;
1549  losthitBylayerStrip = layer + 4;
1550  } else if (hitp.stripTOBHitFilter(hit)) {
1551  losthitBylayer = layer + 14;
1552  losthitBylayerStrip = layer + 7;
1553  } else if (hitp.stripTECHitFilter(hit)) {
1554  losthitBylayer = layer + 20;
1555  losthitBylayerStrip = layer + 13;
1556  }
1557  if (losthitBylayer > -1)
1558  nLostHitByLayerH_->Fill(losthitBylayer, wfac);
1559  if (losthitBylayerPix > -1)
1560  nLostHitByLayerPixH_->Fill(losthitBylayerPix, wfac);
1561  if (losthitBylayerStrip > -1)
1562  nLostHitByLayerStripH_->Fill(losthitBylayerStrip, wfac);
1563  }
1564  }
1565 
1566  if (haveAllHistograms_) {
1567  double etaError = track.etaError();
1568  double thetaError = track.thetaError();
1569  double phiError = track.phiError();
1570  double p = track.p();
1571  double ptError = track.ptError();
1572  double qoverp = track.qoverp();
1573  double qoverpError = track.qoverpError();
1574  double charge = track.charge();
1575 
1576  double dxy = track.dxy(beamSpot->position());
1577  double dxyError = track.dxyError();
1578  double dz = track.dz(beamSpot->position());
1579  double dzError = track.dzError();
1580 
1581  double trkd0 = track.d0();
1582  double chi2 = track.chi2();
1583  double ndof = track.ndof();
1584  double chi2prob = TMath::Prob(track.chi2(), (int)track.ndof());
1585  double chi2oNDF = track.normalizedChi2();
1586  double vx = track.vx();
1587  double vy = track.vy();
1588  double vz = track.vz();
1590  unsigned int track_algo = track.algo();
1591  unsigned int track_origalgo = track.originalAlgo();
1592  // stopping source
1593  int ssmax = trackStoppingSourceH_->getNbinsX();
1594  double stop = track.stopReason() > ssmax ? double(ssmax - 1) : static_cast<double>(track.stopReason());
1595  double distanceOfClosestApproachToPV = track.dxy(pv.position());
1596  double xPointOfClosestApproachwrtPV = track.vx() - pv.position().x();
1597  double yPointOfClosestApproachwrtPV = track.vy() - pv.position().y();
1598  double positionZ0 = track.dz(pv.position());
1599 
1600  reco::TransientTrack transTrack = iSetup.get<TransientTrackRecord>().get(transTrackToken_).build(track);
1601 
1602  double ip3dToPV = 0, iperr3dToPV = 0, sip3dToPV = 0, sip2dToPV = 0;
1603  GlobalVector dir(track.px(), track.py(), track.pz());
1604  std::pair<bool, Measurement1D> ip3d = IPTools::signedImpactParameter3D(transTrack, dir, pv);
1605  std::pair<bool, Measurement1D> ip2d = IPTools::signedTransverseImpactParameter(transTrack, dir, pv);
1606  if (ip3d.first) {
1607  sip3dToPV = ip3d.second.value() / ip3d.second.error();
1608  ip3dToPV = ip3d.second.value();
1609  iperr3dToPV = ip3d.second.error();
1610  }
1611 
1612  double ip3dToBS = 0, iperr3dToBS = 0, sip3dToBS = 0, sip2dToBS = 0;
1613  reco::Vertex beamspotvertex((*beamSpot).position(), (*beamSpot).covariance3D());
1614  std::pair<bool, Measurement1D> ip3dbs = IPTools::signedImpactParameter3D(transTrack, dir, beamspotvertex);
1615  std::pair<bool, Measurement1D> ip2dbs =
1616  IPTools::signedTransverseImpactParameter(transTrack, dir, beamspotvertex);
1617  if (ip3dbs.first) {
1618  sip3dToBS = ip3dbs.second.value() / ip3dbs.second.error();
1619  ip3dToBS = ip3dbs.second.value();
1620  iperr3dToBS = ip3dbs.second.error();
1621  }
1622 
1623  double ip2dToPV = 0, iperr2dToPV = 0;
1624  if (ip2d.first) {
1625  ip2dToPV = ip2d.second.value();
1626  iperr2dToPV = ip2d.second.error();
1627  sip2dToPV = ip2d.second.value() / ip2d.second.error();
1628  }
1629 
1630  double ip2dToBS = 0, iperr2dToBS = 0;
1631  if (ip2dbs.first) {
1632  ip2dToBS = ip2dbs.second.value();
1633  iperr2dToBS = ip2dbs.second.error();
1634  sip2dToBS = ip2d.second.value() / ip2d.second.error();
1635  }
1636 
1637  if (ip2d.first)
1638  sip2dToPV = ip2d.second.value() / ip2d.second.error();
1639  double sipDxyToPV = track.dxy(pv.position()) / track.dxyError();
1640  double sipDzToPV = track.dz(pv.position()) / track.dzError();
1641 
1642  // Fill the histograms
1643  trackDeltaRwrtClosestTrack_->Fill(trackdeltaR, wfac);
1644  trackEtaH_->Fill(eta, wfac);
1645  trackEtaerrH_->Fill(etaError, wfac);
1646  trackCosThetaH_->Fill(std::cos(theta), wfac);
1647  trackThetaerrH_->Fill(thetaError, wfac);
1648  trackPhiH_->Fill(phi, wfac);
1649  trackPhierrH_->Fill(phiError, wfac);
1650  trackPH_->Fill(p, wfac);
1651  trackPtH_->Fill(pt, wfac);
1652  trackPt_ZoomH_->Fill(pt, wfac);
1653  trackPterrH_->Fill(ptError, wfac);
1654  trackqOverpH_->Fill(qoverp, wfac);
1655  trackqOverperrH_->Fill(qoverpError, wfac);
1656  trackChargeH_->Fill(charge, wfac);
1657  trackChi2H_->Fill(chi2, wfac);
1658  trackChi2ProbH_->Fill(chi2prob, wfac);
1659  trackChi2oNDFH_->Fill(chi2oNDF, wfac);
1660  trackd0H_->Fill(trkd0, wfac);
1661  tracknDOFH_->Fill(ndof, wfac);
1662  trackChi2bynDOFH_->Fill(chi2 / ndof, wfac);
1663  trackalgoH_->Fill(track_algo, wfac);
1664  trackorigalgoH_->Fill(track_origalgo, wfac);
1665  trackStoppingSourceH_->Fill(stop, wfac);
1666  trackChi2oNDFVsEtaH_->Fill(eta, chi2oNDF);
1667  trackChi2oNDFVsPhiH_->Fill(phi, chi2oNDF);
1668  trackChi2probVsEtaH_->Fill(eta, chi2prob);
1669  trackChi2probVsPhiH_->Fill(phi, chi2prob);
1670 
1671  nlostHitsH_->Fill(nLostHits, wfac);
1672  nMissingExpectedInnerHitsH_->Fill(nMissingExpectedInnerHits, wfac);
1673  nMissingExpectedOuterHitsH_->Fill(nMissingExpectedOuterHits, wfac);
1674  nlostTrackerHitsH_->Fill(nLostTrackerHits, wfac);
1675 
1676  beamSpotXYposH_->Fill(dxy, wfac);
1677  beamSpotXYposerrH_->Fill(dxyError, wfac);
1678  beamSpotZposH_->Fill(dz, wfac);
1679  beamSpotZposerrH_->Fill(dzError, wfac);
1680 
1681  vertexXposH_->Fill(vx, wfac);
1682  vertexYposH_->Fill(vy, wfac);
1683  vertexZposH_->Fill(vz, wfac);
1684 
1685  int nPixBarrel = 0, nPixEndcap = 0, nStripTIB = 0, nStripTOB = 0, nStripTEC = 0, nStripTID = 0;
1686  if (isRECO_) {
1687  for (auto it = track.recHitsBegin(); it != track.recHitsEnd(); ++it) {
1688  const TrackingRecHit& hit = (**it);
1689  if (hit.isValid()) {
1690  if (hit.geographicalId().det() == DetId::Tracker) {
1691  int subdetId = hit.geographicalId().subdetId();
1692  if (subdetId == PixelSubdetector::PixelBarrel)
1693  ++nPixBarrel;
1694  else if (subdetId == PixelSubdetector::PixelEndcap)
1695  ++nPixEndcap;
1696  else if (subdetId == StripSubdetector::TIB)
1697  ++nStripTIB;
1698  else if (subdetId == StripSubdetector::TOB)
1699  ++nStripTOB;
1700  else if (subdetId == StripSubdetector::TEC)
1701  ++nStripTEC;
1702  else if (subdetId == StripSubdetector::TID)
1703  ++nStripTID;
1704 
1705  // Find on-track clusters
1706  processHit(hit, iSetup, tkGeom, wfac);
1707  }
1708  }
1709  }
1710  }
1711  if (verbose_)
1712  edm::LogInfo("StandaloneTrackMonitor")
1713  << " >>> HITs: nPixBarrel: " << nPixBarrel << " nPixEndcap: " << nPixEndcap << " nStripTIB: " << nStripTIB
1714  << " nStripTOB: " << nStripTOB << " nStripTEC: " << nStripTEC << " nStripTID: " << nStripTID;
1715  if (haveAllHistograms_) {
1716  nPixBarrelH_->Fill(nPixBarrel, wfac);
1717  nPixEndcapH_->Fill(nPixEndcap, wfac);
1718  nStripTIBH_->Fill(nStripTIB, wfac);
1719  nStripTOBH_->Fill(nStripTOB, wfac);
1720  nStripTECH_->Fill(nStripTEC, wfac);
1721  nStripTIDH_->Fill(nStripTID, wfac);
1722  }
1723 
1724  DistanceOfClosestApproachToPVH_->Fill(distanceOfClosestApproachToPV, wfac);
1725  DistanceOfClosestApproachToPVZoomedH_->Fill(distanceOfClosestApproachToPV, wfac);
1726  DistanceOfClosestApproachToPVVsPhiH_->Fill(phi, distanceOfClosestApproachToPV);
1727  xPointOfClosestApproachVsZ0wrtPVH_->Fill(positionZ0, xPointOfClosestApproachwrtPV);
1728  yPointOfClosestApproachVsZ0wrtPVH_->Fill(positionZ0, yPointOfClosestApproachwrtPV);
1729 
1730  ip3dToPVH_->Fill(ip3dToPV, wfac);
1731  iperr3dToPVH_->Fill(iperr3dToPV, wfac);
1732  ip3dToBSH_->Fill(ip3dToBS, wfac);
1733  iperr3dToBSH_->Fill(iperr3dToBS, wfac);
1734  ip2dToPVH_->Fill(ip2dToPV, wfac);
1735  iperr2dToPVH_->Fill(iperr2dToPV, wfac);
1736  iperr2dToPVWtH_->Fill(iperr2dToPV, wfac);
1737  ip2dToBSH_->Fill(ip2dToBS, wfac);
1738  iperr2dToBSH_->Fill(iperr2dToBS, wfac);
1739 
1740  iperr3dToPVWtH_->Fill(iperr3dToPV, wfac);
1741  sip3dToPVH_->Fill(sip3dToPV, wfac);
1742  sip2dToPVH_->Fill(sip2dToPV, wfac);
1743  sip3dToBSH_->Fill(sip3dToBS, wfac);
1744  sip2dToBSH_->Fill(sip2dToBS, wfac);
1745  sip2dToPVWtH_->Fill(sip2dToPV, wfac);
1746  sipDxyToPVH_->Fill(sipDxyToPV, wfac);
1747  sipDzToPVH_->Fill(sipDzToPV, wfac);
1748 
1749  if (nValidPixelHits >= 2) {
1750  ip3dToPV2validpixelhitsH_->Fill(ip3dToPV, wfac);
1751  iperr3dToPV2validpixelhitsH_->Fill(iperr3dToPV, wfac);
1752  ip3dToBS2validpixelhitsH_->Fill(ip3dToBS, wfac);
1753  iperr3dToBS2validpixelhitsH_->Fill(iperr3dToBS, wfac);
1754  sip3dToPV2validpixelhitsH_->Fill(sip3dToPV, wfac);
1755  sip3dToBS2validpixelhitsH_->Fill(sip3dToBS, wfac);
1756  }
1757 
1758  trackIperr3dVsEta2DH_->Fill(eta, iperr3dToPV, wfac);
1759  trackSip2dVsEta2DH_->Fill(eta, sip2dToPV, wfac);
1760  trackIperr3dVsChi2prob2DH_->Fill(chi2prob, iperr3dToPV, wfac);
1761  trackSip2dVsChi2prob2DH_->Fill(chi2prob, sip2dToPV, wfac);
1762 
1763  trackIperr3dVsEtaH_->Fill(eta, iperr3dToPV);
1764  trackSip2dVsEtaH_->Fill(eta, sip2dToPV);
1765 
1766  double trackerLayersWithMeasurement = hitp.trackerLayersWithMeasurement();
1767  double pixelLayersWithMeasurement = hitp.pixelLayersWithMeasurement();
1768  double pixelBLayersWithMeasurement = hitp.pixelBarrelLayersWithMeasurement();
1769  double pixelELayersWithMeasurement = hitp.pixelEndcapLayersWithMeasurement();
1770  double stripLayersWithMeasurement = hitp.stripLayersWithMeasurement();
1771  double stripTIBLayersWithMeasurement = hitp.stripTIBLayersWithMeasurement();
1772  double stripTOBLayersWithMeasurement = hitp.stripTOBLayersWithMeasurement();
1773  double stripTIDLayersWithMeasurement = hitp.stripTIDLayersWithMeasurement();
1774  double stripTECLayersWithMeasurement = hitp.stripTECLayersWithMeasurement();
1775 
1776  trkLayerwithMeasurementH_->Fill(trackerLayersWithMeasurement, wfac);
1777  pixelLayerwithMeasurementH_->Fill(pixelLayersWithMeasurement, wfac);
1778  pixelBLayerwithMeasurementH_->Fill(pixelBLayersWithMeasurement, wfac);
1779  pixelELayerwithMeasurementH_->Fill(pixelELayersWithMeasurement, wfac);
1780  stripLayerwithMeasurementH_->Fill(stripLayersWithMeasurement, wfac);
1781  stripTIBLayerwithMeasurementH_->Fill(stripTIBLayersWithMeasurement, wfac);
1782  stripTOBLayerwithMeasurementH_->Fill(stripTOBLayersWithMeasurement, wfac);
1783  stripTIDLayerwithMeasurementH_->Fill(stripTIDLayersWithMeasurement, wfac);
1784  stripTECLayerwithMeasurementH_->Fill(stripTECLayersWithMeasurement, wfac);
1785 
1786  nallHitsH_->Fill(nAllHits, wfac);
1787  ntrackerHitsH_->Fill(nAllTrackerHits, wfac);
1788  nvalidTrackerHitsH_->Fill(nValidTrackerHits, wfac);
1789  nvalidPixelHitsH_->Fill(nValidPixelHits, wfac);
1790  nvalidPixelBHitsH_->Fill(nValidPixelBHits, wfac);
1791  nvalidPixelEHitsH_->Fill(nValidPixelEHits, wfac);
1792  nvalidStripHitsH_->Fill(nValidStripHits, wfac);
1793  nvalidTIBHitsH_->Fill(nValidTIBHits, wfac);
1794  nvalidTOBHitsH_->Fill(nValidTOBHits, wfac);
1795  nvalidTIDHitsH_->Fill(nValidTIDHits, wfac);
1796  nvalidTECHitsH_->Fill(nValidTECHits, wfac);
1797 
1798  nlostTrackerHitsH_->Fill(nLostTrackerHits, wfac);
1799  nlostPixelHitsH_->Fill(nLostPixHits, wfac);
1800  nlostPixelBHitsH_->Fill(nLostPixBHits, wfac);
1801  nlostPixelEHitsH_->Fill(nLostPixEHits, wfac);
1802  nlostStripHitsH_->Fill(nLostStripHits, wfac);
1803  nlostTIBHitsH_->Fill(nLostStripTIBHits, wfac);
1804  nlostTOBHitsH_->Fill(nLostStripTOBHits, wfac);
1805  nlostTIDHitsH_->Fill(nLostStripTIDHits, wfac);
1806  nlostTECHitsH_->Fill(nLostStripTECHits, wfac);
1807  }
1808  }
1809  } else {
1810  edm::LogError("StandaloneTrackMonitor") << "Error! Failed to get reco::Track collection, " << trackTag_;
1811  }
1812 
1813  if (haveAllHistograms_) {
1814  nTracksH_->Fill(ntracks, wfac);
1815  edm::Handle<std::vector<reco::PFJet> > jetsColl;
1816  iEvent.getByToken(jetsToken_, jetsColl);
1817  nJet_->Fill(jetsColl->size());
1818 
1819  for (auto const& jet : *jetsColl) {
1820  Jet_pt_->Fill(jet.pt(), wfac);
1821  Jet_eta_->Fill(jet.eta(), wfac);
1822  Jet_energy_->Fill(jet.energy(), wfac);
1823  Jet_chargedMultiplicity_->Fill(jet.chargedHadronMultiplicity(), wfac);
1824  }
1825  if (list.size() >= 2) {
1826  TLorentzVector Zp = list[0] + list[1];
1827  Zpt_->Fill(Zp.Pt(), wfac);
1828  ZInvMass_->Fill(Zp.Mag(), wfac);
1829  KalmanVertexFitter kalman(true);
1830  mumuTransientVtx = kalman.vertex(tks);
1831  if (mumuTransientVtx.isValid()) {
1832  const reco::Vertex* theClosestVertex;
1833  edm::Handle<reco::VertexCollection> vertexHandle = iEvent.getHandle(vertexToken_);
1834  if (vertexHandle.isValid()) {
1835  const reco::VertexCollection* vertices = vertexHandle.product();
1836  theClosestVertex = this->findClosestVertex(mumuTransientVtx, vertices);
1837  reco::Vertex theMainVtx;
1838  if (theClosestVertex == nullptr) {
1839  theMainVtx = vertexHandle.product()->front();
1840  } else {
1841  theMainVtx = *theClosestVertex;
1842  }
1843  if (theMainVtx.isValid()) {
1844  const auto& mainVertexPos = theMainVtx.position();
1845  const auto& myVertexPos = mumuTransientVtx.position();
1846 
1847  const double deltaX = myVertexPos.x() - mainVertexPos.x();
1848  const double deltaY = myVertexPos.y() - mainVertexPos.y();
1849  const double deltaZ = myVertexPos.z() - mainVertexPos.z();
1850 
1851  const double deltaNorm = sqrt(deltaX * deltaX + deltaY * deltaY + deltaZ * deltaZ);
1852  const double zpNorm = sqrt(Zp.Px() * Zp.Px() + Zp.Py() * Zp.Py() + Zp.Pz() * Zp.Pz());
1853 
1854  const double cosPhi3D = (Zp.Px() * deltaX + Zp.Py() * deltaY + Zp.Pz() * deltaZ) / (zpNorm * deltaNorm);
1855  cosPhi3DdileptonH_->Fill(cosPhi3D, wfac);
1856  }
1857  }
1858  }
1859  }
1860  }
1861 
1862  // off track cluster properties (only on RECO data-tier)
1863  if (isRECO_) {
1864  processClusters(iEvent, iSetup, tkGeom, wfac);
1865  }
1866 
1867  if (verbose_)
1868  edm::LogInfo("StandaloneTrackMonitor") << "Ends StandaloneTrackMonitor successfully";
1869 }
1870 void StandaloneTrackMonitor::processClusters(edm::Event const& iEvent,
1871  edm::EventSetup const& iSetup,
1872  const TrackerGeometry& tkGeom,
1873  double wfac) {
1874  // SiStripClusters
1875  edm::Handle<edmNew::DetSetVector<SiStripCluster> > clusterHandle;
1876  iEvent.getByToken(clusterToken_, clusterHandle);
1877 
1878  if (clusterHandle.isValid()) {
1879  // Loop on Dets
1880  for (edmNew::DetSetVector<SiStripCluster>::const_iterator dsvit = clusterHandle->begin();
1881  dsvit != clusterHandle->end();
1882  ++dsvit) {
1883  uint32_t detId = dsvit->id();
1884  std::map<uint32_t, std::set<const SiStripCluster*> >::iterator jt = clusterMap_.find(detId);
1885  bool detid_found = (jt != clusterMap_.end()) ? true : false;
1886 
1887  // Loop on Clusters
1888  for (edmNew::DetSet<SiStripCluster>::const_iterator clusit = dsvit->begin(); clusit != dsvit->end(); ++clusit) {
1889  if (detid_found) {
1890  std::set<const SiStripCluster*>& s = jt->second;
1891  if (s.find(&*clusit) != s.end())
1892  continue;
1893  }
1894 
1895  siStripClusterInfo_.setCluster(*clusit, detId);
1896  float charge = siStripClusterInfo_.charge();
1897  float width = siStripClusterInfo_.width();
1898 
1899  const GeomDetUnit* detUnit = tkGeom.idToDetUnit(detId);
1900  float thickness = detUnit->surface().bounds().thickness(); // unit cm
1901  if (thickness > 0.035) {
1902  hOffTrkClusChargeThickH_->Fill(charge, wfac);
1903  hOffTrkClusWidthThickH_->Fill(width, wfac);
1904  } else {
1905  hOffTrkClusChargeThinH_->Fill(charge, wfac);
1906  hOffTrkClusWidthThinH_->Fill(width, wfac);
1907  }
1908  }
1909  }
1910  } else {
1911  edm::LogError("StandaloneTrackMonitor") << "ClusterCollection " << clusterTag_ << " not valid!!" << std::endl;
1912  }
1913 }
1914 void StandaloneTrackMonitor::processHit(const TrackingRecHit& recHit,
1915  edm::EventSetup const& iSetup,
1916  const TrackerGeometry& tkGeom,
1917  double wfac) {
1918  uint32_t detid = recHit.geographicalId();
1919  const GeomDetUnit* detUnit = tkGeom.idToDetUnit(detid);
1920  float thickness = detUnit->surface().bounds().thickness(); // unit cm
1921 
1922  auto const& thit = static_cast<BaseTrackerRecHit const&>(recHit);
1923  if (!thit.isValid())
1924  return;
1925 
1926  auto const& clus = thit.firstClusterRef();
1927  if (!clus.isValid())
1928  return;
1929  if (!clus.isStrip())
1930  return;
1931 
1932  if (thit.isMatched()) {
1933  const SiStripMatchedRecHit2D& matchedHit = dynamic_cast<const SiStripMatchedRecHit2D&>(recHit);
1934 
1935  auto& clusterM = matchedHit.monoCluster();
1936  siStripClusterInfo_.setCluster(clusterM, detid);
1937 
1938  if (thickness > 0.035) {
1939  hOnTrkClusChargeThickH_->Fill(siStripClusterInfo_.charge(), wfac);
1940  hOnTrkClusWidthThickH_->Fill(siStripClusterInfo_.width(), wfac);
1941  } else {
1942  hOnTrkClusChargeThinH_->Fill(siStripClusterInfo_.charge(), wfac);
1943  hOnTrkClusWidthThinH_->Fill(siStripClusterInfo_.width(), wfac);
1944  }
1945  addClusterToMap(detid, &clusterM);
1946 
1947  auto& clusterS = matchedHit.stereoCluster();
1948 
1949  siStripClusterInfo_.setCluster(clusterS, detid);
1950  if (thickness > 0.035) {
1951  hOnTrkClusChargeThickH_->Fill(siStripClusterInfo_.charge(), wfac);
1952  hOnTrkClusWidthThickH_->Fill(siStripClusterInfo_.width(), wfac);
1953  } else {
1954  hOnTrkClusChargeThinH_->Fill(siStripClusterInfo_.charge(), wfac);
1955  hOnTrkClusWidthThinH_->Fill(siStripClusterInfo_.width(), wfac);
1956  }
1957  addClusterToMap(detid, &clusterS);
1958  } else {
1959  auto& cluster = clus.stripCluster();
1960  siStripClusterInfo_.setCluster(cluster, detid);
1961  if (thickness > 0.035) {
1962  hOnTrkClusChargeThickH_->Fill(siStripClusterInfo_.charge(), wfac);
1963  hOnTrkClusWidthThickH_->Fill(siStripClusterInfo_.width(), wfac);
1964  } else {
1965  hOnTrkClusChargeThinH_->Fill(siStripClusterInfo_.charge(), wfac);
1966  hOnTrkClusWidthThinH_->Fill(siStripClusterInfo_.width(), wfac);
1967  }
1968 
1969  addClusterToMap(detid, &cluster);
1970  }
1971 }
1972 void StandaloneTrackMonitor::addClusterToMap(uint32_t detid, const SiStripCluster* cluster) {
1973  std::map<uint32_t, std::set<const SiStripCluster*> >::iterator it = clusterMap_.find(detid);
1974  if (it == clusterMap_.end()) {
1975  std::set<const SiStripCluster*> s;
1976  s.insert(cluster);
1977  clusterMap_.insert(std::pair<uint32_t, std::set<const SiStripCluster*> >(detid, s));
1978  } else {
1979  std::set<const SiStripCluster*>& s = it->second;
1980  s.insert(cluster);
1981  }
1982 }
1983 
1984 const reco::Vertex* StandaloneTrackMonitor::findClosestVertex(
1985  const TransientVertex aTransVtx,
1986  const reco::VertexCollection* vertices) const { //same as in /DQMOffline/Alignment/src/DiMuonVertexSelector.cc
1987  reco::Vertex* defaultVtx = nullptr;
1988 
1989  if (!aTransVtx.isValid())
1990  return defaultVtx;
1991 
1992  // find the closest vertex to the secondary vertex in 3D
1993  VertexDistance3D vertTool3D;
1994  float minD = 9999.;
1995  int closestVtxIndex = 0;
1996  int counter = 0;
1997  for (const auto& vtx : *vertices) {
1998  double dist3D = vertTool3D.distance(aTransVtx, vtx).value();
1999  if (dist3D < minD) {
2000  minD = dist3D;
2001  closestVtxIndex = counter;
2002  }
2003  counter++;
2004  }
2005 
2006  if ((*vertices).at(closestVtxIndex).isValid()) {
2007  return &(vertices->at(closestVtxIndex));
2008  } else {
2009  return defaultVtx;
2010  }
2011 }
2012 
2013 // Define this as a plug-in
2014 #include "FWCore/Framework/interface/MakerMacros.h"
2015 DEFINE_FWK_MODULE(StandaloneTrackMonitor);
StandaloneTrackMonitor::nValidHitsPixEVsPhiH_
MonitorElement * nValidHitsPixEVsPhiH_
Definition: StandaloneTrackMonitor.cc:286
StandaloneTrackMonitor::pixelLayerwithMeasurementH_
MonitorElement * pixelLayerwithMeasurementH_
Definition: StandaloneTrackMonitor.cc:210
reco::HitPattern::numberOfLostTrackerHits
int numberOfLostTrackerHits(HitCategory category) const
Definition: HitPattern.h:893
StandaloneTrackMonitor::tkGeom_
const TrackerGeometry * tkGeom_
Definition: StandaloneTrackMonitor.cc:118
StandaloneTrackMonitor::nLostHitByLayerPixH_
MonitorElement * nLostHitByLayerPixH_
Definition: StandaloneTrackMonitor.cc:249
StandaloneTrackMonitor::nLostHitsPixBVsIterationH_
MonitorElement * nLostHitsPixBVsIterationH_
Definition: StandaloneTrackMonitor.cc:303
rpcPointValidation_cfi.recHit
recHit
Definition: rpcPointValidation_cfi.py:7
reco::HitPattern::numberOfValidPixelHits
int numberOfValidPixelHits() const
Definition: HitPattern.h:831
StandaloneTrackMonitor::hOnTrkClusChargeThickH_
MonitorElement * hOnTrkClusChargeThickH_
Definition: StandaloneTrackMonitor.cc:338
StandaloneTrackMonitor::trackTag_
const edm::InputTag trackTag_
Definition: StandaloneTrackMonitor.cc:83
Calorimetry_cff.thickness
thickness
Definition: Calorimetry_cff.py:115
StripSubdetector::TEC
static constexpr auto TEC
Definition: StripSubdetector.h:19
StandaloneTrackMonitor::nValidHitsPixBVsPhiH_
MonitorElement * nValidHitsPixBVsPhiH_
Definition: StandaloneTrackMonitor.cc:285
reco::HitPattern::stripTECHitFilter
static bool stripTECHitFilter(uint16_t pattern)
Definition: HitPattern.h:632
reco::HitPattern::getLayer
static uint32_t getLayer(uint16_t pattern)
Definition: HitPattern.h:721
StandaloneTrackMonitor::nLostHitsPixEVsIterationH_
MonitorElement * nLostHitsPixEVsIterationH_
Definition: StandaloneTrackMonitor.cc:304
StandaloneTrackMonitor::trackDeltaRwrtClosestTrack_
MonitorElement * trackDeltaRwrtClosestTrack_
Definition: StandaloneTrackMonitor.cc:150
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
deDxTools::esConsumes
ESGetTokenH3DDVariant esConsumes(std::string const &Record, edm::ConsumesCollector &)
Definition: DeDxTools.cc:283
StandaloneTrackMonitor::analyze
void analyze(edm::Event const &iEvent, edm::EventSetup const &iSetup) override
Definition: StandaloneTrackMonitor.cc:1305
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
reco::HitPattern::numberOfLostHits
int numberOfLostHits(HitCategory category) const
Definition: HitPattern.h:891
StandaloneTrackMonitor::hOnTrkClusChargeThinH_
MonitorElement * hOnTrkClusChargeThinH_
Definition: StandaloneTrackMonitor.cc:336
StandaloneTrackMonitor::stripTOBLayerwithMeasurementH_
MonitorElement * stripTOBLayerwithMeasurementH_
Definition: StandaloneTrackMonitor.cc:215
mps_fire.i
i
Definition: mps_fire.py:429
SiStripClusterInfo::setCluster
void setCluster(const SiStripCluster &cluster, int detId)
Definition: SiStripClusterInfo.cc:16
edm::EventSetup::getData
T const & getData(const ESGetToken< T, R > &iToken) const noexcept(false)
Definition: EventSetup.h:119
TrackerGeometry::idToDetUnit
const TrackerGeomDet * idToDetUnit(DetId) const override
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: TrackerGeometry.cc:183
StandaloneTrackMonitor::clusterToken_
const edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > > clusterToken_
Definition: StandaloneTrackMonitor.cc:93
TransientVertex::position
GlobalPoint position() const
Definition: TransientVertex.h:169
StandaloneTrackMonitor::geomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > geomToken_
Definition: StandaloneTrackMonitor.cc:113
muonME0PseudoDigis_cfi.etaError
etaError
Definition: muonME0PseudoDigis_cfi.py:10
StandaloneTrackMonitor::trackStoppingSourceH_
MonitorElement * trackStoppingSourceH_
Definition: StandaloneTrackMonitor.cc:143
StandaloneTrackMonitor::trackChi2oNDFVsEtaH_
MonitorElement * trackChi2oNDFVsEtaH_
Definition: StandaloneTrackMonitor.cc:320
run3scouting_cff.dxyError
dxyError
Definition: run3scouting_cff.py:145
StandaloneTrackMonitor::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: StandaloneTrackMonitor.cc:371
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:36
StandaloneTrackMonitor::StandaloneTrackMonitor
StandaloneTrackMonitor(const edm::ParameterSet &)
Definition: StandaloneTrackMonitor.cc:427
StandaloneTrackMonitor::trackChi2probVsPhiH_
MonitorElement * trackChi2probVsPhiH_
Definition: StandaloneTrackMonitor.cc:323
nano_mu_local_reco_cff.chi2
chi2
Definition: nano_mu_local_reco_cff.py:193
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:150
StandaloneTrackMonitor::hOffTrkClusWidthThickH_
MonitorElement * hOffTrkClusWidthThickH_
Definition: StandaloneTrackMonitor.cc:344
StandaloneTrackMonitor::trackToken_
const edm::EDGetTokenT< reco::TrackCollection > trackToken_
Definition: StandaloneTrackMonitor.cc:89
reco::Vertex::position
const Point & position() const
position
Definition: Vertex.h:128
StandaloneTrackMonitor::nMissingOuterHitEH_
MonitorElement * nMissingOuterHitEH_
Definition: StandaloneTrackMonitor.cc:207
reco::HitPattern::trackerLayersWithMeasurement
int trackerLayersWithMeasurement() const
Definition: HitPattern.cc:535
StandaloneTrackMonitor::nLostHitsTIBVsIterationH_
MonitorElement * nLostHitsTIBVsIterationH_
Definition: StandaloneTrackMonitor.cc:306
reco::HitPattern::numberOfAllHits
int numberOfAllHits(HitCategory category) const
Definition: HitPattern.h:804
StandaloneTrackMonitor::puSummaryToken_
const edm::EDGetTokenT< std::vector< PileupSummaryInfo > > puSummaryToken_
Definition: StandaloneTrackMonitor.cc:92
IPTools::signedTransverseImpactParameter
std::pair< bool, Measurement1D > signedTransverseImpactParameter(const reco::TransientTrack &track, const GlobalVector &direction, const reco::Vertex &vertex)
Definition: IPTools.cc:57
edm::Handle::product
T const * product() const
Definition: Handle.h:70
reco::HitPattern::stripLayersWithMeasurement
int stripLayersWithMeasurement() const
Definition: HitPattern.h:1005
L1BJetProducer_cff.vtx
vtx
Definition: L1BJetProducer_cff.py:10
reco::TrackCollection
std::vector< Track > TrackCollection
collection of Tracks
Definition: TrackFwd.h:14
vertices_cff.ntracks
ntracks
Definition: vertices_cff.py:33
run3scouting_cff.qoverp
qoverp
Definition: run3scouting_cff.py:143
StandaloneTrackMonitor::stripLayerwithMeasurementH_
MonitorElement * stripLayerwithMeasurementH_
Definition: StandaloneTrackMonitor.cc:213
reco::HitPattern::numberOfValidStripTIBHits
int numberOfValidStripTIBHits() const
Definition: HitPattern.h:847
StandaloneTrackMonitor::sip3dToBS2validpixelhitsH_
MonitorElement * sip3dToBS2validpixelhitsH_
Definition: StandaloneTrackMonitor.cc:179
pwdgSkimBPark_cfi.beamSpot
beamSpot
Definition: pwdgSkimBPark_cfi.py:5
IPTools::signedImpactParameter3D
std::pair< bool, Measurement1D > signedImpactParameter3D(const reco::TransientTrack &track, const GlobalVector &direction, const reco::Vertex &vertex)
Definition: IPTools.cc:81
reco::HitPattern::stripTIDLayersWithMeasurement
int stripTIDLayersWithMeasurement() const
Definition: HitPattern.cc:609
reco::HitPattern::pixelEndcapLayersWithMeasurement
int pixelEndcapLayersWithMeasurement() const
Definition: HitPattern.cc:588
reco::HitPattern::numberOfValidStripTECHits
int numberOfValidStripTECHits() const
Definition: HitPattern.h:859
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
StandaloneTrackMonitor::nValidHitsTOBVsPhiH_
MonitorElement * nValidHitsTOBVsPhiH_
Definition: StandaloneTrackMonitor.cc:289
reco::HitPattern::numberOfLostPixelBarrelHits
int numberOfLostPixelBarrelHits(HitCategory category) const
Definition: HitPattern.h:917
l1tTrackerHTMiss_cfi.deltaZ
deltaZ
Definition: l1tTrackerHTMiss_cfi.py:13
StandaloneTrackMonitor::yPointOfClosestApproachVsZ0wrtPVH_
MonitorElement * yPointOfClosestApproachVsZ0wrtPVH_
Definition: StandaloneTrackMonitor.cc:149
edmNew::DetSet::const_iterator
data_type const * const_iterator
Definition: DetSetNew.h:31
nano_mu_digi_cff.bx
bx
Definition: nano_mu_digi_cff.py:41
DeadROC_duringRun.f1
f1
Definition: DeadROC_duringRun.py:219
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
l1t::VertexCollection
std::vector< Vertex > VertexCollection
Definition: Vertex.h:31
electrons_cff.ip3d
ip3d
Definition: electrons_cff.py:319
StandaloneTrackMonitor::hOnTrkClusWidthThinH_
MonitorElement * hOnTrkClusWidthThinH_
Definition: StandaloneTrackMonitor.cc:337
reco::HitPattern::missingHitFilter
static bool missingHitFilter(uint16_t pattern)
Definition: HitPattern.h:798
reco::HitPattern::numberOfAllTrackerHits
int numberOfAllTrackerHits(HitCategory category) const
Definition: HitPattern.h:809
StandaloneTrackMonitor::iperr3dToBS2validpixelhitsH_
MonitorElement * iperr3dToBS2validpixelhitsH_
Definition: StandaloneTrackMonitor.cc:177
StandaloneTrackMonitor::trackChi2probVsEtaH_
MonitorElement * trackChi2probVsEtaH_
Definition: StandaloneTrackMonitor.cc:322
HLT_2024v14_cff.track
track
Definition: HLT_2024v14_cff.py:9565
StandaloneTrackMonitor::Jet_chargedMultiplicity_
MonitorElement * Jet_chargedMultiplicity_
Definition: StandaloneTrackMonitor.cc:357
run3scouting_cff.qoverpError
qoverpError
Definition: run3scouting_cff.py:147
reco::HitPattern::pixelBarrelHitFilter
static bool pixelBarrelHitFilter(uint16_t pattern)
Definition: HitPattern.h:590
StandaloneTrackMonitor::transTrackToken_
const edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord > transTrackToken_
Definition: StandaloneTrackMonitor.cc:114
StandaloneTrackMonitor::hOffTrkClusWidthThinH_
MonitorElement * hOffTrkClusWidthThinH_
Definition: StandaloneTrackMonitor.cc:342
StandaloneTrackMonitor::puScaleFactorFile_
const std::string puScaleFactorFile_
Definition: StandaloneTrackMonitor.cc:101
StandaloneTrackMonitor::hOffTrkClusChargeThinH_
MonitorElement * hOffTrkClusChargeThinH_
Definition: StandaloneTrackMonitor.cc:341
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
reco::HitPattern::stripTECLayersWithMeasurement
int stripTECLayersWithMeasurement() const
Definition: HitPattern.cc:629
StandaloneTrackMonitor::siStripClusterInfo_
SiStripClusterInfo siStripClusterInfo_
Definition: StandaloneTrackMonitor.cc:81
DiDispStaMuonMonitor_cfi.pt
pt
Definition: DiDispStaMuonMonitor_cfi.py:39
DeadROC_duringRun.dir
dir
Definition: DeadROC_duringRun.py:23
StandaloneTrackMonitor::nLostHitsTECVsIterationH_
MonitorElement * nLostHitsTECVsIterationH_
Definition: StandaloneTrackMonitor.cc:308
StandaloneTrackMonitor::mvaTrackTag_
const edm::InputTag mvaTrackTag_
Definition: StandaloneTrackMonitor.cc:104
reco::HitPattern::pixelEndcapHitFilter
static bool pixelEndcapHitFilter(uint16_t pattern)
Definition: HitPattern.h:599
reco::HitPattern::numberOfLostStripHits
int numberOfLostStripHits(HitCategory category) const
Definition: HitPattern.h:925
reco::HitPattern::numberOfValidStripHits
int numberOfValidStripHits() const
Definition: HitPattern.h:843
StandaloneTrackMonitor::jetsTag_
const edm::InputTag jetsTag_
Definition: StandaloneTrackMonitor.cc:88
StandaloneTrackMonitor::trackIperr3dVsEtaH_
MonitorElement * trackIperr3dVsEtaH_
Definition: StandaloneTrackMonitor.cc:325
StandaloneTrackMonitor::nMissingExpectedInnerHitsH_
MonitorElement * nMissingExpectedInnerHitsH_
Definition: StandaloneTrackMonitor.cc:220
KalmanVertexFitter::vertex
CachingVertex< 5 > vertex(const std::vector< reco::TransientTrack > &tracks) const override
Definition: KalmanVertexFitter.h:49
run3scouting_cff.nValidPixelHits
nValidPixelHits
Definition: run3scouting_cff.py:88
StandaloneTrackMonitor::nLostHitsTOBVsIterationH_
MonitorElement * nLostHitsTOBVsIterationH_
Definition: StandaloneTrackMonitor.cc:307
StandaloneTrackMonitor::vertexTag_
const edm::InputTag vertexTag_
Definition: StandaloneTrackMonitor.cc:85
TransientTrackBuilder::build
reco::TransientTrack build(const reco::Track *p) const
Definition: TransientTrackBuilder.cc:16
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:292
AlignmentTracksFromVertexSelector_cfi.vertices
vertices
Definition: AlignmentTracksFromVertexSelector_cfi.py:5
Bounds::thickness
virtual float thickness() const =0
DiMuonV_cfg.tracks
tracks
Definition: DiMuonV_cfg.py:215
StandaloneTrackMonitor::nLostHitsTIDVsEtaH_
MonitorElement * nLostHitsTIDVsEtaH_
Definition: StandaloneTrackMonitor.cc:300
reco::HitPattern::numberOfLostStripTECHits
int numberOfLostStripTECHits(HitCategory category) const
Definition: HitPattern.h:941
quality
string quality
Definition: beamSpotDipStandalone.cc:53
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
StandaloneTrackMonitor::TrackPtHistoPar_
const edm::ParameterSet TrackPtHistoPar_
Definition: StandaloneTrackMonitor.cc:117
dqm::implementation::IBooker::book1DD
MonitorElement * book1DD(TString const &name, TString const &title, int nchX, double lowX, double highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:155
StandaloneTrackMonitor::nLostHitsPixEVsEtaH_
MonitorElement * nLostHitsPixEVsEtaH_
Definition: StandaloneTrackMonitor.cc:295
StandaloneTrackMonitor::nMissingInnerHitEH_
MonitorElement * nMissingInnerHitEH_
Definition: StandaloneTrackMonitor.cc:205
iEvent
int iEvent
Definition: GenABIO.cc:224
StandaloneTrackMonitor::nLostHitsTIBVsEtaH_
MonitorElement * nLostHitsTIBVsEtaH_
Definition: StandaloneTrackMonitor.cc:297
run3scouting_cff.phiError
phiError
Definition: run3scouting_cff.py:149
StandaloneTrackMonitor::nLostHitsTIDVsIterationH_
MonitorElement * nLostHitsTIDVsIterationH_
Definition: StandaloneTrackMonitor.cc:309
MVACollection
std::vector< float > MVACollection
Definition: StandaloneTrackMonitor.cc:51
StandaloneTrackMonitor::nValidHitsPixVsEtaH_
MonitorElement * nValidHitsPixVsEtaH_
Definition: StandaloneTrackMonitor.cc:275
StandaloneTrackMonitor::hOnTrkClusWidthThickH_
MonitorElement * hOnTrkClusWidthThickH_
Definition: StandaloneTrackMonitor.cc:339
StandaloneTrackMonitor::stripTIBLayerwithMeasurementH_
MonitorElement * stripTIBLayerwithMeasurementH_
Definition: StandaloneTrackMonitor.cc:214
StandaloneTrackMonitor::nValidHitsTIDVsEtaH_
MonitorElement * nValidHitsTIDVsEtaH_
Definition: StandaloneTrackMonitor.cc:282
StandaloneTrackMonitor::nValidHitsStripVsPhiH_
MonitorElement * nValidHitsStripVsPhiH_
Definition: StandaloneTrackMonitor.cc:287
SiStripMatchedRecHit2D::monoCluster
SiStripCluster const & monoCluster() const
Definition: SiStripMatchedRecHit2D.h:41
StandaloneTrackMonitor::nValidHitsTECVsEtaH_
MonitorElement * nValidHitsTECVsEtaH_
Definition: StandaloneTrackMonitor.cc:281
TransientVertex::isValid
bool isValid() const
Definition: TransientVertex.h:195
dqm::implementation::IBooker::bookProfile
MonitorElement * bookProfile(TString const &name, TString const &title, int nchX, double lowX, double highX, int, double lowY, double highY, char const *option="s", FUNC onbooking=NOOP())
Definition: DQMStore.h:408
StandaloneTrackMonitor::trackMVAs
std::vector< MonitorElement * > trackMVAs
Definition: StandaloneTrackMonitor.cc:346
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:23
StandaloneTrackMonitor::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: StandaloneTrackMonitor.cc:631
reco::HitPattern::numberOfValidPixelEndcapHits
int numberOfValidPixelEndcapHits() const
Definition: HitPattern.h:839
StandaloneTrackMonitor::trackSip2dVsChi2prob2DH_
MonitorElement * trackSip2dVsChi2prob2DH_
Definition: StandaloneTrackMonitor.cc:334
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
StopReasonName::StopReasonName
const std::string StopReasonName[]
Definition: TrajectoryStopReasons.cc:3
reco::HitPattern::numberOfLostPixelEndcapHits
int numberOfLostPixelEndcapHits(HitCategory category) const
Definition: HitPattern.h:921
StandaloneTrackMonitor::puSummaryTag_
const edm::InputTag puSummaryTag_
Definition: StandaloneTrackMonitor.cc:86
StandaloneTrackMonitor::nValidHitsPixVsPhiH_
MonitorElement * nValidHitsPixVsPhiH_
Definition: StandaloneTrackMonitor.cc:284
StandaloneTrackMonitor::DistanceOfClosestApproachToPVH_
MonitorElement * DistanceOfClosestApproachToPVH_
Definition: StandaloneTrackMonitor.cc:145
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
reco::HitPattern::stripTOBLayersWithMeasurement
int stripTOBLayersWithMeasurement() const
Definition: HitPattern.cc:619
edm::EventSetup::get
T get() const
Definition: EventSetup.h:79
edm::Transition
Transition
Definition: Transition.h:12
StandaloneTrackMonitor::nMissingOuterHitBH_
MonitorElement * nMissingOuterHitBH_
Definition: StandaloneTrackMonitor.cc:206
reco::HitPattern::numberOfValidStripTIDHits
int numberOfValidStripTIDHits() const
Definition: HitPattern.h:851
StandaloneTrackMonitor::TrackPhiHistoPar_
const edm::ParameterSet TrackPhiHistoPar_
Definition: StandaloneTrackMonitor.cc:116
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
StandaloneTrackMonitor::nValidHitsVsCosThetaH_
MonitorElement * nValidHitsVsCosThetaH_
Definition: StandaloneTrackMonitor.cc:272
run3scouting_cff.nValidStripHits
nValidStripHits
Definition: run3scouting_cff.py:89
StandaloneTrackMonitor::nLostHitsTOBVsEtaH_
MonitorElement * nLostHitsTOBVsEtaH_
Definition: StandaloneTrackMonitor.cc:298
StandaloneTrackMonitor::pixelELayerwithMeasurementH_
MonitorElement * pixelELayerwithMeasurementH_
Definition: StandaloneTrackMonitor.cc:212
StripSubdetector::TOB
static constexpr auto TOB
Definition: StripSubdetector.h:18
StandaloneTrackMonitor::ip3dToPV2validpixelhitsH_
MonitorElement * ip3dToPV2validpixelhitsH_
Definition: StandaloneTrackMonitor.cc:174
reco::HitPattern::numberOfLostStripTIDHits
int numberOfLostStripTIDHits(HitCategory category) const
Definition: HitPattern.h:933
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
dqm::impl::MonitorElement::setBinLabel
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:772
StandaloneTrackMonitor::nLostHitsTIDVsPhiH_
MonitorElement * nLostHitsTIDVsPhiH_
Definition: StandaloneTrackMonitor.cc:318
StandaloneTrackMonitor::trackMVAsVsPtProfile
std::vector< MonitorElement * > trackMVAsVsPtProfile
Definition: StandaloneTrackMonitor.cc:348
run3scouting_cff.dzError
dzError
Definition: run3scouting_cff.py:146
reco::HitPattern::getHitPattern
uint16_t getHitPattern(HitCategory category, int position) const
Definition: HitPattern.h:537
hcalRecHitTable_cff.detId
detId
Definition: hcalRecHitTable_cff.py:12
StandaloneTrackMonitor::trackMVAsHP
std::vector< MonitorElement * > trackMVAsHP
Definition: StandaloneTrackMonitor.cc:347
btvMC_cff.vy
vy
Definition: btvMC_cff.py:26
StandaloneTrackMonitor::trackIperr3dVsChi2prob2DH_
MonitorElement * trackIperr3dVsChi2prob2DH_
Definition: StandaloneTrackMonitor.cc:331
StandaloneTrackMonitor::stripTECLayerwithMeasurementH_
MonitorElement * stripTECLayerwithMeasurementH_
Definition: StandaloneTrackMonitor.cc:217
StandaloneTrackMonitor::trackIperr3dVsChi2probH_
MonitorElement * trackIperr3dVsChi2probH_
Definition: StandaloneTrackMonitor.cc:326
QualityMaskCollection
std::vector< unsigned char > QualityMaskCollection
Definition: StandaloneTrackMonitor.cc:52
reco::HitPattern::numberOfLostStripTOBHits
int numberOfLostStripTOBHits(HitCategory category) const
Definition: HitPattern.h:937
StandaloneTrackMonitor::nLostHitsVsCosThetaH_
MonitorElement * nLostHitsVsCosThetaH_
Definition: StandaloneTrackMonitor.cc:253
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:50
StandaloneTrackMonitor::nValidHitsTIBVsEtaH_
MonitorElement * nValidHitsTIBVsEtaH_
Definition: StandaloneTrackMonitor.cc:279
reco::deltaR
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:30
VertexDistance3D::distance
Measurement1D distance(const GlobalPoint &vtx1Position, const GlobalError &vtx1PositionError, const GlobalPoint &vtx2Position, const GlobalError &vtx2PositionError) const override
Definition: VertexDistance3D.cc:17
StandaloneTrackMonitor::hOffTrkClusChargeThickH_
MonitorElement * hOffTrkClusChargeThickH_
Definition: StandaloneTrackMonitor.cc:343
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
StandaloneTrackMonitor::tcProducer_
const edm::InputTag tcProducer_
Definition: StandaloneTrackMonitor.cc:106
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:193
nano_mu_digi_cff.layer
layer
Definition: nano_mu_digi_cff.py:28
StandaloneTrackMonitor::nMissingExpectedOuterHitsH_
MonitorElement * nMissingExpectedOuterHitsH_
Definition: StandaloneTrackMonitor.cc:221
StandaloneTrackMonitor::clusterTag_
const edm::InputTag clusterTag_
Definition: StandaloneTrackMonitor.cc:87
StandaloneTrackMonitor::nvalidTrackerHitsH_
MonitorElement * nvalidTrackerHitsH_
Definition: StandaloneTrackMonitor.cc:184
StandaloneTrackMonitor::nValidHitsPixBVsEtaH_
MonitorElement * nValidHitsPixBVsEtaH_
Definition: StandaloneTrackMonitor.cc:276
StandaloneTrackMonitor::addClusterToMap
void addClusterToMap(uint32_t detid, const SiStripCluster *cluster)
Definition: StandaloneTrackMonitor.cc:1972
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
StandaloneTrackMonitor::nLostHitsTECVsEtaH_
MonitorElement * nLostHitsTECVsEtaH_
Definition: StandaloneTrackMonitor.cc:299
StripSubdetector::TIB
static constexpr auto TIB
Definition: StripSubdetector.h:16
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
StandaloneTrackMonitor::trackChi2oNDFVsPhiH_
MonitorElement * trackChi2oNDFVsPhiH_
Definition: StandaloneTrackMonitor.cc:321
StandaloneTrackMonitor::nLostHitsStripVsEtaH_
MonitorElement * nLostHitsStripVsEtaH_
Definition: StandaloneTrackMonitor.cc:296
StandaloneTrackMonitor::bsToken_
const edm::EDGetTokenT< reco::BeamSpot > bsToken_
Definition: StandaloneTrackMonitor.cc:90
StandaloneTrackMonitor::mvaTrackToken_
const edm::EDGetTokenT< edm::View< reco::Track > > mvaTrackToken_
Definition: StandaloneTrackMonitor.cc:105
StandaloneTrackMonitor::trkLayerwithMeasurementH_
MonitorElement * trkLayerwithMeasurementH_
Definition: StandaloneTrackMonitor.cc:209
ApeEstimator_cff.width
width
Definition: ApeEstimator_cff.py:24
L1TEGammaOffline_cfi.nVertex
nVertex
Definition: L1TEGammaOffline_cfi.py:60
dqm::implementation::IBooker::book2D
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:221
StandaloneTrackMonitor::DistanceOfClosestApproachToPVZoomedH_
MonitorElement * DistanceOfClosestApproachToPVZoomedH_
Definition: StandaloneTrackMonitor.cc:146
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
nano_mu_local_reco_cff.bool
bool
Definition: nano_mu_local_reco_cff.py:14
StandaloneTrackMonitor::clusterMap_
std::map< uint32_t, std::set< const SiStripCluster * > > clusterMap_
Definition: StandaloneTrackMonitor.cc:368
SiStripClusterInfo::width
uint16_t width() const
Definition: SiStripClusterInfo.h:31
SiStripMatchedRecHit2D::stereoCluster
SiStripCluster const & stereoCluster() const
Definition: SiStripMatchedRecHit2D.h:40
StandaloneTrackMonitor::trackSip2dVsEta2DH_
MonitorElement * trackSip2dVsEta2DH_
Definition: StandaloneTrackMonitor.cc:333
reco::TrackBase::algoNames
static const std::string algoNames[]
Definition: TrackBase.h:147
reco::HitPattern::stripTOBHitFilter
static bool stripTOBHitFilter(uint16_t pattern)
Definition: HitPattern.h:628
StandaloneTrackMonitor::xPointOfClosestApproachVsZ0wrtPVH_
MonitorElement * xPointOfClosestApproachVsZ0wrtPVH_
Definition: StandaloneTrackMonitor.cc:148
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
Measurement1D::value
double value() const
Definition: Measurement1D.h:25
StandaloneTrackMonitor::nLostHitByLayerStripH_
MonitorElement * nLostHitByLayerStripH_
Definition: StandaloneTrackMonitor.cc:250
reco::HitPattern::numberOfLostPixelHits
int numberOfLostPixelHits(HitCategory category) const
Definition: HitPattern.h:913
reco::HitPattern::stripTIBHitFilter
static bool stripTIBHitFilter(uint16_t pattern)
Definition: HitPattern.h:620
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it
auto & it
Definition: splitVertices.h:28
StandaloneTrackMonitor::vtrack_
std::vector< float > vtrack_
Definition: StandaloneTrackMonitor.cc:367
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:46
StandaloneTrackMonitor::dqmBeginRun
void dqmBeginRun(const edm::Run &iRun, const edm::EventSetup &iSetup) override
Definition: StandaloneTrackMonitor.cc:627
StandaloneTrackMonitor::iperr3dToPV2validpixelhitsH_
MonitorElement * iperr3dToPV2validpixelhitsH_
Definition: StandaloneTrackMonitor.cc:176
btvMC_cff.vz
vz
Definition: btvMC_cff.py:27
edm
HLT enums.
Definition: AlignableModifier.h:19
StandaloneTrackMonitor::TrackEtaHistoPar_
const edm::ParameterSet TrackEtaHistoPar_
Definition: StandaloneTrackMonitor.cc:115
StandaloneTrackMonitor::vertexToken_
const edm::EDGetTokenT< reco::VertexCollection > vertexToken_
Definition: StandaloneTrackMonitor.cc:91
counter
static std::atomic< unsigned int > counter
Definition: SharedResourceNames.cc:17
StandaloneTrackMonitor::nLostHitsPixVsPhiH_
MonitorElement * nLostHitsPixVsPhiH_
Definition: StandaloneTrackMonitor.cc:311
StandaloneTrackMonitor::sip3dToPV2validpixelhitsH_
MonitorElement * sip3dToPV2validpixelhitsH_
Definition: StandaloneTrackMonitor.cc:178
StandaloneTrackMonitor::processClusters
void processClusters(edm::Event const &iEvent, edm::EventSetup const &iSetup, const TrackerGeometry &tkGeom, double wfac=1)
Definition: StandaloneTrackMonitor.cc:1870
StandaloneTrackMonitor::nLostHitsVsIterationH_
MonitorElement * nLostHitsVsIterationH_
Definition: StandaloneTrackMonitor.cc:255
StandaloneTrackMonitor::nValidHitsStripVsEtaH_
MonitorElement * nValidHitsStripVsEtaH_
Definition: StandaloneTrackMonitor.cc:278
StandaloneTrackMonitor::jetsToken_
const edm::EDGetTokenT< std::vector< reco::PFJet > > jetsToken_
Definition: StandaloneTrackMonitor.cc:94
StandaloneTrackMonitor::trackMVAsHPVsEtaProfile
std::vector< MonitorElement * > trackMVAsHPVsEtaProfile
Definition: StandaloneTrackMonitor.cc:351
reco::HitPattern::numberOfValidPixelBarrelHits
int numberOfValidPixelBarrelHits() const
Definition: HitPattern.h:835
StandaloneTrackMonitor::nValidHitsTOBVsEtaH_
MonitorElement * nValidHitsTOBVsEtaH_
Definition: StandaloneTrackMonitor.cc:280
StandaloneTrackMonitor::nLostHitsTIBVsPhiH_
MonitorElement * nLostHitsTIBVsPhiH_
Definition: StandaloneTrackMonitor.cc:315
btvMC_cff.vx
vx
Definition: btvMC_cff.py:25
reco::HitPattern::pixelBarrelLayersWithMeasurement
int pixelBarrelLayersWithMeasurement() const
Definition: HitPattern.cc:577
StandaloneTrackMonitor::trackIperr3dVsEta2DH_
MonitorElement * trackIperr3dVsEta2DH_
Definition: StandaloneTrackMonitor.cc:330
PixelVertexMonitor_cff.ndof
ndof
Definition: PixelVertexMonitor_cff.py:7
dqm::impl::MonitorElement::getNbinsX
virtual int getNbinsX() const
get # of bins in X-axis
Definition: MonitorElement.cc:575
StandaloneTrackMonitor::nLostHitsPixBVsEtaH_
MonitorElement * nLostHitsPixBVsEtaH_
Definition: StandaloneTrackMonitor.cc:294
reco::HitPattern::pixelLayersWithMeasurement
int pixelLayersWithMeasurement() const
Definition: HitPattern.cc:514
StandaloneTrackMonitor::nLostHitsStripVsIterationH_
MonitorElement * nLostHitsStripVsIterationH_
Definition: StandaloneTrackMonitor.cc:305
SiStripClusterInfo::charge
uint16_t charge() const
Definition: SiStripClusterInfo.h:43
ntrack
#define ntrack
StandaloneTrackMonitor::pixelBLayerwithMeasurementH_
MonitorElement * pixelBLayerwithMeasurementH_
Definition: StandaloneTrackMonitor.cc:211
StandaloneTrackMonitor::ip3dToBS2validpixelhitsH_
MonitorElement * ip3dToBS2validpixelhitsH_
Definition: StandaloneTrackMonitor.cc:175
StandaloneTrackMonitor::nLostHitsTOBVsPhiH_
MonitorElement * nLostHitsTOBVsPhiH_
Definition: StandaloneTrackMonitor.cc:316
StandaloneTrackMonitor::nLostHitsTECVsPhiH_
MonitorElement * nLostHitsTECVsPhiH_
Definition: StandaloneTrackMonitor.cc:317
reco::HitPattern::stripTIDHitFilter
static bool stripTIDHitFilter(uint16_t pattern)
Definition: HitPattern.h:624
StandaloneTrackMonitor::findClosestVertex
const reco::Vertex * findClosestVertex(const TransientVertex aTransVtx, const reco::VertexCollection *vertices) const
Definition: StandaloneTrackMonitor.cc:1984
StandaloneTrackMonitor::processHit
void processHit(const TrackingRecHit &recHit, edm::EventSetup const &iSetup, const TrackerGeometry &tkGeom, double wfac=1)
Definition: StandaloneTrackMonitor.cc:1914
StandaloneTrackMonitor::mvaQualityTokens_
std::vector< std::tuple< edm::EDGetTokenT< MVACollection >, edm::EDGetTokenT< QualityMaskCollection > > > mvaQualityTokens_
Definition: StandaloneTrackMonitor.cc:112
StandaloneTrackMonitor::trackMVAsVsEtaProfile
std::vector< MonitorElement * > trackMVAsVsEtaProfile
Definition: StandaloneTrackMonitor.cc:350
StandaloneTrackMonitor::nLostHitsPixBVsPhiH_
MonitorElement * nLostHitsPixBVsPhiH_
Definition: StandaloneTrackMonitor.cc:312
StandaloneTrackMonitor::nValidHitsTIBVsPhiH_
MonitorElement * nValidHitsTIBVsPhiH_
Definition: StandaloneTrackMonitor.cc:288
StandaloneTrackMonitor::trackScaleFactorFile_
const std::string trackScaleFactorFile_
Definition: StandaloneTrackMonitor.cc:102
StripSubdetector::TID
static constexpr auto TID
Definition: StripSubdetector.h:17
StandaloneTrackMonitor::nLostHitsPixVsIterationH_
MonitorElement * nLostHitsPixVsIterationH_
Definition: StandaloneTrackMonitor.cc:302
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:153
reco::HitPattern::numberOfValidTrackerHits
int numberOfValidTrackerHits() const
Definition: HitPattern.h:819
StandaloneTrackMonitor::stripTIDLayerwithMeasurementH_
MonitorElement * stripTIDLayerwithMeasurementH_
Definition: StandaloneTrackMonitor.cc:216
StandaloneTrackMonitor::nValidHitsPixEVsEtaH_
MonitorElement * nValidHitsPixEVsEtaH_
Definition: StandaloneTrackMonitor.cc:277
StandaloneTrackMonitor::nLostHitsPixVsEtaH_
MonitorElement * nLostHitsPixVsEtaH_
Definition: StandaloneTrackMonitor.cc:293
StandaloneTrackMonitor::nLostHitsPixEVsPhiH_
MonitorElement * nLostHitsPixEVsPhiH_
Definition: StandaloneTrackMonitor.cc:313
reco::HitPattern::numberOfLostStripTIBHits
int numberOfLostStripTIBHits(HitCategory category) const
Definition: HitPattern.h:929
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
StandaloneTrackMonitor::DistanceOfClosestApproachToPVVsPhiH_
MonitorElement * DistanceOfClosestApproachToPVVsPhiH_
Definition: StandaloneTrackMonitor.cc:147
reco::Vertex::isValid
bool isValid() const
Tells whether the vertex is valid.
Definition: Vertex.h:73
StandaloneTrackMonitor::nValidHitsTIDVsPhiH_
MonitorElement * nValidHitsTIDVsPhiH_
Definition: StandaloneTrackMonitor.cc:291
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
edm::Run
Definition: Run.h:45
StandaloneTrackMonitor::nValidHitsTECVsPhiH_
MonitorElement * nValidHitsTECVsPhiH_
Definition: StandaloneTrackMonitor.cc:290
StandaloneTrackMonitor::trackMVAsHPVsPtProfile
std::vector< MonitorElement * > trackMVAsHPVsPtProfile
Definition: StandaloneTrackMonitor.cc:349
StandaloneTrackMonitor::nLostHitsStripVsPhiH_
MonitorElement * nLostHitsStripVsPhiH_
Definition: StandaloneTrackMonitor.cc:314
StandaloneTrackMonitor::mvaProducers_
const std::vector< std::string > mvaProducers_
Definition: StandaloneTrackMonitor.cc:103
reco::HitPattern::numberOfValidStripTOBHits
int numberOfValidStripTOBHits() const
Definition: HitPattern.h:855
reco::HitPattern::stripTIBLayersWithMeasurement
int stripTIBLayersWithMeasurement() const
Definition: HitPattern.cc:599
StandaloneTrackMonitor::nMissingInnerHitBH_
MonitorElement * nMissingInnerHitBH_
Definition: StandaloneTrackMonitor.cc:204
Surface::bounds
const Bounds & bounds() const
Definition: Surface.h:87
Generated for CMSSW Reference Manual by doxygen 1.8.14