|
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
|
#ifndef ZCounting_H
#define ZCounting_H
#include "FWCore/Framework/interface/MakerMacros.h" // definitions for declaring plug-in modules
#include "FWCore/Framework/interface/Frameworkfwd.h" // declaration of EDM types
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h" // Parameters
#include "FWCore/Utilities/interface/InputTag.h"
#include <string> // string class
#include <cassert>
#include "DataFormats/HLTReco/interface/TriggerEvent.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/MuonReco/interface/MuonFwd.h"
#include "DataFormats/MuonReco/interface/MuonSelectors.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMOffline/Lumi/interface/TriggerTools.h"
class ZCounting : public DQMEDAnalyzer {
public:
ZCounting(const edm::ParameterSet& ps);
~ZCounting() override;
enum MuonIDTypes { NoneID, LooseID, MediumID, TightID, CustomTightID };
enum MuonIsoTypes { NoneIso, TrackerIso, PFIso };
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;
private:
//other functions
bool passMuonID(const reco::Muon& muon, const reco::Vertex* vtx);
bool passMuonIso(const reco::Muon& muon);
bool isCustomTightMuon(const reco::Muon& muon);
bool passGlobalMuon(const reco::Muon& muon);
bool passTrack(const reco::Track& track);
// EDM object collection names
const edm::InputTag triggerResultsInputTag_;
edm::EDGetTokenT<reco::VertexCollection> fPVName_token;
// Muons
edm::EDGetTokenT<reco::MuonCollection> fMuonName_token;
std::vector<std::string> fMuonHLTNames;
std::vector<std::string> fMuonHLTObjectNames;
edm::EDGetTokenT<reco::TrackCollection> fStandaloneRegName_token;
edm::EDGetTokenT<reco::TrackCollection> fStandaloneUpdName_token;
// Tracks
edm::EDGetTokenT<reco::TrackCollection> fTrackName_token;
// other input
const double PtCutL1_;
const double PtCutL2_;
const double EtaCutL1_;
const double EtaCutL2_;
const int MassBin_;
const double MassMin_;
const double MassMax_;
const int LumiBin_;
const double LumiMin_;
const double LumiMax_;
const int PVBin_;
const double PVMin_;
const double PVMax_;
const double VtxNTracksFitCut_;
const double VtxNdofCut_;
const double VtxAbsZCut_;
const double VtxRhoCut_;
const std::string IDTypestr_;
const std::string IsoTypestr_;
const double IsoCut_;
// muon ID and ISO parameters
MuonIDTypes IDType_{NoneID};
MuonIsoTypes IsoType_{NoneIso};
// trigger objects
HLTConfigProvider hltConfigProvider_;
TriggerTools* triggers;
// constants
const double DRMAX_HLT = 0.01; // max deltaR^2 matching between muon and hlt object
const double DRMAX_IO = 0.09; // max deltaR^2 matching between inner and outer muon track
const double MIN_PT_TRK = 15; // minimum pT of inner track considered for matching outer to inner track
const double MAX_ETA_TRK = 2.5; // maximum |eta| of inner track considered for matching outer to inner track
const double MIN_PT_STA = 15; // minimum pT of outer track considered for matching inner to outer track
const double MAX_ETA_STA = 2.4; // maximum |eta| of outer track considered for matching inner to outer track
const int N_STA_HITS = 1; // minimum number of valid standalone hits
const double MUON_MASS = 0.105658369;
const double MUON_BOUND = 0.9;
// General Histograms
MonitorElement* h_npv;
// Muon Histograms
MonitorElement* h_mass_2HLT_BB;
MonitorElement* h_mass_2HLT_BE;
MonitorElement* h_mass_2HLT_EE;
MonitorElement* h_mass_1HLT_BB;
MonitorElement* h_mass_1HLT_BE;
MonitorElement* h_mass_1HLT_EE;
MonitorElement* h_mass_ID_fail_BB;
MonitorElement* h_mass_ID_fail_BE;
MonitorElement* h_mass_ID_fail_EE;
MonitorElement* h_mass_Glo_pass_BB;
MonitorElement* h_mass_Glo_pass_BE;
MonitorElement* h_mass_Glo_pass_EE;
MonitorElement* h_mass_Glo_fail_BB;
MonitorElement* h_mass_Glo_fail_BE;
MonitorElement* h_mass_Glo_fail_EE;
MonitorElement* h_mass_Sta_pass_BB;
MonitorElement* h_mass_Sta_pass_BE;
MonitorElement* h_mass_Sta_pass_EE;
MonitorElement* h_mass_Sta_fail_BB;
MonitorElement* h_mass_Sta_fail_BE;
MonitorElement* h_mass_Sta_fail_EE;
};
#endif
|