|
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
|
#include "DataFormats/HcalRecHit/interface/HBHERecHit.h"
#include "DataFormats/HcalRecHit/interface/HBHERecHitAuxSetter.h"
#include "DataFormats/HcalRecHit/interface/CaloRecHitAuxSetter.h"
HcalDetId HBHERecHit::idFront() const {
if (auxPhase1_ & (1U << HBHERecHitAuxSetter::OFF_COMBINED)) {
const HcalDetId myId(id());
return HcalDetId(myId.subdet(), myId.ieta(), myId.iphi(), auxHBHE_ & 0xf);
} else {
return id();
}
}
bool HBHERecHit::isMerged() const { return auxPhase1_ & (1U << HBHERecHitAuxSetter::OFF_COMBINED); }
std::ostream& operator<<(std::ostream& s, const HBHERecHit& hit) {
s << hit.id() << ": " << hit.energy() << " GeV";
if (hit.eraw() > -0.9e21) {
s << ", eraw=" << hit.eraw() << " GeV";
}
if (hit.eaux() > -0.9e21) {
s << ", eaux=" << hit.eaux() << " GeV";
}
if (hit.time() > -998) {
s << ", t= " << hit.time() << " to " << hit.timeFalling() << " ns";
}
return s;
}
void HBHERecHit::getMergedIds(std::vector<HcalDetId>* ids) const {
if (ids) {
ids->clear();
if (auxPhase1_ & (1U << HBHERecHitAuxSetter::OFF_COMBINED)) {
const unsigned nMerged = CaloRecHitAuxSetter::getField(
auxPhase1_, HBHERecHitAuxSetter::MASK_NSAMPLES, HBHERecHitAuxSetter::OFF_NSAMPLES);
ids->reserve(nMerged);
const HcalDetId myId(id());
for (unsigned i = 0; i < nMerged; ++i) {
const unsigned depth = CaloRecHitAuxSetter::getField(auxHBHE_, 0xf, i * 4);
ids->emplace_back(myId.subdet(), myId.ieta(), myId.iphi(), depth);
}
}
}
}
|