CMS 3D CMS Logo
DDHGCalMixLayer.cc
Go to the documentation of this file.
1 // File: DDHGCalMixLayer.cc
3 // Description: Geometry factory class for HGCal (Mix) adopted for DD4HEP
5 
6 #include <cmath>
7 #include <memory>
8 #include <string>
9 #include <unordered_set>
10 #include <vector>
11 
12 #include "Geometry/HGCalCommonData/interface/HGCalGeomTools.h"
13 #include "Geometry/HGCalCommonData/interface/HGCalParameters.h"
14 #include "Geometry/HGCalCommonData/interface/HGCalProperty.h"
15 #include "Geometry/HGCalCommonData/interface/HGCalTileIndex.h"
16 #include "Geometry/HGCalCommonData/interface/HGCalTypes.h"
17 #include "Geometry/HGCalCommonData/interface/HGCalWaferIndex.h"
18 #include "Geometry/HGCalCommonData/interface/HGCalWaferType.h"
19 #include "DD4hep/DetFactoryHelper.h"
20 #include "DataFormats/Math/interface/angle_units.h"
21 #include "DetectorDescription/DDCMS/interface/DDPlugins.h"
22 #include "DetectorDescription/DDCMS/interface/DDutils.h"
23 #include "FWCore/MessageLogger/interface/MessageLogger.h"
24 
25 //#define EDM_ML_DEBUG
26 using namespace angle_units::operators;
27 
28 struct HGCalMixLayer {
29  HGCalMixLayer() { throw cms::Exception("HGCalGeom") << "Wrong initialization to HGCalMixLayer"; }
30  HGCalMixLayer(cms::DDParsingContext& ctxt, xml_h e) {
31  cms::DDNamespace ns(ctxt, e, true);
32  cms::DDAlgoArguments args(ctxt, e);
33 
34 #ifdef EDM_ML_DEBUG
35  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: Creating an instance";
36 #endif
37 
38  static constexpr double tol1 = 0.01 * dd4hep::mm;
39  dd4hep::Volume mother = ns.volume(args.parentName());
40 
41  waferTypes_ = args.value<int>("WaferTypes");
42  facingTypes_ = args.value<int>("FacingTypes");
43  partialTypes_ = args.value<int>("PartialTypes");
44  orientationTypes_ = args.value<int>("OrientationTypes");
45  phiBinsScint_ = args.value<int>("NPhiBinScint");
46 #ifdef EDM_ML_DEBUG
47  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer::Number of types of wafers: " << waferTypes_
48  << " facings: " << facingTypes_ << " partials: " << partialTypes_
49  << " Orientations: " << orientationTypes_ << "; number of cells along phi "
50  << phiBinsScint_;
51 #endif
52  firstLayer_ = args.value<int>("FirstLayer");
53  absorbMode_ = args.value<int>("AbsorberMode");
54  sensitiveMode_ = args.value<int>("SensitiveMode");
55 #ifdef EDM_ML_DEBUG
56  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer::First Layer " << firstLayer_ << " and "
57  << "Absober:Sensitive mode " << absorbMode_ << ":" << sensitiveMode_;
58 #endif
59  zMinBlock_ = args.value<double>("zMinBlock");
60  waferSize_ = args.value<double>("waferSize");
61  waferSepar_ = args.value<double>("SensorSeparation");
62  sectors_ = args.value<int>("Sectors");
63  alpha_ = (1._pi) / sectors_;
64  cosAlpha_ = cos(alpha_);
65 #ifdef EDM_ML_DEBUG
66  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: zStart " << cms::convert2mm(zMinBlock_) << " wafer width "
67  << cms::convert2mm(waferSize_) << " separations " << cms::convert2mm(waferSepar_)
68  << " sectors " << sectors_ << ":" << convertRadToDeg(alpha_) << ":" << cosAlpha_;
69 #endif
70  slopeB_ = args.value<std::vector<double>>("SlopeBottom");
71  zFrontB_ = args.value<std::vector<double>>("ZFrontBottom");
72  rMinFront_ = args.value<std::vector<double>>("RMinFront");
73  slopeT_ = args.value<std::vector<double>>("SlopeTop");
74  zFrontT_ = args.value<std::vector<double>>("ZFrontTop");
75  rMaxFront_ = args.value<std::vector<double>>("RMaxFront");
76 #ifdef EDM_ML_DEBUG
77  for (unsigned int i = 0; i < slopeB_.size(); ++i)
78  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] Zmin " << cms::convert2mm(zFrontB_[i]) << " Rmin "
79  << cms::convert2mm(rMinFront_[i]) << " Slope " << slopeB_[i];
80  for (unsigned int i = 0; i < slopeT_.size(); ++i)
81  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] Zmin " << cms::convert2mm(zFrontT_[i]) << " Rmax "
82  << cms::convert2mm(rMaxFront_[i]) << " Slope " << slopeT_[i];
83 #endif
84 
85  waferFull_ = args.value<std::vector<std::string>>("WaferNamesFull");
86  waferPart_ = args.value<std::vector<std::string>>("WaferNamesPartial");
87 #ifdef EDM_ML_DEBUG
88  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: " << waferFull_.size() << " full and " << waferPart_.size()
89  << " partial modules\nDDHGCalMixLayer:Full Modules:";
90  unsigned int i1max = static_cast<unsigned int>(waferFull_.size());
91  for (unsigned int i1 = 0; i1 < i1max; i1 += 2) {
92  std::ostringstream st1;
93  unsigned int i2 = std::min((i1 + 2), i1max);
94  for (unsigned int i = i1; i < i2; ++i)
95  st1 << " [" << i << "] " << waferFull_[i];
96  edm::LogVerbatim("HGCalGeom") << st1.str() << std::endl;
97  }
98  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: Partial Modules:";
99  i1max = static_cast<unsigned int>(waferPart_.size());
100  for (unsigned int i1 = 0; i1 < i1max; i1 += 2) {
101  std::ostringstream st1;
102  unsigned int i2 = std::min((i1 + 2), i1max);
103  for (unsigned int i = i1; i < i2; ++i)
104  st1 << " [" << i << "] " << waferPart_[i];
105  edm::LogVerbatim("HGCalGeom") << st1.str() << std::endl;
106  }
107 #endif
108 
109  materials_ = args.value<std::vector<std::string>>("MaterialNames");
110  names_ = args.value<std::vector<std::string>>("VolumeNames");
111  thick_ = args.value<std::vector<double>>("Thickness");
112  copyNumber_.resize(materials_.size(), 1);
113 #ifdef EDM_ML_DEBUG
114  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: " << materials_.size() << " types of volumes";
115  for (unsigned int i = 0; i < names_.size(); ++i)
116  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << names_[i] << " of thickness "
117  << cms::convert2mm(thick_[i]) << " filled with " << materials_[i]
118  << " first copy number " << copyNumber_[i];
119 #endif
120  layers_ = args.value<std::vector<int>>("Layers");
121  layerThick_ = args.value<std::vector<double>>("LayerThick");
122 #ifdef EDM_ML_DEBUG
123  edm::LogVerbatim("HGCalGeom") << "There are " << layers_.size() << " blocks";
124  for (unsigned int i = 0; i < layers_.size(); ++i)
125  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] of thickness " << cms::convert2mm(layerThick_[i])
126  << " with " << layers_[i] << " layers";
127 #endif
128  layerType_ = args.value<std::vector<int>>("LayerType");
129  layerSense_ = args.value<std::vector<int>>("LayerSense");
130  layerTypes_ = args.value<std::vector<int>>("LayerTypes");
131 #ifdef EDM_ML_DEBUG
132  for (unsigned int i = 0; i < layerTypes_.size(); ++i)
133  edm::LogVerbatim("HGCalGeom") << "LayerTypes [" << i << "] " << layerTypes_[i];
134 #endif
135  if (firstLayer_ > 0) {
136  for (unsigned int i = 0; i < layerType_.size(); ++i) {
137  if (layerSense_[i] > 0) {
138  int ii = layerType_[i];
139  copyNumber_[ii] = firstLayer_;
140 #ifdef EDM_ML_DEBUG
141  edm::LogVerbatim("HGCalGeom") << "First copy number for layer type " << i << ":" << ii << " with "
142  << materials_[ii] << " changed to " << copyNumber_[ii];
143 #endif
144  break;
145  }
146  }
147  } else {
148  firstLayer_ = 1;
149  }
150 #ifdef EDM_ML_DEBUG
151  edm::LogVerbatim("HGCalGeom") << "There are " << layerType_.size() << " layers";
152  for (unsigned int i = 0; i < layerType_.size(); ++i)
153  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerType_[i] << " sensitive class "
154  << layerSense_[i];
155 #endif
156  materialTop_ = args.value<std::vector<std::string>>("TopMaterialNames");
157  namesTop_ = args.value<std::vector<std::string>>("TopVolumeNames");
158  layerThickTop_ = args.value<std::vector<double>>("TopLayerThickness");
159  layerTypeTop_ = args.value<std::vector<int>>("TopLayerType");
160  copyNumberTop_.resize(materialTop_.size(), firstLayer_);
161 #ifdef EDM_ML_DEBUG
162  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: " << materialTop_.size() << " types of volumes in the top part";
163  for (unsigned int i = 0; i < materialTop_.size(); ++i)
164  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << namesTop_[i] << " of thickness "
165  << cms::convert2mm(layerThickTop_[i]) << " filled with " << materialTop_[i]
166  << " first copy number " << copyNumberTop_[i];
167  edm::LogVerbatim("HGCalGeom") << "There are " << layerTypeTop_.size() << " layers in the top part";
168  for (unsigned int i = 0; i < layerTypeTop_.size(); ++i)
169  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerTypeTop_[i];
170 #endif
171  waferIndex_ = args.value<std::vector<int>>("WaferIndex");
172  waferProperty_ = args.value<std::vector<int>>("WaferProperties");
173  waferLayerStart_ = args.value<std::vector<int>>("WaferLayerStart");
174 #ifdef EDM_ML_DEBUG
175  edm::LogVerbatim("HGCalGeom") << "waferProperties with " << waferIndex_.size() << " entries in "
176  << waferLayerStart_.size() << " layers";
177  for (unsigned int k = 0; k < waferLayerStart_.size(); ++k)
178  edm::LogVerbatim("HGCalGeom") << "LayerStart[" << k << "] " << waferLayerStart_[k];
179  for (unsigned int k = 0; k < waferIndex_.size(); ++k)
180  edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << waferIndex_[k] << " ("
181  << HGCalWaferIndex::waferLayer(waferIndex_[k]) << ", "
182  << HGCalWaferIndex::waferU(waferIndex_[k]) << ", "
183  << HGCalWaferIndex::waferV(waferIndex_[k]) << ") : ("
184  << HGCalProperty::waferThick(waferProperty_[k]) << ":"
185  << HGCalProperty::waferPartial(waferProperty_[k]) << ":"
186  << HGCalProperty::waferOrient(waferProperty_[k]) << ")";
187 #endif
188  tileRMin_ = args.value<std::vector<double>>("TileRMin");
189  tileRMax_ = args.value<std::vector<double>>("TileRMax");
190  tileIndex_ = args.value<std::vector<int>>("TileLayerRings");
191  tilePhis_ = args.value<std::vector<int>>("TilePhiRange");
192  tileLayerStart_ = args.value<std::vector<int>>("TileLayerStart");
193 #ifdef EDM_ML_DEBUG
194  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer:: with " << tileRMin_.size() << " rings";
195  for (unsigned int k = 0; k < tileRMin_.size(); ++k)
196  edm::LogVerbatim("HGCalGeom") << "Ring[" << k << "] " << cms::convert2mm(tileRMin_[k]) << " : "
197  << cms::convert2mm(tileRMax_[k]);
198  edm::LogVerbatim("HGCalGeom") << "TileProperties with " << tileIndex_.size() << " entries in "
199  << tileLayerStart_.size() << " layers";
200  for (unsigned int k = 0; k < tileLayerStart_.size(); ++k)
201  edm::LogVerbatim("HGCalGeom") << "LayerStart[" << k << "] " << tileLayerStart_[k];
202  for (unsigned int k = 0; k < tileIndex_.size(); ++k)
203  edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << tileIndex_[k] << " ("
204  << "Layer " << std::get<0>(HGCalTileIndex::tileUnpack(tileIndex_[k])) << " Ring "
205  << std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[k])) << ":"
206  << std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[k])) << ") Phi "
207  << std::get<1>(HGCalTileIndex::tileUnpack(tilePhis_[k])) << ":"
208  << std::get<2>(HGCalTileIndex::tileUnpack(tilePhis_[k]));
209 #endif
210 
211 #ifdef EDM_ML_DEBUG
212  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: NameSpace " << ns.name();
213 
214  edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalMixLayer...";
215  copies_.clear();
216 #endif
217 
218  double zi(zMinBlock_);
219  int laymin(0);
220  for (unsigned int i = 0; i < layers_.size(); i++) {
221  double zo = zi + layerThick_[i];
222  double routF = HGCalGeomTools::radius(zi, zFrontT_, rMaxFront_, slopeT_);
223  int laymax = laymin + layers_[i];
224  double zz = zi;
225  double thickTot(0);
226  for (int ly = laymin; ly < laymax; ++ly) {
227  int ii = layerType_[ly];
228  int copy = copyNumber_[ii];
229  double hthick = 0.5 * thick_[ii];
230  double rinB = HGCalGeomTools::radius(zo, zFrontB_, rMinFront_, slopeB_);
231  zz += hthick;
232  thickTot += thick_[ii];
233 
234  std::string name = names_[ii] + std::to_string(copy);
235 #ifdef EDM_ML_DEBUG
236  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: Layer " << ly << ":" << ii << " Front "
237  << cms::convert2mm(zi) << ", " << cms::convert2mm(routF) << " Back "
238  << cms::convert2mm(zo) << ", " << cms::convert2mm(rinB)
239  << " superlayer thickness " << cms::convert2mm(layerThick_[i]);
240 #endif
241 
242  dd4hep::Material matter = ns.material(materials_[ii]);
243  dd4hep::Volume glog;
244 
245  if (layerSense_[ly] < 1) {
246  std::vector<double> pgonZ, pgonRin, pgonRout;
247  double rmax =
248  (std::min(routF, HGCalGeomTools::radius(zz + hthick, zFrontT_, rMaxFront_, slopeT_)) * cosAlpha_) - tol1;
249  HGCalGeomTools::radius(zz - hthick,
250  zz + hthick,
251  zFrontB_,
252  rMinFront_,
253  slopeB_,
254  zFrontT_,
255  rMaxFront_,
256  slopeT_,
257  -layerSense_[ly],
258  pgonZ,
259  pgonRin,
260  pgonRout);
261  for (unsigned int isec = 0; isec < pgonZ.size(); ++isec) {
262  pgonZ[isec] -= zz;
263  if (layerSense_[ly] == 0 || absorbMode_ == 0)
264  pgonRout[isec] = rmax;
265  else
266  pgonRout[isec] = pgonRout[isec] * cosAlpha_ - tol1;
267  }
268 
269  dd4hep::Solid solid = dd4hep::Polyhedra(sectors_, -alpha_, 2._pi, pgonZ, pgonRin, pgonRout);
270  ns.addSolidNS(ns.prepend(name), solid);
271  glog = dd4hep::Volume(solid.name(), solid, matter);
272  ns.addVolumeNS(glog);
273 #ifdef EDM_ML_DEBUG
274  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: " << solid.name() << " polyhedra of " << sectors_
275  << " sectors covering " << convertRadToDeg(-alpha_) << ":"
276  << convertRadToDeg(-alpha_ + 2._pi) << " with " << pgonZ.size() << " sections";
277  for (unsigned int k = 0; k < pgonZ.size(); ++k)
278  edm::LogVerbatim("HGCalGeom") << "[" << k << "] z " << cms::convert2mm(pgonZ[k]) << " R "
279  << cms::convert2mm(pgonRin[k]) << ":" << cms::convert2mm(pgonRout[k]);
280 #endif
281  } else {
282  double rins =
283  (sensitiveMode_ < 1) ? rinB : HGCalGeomTools::radius(zz + hthick, zFrontB_, rMinFront_, slopeB_);
284  double routs =
285  (sensitiveMode_ < 1) ? routF : HGCalGeomTools::radius(zz - hthick, zFrontT_, rMaxFront_, slopeT_);
286  dd4hep::Solid solid = dd4hep::Tube(rins, routs, hthick, 0.0, 2._pi);
287  ns.addSolidNS(ns.prepend(name), solid);
288  glog = dd4hep::Volume(solid.name(), solid, matter);
289  ns.addVolumeNS(glog);
290 
291 #ifdef EDM_ML_DEBUG
292  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: " << solid.name() << " Tubs made of " << matter.name()
293  << " of dimensions " << cms::convert2mm(rinB) << ":" << cms::convert2mm(rins)
294  << ", " << cms::convert2mm(routF) << ":" << cms::convert2mm(routs) << ", "
295  << cms::convert2mm(hthick) << ", 0.0, 360.0 and positioned in: " << glog.name()
296  << " number " << copy;
297 #endif
298  positionMix(ctxt, e, glog, name, copy, thick_[ii], matter);
299  }
300 
301  dd4hep::Position r1(0, 0, zz);
302  mother.placeVolume(glog, copy, r1);
303  ++copyNumber_[ii];
304 #ifdef EDM_ML_DEBUG
305  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: " << glog.name() << " number " << copy << " positioned in "
306  << mother.name() << " at (0,0," << cms::convert2mm(zz) << ") with no rotation";
307 #endif
308  zz += hthick;
309  } // End of loop over layers in a block
310  zi = zo;
311  laymin = laymax;
312  if (std::abs(thickTot - layerThick_[i]) > tol2_) {
313  if (thickTot > layerThick_[i]) {
314  edm::LogError("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
315  << " is smaller than " << cms::convert2mm(thickTot)
316  << ": thickness of all its components **** ERROR ****";
317  } else {
318  edm::LogWarning("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
319  << " does not match with " << cms::convert2mm(thickTot) << " of the components";
320  }
321  }
322  } // End of loop over blocks
323 
324 #ifdef EDM_ML_DEBUG
325  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: " << copies_.size() << " different wafer copy numbers";
326  int k(0);
327  for (std::unordered_set<int>::const_iterator itr = copies_.begin(); itr != copies_.end(); ++itr, ++k) {
328  edm::LogVerbatim("HGCalGeom") << "Copy [" << k << "] : " << (*itr);
329  }
330  copies_.clear();
331  edm::LogVerbatim("HGCalGeom") << "<<== End of DDHGCalMixLayer construction...";
332 #endif
333  }
334 
335  void positionMix(cms::DDParsingContext& ctxt,
336  xml_h e,
337  const dd4hep::Volume& glog,
338  const std::string& nameM,
339  int copyM,
340  double thick,
341  const dd4hep::Material& matter) {
342  cms::DDNamespace ns(ctxt, e, true);
343 
344  // Make the top part first
345  for (unsigned int ly = 0; ly < layerTypeTop_.size(); ++ly) {
346  int ii = layerTypeTop_[ly];
347  copyNumberTop_[ii] = copyM;
348  }
349  double hthick = 0.5 * thick;
350  double dphi = (2._pi) / phiBinsScint_;
351  double thickTot(0), zpos(-hthick);
352  for (unsigned int ly = 0; ly < layerTypeTop_.size(); ++ly) {
353  int ii = layerTypeTop_[ly];
354  int copy = copyNumberTop_[ii];
355  int layer = copy - firstLayer_;
356  double hthickl = 0.5 * layerThickTop_[ii];
357  thickTot += layerThickTop_[ii];
358  zpos += hthickl;
359  dd4hep::Material matter1 = ns.material(materialTop_[ii]);
360  unsigned int k = 0;
361  int firstTile = tileLayerStart_[layer];
362  int lastTile = ((layer + 1 < static_cast<int>(tileLayerStart_.size())) ? tileLayerStart_[layer + 1]
363  : static_cast<int>(tileIndex_.size()));
364 #ifdef EDM_ML_DEBUG
365  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: Layer " << ly << ":" << ii << " Copy " << copy << " Tiles "
366  << firstTile << ":" << lastTile;
367 #endif
368  for (int ti = firstTile; ti < lastTile; ++ti) {
369  double r1 = tileRMin_[std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
370  double r2 = tileRMax_[std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
371  int fimin = std::get<1>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
372  int fimax = std::get<2>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
373  double phi1 = dphi * (fimin - 1);
374  double phi2 = dphi * (fimax - fimin + 1);
375 #ifdef EDM_ML_DEBUG
376  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: Layer " << copy << " iR "
377  << std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ly])) << ":"
378  << std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ly])) << " R "
379  << cms::convert2mm(r1) << ":" << cms::convert2mm(r2) << " Thick "
380  << cms::convert2mm((2.0 * hthickl)) << " phi " << fimin << ":" << fimax << ":"
381  << convertRadToDeg(phi1) << ":" << convertRadToDeg(phi2);
382 #endif
383  std::string name = namesTop_[ii] + "L" + std::to_string(copy) + "F" + std::to_string(k);
384  ++k;
385  dd4hep::Solid solid = dd4hep::Tube(r1, r2, hthickl, phi1, phi2);
386  ns.addSolidNS(ns.prepend(name), solid);
387  dd4hep::Volume glog1 = dd4hep::Volume(solid.name(), solid, matter1);
388  ns.addVolumeNS(glog1);
389 #ifdef EDM_ML_DEBUG
390  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: " << glog1.name() << " Tubs made of " << materialTop_[ii]
391  << " of dimensions " << cms::convert2mm(r1) << ", " << cms::convert2mm(r2) << ", "
392  << cms::convert2mm(hthickl) << ", " << convertRadToDeg(phi1) << ", "
393  << convertRadToDeg(phi2);
394 #endif
395  dd4hep::Position tran(0, 0, zpos);
396  glog.placeVolume(glog1, copy, tran);
397 #ifdef EDM_ML_DEBUG
398  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: Position " << glog1.name() << " number " << copy << " in "
399  << glog.name() << " at (0, 0, " << cms::convert2mm(zpos) << ") with no rotation";
400 #endif
401  }
402  ++copyNumberTop_[ii];
403  zpos += hthickl;
404  }
405  if (std::abs(thickTot - thick) > tol2_) {
406  if (thickTot > thick) {
407  edm::LogError("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(thick) << " is smaller than "
408  << cms::convert2mm(thickTot)
409  << ": thickness of all its components in the top part **** ERROR ****";
410  } else {
411  edm::LogWarning("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(thick)
412  << " does not match with " << cms::convert2mm(thickTot)
413  << " of the components in top part";
414  }
415  }
416 
417  // Make the bottom part next
418  int layer = (copyM - firstLayer_);
419  static const double sqrt3 = std::sqrt(3.0);
420  int layercenter = (layerTypes_[layer] == HGCalTypes::CornerCenteredLambda)
421  ? 1
422  : ((layerTypes_[layer] == HGCalTypes::CornerCenteredY) ? 2 : 0);
423  int layerType = (layerTypes_[layer] == HGCalTypes::WaferCenteredBack) ? 1 : 0;
424  int firstWafer = waferLayerStart_[layer];
425  int lastWafer = ((layer + 1 < static_cast<int>(waferLayerStart_.size())) ? waferLayerStart_[layer + 1]
426  : static_cast<int>(waferIndex_.size()));
427  double r = 0.5 * (waferSize_ + waferSepar_);
428  double R = 2.0 * r / sqrt3;
429  double dy = 0.75 * R;
430  const auto& xyoff = geomTools_.shiftXY(layercenter, (waferSize_ + waferSepar_));
431 #ifdef EDM_ML_DEBUG
432  int ium(0), ivm(0), kount(0);
433  std::vector<int> ntype(3, 0);
434  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: " << glog.name() << " r " << cms::convert2mm(r) << " R "
435  << cms::convert2mm(R) << " dy " << cms::convert2mm(dy) << " Shift "
436  << cms::convert2mm(xyoff.first) << ":" << cms::convert2mm(xyoff.second)
437  << " WaferSize " << cms::convert2mm((waferSize_ + waferSepar_)) << " index "
438  << firstWafer << ":" << (lastWafer - 1);
439 #endif
440  for (int k = firstWafer; k < lastWafer; ++k) {
441  int u = HGCalWaferIndex::waferU(waferIndex_[k]);
442  int v = HGCalWaferIndex::waferV(waferIndex_[k]);
443 #ifdef EDM_ML_DEBUG
444  int iu = std::abs(u);
445  int iv = std::abs(v);
446 #endif
447  int nr = 2 * v;
448  int nc = -2 * u + v;
449  double xpos = xyoff.first + nc * r;
450  double ypos = xyoff.second + nr * dy;
451  int type = HGCalProperty::waferThick(waferProperty_[k]);
452  int part = HGCalProperty::waferPartial(waferProperty_[k]);
453  int orien = HGCalProperty::waferOrient(waferProperty_[k]);
454  std::string wafer;
455  int i(999);
456  if (part == HGCalTypes::WaferFull) {
457  i = layerType * waferTypes_ + type;
458  wafer = waferFull_[i];
459  } else {
460  i = (part - 1) * waferTypes_ * facingTypes_ * orientationTypes_ + layerType * waferTypes_ * orientationTypes_ +
461  type * orientationTypes_ + orien;
462 #ifdef EDM_ML_DEBUG
463  edm::LogVerbatim("HGCalGeom") << " layertype:type:part:orien:ind " << layerType << ":" << type << ":" << part
464  << ":" << orien << ":" << i << ":" << waferPart_.size();
465 #endif
466  wafer = waferPart_[i];
467  }
468  int copy = HGCalTypes::packTypeUV(type, u, v);
469 #ifdef EDM_ML_DEBUG
470  edm::LogVerbatim("HGCalGeom") << " DDHGCalMixLayer: Layer " << HGCalWaferIndex::waferLayer(waferIndex_[k])
471  << " Wafer " << wafer << " number " << copy << " type :part:orien:ind " << type
472  << ":" << part << ":" << orien << ":" << i << " layer:u:v " << (layer + firstLayer_)
473  << ":" << u << ":" << v;
474  if (iu > ium)
475  ium = iu;
476  if (iv > ivm)
477  ivm = iv;
478  kount++;
479  if (copies_.count(copy) == 0)
480  copies_.insert(copy);
481 #endif
482  dd4hep::Position tran(xpos, ypos, 0.0);
483  glog.placeVolume(ns.volume(wafer), copy, tran);
484 #ifdef EDM_ML_DEBUG
485  ++ntype[type];
486  edm::LogVerbatim("HGCalGeom") << " DDHGCalMixLayer: " << wafer << " number " << copy << " type " << layerType
487  << ":" << type << " positioned in " << glog.name() << " at ("
488  << cms::convert2mm(xpos) << "," << cms::convert2mm(ypos) << ",0) with no rotation";
489 #endif
490  }
491 
492 #ifdef EDM_ML_DEBUG
493  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixLayer: Maximum # of u " << ium << " # of v " << ivm << " and " << kount
494  << " wafers (" << ntype[0] << ":" << ntype[1] << ":" << ntype[2] << ") for "
495  << glog.name();
496 #endif
497  }
498 
499  //Required data members to cache the values from XML file
500  HGCalGeomTools geomTools_;
501  static constexpr double tol2_ = 0.00001 * dd4hep::mm;
502 
503  int waferTypes_; // Number of wafer types
504  int facingTypes_; // Types of facings of modules toward IP
505  int partialTypes_; // Number of partial wafer types
506  int orientationTypes_; // Number of partial wafer orienations
507  int phiBinsScint_; // Maximum number of cells along phi
508  int firstLayer_; // Copy # of the first sensitive layer
509  int absorbMode_; // Absorber mode
510  int sensitiveMode_; // Sensitive mode
511  double zMinBlock_; // Starting z-value of the block
512  double waferSize_; // Width of the wafer
513  double waferSepar_; // Sensor separation
514  int sectors_; // Sectors
515  std::vector<double> slopeB_; // Slope at the lower R
516  std::vector<double> zFrontB_; // Starting Z values for the slopes
517  std::vector<double> rMinFront_; // Corresponding rMin's
518  std::vector<double> slopeT_; // Slopes at the larger R
519  std::vector<double> zFrontT_; // Starting Z values for the slopes
520  std::vector<double> rMaxFront_; // Corresponding rMax's
521  std::vector<std::string> waferFull_; // Names of full wafer modules
522  std::vector<std::string> waferPart_; // Names of partial wafer modules
523  std::vector<std::string> materials_; // Materials
524  std::vector<std::string> names_; // Names
525  std::vector<double> thick_; // Thickness of the material
526  std::vector<int> copyNumber_; // Initial copy numbers
527  std::vector<int> layers_; // Number of layers in a section
528  std::vector<double> layerThick_; // Thickness of each section
529  std::vector<int> layerType_; // Type of the layer
530  std::vector<int> layerSense_; // Content of a layer (sensitive?)
531  std::vector<std::string> materialTop_; // Materials of top layers
532  std::vector<std::string> namesTop_; // Names of top layers
533  std::vector<double> layerThickTop_; // Thickness of the top sections
534  std::vector<int> layerTypeTop_; // Type of the Top layer
535  std::vector<int> copyNumberTop_; // Initial copy numbers (top section)
536  std::vector<int> layerTypes_; // Layer type of silicon layers
537  std::vector<int> waferIndex_; // Wafer index for the types
538  std::vector<int> waferProperty_; // Wafer property
539  std::vector<int> waferLayerStart_; // Start index of wafers in each layer
540  std::vector<double> tileRMin_; // Minimum radius of each ring
541  std::vector<double> tileRMax_; // Maximum radius of each ring
542  std::vector<int> tileIndex_; // Index of tile (layer/start|end ring)
543  std::vector<int> tilePhis_; // Tile phi range for each index
544  std::vector<int> tileLayerStart_; // Start index of tiles in each layer
545  std::unordered_set<int> copies_; // List of copy #'s
546  double alpha_, cosAlpha_;
547 };
548 
549 static long algorithm(dd4hep::Detector& /* description */, cms::DDParsingContext& ctxt, xml_h e) {
550  HGCalMixLayer healgo(ctxt, e);
551  return cms::s_executed;
552 }
553 
554 DECLARE_DDCMS_DETELEMENT(DDCMS_hgcal_DDHGCalMixLayer, algorithm)
writedatasetfile.args
args
Definition: writedatasetfile.py:18
testProducerWithPsetDescEmpty_cfi.i2
i2
Definition: testProducerWithPsetDescEmpty_cfi.py:46
HGCalMixLayer::copies_
std::unordered_set< int > copies_
Definition: DDHGCalMixLayer.cc:545
HGCalMixLayer::waferLayerStart_
std::vector< int > waferLayerStart_
Definition: DDHGCalMixLayer.cc:539
mps_fire.i
i
Definition: mps_fire.py:428
geometryCSVtoXML.zz
zz
Definition: geometryCSVtoXML.py:19
g4SimHits_cfi.Material
Material
Definition: g4SimHits_cfi.py:594
HGCalMixLayer::tileRMin_
std::vector< double > tileRMin_
Definition: DDHGCalMixLayer.cc:540
MessageLogger.h
HGCalMixLayer::layerTypeTop_
std::vector< int > layerTypeTop_
Definition: DDHGCalMixLayer.cc:534
HGCalMixLayer
Definition: DDHGCalMixLayer.cc:28
HGCalMixLayer::materials_
std::vector< std::string > materials_
Definition: DDHGCalMixLayer.cc:523
HGCalProperty.h
DECLARE_DDCMS_DETELEMENT
#define DECLARE_DDCMS_DETELEMENT(name, func)
Definition: DDPlugins.h:25
HGCalMixLayer::phiBinsScint_
int phiBinsScint_
Definition: DDHGCalMixLayer.cc:507
cms::DDNamespace::material
dd4hep::Material material(const std::string &name) const
Definition: DDNamespace.cc:154
filterCSVwithJSON.copy
copy
Definition: filterCSVwithJSON.py:36
HGCalMixLayer::layerType_
std::vector< int > layerType_
Definition: DDHGCalMixLayer.cc:529
HGCalMixLayer::waferIndex_
std::vector< int > waferIndex_
Definition: DDHGCalMixLayer.cc:537
min
T min(T a, T b)
Definition: MathUtil.h:58
cms::DDParsingContext
Definition: DDParsingContext.h:14
angle_units::operators::convertRadToDeg
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
cms::DDNamespace
Definition: DDNamespace.h:22
gpuVertexFinder::iv
int32_t *__restrict__ iv
Definition: gpuClusterTracksDBSCAN.h:42
HGCalMixLayer::tilePhis_
std::vector< int > tilePhis_
Definition: DDHGCalMixLayer.cc:543
testProducerWithPsetDescEmpty_cfi.i1
i1
Definition: testProducerWithPsetDescEmpty_cfi.py:45
HGCalWaferIndex::waferU
int32_t waferU(const int32_t index)
Definition: HGCalWaferIndex.cc:27
HGCalProperty::waferOrient
int32_t waferOrient(const int32_t property)
Definition: HGCalProperty.cc:17
HGCalMixLayer::namesTop_
std::vector< std::string > namesTop_
Definition: DDHGCalMixLayer.cc:532
findQualityFiles.v
v
Definition: findQualityFiles.py:179
algorithm
static long algorithm(dd4hep::Detector &, cms::DDParsingContext &ctxt, xml_h e)
Definition: DDHGCalMixLayer.cc:549
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
HGCalMixLayer::thick_
std::vector< double > thick_
Definition: DDHGCalMixLayer.cc:525
HGCalMixLayer::tileLayerStart_
std::vector< int > tileLayerStart_
Definition: DDHGCalMixLayer.cc:544
HGCalMixLayer::geomTools_
HGCalGeomTools geomTools_
Definition: DDHGCalMixLayer.cc:500
HGCalMixLayer::partialTypes_
int partialTypes_
Definition: DDHGCalMixLayer.cc:505
HGCalMixLayer::facingTypes_
int facingTypes_
Definition: DDHGCalMixLayer.cc:504
HGCalTileIndex.h
part
part
Definition: HCALResponse.h:20
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
HGCalWaferIndex.h
HGCalMixLayer::layers_
std::vector< int > layers_
Definition: DDHGCalMixLayer.cc:527
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
HGCalMixLayer::HGCalMixLayer
HGCalMixLayer()
Definition: DDHGCalMixLayer.cc:29
cms::DDNamespace::prepend
std::string prepend(const std::string &) const
Definition: DDNamespace.cc:98
HGCalWaferType.h
HGCalTypes::WaferFull
Definition: HGCalTypes.h:58
HGCalWaferIndex::waferV
int32_t waferV(const int32_t index)
Definition: HGCalWaferIndex.cc:32
HGCalMixLayer::orientationTypes_
int orientationTypes_
Definition: DDHGCalMixLayer.cc:506
PixelTestBeamValidation_cfi.Position
Position
Definition: PixelTestBeamValidation_cfi.py:83
dqmdumpme.k
k
Definition: dqmdumpme.py:60
HGCalMixLayer::layerTypes_
std::vector< int > layerTypes_
Definition: DDHGCalMixLayer.cc:536
HGCalMixLayer::layerSense_
std::vector< int > layerSense_
Definition: DDHGCalMixLayer.cc:530
HGCalTypes::packTypeUV
static int32_t packTypeUV(int type, int u, int v)
Definition: HGCalTypes.cc:3
HGCalMixLayer::waferPart_
std::vector< std::string > waferPart_
Definition: DDHGCalMixLayer.cc:522
phase1PixelTopology::layer
constexpr std::array< uint8_t, layerIndexSize > layer
Definition: phase1PixelTopology.h:99
HGCalMixLayer::waferTypes_
int waferTypes_
Definition: DDHGCalMixLayer.cc:503
cms::Volume
dd4hep::Volume Volume
Definition: DDFilteredView.h:47
HGCalMixLayer::positionMix
void positionMix(cms::DDParsingContext &ctxt, xml_h e, const dd4hep::Volume &glog, const std::string &nameM, int copyM, double thick, const dd4hep::Material &matter)
Definition: DDHGCalMixLayer.cc:335
EgHLTOffHistBins_cfi.nr
nr
Definition: EgHLTOffHistBins_cfi.py:4
HGCalTileIndex::tileUnpack
std::tuple< int32_t, int32_t, int32_t > tileUnpack(int32_t index)
Definition: HGCalTileIndex.cc:44
HGCalTypes::CornerCenteredLambda
Definition: HGCalTypes.h:86
HGCalTypes::CornerCenteredY
Definition: HGCalTypes.h:85
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
angle_units::operators
Definition: angle_units.h:11
gainCalibHelper::gainCalibPI::type
type
Definition: SiPixelGainCalibHelper.h:40
DDPlugins.h
HGCalMixLayer::sensitiveMode_
int sensitiveMode_
Definition: DDHGCalMixLayer.cc:510
HGCalMixLayer::rMaxFront_
std::vector< double > rMaxFront_
Definition: DDHGCalMixLayer.cc:520
diffTwoXMLs.r2
r2
Definition: diffTwoXMLs.py:73
HGCalMixLayer::layerThick_
std::vector< double > layerThick_
Definition: DDHGCalMixLayer.cc:528
HGCalMixLayer::layerThickTop_
std::vector< double > layerThickTop_
Definition: DDHGCalMixLayer.cc:533
HGCalProperty::waferThick
int32_t waferThick(const int32_t property)
Definition: HGCalProperty.cc:9
cms::convert2mm
constexpr NumType convert2mm(NumType length)
Definition: DDutils.h:7
HGCalMixLayer::slopeT_
std::vector< double > slopeT_
Definition: DDHGCalMixLayer.cc:518
HGCalGeomTools
Definition: HGCalGeomTools.h:8
cms::DDNamespace::addSolidNS
dd4hep::Solid addSolidNS(const std::string &name, dd4hep::Solid solid) const
Definition: DDNamespace.cc:280
cms::DDAlgoArguments
Definition: DDAlgoArguments.h:28
PVValHelper::dy
Definition: PVValidationHelpers.h:50
HGCalMixLayer::materialTop_
std::vector< std::string > materialTop_
Definition: DDHGCalMixLayer.cc:531
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
cms::DDNamespace::addVolumeNS
dd4hep::Volume addVolumeNS(dd4hep::Volume vol) const
Definition: DDNamespace.cc:190
cms::s_executed
static constexpr long s_executed
Definition: DDAlgoArguments.h:16
HGCalProperty::waferPartial
int32_t waferPartial(const int32_t property)
Definition: HGCalProperty.cc:13
alignCSCRings.r
r
Definition: alignCSCRings.py:93
HGCalMixLayer::waferFull_
std::vector< std::string > waferFull_
Definition: DDHGCalMixLayer.cc:521
align::Detector
Definition: StructureType.h:92
HGCalMixLayer::zFrontB_
std::vector< double > zFrontB_
Definition: DDHGCalMixLayer.cc:516
HGCalGeomTools.h
HGCalWaferIndex::waferLayer
int32_t waferLayer(const int32_t index)
Definition: HGCalWaferIndex.cc:23
HGCalGeomTools::radius
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)
Definition: HGCalGeomTools.cc:12
HGCalMixLayer::slopeB_
std::vector< double > slopeB_
Definition: DDHGCalMixLayer.cc:515
HGCalMixLayer::absorbMode_
int absorbMode_
Definition: DDHGCalMixLayer.cc:509
HGCalMixLayer::HGCalMixLayer
HGCalMixLayer(cms::DDParsingContext &ctxt, xml_h e)
Definition: DDHGCalMixLayer.cc:30
HGCalMixLayer::sectors_
int sectors_
Definition: DDHGCalMixLayer.cc:514
diffTwoXMLs.r1
r1
Definition: diffTwoXMLs.py:53
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
HGCalMixLayer::zFrontT_
std::vector< double > zFrontT_
Definition: DDHGCalMixLayer.cc:519
HGCalGeometryMode::Polyhedra
Definition: HGCalGeometryMode.h:37
HGCalMixLayer::tileRMax_
std::vector< double > tileRMax_
Definition: DDHGCalMixLayer.cc:541
Exception
Definition: hltDiff.cc:245
HGCalMixLayer::copyNumberTop_
std::vector< int > copyNumberTop_
Definition: DDHGCalMixLayer.cc:535
HGCalMixLayer::names_
std::vector< std::string > names_
Definition: DDHGCalMixLayer.cc:524
HGCalMixLayer::waferProperty_
std::vector< int > waferProperty_
Definition: DDHGCalMixLayer.cc:538
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
HGCalTypes::WaferCenteredBack
Definition: HGCalTypes.h:84
HGCalMixLayer::copyNumber_
std::vector< int > copyNumber_
Definition: DDHGCalMixLayer.cc:526
HGCalMixLayer::waferSepar_
double waferSepar_
Definition: DDHGCalMixLayer.cc:513
HGCalMixLayer::tileIndex_
std::vector< int > tileIndex_
Definition: DDHGCalMixLayer.cc:542
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
HGCalMixLayer::waferSize_
double waferSize_
Definition: DDHGCalMixLayer.cc:512
angle_units.h
HGCalTypes.h
HGCalMixLayer::rMinFront_
std::vector< double > rMinFront_
Definition: DDHGCalMixLayer.cc:517
HGCalMixLayer::cosAlpha_
double cosAlpha_
Definition: DDHGCalMixLayer.cc:546
edm::Log
Definition: MessageLogger.h:70
cms::DDNamespace::name
std::string_view name() const
Definition: DDNamespace.h:79
dttmaxenums::R
Definition: DTTMax.h:29
HGCalParameters.h
cuy.ii
ii
Definition: cuy.py:589
cms::DDNamespace::volume
dd4hep::Volume volume(const std::string &name, bool exc=true) const
Definition: DDNamespace.cc:264
DDutils.h
HGCalMixLayer::firstLayer_
int firstLayer_
Definition: DDHGCalMixLayer.cc:508
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
HGCalMixLayer::zMinBlock_
double zMinBlock_
Definition: DDHGCalMixLayer.cc:511
Generated for CMSSW Reference Manual by doxygen 1.8.16