|
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
|
/** \class CastorChannelCoder
Container for ADC<->fQ conversion constants for HCAL/Castor QIE
$Original author: ratnikov
*/
#include <iostream>
#include "CalibFormats/CastorObjects/interface/CastorChannelCoder.h"
#include "CalibFormats/CastorObjects/interface/QieShape.h"
CastorChannelCoder::CastorChannelCoder(const float fOffset[16], const float fSlope[16]) { // [CapId][Range]
for (int range = 0; range < 4; range++) {
for (int capId = 0; capId < 4; capId++) {
mOffset[capId][range] = fOffset[index(capId, range)];
mSlope[capId][range] = fSlope[index(capId, range)];
}
}
}
double CastorChannelCoder::charge(const reco::castor::QieShape& fShape, int fAdc, int fCapId) const {
int range = (fAdc >> 6) & 0x3;
double charge = fShape.linearization(fAdc) / mSlope[fCapId][range] + mOffset[fCapId][range];
// std::cout << "CastorChannelCoder::charge-> " << fAdc << '/' << fCapId
// << " result: " << charge << std::endl;
return charge;
}
int CastorChannelCoder::adc(const reco::castor::QieShape& fShape, double fCharge, int fCapId) const {
int adc = -1; //nothing found yet
// search for the range
for (int range = 0; range < 4; range++) {
double qieCharge = (fCharge - mOffset[fCapId][range]) * mSlope[fCapId][range];
double qieChargeMax = fShape.linearization(32 * range + 31) + 0.5 * fShape.binSize(32 * range + 31);
if (range == 3 && qieCharge > qieChargeMax)
adc = 127; // overflow
if (qieCharge > qieChargeMax)
continue; // next range
for (int bin = 32 * range; bin < 32 * (range + 1); bin++) {
if (qieCharge < fShape.linearization(bin) + 0.5 * fShape.binSize(bin)) {
adc = bin;
break;
}
}
if (adc >= 0)
break; // found
}
if (adc < 0)
adc = 0; // underflow
// std::cout << "CastorChannelCoder::adc-> " << fCharge << '/' << fCapId
// << " result: " << adc << std::endl;
return adc;
}
|