|
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
|
#include "CondFormats/SiStripObjects/interface/SiStripBadStrip.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CondFormats/SiStripObjects/interface/SiStripDetSummary.h"
#include <algorithm>
bool SiStripBadStrip::put(const uint32_t& DetId, Range input) {
// put in SiStripBadStrip::v_badstrips of DetId
Registry::iterator p = std::lower_bound(indexes.begin(), indexes.end(), DetId, SiStripBadStrip::StrictWeakOrdering());
if (p != indexes.end() && p->detid == DetId) {
edm::LogError("SiStripBadStrip") << "[" << __PRETTY_FUNCTION__ << "] SiStripBadStrip for DetID " << DetId
<< " is already stored. Skipping this put" << std::endl;
return false;
}
size_t sd = input.second - input.first;
DetRegistry detregistry;
detregistry.detid = DetId;
detregistry.ibegin = v_badstrips.size();
detregistry.iend = v_badstrips.size() + sd;
indexes.insert(p, detregistry);
v_badstrips.insert(v_badstrips.end(), input.first, input.second);
return true;
}
const SiStripBadStrip::Range SiStripBadStrip::getRange(const uint32_t DetId) const {
// get SiStripBadStrip Range of DetId
RegistryIterator p = std::lower_bound(indexes.begin(), indexes.end(), DetId, SiStripBadStrip::StrictWeakOrdering());
if (p == indexes.end() || p->detid != DetId)
return SiStripBadStrip::Range(v_badstrips.end(), v_badstrips.end());
else {
__builtin_prefetch((&v_badstrips.front()) + p->ibegin);
__builtin_prefetch((&v_badstrips.front()) + p->ibegin + 24);
__builtin_prefetch((&v_badstrips.front()) + p->iend - 24);
return SiStripBadStrip::Range(v_badstrips.begin() + p->ibegin, v_badstrips.begin() + p->iend);
}
}
SiStripBadStrip::Range SiStripBadStrip::getRangeByPos(unsigned short pos) const {
if (pos > indexes.size())
return Range(v_badstrips.end(), v_badstrips.end());
auto p = indexes.begin() + pos;
__builtin_prefetch((&v_badstrips.front()) + p->ibegin);
__builtin_prefetch((&v_badstrips.front()) + p->ibegin + 24);
__builtin_prefetch((&v_badstrips.front()) + p->iend - 24);
return Range(v_badstrips.begin() + p->ibegin, v_badstrips.begin() + p->iend);
}
void SiStripBadStrip::getDetIds(std::vector<uint32_t>& DetIds_) const {
// returns vector of DetIds in map
SiStripBadStrip::RegistryIterator begin = indexes.begin();
SiStripBadStrip::RegistryIterator end = indexes.end();
for (SiStripBadStrip::RegistryIterator p = begin; p != end; ++p) {
DetIds_.push_back(p->detid);
}
}
void SiStripBadStrip::printSummary(std::stringstream& ss, const TrackerTopology* trackerTopo) const {
SiStripDetSummary summaryBadModules{trackerTopo};
SiStripDetSummary summaryBadStrips{trackerTopo};
// Loop on the vector<DetRegistry> and take the bad modules and bad strips
Registry::const_iterator it = indexes.begin();
for (; it != indexes.end(); ++it) {
summaryBadModules.add(it->detid);
summaryBadStrips.add(it->iend - it->ibegin);
}
ss << "Summary of bad modules in detector:" << std::endl;
summaryBadModules.print(ss, false);
ss << "Summary of bad strip in detectors:" << std::endl;
summaryBadStrips.print(ss, false);
}
void SiStripBadStrip::printDebug(std::stringstream& ss, const TrackerTopology* /*trackerTopo*/) const {
ss << "Printing all bad strips for all DetIds" << std::endl;
// Loop on the vector<DetRegistry> and take the bad modules and bad strips
Registry::const_iterator it = indexes.begin();
for (; it != indexes.end(); ++it) {
// ss << "For DetId = " << it->detid << std::endl;
SiStripBadStrip::Range range(getRange(it->detid));
for (std::vector<unsigned int>::const_iterator badStrip = range.first; badStrip != range.second; ++badStrip) {
ss << "DetId=" << it->detid << " Strip=" << decode(*badStrip).firstStrip << ":" << decode(*badStrip).range
<< " flag=" << decode(*badStrip).flag << std::endl;
}
}
}
|