|
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
|
#include "Geometry/MuonNumbering/interface/MuonGeometryNumbering.h"
#include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
#include "Geometry/MuonNumbering/interface/MuonGeometryConstants.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
//#define EDM_ML_DEBUG
MuonGeometryNumbering::MuonGeometryNumbering(const MuonGeometryConstants &muonConstants) {
// Get constant values from muonConstants
theLevelPart = muonConstants.getValue("level");
theSuperPart = muonConstants.getValue("super");
theBasePart = muonConstants.getValue("base");
theStartCopyNo = muonConstants.getValue("xml_starts_with_copyno");
// some consistency checks
if (theBasePart != 1) {
edm::LogWarning("MuonGeom") << "MuonGeometryNumbering finds unusual base constant:" << theBasePart;
}
if (theSuperPart < 100) {
edm::LogWarning("MuonGeom") << "MuonGeometryNumbering finds unusual super constant:" << theSuperPart;
}
if (theLevelPart < 10 * theSuperPart) {
edm::LogWarning("MuonGeom") << "MuonGeometryNumbering finds unusual level constant:" << theLevelPart;
}
if ((theStartCopyNo != 0) && (theStartCopyNo != 1)) {
edm::LogWarning("MuonGeom") << "MuonGeometryNumbering finds unusual start value for copy numbers:"
<< theStartCopyNo;
}
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "MuonGeometryNumbering configured with"
<< " Level = " << theLevelPart << " Super = " << theSuperPart
<< " Base = " << theBasePart << " StartCopyNo = " << theStartCopyNo;
#endif
}
MuonBaseNumber MuonGeometryNumbering::geoHistoryToBaseNumber(const DDGeoHistory &history) const {
MuonBaseNumber num;
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "MuonGeometryNumbering create MuonBaseNumber for " << history;
#endif
//loop over all parents and check
DDGeoHistory::const_iterator cur = history.begin();
DDGeoHistory::const_iterator end = history.end();
while (cur != end) {
const DDLogicalPart &ddlp = cur->logicalPart();
const int tag = getInt("CopyNoTag", ddlp) / theLevelPart;
if (tag > 0) {
const int offset = getInt("CopyNoOffset", ddlp);
const int copyno = (cur->copyno()) + offset % theSuperPart;
const int super = offset / theSuperPart;
num.addBase(tag, super, copyno - theStartCopyNo);
}
cur++;
}
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "MuonGeometryNumbering::" << num.getLevels();
for (int i = 1; i <= num.getLevels(); i++) {
edm::LogVerbatim("MuonGeom") << "[" << i << "] " << num.getSuperNo(i) << " " << num.getBaseNo(i);
}
#endif
return num;
}
MuonBaseNumber MuonGeometryNumbering::geoHistoryToBaseNumber(const cms::ExpandedNodes &nodes) const {
MuonBaseNumber num;
int ctr(0);
for (auto const &it : nodes.tags) {
int tag = it / theLevelPart;
if (tag > 0) {
int offset = nodes.offsets[ctr];
int copyno = nodes.copyNos[ctr] + offset % theSuperPart;
int super = offset / theSuperPart;
num.addBase(tag, super, copyno - theStartCopyNo);
}
++ctr;
}
return num;
}
int MuonGeometryNumbering::getInt(const std::string &s, const DDLogicalPart &part) const {
DDValue val(s);
std::vector<const DDsvalues_type *> result = part.specifics();
std::vector<const DDsvalues_type *>::iterator it = result.begin();
bool foundIt = false;
for (; it != result.end(); ++it) {
foundIt = DDfetch(*it, val);
if (foundIt)
break;
}
if (foundIt) {
std::vector<double> temp = val.doubles();
if (temp.size() != 1) {
edm::LogError("MuonGeom") << "MuonGeometryNumbering:: ERROR: I need only 1 " << s << " in DDLogicalPart "
<< part.name();
abort();
}
return int(temp[0]);
} else
return 0;
}
|