|
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
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
|
#include "DQM/L1TMonitorClient/interface/L1TStage2CaloLayer2DEClientSummary.h"
L1TStage2CaloLayer2DEClientSummary::L1TStage2CaloLayer2DEClientSummary(const edm::ParameterSet &ps)
: monitor_dir_(ps.getUntrackedParameter<std::string>("monitorDir", "")),
hlSummary(nullptr),
jetSummary(nullptr),
egSummary(nullptr),
tauSummary(nullptr),
sumSummary(nullptr) {}
L1TStage2CaloLayer2DEClientSummary::~L1TStage2CaloLayer2DEClientSummary() {}
void L1TStage2CaloLayer2DEClientSummary::dqmEndLuminosityBlock(DQMStore::IBooker &ibooker,
DQMStore::IGetter &igetter,
const edm::LuminosityBlock &lumiSeg,
const edm::EventSetup &c) {
book(ibooker);
processHistograms(igetter);
}
void L1TStage2CaloLayer2DEClientSummary::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {
book(ibooker);
processHistograms(igetter);
}
void L1TStage2CaloLayer2DEClientSummary::book(DQMStore::IBooker &ibooker) {
ibooker.setCurrentFolder(monitor_dir_);
if (hlSummary == nullptr) {
hlSummary = ibooker.book1D("High level summary", "Event by event comparison summary", 5, 1, 6);
hlSummary->setBinLabel(1, "good events");
hlSummary->setBinLabel(2, "good jets");
hlSummary->setBinLabel(3, "good e/gs");
hlSummary->setBinLabel(4, "good taus");
hlSummary->setBinLabel(5, "good sums");
} else {
hlSummary->Reset();
}
if (jetSummary == nullptr) {
jetSummary = ibooker.book1D("Jet Agreement Summary", "Jet Agreement Summary", 4, 1, 5);
jetSummary->setBinLabel(1, "good jets");
jetSummary->setBinLabel(2, "jets pos off only");
jetSummary->setBinLabel(3, "jets Et off only ");
jetSummary->setBinLabel(4, "jets qual off only ");
} else {
jetSummary->Reset();
}
if (egSummary == nullptr) {
egSummary = ibooker.book1D("EG Agreement Summary", "EG Agreement Summary", 7, 1, 8);
egSummary->setBinLabel(1, "good non-iso e/gs");
egSummary->setBinLabel(2, "non-iso e/gs pos off");
egSummary->setBinLabel(3, "non-iso e/gs Et off");
egSummary->setBinLabel(4, "good iso e/gs");
egSummary->setBinLabel(5, "iso e/gs pos off");
egSummary->setBinLabel(6, "iso e/gs Et off");
egSummary->setBinLabel(7, "e/gs Iso off");
} else {
egSummary->Reset();
}
if (tauSummary == nullptr) {
tauSummary = ibooker.book1D("Tau Agreement Summary", "Tau Agremeent Summary", 7, 1, 8);
tauSummary->setBinLabel(1, "good non-iso taus");
tauSummary->setBinLabel(2, "non-iso taus pos off");
tauSummary->setBinLabel(3, "non-iso taus Et off");
tauSummary->setBinLabel(4, "good iso taus");
tauSummary->setBinLabel(5, "iso taus pos off");
tauSummary->setBinLabel(6, "iso taus Et off");
tauSummary->setBinLabel(7, "taus Iso off");
} else {
tauSummary->Reset();
}
if (sumSummary == nullptr) {
sumSummary = ibooker.book1D("Energy Sum Agreement Summary", "Sum Agreement Summary", 9, 1, 10);
sumSummary->setBinLabel(1, "good sums");
sumSummary->setBinLabel(2, "good ETT sums");
sumSummary->setBinLabel(3, "good HTT sums");
sumSummary->setBinLabel(4, "good MET sums");
sumSummary->setBinLabel(5, "good MHT sums");
sumSummary->setBinLabel(6, "good MBHF sums");
sumSummary->setBinLabel(7, "good TowCount sums");
sumSummary->setBinLabel(8, "good AsymCount sums");
sumSummary->setBinLabel(9, "good CentrCount sums");
} else {
sumSummary->Reset();
}
}
void L1TStage2CaloLayer2DEClientSummary::processHistograms(DQMStore::IGetter &igetter) {
// get reference to relevant summary MonitorElement instances
// - high level summary
// - eg agreement summary
// - energy sum agreement summary
// - jet agreement summary
// - tau agreement summary
MonitorElement *hlSummary_ = igetter.get(monitor_dir_ + "/expert/CaloL2 Object Agreement Summary");
MonitorElement *jetSummary_ = igetter.get(monitor_dir_ + "/expert/Jet Agreement Summary");
MonitorElement *egSummary_ = igetter.get(monitor_dir_ + "/expert/EG Agreement Summary");
MonitorElement *tauSummary_ = igetter.get(monitor_dir_ + "/expert/Tau Agreement Summary");
MonitorElement *sumSummary_ = igetter.get(monitor_dir_ + "/expert/Energy Sum Agreement Summary");
// check for existance of object
if (hlSummary_) {
// by default show 0% agreement (for edge case when no objects are found)
double evtRatio = 0, jetRatio = 0, egRatio = 0, tauRatio = 0, sumRatio = 0;
double totalEvents = hlSummary_->getBinContent(1);
double goodEvents = hlSummary_->getBinContent(2);
double totalJets = hlSummary_->getBinContent(3);
double goodJets = hlSummary_->getBinContent(4);
double totalEg = hlSummary_->getBinContent(5);
double goodEg = hlSummary_->getBinContent(6);
double totalTau = hlSummary_->getBinContent(7);
double goodTau = hlSummary_->getBinContent(8);
double totalSums = hlSummary_->getBinContent(9);
double goodSums = hlSummary_->getBinContent(10);
if (totalEvents != 0)
evtRatio = goodEvents / totalEvents;
if (totalJets != 0)
jetRatio = goodJets / totalJets;
if (totalEg != 0)
egRatio = goodEg / totalEg;
if (totalTau != 0)
tauRatio = goodTau / totalTau;
if (totalSums != 0)
sumRatio = goodSums / totalSums;
hlSummary->setBinContent(1, evtRatio);
hlSummary->setBinContent(2, jetRatio);
hlSummary->setBinContent(3, egRatio);
hlSummary->setBinContent(4, tauRatio);
hlSummary->setBinContent(5, sumRatio);
}
if (jetSummary_) {
// by default show 0% agreement (for edge case when no objects are found)
double goodRatio = 0, posOffRatio = 0, etOffRatio = 0, qualOffRatio = 0;
double totalJets = jetSummary_->getBinContent(1);
double goodJets = jetSummary_->getBinContent(2);
double jetPosOff = jetSummary_->getBinContent(3);
double jetEtOff = jetSummary_->getBinContent(4);
double jetQualOff = jetSummary_->getBinContent(5);
if (totalJets != 0) {
goodRatio = goodJets / totalJets;
posOffRatio = jetPosOff / totalJets;
etOffRatio = jetEtOff / totalJets;
qualOffRatio = jetQualOff / totalJets;
}
jetSummary->setBinContent(1, goodRatio);
jetSummary->setBinContent(2, posOffRatio);
jetSummary->setBinContent(3, etOffRatio);
jetSummary->setBinContent(4, qualOffRatio);
}
if (egSummary_) {
// by default show 0% agreement (for edge case when no objects are found)
double goodEgRatio = 0, egPosOffRatio = 0, egEtOffRatio = 0, goodIsoEgRatio = 0, isoEgPosOffRatio = 0,
isoEgEtOffRatio = 0, egIsoOffRatio = 0;
double totalEgs = egSummary_->getBinContent(1);
double goodEgs = egSummary_->getBinContent(2);
double egPosOff = egSummary_->getBinContent(3);
double egEtOff = egSummary_->getBinContent(4);
double totalIsoEgs = egSummary_->getBinContent(5);
double goodIsoEgs = egSummary_->getBinContent(6);
double isoEgPosOff = egSummary_->getBinContent(7);
double isoEgEtOff = egSummary_->getBinContent(8);
double egIsoOff = egSummary_->getBinContent(9);
if (totalEgs != 0) {
goodEgRatio = goodEgs / totalEgs;
egPosOffRatio = egPosOff / totalEgs;
egEtOffRatio = egEtOff / totalEgs;
}
if (totalIsoEgs != 0) {
goodIsoEgRatio = goodIsoEgs / totalIsoEgs;
isoEgPosOffRatio = isoEgPosOff / totalIsoEgs;
isoEgEtOffRatio = isoEgEtOff / totalIsoEgs;
}
if ((totalEgs + totalIsoEgs) > 0)
egIsoOffRatio = egIsoOff / (totalEgs + totalIsoEgs);
egSummary->setBinContent(1, goodEgRatio);
egSummary->setBinContent(2, egPosOffRatio);
egSummary->setBinContent(3, egEtOffRatio);
egSummary->setBinContent(4, goodIsoEgRatio);
egSummary->setBinContent(5, isoEgPosOffRatio);
egSummary->setBinContent(6, isoEgEtOffRatio);
egSummary->setBinContent(7, egIsoOffRatio);
}
if (tauSummary_) {
// by default show 0% agreement (for edge case when no objects are found)
double goodTauRatio = 0, tauPosOffRatio = 0, tauEtOffRatio = 0, goodIsoTauRatio = 0, isoTauPosOffRatio = 0,
isoTauEtOffRatio = 0, tauIsoOffRatio = 0;
double totalTaus = tauSummary_->getBinContent(1);
double goodTaus = tauSummary_->getBinContent(2);
double tauPosOff = tauSummary_->getBinContent(3);
double tauEtOff = tauSummary_->getBinContent(4);
double totalIsoTaus = tauSummary_->getBinContent(5);
double goodIsoTaus = tauSummary_->getBinContent(6);
double isoTauPosOff = tauSummary_->getBinContent(7);
double isoTauEtOff = tauSummary_->getBinContent(8);
double tauIsoOff = tauSummary_->getBinContent(9);
if (totalTaus != 0) {
goodTauRatio = goodTaus / totalTaus;
tauPosOffRatio = tauPosOff / totalTaus;
tauEtOffRatio = tauEtOff / totalTaus;
}
if (totalIsoTaus != 0) {
goodIsoTauRatio = goodIsoTaus / totalIsoTaus;
isoTauPosOffRatio = isoTauPosOff / totalIsoTaus;
isoTauEtOffRatio = isoTauEtOff / totalIsoTaus;
}
if ((totalTaus + totalIsoTaus) > 0)
tauIsoOffRatio = tauIsoOff / (totalTaus + totalIsoTaus);
tauSummary->setBinContent(1, goodTauRatio);
tauSummary->setBinContent(2, tauPosOffRatio);
tauSummary->setBinContent(3, tauEtOffRatio);
tauSummary->setBinContent(4, goodIsoTauRatio);
tauSummary->setBinContent(5, isoTauPosOffRatio);
tauSummary->setBinContent(6, isoTauEtOffRatio);
tauSummary->setBinContent(7, tauIsoOffRatio);
}
if (sumSummary_) {
// by default show 0% agreement (for edge case when no objects are found)
double goodSumRatio = 0, goodETTRatio = 0, goodHTTRatio = 0, goodMETRatio = 0, goodMHTRatio = 0, goodMBHFRatio = 0,
goodTowCountRatio = 0, goodAsymCountRatio = 0, goodCentrCountRatio = 0;
double totalSums = sumSummary_->getBinContent(1);
double goodSums = sumSummary_->getBinContent(2);
double totalETT = sumSummary_->getBinContent(3);
double goodETT = sumSummary_->getBinContent(4);
double totalHTT = sumSummary_->getBinContent(5);
double goodHTT = sumSummary_->getBinContent(6);
double totalMET = sumSummary_->getBinContent(7);
double goodMET = sumSummary_->getBinContent(8);
double totalMHT = sumSummary_->getBinContent(9);
double goodMHT = sumSummary_->getBinContent(10);
double totalMBHF = sumSummary_->getBinContent(11);
double goodMBHF = sumSummary_->getBinContent(12);
double totalTowCount = sumSummary_->getBinContent(13);
double goodTowCount = sumSummary_->getBinContent(14);
double totalAsymCount = sumSummary_->getBinContent(15);
double goodAsymCount = sumSummary_->getBinContent(16);
double totalCentrCount = sumSummary_->getBinContent(17);
double goodCentrCount = sumSummary_->getBinContent(18);
if (totalSums)
goodSumRatio = goodSums / totalSums;
if (totalETT)
goodETTRatio = goodETT / totalETT;
if (totalHTT)
goodHTTRatio = goodHTT / totalHTT;
if (totalMET)
goodMETRatio = goodMET / totalMET;
if (totalMHT)
goodMHTRatio = goodMHT / totalMHT;
if (totalMBHF)
goodMBHFRatio = goodMBHF / totalMBHF;
if (totalTowCount)
goodTowCountRatio = goodTowCount / totalTowCount;
if (totalAsymCount)
goodAsymCountRatio = goodAsymCount / totalAsymCount;
if (totalCentrCount)
goodCentrCountRatio = goodCentrCount / totalCentrCount;
sumSummary->setBinContent(1, goodSumRatio);
sumSummary->setBinContent(2, goodETTRatio);
sumSummary->setBinContent(3, goodHTTRatio);
sumSummary->setBinContent(4, goodMETRatio);
sumSummary->setBinContent(5, goodMHTRatio);
sumSummary->setBinContent(6, goodMBHFRatio);
sumSummary->setBinContent(7, goodTowCountRatio);
sumSummary->setBinContent(8, goodAsymCountRatio);
sumSummary->setBinContent(9, goodCentrCountRatio);
}
}
|