|
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
|
/**
\class CastorQIEData
\author Panos Katsas (UoA)
POOL object to store pedestal values 4xCapId
*/
#include "CondFormats/CastorObjects/interface/CastorQIEShape.h"
namespace {
const float binMin[32] = {-1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14,
16, 18, 20, 22, 24, 26, 28, 31, 34, 37, 40, 44, 48, 52, 57, 62};
}
CastorQIEShape::CastorQIEShape() {
for (int i = 0; i < 32; i++)
mValues[i] = binMin[i];
expand();
}
CastorQIEShape::~CastorQIEShape() {}
void CastorQIEShape::expand() {
int scale = 1;
for (unsigned range = 1; range < 4; range++) {
scale = scale * 5;
unsigned index = range * 32;
mValues[index] = mValues[index - 2]; // link to previous range
for (unsigned i = 1; i < 32; i++) {
mValues[index + i] = mValues[index + i - 1] + scale * (mValues[i] - mValues[i - 1]);
}
}
mValues[128] = 2 * mValues[127] - mValues[126]; // extrapolate
}
float CastorQIEShape::lowEdge(unsigned fAdc) const {
if (fAdc < 128)
return mValues[fAdc];
return 0.;
}
float CastorQIEShape::center(unsigned fAdc) const {
if (fAdc < 128) {
if (fAdc % 32 == 31)
return 0.5 * (3 * mValues[fAdc] - mValues[fAdc - 1]); // extrapolate
else
return 0.5 * (mValues[fAdc] + mValues[fAdc + 1]); // interpolate
}
return 0.;
}
float CastorQIEShape::highEdge(unsigned fAdc) const {
if (fAdc >= 128)
return 0;
if (fAdc == 127)
return mValues[fAdc + 1];
if (fAdc % 32 == 31)
return mValues[fAdc + 3];
return mValues[fAdc + 1];
}
bool CastorQIEShape::setLowEdge(float fValue, unsigned fAdc) {
if (fAdc >= 32)
return false;
mValues[fAdc] = fValue;
return true;
}
bool CastorQIEShape::setLowEdges(const float fValue[32]) {
bool result = true;
for (int adc = 0; adc < 32; adc++)
result = result && setLowEdge(fValue[adc], adc);
expand();
return result;
}
|