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
|
#include "DataFormats/HcalDigi/interface/CastorDataFrame.h"
CastorDataFrame::CastorDataFrame() : id_(0), size_(0), hcalPresamples_(0) {}
CastorDataFrame::CastorDataFrame(const HcalCastorDetId& id) : id_(id), size_(0), hcalPresamples_(0) {
// TODO : test id for CASTOR
}
void CastorDataFrame::setSize(int size) {
if (size > MAXSAMPLES)
size_ = MAXSAMPLES;
else if (size <= 0)
size_ = 0;
else
size_ = size;
}
void CastorDataFrame::setPresamples(int ps) { hcalPresamples_ |= ps & 0xF; }
//void CastorDataFrame::setReadoutIds(const HcalElectronicsId& eid) {
// electronicsId_=eid;
//}
bool CastorDataFrame::validate(int firstSample, int nSamples) const {
int capid = -1;
bool ok = true;
for (int i = 0; ok && i < nSamples && i + firstSample < size_; i++) {
if (data_[i + firstSample].er() || !data_[i + firstSample].dv())
ok = false;
if (i == 0)
capid = data_[i + firstSample].capid();
if (capid != data_[i + firstSample].capid())
ok = false;
capid = (capid + 1) % 4;
}
return ok;
}
void CastorDataFrame::setZSInfo(bool unsuppressed, bool markAndPass, uint32_t crossingMask) {
hcalPresamples_ &= 0x7FC00F0F; // preserve actual presamples and fiber idle offset
if (markAndPass)
hcalPresamples_ |= 0x10;
if (unsuppressed)
hcalPresamples_ |= 0x20;
hcalPresamples_ |= (crossingMask & 0x3FF) << 12;
}
int CastorDataFrame::fiberIdleOffset() const {
int val = (hcalPresamples_ & 0xF00) >> 8;
return (val == 0) ? (-1000) : (((val & 0x8) == 0) ? (-(val & 0x7)) : (val & 0x7));
}
void CastorDataFrame::setFiberIdleOffset(int offset) {
hcalPresamples_ &= 0x7FFFF0FF;
if (offset >= 7)
hcalPresamples_ |= 0xF00;
else if (offset >= 0)
hcalPresamples_ |= (0x800) | (offset << 8);
else if (offset >= -7)
hcalPresamples_ |= ((-offset) << 8);
else
hcalPresamples_ |= 0x700;
}
std::ostream& operator<<(std::ostream& s, const CastorDataFrame& digi) {
s << digi.id() << " " << digi.size() << " samples " << digi.presamples() << " presamples " << std::endl;
if (digi.fiberIdleOffset() != 0) {
if (digi.fiberIdleOffset() == -1000)
s << " nofiberOffset";
else
s << " fiberOffset=" << digi.fiberIdleOffset();
}
for (int i = 0; i < digi.size(); i++)
s << " " << digi.sample(i) << std::endl;
return s;
}
|