|
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
|
#include "Geometry/MuonNumbering/interface/DD4hep_MuonNumbering.h"
#include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
using namespace cms;
using namespace edm;
const MuonBaseNumber MuonNumbering::geoHistoryToBaseNumber(const cms::ExpandedNodes& nodes) const {
MuonBaseNumber num;
int levelPart = get("level");
int superPart = get("super");
int basePart = get("base");
int startCopyNo = get("xml_starts_with_copyno");
// some consistency checks
if (basePart != 1) {
edm::LogError("Geometry") << "MuonNumbering finds unusual base constant: " << basePart;
}
if (superPart < 100) {
edm::LogError("Geometry") << "MuonNumbering finds unusual super part: " << superPart
<< " -- resetting super part to 100";
superPart = 100; // Reset to sensible value so calculations below don't go out of range
}
if (levelPart < 10 * superPart) {
edm::LogError("Geometry") << "MuonNumbering finds unusual level constant: " << levelPart;
}
if ((startCopyNo != 0) && (startCopyNo != 1)) {
edm::LogError("Geometry") << "MuonNumbering finds unusual start value for copy numbers: " << startCopyNo;
}
int ctr(0);
for (auto const& it : nodes.tags) {
int tag = it / levelPart;
if (tag > 0) {
int offset = nodes.offsets[ctr];
int copyno = nodes.copyNos[ctr] + offset % superPart;
int super = offset / superPart;
num.addBase(tag, super, copyno - startCopyNo);
}
++ctr;
}
return num;
}
const int MuonNumbering::get(const char* key) const {
int result(0);
auto const& it = values_.find(key);
if (it != end(values_))
result = it->second;
return result;
}
void MuonNumbering::put(std::string_view str, int num) { values_.emplace(str, num); }
#include "FWCore/Utilities/interface/typelookup.h"
TYPELOOKUP_DATA_REG(MuonNumbering);
|