|
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
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
|
#include "CondFormats/PPSObjects/interface/CTPPSPixelGainCalibration.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <algorithm>
#include <cstring>
//
// Constructors
//
CTPPSPixelGainCalibration::CTPPSPixelGainCalibration() : minPed_(0.), maxPed_(255.), minGain_(0.), maxGain_(255.) {
DetRegistry myreg;
myreg.detid = 0;
myreg.ibegin = 0;
myreg.iend = 0;
myreg.ncols = 0;
myreg.nrows = 0;
indexes = myreg;
edm::LogInfo("CTPPSPixelGainCalibration") << "Default instance of CTPPSPixelGainCalibration";
}
CTPPSPixelGainCalibration::CTPPSPixelGainCalibration(const uint32_t& detId,
const uint32_t& sensorSize = 24960,
const uint32_t& nCols = 156)
: minPed_(0.), maxPed_(255.), minGain_(0.), maxGain_(255.) {
DetRegistry myreg;
myreg.detid = detId;
myreg.ibegin = 0;
myreg.iend = sensorSize;
myreg.ncols = nCols;
myreg.nrows = sensorSize / nCols;
indexes = myreg;
edm::LogInfo("CTPPSPixelGainCalibration")
<< "Instance of CTPPSPixelGainCalibration setting detId = " << detId << " and npix = " << sensorSize
<< " and ncols = " << nCols << " and nrows= " << myreg.nrows;
}
//
CTPPSPixelGainCalibration::CTPPSPixelGainCalibration(float minPed, float maxPed, float minGain, float maxGain)
: minPed_(minPed), maxPed_(maxPed), minGain_(minGain), maxGain_(maxGain) {
DetRegistry myreg;
myreg.detid = 0;
myreg.ibegin = 0;
myreg.iend = 0;
myreg.ncols = 0;
myreg.nrows = 0;
indexes = myreg;
edm::LogInfo("CTPPSPixelGainCalibration") << "Instance of CTPPSPixelGainCalibration setting mininums and maximums";
}
CTPPSPixelGainCalibration::CTPPSPixelGainCalibration(const uint32_t& detid,
const std::vector<float>& peds,
const std::vector<float>& gains,
float minPed,
float maxPed,
float minGain,
float maxGain)
: minPed_(minPed), maxPed_(maxPed), minGain_(minGain), maxGain_(maxGain) {
setGainsPeds(detid, peds, gains);
edm::LogInfo("CTPPSPixelGainCalibration")
<< "Instance of CTPPSPixelGainCalibration setting detId = " << detid << "and vectors of peds and gains";
}
void CTPPSPixelGainCalibration::setGainsPeds(const uint32_t& detid,
const std::vector<float>& peds,
const std::vector<float>& gains) {
int sensorSize = peds.size();
int gainsSize = gains.size();
if (gainsSize != sensorSize)
throw cms::Exception("CTPPSPixelGainCalibration payload configuration error")
<< "[CTPPSPixelGainCalibration::CTPPSPixelGainCalibration] peds and gains vector sizes don't match for detid "
<< detid << " size peds = " << sensorSize << " size gains = " << gainsSize;
DetRegistry myreg;
myreg.detid = detid;
myreg.ibegin = 0;
myreg.iend = sensorSize;
myreg.ncols =
sensorSize /
160; // each ROC is made of 80 rows and 52 columns, each sensor is made of 160 rows and either 104 or 156 columns (2x2 or 2x3 ROCs)
myreg.nrows = 160; // for now this is hardwired.
indexes = myreg;
for (int i = 0; i < sensorSize; ++i)
putData(i, peds[i], gains[i]);
}
void CTPPSPixelGainCalibration::putData(uint32_t ipix, float ped, float gain) {
if (v_pedestals.size() > ipix) //if the vector is too big, this pixel has already been set
{
if (ped >= 0 && gain >= 0)
throw cms::Exception("CTPPSPixelGainCalibration fill error")
<< "[CTPPSPixelGainCalibration::putData] attempting to fill the vectors that are already filled detid = "
<< indexes.detid << " ipix " << ipix;
else // in case it is for setting the noisy or dead flag of already filled pixel you are allowed to reset it
{
edm::LogInfo("CTPPSPixelGainCalibration") << "resetting pixel calibration for noisy or dead flag";
resetPixData(ipix, ped, gain);
}
} else if (v_pedestals.size() < ipix)
throw cms::Exception("CTPPSPixelGainCalibration fill error")
<< "[CTPPSPixelGainCalibration::putData] the order got scrambled detid = " << indexes.detid << " ipix " << ipix;
else { //the size has to match exactly the index, the index - 0 pixel starts the vector, the one with index 1 should be pushed back in a vector of size== 1 (to become size==2) so on and o forth
v_pedestals.push_back(ped);
v_gains.push_back(gain);
}
}
void CTPPSPixelGainCalibration::setIndexes(const uint32_t& detid) {
indexes.detid = detid;
indexes.ibegin = 0;
indexes.iend = v_pedestals.size();
indexes.ncols = indexes.iend / 160; // nrows is 160 in CTPPS
indexes.nrows = 160;
}
void CTPPSPixelGainCalibration::resetPixData(uint32_t ipix, float ped, float gain) {
if (v_pedestals.size() <= ipix) {
edm::LogError("CTPPSPixelGainCalibration") << "this element ipix = " << ipix << " has not been set yet";
return;
}
v_pedestals[ipix] = ped; //modify value to already exisiting element
v_gains[ipix] = gain;
}
float CTPPSPixelGainCalibration::getPed(const int& col, const int& row) const {
int nCols = indexes.ncols;
int nRows = indexes.nrows;
if (col >= nCols || row >= nRows) {
throw cms::Exception("CorruptedData")
<< "[CTPPSPixelGainCalibration::getPed] Pixel out of range: col " << col << " row " << row;
}
int ipix = col + row * nCols;
float tmpped = getPed(ipix /*,isDead,isNoisy*/);
return tmpped;
}
float CTPPSPixelGainCalibration::getGain(const int& col, const int& row) const {
int nCols = indexes.ncols;
int nRows = indexes.nrows;
if (col >= nCols || row >= nRows) {
throw cms::Exception("CorruptedData")
<< "[CTPPSPixelGainCalibration::getGain] Pixel out of range: col " << col << " row " << row;
}
int ipix = col + row * nCols;
float tmpgain = getGain(ipix);
return tmpgain;
}
|