24 #include <unordered_set> 54 std::map<int, HGCalGeomParameters::layerParameters>
layers;
55 std::vector<HGCalParameters::hgtrform> trforms;
56 std::vector<bool> trformUse;
62 int nsiz =
static_cast<int>(
copy.size());
63 int lay = (nsiz > 0) ?
copy[nsiz - 1] : 0;
64 int zp = (nsiz > 2) ?
copy[nsiz - 3] : -1;
68 throw cms::Exception(
"DDException") <<
"Funny layer # " << lay <<
" zp " << zp <<
" in " << nsiz <<
" components";
71 php.
layer_.emplace_back(lay);
72 auto itr =
layers.find(lay);
74 double rin(0), rout(0);
78 const std::vector<double>& rmin = polyhedra.
rMinVec();
79 const std::vector<double>& rmax = polyhedra.
rMaxVec();
92 const CLHEP::HepRep3x3
rotation(
x.X(),
y.X(),
z.X(),
x.Y(),
y.Y(),
z.Y(),
x.Z(),
y.Z(),
z.Z());
93 const CLHEP::HepRotation hr(
rotation);
101 const CLHEP::Hep3Vector h3v(
xx,
yy,
zz);
109 trforms.emplace_back(mytrf);
110 trformUse.emplace_back(
false);
119 std::unordered_map<int32_t, int32_t> copies;
121 std::vector<int32_t> wafer2copy;
122 std::vector<HGCalGeomParameters::cellParameters> wafers;
124 DDValue val1(attribute, sdTag2, 0.0);
129 throw cms::Exception(
"DDException") <<
"Attribute " << val1 <<
" not found but needed.";
132 std::unordered_set<std::string>
names;
137 int nsiz =
static_cast<int>(
copy.size());
138 int wafer = (nsiz > 0) ?
copy[nsiz - 1] : 0;
139 int layer = (nsiz > 1) ?
copy[nsiz - 2] : 0;
141 throw cms::Exception(
"DDException") <<
"Funny wafer # " << wafer <<
" in " << nsiz <<
" components";
142 }
else if (
layer > static_cast<int>(
layers.size())) {
144 <<
" among " << nsiz <<
" components";
146 auto itr = copies.find(wafer);
147 auto cpy = copiesInLayers[
layer].find(wafer);
148 if (itr != copies.end() && cpy == copiesInLayers[
layer].end()) {
149 copiesInLayers[
layer][wafer] = itr->second;
151 if (itr == copies.end()) {
152 copies[wafer] = wafer2copy.size();
153 copiesInLayers[
layer][wafer] = wafer2copy.size();
160 wafer2copy.emplace_back(wafer);
163 wafers.emplace_back(cell);
165 std::vector<double>
zv, rv;
200 std::map<int, int> wafertype;
201 std::map<int, HGCalGeomParameters::cellParameters> cellsf, cellsc;
202 DDValue val2(attribute, sdTag3, 0.0);
207 throw cms::Exception(
"DDException") <<
"Attribute " << val2 <<
" not found but needed.";
214 int nsiz =
static_cast<int>(
copy.size());
215 int cellx = (nsiz > 0) ?
copy[nsiz - 1] : 0;
216 int wafer = (nsiz > 1) ?
copy[nsiz - 2] : 0;
221 <<
"Funny cell # " << cell <<
" type " <<
type <<
" in " << nsiz <<
" components";
223 auto ktr = wafertype.find(wafer);
224 if (ktr == wafertype.end())
225 wafertype[wafer] =
type;
227 std::map<int, HGCalGeomParameters::cellParameters>::iterator itr;
230 itr = cellsf.find(cell);
231 newc = (itr == cellsf.end());
233 itr = cellsc.find(cell);
234 newc = (itr == cellsc.end());
238 bool half = (
name.find(
"Half") != std::string::npos);
251 edm::LogVerbatim(
"HGCalGeom") <<
"Wafer " << wafer <<
" Type " <<
type <<
" Cell " << cellx <<
" local " 252 <<
xx <<
":" <<
yy <<
" new " <<
p1 <<
":" <<
p2;
268 layers, trforms, trformUse, copies, copiesInLayers, wafer2copy, wafers, wafertype, cellsf, cellsc, php);
279 std::map<int, HGCalGeomParameters::layerParameters>
layers;
280 std::vector<HGCalParameters::hgtrform> trforms;
281 std::vector<bool> trformUse;
282 std::vector<std::pair<int, int> > trused;
285 const std::vector<double>& pars = fv.
parameters();
288 int nsiz =
static_cast<int>(
copy.size());
289 int lay = (nsiz > 0) ?
copy[0] : 0;
290 int zp = (nsiz > 2) ?
copy[2] : -1;
294 throw cms::Exception(
"DDException") <<
"Funny layer # " << lay <<
" zp " << zp <<
" in " << nsiz <<
" components";
297 php.
layer_.emplace_back(lay);
298 auto itr =
layers.find(lay);
300 if (itr ==
layers.end()) {
301 double rin(0), rout(0);
302 if (dd4hep::isA<dd4hep::Polyhedra>(fv.
solid())) {
305 }
else if (dd4hep::isA<dd4hep::Tube>(fv.
solid())) {
306 dd4hep::Tube tubeSeg(fv.
solid());
313 std::pair<int, int> layz(lay, zp);
314 if (
std::find(trused.begin(), trused.end(), layz) == trused.end()) {
315 trused.emplace_back(layz);
318 const CLHEP::HepRep3x3
rotation(
x.X(),
y.X(),
z.X(),
x.Y(),
y.Y(),
z.Y(),
x.Z(),
y.Z(),
z.Z());
319 const CLHEP::HepRotation hr(
rotation);
327 const CLHEP::Hep3Vector h3v(
xx,
yy,
zz);
335 trforms.emplace_back(mytrf);
336 trformUse.emplace_back(
false);
345 std::unordered_map<int32_t, int32_t> copies;
347 std::vector<int32_t> wafer2copy;
348 std::vector<HGCalGeomParameters::cellParameters> wafers;
353 throw cms::Exception(
"DDException") <<
"Attribute " << sdTag2 <<
" not found but needed.";
356 std::unordered_set<std::string>
names;
360 int nsiz =
static_cast<int>(
copy.size());
361 int wafer = (nsiz > 0) ?
copy[0] : 0;
364 throw cms::Exception(
"DDException") <<
"Funny wafer # " << wafer <<
" in " << nsiz <<
" components";
365 }
else if (
layer > static_cast<int>(
layers.size())) {
367 <<
" among " << nsiz <<
" components";
369 auto itr = copies.find(wafer);
370 auto cpy = copiesInLayers[
layer].find(wafer);
371 if (itr != copies.end() && cpy == copiesInLayers[
layer].end()) {
372 copiesInLayers[
layer][wafer] = itr->second;
374 if (itr == copies.end()) {
375 copies[wafer] = wafer2copy.size();
376 copiesInLayers[
layer][wafer] = wafer2copy.size();
383 wafer2copy.emplace_back(wafer);
386 wafers.emplace_back(cell);
389 const std::vector<double>& pars = fv1.
parameters();
424 std::map<int, int> wafertype;
425 std::map<int, HGCalGeomParameters::cellParameters> cellsf, cellsc;
430 throw cms::Exception(
"DDException") <<
"Attribute " << sdTag3 <<
" not found but needed.";
436 int nsiz =
static_cast<int>(
copy.size());
437 int cellx = (nsiz > 0) ?
copy[0] : 0;
438 int wafer = (nsiz > 1) ?
copy[1] : 0;
443 <<
"Funny cell # " << cell <<
" type " <<
type <<
" in " << nsiz <<
" components";
445 auto ktr = wafertype.find(wafer);
446 if (ktr == wafertype.end())
447 wafertype[wafer] =
type;
449 std::map<int, HGCalGeomParameters::cellParameters>::iterator itr;
452 itr = cellsf.find(cell);
453 newc = (itr == cellsf.end());
455 itr = cellsc.find(cell);
456 newc = (itr == cellsc.end());
460 bool half = (
name.find(
"Half") != std::string::npos);
473 edm::LogVerbatim(
"HGCalGeom") <<
"Wafer " << wafer <<
" Type " <<
type <<
" Cell " << cellx <<
" local " 474 <<
xx <<
":" <<
yy <<
" new " <<
p1 <<
":" <<
p2;
490 layers, trforms, trformUse, copies, copiesInLayers, wafer2copy, wafers, wafertype, cellsf, cellsc, php);
494 std::vector<HGCalParameters::hgtrform>& trforms,
495 std::vector<bool>& trformUse,
496 const std::unordered_map<int32_t, int32_t>& copies,
498 const std::vector<int32_t>& wafer2copy,
499 const std::vector<HGCalGeomParameters::cellParameters>& wafers,
500 const std::map<int, int>& wafertype,
501 const std::map<int, HGCalGeomParameters::cellParameters>& cellsf,
502 const std::map<int, HGCalGeomParameters::cellParameters>& cellsc,
504 if (((cellsf.size() + cellsc.size()) == 0) || (wafers.empty()) || (
layers.empty())) {
505 throw cms::Exception(
"DDException") <<
"HGCalGeomParameters: mismatch between geometry and specpar: cells " 506 << cellsf.size() <<
":" << cellsc.size() <<
" wafers " << wafers.size()
507 <<
" layers " <<
layers.size();
510 for (
unsigned int i = 0;
i <
layers.size(); ++
i) {
522 for (
unsigned int i = 0;
i < php.
layer_.size(); ++
i) {
523 for (
unsigned int i1 = 0;
i1 < trforms.size(); ++
i1) {
524 if (!trformUse[
i1] && php.
layerGroup_[trforms[
i1].lay - 1] == static_cast<int>(
i + 1)) {
526 trforms[
i1].lay = (
i + 1);
527 trformUse[
i1] =
true;
530 for (
unsigned int i2 =
i1 + 1;
i2 < trforms.size(); ++
i2) {
531 if (!trformUse[
i2] && trforms[
i2].zp == trforms[
i1].zp &&
532 php.
layerGroup_[trforms[
i2].lay - 1] == static_cast<int>(
i + 1)) {
535 trformUse[
i2] =
true;
546 for (
unsigned i = 0;
i < wafer2copy.size(); ++
i) {
550 auto ktr = wafertype.find(wafer2copy[
i]);
551 int typet = (ktr == wafertype.end()) ? 0 : (ktr->second);
553 double r = wafers[
i].xyz.perp();
555 for (
int k = 1;
k < 4; ++
k) {
566 std::vector<HGCalGeomParameters::cellParameters>::const_iterator itrf = wafers.end();
567 for (
unsigned int i = 0;
i < cellsf.size(); ++
i) {
568 auto itr = cellsf.find(
i);
569 if (itr == cellsf.end()) {
570 throw cms::Exception(
"DDException") <<
"HGCalGeomParameters: missing info for fine cell number " <<
i;
572 double xx = (itr->second).xyz.x();
573 double yy = (itr->second).xyz.y();
574 int waf = (itr->second).wafer;
582 for (
unsigned int i = 0;
i < cellsc.size(); ++
i) {
583 auto itr = cellsc.find(
i);
584 if (itr == cellsc.end()) {
585 throw cms::Exception(
"DDException") <<
"HGCalGeomParameters: missing info for coarse cell number " <<
i;
587 double xx = (itr->second).xyz.x();
588 double yy = (itr->second).xyz.y();
589 int waf = (itr->second).wafer;
602 php.
depth_.emplace_back(
i + 1);
642 edm::LogVerbatim(
"HGCalGeom") <<
"HGCalGeomParameters: wafer radius " << php.
waferR_ <<
" and dimensions of the " 656 <<
" coarse cells in a wafer";
675 for (std::unordered_map<int, int>::const_iterator itr = theModules.begin(); itr != theModules.end(); ++itr, ++k2) {
676 edm::LogVerbatim(
"HGCalGeom") <<
"[" << k2 <<
"] " << itr->first <<
":" << itr->second;
685 std::map<int, HGCalGeomParameters::layerParameters>
layers;
696 int nsiz =
static_cast<int>(
copy.size());
697 if (nsiz < levelTop) {
698 int lay =
copy[nsiz - 1];
704 edm::LogVerbatim(
"HGCalGeom") <<
sol.name() <<
" shape " <<
sol.shape() <<
" size " << nsiz <<
":" << levelTop
705 <<
" lay " << lay <<
" z " <<
zside;
709 <<
"Funny layer # " << lay <<
" zp " <<
zside <<
" in " << nsiz <<
" components";
712 php.
layer_.emplace_back(lay);
716 auto itr =
layers.find(lay);
717 if (itr ==
layers.end()) {
722 (itr->second).rmin =
std::min(rin, (itr->second).rmin);
723 (itr->second).rmax =
std::max(rout, (itr->second).rmax);
725 if (trforms.find(std::make_pair(lay,
zside)) == trforms.end()) {
728 const CLHEP::HepRep3x3
rotation(
x.X(),
y.X(),
z.X(),
x.Y(),
y.Y(),
z.Y(),
x.Z(),
y.Z(),
z.Z());
729 const CLHEP::HepRotation hr(
rotation);
744 trforms[std::make_pair(lay,
zside)] = mytrf;
762 std::map<int, HGCalGeomParameters::layerParameters>
layers;
773 int nsiz =
static_cast<int>(fv.
level());
774 if (nsiz < levelTop) {
782 <<
" size " << nsiz <<
":" << levelTop <<
" lay " << lay <<
" z " <<
zside <<
":" 787 <<
"Funny layer # " << lay <<
" zp " <<
zside <<
" in " << nsiz <<
" components";
790 php.
layer_.emplace_back(lay);
791 const std::vector<double>& pars = fv.
parameters();
794 auto itr =
layers.find(lay);
795 if (itr ==
layers.end()) {
800 (itr->second).rmin =
std::min(rin, (itr->second).rmin);
801 (itr->second).rmax =
std::max(rout, (itr->second).rmax);
803 if (trforms.find(std::make_pair(lay,
zside)) == trforms.end()) {
806 const CLHEP::HepRep3x3
rotation(
x.X(),
y.X(),
z.X(),
x.Y(),
y.Y(),
z.Y(),
x.Z(),
y.Z(),
z.Z());
807 const CLHEP::HepRotation hr(
rotation);
822 trforms[std::make_pair(lay,
zside)] = mytrf;
839 edm::LogVerbatim(
"HGCalGeom") <<
"HGCalGeomParameters (DDD)::loadGeometryHexagonModule called with tags " << sdTag1
840 <<
":" << sdTag2 <<
" firstLayer " << firstLayer <<
":" << php.
firstMixedLayer_;
841 int ntot1(0), ntot2(0);
843 std::map<int, HGCalGeomParameters::layerParameters>
layers;
844 std::map<std::pair<int, int>,
double> zvals;
849 DDValue val1(attribute, sdTag2, 0.0);
858 int nsiz =
static_cast<int>(
copy.size());
859 if (levelTop < nsiz) {
860 int lay =
copy[levelTop];
866 <<
"Funny layer # " << lay <<
" zp " <<
zside <<
" in " << nsiz <<
" components";
868 if (zvals.find(std::make_pair(lay,
zside)) == zvals.end()) {
871 std::ostringstream st1;
872 st1 <<
"Name0 " << fv1.
name() <<
" LTop " << levelTop <<
":" << lay <<
" ZSide " <<
zside <<
" # of levels " 874 for (
const auto&
c :
copy)
876 st1 <<
" Z " << zvals[std::make_pair(lay,
zside)];
885 DDValue val2(attribute, sdTag1, 0.0);
894 int nsiz =
static_cast<int>(
copy.size());
895 if (levelTop < nsiz) {
896 int lay =
copy[levelTop];
902 std::ostringstream st2;
903 st2 <<
"Name1 " <<
sol.name() <<
" shape " <<
sol.shape() <<
" LTop " << levelTop <<
":" << lay <<
" ZSide " 905 for (
const auto&
c :
copy)
911 <<
"Funny layer # " << lay <<
" zp " <<
zside <<
" in " << nsiz <<
" components";
913 if (zvals.find(std::make_pair(lay,
zside)) != zvals.end()) {
915 php.
layer_.emplace_back(lay);
916 auto itr =
layers.find(lay);
917 if (itr ==
layers.end()) {
923 double zp = zvals[std::make_pair(lay, 1)];
927 std::ostringstream st3;
928 st3 <<
"Name1 " << fv2.
name() <<
" LTop " << levelTop <<
":" << lay <<
" ZSide " <<
zside <<
" # of levels " 930 for (
const auto&
c :
copy)
932 st3 <<
" R " << rin <<
":" << rout;
937 if (trforms.find(std::make_pair(lay,
zside)) == trforms.end()) {
940 const CLHEP::HepRep3x3
rotation(
x.X(),
y.X(),
z.X(),
x.Y(),
y.Y(),
z.Y(),
x.Z(),
y.Z(),
z.Z());
941 const CLHEP::HepRotation hr(
rotation);
948 const CLHEP::Hep3Vector h3v(
xx,
yy, zvals[std::make_pair(lay,
zside)]);
956 trforms[std::make_pair(lay,
zside)] = mytrf;
967 edm::LogVerbatim(
"HGCalGeom") <<
"Total # of views " << ntot1 <<
":" << ntot2;
978 edm::LogVerbatim(
"HGCalGeom") <<
"HGCalGeomParameters (DD4hep)::loadGeometryHexagonModule called with tags " << sdTag1
979 <<
":" << sdTag2 <<
" firstLayer " << firstLayer;
980 int ntot1(0), ntot2(0);
982 std::map<int, HGCalGeomParameters::layerParameters>
layers;
984 std::map<std::pair<int, int>,
double> zvals;
993 int nsiz =
static_cast<int>(fv1.
level());
994 if (nsiz > levelTop) {
996 int lay =
copy[nsiz - levelTop - 1];
1002 <<
"Funny layer # " << lay <<
" zp " <<
zside <<
" in " << nsiz <<
" components";
1004 if (zvals.find(std::make_pair(lay,
zside)) == zvals.end()) {
1007 std::ostringstream st1;
1008 st1 <<
"Name0 " << fv1.
name() <<
" LTop " << levelTop <<
":" << lay <<
" ZSide " <<
zside <<
" # of levels " 1010 for (
const auto&
c :
copy)
1012 st1 <<
" Z " << zvals[std::make_pair(lay,
zside)];
1024 int nsiz =
static_cast<int>(fv2.
level());
1028 if (nsiz > levelTop) {
1030 int lay =
copy[nsiz - levelTop - 1];
1035 std::ostringstream st2;
1037 << levelTop <<
":" << lay <<
" ZSide " <<
zside <<
":" << php.
levelZSide_ <<
" # of levels " << nsiz;
1044 <<
"Funny layer # " << lay <<
" zp " <<
zside <<
" in " << nsiz <<
" components";
1046 if (zvals.find(std::make_pair(lay,
zside)) != zvals.end()) {
1048 php.
layer_.emplace_back(lay);
1049 auto itr =
layers.find(lay);
1050 if (itr ==
layers.end()) {
1051 const std::vector<double>& pars = fv2.
parameters();
1056 double zp = zvals[std::make_pair(lay, 1)];
1060 std::ostringstream st3;
1061 st3 <<
"Name2 " << fv2.
name() <<
" LTop " << levelTop <<
":" << lay <<
" ZSide " <<
zside <<
" # of levels " 1063 for (
const auto&
c :
copy)
1065 st3 <<
" R " << rin <<
":" << rout;
1070 if (trforms.find(std::make_pair(lay,
zside)) == trforms.end()) {
1073 const CLHEP::HepRep3x3
rotation(
x.X(),
y.X(),
z.X(),
x.Y(),
y.Y(),
z.Y(),
x.Z(),
y.Z(),
z.Z());
1074 const CLHEP::HepRotation hr(
rotation);
1081 const CLHEP::Hep3Vector h3v(
xx,
yy, zvals[std::make_pair(lay,
zside)]);
1089 trforms[std::make_pair(lay,
zside)] = mytrf;
1099 edm::LogVerbatim(
"HGCalGeom") <<
"Total # of views " << ntot1 <<
":" << ntot2;
1106 const int& firstLayer,
1108 double rmin(0), rmax(0);
1109 for (
unsigned int i = 0;
i <
layers.size(); ++
i) {
1111 if (
layer.first == static_cast<int>(
i + firstLayer)) {
1117 rmin =
layer.second.rmin;
1118 rmax =
layer.second.rmax;
1120 if (rmin >
layer.second.rmin)
1121 rmin =
layer.second.rmin;
1122 if (rmax <
layer.second.rmax)
1123 rmax =
layer.second.rmax;
1129 php.
rLimit_.emplace_back(rmin);
1130 php.
rLimit_.emplace_back(rmax);
1135 for (
unsigned int i = 0;
i < php.
layer_.size(); ++
i) {
1136 for (
auto& trform : trforms) {
1137 if (trform.first.first == static_cast<int>(
i + firstLayer)) {
1183 if (!dummy2.empty())
1192 if (fv1.firstChild()) {
1194 const auto&
dummy = getDDDArray(
"WaferSize",
sv, 0);
1195 waferSize_ =
dummy[0];
1201 if (fv2.firstChild()) {
1203 php.cellSize_ = getDDDArray(
"CellSize",
sv, 0);
1206 loadSpecParsHexagon(php);
1215 php.
boundR_ = fv.
get<std::vector<double> >(sdTag4,
"RadiusBound");
1217 php.
rLimit_ = fv.
get<std::vector<double> >(sdTag4,
"RadiusLimits");
1225 php.
slopeMin_ = fv.
get<std::vector<double> >(sdTag4,
"Slope");
1230 const auto&
dummy = fv.
get<std::vector<double> >(sdTag2,
"WaferSize");
1234 php.
cellSize_ = fv.
get<std::vector<double> >(sdTag3,
"CellSize");
1238 const auto& dummy2 = fv.
get<std::vector<double> >(sdTag1,
"LayerOffset");
1239 if (!dummy2.empty()) {
1250 edm::LogVerbatim(
"HGCalGeom") <<
"HGCalGeomParameters: wafer radius ranges" 1251 <<
" for cell grouping " << php.
boundR_[0] <<
":" << php.
boundR_[1] <<
":" 1256 <<
"for the 3 ranges:";
1262 for (
unsigned int k = 0;
k < php.
cellSize_.size(); ++
k)
1316 std::vector<int> layerType,
waferIndex, waferProperties;
1317 std::vector<double> cassetteShift;
1324 cassetteShift = fv.
vector(
"CassetteShiftEE");
1331 cassetteShift = fv.
vector(
"CassetteShiftHE");
1364 const auto&
dummy = fv.
get<std::vector<double> >(sdTag1,
"RadiusCuts");
1365 if (
dummy.size() > 3) {
1374 php.
slopeMin_ = fv.
get<std::vector<double> >(sdTag1,
"SlopeBottom");
1375 php.
zFrontMin_ = fv.
get<std::vector<double> >(sdTag1,
"ZFrontBottom");
1377 php.
rMinFront_ = fv.
get<std::vector<double> >(sdTag1,
"RMinFront");
1380 php.
slopeTop_ = fv.
get<std::vector<double> >(sdTag1,
"SlopeTop");
1381 php.
zFrontTop_ = fv.
get<std::vector<double> >(sdTag1,
"ZFrontTop");
1383 php.
rMaxFront_ = fv.
get<std::vector<double> >(sdTag1,
"RMaxFront");
1386 for (
unsigned int k = 0;
k < kmax; ++
k)
1389 const auto& dummy2 = fv.
get<std::vector<double> >(sdTag1,
"LayerOffset");
1390 if (!dummy2.empty()) {
1396 for (
auto const& it : vmap) {
1397 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"RadiusMixBoundary")) {
1398 for (
const auto&
i : it.second)
1400 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"ZRanges")) {
1401 for (
const auto&
i : it.second)
1403 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"LayerCenter")) {
1404 for (
const auto&
i : it.second)
1413 std::vector<int> layerType,
waferIndex, waferProperties;
1414 std::vector<double> cassetteShift;
1416 for (
auto const& it : vmap) {
1417 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"WaferIndexEE")) {
1418 for (
const auto&
i : it.second)
1420 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"WaferPropertiesEE")) {
1421 for (
const auto&
i : it.second)
1422 waferProperties.emplace_back(std::round(
i));
1426 for (
auto const& it : vmap) {
1427 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"CassetteShiftEE")) {
1428 for (
const auto&
i : it.second)
1429 cassetteShift.emplace_back(
i);
1434 for (
auto const& it : vmap) {
1435 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"WaferIndexHE")) {
1436 for (
const auto&
i : it.second)
1438 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"WaferPropertiesHE")) {
1439 for (
const auto&
i : it.second)
1440 waferProperties.emplace_back(std::round(
i));
1444 for (
auto const& it : vmap) {
1445 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"CassetteShiftHE")) {
1446 for (
const auto&
i : it.second)
1447 cassetteShift.emplace_back(
i);
1454 for (
auto const& it : vmap) {
1455 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"LayerTypesEE")) {
1456 for (
const auto&
i : it.second)
1457 layerType.emplace_back(std::round(
i));
1461 for (
auto const& it : vmap) {
1462 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"LayerTypesHE")) {
1463 for (
const auto&
i : it.second)
1464 layerType.emplace_back(std::round(
i));
1483 <<
"parameters for 120 to 200 micron " 1488 <<
"parameters for 200 to 300 micron " 1513 const std::vector<int>& layerType,
1515 const std::vector<int>& waferProperties) {
1517 for (
unsigned int k = 0;
k < layerType.size(); ++
k) {
1526 php.
layerRotV_.emplace_back(std::make_pair(cth, sth));
1592 std::vector<int> tileHEX1, tileHEX2, tileHEX3, tileHEX4;
1593 std::vector<double> tileRMin, tileRMax;
1594 std::vector<int> tileRingMin, tileRingMax;
1595 std::vector<double> cassetteShift;
1602 tileRMin = fv.
vector(
"TileRMin");
1603 tileRMax = fv.
vector(
"TileRMax");
1611 cassetteShift = fv.
vector(
"CassetteShiftHE");
1634 for (
auto const& it : vmap) {
1635 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"RadiusMixBoundary")) {
1636 for (
const auto&
i : it.second)
1638 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"ZRanges")) {
1639 for (
const auto&
i : it.second)
1641 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"LayerCenter")) {
1642 for (
const auto&
i : it.second)
1649 php.
rMinLayerBH_ = fv.
get<std::vector<double> >(sdTag1,
"RMinLayerBH");
1656 php.
slopeMin_ = fv.
get<std::vector<double> >(sdTag1,
"SlopeBottom");
1657 php.
zFrontMin_ = fv.
get<std::vector<double> >(sdTag1,
"ZFrontBottom");
1659 php.
rMinFront_ = fv.
get<std::vector<double> >(sdTag1,
"RMinFront");
1662 php.
slopeTop_ = fv.
get<std::vector<double> >(sdTag1,
"SlopeTop");
1663 php.
zFrontTop_ = fv.
get<std::vector<double> >(sdTag1,
"ZFrontTop");
1665 php.
rMaxFront_ = fv.
get<std::vector<double> >(sdTag1,
"RMaxFront");
1668 for (
unsigned int k = 0;
k < kmax; ++
k)
1671 const auto& dummy2 = fv.
get<std::vector<double> >(sdTag1,
"LayerOffset");
1679 std::vector<int> tileHEX1, tileHEX2, tileHEX3, tileHEX4;
1680 std::vector<double> tileRMin, tileRMax;
1681 std::vector<int> tileRingMin, tileRingMax;
1682 std::vector<double> cassetteShift;
1683 for (
auto const& it : vmap) {
1684 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"TileIndex")) {
1685 for (
const auto&
i : it.second)
1686 tileIndx.emplace_back(std::round(
i));
1687 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"TileProperty")) {
1688 for (
const auto&
i : it.second)
1690 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"TileHEX1")) {
1691 for (
const auto&
i : it.second)
1692 tileHEX1.emplace_back(std::round(
i));
1693 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"TileHEX2")) {
1694 for (
const auto&
i : it.second)
1695 tileHEX2.emplace_back(std::round(
i));
1696 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"TileHEX3")) {
1697 for (
const auto&
i : it.second)
1698 tileHEX3.emplace_back(std::round(
i));
1699 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"TileHEX4")) {
1700 for (
const auto&
i : it.second)
1701 tileHEX4.emplace_back(std::round(
i));
1702 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"TileRMin")) {
1703 for (
const auto&
i : it.second)
1705 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"TileRMax")) {
1706 for (
const auto&
i : it.second)
1708 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"TileRingMin")) {
1709 for (
const auto&
i : it.second)
1710 tileRingMin.emplace_back(std::round(
i));
1711 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"TileRingMax")) {
1712 for (
const auto&
i : it.second)
1713 tileRingMax.emplace_back(std::round(
i));
1717 for (
auto const& it : vmap) {
1718 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"CassetteShiftHE")) {
1719 for (
const auto&
i : it.second)
1720 cassetteShift.emplace_back(
i);
1773 const std::vector<int>& tileIndx,
1775 const std::vector<int>& tileHEX1,
1776 const std::vector<int>& tileHEX2,
1777 const std::vector<int>& tileHEX3,
1778 const std::vector<int>& tileHEX4,
1779 const std::vector<double>& tileRMin,
1780 const std::vector<double>& tileRMax,
1781 const std::vector<int>& tileRingMin,
1782 const std::vector<int>& tileRingMax) {
1784 for (
unsigned int k = 0;
k < tileIndx.size(); ++
k) {
1795 <<
":" << tileHEX2[
k] <<
":" << tileHEX3[
k] <<
":" << tileHEX4[
k] <<
std::dec;
1799 for (
unsigned int k = 0;
k < tileRMin.size(); ++
k) {
1802 edm::LogVerbatim(
"HGCalGeom") <<
"TileRingR[" <<
k <<
"] " << tileRMin[
k] <<
":" << tileRMax[
k];
1806 for (
unsigned k = 0;
k < tileRingMin.size(); ++
k) {
1809 edm::LogVerbatim(
"HGCalGeom") <<
"TileRingRange[" <<
k <<
"] " << tileRingMin[
k] <<
":" << tileRingMax[
k];
1818 edm::LogVerbatim(
"HGCalGeom") <<
"Input waferWidth " << waferW <<
":" << rmin <<
" R Limits: " << rin <<
":" << rout
1819 <<
" Fine " << rMaxFine;
1827 double dx = 0.5 * waferW;
1828 double dy = 3.0 *
dx *
tan(30._deg);
1829 double rr = 2.0 *
dx *
tan(30._deg);
1830 int ncol =
static_cast<int>(2.0 * rout / waferW) + 1;
1831 int nrow =
static_cast<int>(rout / (waferW *
tan(30._deg))) + 1;
1832 int ns2 = (2 * ncol + 1) * (2 * nrow + 1) * php.
layer_.size();
1833 int incm(0), inrm(0), kount(0),
ntot(0);
1839 for (
int nr = -nrow;
nr <= nrow; ++
nr) {
1840 int inr = (
nr >= 0) ?
nr : -
nr;
1841 for (
int nc = -ncol; nc <= ncol; ++nc) {
1842 int inc = (nc >= 0) ? nc : -nc;
1843 if (inr % 2 == inc % 2) {
1844 double xpos = nc *
dx;
1845 double ypos =
nr *
dy;
1847 double rpos =
std::sqrt(xpos * xpos + ypos * ypos);
1848 int typet = (rpos < rMaxFine) ? 1 : 2;
1850 for (
int k = 1;
k < 4; ++
k) {
1851 if ((rpos + rmin) <= php.
boundR_[
k]) {
1866 <<
" Location " <<
corner.first <<
" Position " << xpos <<
":" << ypos
1867 <<
" Layers " << php.
layer_.size();
1874 for (
unsigned int il = 0; il < php.
layer_.size(); ++il) {
1875 std::pair<int, int>
corner =
1878 auto cpy = copiesInLayers[php.
layer_[il]].find(
copy);
1879 if (cpy == copiesInLayers[php.
layer_[il]].end())
1897 edm::LogVerbatim(
"HGCalGeom") <<
"HGCalWaferHexagon: # of columns " << incm <<
" # of rows " << inrm <<
" and " 1898 << kount <<
":" <<
ntot <<
" wafers; R " << rin <<
":" << rout;
1900 for (
unsigned int k = 0;
k < copiesInLayers.size(); ++
k) {
1901 const auto& theModules = copiesInLayers[
k];
1904 for (std::unordered_map<int, int>::const_iterator itr = theModules.begin(); itr != theModules.end(); ++itr, ++k2) {
1905 edm::LogVerbatim(
"HGCalGeom") <<
"[" << k2 <<
"] " << itr->first <<
":" << itr->second;
1924 edm::LogVerbatim(
"HGCalGeom") <<
"Input waferWidth " << waferW <<
":" << waferS <<
" R Max: " << rout;
1932 double r = 0.5 * (waferW + waferS);
1934 double dy = 0.75 *
R;
1935 double r1 = 0.5 * waferW;
1937 int N = (
r == 0) ? 3 : (static_cast<int>(0.5 * rout /
r) + 3);
1938 int ns1 = (2 *
N + 1) * (2 *
N + 1);
1941 edm::LogVerbatim(
"HGCalGeom") <<
"wafer " << waferW <<
":" << waferS <<
" r " <<
r <<
" dy " <<
dy <<
" N " <<
N 1942 <<
" sizes " << ns1 <<
":" << ns2;
1943 std::vector<int> indtypes(ns1 + 1);
1949 int ipos(0), lpos(0), uvmax(0), nwarn(0);
1950 std::vector<int> uvmx(php.
zLayerHex_.size(), 0);
1951 for (
int v = -
N;
v <=
N; ++
v) {
1952 for (
int u = -
N; u <=
N; ++u) {
1954 int nc = -2 * u +
v;
1955 double xpos = nc *
r;
1956 double ypos =
nr *
dy;
1961 wafersInLayers[indx] = ipos;
1972 double xpos0 = xpos + xyoff.first;
1973 double ypos0 = ypos + xyoff.second;
1985 typesInLayers[kndx] = lpos;
1988 indtypes.emplace_back(kndx);
1990 std::pair<int, int>
corner =
1995 <<
" u " << u <<
" v " <<
v <<
" with " <<
corner.first <<
" corners";
2004 edm::LogVerbatim(
"HGCalGeom") <<
"Layer " << lay <<
" u|v " << u <<
":" <<
v <<
" with corner " 2025 <<
"Layer:u:v " <<
i <<
":" << lay <<
":" << u <<
":" <<
v <<
" Part " << corner0.first <<
":" 2026 <<
part <<
" Orient " << corner0.second <<
":" <<
orient <<
" Position " << xpos0 <<
":" << ypos0
2028 <<
" Compatibiliety Flag " <<
ok;
2034 waferTypes[
wl] = corner0;
2037 <<
"Layer " << lay <<
" u|v " << u <<
":" <<
v <<
" Index " << std::hex <<
wl <<
std::dec <<
" pos " 2038 << xpos0 <<
":" << ypos0 <<
" R " <<
r1 <<
":" << R1 <<
" Range " << php.
rMinLayHex_[
i] <<
":" 2040 <<
corner.second <<
" croner0 " << corner0.first <<
":" << corner0.second;
2045 edm::LogVerbatim(
"HGCalGeom") <<
"Layer " << lay <<
" u|v " << u <<
":" <<
v <<
" with corner " 2054 edm::LogWarning(
"HGCalGeom") <<
"HGCalGeomParameters::loadWafer8: there are " << nwarn
2055 <<
" wafers with non-matching partial- orientation types";
2093 int id = indtypes[
k];
2114 << xyoff.first <<
":" << xyoff.second;
2135 for (
auto const& it : vmap) {
2136 if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"waferFine")) {
2137 for (
const auto&
i : it.second)
2139 }
else if (dd4hep::dd::compareEqual(dd4hep::dd::noNamespace(it.first),
"waferCoarse")) {
2140 for (
const auto&
i : it.second)
2151 for (
unsigned int k = 0;
k < php.
cellFine_.size(); ++
k)
2173 for (
unsigned int k = 0;
k < 2; ++
k) {
2198 for (
unsigned int k = 0;
k < 2; ++
k) {
2208 <<
" dEta = " << php.
cellSize_[
k] <<
"\nOld[0] new R = " << rv <<
" Eta = " <<
eta;
2216 edm::LogVerbatim(
"HGCalGeom") <<
"Old [" << kount <<
"] new R = " << rv <<
" Eta = " <<
eta;
2224 std::vector<double>::iterator
low,
high;
2241 edm::LogVerbatim(
"HGCalGeom") <<
"Modified irlow " << irlow <<
" dr " << drlow;
2260 edm::LogVerbatim(
"HGCalGeom") <<
"Modified irhigh " << irhigh <<
" dr " << drhigh;
2266 edm::LogVerbatim(
"HGCalGeom") <<
"Old Layer " <<
k <<
" Type " <<
kk <<
" Low edge " << irlow <<
":" << drlow
2267 <<
" Top edge " << irhigh <<
":" << drhigh;
2272 for (
unsigned int k = 0;
k < 2; ++
k) {
2301 mytr.
h = 0.5 * (rmax - rmin);
2330 const std::vector<double>& fvec =
value.doubles();
2331 int nval = fvec.size();
2335 <<
"HGCalGeomParameters: # of " <<
str <<
" bins " << nval <<
" < " <<
nmin <<
" ==> illegal";
2338 if (nval < 1 &&
nmin == 0) {
2340 <<
"HGCalGeomParameters: # of " <<
str <<
" bins " << nval <<
" < 1 ==> illegal" 2341 <<
" (nmin=" <<
nmin <<
")";
2347 throw cms::Exception(
"DDException") <<
"HGCalGeomParameters: cannot get array " <<
str;
2349 std::vector<double> fvec;
2355 const std::vector<HGCalGeomParameters::cellParameters>& wafers,
2356 std::vector<HGCalGeomParameters::cellParameters>::const_iterator& itrf,
2360 if (itrf == wafers.end()) {
2361 for (std::vector<HGCalGeomParameters::cellParameters>::const_iterator itr = wafers.begin(); itr != wafers.end();
2363 if (itr->wafer == wafer) {
2369 double dx(0),
dy(0);
2370 if (itrf != wafers.end()) {
2371 dx = (
xx - itrf->xyz.x());
2374 dy = (
yy - itrf->xyz.y());
2378 return std::make_pair(
dx,
dy);
2382 std::for_each(
v.begin(),
v.end(), [
s](
double&
n) {
n *=
s; });
std::vector< int > iradMaxBH_
std::vector< double > waferPosY_
Log< level::Info, true > LogVerbatim
std::vector< int > layer_
static constexpr int scintillatorCassette
std::vector< double > moduleDzR_
std::vector< int > depthLayerF_
std::vector< int > depth_
std::vector< double > zFrontMin_
int32_t tileType(const int32_t property)
HGCalGeomTools geomTools_
nav_type copyNumbers() const
return the stack of copy numbers
hgtrap getModule(unsigned int k, bool reco) const
std::vector< double > moduleHR_
std::vector< double > rMinVec(void) const
layer_map copiesInLayers_
static constexpr int k_OffsetRotation
std::vector< double > const & vector(std::string_view iKey) const
The DDVector information.
std::vector< std::pair< double, double > > layerRotV_
std::vector< bool > cellCoarseHalf_
int scintType(const int layer) const
std::vector< bool > cellFineHalf_
const std::vector< int > copyNos() const
The list of the volume copy numbers.
const cms::DDSolidShape shape() const
std::vector< int > moduleLayR_
static constexpr int siliconFileEE
int32_t waferU(const int32_t index)
std::vector< int > cellFine_
void rescale(std::vector< double > &, const double s)
for(int i=first, nt=offsets[nh];i< nt;i+=gridDim.x *blockDim.x)
int32_t waferLayer(const int32_t index)
static int32_t getUnpackedCell6(int id)
std::vector< double > moduleHS_
std::string name(Mapping a, V value)
std::vector< double > xVec(void) const
std::vector< double > const & vector(std::string_view iKey) const
returns an empty container if not found
std::vector< double > zVec(void) const
Sin< T >::type sin(const T &t)
std::vector< double > trformTranY_
void loadGeometryHexagonModule(const DDCompactView *cpv, HGCalParameters &php, const std::string &sdTag1, const std::string &sdTag2, int firstLayer)
Global3DPoint GlobalPoint
std::vector< double > rMaxVec(void) const
std::vector< double > cellFineY_
void loadSpecParsHexagon(const DDFilteredView &fv, HGCalParameters &php, const DDCompactView *cpv, const std::string &sdTag1, const std::string &sdTag2)
const std::string & name() const
The name of a logical-part of the current node in the filtered-view.
std::vector< double > trformRotZY_
static int getRotation(int zside, int type, int rotn)
void resetZero(std::vector< double > &)
std::vector< uint32_t > trformIndex_
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
std::vector< int > layerGroupM_
Compact representation of the geometrical detector hierarchy.
bool DDfetch(const DDsvalues_type *, DDValue &)
helper for retrieving DDValues from DDsvalues_type *.
HGCalGeometryMode::GeometryMode mode_
const std::string names[nVars_]
int32_t waferOrient(const int32_t property)
std::vector< double > trformRotXX_
std::vector< int > nPhiBinBH_
std::unordered_map< std::string, std::vector< double > > DDVectorsMap
A DDSolid represents the shape of a part.
static constexpr double tol
void fillTrForm(const hgtrform &mytr)
void loadSpecParsHexagon8(const DDFilteredView &fv, HGCalParameters &php)
static constexpr uint32_t k_CornerSize
static bool goodTypeMode(double xpos, double ypos, double delX, double delY, double rin, double rout, int part, int rotn, bool debug=false)
constexpr std::array< uint8_t, layerIndexSize > layer
std::vector< double > zVec(void) const
wafer_map wafersInLayers_
static int32_t getUnpackedCellType6(int id)
std::vector< double > rMinLayerBH_
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DD3Vector
std::vector< double > trformRotZX_
double scintCellSize(const int layer) const
std::string_view name() const
std::vector< double > cellCoarseX_
static constexpr int scintillatorFile
int32_t waferCassette(const int32_t property)
std::vector< std::pair< int, int > > tileRingRange_
static constexpr int siliconCassetteHE
std::vector< int > firstModule_
std::vector< int > cellCoarse_
std::vector< double > trformRotYZ_
std::vector< std::pair< unsigned int, DDValue > > DDsvalues_type
void loadCellTrapezoid(HGCalParameters &php)
std::pair< double, double > cellPosition(const std::vector< cellParameters > &wafers, std::vector< cellParameters >::const_iterator &itrf, int wafer, double xx, double yy)
std::vector< double > boundR_
const std::string & name() const
Returns the name.
std::vector< double > cellSize_
std::vector< int > waferUVMaxLayer_
std::vector< double > moduleDzS_
bool next()
set current node to the next node in the filtered tree
std::vector< int > layerIndex_
std::vector< double > moduleAlphaR_
std::vector< double > yLayerHex_
void loadCellParsHexagon(const DDCompactView *cpv, HGCalParameters &php)
std::vector< double > trformRotXY_
std::vector< double > rMaxFront_
Cos< T >::type cos(const T &t)
std::vector< double > trformRotYX_
static constexpr int siliconCassetteEE
Tan< T >::type tan(const T &t)
Abs< T >::type abs(const T &t)
std::vector< double > getDDDArray(const std::string &str, const DDsvalues_type &sv, const int nmin)
const int level() const
get Iterator level
int32_t tileSiPM(const int32_t property)
std::vector< double > waferThickness_
std::vector< double > slopeTop_
std::vector< int > layerCenter_
static constexpr double k_ScaleFromDD4hep
std::vector< double > moduleBlR_
int32_t waferThick(const int32_t property)
const std::array< const cms::dd::NameValuePair< DDSolidShape >, 21 > DDSolidShapeMap
std::vector< double > rMinLayHex_
Quality *__restrict__ uint16_t nmin
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
void fillModule(const hgtrap &mytr, bool reco)
std::vector< double > moduleTlS_
bool firstChild()
set the current node to the first child
std::vector< double > zLayerHex_
std::vector< int > layerType_
std::vector< int > dbl_to_int(const std::vector< double > &vecdbl)
Converts a std::vector of doubles to a std::vector of int.
static constexpr double tan30deg_
void loadWaferHexagon(HGCalParameters &php)
DDsvalues_type mergedSpecifics() const
std::vector< double > rMaxLayHex_
std::vector< double > trformTranX_
static constexpr double k_ScaleFromDD4hepToG4
std::unordered_map< int32_t, std::pair< int32_t, int32_t > > waferT_map
const DDLogicalPart & logicalPart() const
The logical-part of the current node in the filtered-view.
std::vector< double > zRanges_
std::vector< double > slopeMin_
static constexpr int32_t WaferCenterR
T get(const std::string &)
extract attribute value
std::vector< int > lastModule_
static constexpr double k_ScaleToDDD
std::vector< double > radiusMixBoundary_
std::vector< double > cellThickness_
std::vector< double > trformRotXZ_
ROOT::Math::PositionVector3D< ROOT::Math::Cartesian3D< double > > XYZPointD
point in space with cartesian internal representation
std::vector< int > layerGroup_
std::unordered_map< int32_t, int32_t > wafer_map
static std::pair< int, int > getTypeMode(const double &xpos, const double &ypos, const double &delX, const double &delY, const double &rin, const double &rout, const int &waferType, const int &mode, bool debug=false)
std::vector< double > radius200to300_
std::vector< double > radius100to200_
const DDSolid & solid(void) const
Returns a reference object of the solid being the shape of this LogicalPart.
std::vector< double > rMinFront_
dd4hep::Solid solid() const
std::vector< int > iradMinBH_
std::vector< double > trformRotYY_
std::vector< double > cellFineX_
static constexpr double k_ScaleFromDDDToG4
const DDRotationMatrix & rotation() const
The absolute rotation of the current node.
std::vector< double > trformRotZZ_
int32_t waferPartial(const int32_t property)
const RotationMatrix rotation() const
static constexpr double k_ScaleFromDDD
std::vector< double > moduleAlphaS_
std::vector< int > layerGroupO_
void loadGeometryHexagon(const DDFilteredView &_fv, HGCalParameters &php, const std::string &sdTag1, const DDCompactView *cpv, const std::string &sdTag2, const std::string &sdTag3, HGCalGeometryMode::WaferMode mode)
void loadSpecParsTrapezoid(const DDFilteredView &fv, HGCalParameters &php)
std::vector< double > moduleBlS_
int32_t waferV(const int32_t index)
bool firstChild()
set the current node to the first child ...
std::vector< int > waferCopy_
std::vector< double > cassetteShift_
std::vector< int > depthIndex_
std::vector< int > layerFrontBH_
static int32_t packTypeUV(int type, int u, int v)
std::vector< double > rLimit_
std::vector< double > zFrontTop_
std::vector< double > radiusLayer_[2]
std::vector< int > waferTypeT_
std::vector< int > levelT_
std::vector< double > cellCoarseY_
Log< level::Warning, false > LogWarning
std::vector< std::pair< double, double > > tileRingR_
waferInfo_map waferInfoMap_
static int32_t layerType(int type)
std::vector< int > moduleLayS_
static constexpr int siliconFileHE
std::vector< double > trformTranZ_
const DDTranslation & translation() const
The absolute translation of the current node.
std::vector< double > waferPosX_
int scintCells(const int layer) const
void addTrForm(const CLHEP::Hep3Vector &h3v)
tileInfo_map tileInfoMap_
const std::vector< double > parameters() const
extract shape parameters
std::vector< double > moduleTlR_
std::vector< int > waferTypeL_
std::vector< double > xLayerHex_
int32_t tileProperty(const int32_t type, const int32_t sipm)
void loadWaferHexagon8(HGCalParameters &php)
std::vector< std::unordered_map< int32_t, int32_t > > layer_map
void loadGeometryHexagon8(const DDFilteredView &_fv, HGCalParameters &php, int firstLayer)
def cp(fromDir, toDir, listOfFiles, overwrite=False, smallList=False)
const Translation translation() const