|
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
|
/**
* \class HLTBeamModeFilter
*
*
* Description: see class header.
*
*
* \author: Vasile Mihai Ghete - HEPHY Vienna
*
*
*/
// this class header
#include "HLTBeamModeFilter.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
// system include files
#include <vector>
#include <iostream>
// user include files
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerEvmReadoutRecord.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/MessageLogger/interface/MessageDrop.h"
// constructor(s)
HLTBeamModeFilter::HLTBeamModeFilter(const edm::ParameterSet& parSet)
: HLTFilter(parSet),
m_l1GtEvmReadoutRecordTag(parSet.getParameter<edm::InputTag>("L1GtEvmReadoutRecordTag")),
m_l1GtEvmReadoutRecordToken(consumes<L1GlobalTriggerEvmReadoutRecord>(m_l1GtEvmReadoutRecordTag)),
m_allowedBeamMode(parSet.getParameter<std::vector<unsigned int> >("AllowedBeamMode")),
m_isDebugEnabled(edm::isDebugEnabled()) {
if (m_isDebugEnabled) {
LogDebug("HLTBeamModeFilter") << std::endl;
LogTrace("HLTBeamModeFilter") << "Input tag for L1 GT EVM record: " << m_l1GtEvmReadoutRecordTag
<< "\nAllowed beam modes:" << std::endl;
for (std::vector<unsigned int>::const_iterator itMode = m_allowedBeamMode.begin();
itMode != m_allowedBeamMode.end();
++itMode) {
LogTrace("HLTBeamModeFilter") << " " << (*itMode) << std::endl;
}
LogTrace("HLTBeamModeFilter") << std::endl;
}
}
// destructor
HLTBeamModeFilter::~HLTBeamModeFilter() {
// empty now
}
// member functions
void HLTBeamModeFilter::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
edm::ParameterSetDescription desc;
makeHLTFilterDescription(desc);
// #
// # InputTag for the L1 Global Trigger EVM readout record
// # gtDigis GT Emulator
// # l1GtEvmUnpack GT EVM Unpacker (default module name)
// # gtEvmDigis GT EVM Unpacker in RawToDigi standard sequence
// #
// # cloned GT unpacker in HLT = gtEvmDigis
desc.add<edm::InputTag>("L1GtEvmReadoutRecordTag", edm::InputTag("gtEvmDigis"));
// #
// # vector of allowed beam modes
// # default value: 11 (STABLE)
std::vector<unsigned int> allowedBeamMode(1, 11);
desc.add<std::vector<unsigned int> >("AllowedBeamMode", allowedBeamMode);
descriptions.add("hltBeamModeFilter", desc);
}
bool HLTBeamModeFilter::hltFilter(edm::Event& iEvent,
const edm::EventSetup& evSetup,
trigger::TriggerFilterObjectWithRefs& filterproduct) const {
// for MC samples, return always true (not even checking validity of L1GlobalTriggerEvmReadoutRecord)
// eventually, the BST information will be filled also in MC simulation to spare this check
if (!(iEvent.isRealData())) {
//LogDebug("HLTBeamModeFilter") << "\nRunning on MC sample."
// << "\nAll events are automatically selected, without testing the beam mode"
// << "\n" << std::endl;
return true;
}
// data
// initialize filter result
bool filterResult = false;
// get L1GlobalTriggerEvmReadoutRecord and beam mode
edm::Handle<L1GlobalTriggerEvmReadoutRecord> gtEvmReadoutRecord;
iEvent.getByToken(m_l1GtEvmReadoutRecordToken, gtEvmReadoutRecord);
if (!gtEvmReadoutRecord.isValid()) {
edm::LogWarning("HLTBeamModeFilter") << "\nWarning: L1GlobalTriggerEvmReadoutRecord with input tag "
<< m_l1GtEvmReadoutRecordTag
<< "\nrequested in configuration, but not found in the event." << std::endl;
return false;
}
const uint16_t beamModeValue = (gtEvmReadoutRecord->gtfeWord()).beamMode();
edm::LogInfo("HLTBeamModeFilter") << "Beam mode: " << beamModeValue;
for (unsigned int itMode : m_allowedBeamMode) {
if (beamModeValue == itMode) {
filterResult = true;
//LogTrace("HLTBeamModeFilter") << "Event selected - beam mode: "
// << beamModeValue << "\n" << std::endl;
break;
}
}
//
return filterResult;
}
// register as framework plugin
#include "FWCore/Framework/interface/MakerMacros.h"
#include <cstdint>
DEFINE_FWK_MODULE(HLTBeamModeFilter);
|