|
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
|
#include "DataFormats/HLTReco/interface/TriggerObject.h"
#include "DataFormats/Math/interface/LorentzVector.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_Muon_BJet.h"
SUSY_HLT_Muon_BJet::SUSY_HLT_Muon_BJet(const edm::ParameterSet &ps) {
edm::LogInfo("SUSY_HLT_Muon_BJet") << "Constructor SUSY_HLT_Muon_BJet::SUSY_HLT_Muon_BJet " << std::endl;
// Get parameters from configuration file
theTrigSummary_ = consumes<trigger::TriggerEvent>(ps.getParameter<edm::InputTag>("trigSummary"));
theMuonCollection_ = consumes<reco::MuonCollection>(ps.getParameter<edm::InputTag>("MuonCollection"));
thePfJetCollection_ = consumes<reco::PFJetCollection>(ps.getParameter<edm::InputTag>("pfJetCollection"));
theCaloJetCollection_ = consumes<reco::CaloJetCollection>(ps.getParameter<edm::InputTag>("caloJetCollection"));
triggerResults_ = consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>("TriggerResults"));
HLTProcess_ = ps.getParameter<std::string>("HLTProcess");
triggerPath_ = ps.getParameter<std::string>("TriggerPath");
triggerFilterMuon_ = ps.getParameter<edm::InputTag>("TriggerFilterMuon");
triggerFilterJet_ = ps.getParameter<edm::InputTag>("TriggerFilterJet");
ptThrJet_ = ps.getUntrackedParameter<double>("PtThrJet");
etaThrJet_ = ps.getUntrackedParameter<double>("EtaThrJet");
}
SUSY_HLT_Muon_BJet::~SUSY_HLT_Muon_BJet() {
edm::LogInfo("SUSY_HLT_Muon_BJet") << "Destructor SUSY_HLT_Muon_BJet::~SUSY_HLT_Muon_BJet " << std::endl;
}
void SUSY_HLT_Muon_BJet::dqmBeginRun(edm::Run const &run,
edm::EventSetup const &e) //
{
bool changed;
if (!fHltConfig.init(run, e, HLTProcess_, changed)) {
edm::LogError("SUSY_HLT_Muon_BJet") << "Initialization of HLTConfigProvider failed!!";
return;
}
bool pathFound = false;
const std::vector<std::string> allTrigNames = fHltConfig.triggerNames();
for (size_t j = 0; j < allTrigNames.size(); ++j) {
if (allTrigNames[j].find(triggerPath_) != std::string::npos) {
pathFound = true;
}
}
if (!pathFound) {
LogDebug("SUSY_HLT_Muon_BJet") << "Path not found"
<< "\n";
return;
}
edm::LogInfo("SUSY_HLT_Muon_BJet") << "SUSY_HLT_Muon_BJet::beginRun" << std::endl;
}
void SUSY_HLT_Muon_BJet::bookHistograms(DQMStore::IBooker &ibooker_, edm::Run const &, edm::EventSetup const &) {
edm::LogInfo("SUSY_HLT_Muon_BJet") << "SUSY_HLT_Muon_BJet::bookHistograms" << std::endl;
// book at beginRun
bookHistos(ibooker_);
}
void SUSY_HLT_Muon_BJet::analyze(edm::Event const &e, edm::EventSetup const &eSetup) {
edm::LogInfo("SUSY_HLT_Muon_BJet") << "SUSY_HLT_Muon_BJet::analyze" << std::endl;
//-------------------------------
//--- Trigger
//-------------------------------
edm::Handle<edm::TriggerResults> hltresults;
e.getByToken(triggerResults_, hltresults);
if (!hltresults.isValid()) {
edm::LogError("SUSY_HLT_Muon_BJet") << "invalid collection: TriggerResults"
<< "\n";
return;
}
edm::Handle<trigger::TriggerEvent> triggerSummary;
e.getByToken(theTrigSummary_, triggerSummary);
if (!triggerSummary.isValid()) {
edm::LogError("SUSY_HLT_Muon_BJet") << "invalid collection: TriggerSummary"
<< "\n";
return;
}
// get online objects
// std::vector<float> ptMuon, etaMuon, phiMuon;
size_t filterIndex = triggerSummary->filterIndex(triggerFilterMuon_);
trigger::TriggerObjectCollection triggerObjects = triggerSummary->getObjects();
if (!(filterIndex >= triggerSummary->sizeFilters())) {
const trigger::Keys &keys = triggerSummary->filterKeys(filterIndex);
for (size_t j = 0; j < keys.size(); ++j) {
trigger::TriggerObject foundObject = triggerObjects[keys[j]];
if (fabs(foundObject.id()) == 13) { // It's a muon
h_triggerMuPt->Fill(foundObject.pt());
h_triggerMuEta->Fill(foundObject.eta());
h_triggerMuPhi->Fill(foundObject.phi());
// ptMuon.push_back(foundObject.pt());
// etaMuon.push_back(foundObject.eta());
// phiMuon.push_back(foundObject.phi());
}
}
}
size_t filterIndex2 = triggerSummary->filterIndex(triggerFilterJet_);
if (!(filterIndex2 >= triggerSummary->sizeFilters())) {
const trigger::Keys &keys = triggerSummary->filterKeys(filterIndex2);
for (size_t j = 0; j < keys.size(); ++j) {
trigger::TriggerObject foundObject = triggerObjects[keys[j]];
h_triggerJetPt->Fill(foundObject.pt());
h_triggerJetEta->Fill(foundObject.eta());
h_triggerJetPhi->Fill(foundObject.phi());
}
}
}
void SUSY_HLT_Muon_BJet::bookHistos(DQMStore::IBooker &ibooker_) {
ibooker_.cd();
ibooker_.setCurrentFolder("HLT/SUSYBSM/" + triggerPath_);
// offline quantities
// online quantities
h_triggerMuPt = ibooker_.book1D("triggerMuPt", "Trigger Muon Pt; GeV", 50, 0.0, 500.0);
h_triggerMuEta = ibooker_.book1D("triggerMuEta", "Trigger Muon Eta", 20, -3.0, 3.0);
h_triggerMuPhi = ibooker_.book1D("triggerMuPhi", "Trigger Muon Phi", 20, -3.5, 3.5);
h_triggerJetPt = ibooker_.book1D("triggerJetPt", "Trigger Jet Pt; GeV", 50, 0.0, 500.0);
h_triggerJetEta = ibooker_.book1D("triggerJetEta", "Trigger Jet Eta", 20, -3.0, 3.0);
h_triggerJetPhi = ibooker_.book1D("triggerJetPhi", "Trigger Jet Phi", 20, -3.5, 3.5);
ibooker_.cd();
}
// define this as a plug-in
DEFINE_FWK_MODULE(SUSY_HLT_Muon_BJet);
|