|
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
|
#include "DQM/EcalMonitorTasks/interface/IntegrityTask.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/Exception.h"
namespace ecaldqm {
IntegrityTask::IntegrityTask() : DQWorkerTask() {}
void IntegrityTask::beginEvent(edm::Event const& _evt,
edm::EventSetup const& _es,
bool const& ByLumiResetSwitch,
bool&) {
if (ByLumiResetSwitch) {
MEs_.at("MapByLumi").reset(GetElectronicsMap());
MEs_.at("ByLumi").reset(GetElectronicsMap());
MEs_.at("TTIDByLumi").reset(GetElectronicsMap());
}
}
template <typename IDCollection>
void IntegrityTask::runOnDetIdCollection(IDCollection const& _ids, Collections _collection) {
// Collection is empty if there are no errors
if (_ids.empty())
return;
MESet* set(nullptr);
switch (_collection) {
case kEBGainErrors:
case kEEGainErrors:
set = &MEs_.at("Gain");
break;
case kEBChIdErrors:
case kEEChIdErrors:
set = &MEs_.at("ChId");
break;
case kEBGainSwitchErrors:
case kEEGainSwitchErrors:
set = &MEs_.at("GainSwitch");
break;
default:
return;
}
MESet& meMapByLumi(MEs_.at("MapByLumi"));
MESet& meByLumi(MEs_.at("ByLumi"));
MESet& meTotal(MEs_.at("Total"));
MESet& meTrendNErrors(MEs_.at("TrendNErrors"));
std::for_each(_ids.begin(), _ids.end(), [&](typename IDCollection::value_type const& id) {
set->fill(getEcalDQMSetupObjects(), id);
int dccid(dccId(id, GetElectronicsMap()));
meByLumi.fill(getEcalDQMSetupObjects(), dccid);
meTotal.fill(getEcalDQMSetupObjects(), dccid);
// Fill Integrity Errors Map with channel errors for this lumi
meMapByLumi.fill(getEcalDQMSetupObjects(), id);
meTrendNErrors.fill(getEcalDQMSetupObjects(), double(timestamp_.iLumi), 1.);
});
}
void IntegrityTask::runOnElectronicsIdCollection(EcalElectronicsIdCollection const& _ids, Collections _collection) {
// Collection is empty if there are no errors
if (_ids.empty())
return;
MESet* set(nullptr);
switch (_collection) {
case kTowerIdErrors:
set = &MEs_.at("TowerId");
break;
case kBlockSizeErrors:
set = &MEs_.at("BlockSize");
break;
default:
return;
}
MESet& meMapByLumi(MEs_.at("MapByLumi"));
MESet& meByLumi(MEs_.at("ByLumi"));
MESet& meTotal(MEs_.at("Total"));
MESet& meTrendNErrors(MEs_.at("TrendNErrors"));
MESet& meTTIDTotal(MEs_.at("TTIDTotal"));
MESet& meTTIDByLumi(MEs_.at("TTIDByLumi"));
std::for_each(_ids.begin(), _ids.end(), [&](EcalElectronicsIdCollection::value_type const& id) {
set->fill(getEcalDQMSetupObjects(), id);
int dccid(id.dccId());
double nCrystals(0.);
std::vector<DetId> chIds(GetElectronicsMap()->dccTowerConstituents(dccid, id.towerId()));
if (dccid <= kEEmHigh + 1 || dccid >= kEEpLow + 1)
nCrystals = chIds.size();
else
nCrystals = 25.;
meByLumi.fill(getEcalDQMSetupObjects(), dccid, nCrystals);
meTotal.fill(getEcalDQMSetupObjects(), dccid, nCrystals);
if (_collection == kTowerIdErrors) {
meTTIDByLumi.fill(getEcalDQMSetupObjects(), dccid, nCrystals);
meTTIDTotal.fill(getEcalDQMSetupObjects(), dccid, nCrystals);
}
// Fill Integrity Errors Map with tower errors for this lumi
// Since binned by crystal for compatibility with channel errors,
// fill with constituent channels of tower
for (std::vector<DetId>::iterator chItr(chIds.begin()); chItr != chIds.end(); ++chItr)
meMapByLumi.fill(getEcalDQMSetupObjects(), *chItr);
meTrendNErrors.fill(getEcalDQMSetupObjects(), double(timestamp_.iLumi), nCrystals);
});
}
DEFINE_ECALDQM_WORKER(IntegrityTask);
} // namespace ecaldqm
|