CMS 3D CMS Logo

List of all members | Public Member Functions | Public Attributes | Static Public Attributes
HGCalMixRotatedLayer Struct Reference

Public Member Functions

 HGCalMixRotatedLayer ()
 
 HGCalMixRotatedLayer (cms::DDParsingContext &ctxt, xml_h e)
 
void positionMix (cms::DDParsingContext &ctxt, xml_h e, const dd4hep::Volume &glog, const std::string &nameM, int copyM, double thick, const dd4hep::Material &matter)
 

Public Attributes

int absorbMode_
 
double alpha_
 
HGCalCassette cassette_
 
int cassettes_
 
std::vector< double > cassetteShift_
 
std::unordered_set< int > copies_
 
std::vector< int > copyNumber_
 
std::vector< int > copyNumberTop_
 
double cosAlpha_
 
int facingTypes_
 
int firstLayer_
 
int forFireworks_
 
HGCalGeomTools geomTools_
 
std::vector< int > layerOrient_
 
std::vector< int > layers_
 
std::vector< int > layerSense_
 
std::vector< double > layerThick_
 
std::vector< double > layerThickTop_
 
std::vector< int > layerType_
 
std::vector< int > layerTypeTop_
 
std::vector< std::string > materials_
 
std::vector< std::string > materialTop_
 
std::vector< std::string > names_
 
std::vector< std::string > namesTop_
 
int orientationTypes_
 
int phiBinsScint_
 
int placeOffset_
 
std::vector< double > rMaxFront_
 
std::vector< double > rMinFront_
 
int sectors_
 
int sensitiveMode_
 
std::vector< double > slopeB_
 
std::vector< double > slopeT_
 
std::vector< double > thick_
 
std::vector< int > tileIndex_
 
std::vector< int > tileLayerStart_
 
std::vector< int > tilePhis_
 
std::vector< double > tileRMax_
 
std::vector< double > tileRMin_
 
std::vector< std::string > waferFull_
 
std::vector< int > waferIndex_
 
std::vector< int > waferLayerStart_
 
std::vector< std::string > waferPart_
 
std::vector< int > waferProperty_
 
double waferSepar_
 
double waferSize_
 
int waferTypes_
 
std::vector< double > zFrontB_
 
std::vector< double > zFrontT_
 
double zMinBlock_
 

Static Public Attributes

static constexpr double tol2_ = 0.00001 * dd4hep::mm
 

Detailed Description

Definition at line 30 of file DDHGCalMixRotatedLayer.cc.

Constructor & Destructor Documentation

◆ HGCalMixRotatedLayer() [1/2]

HGCalMixRotatedLayer::HGCalMixRotatedLayer ( )
inline

Definition at line 31 of file DDHGCalMixRotatedLayer.cc.

References Exception.

31 { throw cms::Exception("HGCalGeom") << "Wrong initialization to HGCalMixRotatedLayer"; }

◆ HGCalMixRotatedLayer() [2/2]

HGCalMixRotatedLayer::HGCalMixRotatedLayer ( cms::DDParsingContext ctxt,
xml_h  e 
)
inline

Definition at line 32 of file DDHGCalMixRotatedLayer.cc.

References funct::abs(), cms::DDNamespace::addSolidNS(), cms::DDNamespace::addVolumeNS(), writedatasetfile::args, cms::convert2mm(), angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, funct::cos(), MillePedeFileConverter_cfg::e, mps_fire::i, testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, cuy::ii, dqmiolumiharvest::j, dqmdumpme::k, HGCalTypes::layerType(), cms::DDNamespace::material(), g4SimHits_cfi::Material, SiStripPI::min, Skims_PA_cff::name, cms::DDNamespace::name(), HGCalGeometryMode::Polyhedra, PixelTestBeamValidation_cfi::Position, cms::DDNamespace::prepend(), HGCalGeomTools::radius(), AlCaHLTBitMon_QueryRunRegistry::string, HGCalTileIndex::tileUnpack(), cond::impl::to_string(), cms::DDNamespace::volume(), HGCalWaferIndex::waferLayer(), HGCalProperty::waferOrient(), HGCalProperty::waferPartial(), HGCalProperty::waferThick(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), and geometryCSVtoXML::zz.

