CMS 3D CMS Logo
TrackSplittingMonitor.cc
Go to the documentation of this file.
1 /*
2  * See header file for a description of this class.
3  *
4  * \author Suchandra Dutta , Giorgia Mila
5  */
6 
7 #include "DQM/TrackingMonitor/interface/TrackSplittingMonitor.h"
8 #include "DQMServices/Core/interface/DQMStore.h"
9 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
10 #include "DataFormats/TrackCandidate/interface/TrackCandidateCollection.h"
11 #include "DataFormats/TrackReco/interface/Track.h"
12 #include "DataFormats/TrackReco/interface/TrackFwd.h"
13 #include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"
14 #include "FWCore/MessageLogger/interface/MessageLogger.h"
15 #include "FWCore/Utilities/interface/InputTag.h"
16 #include "TrackingTools/PatternTools/interface/TSCBLBuilderNoMaterial.h"
17 #include "TrackingTools/PatternTools/interface/TSCPBuilderNoMaterial.h"
18 #include "TrackingTools/Records/interface/TransientRecHitRecord.h"
19 #include "TrackingTools/Records/interface/TransientTrackRecord.h"
20 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"
21 #include "TrackingTools/TransientTrack/interface/TransientTrack.h"
22 #include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
23 //#include "DQM/TrackingMonitor/interface/TrackAnalyzer.h"
24 #include <string>
25 
26 TrackSplittingMonitor::TrackSplittingMonitor(const edm::ParameterSet& iConfig)
27  : dqmStore_(edm::Service<DQMStore>().operator->()),
28  conf_(iConfig),
29  mfToken_(esConsumes()),
30  tkGeomToken_(esConsumes()),
31  dtGeomToken_(esConsumes()),
32  cscGeomToken_(esConsumes()),
33  rpcGeomToken_(esConsumes()) {
34  splitTracks_ = conf_.getParameter<edm::InputTag>("splitTrackCollection");
35  splitMuons_ = conf_.getParameter<edm::InputTag>("splitMuonCollection");
36  splitTracksToken_ = consumes<std::vector<reco::Track> >(splitTracks_);
37  splitMuonsToken_ = mayConsume<std::vector<reco::Muon> >(splitMuons_);
38 
39  plotMuons_ = conf_.getParameter<bool>("ifPlotMuons");
40 
41  // cuts
42  pixelHitsPerLeg_ = conf_.getParameter<int>("pixelHitsPerLeg");
43  totalHitsPerLeg_ = conf_.getParameter<int>("totalHitsPerLeg");
44  d0Cut_ = conf_.getParameter<double>("d0Cut");
45  dzCut_ = conf_.getParameter<double>("dzCut");
46  ptCut_ = conf_.getParameter<double>("ptCut");
47  norchiCut_ = conf_.getParameter<double>("norchiCut");
48 }
49 
50 TrackSplittingMonitor::~TrackSplittingMonitor() {}
51 
52 void TrackSplittingMonitor::bookHistograms(DQMStore::IBooker& ibooker,
53  edm::Run const& /* iRun */,
54  edm::EventSetup const& /* iSetup */)
55 
56 {
57  std::string MEFolderName = conf_.getParameter<std::string>("FolderName");
58  ibooker.setCurrentFolder(MEFolderName);
59 
60  // bin declarations
61  int ddxyBin = conf_.getParameter<int>("ddxyBin");
62  double ddxyMin = conf_.getParameter<double>("ddxyMin");
63  double ddxyMax = conf_.getParameter<double>("ddxyMax");
64 
65  int ddzBin = conf_.getParameter<int>("ddzBin");
66  double ddzMin = conf_.getParameter<double>("ddzMin");
67  double ddzMax = conf_.getParameter<double>("ddzMax");
68 
69  int dphiBin = conf_.getParameter<int>("dphiBin");
70  double dphiMin = conf_.getParameter<double>("dphiMin");
71  double dphiMax = conf_.getParameter<double>("dphiMax");
72 
73  int dthetaBin = conf_.getParameter<int>("dthetaBin");
74  double dthetaMin = conf_.getParameter<double>("dthetaMin");
75  double dthetaMax = conf_.getParameter<double>("dthetaMax");
76 
77  int dptBin = conf_.getParameter<int>("dptBin");
78  double dptMin = conf_.getParameter<double>("dptMin");
79  double dptMax = conf_.getParameter<double>("dptMax");
80 
81  int dcurvBin = conf_.getParameter<int>("dcurvBin");
82  double dcurvMin = conf_.getParameter<double>("dcurvMin");
83  double dcurvMax = conf_.getParameter<double>("dcurvMax");
84 
85  int normBin = conf_.getParameter<int>("normBin");
86  double normMin = conf_.getParameter<double>("normMin");
87  double normMax = conf_.getParameter<double>("normMax");
88 
89  // declare histogram
90  ddxyAbsoluteResiduals_tracker_ =
91  ibooker.book1D("ddxyAbsoluteResiduals_tracker", "ddxyAbsoluteResiduals_tracker", ddxyBin, ddxyMin, ddxyMax);
92  ddzAbsoluteResiduals_tracker_ =
93  ibooker.book1D("ddzAbsoluteResiduals_tracker", "ddzAbsoluteResiduals_tracker", ddzBin, ddzMin, ddzMax);
94  dphiAbsoluteResiduals_tracker_ =
95  ibooker.book1D("dphiAbsoluteResiduals_tracker", "dphiAbsoluteResiduals_tracker", dphiBin, dphiMin, dphiMax);
96  dthetaAbsoluteResiduals_tracker_ = ibooker.book1D(
97  "dthetaAbsoluteResiduals_tracker", "dthetaAbsoluteResiduals_tracker", dthetaBin, dthetaMin, dthetaMax);
98  dptAbsoluteResiduals_tracker_ =
99  ibooker.book1D("dptAbsoluteResiduals_tracker", "dptAbsoluteResiduals_tracker", dptBin, dptMin, dptMax);
100  dcurvAbsoluteResiduals_tracker_ =
101  ibooker.book1D("dcurvAbsoluteResiduals_tracker", "dcurvAbsoluteResiduals_tracker", dcurvBin, dcurvMin, dcurvMax);
102 
103  ddxyNormalizedResiduals_tracker_ =
104  ibooker.book1D("ddxyNormalizedResiduals_tracker", "ddxyNormalizedResiduals_tracker", normBin, normMin, normMax);
105  ddzNormalizedResiduals_tracker_ =
106  ibooker.book1D("ddzNormalizedResiduals_tracker", "ddzNormalizedResiduals_tracker", normBin, normMin, normMax);
107  dphiNormalizedResiduals_tracker_ =
108  ibooker.book1D("dphiNormalizedResiduals_tracker", "dphiNormalizedResiduals_tracker", normBin, normMin, normMax);
109  dthetaNormalizedResiduals_tracker_ = ibooker.book1D(
110  "dthetaNormalizedResiduals_tracker", "dthetaNormalizedResiduals_tracker", normBin, normMin, normMax);
111  dptNormalizedResiduals_tracker_ =
112  ibooker.book1D("dptNormalizedResiduals_tracker", "dptNormalizedResiduals_tracker", normBin, normMin, normMax);
113  dcurvNormalizedResiduals_tracker_ =
114  ibooker.book1D("dcurvNormalizedResiduals_tracker", "dcurvNormalizedResiduals_tracker", normBin, normMin, normMax);
115 
116  if (plotMuons_) {
117  ddxyAbsoluteResiduals_global_ =
118  ibooker.book1D("ddxyAbsoluteResiduals_global", "ddxyAbsoluteResiduals_global", ddxyBin, ddxyMin, ddxyMax);
119  ddzAbsoluteResiduals_global_ =
120  ibooker.book1D("ddzAbsoluteResiduals_global", "ddzAbsoluteResiduals_global", ddzBin, ddzMin, ddzMax);
121  dphiAbsoluteResiduals_global_ =
122  ibooker.book1D("dphiAbsoluteResiduals_global", "dphiAbsoluteResiduals_global", dphiBin, dphiMin, dphiMax);
123  dthetaAbsoluteResiduals_global_ = ibooker.book1D(
124  "dthetaAbsoluteResiduals_global", "dthetaAbsoluteResiduals_global", dthetaBin, dthetaMin, dthetaMax);
125  dptAbsoluteResiduals_global_ =
126  ibooker.book1D("dptAbsoluteResiduals_global", "dptAbsoluteResiduals_global", dptBin, dptMin, dptMax);
127  dcurvAbsoluteResiduals_global_ =
128  ibooker.book1D("dcurvAbsoluteResiduals_global", "dcurvAbsoluteResiduals_global", dcurvBin, dcurvMin, dcurvMax);
129 
130  ddxyNormalizedResiduals_global_ =
131  ibooker.book1D("ddxyNormalizedResiduals_global", "ddxyNormalizedResiduals_global", normBin, normMin, normMax);
132  ddzNormalizedResiduals_global_ =
133  ibooker.book1D("ddzNormalizedResiduals_global", "ddzNormalizedResiduals_global", normBin, normMin, normMax);
134  dphiNormalizedResiduals_global_ =
135  ibooker.book1D("dphiNormalizedResiduals_global", "dphiNormalizedResiduals_global", normBin, normMin, normMax);
136  dthetaNormalizedResiduals_global_ = ibooker.book1D(
137  "dthetaNormalizedResiduals_global", "dthetaNormalizedResiduals_global", normBin, normMin, normMax);
138  dptNormalizedResiduals_global_ =
139  ibooker.book1D("dptNormalizedResiduals_global", "dptNormalizedResiduals_global", normBin, normMin, normMax);
140  dcurvNormalizedResiduals_global_ =
141  ibooker.book1D("dcurvNormalizedResiduals_global", "dcurvNormalizedResiduals_global", normBin, normMin, normMax);
142  }
143 
144  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta d_{xy})/#sqrt{2} [#mum]");
145  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta d_{z})/#sqrt{2} [#mum]");
146  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta #phi)/#sqrt{2} [mrad]");
147  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta #theta)/#sqrt{2} [mrad]");
148  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta pT)/#sqrt{2} [GeV]");
149  ddxyAbsoluteResiduals_tracker_->setAxisTitle("(#delta (1/pT))/#sqrt{2} [GeV^{-1}]");
150 
151  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta d_{xy}/#sigma(d_{xy}");
152  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta d_{z}/#sigma(d_{z})");
153  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta #phi/#sigma(d_{#phi})");
154  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta #theta/#sigma(d_{#theta})");
155  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta p_{T}/#sigma(p_{T})");
156  ddxyNormalizedResiduals_tracker_->setAxisTitle("#delta 1/p_{T}/#sigma(1/p_{T})");
157 
158  if (plotMuons_) {
159  ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta d_{xy})/#sqrt{2} [#mum]");
160  ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta d_{z})/#sqrt{2} [#mum]");
161  ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta #phi)/#sqrt{2} [mrad]");
162  ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta #theta)/#sqrt{2} [mrad]");
163  ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta pT)/#sqrt{2} [GeV]");
164  ddxyAbsoluteResiduals_global_->setAxisTitle("(#delta (1/pT))/#sqrt{2} [GeV^{-1}]");
165 
166  ddxyNormalizedResiduals_global_->setAxisTitle("#delta d_{xy}/#sigma(d_{xy}");
167  ddxyNormalizedResiduals_global_->setAxisTitle("#delta d_{z}/#sigma(d_{z})");
168  ddxyNormalizedResiduals_global_->setAxisTitle("#delta #phi/#sigma(d_{#phi})");
169  ddxyNormalizedResiduals_global_->setAxisTitle("#delta #theta/#sigma(d_{#theta})");
170  ddxyNormalizedResiduals_global_->setAxisTitle("#delta p_{T}/#sigma(p_{T})");
171  ddxyNormalizedResiduals_global_->setAxisTitle("#delta 1/p_{T}/#sigma(1/p_{T})");
172  }
173 }
174 
175 //
176 // -- Analyse
177 //
178 void TrackSplittingMonitor::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
179  theMagField = &iSetup.getData(mfToken_);
180  theGeometry = &iSetup.getData(tkGeomToken_);
181  dtGeometry = &iSetup.getData(dtGeomToken_);
182  cscGeometry = &iSetup.getData(cscGeomToken_);
183  rpcGeometry = &iSetup.getData(rpcGeomToken_);
184 
185  edm::Handle<std::vector<reco::Track> > splitTracks = iEvent.getHandle(splitTracksToken_);
186  if (!splitTracks.isValid())
187  return;
188 
189  edm::Handle<std::vector<reco::Muon> > splitMuons;
190  if (plotMuons_) {
191  splitMuons = iEvent.getHandle(splitMuonsToken_);
192  }
193 
194  if (splitTracks->size() == 2) {
195  // check that there are 2 tracks in split track collection
196  edm::LogInfo("TrackSplittingMonitor") << "Split Track size: " << splitTracks->size();
197 
198  // split tracks calculations
199  reco::Track track1 = splitTracks->at(0);
200  reco::Track track2 = splitTracks->at(1);
201 
202  // -------------------------- basic selection ---------------------------
203 
204  // hit counting
205  // looping through the hits for track 1
206  double nRechits1 = 0;
207  double nRechitinBPIX1 = 0;
208  for (trackingRecHit_iterator iHit = track1.recHitsBegin(); iHit != track1.recHitsEnd(); ++iHit) {
209  if ((*iHit)->isValid()) {
210  nRechits1++;
211  int type = (*iHit)->geographicalId().subdetId();
212  if (type == int(PixelSubdetector::PixelBarrel)) {
213  ++nRechitinBPIX1;
214  }
215  }
216  }
217  // looping through the hits for track 2
218  double nRechits2 = 0;
219  double nRechitinBPIX2 = 0;
220  for (trackingRecHit_iterator iHit = track2.recHitsBegin(); iHit != track2.recHitsEnd(); ++iHit) {
221  if ((*iHit)->isValid()) {
222  nRechits2++;
223  int type = (*iHit)->geographicalId().subdetId();
224  if (type == int(PixelSubdetector::PixelBarrel)) {
225  ++nRechitinBPIX2;
226  }
227  }
228  }
229 
230  // DCA of each track
231  double d01 = track1.d0();
232  double dz1 = track1.dz();
233  double d02 = track2.d0();
234  double dz2 = track2.dz();
235 
236  // pT of each track
237  double pt1 = track1.pt();
238  double pt2 = track2.pt();
239 
240  // chi2 of each track
241  double norchi1 = track1.normalizedChi2();
242  double norchi2 = track2.normalizedChi2();
243 
244  // basic selection
245  // pixel hits and total hits
246  if ((nRechitinBPIX1 >= pixelHitsPerLeg_) && (nRechitinBPIX1 >= pixelHitsPerLeg_) &&
247  (nRechits1 >= totalHitsPerLeg_) && (nRechits2 >= totalHitsPerLeg_)) {
248  // dca cut
249  if (((fabs(d01) < d0Cut_)) && (fabs(d02) < d0Cut_) && (fabs(dz2) < dzCut_) && (fabs(dz2) < dzCut_)) {
250  // pt cut
251  if ((pt1 + pt2) / 2 < ptCut_) {
252  // chi2 cut
253  if ((norchi1 < norchiCut_) && (norchi2 < norchiCut_)) {
254  // passed all cuts...
255  edm::LogInfo("TrackSplittingMonitor") << " Setected after all cuts ?";
256 
257  double ddxyVal = d01 - d02;
258  double ddzVal = dz1 - dz2;
259  double dphiVal = track1.phi() - track2.phi();
260  double dthetaVal = track1.theta() - track2.theta();
261  double dptVal = pt1 - pt2;
262  double dcurvVal = (1 / pt1) - (1 / pt2);
263 
264  double d01ErrVal = track1.d0Error();
265  double d02ErrVal = track2.d0Error();
266  double dz1ErrVal = track1.dzError();
267  double dz2ErrVal = track2.dzError();
268  double phi1ErrVal = track1.phiError();
269  double phi2ErrVal = track2.phiError();
270  double theta1ErrVal = track1.thetaError();
271  double theta2ErrVal = track2.thetaError();
272  double pt1ErrVal = track1.ptError();
273  double pt2ErrVal = track2.ptError();
274 
275  ddxyAbsoluteResiduals_tracker_->Fill(10000.0 * ddxyVal / sqrt(2.0));
276  ddxyAbsoluteResiduals_tracker_->Fill(10000.0 * ddzVal / sqrt(2.0));
277  ddxyAbsoluteResiduals_tracker_->Fill(1000.0 * dphiVal / sqrt(2.0));
278  ddxyAbsoluteResiduals_tracker_->Fill(1000.0 * dthetaVal / sqrt(2.0));
279  ddxyAbsoluteResiduals_tracker_->Fill(dptVal / sqrt(2.0));
280  ddxyAbsoluteResiduals_tracker_->Fill(dcurvVal / sqrt(2.0));
281 
282  ddxyNormalizedResiduals_tracker_->Fill(ddxyVal / sqrt(d01ErrVal * d01ErrVal + d02ErrVal * d02ErrVal));
283  ddxyNormalizedResiduals_tracker_->Fill(ddzVal / sqrt(dz1ErrVal * dz1ErrVal + dz2ErrVal * dz2ErrVal));
284  ddxyNormalizedResiduals_tracker_->Fill(dphiVal / sqrt(phi1ErrVal * phi1ErrVal + phi2ErrVal * phi2ErrVal));
285  ddxyNormalizedResiduals_tracker_->Fill(dthetaVal /
286  sqrt(theta1ErrVal * theta1ErrVal + theta2ErrVal * theta2ErrVal));
287  ddxyNormalizedResiduals_tracker_->Fill(dptVal / sqrt(pt1ErrVal * pt1ErrVal + pt2ErrVal * pt2ErrVal));
288  ddxyNormalizedResiduals_tracker_->Fill(
289  dcurvVal / sqrt(pow(pt1ErrVal, 2) / pow(pt1, 4) + pow(pt2ErrVal, 2) / pow(pt2, 4)));
290 
291  // if do the same for split muons
292  if (plotMuons_ && splitMuons.isValid()) {
293  int gmCtr = 0;
294  bool topGlobalMuonFlag = false;
295  bool bottomGlobalMuonFlag = false;
296  int topGlobalMuon = -1;
297  int bottomGlobalMuon = -1;
298  double topGlobalMuonNorchi2 = 1e10;
299  double bottomGlobalMuonNorchi2 = 1e10;
300 
301  // check if usable split global muons
302  for (std::vector<reco::Muon>::const_iterator gmI = splitMuons->begin(); gmI != splitMuons->end(); gmI++) {
303  if (gmI->isTrackerMuon() && gmI->isStandAloneMuon() && gmI->isGlobalMuon()) {
304  reco::TrackRef trackerTrackRef1(splitTracks, 0);
305  reco::TrackRef trackerTrackRef2(splitTracks, 1);
306 
307  if (gmI->innerTrack() == trackerTrackRef1) {
308  if (gmI->globalTrack()->normalizedChi2() < topGlobalMuonNorchi2) {
309  topGlobalMuonFlag = true;
310  topGlobalMuonNorchi2 = gmI->globalTrack()->normalizedChi2();
311  topGlobalMuon = gmCtr;
312  }
313  }
314  if (gmI->innerTrack() == trackerTrackRef2) {
315  if (gmI->globalTrack()->normalizedChi2() < bottomGlobalMuonNorchi2) {
316  bottomGlobalMuonFlag = true;
317  bottomGlobalMuonNorchi2 = gmI->globalTrack()->normalizedChi2();
318  bottomGlobalMuon = gmCtr;
319  }
320  }
321  }
322  gmCtr++;
323  }
324 
325  if (bottomGlobalMuonFlag && topGlobalMuonFlag) {
326  reco::Muon muonTop = splitMuons->at(topGlobalMuon);
327  reco::Muon muonBottom = splitMuons->at(bottomGlobalMuon);
328 
329  reco::TrackRef glb1 = muonTop.globalTrack();
330  reco::TrackRef glb2 = muonBottom.globalTrack();
331 
332  double ddxyValGlb = glb1->d0() - glb2->d0();
333  double ddzValGlb = glb1->dz() - glb2->dz();
334  double dphiValGlb = glb1->phi() - glb2->phi();
335  double dthetaValGlb = glb1->theta() - glb2->theta();
336  double dptValGlb = glb1->pt() - glb2->pt();
337  double dcurvValGlb = (1 / glb1->pt()) - (1 / glb2->pt());
338 
339  double d01ErrValGlb = glb1->d0Error();
340  double d02ErrValGlb = glb2->d0Error();
341  double dz1ErrValGlb = glb1->dzError();
342  double dz2ErrValGlb = glb2->dzError();
343  double phi1ErrValGlb = glb1->phiError();
344  double phi2ErrValGlb = glb2->phiError();
345  double theta1ErrValGlb = glb1->thetaError();
346  double theta2ErrValGlb = glb2->thetaError();
347  double pt1ErrValGlb = glb1->ptError();
348  double pt2ErrValGlb = glb2->ptError();
349 
350  ddxyAbsoluteResiduals_global_->Fill(10000.0 * ddxyValGlb / sqrt(2.0));
351  ddxyAbsoluteResiduals_global_->Fill(10000.0 * ddzValGlb / sqrt(2.0));
352  ddxyAbsoluteResiduals_global_->Fill(1000.0 * dphiValGlb / sqrt(2.0));
353  ddxyAbsoluteResiduals_global_->Fill(1000.0 * dthetaValGlb / sqrt(2.0));
354  ddxyAbsoluteResiduals_global_->Fill(dptValGlb / sqrt(2.0));
355  ddxyAbsoluteResiduals_global_->Fill(dcurvValGlb / sqrt(2.0));
356 
357  ddxyNormalizedResiduals_global_->Fill(ddxyValGlb /
358  sqrt(d01ErrValGlb * d01ErrValGlb + d02ErrValGlb * d02ErrValGlb));
359  ddxyNormalizedResiduals_global_->Fill(ddzValGlb /
360  sqrt(dz1ErrValGlb * dz1ErrValGlb + dz2ErrValGlb * dz2ErrValGlb));
361  ddxyNormalizedResiduals_global_->Fill(
362  dphiValGlb / sqrt(phi1ErrValGlb * phi1ErrValGlb + phi2ErrValGlb * phi2ErrValGlb));
363  ddxyNormalizedResiduals_global_->Fill(
364  dthetaValGlb / sqrt(theta1ErrValGlb * theta1ErrValGlb + theta2ErrValGlb * theta2ErrValGlb));
365  ddxyNormalizedResiduals_global_->Fill(dptValGlb /
366  sqrt(pt1ErrValGlb * pt1ErrValGlb + pt2ErrValGlb * pt2ErrValGlb));
367  ddxyNormalizedResiduals_global_->Fill(
368  dcurvValGlb / sqrt(pow(pt1ErrValGlb, 2) / pow(pt1, 4) + pow(pt2ErrValGlb, 2) / pow(pt2, 4)));
369  }
370 
371  } // end of split muons loop
372  }
373  }
374  }
375  }
376  }
377 }
378 
379 DEFINE_FWK_MODULE(TrackSplittingMonitor);
TrackSplittingMonitor::cscGeomToken_
const edm::ESGetToken< CSCGeometry, MuonGeometryRecord > cscGeomToken_
Definition: TrackSplittingMonitor.h:65
TSCBLBuilderNoMaterial.h
TrackSplittingMonitor_cfi.dphiMax
dphiMax
Definition: TrackSplittingMonitor_cfi.py:34
TrackSplittingMonitor::dthetaNormalizedResiduals_global_
MonitorElement * dthetaNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:112
TrackSplittingMonitor::dcurvAbsoluteResiduals_global_
MonitorElement * dcurvAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:107
reco::TrackBase::ptError
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:754
MessageLogger.h
PixelSubdetector::PixelBarrel
Definition: PixelSubdetector.h:11
TrackSplittingMonitor::dphiAbsoluteResiduals_tracker_
MonitorElement * dphiAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:90
TrackCandidateCollection.h
TrackSplittingMonitor::ptCut_
double ptCut_
Definition: TrackSplittingMonitor.h:84
TrackSplittingMonitor::d0Cut_
double d0Cut_
Definition: TrackSplittingMonitor.h:82
reco::Track::recHitsBegin
trackingRecHit_iterator recHitsBegin() const
Iterator to first hit on the track.
Definition: Track.h:88
TrackSplittingMonitor_cfi.ddzMin
ddzMin
Definition: TrackSplittingMonitor_cfi.py:29
TrackSplittingMonitor_cfi.normMax
normMax
Definition: TrackSplittingMonitor_cfi.py:50
TrackSplittingMonitor_cfi.ddzMax
ddzMax
Definition: TrackSplittingMonitor_cfi.py:30
TrackSplittingMonitor::mfToken_
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > mfToken_
Definition: TrackSplittingMonitor.h:62
TrackSplittingMonitor_cfi.dcurvBin
dcurvBin
Definition: TrackSplittingMonitor_cfi.py:44
edm::Run
Definition: Run.h:45
TrackSplittingMonitor::dptAbsoluteResiduals_tracker_
MonitorElement * dptAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:92
TrackSplittingMonitor::TrackSplittingMonitor
TrackSplittingMonitor(const edm::ParameterSet &)
Definition: TrackSplittingMonitor.cc:26
TrackSplittingMonitor::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: TrackSplittingMonitor.cc:178
edm
HLT enums.
Definition: AlignableModifier.h:19
TransientRecHitRecord.h
TrackSplittingMonitor::dtGeomToken_
const edm::ESGetToken< DTGeometry, MuonGeometryRecord > dtGeomToken_
Definition: TrackSplittingMonitor.h:64
reco::TrackBase::thetaError
double thetaError() const
error on theta
Definition: TrackBase.h:757
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
TrackSplittingMonitor_cfi.dthetaMin
dthetaMin
Definition: TrackSplittingMonitor_cfi.py:37
TrackSplittingMonitor::dphiNormalizedResiduals_global_
MonitorElement * dphiNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:111
DQMStore.h
TrackSplittingMonitor::ddzNormalizedResiduals_tracker_
MonitorElement * ddzNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:96
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:125
TrackSplittingMonitor::dptAbsoluteResiduals_global_
MonitorElement * dptAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:106
TransientTrack.h
edm::Handle
Definition: AssociativeIterator.h:50
TrackSplittingMonitor_cfi.dcurvMin
dcurvMin
Definition: TrackSplittingMonitor_cfi.py:45
TrackSplittingMonitor_cfi.dcurvMax
dcurvMax
Definition: TrackSplittingMonitor_cfi.py:46
reco::Muon
Definition: Muon.h:27
edm::Ref< TrackCollection >
Utilities.operator
operator
Definition: Utilities.py:24
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:637
TrackSplittingMonitor::dphiNormalizedResiduals_tracker_
MonitorElement * dphiNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:97
TrackSplittingMonitor::theGeometry
const TrackerGeometry * theGeometry
Definition: TrackSplittingMonitor.h:68
TrackSplittingMonitor::dcurvAbsoluteResiduals_tracker_
MonitorElement * dcurvAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:93
TrackSplittingMonitor::conf_
edm::ParameterSet conf_
Definition: TrackSplittingMonitor.h:60
Track.h
TrackSplittingMonitor::ddxyAbsoluteResiduals_global_
MonitorElement * ddxyAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:102
TrackFwd.h
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
dqm::legacy::DQMStore
Definition: DQMStore.h:728
BeamSpot.h
TrackSplittingMonitor::splitTracksToken_
edm::EDGetTokenT< std::vector< reco::Track > > splitTracksToken_
Definition: TrackSplittingMonitor.h:76
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
TrackSplittingMonitor.h
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
reco::Track::recHitsEnd
trackingRecHit_iterator recHitsEnd() const
Iterator to last hit on the track.
Definition: Track.h:91
reco::Track
Definition: Track.h:27
TrackSplittingMonitor::dptNormalizedResiduals_tracker_
MonitorElement * dptNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:99
Service
HLT_FULL_cff.pt1
pt1
Definition: HLT_FULL_cff.py:9893
reco::TrackBase::dz
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:622
TrackSplittingMonitor::dthetaAbsoluteResiduals_global_
MonitorElement * dthetaAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:105
TrackSplittingMonitor_cfi.ddxyMin
ddxyMin
Definition: TrackSplittingMonitor_cfi.py:25
TrackSplittingMonitor::totalHitsPerLeg_
int totalHitsPerLeg_
Definition: TrackSplittingMonitor.h:81
TrackSplittingMonitor::theMagField
const MagneticField * theMagField
Definition: TrackSplittingMonitor.h:69
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:649
edm::OwnVector::const_iterator
Definition: OwnVector.h:41
TrackSplittingMonitor::ddzAbsoluteResiduals_global_
MonitorElement * ddzAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:103
TrackSplittingMonitor::~TrackSplittingMonitor
~TrackSplittingMonitor() override
Definition: TrackSplittingMonitor.cc:50
TrackSplittingMonitor::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: TrackSplittingMonitor.cc:52
TrackSplittingMonitor::ddxyAbsoluteResiduals_tracker_
MonitorElement * ddxyAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:88
TransientTrackBuilder.h
edm::ParameterSet
Definition: ParameterSet.h:47
TrackSplittingMonitor_cfi.dthetaMax
dthetaMax
Definition: TrackSplittingMonitor_cfi.py:38
TrackSplittingMonitor::cscGeometry
const CSCGeometry * cscGeometry
Definition: TrackSplittingMonitor.h:71
reco::Muon::globalTrack
virtual TrackRef globalTrack() const
reference to Track reconstructed in both tracked and muon detector
Definition: Muon.h:51
TrackSplittingMonitor::dphiAbsoluteResiduals_global_
MonitorElement * dphiAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:104
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
TrackSplittingMonitor::pixelHitsPerLeg_
int pixelHitsPerLeg_
Definition: TrackSplittingMonitor.h:80
TrackSplittingMonitor
Definition: TrackSplittingMonitor.h:41
TrackSplittingMonitor::splitMuonsToken_
edm::EDGetTokenT< std::vector< reco::Muon > > splitMuonsToken_
Definition: TrackSplittingMonitor.h:77
reco::TrackBase::dzError
double dzError() const
error on dz
Definition: TrackBase.h:778
TrackSplittingMonitor::dtGeometry
const DTGeometry * dtGeometry
Definition: TrackSplittingMonitor.h:70
iEvent
int iEvent
Definition: GenABIO.cc:224
reco::TrackBase::d0
double d0() const
dxy parameter in perigee convention (d0 = -dxy)
Definition: TrackBase.h:611
edm::EventSetup
Definition: EventSetup.h:58
HLT_FULL_cff.pt2
pt2
Definition: HLT_FULL_cff.py:9895
TrackSplittingMonitor_cfi.ddzBin
ddzBin
Definition: TrackSplittingMonitor_cfi.py:28
reco::TrackBase::normalizedChi2
double normalizedChi2() const
chi-squared divided by n.d.o.f. (or chi-squared * 1e6 if n.d.o.f. is zero)
Definition: TrackBase.h:593
TrackSplittingMonitor_cfi.normMin
normMin
Definition: TrackSplittingMonitor_cfi.py:49
TransientTrackRecord.h
TSCPBuilderNoMaterial.h
TrajectorySeedCollection.h
TrackSplittingMonitor::plotMuons_
bool plotMuons_
Definition: TrackSplittingMonitor.h:79
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
InputTag.h
TrackSplittingMonitor::dcurvNormalizedResiduals_tracker_
MonitorElement * dcurvNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:100
edm::EventSetup::getData
bool getData(T &iHolder) const
Definition: EventSetup.h:127
TrackSplittingMonitor_cfi.normBin
normBin
Definition: TrackSplittingMonitor_cfi.py:48
TrackSplittingMonitor::dthetaAbsoluteResiduals_tracker_
MonitorElement * dthetaAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:91
TrackSplittingMonitor_cfi.ddxyMax
ddxyMax
Definition: TrackSplittingMonitor_cfi.py:26
TrackSplittingMonitor::ddzAbsoluteResiduals_tracker_
MonitorElement * ddzAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:89
TrackSplittingMonitor::ddzNormalizedResiduals_global_
MonitorElement * ddzNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:110
TrackSplittingMonitor::ddxyNormalizedResiduals_global_
MonitorElement * ddxyNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:109
reco::TrackBase::d0Error
double d0Error() const
error on d0
Definition: TrackBase.h:772
TrackSplittingMonitor::tkGeomToken_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord > tkGeomToken_
Definition: TrackSplittingMonitor.h:63
TrackSplittingMonitor::rpcGeometry
const RPCGeometry * rpcGeometry
Definition: TrackSplittingMonitor.h:72
TrackSplittingMonitor_cfi.dthetaBin
dthetaBin
Definition: TrackSplittingMonitor_cfi.py:36
TrackSplittingMonitor_cfi.dphiBin
dphiBin
Definition: TrackSplittingMonitor_cfi.py:32
TrackSplittingMonitor::rpcGeomToken_
const edm::ESGetToken< RPCGeometry, MuonGeometryRecord > rpcGeomToken_
Definition: TrackSplittingMonitor.h:66
TrackSplittingMonitor::splitMuons_
edm::InputTag splitMuons_
Definition: TrackSplittingMonitor.h:75
TrackSplittingMonitor::splitTracks_
edm::InputTag splitTracks_
Definition: TrackSplittingMonitor.h:74
RecoMuonCosmics_cff.splitMuons
splitMuons
Definition: RecoMuonCosmics_cff.py:285
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
TrackSplittingMonitor::dcurvNormalizedResiduals_global_
MonitorElement * dcurvNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:114
TrackSplittingMonitor::dptNormalizedResiduals_global_
MonitorElement * dptNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:113
TrajectoryStateTransform.h
TrackSplittingMonitor_cfi.dphiMin
dphiMin
Definition: TrackSplittingMonitor_cfi.py:33
dqm::implementation::IBooker
Definition: DQMStore.h:43
reco::TrackBase::phiError
double phiError() const
error on phi
Definition: TrackBase.h:766
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
TrackSplittingMonitor_cfi.ddxyBin
ddxyBin
Definition: TrackSplittingMonitor_cfi.py:24
TrackSplittingMonitor_cfi.dptMax
dptMax
Definition: TrackSplittingMonitor_cfi.py:42
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
DeDxTools::esConsumes
ESGetTokenH3DDVariant esConsumes(std::string const &Reccord, edm::ConsumesCollector &)
Definition: DeDxTools.cc:283
edm::Event
Definition: Event.h:73
reco::TrackBase::theta
double theta() const
polar angle
Definition: TrackBase.h:602
dqm::impl::MonitorElement::setAxisTitle
virtual void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:791
edm::InputTag
Definition: InputTag.h:15
TrackSplittingMonitor::norchiCut_
double norchiCut_
Definition: TrackSplittingMonitor.h:85
TrackSplittingMonitor::ddxyNormalizedResiduals_tracker_
MonitorElement * ddxyNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:95
dqm::implementation::IBooker::book1D
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
Definition: DQMStore.h:98
TrackSplittingMonitor_cfi.dptBin
dptBin
Definition: TrackSplittingMonitor_cfi.py:40
TrackSplittingMonitor_cfi.dptMin
dptMin
Definition: TrackSplittingMonitor_cfi.py:41
TrackSplittingMonitor::dzCut_
double dzCut_
Definition: TrackSplittingMonitor.h:83
TrackSplittingMonitor::dthetaNormalizedResiduals_tracker_
MonitorElement * dthetaNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:98
Generated for CMSSW Reference Manual by doxygen 1.8.16