8 #include "DD4hep/DetFactoryHelper.h"
21 const auto& waferSize =
args.value<
double>(
"WaferSize");
22 const auto& waferT =
args.value<
double>(
"WaferThick");
23 const auto& waferSepar =
args.value<
double>(
"SensorSeparation");
24 const auto&
nCells =
args.value<
int>(
"NCells");
25 const auto& cellType =
args.value<
int>(
"CellType");
27 const auto& cellNames =
args.value<std::vector<std::string>>(
"CellNames");
32 <<
" Cells/Wafer " <<
nCells <<
" Cell Type " << cellType <<
" Material " << material
33 <<
" Names " << motherName <<
" NameSpace " << ns.
name() <<
" # of cells "
35 for (
unsigned int k = 0;
k < cellNames.size(); ++
k)
39 static const double sqrt3 =
std::sqrt(3.0);
40 double rM = 0.5 * (waferSize + waferSepar);
41 double RM2 = rM / sqrt3;
42 double R = waferSize / (3.0 *
nCells);
43 double r = 0.5 *
R * sqrt3;
46 std::vector<double> xM = {rM, 0, -rM, -rM, 0, rM};
47 std::vector<double> yM = {RM2, 2 * RM2, RM2, -RM2, -2 * RM2, -RM2};
48 std::vector<double>
zw = {-0.5 * waferT, 0.5 * waferT};
49 std::vector<double> zx(2, 0), zy(2, 0),
scale(2, 1.0);
58 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalWafer8: " << solid.name() <<
" extruded polygon made of " << material
63 for (
unsigned int k = 0;
k < xM.size(); ++
k)
71 for (
int u = 0; u < 2 *
nCells; ++u) {
78 double yp = (u - 0.5 *
v - n2) * 2 *
r;
79 double xp = (1.5 * (
v -
nCells) + 1.0) *
R;
81 if ((u == 0) && (
v == 0))
83 else if ((u == 0) && (
v ==
nCells - 1))
91 else if ((u ==
nCells) && (
v == 0))
95 else if ((
v - u) == (
nCells - 1))
99 else if (u == (2 *
nCells - 1))
107 glog.placeVolume(ns.
volume(cellNames[cell]),
copy, dd4hep::Transform3D(
rotation, tran));
109 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalWafer8: " << cellNames[cell] <<
" number " <<
copy <<
" position in "
111 <<
",0) with no rotation";
117 edm::LogVerbatim(
"HGCalGeom") <<
"\nDDHGCalWafer8::Counter : " <<
counter <<
"\n===============================\n";