|
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
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
|
// -*-C++-*-
#ifndef L1TGCT_H
#define L1TGCT_H
/*
* \file L1TGCT.h
*
* \author J. Berryhill
*
* Revision 1.23 2010/05/30 10:01:58 tapper
* Added one histogram, correlation of sum ET and HT and changed a few labels for the better.
*
* Revision 1.22 2009/11/19 14:33:13 puigh
* modify beginJob
*
* Revision 1.21 2009/11/02 17:00:04 tapper
* Changes to L1TdeGCT (to include energy sums), to L1TDEMON (should not make any difference now) and L1TGCT to add multiple BXs.
*
* Revision 1.20 2009/06/23 09:48:55 tapper
* Added missing occupancy plot for central and forward jets.
*
* Revision 1.19 2009/06/22 15:58:20 tapper
* Added MET vs MHT correlation plots (both for magnitude and phi). Still untested!
*
* Revision 1.18 2009/06/22 15:47:04 tapper
* Removed rank difference histograms and added MHT. Untested so far!
*
* Revision 1.17 2009/05/27 21:49:26 jad
* updated Total and Missing Energy histograms and added Overlow plots
*
* Revision 1.16 2008/11/11 13:20:31 tapper
* A whole list of house keeping:
* 1. New shifter histogram with central and forward jets together.
* 2. Relabelled Ring 0 and Ring 1 to Ring 1 and Ring 2 for HF rings.
* 3. Tidied up some histograms names to make all consistent.
* 4. Switched eta and phi in 2D plots to match RCT.
* 5. Removed 1D eta and phi plots. Will not be needed for Qtests in future.
*
* Revision 1.15 2008/09/21 14:33:12 jad
* updated HF Sums & Counts and added individual Jet Candidates and differences
*
* Revision 1.14 2008/06/09 11:08:05 tapper
* Removed electron sub-folders with histograms per eta and phi bin.
*
* Revision 1.13 2008/06/02 11:08:58 tapper
* Added HF ring histograms....
*
* Revision 1.12 2008/04/28 09:23:07 tapper
* Added 1D eta and phi histograms for electrons and jets as input to Q tests.
*
* Revision 1.11 2008/04/25 15:40:21 tapper
* Added histograms to EventInfo//errorSummarySegments.
*
* Revision 1.10 2008/03/01 00:40:00 lat
* DQM core migration.
*
* Revision 1.9 2008/02/20 19:24:24 tapper
* Removed noisy include.
*
* Revision 1.8 2008/02/20 18:59:29 tapper
* Ported GCTMonitor histograms into L1TGCT
*
* Revision 1.7 2007/09/04 02:54:21 wittich
* - fix dupe ME in RCT
* - put in rank>0 req in GCT
* - various small other fixes
*
* Revision 1.6 2007/08/31 18:14:20 wittich
* update GCT packages to reflect GctRawToDigi, and move to raw plots
*
* Revision 1.5 2007/08/31 11:02:55 wittich
* cerr -> LogInfo
*
* Revision 1.4 2007/02/22 19:43:52 berryhil
*
*
*
* InputTag parameters added for all modules
*
* Revision 1.3 2007/02/19 22:49:53 wittich
* - Add RCT monitor
*
* Revision 1.2 2007/02/19 21:11:23 wittich
* - Updates for integrating GCT monitor.
* + Adapted right now only the L1E elements thereof.
* + added DataFormats/L1Trigger to build file.
*
*
*/
// system include files
#include <memory>
#include <unistd.h>
#include <iostream>
#include <fstream>
#include <vector>
// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
// DQM
#include "DQMServices/Core/interface/DQMStore.h"
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctCollections.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
//
// class declaration
//
class L1TGCT : public DQMEDAnalyzer {
public:
// Constructor
L1TGCT(const edm::ParameterSet& ps);
// Destructor
~L1TGCT() override;
protected:
// Analyze
void analyze(const edm::Event& e, const edm::EventSetup& c) override;
void bookHistograms(DQMStore::IBooker& ibooker, edm::Run const&, edm::EventSetup const&) override;
private:
// ----------member data ---------------------------
std::string monitorDir_;
// trigger type information
MonitorElement* triggerType_;
// Jet and EM stuff
MonitorElement* l1GctAllJetsEtEtaPhi_;
MonitorElement* l1GctCenJetsEtEtaPhi_;
MonitorElement* l1GctForJetsEtEtaPhi_;
MonitorElement* l1GctTauJetsEtEtaPhi_;
MonitorElement* l1GctIsoTauJetsEtEtaPhi_;
MonitorElement* l1GctIsoEmRankEtaPhi_;
MonitorElement* l1GctNonIsoEmRankEtaPhi_;
MonitorElement* l1GctAllJetsOccEtaPhi_;
MonitorElement* l1GctCenJetsOccEtaPhi_;
MonitorElement* l1GctForJetsOccEtaPhi_;
MonitorElement* l1GctTauJetsOccEtaPhi_;
MonitorElement* l1GctIsoTauJetsOccEtaPhi_;
MonitorElement* l1GctIsoEmOccEtaPhi_;
MonitorElement* l1GctNonIsoEmOccEtaPhi_;
MonitorElement* l1GctCenJetsRank_;
MonitorElement* l1GctForJetsRank_;
MonitorElement* l1GctTauJetsRank_;
MonitorElement* l1GctIsoTauJetsRank_;
MonitorElement* l1GctIsoEmRank_;
MonitorElement* l1GctNonIsoEmRank_;
MonitorElement* l1GctAllJetsOccRankBx_;
MonitorElement* l1GctAllEmOccRankBx_;
// Energy sums stuff
MonitorElement* l1GctEtMiss_;
MonitorElement* l1GctEtMissPhi_;
MonitorElement* l1GctEtMissOf_;
MonitorElement* l1GctEtMissOccBx_;
MonitorElement* l1GctHtMiss_;
MonitorElement* l1GctHtMissPhi_;
MonitorElement* l1GctHtMissOf_;
MonitorElement* l1GctHtMissOccBx_;
MonitorElement* l1GctEtMissHtMissCorr_;
MonitorElement* l1GctEtMissHtMissCorrPhi_;
MonitorElement* l1GctEtTotal_;
MonitorElement* l1GctEtTotalOf_;
MonitorElement* l1GctEtTotalOccBx_;
MonitorElement* l1GctEtHad_;
MonitorElement* l1GctEtHadOf_;
MonitorElement* l1GctEtHadOccBx_;
MonitorElement* l1GctEtTotalEtHadCorr_;
// HF Rings stuff
MonitorElement* l1GctHFRing1PosEtaNegEta_;
MonitorElement* l1GctHFRing2PosEtaNegEta_;
MonitorElement* l1GctHFRing1TowerCountPosEtaNegEta_;
MonitorElement* l1GctHFRing2TowerCountPosEtaNegEta_;
MonitorElement* l1GctHFRing1TowerCountPosEta_;
MonitorElement* l1GctHFRing1TowerCountNegEta_;
MonitorElement* l1GctHFRing2TowerCountPosEta_;
MonitorElement* l1GctHFRing2TowerCountNegEta_;
MonitorElement* l1GctHFRing1ETSumPosEta_;
MonitorElement* l1GctHFRing1ETSumNegEta_;
MonitorElement* l1GctHFRing2ETSumPosEta_;
MonitorElement* l1GctHFRing2ETSumNegEta_;
MonitorElement* l1GctHFRingRatioPosEta_;
MonitorElement* l1GctHFRingRatioNegEta_;
MonitorElement* l1GctHFRingETSumOccBx_;
MonitorElement* l1GctHFRingTowerCountOccBx_;
int nev_; // Number of events processed
std::string outputFile_; //file name for ROOT ouput
bool verbose_;
bool monitorDaemon_;
std::ofstream logFile_;
edm::InputTag gctCenJetsSource_;
edm::InputTag gctForJetsSource_;
edm::InputTag gctTauJetsSource_;
edm::InputTag gctIsoTauJetsSource_;
edm::InputTag gctEnergySumsSource_;
edm::InputTag gctIsoEmSource_;
edm::InputTag gctNonIsoEmSource_;
/// filter TriggerType
bool m_stage1_layer2_;
int filterTriggerType_;
//define Token(-s)
edm::EDGetTokenT<L1GctEmCandCollection> gctIsoEmSourceToken_;
edm::EDGetTokenT<L1GctEmCandCollection> gctNonIsoEmSourceToken_;
edm::EDGetTokenT<L1GctJetCandCollection> gctCenJetsSourceToken_;
edm::EDGetTokenT<L1GctJetCandCollection> gctForJetsSourceToken_;
edm::EDGetTokenT<L1GctJetCandCollection> gctTauJetsSourceToken_;
edm::EDGetTokenT<L1GctJetCandCollection> gctIsoTauJetsSourceToken_;
edm::EDGetTokenT<L1GctHFRingEtSumsCollection> gctEnergySumsSourceToken_;
edm::EDGetTokenT<L1GctHFBitCountsCollection> l1HFCountsToken_;
edm::EDGetTokenT<L1GctEtMissCollection> l1EtMissToken_;
edm::EDGetTokenT<L1GctHtMissCollection> l1HtMissToken_;
edm::EDGetTokenT<L1GctEtHadCollection> l1EtHadToken_;
edm::EDGetTokenT<L1GctEtTotalCollection> l1EtTotalToken_;
};
#endif
|