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 
42 }
43 
44 TrackSplittingMonitor::~TrackSplittingMonitor() {
45 }
46 
47 void TrackSplittingMonitor::beginJob(void) {
48 
49  std::string MEFolderName = conf_.getParameter<std::string>("FolderName");
50  dqmStore_->setCurrentFolder(MEFolderName);
51 
52  //get input tags
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  iEvent.getByToken(splitTracksToken_, splitTracks);
171  if (!splitTracks.isValid()) return;
172 
173  edm::Handle<std::vector<reco::Muon> > splitMuons;
174  if (plotMuons_){
175  // iEvent.getByLabel(splitMuons_, splitMuons);
176  iEvent.getByToken(splitMuonsToken_, splitMuons);
177  }
178 
179  if (splitTracks->size() == 2){
180  // check that there are 2 tracks in split track collection
181  edm::LogInfo("TrackSplittingMonitor") << "Split Track size: " << splitTracks->size();
182 
183  // split tracks calculations
184  reco::Track track1 = splitTracks->at(0);
185  reco::Track track2 = splitTracks->at(1);
186 
187 
188  // -------------------------- basic selection ---------------------------
189 
190  // hit counting
191  // looping through the hits for track 1
192  double nRechits1 =0;
193  double nRechitinBPIX1 =0;
194  for (trackingRecHit_iterator iHit = track1.recHitsBegin(); iHit != track1.recHitsEnd(); ++iHit) {
195  if((*iHit)->isValid()) {
196  nRechits1++;
197  int type =(*iHit)->geographicalId().subdetId();
198  if(type==int(PixelSubdetector::PixelBarrel)){++nRechitinBPIX1;}
199  }
200  }
201  // looping through the hits for track 2
202  double nRechits2 =0;
203  double nRechitinBPIX2 =0;
204  for (trackingRecHit_iterator iHit = track2.recHitsBegin(); iHit != track2.recHitsEnd(); ++iHit) {
205  if((*iHit)->isValid()) {
206  nRechits2++;
207  int type =(*iHit)->geographicalId().subdetId();
208  if(type==int(PixelSubdetector::PixelBarrel)){++nRechitinBPIX2;}
209  }
210  }
211 
212  // DCA of each track
213  double d01 = track1.d0();
214  double dz1 = track1.dz();
215  double d02 = track2.d0();
216  double dz2 = track2.dz();
217 
218  // pT of each track
219  double pt1 = track1.pt();
220  double pt2 = track2.pt();
221 
222  // chi2 of each track
223  double norchi1 = track1.normalizedChi2();
224  double norchi2 = track2.normalizedChi2();
225 
226  // basic selection
227  // pixel hits and total hits
228  if ((nRechitinBPIX1 >= pixelHitsPerLeg_)&&(nRechitinBPIX1 >= pixelHitsPerLeg_)&&(nRechits1 >= totalHitsPerLeg_)&&(nRechits2 >= totalHitsPerLeg_)){
229  // dca cut
230  if ( ((fabs(d01) < d0Cut_))&&(fabs(d02) < d0Cut_)&&(fabs(dz2) < dzCut_)&&(fabs(dz2) < dzCut_) ){
231  // pt cut
232  if ( (pt1+pt2)/2 < ptCut_){
233  // chi2 cut
234  if ( (norchi1 < norchiCut_)&&(norchi2 < norchiCut_) ){
235 
236  // passed all cuts...
237  edm::LogInfo("TrackSplittingMonitor") << " Setected after all cuts ?";
238 
239  double ddxyVal = d01 - d02;
240  double ddzVal = dz1 - dz2;
241  double dphiVal = track1.phi() - track2.phi();
242  double dthetaVal = track1.theta() - track2.theta();
243  double dptVal = pt1 - pt2;
244  double dcurvVal = (1/pt1) - (1/pt2);
245 
246  double d01ErrVal = track1.d0Error();
247  double d02ErrVal = track2.d0Error();
248  double dz1ErrVal = track1.dzError();
249  double dz2ErrVal = track2.dzError();
250  double phi1ErrVal = track1.phiError();
251  double phi2ErrVal = track2.phiError();
252  double theta1ErrVal = track1.thetaError();
253  double theta2ErrVal = track2.thetaError();
254  double pt1ErrVal = track1.ptError();
255  double pt2ErrVal = track2.ptError();
256 
257  ddxyAbsoluteResiduals_tracker_->Fill( 10000.0*ddxyVal/sqrt(2.0) );
258  ddxyAbsoluteResiduals_tracker_->Fill( 10000.0*ddzVal/sqrt(2.0) );
259  ddxyAbsoluteResiduals_tracker_->Fill( 1000.0*dphiVal/sqrt(2.0) );
260  ddxyAbsoluteResiduals_tracker_->Fill( 1000.0*dthetaVal/sqrt(2.0) );
261  ddxyAbsoluteResiduals_tracker_->Fill( dptVal/sqrt(2.0) );
262  ddxyAbsoluteResiduals_tracker_->Fill( dcurvVal/sqrt(2.0) );
263 
264  ddxyNormalizedResiduals_tracker_->Fill( ddxyVal/sqrt( d01ErrVal*d01ErrVal + d02ErrVal*d02ErrVal ) );
265  ddxyNormalizedResiduals_tracker_->Fill( ddzVal/sqrt( dz1ErrVal*dz1ErrVal + dz2ErrVal*dz2ErrVal ) );
266  ddxyNormalizedResiduals_tracker_->Fill( dphiVal/sqrt( phi1ErrVal*phi1ErrVal + phi2ErrVal*phi2ErrVal ) );
267  ddxyNormalizedResiduals_tracker_->Fill( dthetaVal/sqrt( theta1ErrVal*theta1ErrVal + theta2ErrVal*theta2ErrVal ) );
268  ddxyNormalizedResiduals_tracker_->Fill( dptVal/sqrt( pt1ErrVal*pt1ErrVal + pt2ErrVal*pt2ErrVal ) );
269  ddxyNormalizedResiduals_tracker_->Fill( dcurvVal/sqrt( pow(pt1ErrVal,2)/pow(pt1,4) + pow(pt2ErrVal,2)/pow(pt2,4) ) );
270 
271  // if do the same for split muons
272  if (plotMuons_ && splitMuons.isValid()){
273 
274  int gmCtr = 0;
275  bool topGlobalMuonFlag = false;
276  bool bottomGlobalMuonFlag = false;
277  int topGlobalMuon = -1;
278  int bottomGlobalMuon = -1;
279  double topGlobalMuonNorchi2 = 1e10;
280  double bottomGlobalMuonNorchi2 = 1e10;
281 
282  // check if usable split global muons
283  for (std::vector<reco::Muon>::const_iterator gmI = splitMuons->begin(); gmI != splitMuons->end(); gmI++){
284  if ( gmI->isTrackerMuon() && gmI->isStandAloneMuon() && gmI->isGlobalMuon() ){
285 
286  reco::TrackRef trackerTrackRef1( splitTracks, 0 );
287  reco::TrackRef trackerTrackRef2( splitTracks, 1 );
288 
289  if (gmI->innerTrack() == trackerTrackRef1){
290  if (gmI->globalTrack()->normalizedChi2() < topGlobalMuonNorchi2){
291  topGlobalMuonFlag = true;
292  topGlobalMuonNorchi2 = gmI->globalTrack()->normalizedChi2();
293  topGlobalMuon = gmCtr;
294  }
295  }
296  if (gmI->innerTrack() == trackerTrackRef2){
297  if (gmI->globalTrack()->normalizedChi2() < bottomGlobalMuonNorchi2){
298  bottomGlobalMuonFlag = true;
299  bottomGlobalMuonNorchi2 = gmI->globalTrack()->normalizedChi2();
300  bottomGlobalMuon = gmCtr;
301  }
302  }
303  }
304  gmCtr++;
305  }
306 
307  if (bottomGlobalMuonFlag && topGlobalMuonFlag) {
308 
309  reco::Muon muonTop = splitMuons->at( topGlobalMuon );
310  reco::Muon muonBottom = splitMuons->at( bottomGlobalMuon );
311 
312  reco::TrackRef glb1 = muonTop.globalTrack();
313  reco::TrackRef glb2 = muonBottom.globalTrack();
314 
315  double ddxyValGlb = glb1->d0() - glb2->d0();
316  double ddzValGlb = glb1->dz() - glb2->dz();
317  double dphiValGlb = glb1->phi() - glb2->phi();
318  double dthetaValGlb = glb1->theta() - glb2->theta();
319  double dptValGlb = glb1->pt() - glb2->pt();
320  double dcurvValGlb = (1/glb1->pt()) - (1/glb2->pt());
321 
322  double d01ErrValGlb = glb1->d0Error();
323  double d02ErrValGlb = glb2->d0Error();
324  double dz1ErrValGlb = glb1->dzError();
325  double dz2ErrValGlb = glb2->dzError();
326  double phi1ErrValGlb = glb1->phiError();
327  double phi2ErrValGlb = glb2->phiError();
328  double theta1ErrValGlb = glb1->thetaError();
329  double theta2ErrValGlb = glb2->thetaError();
330  double pt1ErrValGlb = glb1->ptError();
331  double pt2ErrValGlb = glb2->ptError();
332 
333  ddxyAbsoluteResiduals_global_->Fill( 10000.0*ddxyValGlb/sqrt(2.0) );
334  ddxyAbsoluteResiduals_global_->Fill( 10000.0*ddzValGlb/sqrt(2.0) );
335  ddxyAbsoluteResiduals_global_->Fill( 1000.0*dphiValGlb/sqrt(2.0) );
336  ddxyAbsoluteResiduals_global_->Fill( 1000.0*dthetaValGlb/sqrt(2.0) );
337  ddxyAbsoluteResiduals_global_->Fill( dptValGlb/sqrt(2.0) );
338  ddxyAbsoluteResiduals_global_->Fill( dcurvValGlb/sqrt(2.0) );
339 
340  ddxyNormalizedResiduals_global_->Fill( ddxyValGlb/sqrt( d01ErrValGlb*d01ErrValGlb + d02ErrValGlb*d02ErrValGlb ) );
341  ddxyNormalizedResiduals_global_->Fill( ddzValGlb/sqrt( dz1ErrValGlb*dz1ErrValGlb + dz2ErrValGlb*dz2ErrValGlb ) );
342  ddxyNormalizedResiduals_global_->Fill( dphiValGlb/sqrt( phi1ErrValGlb*phi1ErrValGlb + phi2ErrValGlb*phi2ErrValGlb ) );
343  ddxyNormalizedResiduals_global_->Fill( dthetaValGlb/sqrt( theta1ErrValGlb*theta1ErrValGlb + theta2ErrValGlb*theta2ErrValGlb ) );
344  ddxyNormalizedResiduals_global_->Fill( dptValGlb/sqrt( pt1ErrValGlb*pt1ErrValGlb + pt2ErrValGlb*pt2ErrValGlb ) );
345  ddxyNormalizedResiduals_global_->Fill( dcurvValGlb/sqrt( pow(pt1ErrValGlb,2)/pow(pt1,4) + pow(pt2ErrValGlb,2)/pow(pt2,4) ) );
346 
347  }
348 
349 
350  } // end of split muons loop
351  }
352  }
353  }
354  }
355  }
356 }
357 
358 
359 void TrackSplittingMonitor::endJob(void) {
360  bool outputMEsInRootFile = conf_.getParameter<bool>("OutputMEsInRootFile");
361  std::string outputFileName = conf_.getParameter<std::string>("OutputFileName");
362  if(outputMEsInRootFile){
363  dqmStore_->showDirStructure();
364  dqmStore_->save(outputFileName);
365  }
366 }
367 
368 DEFINE_FWK_MODULE(TrackSplittingMonitor);
type
type
Definition: HCALResponse.h:21
TrackSplittingMonitor::theGeometry
edm::ESHandle< TrackerGeometry > theGeometry
Definition: TrackSplittingMonitor.h:64
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
reco::TrackBase::d0Error
double d0Error() const
error on d0
Definition: TrackBase.h:209
TrackSplittingMonitor::dthetaAbsoluteResiduals_global_
MonitorElement * dthetaAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:102
TrackSplittingMonitor::dcurvNormalizedResiduals_global_
MonitorElement * dcurvNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:111
reco::TrackBase::d0
double d0() const
dxy parameter in perigee convention (d0 = - dxy)
Definition: TrackBase.h:121
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:872
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:109
TrackSplittingMonitor::dthetaAbsoluteResiduals_tracker_
MonitorElement * dthetaAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:88
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
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:115
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
TrackSplittingMonitor::ddzNormalizedResiduals_global_
MonitorElement * ddzNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:107
TrackSplittingMonitor::ddxyNormalizedResiduals_tracker_
MonitorElement * ddxyNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:92
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:137
TrackSplittingMonitor::dcurvAbsoluteResiduals_tracker_
MonitorElement * dcurvAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:90
TrackSplittingMonitor::rpcGeometry
edm::ESHandle< RPCGeometry > rpcGeometry
Definition: TrackSplittingMonitor.h:68
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
TrackSplittingMonitor::ddzAbsoluteResiduals_global_
MonitorElement * ddzAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:100
TrackSplittingMonitor::splitMuonsToken_
edm::EDGetTokenT< std::vector< reco::Muon > > splitMuonsToken_
Definition: TrackSplittingMonitor.h:73
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:161
TrackSplittingMonitor::dphiNormalizedResiduals_global_
MonitorElement * dphiNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:108
dumpDBToFile_GT_ttrig_cfg.outputFileName
tuple outputFileName
Definition: dumpDBToFile_GT_ttrig_cfg.py:31
TrackSplittingMonitor::dthetaNormalizedResiduals_tracker_
MonitorElement * dthetaNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:95
iEvent
int iEvent
Definition: GenABIO.cc:243
TrackSplittingMonitor::ddxyNormalizedResiduals_global_
MonitorElement * ddxyNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:106
TrackSplittingMonitor::dthetaNormalizedResiduals_global_
MonitorElement * dthetaNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:109
TrackSplittingMonitor::dphiNormalizedResiduals_tracker_
MonitorElement * dphiNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:94
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
TrackSplittingMonitor::beginJob
virtual void beginJob(void)
Definition: TrackSplittingMonitor.cc:47
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:129
TrackSplittingMonitor::dcurvAbsoluteResiduals_global_
MonitorElement * dcurvAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:104
reco::TrackBase::ptError
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:192
TrackSplittingMonitor::dphiAbsoluteResiduals_tracker_
MonitorElement * dphiAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:87
reco::TrackBase::phiError
double phiError() const
error on phi
Definition: TrackBase.h:205
TrackSplittingMonitor::theMagField
edm::ESHandle< MagneticField > theMagField
Definition: TrackSplittingMonitor.h:65
TrackSplittingMonitor::ddzNormalizedResiduals_tracker_
MonitorElement * ddzNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:93
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", const uint32_t run=0, SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:2296
TrackSplittingMonitor::dptAbsoluteResiduals_global_
MonitorElement * dptAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:103
reco::Track::recHitsBegin
trackingRecHit_iterator recHitsBegin() const
Iterator to first hit on the track.
Definition: Track.h:62
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
TrackSplittingMonitor::dptAbsoluteResiduals_tracker_
MonitorElement * dptAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:89
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:125
reco::TrackBase::dzError
double dzError() const
error on dz
Definition: TrackBase.h:213
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
TrackSplittingMonitor::TrackSplittingMonitor
TrackSplittingMonitor(const edm::ParameterSet &)
Definition: TrackSplittingMonitor.cc:32
TrackSplittingMonitor::ddzAbsoluteResiduals_tracker_
MonitorElement * ddzAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:86
TrackSplittingMonitor::ddxyAbsoluteResiduals_global_
MonitorElement * ddxyAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:99
TrackSplittingMonitor::dtGeometry
edm::ESHandle< DTGeometry > dtGeometry
Definition: TrackSplittingMonitor.h:66
TrackSplittingMonitor::dcurvNormalizedResiduals_tracker_
MonitorElement * dcurvNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:97
TrackSplittingMonitor::dphiAbsoluteResiduals_global_
MonitorElement * dphiAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:101
TrackSplittingMonitor::splitTracksToken_
edm::EDGetTokenT< std::vector< reco::Track > > splitTracksToken_
Definition: TrackSplittingMonitor.h:72
TrackSplittingMonitor::dptNormalizedResiduals_tracker_
MonitorElement * dptNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:96
TrackSplittingMonitor::dptNormalizedResiduals_global_
MonitorElement * dptNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:110
TrackSplittingMonitor::cscGeometry
edm::ESHandle< CSCGeometry > cscGeometry
Definition: TrackSplittingMonitor.h:67
DQMStore::showDirStructure
void showDirStructure(void) const
Definition: DQMStore.cc:2961
RecoMuonCosmics_cff.splitMuons
tuple splitMuons
Definition: RecoMuonCosmics_cff.py:240
TrackSplittingMonitor::ddxyAbsoluteResiduals_tracker_
MonitorElement * ddxyAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:85
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:928
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:584
reco::TrackBase::thetaError
double thetaError() const
error on theta
Definition: TrackBase.h:199
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:64
Generated for CMSSW Reference Manual by doxygen 1.8.5