|
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
|
#include "DQMOffline/Trigger/interface/HLTTauDQMPathSummaryPlotter.h"
#include "DQMOffline/Trigger/interface/HLTTauDQMPath.h"
HLTTauDQMPathSummaryPlotter::HLTTauDQMPathSummaryPlotter(const edm::ParameterSet& pset,
bool doRefAnalysis,
const std::string& dqmBaseFolder,
double hltMatchDr)
: HLTTauDQMPlotter(pset, dqmBaseFolder), hltMatchDr_(hltMatchDr), doRefAnalysis_(doRefAnalysis) {}
HLTTauDQMPathSummaryPlotter::~HLTTauDQMPathSummaryPlotter() = default;
void HLTTauDQMPathSummaryPlotter::bookHistograms(HistoWrapper& iWrapper, DQMStore::IBooker& iBooker) {
if (!isValid() || pathObjects_.empty())
return;
//Create the histograms
iBooker.setCurrentFolder(triggerTag() + "/helpers");
all_events = iWrapper.book1D(iBooker, "RefEvents", "All events", pathObjects_.size(), 0, pathObjects_.size(), kVital);
accepted_events = iWrapper.book1D(iBooker,
"PathTriggerBits",
"Accepted Events per Path;;entries",
pathObjects_.size(),
0,
pathObjects_.size(),
kVital);
for (size_t i = 0; i < pathObjects_.size(); ++i) {
if (all_events)
all_events->setBinLabel(i + 1, pathObjects_[i]->getPathName());
if (accepted_events)
accepted_events->setBinLabel(i + 1, pathObjects_[i]->getPathName());
}
iBooker.setCurrentFolder(triggerTag());
}
void HLTTauDQMPathSummaryPlotter::analyze(const edm::TriggerResults& triggerResults,
const trigger::TriggerEvent& triggerEvent,
const HLTTauDQMOfflineObjects& refCollection) {
if (doRefAnalysis_) {
std::vector<HLTTauDQMPath::Object> triggerObjs;
std::vector<HLTTauDQMPath::Object> matchedTriggerObjs;
HLTTauDQMOfflineObjects matchedOfflineObjs;
for (size_t i = 0; i < pathObjects_.size(); ++i) {
const HLTTauDQMPath* path = pathObjects_[i];
const int lastFilter = path->filtersSize() - 1;
if (path->goodOfflineEvent(lastFilter, refCollection)) {
if (all_events)
all_events->Fill(i + 0.5);
}
if (path->fired(triggerResults)) {
triggerObjs.clear();
matchedTriggerObjs.clear();
matchedOfflineObjs.clear();
path->getFilterObjects(triggerEvent, lastFilter, triggerObjs);
if (path->offlineMatching(
lastFilter, triggerObjs, refCollection, hltMatchDr_, matchedTriggerObjs, matchedOfflineObjs)) {
if (accepted_events)
accepted_events->Fill(i + 0.5);
}
}
}
} else {
for (size_t i = 0; i < pathObjects_.size(); ++i) {
const HLTTauDQMPath* path = pathObjects_[i];
if (all_events)
all_events->Fill(i + 0.5);
if (path->fired(triggerResults)) {
if (accepted_events)
accepted_events->Fill(i + 0.5);
}
}
}
}
|