|
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
|
#include "CondFormats/SiPixelObjects/interface/SiPixelCalibConfiguration.h"
#include "CalibFormats/SiPixelObjects/interface/PixelCalibConfiguration.h"
#include <cstdlib>
short SiPixelCalibConfiguration::vcalIndexForEvent(const uint32_t &eventnumber) const {
uint32_t relative_event = std::abs((int32_t)eventnumber - 1) % patternSize();
short relative_pattern = relative_event / getNTriggers();
return relative_pattern;
}
short SiPixelCalibConfiguration::vcalForEvent(const uint32_t &eventnumber) const {
short result = fVCalValues[vcalIndexForEvent(eventnumber)];
return result;
}
std::vector<short> SiPixelCalibConfiguration::columnPatternForEvent(const uint32_t &eventnumber) const {
std::vector<short> result;
uint32_t patternnumber = eventnumber / patternSize();
uint32_t colpatternnumber = patternnumber % nColumnPatterns();
uint32_t nminuscol = 0;
for (size_t icol = 0; icol < fColumnPattern.size(); icol++) {
if (fColumnPattern[icol] == -1)
nminuscol++;
else if (nminuscol > colpatternnumber)
break;
else if (nminuscol == colpatternnumber) {
short val = fColumnPattern[icol];
result.push_back(val);
}
}
return result;
}
std::vector<short> SiPixelCalibConfiguration::rowPatternForEvent(const uint32_t &eventnumber) const {
std::vector<short> result;
uint32_t patternnumber = eventnumber / patternSize();
uint32_t rowpatternnumber = patternnumber / nColumnPatterns();
uint32_t nminusrow = 0;
for (size_t irow = 0; irow < fRowPattern.size(); irow++) {
if (fRowPattern[irow] == -1)
nminusrow++;
else if (nminusrow > rowpatternnumber)
break;
else if (nminusrow == rowpatternnumber) {
short val = fRowPattern[irow];
result.push_back(val);
}
}
return result;
}
uint32_t SiPixelCalibConfiguration::nextPatternChangeForEvent(const uint32_t &eventnumber) const {
uint32_t relative_event = eventnumber / patternSize();
relative_event += 1;
return relative_event * patternSize();
}
uint32_t SiPixelCalibConfiguration::expectedTotalEvents() const { return patternSize() * nPatterns(); }
SiPixelCalibConfiguration::SiPixelCalibConfiguration(const pos::PixelCalibConfiguration &fancyConfig)
: fNTriggers(0),
fRowPattern(std::vector<short>(0)),
fColumnPattern(std::vector<short>(0)),
fVCalValues(std::vector<short>(0)),
fMode("unknown")
{ // copy constructor that uses the complex object
fNTriggers = fancyConfig.nTriggersPerPattern();
std::vector<int> vcalpoints(0);
std::cout << "scan name = " << fancyConfig.scanName(0) << std::endl;
std::vector<uint32_t> vcalpointsuint32 = fancyConfig.scanValues(fancyConfig.scanName(0));
for (size_t ical = 0; ical < vcalpointsuint32.size(); ++ical) {
short vcalinput = vcalpointsuint32[ical];
std::cout << "Vcal value " << ical << " = " << vcalinput << std::endl;
fVCalValues.push_back(vcalinput);
}
// copy row and column patterns
std::vector<std::vector<uint32_t> > cols = fancyConfig.columnList();
std::vector<std::vector<uint32_t> > rows = fancyConfig.rowList();
for (uint32_t i = 0; i < cols.size(); ++i) {
for (uint32_t j = 0; j < cols[i].size(); ++j) {
short colval = cols[i][j];
fColumnPattern.push_back(colval);
}
fColumnPattern.push_back(-1);
}
for (uint32_t i = 0; i < rows.size(); ++i) {
for (uint32_t j = 0; j < rows[i].size(); ++j) {
short rowval = rows[i][j];
fRowPattern.push_back(rowval);
}
fRowPattern.push_back(-1);
}
fMode = fancyConfig.mode();
}
uint32_t SiPixelCalibConfiguration::nRowPatterns() const {
uint32_t nrows = 0;
for (std::vector<short>::const_iterator i = fRowPattern.begin(); i != fRowPattern.end(); ++i) {
if (*i == -1)
nrows++;
}
return nrows;
}
uint32_t SiPixelCalibConfiguration::nColumnPatterns() const {
uint32_t ncols = 0;
for (std::vector<short>::const_iterator i = fColumnPattern.begin(); i != fColumnPattern.end(); ++i) {
if (*i == -1)
ncols++;
}
return ncols;
}
|