21 edm::LogVerbatim(
"HGCalGeom") <<
"HGCalParametersFromDD (DDD)::build called with " 22 <<
"names " <<
name <<
":" << namew <<
":" << namec <<
":" << namet;
55 std::unique_ptr<HGCalGeomParameters>
geom = std::make_unique<HGCalGeomParameters>();
57 attribute =
"OnlyForHGCalNumbering";
79 php.
levelZSide_ =
static_cast<int>(getDDDValue(
"LevelZSide",
sv));
86 php.
waferZSide_ =
static_cast<int>(getDDDValue(
"WaferZside",
sv));
88 php.
useSimWt_ =
static_cast<int>(getDDDValue(
"UseSimWt",
sv));
102 attribute =
"OnlyForHGCalNumbering";
111 php.
nCellsFine_ =
static_cast<int>(getDDDValue(
"NumberOfCellsFine", sv2));
112 php.
nCellsCoarse_ =
static_cast<int>(getDDDValue(
"NumberOfCellsCoarse", sv2));
119 php.
useOffset_ =
static_cast<int>(getDDDValue(
"UseOffset", sv2));
132 for (
int k = 0;
k < 2; ++
k)
133 getCellPosition(php,
k);
139 geom->loadSpecParsHexagon(fv, php, cpv, namew, namec);
141 geom->loadGeometryHexagon(fv, php,
name, cpv, namew, namec,
mode);
143 geom->loadCellParsHexagon(cpv, php);
149 geom->loadSpecParsHexagon(fv, php, cpv, namew, namec);
151 geom->loadGeometryHexagon(fv, php,
name, cpv, namew, namec,
mode);
153 geom->loadWaferHexagon(php);
155 geom->loadCellParsHexagon(cpv, php);
160 geom->loadSpecParsHexagon8(fv, php);
162 geom->loadGeometryHexagon8(fv, php, 1);
166 geom->loadWaferHexagon8(php);
169 geom->loadSpecParsHexagon8(fv, php);
171 geom->loadGeometryHexagon8(fv, php, 1);
175 geom->loadWaferHexagon8(php);
178 geom->loadSpecParsHexagon8(fv, php);
180 geom->loadGeometryHexagonModule(cpv, php,
name, namec, 1);
184 geom->loadWaferHexagon8(php);
199 php.
waferZSide_ =
static_cast<int>(getDDDValue(
"WaferZside",
sv));
201 php.
useSimWt_ =
static_cast<int>(getDDDValue(
"UseSimWt",
sv));
212 geom->loadSpecParsTrapezoid(fv, php);
216 geom->loadCellTrapezoid(php);
218 edm::LogError(
"HGCalGeom") <<
"Unknown Geometry type " << php.
mode_ <<
" for HGCal " <<
name <<
":" << namew
221 <<
"Unknown Geometry type " << php.
mode_ <<
" for HGCal " <<
name <<
":" << namew <<
":" << namec;
224 edm::LogError(
"HGCalGeom") <<
" Attribute " <<
val <<
" not found but needed.";
225 throw cms::Exception(
"DDException") <<
"Attribute " <<
val <<
" not found but needed.";
229 <<
" with flag " <<
ok;
242 edm::LogVerbatim(
"HGCalGeom") <<
"HGCalParametersFromDD (DD4hep)::build called with " 243 <<
"names " <<
name <<
":" << namew <<
":" << namec <<
":" << namet <<
":" <<
name2;
248 std::vector<std::string> tempS;
249 std::vector<double> tempD;
251 tempS = fv.
get<std::vector<std::string> >(
name2,
"GeometryMode");
253 tempS = fv.
get<std::vector<std::string> >(
name,
"GeometryMode");
255 std::string sv = (!tempS.empty()) ? tempS[0] :
"HGCalGeometryMode::Hexagon8Full";
278 std::unique_ptr<HGCalGeomParameters>
geom = std::make_unique<HGCalGeomParameters>();
280 tempS = fv.
get<std::vector<std::string> >(namet,
"WaferMode");
281 std::string sv2 = (!tempS.empty()) ? tempS[0] :
"HGCalGeometryMode::Polyhedra";
295 tempD = fv.
get<std::vector<double> >(
name,
"LevelZSide");
299 tempD = fv.
get<std::vector<double> >(
name,
"FirstMixedLayer");
301 tempD = fv.
get<std::vector<double> >(
name,
"DetectorType");
304 tempD = fv.
get<std::vector<double> >(
name,
"WaferMaskMode");
306 tempD = fv.
get<std::vector<double> >(
name,
"WaferZside");
309 tempD = fv.
get<std::vector<double> >(
name,
"LayerRotation");
311 tempD = fv.
get<std::vector<double> >(
name,
"UseSimWt");
316 tempD = fv.
get<std::vector<double> >(
name,
"Cassettes");
328 tempS = fv.
get<std::vector<std::string> >(namet,
"WaferMode");
329 std::string sv2 = (!tempS.empty()) ? tempS[0] :
"HGCalGeometryMode::ExtrudedPolygon";
331 tempD = fv.
get<std::vector<double> >(namet,
"NumberOfCellsFine");
333 tempD = fv.
get<std::vector<double> >(namet,
"NumberOfCellsCoarse");
335 tempD = fv.
get<std::vector<double> >(namet,
"WaferSize");
337 tempD = fv.
get<std::vector<double> >(namet,
"WaferThickness");
339 tempD = fv.
get<std::vector<double> >(namet,
"SensorSeparation");
341 tempD = fv.
get<std::vector<double> >(namet,
"SensorSizeOffset");
343 tempD = fv.
get<std::vector<double> >(namet,
"GuardRingOffset");
345 tempD = fv.
get<std::vector<double> >(namet,
"MouseBite");
347 tempD = fv.
get<std::vector<double> >(namet,
"UseOffset");
361 for (
int k = 0;
k < 2; ++
k)
362 getCellPosition(php,
k);
367 geom->loadSpecParsHexagon(fv, php,
name, namew, namec,
name2);
369 geom->loadGeometryHexagon(cpv, php,
name, namew, namec,
mode);
371 geom->loadCellParsHexagon(vmap, php);
377 geom->loadSpecParsHexagon(fv, php,
name, namew, namec,
name2);
379 geom->loadGeometryHexagon(cpv, php,
name, namew, namec,
mode);
381 geom->loadWaferHexagon(php);
383 geom->loadCellParsHexagon(vmap, php);
388 geom->loadSpecParsHexagon8(fv, vmap, php,
name);
390 geom->loadGeometryHexagon8(cpv, php,
name, 1);
394 geom->loadWaferHexagon8(php);
397 geom->loadSpecParsHexagon8(fv, vmap, php,
name);
399 geom->loadGeometryHexagon8(cpv, php,
name, 1);
403 geom->loadWaferHexagon8(php);
406 geom->loadSpecParsHexagon8(fv, vmap, php,
name);
408 geom->loadGeometryHexagonModule(cpv, php,
name, namec, 1);
412 geom->loadWaferHexagon8(php);
418 tempD = fv.
get<std::vector<double> >(
name,
"LevelZSide");
421 tempD = fv.
get<std::vector<double> >(
name,
"FirstLayer");
423 tempD = fv.
get<std::vector<double> >(
name,
"FirstMixedLayer");
425 tempD = fv.
get<std::vector<double> >(
name,
"DetectorType");
427 tempD = fv.
get<std::vector<double> >(
name,
"WaferThickness");
429 tempD = fv.
get<std::vector<double> >(
name,
"MinimumTileSize");
434 tempD = fv.
get<std::vector<double> >(
name,
"WaferMaskMode");
436 tempD = fv.
get<std::vector<double> >(
name,
"WaferZside");
439 tempD = fv.
get<std::vector<double> >(
name,
"UseSimWt");
443 tempD = fv.
get<std::vector<double> >(
name,
"Cassettes");
455 geom->loadSpecParsTrapezoid(fv, vmap, php,
name);
459 geom->loadCellTrapezoid(php);
461 edm::LogError(
"HGCalGeom") <<
"Unknown Geometry type " << php.
mode_ <<
" for HGCal " <<
name <<
":" << namew
464 <<
"Unknown Geometry type " << php.
mode_ <<
" for HGCal " <<
name <<
":" << namew <<
":" << namec;
467 edm::LogError(
"HGCalGeom") <<
" Attribute Volume:" <<
name <<
" not found but needed.";
468 throw cms::Exception(
"DDException") <<
"Attribute Volume:" <<
name <<
" not found but needed.";
472 <<
" with flag " <<
ok;
487 std::vector<int> indtypes;
491 double r = 0.5 *
R *
sqrt(3.0);
494 for (
int u = 0; u < 2 *
N; ++u) {
495 for (
int v = 0;
v < 2 *
N; ++
v) {
496 if (((
v - u) <
N) && (u -
v) <=
N) {
497 double yp = (u - 0.5 *
v - n2) * 2 * r;
498 double xp = (1.5 * (
v -
N) + 1.0) *
R;
499 int id =
v * 100 + u;
501 indtypes.emplace_back(
id);
510 cellIndex[
id] = ipos;
525 int id = indtypes[
k];
532 int id = indtypes[
k];
543 const std::vector<double>& fvec =
val.doubles();
545 throw cms::Exception(
"HGCalGeom") <<
"getDDDValue::Failed to get " <<
s <<
" tag.";
549 throw cms::Exception(
"HGCalGeom") <<
"getDDDValue::Failed to fetch " <<
s <<
" tag";
556 const std::vector<double>& fvec =
val.doubles();
558 throw cms::Exception(
"HGCalGeom") <<
"getDDDArray::Failed to get " <<
s <<
" tag.";
562 throw cms::Exception(
"HGCalGeom") <<
"getDDDArray:Failed to fetch " <<
s <<
" tag";
Log< level::Info, true > LogVerbatim
static constexpr int scintillatorCassette
constexpr NumType convertRadToDeg(NumType radians)
std::vector< double > cellFineY_
const cms::DDDetector * detector() const
Log< level::Error, false > LogError
bool build(const DDCompactView *cpv, HGCalParameters &php, const std::string &name, const std::string &namew, const std::string &namec, const std::string &namet)
Compact representation of the geometrical detector hierarchy.
bool DDfetch(const DDsvalues_type *, DDValue &)
helper for retrieving DDValues from DDsvalues_type *.
HGCalGeometryMode::GeometryMode mode_
std::unordered_map< std::string, std::vector< double > > DDVectorsMap
WaferMode getGeometryWaferMode(const char *s, const DDsvalues_type &sv)
std::vector< double > cellCoarseX_
static constexpr int siliconCassetteHE
std::vector< std::pair< unsigned int, DDValue > > DDsvalues_type
std::vector< double > cellSize_
static constexpr int siliconCassetteEE
GeometryMode getGeometryMode(const char *s, const DDsvalues_type &sv)
static constexpr double k_ScaleFromDD4hep
bool firstChild()
set the current node to the first child
std::vector< int > dbl_to_int(const std::vector< double > &vecdbl)
Converts a std::vector of doubles to a std::vector of int.
DDsvalues_type mergedSpecifics() const
T get(const std::string &)
extract attribute value
static constexpr double k_ScaleToDDD
std::unordered_map< int32_t, int32_t > wafer_map
double getDDDValue(const char *s, const DDsvalues_type &sv)
wafer_map cellCoarseIndex_
std::vector< double > cellFineX_
static constexpr double k_ScaleFromDDD
bool firstChild()
set the current node to the first child ...
std::vector< int > levelT_
std::vector< double > cellCoarseY_
cms::DDVectorsMap const & vectors() const
std::vector< double > getDDDArray(const char *s, const DDsvalues_type &sv)
void getCellPosition(HGCalParameters &php, int type)