32  {
33  cms::DDNamespace ns(ctxt, e, true);
35 
36 #ifdef EDM_ML_DEBUG
37  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: Creating an instance";
38 #endif
39 
40  static constexpr double tol1 = 0.01 * dd4hep::mm;
41  dd4hep::Volume mother = ns.volume(args.parentName());
42 
43  waferTypes_ = args.value<int>("WaferTypes");
44  facingTypes_ = args.value<int>("FacingTypes");
45  orientationTypes_ = args.value<int>("OrientationTypes");
46  placeOffset_ = args.value<int>("PlaceOffset");
47  phiBinsScint_ = args.value<int>("NPhiBinScint");
48  forFireworks_ = args.value<int>("ForFireWorks");
49 #ifdef EDM_ML_DEBUG
50  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer::Number of types of wafers: " << waferTypes_
51  << " facings: " << facingTypes_ << " Orientations: " << orientationTypes_
52  << " PlaceOffset: " << placeOffset_ << "; number of cells along phi " << phiBinsScint_
53  << " forFireworks_: " << forFireworks_;
54 #endif
55  firstLayer_ = args.value<int>("FirstLayer");
56  absorbMode_ = args.value<int>("AbsorberMode");
57  sensitiveMode_ = args.value<int>("SensitiveMode");
58 #ifdef EDM_ML_DEBUG
59  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer::First Layer " << firstLayer_ << " and "
60  << "Absober:Sensitive mode " << absorbMode_ << ":" << sensitiveMode_;
61 #endif
62  zMinBlock_ = args.value<double>("zMinBlock");
63  waferSize_ = args.value<double>("waferSize");
64  waferSepar_ = args.value<double>("SensorSeparation");
65  sectors_ = args.value<int>("Sectors");
66  cassettes_ = args.value<int>("Cassettes");
67  alpha_ = (1._pi) / sectors_;
68  cosAlpha_ = cos(alpha_);
69 #ifdef EDM_ML_DEBUG
70  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: zStart " << cms::convert2mm(zMinBlock_) << " wafer width "
71  << cms::convert2mm(waferSize_) << " separations " << cms::convert2mm(waferSepar_)
72  << " sectors " << sectors_ << ":" << convertRadToDeg(alpha_) << ":" << cosAlpha_
73  << " with " << cassettes_ << " cassettes";
74 #endif
75  slopeB_ = args.value<std::vector<double>>("SlopeBottom");
76  zFrontB_ = args.value<std::vector<double>>("ZFrontBottom");
77  rMinFront_ = args.value<std::vector<double>>("RMinFront");
78  slopeT_ = args.value<std::vector<double>>("SlopeTop");
79  zFrontT_ = args.value<std::vector<double>>("ZFrontTop");
80  rMaxFront_ = args.value<std::vector<double>>("RMaxFront");
81 #ifdef EDM_ML_DEBUG
82  for (unsigned int i = 0; i < slopeB_.size(); ++i)
83  edm::LogVerbatim("HGCalGeom") << "Bottom Block [" << i << "] Zmin " << cms::convert2mm(zFrontB_[i]) << " Rmin "
84  << cms::convert2mm(rMinFront_[i]) << " Slope " << slopeB_[i];
85  for (unsigned int i = 0; i < slopeT_.size(); ++i)
86  edm::LogVerbatim("HGCalGeom") << "Top Block [" << i << "] Zmin " << cms::convert2mm(zFrontT_[i]) << " Rmax "
87  << cms::convert2mm(rMaxFront_[i]) << " Slope " << slopeT_[i];
88 #endif
89 
90  waferFull_ = args.value<std::vector<std::string>>("WaferNamesFull");
91  waferPart_ = args.value<std::vector<std::string>>("WaferNamesPartial");
92 #ifdef EDM_ML_DEBUG
93  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: " << waferFull_.size() << " full and "
94  << waferPart_.size() << " partial modules\nDDHGCalMixRotatedLayer:Full Modules:";
95  unsigned int i1max = static_cast<unsigned int>(waferFull_.size());
96  for (unsigned int i1 = 0; i1 < i1max; i1 += 2) {
97  std::ostringstream st1;
98  unsigned int i2 = std::min((i1 + 2), i1max);
99  for (unsigned int i = i1; i < i2; ++i)
100  st1 << " [" << i << "] " << waferFull_[i];
101  edm::LogVerbatim("HGCalGeom") << st1.str();
102  }
103  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: Partial Modules:";
104  i1max = static_cast<unsigned int>(waferPart_.size());
105  for (unsigned int i1 = 0; i1 < i1max; i1 += 2) {
106  std::ostringstream st1;
107  unsigned int i2 = std::min((i1 + 2), i1max);
108  for (unsigned int i = i1; i < i2; ++i)
109  st1 << " [" << i << "] " << waferPart_[i];
110  edm::LogVerbatim("HGCalGeom") << st1.str();
111  }
112 #endif
113 
114  materials_ = args.value<std::vector<std::string>>("MaterialNames");
115  names_ = args.value<std::vector<std::string>>("VolumeNames");
116  thick_ = args.value<std::vector<double>>("Thickness");
117  copyNumber_.resize(materials_.size(), 1);
118 #ifdef EDM_ML_DEBUG
119  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: " << materials_.size() << " types of volumes";
120  for (unsigned int i = 0; i < names_.size(); ++i)
121  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << names_[i] << " of thickness "
122  << cms::convert2mm(thick_[i]) << " filled with " << materials_[i]
123  << " first copy number " << copyNumber_[i];
124 #endif
125  layers_ = args.value<std::vector<int>>("Layers");
126  layerThick_ = args.value<std::vector<double>>("LayerThick");
127 #ifdef EDM_ML_DEBUG
128  edm::LogVerbatim("HGCalGeom") << "There are " << layers_.size() << " blocks";
129  for (unsigned int i = 0; i < layers_.size(); ++i)
130  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] of thickness " << cms::convert2mm(layerThick_[i])
131  << " with " << layers_[i] << " layers";
132 #endif
133  layerType_ = args.value<std::vector<int>>("LayerType");
134  layerSense_ = args.value<std::vector<int>>("LayerSense");
135  layerOrient_ = args.value<std::vector<int>>("LayerTypes");
136  for (unsigned int k = 0; k < layerOrient_.size(); ++k)
138 #ifdef EDM_ML_DEBUG
139  for (unsigned int i = 0; i < layerOrient_.size(); ++i)
140  edm::LogVerbatim("HGCalGeom") << "LayerTypes [" << i << "] " << layerOrient_[i];
141 #endif
142  if (firstLayer_ > 0) {
143  for (unsigned int i = 0; i < layerType_.size(); ++i) {
144  if (layerSense_[i] > 0) {
145  int ii = layerType_[i];
147 #ifdef EDM_ML_DEBUG
148  edm::LogVerbatim("HGCalGeom") << "First copy number for layer type " << i << ":" << ii << " with "
149  << materials_[ii] << " changed to " << copyNumber_[ii];
150 #endif
151  break;
152  }
153  }
154  } else {
155  firstLayer_ = 1;
156  }
157 #ifdef EDM_ML_DEBUG
158  edm::LogVerbatim("HGCalGeom") << "There are " << layerType_.size() << " layers";
159  for (unsigned int i = 0; i < layerType_.size(); ++i)
160  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerType_[i] << " sensitive class "
161  << layerSense_[i];
162 #endif
163  materialTop_ = args.value<std::vector<std::string>>("TopMaterialNames");
164  namesTop_ = args.value<std::vector<std::string>>("TopVolumeNames");
165  layerThickTop_ = args.value<std::vector<double>>("TopLayerThickness");
166  layerTypeTop_ = args.value<std::vector<int>>("TopLayerType");
167  copyNumberTop_.resize(materialTop_.size(), firstLayer_);
168 #ifdef EDM_ML_DEBUG
169  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: " << materialTop_.size()
170  << " types of volumes in the top part";
171  for (unsigned int i = 0; i < materialTop_.size(); ++i)
172  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << namesTop_[i] << " of thickness "
173  << cms::convert2mm(layerThickTop_[i]) << " filled with " << materialTop_[i]
174  << " first copy number " << copyNumberTop_[i];
175  edm::LogVerbatim("HGCalGeom") << "There are " << layerTypeTop_.size() << " layers in the top part";
176  for (unsigned int i = 0; i < layerTypeTop_.size(); ++i)
177  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerTypeTop_[i];
178 #endif
179  waferIndex_ = args.value<std::vector<int>>("WaferIndex");
180  waferProperty_ = args.value<std::vector<int>>("WaferProperties");
181  waferLayerStart_ = args.value<std::vector<int>>("WaferLayerStart");
182  cassetteShift_ = args.value<std::vector<double>>("CassetteShift");
183 #ifdef EDM_ML_DEBUG
184  edm::LogVerbatim("HGCalGeom") << "waferProperties with " << waferIndex_.size() << " entries in "
185  << waferLayerStart_.size() << " layers";
186  for (unsigned int k = 0; k < waferLayerStart_.size(); ++k)
187  edm::LogVerbatim("HGCalGeom") << "LayerStart[" << k << "] " << waferLayerStart_[k];
188  for (unsigned int k = 0; k < waferIndex_.size(); ++k)
189  edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << waferIndex_[k] << " ("
192  << HGCalWaferIndex::waferV(waferIndex_[k]) << ") : ("
196  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: " << cassetteShift_.size()
197  << " elements for cassette shifts";
198  unsigned int j1max = cassetteShift_.size();
199  for (unsigned int j1 = 0; j1 < j1max; j1 += 6) {
200  std::ostringstream st1;
201  unsigned int j2 = std::min((j1 + 6), j1max);
202  for (unsigned int j = j1; j < j2; ++j)
203  st1 << " [" << j << "] " << std::setw(9) << cassetteShift_[j];
204  edm::LogVerbatim("HGCalGeom") << st1.str();
205  }
206 #endif
207  tileRMin_ = args.value<std::vector<double>>("TileRMin");
208  tileRMax_ = args.value<std::vector<double>>("TileRMax");
209  tileIndex_ = args.value<std::vector<int>>("TileLayerRings");
210  tilePhis_ = args.value<std::vector<int>>("TilePhiRange");
211  tileLayerStart_ = args.value<std::vector<int>>("TileLayerStart");
212 #ifdef EDM_ML_DEBUG
213  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer:: with " << tileRMin_.size() << " rings";
214  for (unsigned int k = 0; k < tileRMin_.size(); ++k)
215  edm::LogVerbatim("HGCalGeom") << "Ring[" << k << "] " << cms::convert2mm(tileRMin_[k]) << " : "
217  edm::LogVerbatim("HGCalGeom") << "TileProperties with " << tileIndex_.size() << " entries in "
218  << tileLayerStart_.size() << " layers";
219  for (unsigned int k = 0; k < tileLayerStart_.size(); ++k)
220  edm::LogVerbatim("HGCalGeom") << "LayerStart[" << k << "] " << tileLayerStart_[k];
221  for (unsigned int k = 0; k < tileIndex_.size(); ++k)
222  edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << tileIndex_[k] << " ("
223  << "Layer " << std::get<0>(HGCalTileIndex::tileUnpack(tileIndex_[k])) << " Ring "
224  << std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[k])) << ":"
225  << std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[k])) << ") Phi "
226  << std::get<1>(HGCalTileIndex::tileUnpack(tilePhis_[k])) << ":"
227  << std::get<2>(HGCalTileIndex::tileUnpack(tilePhis_[k]));
228 #endif
230 
231 #ifdef EDM_ML_DEBUG
232  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: NameSpace " << ns.name();
233 
234  edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalMixRotatedLayer...";
235  copies_.clear();
236 #endif
237 
238  double zi(zMinBlock_);
239  int laymin(0);
240  for (unsigned int i = 0; i < layers_.size(); i++) {
241  double zo = zi + layerThick_[i];
242  double routF = HGCalGeomTools::radius(zi, zFrontT_, rMaxFront_, slopeT_);
243  int laymax = laymin + layers_[i];
244  double zz = zi;
245  double thickTot(0);
246  for (int ly = laymin; ly < laymax; ++ly) {
247  int ii = layerType_[ly];
248  int copy = copyNumber_[ii];
249  double hthick = 0.5 * thick_[ii];
250  double rinB = HGCalGeomTools::radius(zo, zFrontB_, rMinFront_, slopeB_);
251  zz += hthick;
252  thickTot += thick_[ii];
253 
255 #ifdef EDM_ML_DEBUG
256  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: Layer " << ly << ":" << ii << " Front "
257  << cms::convert2mm(zi) << ", " << cms::convert2mm(routF) << " Back "
258  << cms::convert2mm(zo) << ", " << cms::convert2mm(rinB)
259  << " superlayer thickness " << cms::convert2mm(layerThick_[i]);
260 #endif
261 
262  dd4hep::Material matter = ns.material(materials_[ii]);
263  dd4hep::Volume glog;
264 
265  if (layerSense_[ly] < 1) {
266  std::vector<double> pgonZ, pgonRin, pgonRout;
267  double rmax =
268  (std::min(routF, HGCalGeomTools::radius(zz + hthick, zFrontT_, rMaxFront_, slopeT_)) * cosAlpha_) - tol1;
269  HGCalGeomTools::radius(zz - hthick,
270  zz + hthick,
271  zFrontB_,
272  rMinFront_,
273  slopeB_,
274  zFrontT_,
275  rMaxFront_,
276  slopeT_,
277  -layerSense_[ly],
278  pgonZ,
279  pgonRin,
280  pgonRout);
281  for (unsigned int isec = 0; isec < pgonZ.size(); ++isec) {
282  pgonZ[isec] -= zz;
283  if (layerSense_[ly] == 0 || absorbMode_ == 0)
284  pgonRout[isec] = rmax;
285  else
286  pgonRout[isec] = pgonRout[isec] * cosAlpha_ - tol1;
287  }
288 
289  dd4hep::Solid solid = dd4hep::Polyhedra(sectors_, -alpha_, 2._pi, pgonZ, pgonRin, pgonRout);
290  ns.addSolidNS(ns.prepend(name), solid);
291  glog = dd4hep::Volume(solid.name(), solid, matter);
292  ns.addVolumeNS(glog);
293 #ifdef EDM_ML_DEBUG
294  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: " << solid.name() << " polyhedra of " << sectors_
295  << " sectors covering " << convertRadToDeg(-alpha_) << ":"
296  << convertRadToDeg(-alpha_ + 2._pi) << " with " << pgonZ.size() << " sections";
297  for (unsigned int k = 0; k < pgonZ.size(); ++k)
298  edm::LogVerbatim("HGCalGeom") << "[" << k << "] z " << cms::convert2mm(pgonZ[k]) << " R "
299  << cms::convert2mm(pgonRin[k]) << ":" << cms::convert2mm(pgonRout[k]);
300 #endif
301  } else {
302  double rins =
304  double routs =
305  (sensitiveMode_ < 1) ? routF : HGCalGeomTools::radius(zz - hthick, zFrontT_, rMaxFront_, slopeT_);
306  dd4hep::Solid solid = dd4hep::Tube(rins, routs, hthick, 0.0, 2._pi);
307  ns.addSolidNS(ns.prepend(name), solid);
308  glog = dd4hep::Volume(solid.name(), solid, matter);
309  ns.addVolumeNS(glog);
310 
311 #ifdef EDM_ML_DEBUG
312  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: " << solid.name() << " Tubs made of "
313  << matter.name() << " of dimensions " << cms::convert2mm(rinB) << ":"
314  << cms::convert2mm(rins) << ", " << cms::convert2mm(routF) << ":"
315  << cms::convert2mm(routs) << ", " << cms::convert2mm(hthick)
316  << ", 0.0, 360.0 and positioned in: " << glog.name() << " number " << copy;
317 #endif
318  positionMix(ctxt, e, glog, name, copy, thick_[ii], matter);
319  }
320 
321  dd4hep::Position r1(0, 0, zz);
322  mother.placeVolume(glog, copy, r1);
323  ++copyNumber_[ii];
324 #ifdef EDM_ML_DEBUG
325  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: " << glog.name() << " number " << copy
326  << " positioned in " << mother.name() << " at (0,0," << cms::convert2mm(zz)
327  << ") with no rotation";
328 #endif
329  zz += hthick;
330  } // End of loop over layers in a block
331  zi = zo;
332  laymin = laymax;
333  if (std::abs(thickTot - layerThick_[i]) > tol2_) {
334  if (thickTot > layerThick_[i]) {
335  edm::LogError("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
336  << " is smaller than " << cms::convert2mm(thickTot)
337  << ": thickness of all its components **** ERROR ****";
338  } else {
339  edm::LogWarning("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
340  << " does not match with " << cms::convert2mm(thickTot) << " of the components";
341  }
342  }
343  } // End of loop over blocks
344 
345 #ifdef EDM_ML_DEBUG
346  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: " << copies_.size() << " different wafer copy numbers";
347  int k(0);
348  for (std::unordered_set<int>::const_iterator itr = copies_.begin(); itr != copies_.end(); ++itr, ++k) {
349  edm::LogVerbatim("HGCalGeom") << "Copy [" << k << "] : " << (*itr);
350  }
351  copies_.clear();
352  edm::LogVerbatim("HGCalGeom") << "<<== End of DDHGCalMixRotatedLayer construction...";
353 #endif
354  }
std::vector< int > tileIndex_
Log< level::Info, true > LogVerbatim
std::vector< int > layerTypeTop_
std::vector< double > thick_
static void radius(double zf, double zb, std::vector< double > const &zFront1, std::vector< double > const &rFront1, std::vector< double > const &slope1, std::vector< double > const &zFront2, std::vector< double > const &rFront2, std::vector< double > const &slope2, int flag, std::vector< double > &zz, std::vector< double > &rin, std::vector< double > &rout)
std::vector< int > waferIndex_
int32_t waferU(const int32_t index)
int32_t waferLayer(const int32_t index)
std::vector< int > layerSense_
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
constexpr NumType convert2mm(NumType length)
Definition: DDutils.h:7
std::vector< int > tileLayerStart_
std::string to_string(const V &value)
Definition: OMSAccess.h:77
std::vector< int > copyNumber_
std::vector< int > copyNumberTop_
std::vector< double > cassetteShift_
Log< level::Error, false > LogError
static constexpr double tol2_
std::vector< double > zFrontB_
int32_t waferOrient(const int32_t property)
std::vector< int > waferProperty_
std::vector< double > layerThickTop_
std::vector< std::string > materialTop_
std::vector< double > rMaxFront_
std::vector< std::string > namesTop_
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
std::vector< int > tilePhis_
std::vector< double > tileRMax_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::vector< int > waferLayerStart_
std::vector< double > layerThick_
int32_t waferThick(const int32_t property)
dd4hep::Volume Volume
std::vector< double > slopeT_
ii
Definition: cuy.py:589
std::unordered_set< int > copies_
std::tuple< int32_t, int32_t, int32_t > tileUnpack(int32_t index)
std::vector< double > zFrontT_
std::vector< int > layerOrient_
std::vector< std::string > waferPart_
void positionMix(cms::DDParsingContext &ctxt, xml_h e, const dd4hep::Volume &glog, const std::string &nameM, int copyM, double thick, const dd4hep::Material &matter)
std::vector< int > layerType_
std::vector< double > tileRMin_
int32_t waferPartial(const int32_t property)
std::vector< double > slopeB_
int32_t waferV(const int32_t index)
std::vector< std::string > waferFull_
Log< level::Warning, false > LogWarning
static int32_t layerType(int type)
Definition: HGCalTypes.cc:42
std::vector< std::string > names_
std::vector< double > rMinFront_
std::vector< std::string > materials_
void setParameter(int cassette, const std::vector< double > &shifts)
Definition: HGCalCassette.cc:8

Member Function Documentation

◆ positionMix()

void HGCalMixRotatedLayer::positionMix ( cms::DDParsingContext ctxt,
xml_h  e,
const dd4hep::Volume &  glog,
const std::string &  nameM,
int  copyM,
double  thick,
const dd4hep::Material &  matter 
)
inline

Definition at line 356 of file DDHGCalMixRotatedLayer.cc.

References funct::abs(), cms::DDNamespace::addSolidNS(), cms::DDNamespace::addVolumeNS(), angle(), HGCalCell::cellPlacementIndex(), cms::convert2mm(), angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, PVValHelper::dy, MillePedeFileConverter_cfg::e, mps_fire::i, cuy::ii, gpuVertexFinder::iv, dqmdumpme::k, HGCalTypes::layerFrontBack(), cms::DDNamespace::material(), g4SimHits_cfi::Material, Skims_PA_cff::name, EgHLTOffHistBins_cfi::nr, HGCalTypes::packTypeUV(), PixelTestBeamValidation_cfi::Position, cms::DDNamespace::prepend(), diffTwoXMLs::r2, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, HGCalTileIndex::tileUnpack(), cond::impl::to_string(), findQualityFiles::v, cms::DDNamespace::volume(), HGCalProperty::waferCassette(), HGCalTypes::WaferFull, HGCalTypes::WaferHDTop, HGCalWaferIndex::waferLayer(), HGCalProperty::waferOrient(), HGCalTypes::WaferPartHDOffset, HGCalProperty::waferPartial(), HGCalTypes::WaferPartLDOffset, HGCalProperty::waferThick(), HGCalTypes::WaferTypeOffset, HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

362  {
363  cms::DDNamespace ns(ctxt, e, true);
364 
365  // Make the top part first
366  for (unsigned int ly = 0; ly < layerTypeTop_.size(); ++ly) {
367  int ii = layerTypeTop_[ly];
368  copyNumberTop_[ii] = copyM;
369  }
370  double hthick = 0.5 * thick;
371  double dphi = (2._pi) / phiBinsScint_;
372  double thickTot(0), zpos(-hthick);
373  for (unsigned int ly = 0; ly < layerTypeTop_.size(); ++ly) {
374  int ii = layerTypeTop_[ly];
375  int copy = copyNumberTop_[ii];
376  int layer = copy - firstLayer_;
377  double hthickl = 0.5 * layerThickTop_[ii];
378  thickTot += layerThickTop_[ii];
379  zpos += hthickl;
380  dd4hep::Material matter1 = ns.material(materialTop_[ii]);
381  unsigned int k = 0;
382  int firstTile = tileLayerStart_[layer];
383  int lastTile = ((layer + 1 < static_cast<int>(tileLayerStart_.size())) ? tileLayerStart_[layer + 1]
384  : static_cast<int>(tileIndex_.size()));
385 #ifdef EDM_ML_DEBUG
386  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: Layer " << ly << ":" << ii << " Copy " << copy
387  << " Tiles " << firstTile << ":" << lastTile;
388 #endif
389  for (int ti = firstTile; ti < lastTile; ++ti) {
390  double r1 = tileRMin_[std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
391  double r2 = tileRMax_[std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
392  int cassette = std::get<0>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
393  int fimin = std::get<1>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
394  int fimax = std::get<2>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
395  double phi1 = dphi * (fimin - 1);
396  double phi2 = (forFireworks_ == 1) ? (dphi * (fimax - fimin + 1)) : (dphi * fimax);
397  auto cshift = cassette_.getShift(layer + 1, 1, cassette);
398 #ifdef EDM_ML_DEBUG
399  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: Layer " << copy << ":" << (layer + 1) << " iR "
400  << std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) << ":"
401  << std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) << " R "
402  << cms::convert2mm(r1) << ":" << cms::convert2mm(r2) << " Thick "
403  << cms::convert2mm((2.0 * hthickl)) << " phi " << fimin << ":" << fimax << ":"
404  << convertRadToDeg(phi1) << ":" << convertRadToDeg(phi2) << " cassette "
405  << cassette << " Shift " << cshift.first << ":" << cshift.second;
406 #endif
408  ++k;
409  dd4hep::Solid solid = dd4hep::Tube(r1, r2, hthickl, phi1, phi2);
410  ns.addSolidNS(ns.prepend(name), solid);
411  dd4hep::Volume glog1 = dd4hep::Volume(solid.name(), solid, matter1);
412  ns.addVolumeNS(glog1);
413 #ifdef EDM_ML_DEBUG
414  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: " << glog1.name() << " Tubs made of "
415  << materialTop_[ii] << " of dimensions " << cms::convert2mm(r1) << ", "
416  << cms::convert2mm(r2) << ", " << cms::convert2mm(hthickl) << ", "
417  << convertRadToDeg(phi1) << ", " << convertRadToDeg(phi2);
418 #endif
419  dd4hep::Position tran(-cshift.first, cshift.second, zpos);
420  glog.placeVolume(glog1, copy, tran);
421 #ifdef EDM_ML_DEBUG
422  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: Position " << glog1.name() << " number " << copy
423  << " in " << glog.name() << " at (" << cms::convert2mm(cshift.first) << ", "
424  << cms::convert2mm(cshift.second) << ", " << cms::convert2mm(zpos)
425  << ") with no rotation";
426 #endif
427  }
428  ++copyNumberTop_[ii];
429  zpos += hthickl;
430  }
431  if (std::abs(thickTot - thick) > tol2_) {
432  if (thickTot > thick) {
433  edm::LogError("HGCalGeom") << "DDHGCalMixRotatedLayer: Thickness of the partition " << cms::convert2mm(thick)
434  << " is smaller than " << cms::convert2mm(thickTot)
435  << ": thickness of all its components in the top part **** ERROR ****";
436  } else {
437  edm::LogWarning("HGCalGeom") << "DDHGCalMixRotatedLayer: Thickness of the partition " << cms::convert2mm(thick)
438  << " does not match with " << cms::convert2mm(thickTot)
439  << " of the components in top part";
440  }
441  }
442 
443  // Make the bottom part next
444  int layer = (copyM - firstLayer_);
445  static const double sqrt3 = std::sqrt(3.0);
446  int layercenter = layerOrient_[layer];
448  int firstWafer = waferLayerStart_[layer];
449  int lastWafer = ((layer + 1 < static_cast<int>(waferLayerStart_.size())) ? waferLayerStart_[layer + 1]
450  : static_cast<int>(waferIndex_.size()));
451  double delx = 0.5 * (waferSize_ + waferSepar_);
452  double dely = 2.0 * delx / sqrt3;
453  double dy = 0.75 * dely;
454  const auto& xyoff = geomTools_.shiftXY(layercenter, (waferSize_ + waferSepar_));
455 #ifdef EDM_ML_DEBUG
456  int ium(0), ivm(0), kount(0);
457  std::vector<int> ntype(3, 0);
458  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: " << glog.name() << " r " << cms::convert2mm(delx)
459  << " R " << cms::convert2mm(dely) << " dy " << cms::convert2mm(dy) << " Shift "
460  << cms::convert2mm(xyoff.first) << ":" << cms::convert2mm(xyoff.second)
461  << " WaferSize " << cms::convert2mm((waferSize_ + waferSepar_)) << " index "
462  << firstWafer << ":" << (lastWafer - 1);
463 #endif
464  for (int k = firstWafer; k < lastWafer; ++k) {
467 #ifdef EDM_ML_DEBUG
468  int iu = std::abs(u);
469  int iv = std::abs(v);
470 #endif
471  int nr = 2 * v;
472  int nc = -2 * u + v;
477  int place = HGCalCell::cellPlacementIndex(1, layertype, orien);
478 #ifdef EDM_ML_DEBUG
479  edm::LogVerbatim("HGCalGeom")
480  << "DDHGCalMixRotatedLayer::index:Property:layertype:type:part:orien:cassette:place:offsets:ind " << k << ":"
481  << waferProperty_[k] << ":" << layertype << ":" << type << ":" << part << ":" << orien << ":" << cassette
482  << ":" << place;
483 #endif
484  auto cshift = cassette_.getShift(layer + 1, -1, cassette);
485  double xpos = xyoff.first - cshift.first + nc * delx;
486  double ypos = xyoff.second + cshift.second + nr * dy;
487 #ifdef EDM_ML_DEBUG
488  double xorig = xyoff.first + nc * delx;
489  double yorig = xyoff.second + nr * dy;
490  double angle = std::atan2(yorig, xorig);
491  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer::Wafer: layer " << layer + 1 << " cassette " << cassette
492  << " Shift " << cshift.first << ":" << cshift.second << " Original " << xorig << ":"
493  << yorig << ":" << convertRadToDeg(angle) << " Final " << xpos << ":" << ypos;
494 #endif
495  std::string wafer;
496  int i(999);
497  if (part == HGCalTypes::WaferFull) {
499 #ifdef EDM_ML_DEBUG
500  edm::LogVerbatim("HGCalGeom") << " FullWafer type:place:ind " << type << ":" << place << ":" << i << ":"
501  << waferFull_.size();
502 #endif
503  wafer = waferFull_[i];
504  } else {
505  int partoffset =
507  i = (part - partoffset) * facingTypes_ * orientationTypes_ +
509 #ifdef EDM_ML_DEBUG
510  edm::LogVerbatim("HGCalGeom") << " layertype:type:part:orien:cassette:place:offsets:ind " << layertype << ":"
511  << type << ":" << part << ":" << orien << ":" << cassette << ":" << place << ":"
512  << partoffset << ":" << HGCalTypes::WaferTypeOffset[type] << ":" << i << ":"
513  << waferPart_.size();
514 #endif
515  wafer = waferPart_[i];
516  }
517  int copy = HGCalTypes::packTypeUV(type, u, v);
518 #ifdef EDM_ML_DEBUG
519  edm::LogVerbatim("HGCalGeom") << " DDHGCalMixRotatedLayer: Layer " << HGCalWaferIndex::waferLayer(waferIndex_[k])
520  << " Wafer " << wafer << " number " << copy << " type :part:orien:ind " << type
521  << ":" << part << ":" << orien << ":" << i << " layer:u:v " << (layer + firstLayer_)
522  << ":" << u << ":" << v;
523  if (iu > ium)
524  ium = iu;
525  if (iv > ivm)
526  ivm = iv;
527  kount++;
528  if (copies_.count(copy) == 0)
529  copies_.insert(copy);
530 #endif
531  dd4hep::Position tran(xpos, ypos, 0.0);
532  glog.placeVolume(ns.volume(wafer), copy, tran);
533 #ifdef EDM_ML_DEBUG
534  ++ntype[type];
535  edm::LogVerbatim("HGCalGeom") << " DDHGCalMixRotatedLayer: " << wafer << " number " << copy << " type "
536  << layertype << ":" << type << " positioned in " << glog.name() << " at ("
537  << cms::convert2mm(xpos) << "," << cms::convert2mm(ypos) << ",0) with no rotation";
538 #endif
539  }
540 
541 #ifdef EDM_ML_DEBUG
542  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedLayer: Maximum # of u " << ium << " # of v " << ivm << " and "
543  << kount << " wafers (" << ntype[0] << ":" << ntype[1] << ":" << ntype[2] << ") for "
544  << glog.name();
545 #endif
546  }
std::vector< int > tileIndex_
Log< level::Info, true > LogVerbatim
static int32_t cellPlacementIndex(int32_t iz, int32_t frontBack, int32_t orient)
Definition: HGCalCell.cc:233
std::vector< int > layerTypeTop_
static constexpr int32_t WaferPartLDOffset
Definition: HGCalTypes.h:57
int32_t *__restrict__ iv
static constexpr int32_t WaferTypeOffset[3]
Definition: HGCalTypes.h:59
std::vector< int > waferIndex_
int32_t waferU(const int32_t index)
int32_t waferLayer(const int32_t index)
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
constexpr NumType convert2mm(NumType length)
Definition: DDutils.h:7
std::vector< int > tileLayerStart_
std::pair< double, double > getShift(int layer, int zside, int cassette) const
std::string to_string(const V &value)
Definition: OMSAccess.h:77
std::vector< int > copyNumberTop_
Log< level::Error, false > LogError
static constexpr double tol2_
int32_t waferOrient(const int32_t property)
std::pair< double, double > shiftXY(int waferPosition, double waferSize) const
std::vector< int > waferProperty_
std::vector< double > layerThickTop_
constexpr std::array< uint8_t, layerIndexSize< TrackerTraits > > layer
std::vector< std::string > materialTop_
int32_t waferCassette(const int32_t property)
std::vector< std::string > namesTop_
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:35
T sqrt(T t)
Definition: SSEVec.h:19
std::vector< int > tilePhis_
std::vector< double > tileRMax_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::vector< int > waferLayerStart_
int32_t waferThick(const int32_t property)
dd4hep::Volume Volume
ii
Definition: cuy.py:589
std::unordered_set< int > copies_
std::tuple< int32_t, int32_t, int32_t > tileUnpack(int32_t index)
static constexpr int32_t WaferPartHDOffset
Definition: HGCalTypes.h:58
std::vector< int > layerOrient_
std::vector< std::string > waferPart_
part
Definition: HCALResponse.h:20
std::vector< double > tileRMin_
int32_t waferPartial(const int32_t property)
int32_t waferV(const int32_t index)
std::vector< std::string > waferFull_
static int32_t packTypeUV(int type, int u, int v)
Definition: HGCalTypes.cc:3
static constexpr int32_t WaferHDTop
Definition: HGCalTypes.h:51
Log< level::Warning, false > LogWarning
static int32_t layerFrontBack(int32_t layerOrient)
Definition: HGCalTypes.h:125
T angle(T x1, T y1, T z1, T x2, T y2, T z2)
Definition: angle.h:11

Member Data Documentation

◆ absorbMode_

int HGCalMixRotatedLayer::absorbMode_

Definition at line 561 of file DDHGCalMixRotatedLayer.cc.

◆ alpha_

double HGCalMixRotatedLayer::alpha_

Definition at line 600 of file DDHGCalMixRotatedLayer.cc.

◆ cassette_

HGCalCassette HGCalMixRotatedLayer::cassette_

Definition at line 550 of file DDHGCalMixRotatedLayer.cc.

◆ cassettes_

int HGCalMixRotatedLayer::cassettes_

Definition at line 567 of file DDHGCalMixRotatedLayer.cc.

◆ cassetteShift_

std::vector<double> HGCalMixRotatedLayer::cassetteShift_

Definition at line 593 of file DDHGCalMixRotatedLayer.cc.

◆ copies_

std::unordered_set<int> HGCalMixRotatedLayer::copies_

Definition at line 599 of file DDHGCalMixRotatedLayer.cc.

◆ copyNumber_

std::vector<int> HGCalMixRotatedLayer::copyNumber_

Definition at line 579 of file DDHGCalMixRotatedLayer.cc.

◆ copyNumberTop_

std::vector<int> HGCalMixRotatedLayer::copyNumberTop_

Definition at line 588 of file DDHGCalMixRotatedLayer.cc.

◆ cosAlpha_

double HGCalMixRotatedLayer::cosAlpha_

Definition at line 600 of file DDHGCalMixRotatedLayer.cc.

◆ facingTypes_

int HGCalMixRotatedLayer::facingTypes_

Definition at line 555 of file DDHGCalMixRotatedLayer.cc.

◆ firstLayer_

int HGCalMixRotatedLayer::firstLayer_

Definition at line 560 of file DDHGCalMixRotatedLayer.cc.

◆ forFireworks_

int HGCalMixRotatedLayer::forFireworks_

Definition at line 559 of file DDHGCalMixRotatedLayer.cc.

◆ geomTools_

HGCalGeomTools HGCalMixRotatedLayer::geomTools_

Definition at line 549 of file DDHGCalMixRotatedLayer.cc.

◆ layerOrient_

std::vector<int> HGCalMixRotatedLayer::layerOrient_

Definition at line 589 of file DDHGCalMixRotatedLayer.cc.

◆ layers_

std::vector<int> HGCalMixRotatedLayer::layers_

Definition at line 580 of file DDHGCalMixRotatedLayer.cc.

◆ layerSense_

std::vector<int> HGCalMixRotatedLayer::layerSense_

Definition at line 583 of file DDHGCalMixRotatedLayer.cc.

◆ layerThick_

std::vector<double> HGCalMixRotatedLayer::layerThick_

Definition at line 581 of file DDHGCalMixRotatedLayer.cc.

◆ layerThickTop_

std::vector<double> HGCalMixRotatedLayer::layerThickTop_

Definition at line 586 of file DDHGCalMixRotatedLayer.cc.

◆ layerType_

std::vector<int> HGCalMixRotatedLayer::layerType_

Definition at line 582 of file DDHGCalMixRotatedLayer.cc.

◆ layerTypeTop_

std::vector<int> HGCalMixRotatedLayer::layerTypeTop_

Definition at line 587 of file DDHGCalMixRotatedLayer.cc.

◆ materials_

std::vector<std::string> HGCalMixRotatedLayer::materials_

Definition at line 576 of file DDHGCalMixRotatedLayer.cc.

◆ materialTop_

std::vector<std::string> HGCalMixRotatedLayer::materialTop_

Definition at line 584 of file DDHGCalMixRotatedLayer.cc.

◆ names_

std::vector<std::string> HGCalMixRotatedLayer::names_

Definition at line 577 of file DDHGCalMixRotatedLayer.cc.

◆ namesTop_

std::vector<std::string> HGCalMixRotatedLayer::namesTop_

Definition at line 585 of file DDHGCalMixRotatedLayer.cc.

◆ orientationTypes_

int HGCalMixRotatedLayer::orientationTypes_

Definition at line 556 of file DDHGCalMixRotatedLayer.cc.

◆ phiBinsScint_

int HGCalMixRotatedLayer::phiBinsScint_

Definition at line 558 of file DDHGCalMixRotatedLayer.cc.

◆ placeOffset_

int HGCalMixRotatedLayer::placeOffset_

Definition at line 557 of file DDHGCalMixRotatedLayer.cc.

◆ rMaxFront_

std::vector<double> HGCalMixRotatedLayer::rMaxFront_

Definition at line 573 of file DDHGCalMixRotatedLayer.cc.

◆ rMinFront_

std::vector<double> HGCalMixRotatedLayer::rMinFront_

Definition at line 570 of file DDHGCalMixRotatedLayer.cc.

◆ sectors_

int HGCalMixRotatedLayer::sectors_

Definition at line 566 of file DDHGCalMixRotatedLayer.cc.

◆ sensitiveMode_

int HGCalMixRotatedLayer::sensitiveMode_

Definition at line 562 of file DDHGCalMixRotatedLayer.cc.

◆ slopeB_

std::vector<double> HGCalMixRotatedLayer::slopeB_

Definition at line 568 of file DDHGCalMixRotatedLayer.cc.

◆ slopeT_

std::vector<double> HGCalMixRotatedLayer::slopeT_

Definition at line 571 of file DDHGCalMixRotatedLayer.cc.

◆ thick_

std::vector<double> HGCalMixRotatedLayer::thick_

Definition at line 578 of file DDHGCalMixRotatedLayer.cc.

◆ tileIndex_

std::vector<int> HGCalMixRotatedLayer::tileIndex_

Definition at line 596 of file DDHGCalMixRotatedLayer.cc.

◆ tileLayerStart_

std::vector<int> HGCalMixRotatedLayer::tileLayerStart_

Definition at line 598 of file DDHGCalMixRotatedLayer.cc.

◆ tilePhis_

std::vector<int> HGCalMixRotatedLayer::tilePhis_

Definition at line 597 of file DDHGCalMixRotatedLayer.cc.

◆ tileRMax_

std::vector<double> HGCalMixRotatedLayer::tileRMax_

Definition at line 595 of file DDHGCalMixRotatedLayer.cc.

◆ tileRMin_

std::vector<double> HGCalMixRotatedLayer::tileRMin_

Definition at line 594 of file DDHGCalMixRotatedLayer.cc.

◆ tol2_

constexpr double HGCalMixRotatedLayer::tol2_ = 0.00001 * dd4hep::mm
static

Definition at line 552 of file DDHGCalMixRotatedLayer.cc.

◆ waferFull_

std::vector<std::string> HGCalMixRotatedLayer::waferFull_

Definition at line 574 of file DDHGCalMixRotatedLayer.cc.

◆ waferIndex_

std::vector<int> HGCalMixRotatedLayer::waferIndex_

Definition at line 590 of file DDHGCalMixRotatedLayer.cc.

◆ waferLayerStart_

std::vector<int> HGCalMixRotatedLayer::waferLayerStart_

Definition at line 592 of file DDHGCalMixRotatedLayer.cc.

◆ waferPart_

std::vector<std::string> HGCalMixRotatedLayer::waferPart_

Definition at line 575 of file DDHGCalMixRotatedLayer.cc.

◆ waferProperty_

std::vector<int> HGCalMixRotatedLayer::waferProperty_

Definition at line 591 of file DDHGCalMixRotatedLayer.cc.

◆ waferSepar_

double HGCalMixRotatedLayer::waferSepar_

Definition at line 565 of file DDHGCalMixRotatedLayer.cc.

◆ waferSize_

double HGCalMixRotatedLayer::waferSize_

Definition at line 564 of file DDHGCalMixRotatedLayer.cc.

◆ waferTypes_

int HGCalMixRotatedLayer::waferTypes_

Definition at line 554 of file DDHGCalMixRotatedLayer.cc.

◆ zFrontB_

std::vector<double> HGCalMixRotatedLayer::zFrontB_

Definition at line 569 of file DDHGCalMixRotatedLayer.cc.

◆ zFrontT_

std::vector<double> HGCalMixRotatedLayer::zFrontT_

Definition at line 572 of file DDHGCalMixRotatedLayer.cc.

◆ zMinBlock_

double HGCalMixRotatedLayer::zMinBlock_

Definition at line 563 of file DDHGCalMixRotatedLayer.cc.