|
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
|
#ifndef DQMOFFLINE_L1TRIGGER_L1TPHASE2MUONOFFLINE_H
#define DQMOFFLINE_L1TRIGGER_L1TPHASE2MUONOFFLINE_H
/**
* \file L1TPhase2MuonOffline.h
*
* \author S. Folgueras
*
*/
// DataFormats
#include "DataFormats/L1Trigger/interface/Muon.h"
#include "DataFormats/L1TMuonPhase2/interface/SAMuon.h"
#include "DataFormats/L1TMuonPhase2/interface/TrackerMuon.h"
#include "DataFormats/HepMCCandidate/interface/GenParticle.h"
#include "DataFormats/L1Trigger/interface/L1MuonParticleFwd.h"
#include "DataFormats/L1Trigger/interface/L1MuonParticle.h"
#include "DataFormats/Math/interface/deltaPhi.h"
// FWCore
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
// DQMServices
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include <memory>
#include "TRegexp.h"
#include <utility>
#include <vector>
class GenMuonGMTPair;
//
// DQM class declaration
//
class L1TPhase2MuonOffline : public DQMEDAnalyzer {
public:
L1TPhase2MuonOffline(const edm::ParameterSet& ps);
~L1TPhase2MuonOffline() override = default;
enum MuType { kSAMuon, kTkMuon, kNMuTypes };
enum VarType { kPt, kEta, kPhi, kIso, kQual, kZ0, kD0, kNVarTypes };
enum EffType { kEffPt, kEffPhi, kEffEta, kEffTypes };
enum ResType { kResPt, kRes1OverPt, kResQOverPt, kResPhi, kResEta, kResCh, kNResTypes };
enum EtaRegion { kEtaRegionAll, kEtaRegionBmtf, kEtaRegionOmtf, kEtaRegionEmtf, kNEtaRegions };
enum QualLevel { kQualOpen, kQualDouble, kQualSingle, kNQualLevels };
protected:
void dqmBeginRun(const edm::Run& run, const edm::EventSetup& iSetup) override;
void analyze(const edm::Event& e, const edm::EventSetup& c) override;
void bookControlHistos(DQMStore::IBooker&, MuType type);
void bookEfficiencyHistos(DQMStore::IBooker& ibooker, MuType type);
void bookResolutionHistos(DQMStore::IBooker& ibooker, MuType type);
void bookHistograms(DQMStore::IBooker& ibooker, const edm::Run& run, const edm::EventSetup& iSetup) override;
//Fill Histos
void fillControlHistos();
void fillEfficiencyHistos();
void fillResolutionHistos();
private:
// Cut and Matching
void getMuonGmtPairs(edm::Handle<l1t::MuonBxCollection>& gmtCands);
// Handles and Tokens
edm::EDGetTokenT<l1t::SAMuonCollection> gmtMuonToken_;
edm::EDGetTokenT<l1t::TrackerMuonCollection> gmtTkMuonToken_;
edm::EDGetTokenT<std::vector<reco::GenParticle>> genParticleToken_;
edm::Handle<l1t::SAMuonCollection> gmtSAMuon_;
edm::Handle<l1t::TrackerMuonCollection> gmtTkMuon_;
edm::Handle<std::vector<reco::GenParticle>> genparticles_;
// PropagateToMuon muonpropagator_;
// vectors of enum values to loop over (and store quantities)
const std::vector<MuType> muonTypes_;
const std::vector<EffType> effTypes_;
const std::vector<ResType> resTypes_;
// const std::vector<VarType> varTypes_;
const std::vector<EtaRegion> etaRegions_;
const std::vector<QualLevel> qualLevels_;
// maps with histogram name bits
// std::map<EffType, std::string> effNames_;
// std::map<EffType, std::string> effLabels_;
std::map<ResType, std::string> resNames_;
std::map<ResType, std::string> resLabels_;
std::map<EtaRegion, std::string> etaNames_;
std::map<QualLevel, std::string> qualNames_;
std::map<MuType, std::string> muonNames_;
// config params
const std::string histFolder_;
const std::vector<edm::ParameterSet> cutsVPSet_;
const std::vector<double> effVsPtBins_;
const std::vector<double> effVsPhiBins_;
const std::vector<double> effVsEtaBins_;
const double maxGmtMuonDR_;
// Helper methods
void matchMuonsToGen(std::vector<const reco::GenParticle*> genmus);
std::vector<float> getHistBinsEff(EffType eff);
std::tuple<int, double, double> getHistBinsRes(ResType res);
MonitorElement* efficiencyNum_[kNMuTypes][kNEtaRegions][kNQualLevels][kEffTypes];
MonitorElement* efficiencyDen_[kNMuTypes][kNEtaRegions][kNQualLevels][kEffTypes];
MonitorElement* resolutionHistos_[kNMuTypes][kNEtaRegions][kNQualLevels][kNResTypes];
MonitorElement* controlHistos_[kNMuTypes][kNVarTypes];
// helper variables
std::vector<GenMuonGMTPair> gmtSAMuonPairs_;
std::vector<GenMuonGMTPair> gmtTkMuonPairs_;
std::vector<std::pair<int, QualLevel>> cuts_;
const float lsb_pt = Phase2L1GMT::LSBpt;
const float lsb_phi = Phase2L1GMT::LSBphi;
const float lsb_eta = Phase2L1GMT::LSBeta;
const float lsb_z0 = Phase2L1GMT::LSBSAz0;
const float lsb_d0 = Phase2L1GMT::LSBSAd0;
};
//
// helper class to manage GMT-GenMuon pairing
//
class GenMuonGMTPair {
public:
GenMuonGMTPair(const reco::GenParticle* mu, const l1t::L1Candidate* gmtmu);
GenMuonGMTPair(const GenMuonGMTPair& muongmtPair);
GenMuonGMTPair& operator=(const GenMuonGMTPair& muongmtPair) = default;
~GenMuonGMTPair() {}
float dR2();
float pt() const { return mu_->pt(); };
float eta() const { return mu_->eta(); };
float phi() const { return mu_->phi(); };
int charge() const { return mu_->charge(); };
// Now properties of the L1 candidate:
float gmtPt() const { return gmtmu_ ? gmtmu_->pt() : -1.; };
float gmtEta() const { return gmtmu_ ? gmtEta_ : -5.; };
float gmtPhi() const { return gmtmu_ ? gmtPhi_ : -5.; };
int gmtCharge() const { return gmtmu_ ? gmtmu_->charge() : -5; };
int gmtQual() const { return gmtmu_ ? gmtmu_->hwQual() : -1; };
L1TPhase2MuonOffline::EtaRegion etaRegion() const;
double getDeltaVar(const L1TPhase2MuonOffline::ResType) const;
double getVar(const L1TPhase2MuonOffline::EffType) const;
private:
const reco::GenParticle* mu_;
const l1t::L1Candidate* gmtmu_;
// L1T muon eta and phi coordinates to be used
// Can be the coordinates from the 2nd muon station or from the vertex
float gmtEta_;
float gmtPhi_;
float muEta_;
float muPhi_;
};
#endif
|