|
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
|
#include "HLTriggerOffline/Btag/interface/HLTVertexPerformanceAnalyzer.h"
using namespace edm;
using namespace reco;
HLTVertexPerformanceAnalyzer::HLTVertexPerformanceAnalyzer(const edm::ParameterSet &iConfig) {
mainFolder_ = iConfig.getParameter<std::string>("mainFolder");
hlTriggerResults_ = consumes<TriggerResults>(iConfig.getParameter<InputTag>("TriggerResults"));
VertexCollection_ =
edm::vector_transform(iConfig.getParameter<std::vector<edm::InputTag>>("Vertex"),
[this](edm::InputTag const &tag) { return mayConsume<reco::VertexCollection>(tag); });
hltPathNames_ = iConfig.getParameter<std::vector<std::string>>("HLTPathNames");
simVertexCollection_ = consumes<std::vector<SimVertex>>(iConfig.getParameter<edm::InputTag>("SimVertexCollection"));
EDConsumerBase::labelsForToken(hlTriggerResults_, label);
hlTriggerResults_Label = label.module;
for (unsigned int i = 0; i < VertexCollection_.size(); i++) {
EDConsumerBase::labelsForToken(VertexCollection_[i], label);
VertexCollection_Label.push_back(label.module);
}
}
HLTVertexPerformanceAnalyzer::~HLTVertexPerformanceAnalyzer() {}
void HLTVertexPerformanceAnalyzer::dqmBeginRun(const edm::Run &iRun, const edm::EventSetup &iSetup) {
triggerConfChanged_ = true;
EDConsumerBase::labelsForToken(hlTriggerResults_, label);
hltConfigProvider_.init(iRun, iSetup, label.process, triggerConfChanged_);
const std::vector<std::string> &allHltPathNames = hltConfigProvider_.triggerNames();
// fill hltPathIndexs_ with the trigger number of each hltPathNames_
for (size_t trgs = 0; trgs < hltPathNames_.size(); trgs++) {
unsigned int found = 1;
int it_mem = -1;
for (size_t it = 0; it < allHltPathNames.size(); ++it) {
found = allHltPathNames.at(it).find(hltPathNames_[trgs]);
if (found == 0) {
it_mem = (int)it;
}
} // for allallHltPathNames
hltPathIndexs_.push_back(it_mem);
} // for hltPathNames_
// fill _isfoundHLTs for each hltPathNames_
for (size_t trgs = 0; trgs < hltPathNames_.size(); trgs++) {
if (hltPathIndexs_[trgs] < 0) {
_isfoundHLTs.push_back(false);
} else {
_isfoundHLTs.push_back(true);
}
}
}
void HLTVertexPerformanceAnalyzer::analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup) {
bool trigRes = false;
using namespace edm;
// get triggerResults
Handle<TriggerResults> TriggerResulsHandler;
if (hlTriggerResults_Label.empty() || hlTriggerResults_Label == "NULL") {
edm::LogInfo("NoTriggerResults") << "TriggerResults ==> Empty";
return;
}
iEvent.getByToken(hlTriggerResults_, TriggerResulsHandler);
if (TriggerResulsHandler.isValid())
trigRes = true;
if (!trigRes) {
edm::LogInfo("NoTriggerResults") << "TriggerResults ==> not readable";
return;
}
const TriggerResults &triggerResults = *(TriggerResulsHandler.product());
// get simVertex
float simPV = 0;
Handle<std::vector<SimVertex>> simVertexCollection;
iEvent.getByToken(simVertexCollection_, simVertexCollection);
if (!simVertexCollection.isValid()) {
edm::LogInfo("NoSimVertex") << "SimVertex collection is invalid";
return;
}
if (simVertexCollection->empty()) {
edm::LogInfo("NoSimVertex") << "SimVertex collection is empty";
return;
}
const SimVertex simPVh = *(simVertexCollection->begin());
simPV = simPVh.position().z();
// fill the DQM plot
Handle<VertexCollection> VertexHandler;
for (unsigned int ind = 0; ind < hltPathNames_.size(); ind++) {
for (unsigned int coll = 0; coll < VertexCollection_.size(); coll++) {
bool VertexOK = false;
if (!_isfoundHLTs[ind])
continue; // if the hltPath is not in the event, skip the event
if (!triggerResults.accept(hltPathIndexs_[ind]))
continue; // if the hltPath was not accepted skip the event
// get the recoVertex
if (!VertexCollection_Label.at(coll).empty() && VertexCollection_Label.at(coll) != "NULL") {
iEvent.getByToken(VertexCollection_.at(coll), VertexHandler);
if (VertexHandler.isValid() > 0)
VertexOK = true;
}
if (VertexOK) {
// calculate the variable (RecoVertex - SimVertex)
float value = VertexHandler->begin()->z() - simPV;
// if value is over/under flow, assign the extreme value
float maxValue = H1_.at(ind)["Vertex_" + VertexCollection_Label.at(coll)]->getTH1F()->GetXaxis()->GetXmax();
if (value > maxValue)
value = maxValue - 0.0001;
if (value < -maxValue)
value = -maxValue + 0.0001;
// fill the histo
H1_.at(ind)["Vertex_" + VertexCollection_Label.at(coll)]->Fill(value);
}
} // for on VertexCollection_
} // for on hltPathNames_
}
void HLTVertexPerformanceAnalyzer::bookHistograms(DQMStore::IBooker &ibooker,
edm::Run const &iRun,
edm::EventSetup const &iSetup) {
// book the DQM plots
using namespace std;
std::string dqmFolder;
for (unsigned int ind = 0; ind < hltPathNames_.size(); ind++) {
dqmFolder = Form("%s/Vertex/%s", mainFolder_.c_str(), hltPathNames_[ind].c_str());
H1_.push_back(std::map<std::string, MonitorElement *>());
ibooker.setCurrentFolder(dqmFolder);
for (unsigned int coll = 0; coll < VertexCollection_.size(); coll++) {
float maxValue = 0.02;
if (VertexCollection_Label.at(coll) == ("hltFastPrimaryVertex"))
maxValue = 2.; // for the hltFastPrimaryVertex use a larger scale (res ~ 1 cm)
float vertexU = maxValue;
float vertexL = -maxValue;
int vertexBins = 100;
if (!VertexCollection_Label.at(coll).empty() && VertexCollection_Label.at(coll) != "NULL") {
H1_.back()["Vertex_" + VertexCollection_Label.at(coll)] =
ibooker.book1D("Vertex_" + VertexCollection_Label.at(coll),
VertexCollection_Label.at(coll).c_str(),
vertexBins,
vertexL,
vertexU);
H1_.back()["Vertex_" + VertexCollection_Label.at(coll)]->setAxisTitle("vertex error (cm)", 1);
}
}
}
}
DEFINE_FWK_MODULE(HLTVertexPerformanceAnalyzer);
|