|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
|
#ifndef TrackEfficiencyMonitor_H
#define TrackEfficiencyMonitor_H
// -*- C++ -*-
//
// Package: TrackEfficiencyMonitor
// Class: TrackEfficiencyMonitor
//
/**\class TrackEfficiencyMonitor TrackEfficiencyMonitor.cc DQM/TrackerMonitorTrack/src/TrackEfficiencyMonitor.cc
Monitoring source to measure the track efficiency
*/
// Original Author: Jeremy Andrea
// Insertion in DQM: Anne-Catherine Le Bihan
// Created: Thu 28 22:45:30 CEST 2008
// system includes
#include <memory>
#include <fstream>
// user includes
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/MuonReco/interface/MuonSelectors.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
#include "RecoMuon/GlobalTrackingTools/interface/DirectTrackerNavigation.h"
#include "RecoMuon/TrackingTools/interface/MuonServiceProxy.h"
#include "RecoTracker/MeasurementDet/interface/MeasurementTracker.h"
#include "RecoTracker/Record/interface/CkfComponentsRecord.h"
#include "RecoTracker/Record/interface/NavigationSchoolRecord.h"
#include "RecoTracker/Record/interface/TrackerRecoGeometryRecord.h"
#include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
namespace reco {
class TransientTrack;
}
class NavigationSchool;
class TrackEfficiencyMonitor : public DQMEDAnalyzer {
public:
typedef reco::Track Track;
typedef reco::TrackCollection TrackCollection;
explicit TrackEfficiencyMonitor(const edm::ParameterSet&);
~TrackEfficiencyMonitor() override;
void analyze(const edm::Event&, const edm::EventSetup&) override;
void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
enum SemiCylinder { Up, Down };
std::pair<TrajectoryStateOnSurface, const DetLayer*> findNextLayer(
TrajectoryStateOnSurface startTSOS, const std::vector<const DetLayer*>& trackCompatibleLayers, bool isUpMuon);
SemiCylinder checkSemiCylinder(const Track&);
void testTrackerTracks(edm::Handle<TrackCollection> tkTracks,
edm::Handle<TrackCollection> staTracks,
const NavigationSchool& navigationSchool);
void testSTATracks(edm::Handle<TrackCollection> tkTracks, edm::Handle<TrackCollection> staTracks);
bool trackerAcceptance(TrajectoryStateOnSurface theTSOS, double theRadius, double theMaxZ);
int compatibleLayers(const NavigationSchool& navigationSchool, TrajectoryStateOnSurface theTSOS);
private:
// ----------member data ---------------------------
std::string histname; //for naming the histograms
DQMStore* dqmStore_;
edm::ParameterSet conf_;
double theRadius_;
double theMaxZ_;
bool isBFieldOff_;
bool trackEfficiency_; //1 if one wants to measure the tracking efficiency
//0 if one wants to measure the muon reco efficiency
edm::InputTag theTKTracksLabel_;
edm::InputTag theSTATracksLabel_;
edm::EDGetTokenT<reco::TrackCollection> theTKTracksToken_;
edm::EDGetTokenT<reco::TrackCollection> theSTATracksToken_;
int failedToPropagate;
int nCompatibleLayers;
bool findDetLayer;
MonitorElement* muonX;
MonitorElement* muonY;
MonitorElement* muonZ;
MonitorElement* muonEta;
MonitorElement* muonPhi;
MonitorElement* muonD0;
MonitorElement* muonCompatibleLayers;
MonitorElement* trackX;
MonitorElement* trackY;
MonitorElement* trackZ;
MonitorElement* trackEta;
MonitorElement* trackPhi;
MonitorElement* trackD0;
MonitorElement* trackCompatibleLayers;
MonitorElement* deltaX;
MonitorElement* deltaY;
MonitorElement* signDeltaX;
MonitorElement* signDeltaY;
MonitorElement* GlobalMuonPtEtaPhiLowPt;
MonitorElement* StandaloneMuonPtEtaPhiLowPt;
MonitorElement* GlobalMuonPtEtaPhiHighPt;
MonitorElement* StandaloneMuonPtEtaPhiHighPt;
const DirectTrackerNavigation* theNavigation;
MuonServiceProxy* theMuonServiceProxy;
const edm::ESGetToken<TransientTrackBuilder, TransientTrackRecord> ttbToken_;
const edm::ESGetToken<Propagator, TrackingComponentsRecord> propToken_;
const edm::ESGetToken<NavigationSchool, NavigationSchoolRecord> navToken_;
const edm::ESGetToken<GeometricSearchTracker, TrackerRecoGeometryRecord> gstToken_;
const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> mfToken_;
const edm::ESGetToken<MeasurementTracker, CkfComponentsRecord> mtToken_;
edm::EDGetTokenT<edm::View<reco::Muon> > muonToken_;
const Propagator* thePropagator;
const TransientTrackBuilder* theTTrackBuilder;
const GeometricSearchTracker* theTracker;
const MagneticField* bField;
const MeasurementTracker* measurementTracker;
};
#endif
|