107 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << matName <<
" initialized at " << &matter;
111 static const double sqrt3 =
std::sqrt(3.0);
113 double r = 0.5 *
R * sqrt3;
115 double dx2 = 0.5 * dx1;
116 double dx3 = 2.5 * dx2;
117 double dx4 = 0.5 * dx2;
119 double dy2 = 0.5 * dy1;
120 double dy3 = 1.5 * dy1;
121 std::vector<double>
xx = {dx1, dx2, -dx2, -dx1, -dx2, dx2, dx3, -dx4, -dx1, -dx1, -dx4, dx3};
122 std::vector<double>
yy = {0, dy1, dy1, 0, -dy1, -dy1, dy2, dy3, dy1, -dy1, -dy3, -dy2};
127 std::vector<double> xw = {
xx[0],
xx[1],
xx[2],
xx[3],
xx[4],
xx[5]};
128 std::vector<double> yw = {
yy[0],
yy[1],
yy[2],
yy[3],
yy[4],
yy[5]};
130 std::vector<double> zx(2, 0), zy(2, 0),
scale(2, 1.0);
134 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
135 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
136 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
137 <<
" and " << xw.size() <<
" edges";
138 for (
unsigned int k = 0;
k < xw.size(); ++
k)
145 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
146 <<
" z|x|y|s (0) " << zc[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
147 <<
" z|x|y|s (1) " << zc[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
148 <<
" and " << xw.size() <<
" edges";
149 for (
unsigned int k = 0;
k < xw.size(); ++
k)
155 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << glog2.
name() <<
" number 1 position in " << glog1.
name() <<
" at "
156 << tran <<
" with " <<
rot;
159 static const int ir0[] = {0, 1, 0};
160 static const int ir1[] = {1, 2, 1};
161 static const int ir2[] = {2, 3, 3};
162 static const int ir3[] = {3, 4, 4};
163 static const int ir4[] = {5, 5, 5};
164 for (
unsigned int i = 0;
i <
truncCN_.size(); ++
i) {
165 std::vector<double> xw = {
xx[ir0[
i]],
xx[ir1[
i]],
xx[ir2[
i]],
xx[ir3[
i]],
xx[ir4[
i]]};
166 std::vector<double> yw = {
yy[ir0[
i]],
yy[ir1[
i]],
yy[ir2[
i]],
yy[ir3[
i]],
yy[ir4[
i]]};
170 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
171 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
172 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
173 <<
" and " << xw.size() <<
" edges";
174 for (
unsigned int k = 0;
k < xw.size(); ++
k)
180 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
181 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
182 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
183 <<
" and " << xw.size() <<
" edges";
184 for (
unsigned int k = 0;
k < xw.size(); ++
k)
190 <<
" at " << tran <<
" with " <<
rot;
194 static const int ie0[] = {1, 5, 0};
195 static const int ie1[] = {2, 6, 1};
196 static const int ie2[] = {3, 7, 8};
197 static const int ie3[] = {10, 3, 9};
198 static const int ie4[] = {11, 4, 5};
199 for (
unsigned int i = 0;
i <
extenCN_.size(); ++
i) {
200 std::vector<double> xw = {
xx[ie0[
i]],
xx[ie1[
i]],
xx[ie2[
i]],
xx[ie3[
i]],
xx[ie4[
i]]};
201 std::vector<double> yw = {
yy[ie0[
i]],
yy[ie1[
i]],
yy[ie2[
i]],
yy[ie3[
i]],
yy[ie4[
i]]};
205 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
206 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
207 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
208 <<
" and " << xw.size() <<
" edges";
209 for (
unsigned int k = 0;
k < xw.size(); ++
k)
215 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
216 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
217 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
218 <<
" and " << xw.size() <<
" edges";
219 for (
unsigned int k = 0;
k < xw.size(); ++
k)
225 <<
" at " << tran <<
" with " <<
rot;
229 static const int ic0[] = {0, 1, 1, 1, 1, 0};
230 static const int ic1[] = {1, 2, 2, 7, 3, 1};
231 static const int ic2[] = {8, 3, 3, 3, 4, 3};
232 static const int ic3[] = {3, 5, 10, 4, 5, 9};
233 static const int ic4[] = {5, 11, 5, 5, 6, 5};
234 for (
unsigned int i = 0;
i <
cornrCN_.size(); ++
i) {
235 std::vector<double> xw = {
xx[ic0[
i]],
xx[ic1[
i]],
xx[ic2[
i]],
xx[ic3[
i]],
xx[ic4[
i]]};
236 std::vector<double> yw = {
yy[ic0[
i]],
yy[ic1[
i]],
yy[ic2[
i]],
yy[ic3[
i]],
yy[ic4[
i]]};
240 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
241 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
242 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
243 <<
" and " << xw.size() <<
" edges";
244 for (
unsigned int k = 0;
k < xw.size(); ++
k)
250 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalCell: " << solid.
name() <<
" extruded polygon made of " << matName
251 <<
" z|x|y|s (0) " <<
zw[0] <<
":" << zx[0] <<
":" << zy[0] <<
":" <<
scale[0]
252 <<
" z|x|y|s (1) " <<
zw[1] <<
":" << zx[1] <<
":" << zy[1] <<
":" <<
scale[1]
253 <<
" and " << xw.size() <<
" edges";
254 for (
unsigned int k = 0;
k < xw.size(); ++
k)
260 <<
" at " << tran <<
" with " <<
rot;