|
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
|
#include "DQM/HcalTasks/interface/FCDTask.h"
bool operator==(const FCDTask::FCDChannel& lhs, const FCDTask::FCDChannel& rhs) {
return ((lhs.crate == rhs.crate) && (lhs.slot == rhs.slot) && (lhs.fiber == rhs.fiber) &&
(lhs.fiberChannel == rhs.fiberChannel));
}
FCDTask::FCDTask(edm::ParameterSet const& ps)
: hcalDbServiceToken_(esConsumes<HcalDbService, HcalDbRecord, edm::Transition::BeginRun>()) {
// tags
_tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagQIE10", edm::InputTag("hcalDigis", "ZDC"));
_tokQIE10 = consumes<QIE10DigiCollection>(_tagQIE10);
// channels
edm::ParameterSet channelPSet = ps.getParameter<edm::ParameterSet>("fcdChannels");
std::vector<int32_t> crates = channelPSet.getUntrackedParameter<std::vector<int32_t> >("crate");
std::vector<int32_t> slots = channelPSet.getUntrackedParameter<std::vector<int32_t> >("slot");
std::vector<int32_t> fibers = channelPSet.getUntrackedParameter<std::vector<int32_t> >("fiber");
std::vector<int32_t> fiberChannels = channelPSet.getUntrackedParameter<std::vector<int32_t> >("fiber_channel");
for (unsigned int i = 0; i < crates.size(); ++i) {
_channels.push_back({crates[i], slots[i], fibers[i], fiberChannels[i]});
}
}
/* virtual */ void FCDTask::bookHistograms(DQMStore::IBooker& ib, edm::Run const& r, edm::EventSetup const& es) {
edm::ESHandle<HcalDbService> dbService = es.getHandle(hcalDbServiceToken_);
_emap = dbService->getHcalMapping();
_ehashmap.initialize(_emap, hcaldqm::electronicsmap::fD2EHashMap);
ib.cd();
//book histos per channel
std::string histoname;
std::vector<HcalGenericDetId> gids = _emap->allPrecisionId();
for (auto& it_gid : gids) {
if (it_gid.genericSubdet() != HcalGenericDetId::HcalGenZDC) {
continue;
}
HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(it_gid));
for (auto& it_channel : _channels) {
if ((eid.crateId() == it_channel.crate) && (eid.slot() == it_channel.slot) &&
(eid.fiberIndex() == it_channel.fiber) && (eid.fiberChanId() == it_channel.fiberChannel)) {
_fcd_eids.push_back(eid);
}
}
}
for (auto& it_eid : _fcd_eids) {
// EM Pos
histoname = std::to_string(it_eid.crateId()) + "-" + std::to_string(it_eid.slot()) + "-" +
std::to_string(it_eid.fiberIndex()) + "-" + std::to_string(it_eid.fiberChanId());
ib.setCurrentFolder("Hcal/FCDTask/ADC");
_cADC[it_eid] = ib.book1DD(histoname.c_str(), histoname.c_str(), 256, 0, 256);
_cADC[it_eid]->setAxisTitle("ADC", 1);
_cADC[it_eid]->setAxisTitle("N", 2);
ib.setCurrentFolder("Hcal/FCDTask/ADC_vs_TS"),
_cADC_vs_TS[it_eid] = ib.book2D(histoname.c_str(), histoname.c_str(), 10, 0, 10, 64, 0, 256);
_cADC_vs_TS[it_eid]->setAxisTitle("TS", 1);
_cADC_vs_TS[it_eid]->setAxisTitle("ADC", 2);
ib.setCurrentFolder("Hcal/FCDTask/TDCTime");
_cTDCTime[it_eid] = ib.book1DD(histoname.c_str(), histoname.c_str(), 500, 0., 250.);
_cTDCTime[it_eid]->setAxisTitle("TDC time [ns]", 1);
ib.setCurrentFolder("Hcal/FCDTask/TDC");
_cTDC[it_eid] = ib.book1DD(histoname.c_str(), histoname.c_str(), 64, -0.5, 63.5);
_cTDC[it_eid]->setAxisTitle("TDC", 1);
}
}
/* virtual */ void FCDTask::analyze(edm::Event const& e, edm::EventSetup const&) {
edm::Handle<QIE10DigiCollection> digis;
if (!e.getByToken(_tokQIE10, digis))
edm::LogError("Collection QIE10DigiCollection for ZDC isn't available" + _tagQIE10.label() + " " +
_tagQIE10.instance());
for (auto it = digis->begin(); it != digis->end(); it++) {
const QIE10DataFrame digi = static_cast<const QIE10DataFrame>(*it);
HcalGenericDetId const& gdid = digi.detid();
HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(gdid));
if (std::find(_fcd_eids.begin(), _fcd_eids.end(), eid) == _fcd_eids.end()) {
continue;
}
for (int i = 0; i < digi.samples(); i++) {
// iter over all samples
_cADC[eid]->Fill(digi[i].adc());
_cADC_vs_TS[eid]->Fill(i, digi[i].adc());
_cTDC[eid]->Fill(digi[i].le_tdc());
if (digi[i].le_tdc() <= 50.) {
double tdctime = 25. * i + 0.5 * digi[i].le_tdc();
_cTDCTime[eid]->Fill(tdctime);
}
}
}
}
DEFINE_FWK_MODULE(FCDTask);
|