|
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
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
|
#ifndef L1TEGammaOffline_H
#define L1TEGammaOffline_H
//Framework
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
//event
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
//DQM
#include "DQMServices/Core/interface/DQMOneEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMOffline/L1Trigger/interface/HistDefinition.h"
//Candidate handling
#include "DataFormats/Candidate/interface/Candidate.h"
// Electron & photon collections
#include "DataFormats/EgammaCandidates/interface/GsfElectron.h"
#include "DataFormats/EgammaCandidates/interface/Photon.h"
// Vertex utilities
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
// Trigger
#include "DataFormats/Common/interface/TriggerResults.h"
#include "DataFormats/HLTReco/interface/TriggerObject.h"
#include "DataFormats/HLTReco/interface/TriggerEvent.h"
// stage2 collections:
#include "DataFormats/L1Trigger/interface/EGamma.h"
class L1TEGammaOffline : public DQMOneEDAnalyzer<> {
public:
L1TEGammaOffline(const edm::ParameterSet& ps);
~L1TEGammaOffline() override;
enum PlotConfig { nVertex, ETvsET, PHIvsPHI };
static const std::map<std::string, unsigned int> PlotConfigNames;
protected:
void dqmBeginRun(edm::Run const&, edm::EventSetup const&) override;
void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
void analyze(edm::Event const& e, edm::EventSetup const& eSetup) override;
void dqmEndRun(edm::Run const& run, edm::EventSetup const& eSetup) override;
private:
bool passesLooseEleId(reco::GsfElectron const& electron) const;
bool passesMediumEleId(reco::GsfElectron const& electron) const;
void bookElectronHistos(DQMStore::IBooker&);
void bookPhotonHistos(DQMStore::IBooker&);
//other functions
double Distance(const reco::Candidate& c1, const reco::Candidate& c2);
double DistancePhi(const reco::Candidate& c1, const reco::Candidate& c2);
double calcDeltaPhi(double phi1, double phi2);
void fillElectrons(edm::Event const& e, const unsigned int nVertex);
void fillPhotons(edm::Event const& e, const unsigned int nVertex);
bool findTagAndProbePair(edm::Handle<reco::GsfElectronCollection> const& electrons);
bool matchesAnHLTObject(double eta, double phi) const;
void normalise2DHistogramsToBinArea();
math::XYZPoint PVPoint_;
//variables from config file
edm::EDGetTokenT<reco::GsfElectronCollection> theGsfElectronCollection_;
edm::EDGetTokenT<std::vector<reco::Photon> > thePhotonCollection_;
edm::EDGetTokenT<reco::VertexCollection> thePVCollection_;
edm::EDGetTokenT<reco::BeamSpot> theBSCollection_;
edm::EDGetTokenT<trigger::TriggerEvent> triggerInputTag_;
edm::EDGetTokenT<edm::TriggerResults> triggerResultsInputTag_;
std::string triggerProcess_;
std::vector<std::string> triggerNames_;
std::string histFolder_;
std::string efficiencyFolder_;
edm::EDGetTokenT<l1t::EGammaBxCollection> stage2CaloLayer2EGammaToken_;
std::vector<double> electronEfficiencyThresholds_;
std::vector<double> electronEfficiencyBins_;
double probeToL1Offset_;
std::vector<double> deepInspectionElectronThresholds_;
std::vector<double> photonEfficiencyThresholds_;
std::vector<double> photonEfficiencyBins_;
double maxDeltaRForL1Matching_;
double maxDeltaRForHLTMatching_;
double recoToL1TThresholdFactor_;
reco::GsfElectron tagElectron_;
reco::GsfElectron probeElectron_;
double tagAndProbleInvariantMass_;
HLTConfigProvider hltConfig_;
std::vector<unsigned int> triggerIndices_;
edm::TriggerResults triggerResults_;
trigger::TriggerEvent triggerEvent_;
dqmoffline::l1t::HistDefinitions histDefinitions_;
// TODO: add turn-on cuts (vectors of doubles)
// Histograms
MonitorElement* h_nVertex_;
MonitorElement* h_tagAndProbeMass_;
// electron reco vs L1
MonitorElement* h_L1EGammaETvsElectronET_EB_;
MonitorElement* h_L1EGammaETvsElectronET_EE_;
MonitorElement* h_L1EGammaETvsElectronET_EB_EE_;
MonitorElement* h_L1EGammaPhivsElectronPhi_EB_;
MonitorElement* h_L1EGammaPhivsElectronPhi_EE_;
MonitorElement* h_L1EGammaPhivsElectronPhi_EB_EE_;
MonitorElement* h_L1EGammaEtavsElectronEta_;
// electron resolutions
MonitorElement* h_resolutionElectronET_EB_;
MonitorElement* h_resolutionElectronET_EE_;
MonitorElement* h_resolutionElectronET_EB_EE_;
MonitorElement* h_resolutionElectronPhi_EB_;
MonitorElement* h_resolutionElectronPhi_EE_;
MonitorElement* h_resolutionElectronPhi_EB_EE_;
MonitorElement* h_resolutionElectronEta_;
// electron turn-ons
std::map<double, MonitorElement*> h_efficiencyElectronET_EB_pass_;
std::map<double, MonitorElement*> h_efficiencyElectronET_EE_pass_;
std::map<double, MonitorElement*> h_efficiencyElectronET_EB_EE_pass_;
std::map<double, MonitorElement*> h_efficiencyElectronPhi_vs_Eta_pass_;
// for deep inspection only
std::map<double, MonitorElement*> h_efficiencyElectronEta_pass_;
std::map<double, MonitorElement*> h_efficiencyElectronPhi_pass_;
std::map<double, MonitorElement*> h_efficiencyElectronNVertex_pass_;
// we could drop the map here, but L1TEfficiency_Harvesting expects
// identical names except for the suffix
std::map<double, MonitorElement*> h_efficiencyElectronET_EB_total_;
std::map<double, MonitorElement*> h_efficiencyElectronET_EE_total_;
std::map<double, MonitorElement*> h_efficiencyElectronET_EB_EE_total_;
std::map<double, MonitorElement*> h_efficiencyElectronPhi_vs_Eta_total_;
// for deep inspection only
std::map<double, MonitorElement*> h_efficiencyElectronEta_total_;
std::map<double, MonitorElement*> h_efficiencyElectronPhi_total_;
std::map<double, MonitorElement*> h_efficiencyElectronNVertex_total_;
// photons
MonitorElement* h_L1EGammaETvsPhotonET_EB_;
MonitorElement* h_L1EGammaETvsPhotonET_EE_;
MonitorElement* h_L1EGammaETvsPhotonET_EB_EE_;
MonitorElement* h_L1EGammaPhivsPhotonPhi_EB_;
MonitorElement* h_L1EGammaPhivsPhotonPhi_EE_;
MonitorElement* h_L1EGammaPhivsPhotonPhi_EB_EE_;
MonitorElement* h_L1EGammaEtavsPhotonEta_;
// electron resolutions
MonitorElement* h_resolutionPhotonET_EB_;
MonitorElement* h_resolutionPhotonET_EE_;
MonitorElement* h_resolutionPhotonET_EB_EE_;
MonitorElement* h_resolutionPhotonPhi_EB_;
MonitorElement* h_resolutionPhotonPhi_EE_;
MonitorElement* h_resolutionPhotonPhi_EB_EE_;
MonitorElement* h_resolutionPhotonEta_;
// photon turn-ons
std::map<double, MonitorElement*> h_efficiencyPhotonET_EB_pass_;
std::map<double, MonitorElement*> h_efficiencyPhotonET_EE_pass_;
std::map<double, MonitorElement*> h_efficiencyPhotonET_EB_EE_pass_;
// we could drop the map here, but L1TEfficiency_Harvesting expects
// identical names except for the suffix
std::map<double, MonitorElement*> h_efficiencyPhotonET_EB_total_;
std::map<double, MonitorElement*> h_efficiencyPhotonET_EE_total_;
std::map<double, MonitorElement*> h_efficiencyPhotonET_EB_EE_total_;
};
#endif
|