|
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
|
#include "Fireworks/Core/interface/FWSimpleProxyBuilderTemplate.h"
// #include "Fireworks/Calo/interface/FWHeatmapProxyBuilderTemplate.h"
#include "Fireworks/Core/interface/FWEventItem.h"
#include "Fireworks/Core/interface/FWProxyBuilderConfiguration.h"
#include "DataFormats/ForwardDetId/interface/HGCScintillatorDetId.h"
#include "DataFormats/ForwardDetId/interface/HGCSiliconDetId.h"
#include "DataFormats/ForwardDetId/interface/HGCalTriggerDetId.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"
template <typename T>
class FWL1THGCalProxyTemplate : public FWSimpleProxyBuilderTemplate<T> {
protected:
void setItem(const FWEventItem *iItem) override {
FWProxyBuilderBase::setItem(iItem);
if (iItem) {
iItem->getConfig()->keepEntries(true);
iItem->getConfig()->assertParam("Layer", 0L, 0L, 52L);
iItem->getConfig()->assertParam("EnergyCutOff", 0.5, 0.2, 5.0);
iItem->getConfig()->assertParam("Heatmap", true);
iItem->getConfig()->assertParam("Z+", true);
iItem->getConfig()->assertParam("Z-", true);
}
}
std::unordered_set<unsigned> getCellsFromTriggerCell(const unsigned trigger_cell_id) const {
constexpr unsigned hSc_triggercell_size_ = 2;
// constexpr unsigned hSc_module_size_ = 12; // in TC units (144 TC / panel = 36 e-links)
DetId trigger_cell_det_id(trigger_cell_id);
unsigned det = trigger_cell_det_id.det();
std::unordered_set<unsigned> cell_det_ids;
// Scintillator
if (det == DetId::HGCalHSc) {
HGCScintillatorDetId trigger_cell_sc_id(trigger_cell_id);
int ieta0 = (trigger_cell_sc_id.ietaAbs() - 1) * hSc_triggercell_size_ + 1;
int iphi0 = (trigger_cell_sc_id.iphi() - 1) * hSc_triggercell_size_ + 1;
for (int ietaAbs = ieta0; ietaAbs < ieta0 + (int)hSc_triggercell_size_; ietaAbs++) {
int ieta = ietaAbs * trigger_cell_sc_id.zside();
for (int iphi = iphi0; iphi < iphi0 + (int)hSc_triggercell_size_; iphi++) {
unsigned cell_id = HGCScintillatorDetId(trigger_cell_sc_id.type(), trigger_cell_sc_id.layer(), ieta, iphi);
#if 0
if (validCellId(DetId::HGCalHSc, cell_id))
#endif
cell_det_ids.emplace(cell_id);
}
}
}
// Silicon
else {
HGCalTriggerDetId trigger_cell_trig_id(trigger_cell_id);
unsigned subdet = trigger_cell_trig_id.subdet();
if (subdet == HGCalTriggerSubdetector::HGCalEETrigger || subdet == HGCalTriggerSubdetector::HGCalHSiTrigger) {
DetId::Detector cell_det =
(subdet == HGCalTriggerSubdetector::HGCalEETrigger ? DetId::HGCalEE : DetId::HGCalHSi);
int layer = trigger_cell_trig_id.layer();
int zside = trigger_cell_trig_id.zside();
int type = trigger_cell_trig_id.type();
int waferu = trigger_cell_trig_id.waferU();
int waferv = trigger_cell_trig_id.waferV();
std::vector<int> cellus = trigger_cell_trig_id.cellU();
std::vector<int> cellvs = trigger_cell_trig_id.cellV();
for (unsigned ic = 0; ic < cellus.size(); ic++) {
HGCSiliconDetId cell_det_id(cell_det, zside, type, layer, waferu, waferv, cellus[ic], cellvs[ic]);
cell_det_ids.emplace(cell_det_id.rawId());
}
}
}
return cell_det_ids;
}
};
|