89 int section, layer, sector, plane, strip, stripGlobal;
90 std::map<int, StripData>::iterator it;
91 std::vector<StripData> stripDataVector;
92 std::vector<StripData>::iterator it2, it3, it4;
94 stripDataVector.push_back(it->second);
95 sort(stripDataVector.begin(), stripDataVector.end(), compareBadDigitRate);
96 printf(
"EKLM data checker report.\n"
97 "Strips with readout errors sorted by error rate:\n");
98 it2 = stripDataVector.begin();
99 while (it2 != stripDataVector.end()) {
100 if (it2->nBadDigits == 0)
103 while (it3 != stripDataVector.end()) {
104 if (it3->nBadDigits != it2->nBadDigits || it3->nDigits != it2->nDigits)
108 sort(it2, it3, compareStripNumber);
109 for (it4 = it2; it4 != it3; ++it4) {
111 it4->strip, §ion, &layer, §or, &plane, &strip);
112 printf(
"Section %d, layer %d, sector %d, plane %d, strip %d: %.1f%% "
114 section, layer, sector, plane, strip,
115 float(it4->nBadDigits) / it4->nDigits * 100,
116 it4->nBadDigits, it4->nDigits);
120 printf(
"Strips with no data collected:\n");
121 for (section = 1; section <=
m_GeoDat->getNSections(); section++) {
122 for (layer = 1; layer <=
m_GeoDat->getNDetectorLayers(section); layer++) {
123 for (sector = 1; sector <=
m_GeoDat->getNSectors(); sector++) {
124 for (plane = 1; plane <=
m_GeoDat->getNPlanes(); plane++) {
125 for (strip = 1; strip <=
m_GeoDat->getNStrips(); strip++) {
127 section, layer, sector, plane, strip);
130 printf(
"Section %d, layer %d, sector %d, plane %d, strip %d.\n",
131 section, layer, sector, plane, strip);