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  * $Date: 2010/02/11 00:11:01 $
5  * $Revision: 1.4 $
6  * \author Suchandra Dutta , Giorgia Mila
7  */
8 
9 #include "DataFormats/TrackReco/interface/Track.h"
10 #include "DataFormats/TrackReco/interface/TrackFwd.h"
11 #include "DataFormats/TrajectorySeed/interface/TrajectorySeedCollection.h"
12 #include "DataFormats/TrackCandidate/interface/TrackCandidateCollection.h"
13 #include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
14 #include "TrackingTools/Records/interface/TransientTrackRecord.h"
15 #include "TrackingTools/TransientTrack/interface/TransientTrack.h"
16 #include "MagneticField/Engine/interface/MagneticField.h"
17 
18 #include "FWCore/MessageLogger/interface/MessageLogger.h"
19 #include "FWCore/Utilities/interface/InputTag.h"
20 #include "DQMServices/Core/interface/DQMStore.h"
21 #include "DQM/TrackingMonitor/interface/TrackAnalyzer.h"
22 #include "DQM/TrackingMonitor/plugins/TrackSplittingMonitor.h"
23 
24 #include "FWCore/Framework/interface/ESHandle.h"
25 #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
26 
27 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"
28 #include "TrackingTools/Records/interface/TransientRecHitRecord.h"
29 #include "TrackingTools/PatternTools/interface/TSCBLBuilderNoMaterial.h"
30 #include "TrackingTools/PatternTools/interface/TSCPBuilderNoMaterial.h"
31 #include "DataFormats/BeamSpot/interface/BeamSpot.h"
32 #include <string>
33 
34 #include "DataFormats/MuonReco/interface/Muon.h"
35 #include "DataFormats/MuonReco/interface/MuonFwd.h"
36 
37 TrackSplittingMonitor::TrackSplittingMonitor(const edm::ParameterSet& iConfig) {
38  dqmStore_ = edm::Service<DQMStore>().operator->();
39  conf_ = iConfig;
40 }
41 
42 TrackSplittingMonitor::~TrackSplittingMonitor() {
43 }
44 
45 void TrackSplittingMonitor::beginJob(void) {
46 
47  std::string MEFolderName = conf_.getParameter<std::string>("FolderName");
48  dqmStore_->setCurrentFolder(MEFolderName);
49 
50  //get input tags
51  splitTracks_ = conf_.getParameter< edm::InputTag >("splitTrackCollection");
52  splitMuons_ = conf_.getParameter< edm::InputTag >("splitMuonCollection");
53  plotMuons_ = conf_.getParameter<bool>("ifPlotMuons");
54 
55  // cuts
56  pixelHitsPerLeg_ = conf_.getParameter<int>("pixelHitsPerLeg");
57  totalHitsPerLeg_ = conf_.getParameter<int>("totalHitsPerLeg");
58  d0Cut_ = conf_.getParameter<double>("d0Cut");
59  dzCut_ = conf_.getParameter<double>("dzCut");
60  ptCut_ = conf_.getParameter<double>("ptCut");
61  norchiCut_ = conf_.getParameter<double>("norchiCut");
62 
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_ = dqmStore_->book1D( "ddxyAbsoluteResiduals_tracker", "ddxyAbsoluteResiduals_tracker", ddxyBin, ddxyMin, ddxyMax );
95  ddzAbsoluteResiduals_tracker_ = dqmStore_->book1D( "ddzAbsoluteResiduals_tracker", "ddzAbsoluteResiduals_tracker", ddzBin, ddzMin, ddzMax );
96  dphiAbsoluteResiduals_tracker_ = dqmStore_->book1D( "dphiAbsoluteResiduals_tracker", "dphiAbsoluteResiduals_tracker", dphiBin, dphiMin, dphiMax );
97  dthetaAbsoluteResiduals_tracker_ = dqmStore_->book1D( "dthetaAbsoluteResiduals_tracker", "dthetaAbsoluteResiduals_tracker", dthetaBin, dthetaMin, dthetaMax );
98  dptAbsoluteResiduals_tracker_ = dqmStore_->book1D( "dptAbsoluteResiduals_tracker", "dptAbsoluteResiduals_tracker", dptBin, dptMin, dptMax );
99  dcurvAbsoluteResiduals_tracker_ = dqmStore_->book1D( "dcurvAbsoluteResiduals_tracker", "dcurvAbsoluteResiduals_tracker", dcurvBin, dcurvMin, dcurvMax );
100 
101  ddxyNormalizedResiduals_tracker_ = dqmStore_->book1D( "ddxyNormalizedResiduals_tracker", "ddxyNormalizedResiduals_tracker", normBin, normMin, normMax );
102  ddzNormalizedResiduals_tracker_ = dqmStore_->book1D( "ddzNormalizedResiduals_tracker", "ddzNormalizedResiduals_tracker", normBin, normMin, normMax );
103  dphiNormalizedResiduals_tracker_ = dqmStore_->book1D( "dphiNormalizedResiduals_tracker", "dphiNormalizedResiduals_tracker", normBin, normMin, normMax );
104  dthetaNormalizedResiduals_tracker_ = dqmStore_->book1D( "dthetaNormalizedResiduals_tracker", "dthetaNormalizedResiduals_tracker", normBin, normMin, normMax );
105  dptNormalizedResiduals_tracker_ = dqmStore_->book1D( "dptNormalizedResiduals_tracker", "dptNormalizedResiduals_tracker", normBin, normMin, normMax );
106  dcurvNormalizedResiduals_tracker_ = dqmStore_->book1D( "dcurvNormalizedResiduals_tracker", "dcurvNormalizedResiduals_tracker", normBin, normMin, normMax );
107 
108  if (plotMuons_){
109  ddxyAbsoluteResiduals_global_ = dqmStore_->book1D( "ddxyAbsoluteResiduals_global", "ddxyAbsoluteResiduals_global", ddxyBin, ddxyMin, ddxyMax );
110  ddzAbsoluteResiduals_global_ = dqmStore_->book1D( "ddzAbsoluteResiduals_global", "ddzAbsoluteResiduals_global", ddzBin, ddzMin, ddzMax );
111  dphiAbsoluteResiduals_global_ = dqmStore_->book1D( "dphiAbsoluteResiduals_global", "dphiAbsoluteResiduals_global", dphiBin, dphiMin, dphiMax );
112  dthetaAbsoluteResiduals_global_ = dqmStore_->book1D( "dthetaAbsoluteResiduals_global", "dthetaAbsoluteResiduals_global", dthetaBin, dthetaMin, dthetaMax );
113  dptAbsoluteResiduals_global_ = dqmStore_->book1D( "dptAbsoluteResiduals_global", "dptAbsoluteResiduals_global", dptBin, dptMin, dptMax );
114  dcurvAbsoluteResiduals_global_ = dqmStore_->book1D( "dcurvAbsoluteResiduals_global", "dcurvAbsoluteResiduals_global", dcurvBin, dcurvMin, dcurvMax );
115 
116  ddxyNormalizedResiduals_global_ = dqmStore_->book1D( "ddxyNormalizedResiduals_global", "ddxyNormalizedResiduals_global", normBin, normMin, normMax );
117  ddzNormalizedResiduals_global_ = dqmStore_->book1D( "ddzNormalizedResiduals_global", "ddzNormalizedResiduals_global", normBin, normMin, normMax );
118  dphiNormalizedResiduals_global_ = dqmStore_->book1D( "dphiNormalizedResiduals_global", "dphiNormalizedResiduals_global", normBin, normMin, normMax );
119  dthetaNormalizedResiduals_global_ = dqmStore_->book1D( "dthetaNormalizedResiduals_global", "dthetaNormalizedResiduals_global", normBin, normMin, normMax );
120  dptNormalizedResiduals_global_ = dqmStore_->book1D( "dptNormalizedResiduals_global", "dptNormalizedResiduals_global", normBin, normMin, normMax );
121  dcurvNormalizedResiduals_global_ = dqmStore_->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 // -- Analyse
158 //
159 void TrackSplittingMonitor::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
160 
161 
162  iSetup.get<IdealMagneticFieldRecord>().get(theMagField);
163  iSetup.get<TrackerDigiGeometryRecord>().get(theGeometry);
164  iSetup.get<MuonGeometryRecord>().get(dtGeometry);
165  iSetup.get<MuonGeometryRecord>().get(cscGeometry);
166  iSetup.get<MuonGeometryRecord>().get(rpcGeometry);
167 
168  edm::Handle<std::vector<reco::Track> > splitTracks;
169  iEvent.getByLabel(splitTracks_, splitTracks);
170  if (!splitTracks.isValid()) return;
171 
172  edm::Handle<std::vector<reco::Muon> > splitMuons;
173  if (plotMuons_){
174  iEvent.getByLabel(splitMuons_, splitMuons);
175  }
176 
177  if (splitTracks->size() == 2){
178  // check that there are 2 tracks in split track collection
179  edm::LogInfo("TrackSplittingMonitor") << "Split Track size: " << splitTracks->size();
180 
181  // split tracks calculations
182  reco::Track track1 = splitTracks->at(0);
183  reco::Track track2 = splitTracks->at(1);
184 
185 
186  // -------------------------- basic selection ---------------------------
187 
188  // hit counting
189  // looping through the hits for track 1
190  double nRechits1 =0;
191  double nRechitinBPIX1 =0;
192  for (trackingRecHit_iterator iHit = track1.recHitsBegin(); iHit != track1.recHitsEnd(); ++iHit) {
193  if((*iHit)->isValid()) {
194  nRechits1++;
195  int type =(*iHit)->geographicalId().subdetId();
196  if(type==int(PixelSubdetector::PixelBarrel)){++nRechitinBPIX1;}
197  }
198  }
199  // looping through the hits for track 2
200  double nRechits2 =0;
201  double nRechitinBPIX2 =0;
202  for (trackingRecHit_iterator iHit = track2.recHitsBegin(); iHit != track2.recHitsEnd(); ++iHit) {
203  if((*iHit)->isValid()) {
204  nRechits2++;
205  int type =(*iHit)->geographicalId().subdetId();
206  if(type==int(PixelSubdetector::PixelBarrel)){++nRechitinBPIX2;}
207  }
208  }
209 
210  // DCA of each track
211  double d01 = track1.d0();
212  double dz1 = track1.dz();
213  double d02 = track2.d0();
214  double dz2 = track2.dz();
215 
216  // pT of each track
217  double pt1 = track1.pt();
218  double pt2 = track2.pt();
219 
220  // chi2 of each track
221  double norchi1 = track1.normalizedChi2();
222  double norchi2 = track2.normalizedChi2();
223 
224  // basic selection
225  // pixel hits and total hits
226  if ((nRechitinBPIX1 >= pixelHitsPerLeg_)&&(nRechitinBPIX1 >= pixelHitsPerLeg_)&&(nRechits1 >= totalHitsPerLeg_)&&(nRechits2 >= totalHitsPerLeg_)){
227  // dca cut
228  if ( ((fabs(d01) < d0Cut_))&&(fabs(d02) < d0Cut_)&&(fabs(dz2) < dzCut_)&&(fabs(dz2) < dzCut_) ){
229  // pt cut
230  if ( (pt1+pt2)/2 < ptCut_){
231  // chi2 cut
232  if ( (norchi1 < norchiCut_)&&(norchi2 < norchiCut_) ){
233 
234  // passed all cuts...
235  edm::LogInfo("TrackSplittingMonitor") << " Setected after all cuts ?";
236 
237  double ddxyVal = d01 - d02;
238  double ddzVal = dz1 - dz2;
239  double dphiVal = track1.phi() - track2.phi();
240  double dthetaVal = track1.theta() - track2.theta();
241  double dptVal = pt1 - pt2;
242  double dcurvVal = (1/pt1) - (1/pt2);
243 
244  double d01ErrVal = track1.d0Error();
245  double d02ErrVal = track2.d0Error();
246  double dz1ErrVal = track1.dzError();
247  double dz2ErrVal = track2.dzError();
248  double phi1ErrVal = track1.phiError();
249  double phi2ErrVal = track2.phiError();
250  double theta1ErrVal = track1.thetaError();
251  double theta2ErrVal = track2.thetaError();
252  double pt1ErrVal = track1.ptError();
253  double pt2ErrVal = track2.ptError();
254 
255  ddxyAbsoluteResiduals_tracker_->Fill( 10000.0*ddxyVal/sqrt(2.0) );
256  ddxyAbsoluteResiduals_tracker_->Fill( 10000.0*ddzVal/sqrt(2.0) );
257  ddxyAbsoluteResiduals_tracker_->Fill( 1000.0*dphiVal/sqrt(2.0) );
258  ddxyAbsoluteResiduals_tracker_->Fill( 1000.0*dthetaVal/sqrt(2.0) );
259  ddxyAbsoluteResiduals_tracker_->Fill( dptVal/sqrt(2.0) );
260  ddxyAbsoluteResiduals_tracker_->Fill( dcurvVal/sqrt(2.0) );
261 
262  ddxyNormalizedResiduals_tracker_->Fill( ddxyVal/sqrt( d01ErrVal*d01ErrVal + d02ErrVal*d02ErrVal ) );
263  ddxyNormalizedResiduals_tracker_->Fill( ddzVal/sqrt( dz1ErrVal*dz1ErrVal + dz2ErrVal*dz2ErrVal ) );
264  ddxyNormalizedResiduals_tracker_->Fill( dphiVal/sqrt( phi1ErrVal*phi1ErrVal + phi2ErrVal*phi2ErrVal ) );
265  ddxyNormalizedResiduals_tracker_->Fill( dthetaVal/sqrt( theta1ErrVal*theta1ErrVal + theta2ErrVal*theta2ErrVal ) );
266  ddxyNormalizedResiduals_tracker_->Fill( dptVal/sqrt( pt1ErrVal*pt1ErrVal + pt2ErrVal*pt2ErrVal ) );
267  ddxyNormalizedResiduals_tracker_->Fill( dcurvVal/sqrt( pow(pt1ErrVal,2)/pow(pt1,4) + pow(pt2ErrVal,2)/pow(pt2,4) ) );
268 
269  // if do the same for split muons
270  if (plotMuons_ && splitMuons.isValid()){
271 
272  int gmCtr = 0;
273  bool topGlobalMuonFlag = false;
274  bool bottomGlobalMuonFlag = false;
275  int topGlobalMuon = -1;
276  int bottomGlobalMuon = -1;
277  double topGlobalMuonNorchi2 = 1e10;
278  double bottomGlobalMuonNorchi2 = 1e10;
279 
280  // check if usable split global muons
281  for (std::vector<reco::Muon>::const_iterator gmI = splitMuons->begin(); gmI != splitMuons->end(); gmI++){
282  if ( gmI->isTrackerMuon() && gmI->isStandAloneMuon() && gmI->isGlobalMuon() ){
283 
284  reco::TrackRef trackerTrackRef1( splitTracks, 0 );
285  reco::TrackRef trackerTrackRef2( splitTracks, 1 );
286 
287  if (gmI->innerTrack() == trackerTrackRef1){
288  if (gmI->globalTrack()->normalizedChi2() < topGlobalMuonNorchi2){
289  topGlobalMuonFlag = true;
290  topGlobalMuonNorchi2 = gmI->globalTrack()->normalizedChi2();
291  topGlobalMuon = gmCtr;
292  }
293  }
294  if (gmI->innerTrack() == trackerTrackRef2){
295  if (gmI->globalTrack()->normalizedChi2() < bottomGlobalMuonNorchi2){
296  bottomGlobalMuonFlag = true;
297  bottomGlobalMuonNorchi2 = gmI->globalTrack()->normalizedChi2();
298  bottomGlobalMuon = gmCtr;
299  }
300  }
301  }
302  gmCtr++;
303  }
304 
305  if (bottomGlobalMuonFlag && topGlobalMuonFlag) {
306 
307  reco::Muon muonTop = splitMuons->at( topGlobalMuon );
308  reco::Muon muonBottom = splitMuons->at( bottomGlobalMuon );
309 
310  reco::TrackRef glb1 = muonTop.globalTrack();
311  reco::TrackRef glb2 = muonBottom.globalTrack();
312 
313  double ddxyValGlb = glb1->d0() - glb2->d0();
314  double ddzValGlb = glb1->dz() - glb2->dz();
315  double dphiValGlb = glb1->phi() - glb2->phi();
316  double dthetaValGlb = glb1->theta() - glb2->theta();
317  double dptValGlb = glb1->pt() - glb2->pt();
318  double dcurvValGlb = (1/glb1->pt()) - (1/glb2->pt());
319 
320  double d01ErrValGlb = glb1->d0Error();
321  double d02ErrValGlb = glb2->d0Error();
322  double dz1ErrValGlb = glb1->dzError();
323  double dz2ErrValGlb = glb2->dzError();
324  double phi1ErrValGlb = glb1->phiError();
325  double phi2ErrValGlb = glb2->phiError();
326  double theta1ErrValGlb = glb1->thetaError();
327  double theta2ErrValGlb = glb2->thetaError();
328  double pt1ErrValGlb = glb1->ptError();
329  double pt2ErrValGlb = glb2->ptError();
330 
331  ddxyAbsoluteResiduals_global_->Fill( 10000.0*ddxyValGlb/sqrt(2.0) );
332  ddxyAbsoluteResiduals_global_->Fill( 10000.0*ddzValGlb/sqrt(2.0) );
333  ddxyAbsoluteResiduals_global_->Fill( 1000.0*dphiValGlb/sqrt(2.0) );
334  ddxyAbsoluteResiduals_global_->Fill( 1000.0*dthetaValGlb/sqrt(2.0) );
335  ddxyAbsoluteResiduals_global_->Fill( dptValGlb/sqrt(2.0) );
336  ddxyAbsoluteResiduals_global_->Fill( dcurvValGlb/sqrt(2.0) );
337 
338  ddxyNormalizedResiduals_global_->Fill( ddxyValGlb/sqrt( d01ErrValGlb*d01ErrValGlb + d02ErrValGlb*d02ErrValGlb ) );
339  ddxyNormalizedResiduals_global_->Fill( ddzValGlb/sqrt( dz1ErrValGlb*dz1ErrValGlb + dz2ErrValGlb*dz2ErrValGlb ) );
340  ddxyNormalizedResiduals_global_->Fill( dphiValGlb/sqrt( phi1ErrValGlb*phi1ErrValGlb + phi2ErrValGlb*phi2ErrValGlb ) );
341  ddxyNormalizedResiduals_global_->Fill( dthetaValGlb/sqrt( theta1ErrValGlb*theta1ErrValGlb + theta2ErrValGlb*theta2ErrValGlb ) );
342  ddxyNormalizedResiduals_global_->Fill( dptValGlb/sqrt( pt1ErrValGlb*pt1ErrValGlb + pt2ErrValGlb*pt2ErrValGlb ) );
343  ddxyNormalizedResiduals_global_->Fill( dcurvValGlb/sqrt( pow(pt1ErrValGlb,2)/pow(pt1,4) + pow(pt2ErrValGlb,2)/pow(pt2,4) ) );
344 
345  }
346 
347 
348  } // end of split muons loop
349  }
350  }
351  }
352  }
353  }
354 }
355 
356 
357 void TrackSplittingMonitor::endJob(void) {
358  bool outputMEsInRootFile = conf_.getParameter<bool>("OutputMEsInRootFile");
359  std::string outputFileName = conf_.getParameter<std::string>("OutputFileName");
360  if(outputMEsInRootFile){
361  dqmStore_->showDirStructure();
362  dqmStore_->save(outputFileName);
363  }
364 }
365 
366 DEFINE_FWK_MODULE(TrackSplittingMonitor);
type
type
Definition: HCALResponse.h:22
TrackSplittingMonitor::theGeometry
edm::ESHandle< TrackerGeometry > theGeometry
Definition: TrackSplittingMonitor.h:59
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
reco::TrackBase::d0Error
double d0Error() const
error on d0
Definition: TrackBase.h:204
TrackSplittingMonitor::dthetaAbsoluteResiduals_global_
MonitorElement * dthetaAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:96
TrackSplittingMonitor::dcurvNormalizedResiduals_global_
MonitorElement * dcurvNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:105
reco::TrackBase::d0
double d0() const
dxy parameter in perigee convention (d0 = - dxy)
Definition: TrackBase.h:122
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:519
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:110
TrackSplittingMonitor::dthetaAbsoluteResiduals_tracker_
MonitorElement * dthetaAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:82
TrackSplittingMonitor::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: TrackSplittingMonitor.cc:159
reco::TrackBase::theta
double theta() const
polar angle
Definition: TrackBase.h:116
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:1883
TrackSplittingMonitor::ddzNormalizedResiduals_global_
MonitorElement * ddzNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:101
TrackSplittingMonitor::ddxyNormalizedResiduals_tracker_
MonitorElement * ddxyNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:86
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:138
TrackSplittingMonitor::dcurvAbsoluteResiduals_tracker_
MonitorElement * dcurvAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:84
TrackSplittingMonitor::rpcGeometry
edm::ESHandle< RPCGeometry > rpcGeometry
Definition: TrackSplittingMonitor.h:63
TrackSplittingMonitor::ddzAbsoluteResiduals_global_
MonitorElement * ddzAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:94
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:147
TrackSplittingMonitor::dphiNormalizedResiduals_global_
MonitorElement * dphiNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:102
TrackSplittingMonitor::dthetaNormalizedResiduals_tracker_
MonitorElement * dthetaNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:89
iEvent
int iEvent
Definition: GenABIO.cc:243
TrackSplittingMonitor::ddxyNormalizedResiduals_global_
MonitorElement * ddxyNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:100
TrackSplittingMonitor::dthetaNormalizedResiduals_global_
MonitorElement * dthetaNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:103
TrackSplittingMonitor::dphiNormalizedResiduals_tracker_
MonitorElement * dphiNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:88
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
TrackSplittingMonitor::beginJob
virtual void beginJob(void)
Definition: TrackSplittingMonitor.cc:45
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:130
TrackSplittingMonitor::dcurvAbsoluteResiduals_global_
MonitorElement * dcurvAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:98
reco::TrackBase::ptError
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:187
TrackSplittingMonitor::dphiAbsoluteResiduals_tracker_
MonitorElement * dphiAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:81
reco::TrackBase::phiError
double phiError() const
error on phi
Definition: TrackBase.h:200
TrackSplittingMonitor::theMagField
edm::ESHandle< MagneticField > theMagField
Definition: TrackSplittingMonitor.h:60
TrackSplittingMonitor::ddzNormalizedResiduals_tracker_
MonitorElement * ddzNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:87
TrackSplittingMonitor::dptAbsoluteResiduals_global_
MonitorElement * dptAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:97
reco::Track::recHitsBegin
trackingRecHit_iterator recHitsBegin() const
Iterator to first hit on the track.
Definition: Track.h:63
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
TrackSplittingMonitor::dptAbsoluteResiduals_tracker_
MonitorElement * dptAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:83
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:359
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:126
reco::TrackBase::dzError
double dzError() const
error on dz
Definition: TrackBase.h:208
ExpressReco_HICollisions_FallBack.outputFileName
tuple outputFileName
Definition: ExpressReco_HICollisions_FallBack.py:11511
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
TrackSplittingMonitor::TrackSplittingMonitor
TrackSplittingMonitor(const edm::ParameterSet &)
Definition: TrackSplittingMonitor.cc:37
TrackSplittingMonitor::ddzAbsoluteResiduals_tracker_
MonitorElement * ddzAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:80
TrackSplittingMonitor::ddxyAbsoluteResiduals_global_
MonitorElement * ddxyAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:93
TrackSplittingMonitor::dtGeometry
edm::ESHandle< DTGeometry > dtGeometry
Definition: TrackSplittingMonitor.h:61
TrackSplittingMonitor::dcurvNormalizedResiduals_tracker_
MonitorElement * dcurvNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:91
TrackSplittingMonitor::dphiAbsoluteResiduals_global_
MonitorElement * dphiAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:95
TrackSplittingMonitor::dptNormalizedResiduals_tracker_
MonitorElement * dptNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:90
TrackSplittingMonitor::dptNormalizedResiduals_global_
MonitorElement * dptNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:104
TrackSplittingMonitor::cscGeometry
edm::ESHandle< CSCGeometry > cscGeometry
Definition: TrackSplittingMonitor.h:62
DQMStore::showDirStructure
void showDirStructure(void) const
Definition: DQMStore.cc:2539
RecoMuonCosmics_cff.splitMuons
tuple splitMuons
Definition: RecoMuonCosmics_cff.py:310
TrackSplittingMonitor::ddxyAbsoluteResiduals_tracker_
MonitorElement * ddxyAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:79
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:903
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:237
reco::TrackBase::thetaError
double thetaError() const
error on theta
Definition: TrackBase.h:194
reco::Muon::globalTrack
virtual TrackRef globalTrack() const
reference to Track reconstructed in both tracked and muon detector
Definition: Muon.h:44
reco::Track::recHitsEnd
trackingRecHit_iterator recHitsEnd() const
Iterator to last hit on the track.
Definition: Track.h:65
Generated for CMSSW Reference Manual by doxygen 1.8.5