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