|
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
|
#include "Geometry/MuonNumbering/interface/MuonOffsetFromDD.h"
#include "CondFormats/GeometryObjects/interface/MuonOffsetMap.h"
#include "DetectorDescription/Core/interface/DDFilteredView.h"
#include "DetectorDescription/Core/interface/DDFilter.h"
#include "DetectorDescription/Core/interface/DDValue.h"
#include "DetectorDescription/Core/interface/DDutils.h"
#include "DetectorDescription/DDCMS/interface/DDFilteredView.h"
#include <iostream>
#include <iomanip>
#include <sstream>
//#define EDM_ML_DEBUG
MuonOffsetFromDD::MuonOffsetFromDD(std::vector<std::string> name) : specpars_(name), nset_(name.size()) {
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "MuonOffsetFromDD initialized with " << nset_ << " specpars";
std::ostringstream st1;
for (const auto& name : specpars_)
st1 << " " << name;
edm::LogVerbatim("MuonGeom") << "MuonOffsetFromDD::SpecPars:" << st1.str();
#endif
}
bool MuonOffsetFromDD::build(const DDCompactView* cpv, MuonOffsetMap& php) {
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "Inside MuonOffsetFromDD::build(const DDCompactView*, MuonOffsetMap&)";
#endif
// Loop over all the sets
std::string attribute = "OnlyForMuonNumbering";
std::string name;
for (unsigned int k = 0; k < nset_; ++k) {
name = "muonstep" + std::to_string(k);
DDSpecificsMatchesValueFilter filter{DDValue(attribute, name, 0)};
DDFilteredView fv(*cpv, filter);
bool dodet = fv.firstChild();
DDsvalues_type sv(fv.mergedSpecifics());
int offset = getNumber("CopyNoOffset", sv);
int tag = getNumber("CopyNoTag", sv);
while (dodet) {
name = fv.logicalPart().name().name();
php.muonMap_[name] = std::make_pair(offset, tag);
dodet = fv.next();
}
}
return this->debugParameters(php);
}
bool MuonOffsetFromDD::build(const cms::DDCompactView* cpv, MuonOffsetMap& php) {
edm::LogVerbatim("MuonGeom") << "Inside MuonOffsetFromDD::build(const cms::DDCompactView*, MuonOffsetMap&)";
// Get the offsets and tags first
int offsets[nset_], tags[nset_];
cms::DDFilteredView fv(cpv->detector(), cpv->detector()->worldVolume());
for (unsigned int k = 0; k < nset_; ++k) {
std::vector<int> off = fv.get<std::vector<int>>(specpars_[k], "CopyNoOffset");
offsets[k] = (!off.empty()) ? off[0] : 0;
std::vector<int> tag = fv.get<std::vector<int>>(specpars_[k], "CopyNoTag");
tags[k] = (!tag.empty()) ? tag[0] : 0;
}
// Now loop over the detectors
std::string attribute = "OnlyForMuonNumbering";
std::string name;
for (unsigned int k = 0; k < nset_; ++k) {
name = "muonstep" + std::to_string(k);
const cms::DDFilter filter(attribute, name);
cms::DDFilteredView fv((*cpv), filter);
while (fv.firstChild()) {
name = static_cast<std::string>(fv.name());
php.muonMap_[name] = std::make_pair(offsets[k], tags[k]);
}
}
return this->debugParameters(php);
}
bool MuonOffsetFromDD::debugParameters(const MuonOffsetMap& php) {
edm::LogVerbatim("MuonGeom") << "MuonOffsetFromDD: Finds " << php.muonMap_.size() << " entries in the map";
#ifdef EDM_ML_DEBUG
unsigned int k(0);
for (auto itr = php.muonMap_.begin(); itr != php.muonMap_.end(); ++itr, ++k) {
edm::LogVerbatim("MuonGeom") << "[" << k << "] " << itr->first << ": (" << (itr->second).first << ", "
<< (itr->second).second << ")";
}
#endif
return true;
}
int MuonOffsetFromDD::getNumber(const std::string& str, const DDsvalues_type& sv) {
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "MuonOffsetFromDD::getNumbers called for " << str;
#endif
DDValue value(str);
if (DDfetch(&sv, value)) {
const std::vector<double>& fvec = value.doubles();
int nval = (!fvec.empty()) ? static_cast<int>(fvec[0]) : 0;
return nval;
} else {
return 0;
}
}
|