|
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
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
|
#include "CondFormats/SiPixelObjects/interface/SiPixelGainCalibration.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <algorithm>
#include <cstring>
//
// Constructors
//
SiPixelGainCalibration::SiPixelGainCalibration()
: minPed_(0.),
maxPed_(255.),
minGain_(0.),
maxGain_(255.),
numberOfRowsToAverageOver_(1),
nBinsToUseForEncoding_(253),
deadFlag_(255),
noisyFlag_(254) {
if (deadFlag_ > 0xFF)
throw cms::Exception("GainCalibration Payload configuration error")
<< "[SiPixelGainCalibration::SiPixelGainCalibration] Dead flag was set to " << deadFlag_
<< ", and it must be set less than or equal to 255";
}
//
SiPixelGainCalibration::SiPixelGainCalibration(float minPed, float maxPed, float minGain, float maxGain)
: minPed_(minPed),
maxPed_(maxPed),
minGain_(minGain),
maxGain_(maxGain),
numberOfRowsToAverageOver_(1),
nBinsToUseForEncoding_(253),
deadFlag_(255),
noisyFlag_(254) {
if (deadFlag_ > 0xFF)
throw cms::Exception("GainCalibration Payload configuration error")
<< "[SiPixelGainCalibration::SiPixelGainCalibration] Dead flag was set to " << deadFlag_
<< ", and it must be set less than or equal to 255";
}
bool SiPixelGainCalibration::put(const uint32_t& DetId, Range input, const int& nCols) {
// put in SiPixelGainCalibration of DetId
Registry::iterator p =
std::lower_bound(indexes.begin(), indexes.end(), DetId, SiPixelGainCalibration::StrictWeakOrdering());
if (p != indexes.end() && p->detid == DetId)
return false;
size_t sd = input.second - input.first;
DetRegistry detregistry;
detregistry.detid = DetId;
detregistry.ncols = nCols;
detregistry.ibegin = v_pedestals.size();
detregistry.iend = v_pedestals.size() + sd;
indexes.insert(p, detregistry);
v_pedestals.insert(v_pedestals.end(), input.first, input.second);
return true;
}
const int SiPixelGainCalibration::getNCols(const uint32_t& DetId) const {
// get number of columns of DetId
RegistryIterator p =
std::lower_bound(indexes.begin(), indexes.end(), DetId, SiPixelGainCalibration::StrictWeakOrdering());
if (p == indexes.end() || p->detid != DetId)
return 0;
else
return p->ncols;
}
const SiPixelGainCalibration::Range SiPixelGainCalibration::getRange(const uint32_t& DetId) const {
// get SiPixelGainCalibration Range of DetId
RegistryIterator p =
std::lower_bound(indexes.begin(), indexes.end(), DetId, SiPixelGainCalibration::StrictWeakOrdering());
if (p == indexes.end() || p->detid != DetId)
return SiPixelGainCalibration::Range(v_pedestals.end(), v_pedestals.end());
else
return SiPixelGainCalibration::Range(v_pedestals.begin() + p->ibegin, v_pedestals.begin() + p->iend);
}
const std::pair<const SiPixelGainCalibration::Range, const int> SiPixelGainCalibration::getRangeAndNCols(
const uint32_t& DetId) const {
RegistryIterator p =
std::lower_bound(indexes.begin(), indexes.end(), DetId, SiPixelGainCalibration::StrictWeakOrdering());
if (p == indexes.end() || p->detid != DetId)
return std::make_pair(SiPixelGainCalibration::Range(v_pedestals.end(), v_pedestals.end()), 0);
else
return std::make_pair(SiPixelGainCalibration::Range(v_pedestals.begin() + p->ibegin, v_pedestals.begin() + p->iend),
p->ncols);
}
void SiPixelGainCalibration::getDetIds(std::vector<uint32_t>& DetIds_) const {
// returns vector of DetIds in map
SiPixelGainCalibration::RegistryIterator begin = indexes.begin();
SiPixelGainCalibration::RegistryIterator end = indexes.end();
for (SiPixelGainCalibration::RegistryIterator p = begin; p != end; ++p) {
DetIds_.push_back(p->detid);
}
}
void SiPixelGainCalibration::setData(
float ped, float gain, std::vector<char>& vped, bool isDeadPixel, bool isNoisyPixel) {
float theEncodedGain = 0;
float theEncodedPed = 0;
if (!isDeadPixel && !isNoisyPixel) {
theEncodedGain = encodeGain(gain);
theEncodedPed = encodePed(ped);
}
unsigned int ped_ = (static_cast<unsigned int>(theEncodedPed)) & 0xFF;
unsigned int gain_ = (static_cast<unsigned int>(theEncodedGain)) & 0xFF;
if (isDeadPixel) {
ped_ = deadFlag_ & 0xFF;
gain_ = deadFlag_ & 0xFF;
}
if (isNoisyPixel) {
ped_ = noisyFlag_ & 0xFF;
gain_ = noisyFlag_ & 0xFF;
}
unsigned int data = (ped_ << 8) | gain_;
vped.resize(vped.size() + 2);
// insert in vector of char
::memcpy((void*)(&vped[vped.size() - 2]), (void*)(&data), 2);
}
float SiPixelGainCalibration::getPed(
const int& col, const int& row, const Range& range, const int& nCols, bool& isDead, bool& isNoisy) const {
int nRows = (range.second - range.first) / 2 / nCols;
const unsigned int ped = *(range.first + 1 + (col * nRows + row) * 2) & 0xFF;
if (col >= nCols || row >= nRows) {
throw cms::Exception("CorruptedData")
<< "[SiPixelGainCalibration::getPed] Pixel out of range: col " << col << " row " << row;
}
if (ped == deadFlag_)
isDead = true;
if (ped == noisyFlag_)
isNoisy = true;
return decodePed(ped);
}
float SiPixelGainCalibration::getGain(
const int& col, const int& row, const Range& range, const int& nCols, bool& isDead, bool& isNoisy) const {
int nRows = (range.second - range.first) / 2 / nCols;
const unsigned int gain = *(range.first + (col * nRows + row) * 2) & 0xFF;
if (col >= nCols || row >= nRows) {
throw cms::Exception("CorruptedData")
<< "[SiPixelGainCalibration::getPed] Pixel out of range: col " << col << " row " << row;
}
if (gain == deadFlag_)
isDead = true;
if (gain == noisyFlag_)
isNoisy = true;
return decodeGain(gain);
}
float SiPixelGainCalibration::encodeGain(const float& gain) {
if (gain < minGain_ || gain > maxGain_) {
throw cms::Exception("InsertFailure") << "[SiPixelGainCalibration::encodeGain] Trying to encode gain (" << gain
<< ") out of range [" << minGain_ << "," << maxGain_ << "]\n";
} else {
double precision = (maxGain_ - minGain_) / static_cast<float>(nBinsToUseForEncoding_);
float encodedGain = (float)((gain - minGain_) / precision);
return encodedGain;
}
}
float SiPixelGainCalibration::encodePed(const float& ped) {
if (ped < minPed_ || ped > maxPed_) {
throw cms::Exception("InsertFailure") << "[SiPixelGainCalibration::encodePed] Trying to encode pedestal (" << ped
<< ") out of range [" << minPed_ << "," << maxPed_ << "]\n";
} else {
double precision = (maxPed_ - minPed_) / static_cast<float>(nBinsToUseForEncoding_);
float encodedPed = (float)((ped - minPed_) / precision);
return encodedPed;
}
}
float SiPixelGainCalibration::decodePed(unsigned int ped) const {
double precision = (maxPed_ - minPed_) / static_cast<float>(nBinsToUseForEncoding_);
float decodedPed = (float)(ped * precision + minPed_);
return decodedPed;
}
float SiPixelGainCalibration::decodeGain(unsigned int gain) const {
double precision = (maxGain_ - minGain_) / static_cast<float>(nBinsToUseForEncoding_);
float decodedGain = (float)(gain * precision + minGain_);
return decodedGain;
}
|