81 static const double sqrt3 =
std::sqrt(3.0);
83 double RM2 = rM / sqrt3;
85 double r = 0.5 * R * sqrt3;
88 std::vector<double> xM = {rM, 0, -rM, -rM, 0, rM};
89 std::vector<double> yM = {RM2, 2 * RM2, RM2, -RM2, -2 * RM2, -RM2};
91 std::vector<double> zx(2, 0), zy(2, 0),
scale(2, 1.0);
98 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalWafer8: " << solid.
name() <<
" extruded polygon made of " << matName
99 <<
" z|x|y|s (0) " << zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
100 <<
" z|x|y|s (1) " << zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
101 <<
" and " << xM.size() <<
" edges";
102 for (
unsigned int k = 0;
k < xM.size(); ++
k)
107 for (
int u = 0; u < 2 *
nCells_; ++u) {
109 if (((
v - u) < nCells_) && (u -
v) <= nCells_) {
110 int n2 = nCells_ / 2;
111 double yp = (u - 0.5 *
v - n2) * 2 * r;
112 double xp = (1.5 * (
v -
nCells_) + 1.0) *
R;
114 if ((u == 0) && (
v == 0))
116 else if ((u == 0) && (
v == nCells_ - 1))
118 else if ((u == nCells_) && (
v == 2 * nCells_ - 1))
120 else if ((u == 2 * nCells_ - 1) && (
v == 2 * nCells_ - 1))
122 else if ((u == 2 * nCells_ - 1) && (
v == nCells_ - 1))
124 else if ((u == nCells_) && (
v == 0))
128 else if ((
v - u) == (nCells_ - 1))
130 else if (
v == (2 * nCells_ - 1))
132 else if (u == (2 * nCells_ - 1))
134 else if ((u -
v) == nCells_)
143 << glog.
name() <<
" at " << tran <<
" with " <<
rot;
DDMaterial is used to define and access material information.
DDName is used to identify DDD entities uniquely.
A DDSolid represents the shape of a part.
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DDTranslation
Represents a uniquely identifyable rotation matrix.
U second(std::pair< T, U > const &p)
std::vector< std::string > cellNames_
A DDLogicalPart aggregates information concerning material, solid and sensitveness ...
static DDSolid extrudedpolygon(const DDName &name, const std::vector< double > &x, const std::vector< double > &y, const std::vector< double > &z, const std::vector< double > &zx, const std::vector< double > &zy, const std::vector< double > &zscale)
void position(const DDLogicalPart &self, const DDLogicalPart &parent, const std::string ©no, const DDTranslation &trans, const DDRotation &rot, const DDDivision *div=0)
auto zw(V v) -> Vec2< typename std::remove_reference< decltype(v[0])>::type >
std::pair< std::string, std::string > DDSplit(const std::string &n)
split into (name,namespace), separator = ':'