|
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
|
/** \class QieShape
basic linearization function for HCAL QIE
$Author: ratnikov
*/
#include <iostream>
#include "CalibFormats/CastorObjects/interface/QieShape.h"
namespace reco {
namespace castor {
QieShape::QieShape(const double fAdcShape[32], const double fAdcBin[32]) {
for (int i = 0; i < 32; i++) { // initial settings
mLinearization[i] = fAdcShape[i];
mBinSize[i] = fAdcBin[i];
// std::cout << "QieShape::QieShape-> #/adc/bin: " << i << '/' << fAdcShape [i] << '/' << fAdcBin [i] << std::endl;
}
double factor = 1;
for (int range = 1; range < 4; range++) {
factor = factor * 5;
int offset = 32 * range;
mLinearization[offset] = mLinearization[offset - 2]; // initial overlap
for (int bin = 1; bin < 32; bin++) {
mLinearization[offset + bin] =
mLinearization[offset + bin - 1] +
factor * (mLinearization[bin] - mLinearization[bin - 1]); // scale initial curve
mBinSize[offset + bin] = factor * mBinSize[bin];
}
}
}
} // namespace castor
} // namespace reco
|