|
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
|
#include "DataFormats/HLTReco/interface/TriggerObject.h"
#include "FWCore/Common/interface/TriggerNames.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "HLTriggerOffline/SUSYBSM/interface/SUSY_HLT_PhotonHT.h"
SUSY_HLT_PhotonHT::SUSY_HLT_PhotonHT(const edm::ParameterSet &ps) {
edm::LogInfo("SUSY_HLT_PhotonHT") << "Constructor SUSY_HLT_PhotonHT::SUSY_HLT_PhotonHT " << std::endl;
// Get parameters from configuration file
theTrigSummary_ = consumes<trigger::TriggerEvent>(ps.getParameter<edm::InputTag>("trigSummary"));
thePfMETCollection_ = consumes<reco::PFMETCollection>(ps.getParameter<edm::InputTag>("pfMETCollection"));
thePhotonCollection_ = consumes<reco::PhotonCollection>(ps.getParameter<edm::InputTag>("photonCollection"));
triggerResults_ = consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>("TriggerResults"));
triggerPath_ = ps.getParameter<std::string>("TriggerPath");
triggerPathAuxiliaryForHadronic_ = ps.getParameter<std::string>("TriggerPathAuxiliaryForHadronic");
triggerFilterPhoton_ = ps.getParameter<edm::InputTag>("TriggerFilterPhoton");
triggerFilterHt_ = ps.getParameter<edm::InputTag>("TriggerFilterHt");
ptThrOffline_ = ps.getUntrackedParameter<double>("ptThrOffline");
htThrOffline_ = ps.getUntrackedParameter<double>("htThrOffline");
}
SUSY_HLT_PhotonHT::~SUSY_HLT_PhotonHT() {
edm::LogInfo("SUSY_HLT_PhotonHT") << "Destructor SUSY_HLT_PhotonHT::~SUSY_HLT_PhotonHT " << std::endl;
}
void SUSY_HLT_PhotonHT::bookHistograms(DQMStore::IBooker &ibooker_, edm::Run const &, edm::EventSetup const &) {
edm::LogInfo("SUSY_HLT_PhotonHT") << "SUSY_HLT_PhotonHT::bookHistograms" << std::endl;
// book at beginRun
bookHistos(ibooker_);
}
void SUSY_HLT_PhotonHT::analyze(edm::Event const &e, edm::EventSetup const &eSetup) {
edm::LogInfo("SUSY_HLT_PhotonHT") << "SUSY_HLT_PhotonHT::analyze" << std::endl;
//-------------------------------
//--- MET / HT
//-------------------------------
edm::Handle<reco::PFMETCollection> pfMETCollection;
e.getByToken(thePfMETCollection_, pfMETCollection);
if (!pfMETCollection.isValid()) {
edm::LogError("SUSY_HLT_PhotonHT") << "invalid met collection"
<< "\n";
return;
}
//-------------------------------
//--- Photon
//-------------------------------
edm::Handle<reco::PhotonCollection> photonCollection;
e.getByToken(thePhotonCollection_, photonCollection);
if (!photonCollection.isValid()) {
edm::LogError("SUSY_HLT_PhotonHT") << "invalid egamma collection"
<< "\n";
return;
}
// check what is in the menu
edm::Handle<edm::TriggerResults> hltresults;
e.getByToken(triggerResults_, hltresults);
if (!hltresults.isValid()) {
edm::LogError("SUSY_HLT_PhotonHT") << "invalid collection: TriggerResults"
<< "\n";
return;
}
//-------------------------------
//--- Trigger
//-------------------------------
edm::Handle<trigger::TriggerEvent> triggerSummary;
e.getByToken(theTrigSummary_, triggerSummary);
if (!triggerSummary.isValid()) {
edm::LogError("SUSY_HLT_PhotonHT") << "invalid collection: TriggerSummary"
<< "\n";
return;
}
// get online objects
trigger::TriggerObjectCollection triggerObjects = triggerSummary->getObjects();
// get the photon object
size_t filterIndexPhoton = triggerSummary->filterIndex(triggerFilterPhoton_);
if (filterIndexPhoton < triggerSummary->sizeFilters()) {
const trigger::Keys &keys = triggerSummary->filterKeys(filterIndexPhoton);
if (!keys.empty()) {
// take the leading photon
float pt = triggerObjects[keys[0]].pt();
h_photonPt->Fill(pt);
}
}
// get ht
size_t filterIndexHt = triggerSummary->filterIndex(triggerFilterHt_);
if (filterIndexHt < triggerSummary->sizeFilters()) {
const trigger::Keys &keys = triggerSummary->filterKeys(filterIndexHt);
if (!keys.empty()) {
float ht = triggerObjects[keys[0]].pt();
h_ht->Fill(ht);
}
}
bool hasFired = false, hasFiredAuxiliaryForHadronicLeg = false;
const edm::TriggerNames &trigNames = e.triggerNames(*hltresults);
unsigned int numTriggers = trigNames.size();
for (unsigned int hltIndex = 0; hltIndex < numTriggers; ++hltIndex) {
if (trigNames.triggerName(hltIndex).find(triggerPath_) != std::string::npos && hltresults->wasrun(hltIndex) &&
hltresults->accept(hltIndex))
hasFired = true;
if (trigNames.triggerName(hltIndex).find(triggerPathAuxiliaryForHadronic_) != std::string::npos &&
hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex))
hasFiredAuxiliaryForHadronicLeg = true;
}
if (hasFiredAuxiliaryForHadronicLeg || !e.isRealData()) {
float recoPhotonPt = !photonCollection->empty() ? photonCollection->begin()->et() : 0;
float recoHt = !pfMETCollection->empty() ? pfMETCollection->begin()->et() : 0;
if (hasFired) {
if (!photonCollection->empty() && recoHt >= htThrOffline_)
h_photonTurnOn_num->Fill(recoPhotonPt);
if (!pfMETCollection->empty() && recoPhotonPt >= ptThrOffline_)
h_htTurnOn_num->Fill(recoHt);
}
if (!photonCollection->empty() && recoHt >= htThrOffline_)
h_photonTurnOn_den->Fill(recoPhotonPt);
if (!pfMETCollection->empty() && recoPhotonPt >= ptThrOffline_)
h_htTurnOn_den->Fill(recoPhotonPt);
}
}
void SUSY_HLT_PhotonHT::bookHistos(DQMStore::IBooker &ibooker_) {
ibooker_.cd();
ibooker_.setCurrentFolder("HLT/SUSYBSM/" + triggerPath_);
// offline quantities
h_photonPt = ibooker_.book1D("photonPt", "Photon transverse momentum; p_{T} (GeV)", 20, 0, 500);
h_ht = ibooker_.book1D("ht", "Hadronic activity;H_{T} (GeV)", 20, 0, 2000);
h_htTurnOn_num = ibooker_.book1D("pfHtTurnOn_num", "PF HT Turn On Numerator", 20, 300, 800);
h_htTurnOn_den = ibooker_.book1D("pfHtTurnOn_den", "PF HT Turn On Denominator", 20, 300, 800);
h_photonTurnOn_num = ibooker_.book1D("photonTurnOn_num", "Photon Turn On Numerator", 20, 70, 130);
h_photonTurnOn_den = ibooker_.book1D("photonTurnOn_den", "Photon Turn On Denominator", 20, 70, 130);
ibooker_.cd();
}
// define this as a plug-in
DEFINE_FWK_MODULE(SUSY_HLT_PhotonHT);
|