|
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
|
/**
\class CastorQIEData
\author Panos Katsas (UoA)
POOL object to store QIE coder parameters for one channel
*/
#include <iostream>
#include "CondFormats/CastorObjects/interface/CastorQIEShape.h"
#include "CondFormats/CastorObjects/interface/CastorQIECoder.h"
namespace {
// pack range/capId in the plain index
unsigned index(unsigned fRange, unsigned fCapId) { return fCapId * 4 + fRange; }
} // namespace
float CastorQIECoder::charge(const CastorQIEShape& fShape, unsigned fAdc, unsigned fCapId) const {
unsigned range = fShape.range(fAdc);
return (fShape.center(fAdc) - offset(fCapId, range)) / slope(fCapId, range);
}
unsigned CastorQIECoder::adc(const CastorQIEShape& fShape, float fCharge, unsigned fCapId) const {
// search for the range
for (unsigned range = 0; range < 4; range++) {
float qieCharge = fCharge * slope(fCapId, range) + offset(fCapId, range);
unsigned minBin = 32 * range;
float qieChargeMax = fShape.highEdge(minBin + 31);
if (qieCharge <= qieChargeMax) {
for (unsigned bin = minBin; bin <= minBin + 31; bin++) {
if (qieCharge < fShape.highEdge(bin)) {
return bin;
}
}
return minBin; // underflow
} else if (range == 3) {
return 127; // overflow
}
}
return 0; //should never get here
}
float CastorQIECoder::offset(unsigned fCapId, unsigned fRange) const { return *((&mOffset00) + index(fRange, fCapId)); }
float CastorQIECoder::slope(unsigned fCapId, unsigned fRange) const { return *((&mSlope00) + index(fRange, fCapId)); }
void CastorQIECoder::setOffset(unsigned fCapId, unsigned fRange, float fValue) {
if (fCapId < 4U && fRange < 4U) { // fCapId >= 0 and fRange >= 0, since fCapId and fRange are unsigned
*((&mOffset00) + index(fRange, fCapId)) = fValue;
} else {
std::cerr << "CastorQIECoder::setOffset-> Wrong parameters capid/range: " << fCapId << '/' << fRange << std::endl;
}
}
void CastorQIECoder::setSlope(unsigned fCapId, unsigned fRange, float fValue) {
if (fCapId < 4U && fRange < 4U) { // fCapId >= 0 and fRange >= 0, since fCapId and fRange are unsigned
*((&mSlope00) + index(fRange, fCapId)) = fValue;
} else {
std::cerr << "CastorQIECoder::setSlope-> Wrong parameters capid/range: " << fCapId << '/' << fRange << std::endl;
}
}
|