|
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
|
#include "CondFormats/SiPhase2TrackerObjects/interface/TrackerDetToDTCELinkCablingMap.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <utility>
#include <algorithm>
#include <iostream>
TrackerDetToDTCELinkCablingMap::TrackerDetToDTCELinkCablingMap() {}
TrackerDetToDTCELinkCablingMap::~TrackerDetToDTCELinkCablingMap() {}
std::unordered_map<DTCELinkId, uint32_t>::const_iterator TrackerDetToDTCELinkCablingMap::dtcELinkIdToDetId(
DTCELinkId const& key) const {
if (cablingMapDTCELinkIdToDetId_.find(key) == cablingMapDTCELinkIdToDetId_.end()) {
throw cms::Exception(
"TrackerDetToDTCELinkCablingMap has been asked to return a DetId associated to a DTCELinkId, but the latter is "
"unknown to the map. ")
<< " (DTC, GBT, Elink) numbers = (" << key.dtc_id() << "," << key.gbtlink_id() << "," << key.elink_id() << ")"
<< std::endl;
}
return cablingMapDTCELinkIdToDetId_.find(key);
}
std::pair<std::unordered_multimap<uint32_t, DTCELinkId>::const_iterator,
std::unordered_multimap<uint32_t, DTCELinkId>::const_iterator>
TrackerDetToDTCELinkCablingMap::detIdToDTCELinkId(uint32_t const key) const {
auto const DTCELinkId_itpair = cablingMapDetIdToDTCELinkId_.equal_range(key);
if (DTCELinkId_itpair.first == cablingMapDetIdToDTCELinkId_.end()) {
throw cms::Exception(
"TrackerDetToDTCELinkCablingMap has been asked to return a DTCELinkId associated to a DetId, but the latter is "
"unknown to the map. ")
<< " DetId = " << key << std::endl;
}
return DTCELinkId_itpair;
}
bool TrackerDetToDTCELinkCablingMap::knowsDTCELinkId(DTCELinkId const& key) const {
return cablingMapDTCELinkIdToDetId_.find(key) != cablingMapDTCELinkIdToDetId_.end();
}
bool TrackerDetToDTCELinkCablingMap::knowsDetId(uint32_t key) const {
return cablingMapDetIdToDTCELinkId_.find(key) != cablingMapDetIdToDTCELinkId_.end();
}
std::vector<DTCELinkId> TrackerDetToDTCELinkCablingMap::getKnownDTCELinkIds() const {
std::vector<DTCELinkId> knownDTCELinkIds(cablingMapDTCELinkIdToDetId_.size());
// Unzip the map into a vector of DTCELinkId, discarding the DetIds
std::transform(cablingMapDTCELinkIdToDetId_.begin(),
cablingMapDTCELinkIdToDetId_.end(),
knownDTCELinkIds.begin(),
[=](auto pair) { return pair.first; });
return knownDTCELinkIds;
}
std::vector<uint32_t> TrackerDetToDTCELinkCablingMap::getKnownDetIds() const {
std::vector<uint32_t> knownDetId;
// To get the list of unique DetIds we need to iterate over the various equal_ranges
// in the map associated to each unique key, and count them only once.
for (auto allpairs_it = cablingMapDetIdToDTCELinkId_.begin(), allpairs_end = cablingMapDetIdToDTCELinkId_.end();
allpairs_it != allpairs_end;) {
// ***Store the first instance of the key***
knownDetId.push_back(uint32_t(allpairs_it->first));
// *** Skip to the end of the equal range ***
// The following is just explicative, the bottom expression is equivalent
//auto const current_key = allpairs_it->first;
//auto const current_key_equal_range = cablingMapDetIdToDTCELinkId_.equal_range(current_key);
//auto const current_key_range_end = current_key_equal_range.second;
auto const current_key_range_end = cablingMapDetIdToDTCELinkId_.equal_range(allpairs_it->first).second;
while (allpairs_it != current_key_range_end)
++allpairs_it;
}
return knownDetId;
}
void TrackerDetToDTCELinkCablingMap::insert(DTCELinkId const& dtcELinkId, uint32_t const detId) {
cablingMapDTCELinkIdToDetId_.insert(std::make_pair(DTCELinkId(dtcELinkId), uint32_t(detId)));
cablingMapDetIdToDTCELinkId_.insert(std::make_pair(uint32_t(detId), DTCELinkId(dtcELinkId)));
}
void TrackerDetToDTCELinkCablingMap::clear() {
cablingMapDTCELinkIdToDetId_.clear();
cablingMapDetIdToDTCELinkId_.clear();
}
|