|
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
|
// -*- C++ -*-
//
// Package: Core
// Class : FWHLTTriggerTableView
//
// Implementation:
// [Notes on implementation]
//
// Original Author:
// Created: Tue Jan 25 16:02:03 CET 2011
//
#include <boost/regex.hpp>
#include "Fireworks/Core/interface/FWHLTTriggerTableView.h"
#include "Fireworks/Core/interface/fwLog.h"
#include "DataFormats/FWLite/interface/Handle.h"
#include "DataFormats/FWLite/interface/Event.h"
#include "FWCore/Common/interface/TriggerNames.h"
#include "DataFormats/Common/interface/TriggerResults.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "Fireworks/Core/interface/FWTriggerTableViewTableManager.h"
#include "Fireworks/Core/interface/CmsShowViewPopup.h"
FWHLTTriggerTableView::FWHLTTriggerTableView(TEveWindowSlot* iParent)
: FWTriggerTableView(iParent, FWViewType::kTableHLT), m_event(nullptr) {
m_columns[0].title = "Filter Name";
m_columns.push_back(Column("Accept"));
m_columns.push_back(Column("Average Accept"));
dataChanged();
}
void FWHLTTriggerTableView::fillTable(fwlite::Event* event) {
if (event != m_event) {
m_event = event;
fillAverageAcceptFractions();
}
fwlite::Handle<edm::TriggerResults> hTriggerResults;
edm::TriggerNames const* triggerNames(nullptr);
try {
hTriggerResults.getByLabel(*event, "TriggerResults", "", m_process.value().c_str());
triggerNames = &event->triggerNames(*hTriggerResults);
} catch (cms::Exception&) {
fwLog(fwlog::kWarning) << " no trigger results with process name HLT is available" << std::endl;
m_tableManager->dataChanged();
return;
}
boost::regex filter(m_regex.value());
for (unsigned int i = 0; i < triggerNames->size(); ++i) {
if (!boost::regex_search(triggerNames->triggerName(i), filter))
continue;
m_columns.at(0).values.push_back(triggerNames->triggerName(i));
m_columns.at(1).values.push_back(Form("%d", hTriggerResults->accept(i)));
m_columns.at(2).values.push_back(Form("%6.1f%%", m_averageAccept[triggerNames->triggerName(i)] * 100));
}
}
void FWHLTTriggerTableView::fillAverageAcceptFractions() {
edm::EventID currentEvent = m_event->id();
// better to keep the keys and just set to zero the values
for (acceptmap_t::iterator it = m_averageAccept.begin(), ed = m_averageAccept.end(); it != ed; ++it) {
it->second = 0;
}
// loop over events
fwlite::Handle<edm::TriggerResults> hTriggerResults;
for (m_event->toBegin(); !m_event->atEnd(); ++(*m_event)) {
hTriggerResults.getByLabel(*m_event, "TriggerResults", "", "HLT");
edm::TriggerNames const* triggerNames(nullptr);
try {
triggerNames = &m_event->triggerNames(*hTriggerResults);
} catch (cms::Exception&) {
fwLog(fwlog::kError) << " exception caught while trying to get trigger info" << std::endl;
break;
}
for (unsigned int i = 0; i < triggerNames->size(); ++i) {
if (hTriggerResults->accept(i)) {
m_averageAccept[triggerNames->triggerName(i)]++;
}
}
}
m_event->to(currentEvent);
double denominator = 1.0 / m_event->size();
for (acceptmap_t::iterator it = m_averageAccept.begin(), ed = m_averageAccept.end(); it != ed; ++it) {
it->second *= denominator;
}
}
|