|
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
|
#include "DQM/EcalMonitorTasks/interface/PNDiodeTask.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
namespace ecaldqm {
PNDiodeTask::PNDiodeTask() : DQWorkerTask() { std::fill_n(enable_, nDCC, false); }
bool PNDiodeTask::filterRunType(short const* _runType) {
bool enable(false);
for (int iDCC(0); iDCC < 54; iDCC++) {
if (_runType[iDCC] == EcalDCCHeaderBlock::LASER_STD || _runType[iDCC] == EcalDCCHeaderBlock::LASER_GAP ||
_runType[iDCC] == EcalDCCHeaderBlock::LED_STD || _runType[iDCC] == EcalDCCHeaderBlock::LED_GAP ||
_runType[iDCC] == EcalDCCHeaderBlock::TESTPULSE_MGPA || _runType[iDCC] == EcalDCCHeaderBlock::TESTPULSE_GAP ||
_runType[iDCC] == EcalDCCHeaderBlock::PEDESTAL_STD || _runType[iDCC] == EcalDCCHeaderBlock::PEDESTAL_GAP) {
enable = true;
enable_[iDCC] = true;
} else
enable_[iDCC] = false;
}
return enable;
}
void PNDiodeTask::runOnErrors(EcalElectronicsIdCollection const& _ids, Collections _collection) {
if (_ids.empty())
return;
MESet* meMEMErrors = &MEs_.at("MEMErrors");
// MEM Box Integrity Errors (TowerIds 69 and 70)
// errorType matches to the following labels in DQM plot
// 0 = TOWERID
// 1 = BLOCKSIZE
// 2 = CHID
// 3 = GAIN
int errorType(-1);
switch (_collection) {
case kMEMTowerIdErrors:
errorType = 0;
break;
case kMEMBlockSizeErrors:
errorType = 1;
break;
case kMEMChIdErrors:
errorType = 2;
break;
case kMEMGainErrors:
errorType = 3;
break;
default:
return;
}
// Integrity errors for MEM boxes (towerIds 69/70)
// Plot contains two bins per dccId. Integer number
// bins correspond to towerId 69 and half integer
// number bins correspond to towerId 70.
std::for_each(_ids.begin(),
_ids.end(),
[&](EcalElectronicsIdCollection::value_type const& id) {
if (id.towerId() == 69)
meMEMErrors->fill(getEcalDQMSetupObjects(), id.dccId() + 0.0, errorType);
else if (id.towerId() == 70)
meMEMErrors->fill(getEcalDQMSetupObjects(), id.dccId() + 0.5, errorType);
else {
edm::LogWarning("EcalDQM")
<< "PNDiodeTask::runOnErrors : one of the ids in the electronics ID collection does not "
<< "correspond to one of the MEM box towerIds (69/70) in lumi number " << timestamp_.iLumi
<< ", event number " << timestamp_.iEvt;
}
});
}
void PNDiodeTask::runOnPnDigis(EcalPnDiodeDigiCollection const& _digis) {
MESet& meOccupancy(MEs_.at("Occupancy"));
MESet& meOccupancySummary(MEs_.at("OccupancySummary"));
MESet& mePedestal(MEs_.at("Pedestal"));
std::for_each(_digis.begin(), _digis.end(), [&](EcalPnDiodeDigiCollection::value_type const& digi) {
const EcalPnDiodeDetId& id(digi.id());
if (!enable_[dccId(id, GetElectronicsMap()) - 1])
return;
meOccupancy.fill(getEcalDQMSetupObjects(), id);
meOccupancySummary.fill(getEcalDQMSetupObjects(), id);
for (int iSample(0); iSample < 4; iSample++) {
if (digi.sample(iSample).gainId() != 1)
break;
mePedestal.fill(getEcalDQMSetupObjects(), id, double(digi.sample(iSample).adc()));
}
});
}
DEFINE_ECALDQM_WORKER(PNDiodeTask);
} // namespace ecaldqm
|