94 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << matName <<
" initialized";
98 static const double sqrt3 =
std::sqrt(3.0);
100 double r = 0.5 *
R * sqrt3;
102 double dx2 = 0.5 * dx1;
103 double dx3 = 2.5 * dx2;
104 double dx4 = 0.5 * dx2;
106 double dy2 = 0.5 * dy1;
107 double dy3 = 1.5 * dy1;
108 std::vector<double>
xx = {dx1, dx2, -dx2, -dx1, -dx2, dx2, dx3, -dx4, -dx1, -dx1, -dx4, dx3};
109 std::vector<double>
yy = {0, dy1, dy1, 0, -dy1, -dy1, dy2, dy3, dy1, -dy1, -dy3, -dy2};
114 std::vector<double> xw = {
xx[0],
xx[1],
xx[2],
xx[3],
xx[4],
xx[5]};
115 std::vector<double> yw = {
yy[0],
yy[1],
yy[2],
yy[3],
yy[4],
yy[5]};
117 std::vector<double> zx(2, 0), zy(2, 0),
scale(2, 1.0);
121 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
122 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
123 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
124 <<
" and " << xw.size() <<
" edges";
125 for (
unsigned int k = 0;
k < xw.size(); ++
k)
132 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
133 <<
" z|x|y|s (0) " << zc[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
134 <<
" z|x|y|s (1) " << zc[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
135 <<
" and " << xw.size() <<
" edges";
136 for (
unsigned int k = 0;
k < xw.size(); ++
k)
142 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << glog2.
name() <<
" number 1 position in " << glog1.
name() <<
" at "
143 << tran <<
" with no rotation";
146 static constexpr
int ir0[] = {0, 1, 0};
147 static constexpr
int ir1[] = {1, 2, 1};
148 static constexpr
int ir2[] = {2, 3, 3};
149 static constexpr
int ir3[] = {3, 4, 4};
150 static constexpr
int ir4[] = {5, 5, 5};
151 for (
unsigned int i = 0;
i <
truncCN_.size(); ++
i) {
152 std::vector<double> xw = {
xx[ir0[
i]],
xx[ir1[
i]],
xx[ir2[
i]],
xx[ir3[
i]],
xx[ir4[
i]]};
153 std::vector<double> yw = {
yy[ir0[
i]],
yy[ir1[
i]],
yy[ir2[
i]],
yy[ir3[
i]],
yy[ir4[
i]]};
157 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
158 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
159 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
160 <<
" and " << xw.size() <<
" edges";
161 for (
unsigned int k = 0;
k < xw.size(); ++
k)
167 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
168 <<
" z|x|y|s (0) " << zc[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
169 <<
" z|x|y|s (1) " << zc[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
170 <<
" and " << xw.size() <<
" edges";
171 for (
unsigned int k = 0;
k < xw.size(); ++
k)
177 <<
" at " << tran <<
" with no rotation";
181 static constexpr
int ie0[] = {1, 5, 0};
182 static constexpr
int ie1[] = {2, 6, 1};
183 static constexpr
int ie2[] = {3, 7, 8};
184 static constexpr
int ie3[] = {10, 3, 9};
185 static constexpr
int ie4[] = {11, 4, 5};
186 for (
unsigned int i = 0;
i <
extenCN_.size(); ++
i) {
187 std::vector<double> xw = {
xx[ie0[
i]],
xx[ie1[
i]],
xx[ie2[
i]],
xx[ie3[
i]],
xx[ie4[
i]]};
188 std::vector<double> yw = {
yy[ie0[
i]],
yy[ie1[
i]],
yy[ie2[
i]],
yy[ie3[
i]],
yy[ie4[
i]]};
192 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
193 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
194 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
195 <<
" and " << xw.size() <<
" edges";
196 for (
unsigned int k = 0;
k < xw.size(); ++
k)
202 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
203 <<
" z|x|y|s (0) " << zc[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
204 <<
" z|x|y|s (1) " << zc[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
205 <<
" and " << xw.size() <<
" edges";
206 for (
unsigned int k = 0;
k < xw.size(); ++
k)
212 <<
" at " << tran <<
" with no rotation";
216 static constexpr
int ic0[] = {0, 1, 1, 1, 1, 0};
217 static constexpr
int ic1[] = {1, 2, 2, 7, 3, 1};
218 static constexpr
int ic2[] = {8, 3, 3, 3, 4, 3};
219 static constexpr
int ic3[] = {3, 5, 10, 4, 5, 9};
220 static constexpr
int ic4[] = {5, 11, 5, 5, 6, 5};
221 for (
unsigned int i = 0;
i <
cornrCN_.size(); ++
i) {
222 std::vector<double> xw = {
xx[ic0[
i]],
xx[ic1[
i]],
xx[ic2[
i]],
xx[ic3[
i]],
xx[ic4[
i]]};
223 std::vector<double> yw = {
yy[ic0[
i]],
yy[ic1[
i]],
yy[ic2[
i]],
yy[ic3[
i]],
yy[ic4[
i]]};
227 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
228 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
229 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
230 <<
" and " << xw.size() <<
" edges";
231 for (
unsigned int k = 0;
k < xw.size(); ++
k)
237 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
238 <<
" z|x|y|s (0) " << zc[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
239 <<
" z|x|y|s (1) " << zc[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
240 <<
" and " << xw.size() <<
" edges";
241 for (
unsigned int k = 0;
k < xw.size(); ++
k)
247 <<
" at " << tran <<
" with no rotation";