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::bookHistograms(DQMStore::IBooker & ibooker,
48  edm::Run const & /* iRun */,
49  edm::EventSetup const & /* iSetup */)
50 
51 {
52 
53  std::string MEFolderName = conf_.getParameter<std::string>("FolderName");
54  ibooker.setCurrentFolder(MEFolderName);
55 
56  // bin declarations
57  int ddxyBin = conf_.getParameter<int>("ddxyBin");
58  double ddxyMin = conf_.getParameter<double>("ddxyMin");
59  double ddxyMax = conf_.getParameter<double>("ddxyMax");
60 
61  int ddzBin = conf_.getParameter<int>("ddzBin");
62  double ddzMin = conf_.getParameter<double>("ddzMin");
63  double ddzMax = conf_.getParameter<double>("ddzMax");
64 
65  int dphiBin = conf_.getParameter<int>("dphiBin");
66  double dphiMin = conf_.getParameter<double>("dphiMin");
67  double dphiMax = conf_.getParameter<double>("dphiMax");
68 
69  int dthetaBin = conf_.getParameter<int>("dthetaBin");
70  double dthetaMin = conf_.getParameter<double>("dthetaMin");
71  double dthetaMax = conf_.getParameter<double>("dthetaMax");
72 
73  int dptBin = conf_.getParameter<int>("dptBin");
74  double dptMin = conf_.getParameter<double>("dptMin");
75  double dptMax = conf_.getParameter<double>("dptMax");
76 
77  int dcurvBin = conf_.getParameter<int>("dcurvBin");
78  double dcurvMin = conf_.getParameter<double>("dcurvMin");
79  double dcurvMax = conf_.getParameter<double>("dcurvMax");
80 
81  int normBin = conf_.getParameter<int>("normBin");
82  double normMin = conf_.getParameter<double>("normMin");
83  double normMax = conf_.getParameter<double>("normMax");
84 
85  // declare histogram
86  ddxyAbsoluteResiduals_tracker_ = ibooker.book1D( "ddxyAbsoluteResiduals_tracker", "ddxyAbsoluteResiduals_tracker", ddxyBin, ddxyMin, ddxyMax );
87  ddzAbsoluteResiduals_tracker_ = ibooker.book1D( "ddzAbsoluteResiduals_tracker", "ddzAbsoluteResiduals_tracker", ddzBin, ddzMin, ddzMax );
88  dphiAbsoluteResiduals_tracker_ = ibooker.book1D( "dphiAbsoluteResiduals_tracker", "dphiAbsoluteResiduals_tracker", dphiBin, dphiMin, dphiMax );
89  dthetaAbsoluteResiduals_tracker_ = ibooker.book1D( "dthetaAbsoluteResiduals_tracker", "dthetaAbsoluteResiduals_tracker", dthetaBin, dthetaMin, dthetaMax );
90  dptAbsoluteResiduals_tracker_ = ibooker.book1D( "dptAbsoluteResiduals_tracker", "dptAbsoluteResiduals_tracker", dptBin, dptMin, dptMax );
91  dcurvAbsoluteResiduals_tracker_ = ibooker.book1D( "dcurvAbsoluteResiduals_tracker", "dcurvAbsoluteResiduals_tracker", dcurvBin, dcurvMin, dcurvMax );
92 
93  ddxyNormalizedResiduals_tracker_ = ibooker.book1D( "ddxyNormalizedResiduals_tracker", "ddxyNormalizedResiduals_tracker", normBin, normMin, normMax );
94  ddzNormalizedResiduals_tracker_ = ibooker.book1D( "ddzNormalizedResiduals_tracker", "ddzNormalizedResiduals_tracker", normBin, normMin, normMax );
95  dphiNormalizedResiduals_tracker_ = ibooker.book1D( "dphiNormalizedResiduals_tracker", "dphiNormalizedResiduals_tracker", normBin, normMin, normMax );
96  dthetaNormalizedResiduals_tracker_ = ibooker.book1D( "dthetaNormalizedResiduals_tracker", "dthetaNormalizedResiduals_tracker", normBin, normMin, normMax );
97  dptNormalizedResiduals_tracker_ = ibooker.book1D( "dptNormalizedResiduals_tracker", "dptNormalizedResiduals_tracker", normBin, normMin, normMax );
98  dcurvNormalizedResiduals_tracker_ = ibooker.book1D( "dcurvNormalizedResiduals_tracker", "dcurvNormalizedResiduals_tracker", normBin, normMin, normMax );
99 
100  if (plotMuons_){
101  ddxyAbsoluteResiduals_global_ = ibooker.book1D( "ddxyAbsoluteResiduals_global", "ddxyAbsoluteResiduals_global", ddxyBin, ddxyMin, ddxyMax );
102  ddzAbsoluteResiduals_global_ = ibooker.book1D( "ddzAbsoluteResiduals_global", "ddzAbsoluteResiduals_global", ddzBin, ddzMin, ddzMax );
103  dphiAbsoluteResiduals_global_ = ibooker.book1D( "dphiAbsoluteResiduals_global", "dphiAbsoluteResiduals_global", dphiBin, dphiMin, dphiMax );
104  dthetaAbsoluteResiduals_global_ = ibooker.book1D( "dthetaAbsoluteResiduals_global", "dthetaAbsoluteResiduals_global", dthetaBin, dthetaMin, dthetaMax );
105  dptAbsoluteResiduals_global_ = ibooker.book1D( "dptAbsoluteResiduals_global", "dptAbsoluteResiduals_global", dptBin, dptMin, dptMax );
106  dcurvAbsoluteResiduals_global_ = ibooker.book1D( "dcurvAbsoluteResiduals_global", "dcurvAbsoluteResiduals_global", dcurvBin, dcurvMin, dcurvMax );
107 
108  ddxyNormalizedResiduals_global_ = ibooker.book1D( "ddxyNormalizedResiduals_global", "ddxyNormalizedResiduals_global", normBin, normMin, normMax );
109  ddzNormalizedResiduals_global_ = ibooker.book1D( "ddzNormalizedResiduals_global", "ddzNormalizedResiduals_global", normBin, normMin, normMax );
110  dphiNormalizedResiduals_global_ = ibooker.book1D( "dphiNormalizedResiduals_global", "dphiNormalizedResiduals_global", normBin, normMin, normMax );
111  dthetaNormalizedResiduals_global_ = ibooker.book1D( "dthetaNormalizedResiduals_global", "dthetaNormalizedResiduals_global", normBin, normMin, normMax );
112  dptNormalizedResiduals_global_ = ibooker.book1D( "dptNormalizedResiduals_global", "dptNormalizedResiduals_global", normBin, normMin, normMax );
113  dcurvNormalizedResiduals_global_ = ibooker.book1D( "dcurvNormalizedResiduals_global", "dcurvNormalizedResiduals_global", normBin, normMin, normMax );
114  }
115 
116  ddxyAbsoluteResiduals_tracker_->setAxisTitle( "(#delta d_{xy})/#sqrt{2} [#mum]" );
117  ddxyAbsoluteResiduals_tracker_->setAxisTitle( "(#delta d_{z})/#sqrt{2} [#mum]" );
118  ddxyAbsoluteResiduals_tracker_->setAxisTitle( "(#delta #phi)/#sqrt{2} [mrad]" );
119  ddxyAbsoluteResiduals_tracker_->setAxisTitle( "(#delta #theta)/#sqrt{2} [mrad]" );
120  ddxyAbsoluteResiduals_tracker_->setAxisTitle( "(#delta pT)/#sqrt{2} [GeV]" );
121  ddxyAbsoluteResiduals_tracker_->setAxisTitle( "(#delta (1/pT))/#sqrt{2} [GeV^{-1}]" );
122 
123  ddxyNormalizedResiduals_tracker_->setAxisTitle( "#delta d_{xy}/#sigma(d_{xy}" );
124  ddxyNormalizedResiduals_tracker_->setAxisTitle( "#delta d_{z}/#sigma(d_{z})" );
125  ddxyNormalizedResiduals_tracker_->setAxisTitle( "#delta #phi/#sigma(d_{#phi})" );
126  ddxyNormalizedResiduals_tracker_->setAxisTitle( "#delta #theta/#sigma(d_{#theta})" );
127  ddxyNormalizedResiduals_tracker_->setAxisTitle( "#delta p_{T}/#sigma(p_{T})" );
128  ddxyNormalizedResiduals_tracker_->setAxisTitle( "#delta 1/p_{T}/#sigma(1/p_{T})" );
129 
130  if (plotMuons_){
131  ddxyAbsoluteResiduals_global_->setAxisTitle( "(#delta d_{xy})/#sqrt{2} [#mum]" );
132  ddxyAbsoluteResiduals_global_->setAxisTitle( "(#delta d_{z})/#sqrt{2} [#mum]" );
133  ddxyAbsoluteResiduals_global_->setAxisTitle( "(#delta #phi)/#sqrt{2} [mrad]" );
134  ddxyAbsoluteResiduals_global_->setAxisTitle( "(#delta #theta)/#sqrt{2} [mrad]" );
135  ddxyAbsoluteResiduals_global_->setAxisTitle( "(#delta pT)/#sqrt{2} [GeV]" );
136  ddxyAbsoluteResiduals_global_->setAxisTitle( "(#delta (1/pT))/#sqrt{2} [GeV^{-1}]" );
137 
138  ddxyNormalizedResiduals_global_->setAxisTitle( "#delta d_{xy}/#sigma(d_{xy}" );
139  ddxyNormalizedResiduals_global_->setAxisTitle( "#delta d_{z}/#sigma(d_{z})" );
140  ddxyNormalizedResiduals_global_->setAxisTitle( "#delta #phi/#sigma(d_{#phi})" );
141  ddxyNormalizedResiduals_global_->setAxisTitle( "#delta #theta/#sigma(d_{#theta})" );
142  ddxyNormalizedResiduals_global_->setAxisTitle( "#delta p_{T}/#sigma(p_{T})" );
143  ddxyNormalizedResiduals_global_->setAxisTitle( "#delta 1/p_{T}/#sigma(1/p_{T})" );
144  }
145 
146 }
147 
148 void TrackSplittingMonitor::beginJob(void)
149 {
150  //get input tags
151  plotMuons_ = conf_.getParameter<bool>("ifPlotMuons");
152 
153  // cuts
154  pixelHitsPerLeg_ = conf_.getParameter<int>("pixelHitsPerLeg");
155  totalHitsPerLeg_ = conf_.getParameter<int>("totalHitsPerLeg");
156  d0Cut_ = conf_.getParameter<double>("d0Cut");
157  dzCut_ = conf_.getParameter<double>("dzCut");
158  ptCut_ = conf_.getParameter<double>("ptCut");
159  norchiCut_ = conf_.getParameter<double>("norchiCut");
160 
161 }
162 
163 //
164 // -- Analyse
165 //
166 void TrackSplittingMonitor::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
167 
168 
169  iSetup.get<IdealMagneticFieldRecord>().get(theMagField);
170  iSetup.get<TrackerDigiGeometryRecord>().get(theGeometry);
171  iSetup.get<MuonGeometryRecord>().get(dtGeometry);
172  iSetup.get<MuonGeometryRecord>().get(cscGeometry);
173  iSetup.get<MuonGeometryRecord>().get(rpcGeometry);
174 
175  edm::Handle<std::vector<reco::Track> > splitTracks;
176  iEvent.getByToken(splitTracksToken_, splitTracks);
177  if (!splitTracks.isValid()) return;
178 
179  edm::Handle<std::vector<reco::Muon> > splitMuons;
180  if (plotMuons_){
181  iEvent.getByToken(splitMuonsToken_, splitMuons);
182  }
183 
184  if (splitTracks->size() == 2){
185  // check that there are 2 tracks in split track collection
186  edm::LogInfo("TrackSplittingMonitor") << "Split Track size: " << splitTracks->size();
187 
188  // split tracks calculations
189  reco::Track track1 = splitTracks->at(0);
190  reco::Track track2 = splitTracks->at(1);
191 
192 
193  // -------------------------- basic selection ---------------------------
194 
195  // hit counting
196  // looping through the hits for track 1
197  double nRechits1 =0;
198  double nRechitinBPIX1 =0;
199  for (trackingRecHit_iterator iHit = track1.recHitsBegin(); iHit != track1.recHitsEnd(); ++iHit) {
200  if((*iHit)->isValid()) {
201  nRechits1++;
202  int type =(*iHit)->geographicalId().subdetId();
203  if(type==int(PixelSubdetector::PixelBarrel)){++nRechitinBPIX1;}
204  }
205  }
206  // looping through the hits for track 2
207  double nRechits2 =0;
208  double nRechitinBPIX2 =0;
209  for (trackingRecHit_iterator iHit = track2.recHitsBegin(); iHit != track2.recHitsEnd(); ++iHit) {
210  if((*iHit)->isValid()) {
211  nRechits2++;
212  int type =(*iHit)->geographicalId().subdetId();
213  if(type==int(PixelSubdetector::PixelBarrel)){++nRechitinBPIX2;}
214  }
215  }
216 
217  // DCA of each track
218  double d01 = track1.d0();
219  double dz1 = track1.dz();
220  double d02 = track2.d0();
221  double dz2 = track2.dz();
222 
223  // pT of each track
224  double pt1 = track1.pt();
225  double pt2 = track2.pt();
226 
227  // chi2 of each track
228  double norchi1 = track1.normalizedChi2();
229  double norchi2 = track2.normalizedChi2();
230 
231  // basic selection
232  // pixel hits and total hits
233  if ((nRechitinBPIX1 >= pixelHitsPerLeg_)&&(nRechitinBPIX1 >= pixelHitsPerLeg_)&&(nRechits1 >= totalHitsPerLeg_)&&(nRechits2 >= totalHitsPerLeg_)){
234  // dca cut
235  if ( ((fabs(d01) < d0Cut_))&&(fabs(d02) < d0Cut_)&&(fabs(dz2) < dzCut_)&&(fabs(dz2) < dzCut_) ){
236  // pt cut
237  if ( (pt1+pt2)/2 < ptCut_){
238  // chi2 cut
239  if ( (norchi1 < norchiCut_)&&(norchi2 < norchiCut_) ){
240 
241  // passed all cuts...
242  edm::LogInfo("TrackSplittingMonitor") << " Setected after all cuts ?";
243 
244  double ddxyVal = d01 - d02;
245  double ddzVal = dz1 - dz2;
246  double dphiVal = track1.phi() - track2.phi();
247  double dthetaVal = track1.theta() - track2.theta();
248  double dptVal = pt1 - pt2;
249  double dcurvVal = (1/pt1) - (1/pt2);
250 
251  double d01ErrVal = track1.d0Error();
252  double d02ErrVal = track2.d0Error();
253  double dz1ErrVal = track1.dzError();
254  double dz2ErrVal = track2.dzError();
255  double phi1ErrVal = track1.phiError();
256  double phi2ErrVal = track2.phiError();
257  double theta1ErrVal = track1.thetaError();
258  double theta2ErrVal = track2.thetaError();
259  double pt1ErrVal = track1.ptError();
260  double pt2ErrVal = track2.ptError();
261 
262  ddxyAbsoluteResiduals_tracker_->Fill( 10000.0*ddxyVal/sqrt(2.0) );
263  ddxyAbsoluteResiduals_tracker_->Fill( 10000.0*ddzVal/sqrt(2.0) );
264  ddxyAbsoluteResiduals_tracker_->Fill( 1000.0*dphiVal/sqrt(2.0) );
265  ddxyAbsoluteResiduals_tracker_->Fill( 1000.0*dthetaVal/sqrt(2.0) );
266  ddxyAbsoluteResiduals_tracker_->Fill( dptVal/sqrt(2.0) );
267  ddxyAbsoluteResiduals_tracker_->Fill( dcurvVal/sqrt(2.0) );
268 
269  ddxyNormalizedResiduals_tracker_->Fill( ddxyVal/sqrt( d01ErrVal*d01ErrVal + d02ErrVal*d02ErrVal ) );
270  ddxyNormalizedResiduals_tracker_->Fill( ddzVal/sqrt( dz1ErrVal*dz1ErrVal + dz2ErrVal*dz2ErrVal ) );
271  ddxyNormalizedResiduals_tracker_->Fill( dphiVal/sqrt( phi1ErrVal*phi1ErrVal + phi2ErrVal*phi2ErrVal ) );
272  ddxyNormalizedResiduals_tracker_->Fill( dthetaVal/sqrt( theta1ErrVal*theta1ErrVal + theta2ErrVal*theta2ErrVal ) );
273  ddxyNormalizedResiduals_tracker_->Fill( dptVal/sqrt( pt1ErrVal*pt1ErrVal + pt2ErrVal*pt2ErrVal ) );
274  ddxyNormalizedResiduals_tracker_->Fill( dcurvVal/sqrt( pow(pt1ErrVal,2)/pow(pt1,4) + pow(pt2ErrVal,2)/pow(pt2,4) ) );
275 
276  // if do the same for split muons
277  if (plotMuons_ && splitMuons.isValid()){
278 
279  int gmCtr = 0;
280  bool topGlobalMuonFlag = false;
281  bool bottomGlobalMuonFlag = false;
282  int topGlobalMuon = -1;
283  int bottomGlobalMuon = -1;
284  double topGlobalMuonNorchi2 = 1e10;
285  double bottomGlobalMuonNorchi2 = 1e10;
286 
287  // check if usable split global muons
288  for (std::vector<reco::Muon>::const_iterator gmI = splitMuons->begin(); gmI != splitMuons->end(); gmI++){
289  if ( gmI->isTrackerMuon() && gmI->isStandAloneMuon() && gmI->isGlobalMuon() ){
290 
291  reco::TrackRef trackerTrackRef1( splitTracks, 0 );
292  reco::TrackRef trackerTrackRef2( splitTracks, 1 );
293 
294  if (gmI->innerTrack() == trackerTrackRef1){
295  if (gmI->globalTrack()->normalizedChi2() < topGlobalMuonNorchi2){
296  topGlobalMuonFlag = true;
297  topGlobalMuonNorchi2 = gmI->globalTrack()->normalizedChi2();
298  topGlobalMuon = gmCtr;
299  }
300  }
301  if (gmI->innerTrack() == trackerTrackRef2){
302  if (gmI->globalTrack()->normalizedChi2() < bottomGlobalMuonNorchi2){
303  bottomGlobalMuonFlag = true;
304  bottomGlobalMuonNorchi2 = gmI->globalTrack()->normalizedChi2();
305  bottomGlobalMuon = gmCtr;
306  }
307  }
308  }
309  gmCtr++;
310  }
311 
312  if (bottomGlobalMuonFlag && topGlobalMuonFlag) {
313 
314  reco::Muon muonTop = splitMuons->at( topGlobalMuon );
315  reco::Muon muonBottom = splitMuons->at( bottomGlobalMuon );
316 
317  reco::TrackRef glb1 = muonTop.globalTrack();
318  reco::TrackRef glb2 = muonBottom.globalTrack();
319 
320  double ddxyValGlb = glb1->d0() - glb2->d0();
321  double ddzValGlb = glb1->dz() - glb2->dz();
322  double dphiValGlb = glb1->phi() - glb2->phi();
323  double dthetaValGlb = glb1->theta() - glb2->theta();
324  double dptValGlb = glb1->pt() - glb2->pt();
325  double dcurvValGlb = (1/glb1->pt()) - (1/glb2->pt());
326 
327  double d01ErrValGlb = glb1->d0Error();
328  double d02ErrValGlb = glb2->d0Error();
329  double dz1ErrValGlb = glb1->dzError();
330  double dz2ErrValGlb = glb2->dzError();
331  double phi1ErrValGlb = glb1->phiError();
332  double phi2ErrValGlb = glb2->phiError();
333  double theta1ErrValGlb = glb1->thetaError();
334  double theta2ErrValGlb = glb2->thetaError();
335  double pt1ErrValGlb = glb1->ptError();
336  double pt2ErrValGlb = glb2->ptError();
337 
338  ddxyAbsoluteResiduals_global_->Fill( 10000.0*ddxyValGlb/sqrt(2.0) );
339  ddxyAbsoluteResiduals_global_->Fill( 10000.0*ddzValGlb/sqrt(2.0) );
340  ddxyAbsoluteResiduals_global_->Fill( 1000.0*dphiValGlb/sqrt(2.0) );
341  ddxyAbsoluteResiduals_global_->Fill( 1000.0*dthetaValGlb/sqrt(2.0) );
342  ddxyAbsoluteResiduals_global_->Fill( dptValGlb/sqrt(2.0) );
343  ddxyAbsoluteResiduals_global_->Fill( dcurvValGlb/sqrt(2.0) );
344 
345  ddxyNormalizedResiduals_global_->Fill( ddxyValGlb/sqrt( d01ErrValGlb*d01ErrValGlb + d02ErrValGlb*d02ErrValGlb ) );
346  ddxyNormalizedResiduals_global_->Fill( ddzValGlb/sqrt( dz1ErrValGlb*dz1ErrValGlb + dz2ErrValGlb*dz2ErrValGlb ) );
347  ddxyNormalizedResiduals_global_->Fill( dphiValGlb/sqrt( phi1ErrValGlb*phi1ErrValGlb + phi2ErrValGlb*phi2ErrValGlb ) );
348  ddxyNormalizedResiduals_global_->Fill( dthetaValGlb/sqrt( theta1ErrValGlb*theta1ErrValGlb + theta2ErrValGlb*theta2ErrValGlb ) );
349  ddxyNormalizedResiduals_global_->Fill( dptValGlb/sqrt( pt1ErrValGlb*pt1ErrValGlb + pt2ErrValGlb*pt2ErrValGlb ) );
350  ddxyNormalizedResiduals_global_->Fill( dcurvValGlb/sqrt( pow(pt1ErrValGlb,2)/pow(pt1,4) + pow(pt2ErrValGlb,2)/pow(pt2,4) ) );
351 
352  }
353 
354 
355  } // end of split muons loop
356  }
357  }
358  }
359  }
360  }
361 }
362 
363 
364 void TrackSplittingMonitor::endJob(void) {
365  bool outputMEsInRootFile = conf_.getParameter<bool>("OutputMEsInRootFile");
366  std::string outputFileName = conf_.getParameter<std::string>("OutputFileName");
367  if(outputMEsInRootFile){
368  dqmStore_->showDirStructure();
369  dqmStore_->save(outputFileName);
370  }
371 }
372 
373 DEFINE_FWK_MODULE(TrackSplittingMonitor);
type
type
Definition: HCALResponse.h:21
TrackSplittingMonitor::theGeometry
edm::ESHandle< TrackerGeometry > theGeometry
Definition: TrackSplittingMonitor.h:69
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:107
TrackSplittingMonitor::dcurvNormalizedResiduals_global_
MonitorElement * dcurvNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:116
reco::TrackBase::d0
double d0() const
dxy parameter in perigee convention (d0 = - dxy)
Definition: TrackBase.h:121
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:93
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:166
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:112
TrackSplittingMonitor::ddxyNormalizedResiduals_tracker_
MonitorElement * ddxyNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:97
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:137
TrackSplittingMonitor::dcurvAbsoluteResiduals_tracker_
MonitorElement * dcurvAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:95
TrackSplittingMonitor::rpcGeometry
edm::ESHandle< RPCGeometry > rpcGeometry
Definition: TrackSplittingMonitor.h:73
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
TrackSplittingMonitor::ddzAbsoluteResiduals_global_
MonitorElement * ddzAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:105
TrackSplittingMonitor::splitMuonsToken_
edm::EDGetTokenT< std::vector< reco::Muon > > splitMuonsToken_
Definition: TrackSplittingMonitor.h:78
TrackSplittingMonitor::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: TrackSplittingMonitor.cc:47
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:164
TrackSplittingMonitor::dphiNormalizedResiduals_global_
MonitorElement * dphiNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:113
dumpDBToFile_GT_ttrig_cfg.outputFileName
tuple outputFileName
Definition: dumpDBToFile_GT_ttrig_cfg.py:31
TrackSplittingMonitor::dthetaNormalizedResiduals_tracker_
MonitorElement * dthetaNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:100
iEvent
int iEvent
Definition: GenABIO.cc:230
TrackSplittingMonitor::ddxyNormalizedResiduals_global_
MonitorElement * ddxyNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:111
TrackSplittingMonitor::dthetaNormalizedResiduals_global_
MonitorElement * dthetaNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:114
TrackSplittingMonitor::dphiNormalizedResiduals_tracker_
MonitorElement * dphiNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:99
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:129
TrackSplittingMonitor::dcurvAbsoluteResiduals_global_
MonitorElement * dcurvAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:109
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:92
reco::TrackBase::phiError
double phiError() const
error on phi
Definition: TrackBase.h:205
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:113
TrackSplittingMonitor::theMagField
edm::ESHandle< MagneticField > theMagField
Definition: TrackSplittingMonitor.h:70
TrackSplittingMonitor::ddzNormalizedResiduals_tracker_
MonitorElement * ddzNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:98
TrackSplittingMonitor::dptAbsoluteResiduals_global_
MonitorElement * dptAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:108
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:94
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
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::ddzAbsoluteResiduals_tracker_
MonitorElement * ddzAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:91
TrackSplittingMonitor::ddxyAbsoluteResiduals_global_
MonitorElement * ddxyAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:104
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, const uint32_t lumi=0, SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE", const bool resetMEsAfterWriting=false)
Definition: DQMStore.cc:2490
TrackSplittingMonitor::dtGeometry
edm::ESHandle< DTGeometry > dtGeometry
Definition: TrackSplittingMonitor.h:71
TrackSplittingMonitor::dcurvNormalizedResiduals_tracker_
MonitorElement * dcurvNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:102
TrackSplittingMonitor::dphiAbsoluteResiduals_global_
MonitorElement * dphiAbsoluteResiduals_global_
Definition: TrackSplittingMonitor.h:106
TrackSplittingMonitor::splitTracksToken_
edm::EDGetTokenT< std::vector< reco::Track > > splitTracksToken_
Definition: TrackSplittingMonitor.h:77
TrackSplittingMonitor::dptNormalizedResiduals_tracker_
MonitorElement * dptNormalizedResiduals_tracker_
Definition: TrackSplittingMonitor.h:101
TrackSplittingMonitor::dptNormalizedResiduals_global_
MonitorElement * dptNormalizedResiduals_global_
Definition: TrackSplittingMonitor.h:115
TrackSplittingMonitor::cscGeometry
edm::ESHandle< CSCGeometry > cscGeometry
Definition: TrackSplittingMonitor.h:72
DQMStore::showDirStructure
void showDirStructure(void) const
Definition: DQMStore.cc:3282
RecoMuonCosmics_cff.splitMuons
tuple splitMuons
Definition: RecoMuonCosmics_cff.py:240
TrackSplittingMonitor::ddxyAbsoluteResiduals_tracker_
MonitorElement * ddxyAbsoluteResiduals_tracker_
Definition: TrackSplittingMonitor.h:90
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:957
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
edm::Run
Definition: Run.h:41
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:65
Generated for CMSSW Reference Manual by doxygen 1.8.5