CMS 3D CMS Logo
List of all members | Public Member Functions | Public Attributes | Static Public Attributes
HGCalMixRotatedCassette Struct Reference

Public Member Functions

 HGCalMixRotatedCassette ()
 
 HGCalMixRotatedCassette (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, int absType)
 

Public Attributes

int absorbMode_
 
double alpha_
 
HGCalCassette cassette_
 
int cassettes_
 
std::vector< double > cassetteShift_
 
std::unordered_set< int > copies_
 
std::vector< int > copyNumber_
 
int copyNumberCoverTop_
 
std::vector< int > copyNumberTop_
 
double cosAlpha_
 
int coverTypeTop_
 
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 partialTypes_
 
std::vector< std::string > passiveFull_
 
std::vector< std::string > passivePart_
 
int passiveTypes_
 
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 DDHGCalMixRotatedCassette.cc.

Constructor & Destructor Documentation

◆ HGCalMixRotatedCassette() [1/2]

HGCalMixRotatedCassette::HGCalMixRotatedCassette ( )
inline

Definition at line 31 of file DDHGCalMixRotatedCassette.cc.

References Exception.

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

◆ HGCalMixRotatedCassette() [2/2]

HGCalMixRotatedCassette::HGCalMixRotatedCassette ( cms::DDParsingContext &  ctxt,
xml_h  e 
)
inline

Definition at line 32 of file DDHGCalMixRotatedCassette.cc.

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

Member Function Documentation

◆ positionMix()

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

Definition at line 384 of file DDHGCalMixRotatedCassette.cc.

References funct::abs(), cms::DDNamespace::addSolidNS(), cms::DDNamespace::addVolumeNS(), angle(), Matriplex::atan2(), HGCalCassette::cassetteType(), HGCalCell::cellPlacementIndex(), cms::convert2mm(), angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, PVValHelper::dy, MillePedeFileConverter_cfg::e, mps_fire::i, cuy::ii, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, isotrackApplyRegressor::k, HGCalTypes::layerFrontBack(), cms::DDNamespace::material(), g4SimHits_cfi::Material, mergeVDriftHistosByStation::name, EgHLTOffHistBins_cfi::nr, HGCalTypes::packTypeUV(), PixelTestBeamValidation_cfi::Position, cms::DDNamespace::prepend(), diffTwoXMLs::r2, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, HGCalTileIndex::tileUnpack(), 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().

391  {
392  cms::DDNamespace ns(ctxt, e, true);
393 
394  // Make the top part first
395  for (unsigned int ly = 0; ly < layerTypeTop_.size(); ++ly) {
396  int ii = layerTypeTop_[ly];
397  copyNumberTop_[ii] = copyM;
398  }
399  double hthick = 0.5 * thick;
400  double dphi = (2._pi) / phiBinsScint_;
401  double thickTot(0), zpos(-hthick);
402  if (absType < 0) {
403  for (unsigned int ly = 0; ly < layerTypeTop_.size(); ++ly) {
404  int ii = layerTypeTop_[ly];
405  int copy = copyNumberTop_[ii];
406  int layer = copy - firstLayer_;
407  double hthickl = 0.5 * layerThickTop_[ii];
408  thickTot += layerThickTop_[ii];
409  zpos += hthickl;
410  dd4hep::Material matter1 = ns.material(materialTop_[ii]);
411  unsigned int k = 0;
412  int firstTile = tileLayerStart_[layer];
413  int lastTile = ((layer + 1 < static_cast<int>(tileLayerStart_.size())) ? tileLayerStart_[layer + 1]
414  : static_cast<int>(tileIndex_.size()));
415 #ifdef EDM_ML_DEBUG
416  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Layer " << ly << ":" << ii << " Copy " << copy
417  << " Tiles " << firstTile << ":" << lastTile;
418 #endif
419  for (int ti = firstTile; ti < lastTile; ++ti) {
420  double r1 = tileRMin_[std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
421  double r2 = tileRMax_[std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
422  int cassette = std::get<0>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
423  int fimin = std::get<1>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
424  int fimax = std::get<2>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
425  double phi1 = dphi * (fimin - 1);
426  double phi2 = (forFireworks_ == 1) ? (dphi * (fimax - fimin + 1)) : (dphi * fimax);
427  auto cshift = cassette_.getShift(layer + 1, 1, cassette);
428 #ifdef EDM_ML_DEBUG
429  int cassette0 = HGCalCassette::cassetteType(2, 1, cassette); //
430  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Layer " << copy << ":" << (layer + 1) << " iR "
431  << std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) << ":"
432  << std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) << " R "
433  << cms::convert2mm(r1) << ":" << cms::convert2mm(r2) << " Thick "
434  << cms::convert2mm((2.0 * hthickl)) << " phi " << fimin << ":" << fimax << ":"
435  << convertRadToDeg(phi1) << ":" << convertRadToDeg(phi2) << " cassette "
436  << cassette << ":" << cassette0 << " Shift " << cms::convert2mm(cshift.first)
437  << ":" << cms::convert2mm(cshift.second);
438 #endif
439  std::string name = namesTop_[ii] + "L" + std::to_string(copy) + "F" + std::to_string(k);
440  ++k;
441  dd4hep::Solid solid = dd4hep::Tube(r1, r2, hthickl, phi1, phi2);
442  ns.addSolidNS(ns.prepend(name), solid);
443  dd4hep::Volume glog1 = dd4hep::Volume(solid.name(), solid, matter1);
444  ns.addVolumeNS(glog1);
445 #ifdef EDM_ML_DEBUG
446  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << glog1.name() << " Tubs made of "
447  << materialTop_[ii] << " of dimensions " << cms::convert2mm(r1) << ", "
448  << cms::convert2mm(r2) << ", " << cms::convert2mm(hthickl) << ", "
449  << convertRadToDeg(phi1) << ", " << convertRadToDeg(phi2);
450 #endif
451  dd4hep::Position tran(-cshift.first, cshift.second, zpos);
452  glog.placeVolume(glog1, copy, tran);
453 #ifdef EDM_ML_DEBUG
454  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Position " << glog1.name() << " number " << copy
455  << " in " << glog.name() << " at (" << cms::convert2mm(cshift.first) << ", "
456  << cms::convert2mm(cshift.second) << ", " << cms::convert2mm(zpos)
457  << ") with no rotation";
458 #endif
459  }
460  ++copyNumberTop_[ii];
461  zpos += hthickl;
462  }
463  if (std::abs(thickTot - thick) > tol2_) {
464  if (thickTot > thick) {
465  edm::LogError("HGCalGeom") << "DDHGCalMixRotatedCassette: Thickness of the partition "
466  << cms::convert2mm(thick) << " is smaller than " << cms::convert2mm(thickTot)
467  << ": thickness of all its components in the top part **** ERROR ****";
468  } else {
469  edm::LogWarning("HGCalGeom") << "DDHGCalMixRotatedCassette: Thickness of the partition "
470  << cms::convert2mm(thick) << " does not match with " << cms::convert2mm(thickTot)
471  << " of the components in top part";
472  }
473  }
474  } else {
475  int ii = coverTypeTop_;
476  int copy = copyNumberCoverTop_;
477  int layer = copy - firstLayer_;
478  double hthickl = 0.5 * layerThickTop_[ii];
479  zpos += hthickl;
480  dd4hep::Material matter1 = ns.material(materialTop_[ii]);
481  unsigned int k = 0;
482  int firstTile = tileLayerStart_[layer];
483  int lastTile = ((layer + 1 < static_cast<int>(tileLayerStart_.size())) ? tileLayerStart_[layer + 1]
484  : static_cast<int>(tileIndex_.size()));
485 #ifdef EDM_ML_DEBUG
486  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Layer " << layer << ":" << ii << " Copy " << copy
487  << " Tiles " << firstTile << ":" << lastTile;
488 #endif
489  for (int ti = firstTile; ti < lastTile; ++ti) {
490  double r1 = tileRMin_[std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
491  double r2 = tileRMax_[std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
492  int cassette = std::get<0>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
493  int fimin = std::get<1>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
494  int fimax = std::get<2>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
495  double phi1 = dphi * (fimin - 1);
496  double phi2 = (forFireworks_ == 1) ? (dphi * (fimax - fimin + 1)) : (dphi * fimax);
497  auto cshift = cassette_.getShift(layer + 1, 1, cassette);
498 #ifdef EDM_ML_DEBUG
499  int cassette0 = HGCalCassette::cassetteType(2, 1, cassette); //
500  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Layer " << copy << ":" << (layer + 1) << " iR "
501  << std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) << ":"
502  << std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) << " R "
503  << cms::convert2mm(r1) << ":" << cms::convert2mm(r2) << " Thick "
504  << cms::convert2mm(2.0 * hthickl) << " phi " << fimin << ":" << fimax << ":"
505  << convertRadToDeg(phi1) << ":" << convertRadToDeg(phi2) << " cassette "
506  << cassette << ":" << cassette0 << " Shift " << cms::convert2mm(cshift.first)
507  << ":" << cms::convert2mm(cshift.second);
508 #endif
509  std::string name = namesTop_[ii] + "L" + std::to_string(copy) + "F" + std::to_string(k);
510  ++k;
511  dd4hep::Solid solid = dd4hep::Tube(r1, r2, hthickl, phi1, phi2);
512  ns.addSolidNS(ns.prepend(name), solid);
513  dd4hep::Volume glog1 = dd4hep::Volume(solid.name(), solid, matter1);
514  ns.addVolumeNS(glog1);
515 #ifdef EDM_ML_DEBUG
516  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << glog1.name() << " Tubs made of "
517  << matter1.name() << " of dimensions " << cms::convert2mm(r1) << ", "
518  << cms::convert2mm(r2) << ", " << cms::convert2mm(hthickl) << ", "
519  << convertRadToDeg(phi1) << ", " << convertRadToDeg(phi2);
520 #endif
521  dd4hep::Position tran(-cshift.first, cshift.second, zpos);
522  glog.placeVolume(glog1, copy, tran);
523 #ifdef EDM_ML_DEBUG
524  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Position " << glog1.name() << " number " << copy
525  << " in " << glog.name() << " at (" << cms::convert2mm(cshift.first) << ", "
526  << cms::convert2mm(cshift.second) << ", " << cms::convert2mm(zpos)
527  << ") with no rotation";
528 #endif
529  }
530  ++copyNumberCoverTop_;
531  }
532 
533  // Make the bottom part next
534  int layer = (copyM - firstLayer_);
535  static const double sqrt3 = std::sqrt(3.0);
536  int layercenter = layerOrient_[layer];
537  int layertype = HGCalTypes::layerFrontBack(layerOrient_[layer]);
538  int firstWafer = waferLayerStart_[layer];
539  int lastWafer = ((layer + 1 < static_cast<int>(waferLayerStart_.size())) ? waferLayerStart_[layer + 1]
540  : static_cast<int>(waferIndex_.size()));
541  double delx = 0.5 * (waferSize_ + waferSepar_);
542  double dely = 2.0 * delx / sqrt3;
543  double dy = 0.75 * dely;
544  const auto& xyoff = geomTools_.shiftXY(layercenter, (waferSize_ + waferSepar_));
545 #ifdef EDM_ML_DEBUG
546  int ium(0), ivm(0), kount(0);
547  std::vector<int> ntype(3, 0);
548  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << glog.name() << " r " << cms::convert2mm(delx)
549  << " R " << cms::convert2mm(dely) << " dy " << cms::convert2mm(dy) << " Shift "
550  << cms::convert2mm(xyoff.first) << ":" << cms::convert2mm(xyoff.second)
551  << " WaferSize " << cms::convert2mm((waferSize_ + waferSepar_)) << " index "
552  << firstWafer << ":" << (lastWafer - 1) << " Copy " << copyM << ":" << layer;
553 #endif
554  for (int k = firstWafer; k < lastWafer; ++k) {
555  int u = HGCalWaferIndex::waferU(waferIndex_[k]);
556  int v = HGCalWaferIndex::waferV(waferIndex_[k]);
557 #ifdef EDM_ML_DEBUG
558  int iu = std::abs(u);
559  int iv = std::abs(v);
560 #endif
561  int nr = 2 * v;
562  int nc = -2 * u + v;
563  int type = HGCalProperty::waferThick(waferProperty_[k]);
564  int part = HGCalProperty::waferPartial(waferProperty_[k]);
565  int orien = HGCalProperty::waferOrient(waferProperty_[k]);
566  int cassette = HGCalProperty::waferCassette(waferProperty_[k]);
567  int place = HGCalCell::cellPlacementIndex(1, layertype, orien);
568 #ifdef EDM_ML_DEBUG
569  edm::LogVerbatim("HGCalGeom")
570  << "DDHGCalMixRotatedCassette::index:Property:layertype:type:part:orien:cassette:place:offsets:ind " << k
571  << ":" << waferProperty_[k] << ":" << layertype << ":" << type << ":" << part << ":" << orien << ":"
572  << cassette << ":" << place;
573 #endif
574  auto cshift = cassette_.getShift(layer + 1, -1, cassette);
575  double xpos = xyoff.first - cshift.first + nc * delx;
576  double ypos = xyoff.second + cshift.second + nr * dy;
577 #ifdef EDM_ML_DEBUG
578  double xorig = xyoff.first + nc * delx;
579  double yorig = xyoff.second + nr * dy;
580  double angle = std::atan2(yorig, xorig);
581  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette::Wafer: layer " << layer + 1 << " cassette "
582  << cassette << " Shift " << cms::convert2mm(cshift.first) << ":"
583  << cms::convert2mm(cshift.second) << " Original " << cms::convert2mm(xorig) << ":"
584  << cms::convert2mm(yorig) << ":" << convertRadToDeg(angle) << " Final "
585  << cms::convert2mm(xpos) << ":" << cms::convert2mm(ypos);
586 #endif
587  std::string wafer;
588  int i(999);
589  if (absType < 0) {
590  if (part == HGCalTypes::WaferFull) {
591  i = type * facingTypes_ * orientationTypes_ + place - placeOffset_;
592 #ifdef EDM_ML_DEBUG
593  edm::LogVerbatim("HGCalGeom") << " FullWafer type:place:ind " << type << ":" << place << ":" << i << ":"
594  << waferFull_.size();
595 #endif
596  wafer = waferFull_[i];
597  } else {
598  int partoffset =
599  (part >= HGCalTypes::WaferHDTop) ? HGCalTypes::WaferPartHDOffset : HGCalTypes::WaferPartLDOffset;
600  i = (part - partoffset) * facingTypes_ * orientationTypes_ +
601  HGCalTypes::WaferTypeOffset[type] * facingTypes_ * orientationTypes_ + place - placeOffset_;
602 
603 #ifdef EDM_ML_DEBUG
604  edm::LogVerbatim("HGCalGeom") << " layertype:type:part:orien:cassette:place:offsets:ind " << layertype << ":"
605  << type << ":" << part << ":" << orien << ":" << cassette << ":" << place << ":"
606  << partoffset << ":" << HGCalTypes::WaferTypeOffset[type] << ":" << i << ":"
607  << waferPart_.size();
608 #endif
609  wafer = waferPart_[i];
610  }
611  } else {
612  type = absType;
613  if (part == HGCalTypes::WaferFull) {
614  i = absType - 1;
615  wafer = passiveFull_[i];
616 #ifdef EDM_ML_DEBUG
617  edm::LogVerbatim("HGCalGeom") << " layertype:abstype:part:orien:cassette:offsets:ind " << layertype << ":"
618  << absType << ":" << part << ":" << orien << ":" << cassette << ":"
619  << ":" << partialTypes_ << ":" << orientationTypes_ << " passive " << i << ":"
620  << wafer;
621 #endif
622  } else {
623  int partoffset = (part >= HGCalTypes::WaferHDTop)
624  ? HGCalTypes::WaferPartHDOffset
625  : (HGCalTypes::WaferPartLDOffset - HGCalTypes::WaferTypeOffset[1]);
626  i = (part - partoffset) * facingTypes_ * orientationTypes_ +
627  (absType - 1) * facingTypes_ * orientationTypes_ * partialTypes_ + place - placeOffset_;
628 #ifdef EDM_ML_DEBUG
629  edm::LogVerbatim("HGCalGeom") << " layertype:abstype:part:orien:cassette:3Types:offset:ind " << layertype
630  << ":" << absType << ":" << part << ":" << orien << ":" << cassette << ":"
631  << partialTypes_ << ":" << facingTypes_ << ":" << orientationTypes_ << ":"
632  << partoffset << ":" << i << ":" << passivePart_.size();
633 #endif
634  wafer = passivePart_[i];
635  }
636  }
637 
638  int copy = HGCalTypes::packTypeUV(type, u, v);
639 #ifdef EDM_ML_DEBUG
640  edm::LogVerbatim("HGCalGeom") << " DDHGCalMixRotatedCassette: Layer "
641  << HGCalWaferIndex::waferLayer(waferIndex_[k]) << " Wafer " << wafer << " number "
642  << copy << " type :part:orien:ind " << type << ":" << part << ":" << orien << ":"
643  << i << " layer:u:v " << (layer + firstLayer_) << ":" << u << ":" << v;
644  if (iu > ium)
645  ium = iu;
646  if (iv > ivm)
647  ivm = iv;
648  kount++;
649  if (copies_.count(copy) == 0)
650  copies_.insert(copy);
651 #endif
652  dd4hep::Position tran(xpos, ypos, 0.0);
653  glog.placeVolume(ns.volume(wafer), copy, tran);
654 #ifdef EDM_ML_DEBUG
655  ++ntype[type];
656  edm::LogVerbatim("HGCalGeom") << " DDHGCalMixRotatedCassette: " << wafer << " number " << copy << " type "
657  << layertype << ":" << type << " positioned in " << glog.name() << " at ("
658  << cms::convert2mm(xpos) << "," << cms::convert2mm(ypos) << ",0) with no rotation";
659 #endif
660  }
661 
662 #ifdef EDM_ML_DEBUG
663  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Maximum # of u " << ium << " # of v " << ivm << " and "
664  << kount << " wafers (" << ntype[0] << ":" << ntype[1] << ":" << ntype[2] << ") for "
665  << glog.name();
666 #endif
667  }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:134
HGCalCell::cellPlacementIndex
static int32_t cellPlacementIndex(int32_t iz, int32_t frontBack, int32_t orient)
Definition: HGCalCell.cc:239
EgHLTOffHistBins_cfi.nr
nr
Definition: EgHLTOffHistBins_cfi.py:4
HGCalTypes::WaferPartLDOffset
static constexpr int32_t WaferPartLDOffset
Definition: HGCalTypes.h:57
HGCalMixRotatedCassette::copies_
std::unordered_set< int > copies_
Definition: DDHGCalMixRotatedCassette.cc:726
mps_fire.i
i
Definition: mps_fire.py:429
HGCalTypes::WaferTypeOffset
static constexpr int32_t WaferTypeOffset[3]
Definition: HGCalTypes.h:61
HGCalMixRotatedCassette::namesTop_
std::vector< std::string > namesTop_
Definition: DDHGCalMixRotatedCassette.cc:710
HGCalWaferIndex::waferU
int32_t waferU(const int32_t index)
Definition: HGCalWaferIndex.cc:27
HGCalWaferIndex::waferLayer
int32_t waferLayer(const int32_t index)
Definition: HGCalWaferIndex.cc:23
angle_units::operators::convertRadToDeg
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
cms::convert2mm
constexpr NumType convert2mm(NumType length)
Definition: DDutils.h:7
HGCalCassette::getShift
std::pair< double, double > getShift(int layer, int zside, int cassette, bool scnt=false) const
Definition: HGCalCassette.cc:55
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv
int32_t *__restrict__ iv
Definition: fitVertices.h:39
HGCalMixRotatedCassette::waferPart_
std::vector< std::string > waferPart_
Definition: DDHGCalMixRotatedCassette.cc:698
PixelTestBeamValidation_cfi.Position
Position
Definition: PixelTestBeamValidation_cfi.py:83
HGCalMixRotatedCassette::layerThickTop_
std::vector< double > layerThickTop_
Definition: DDHGCalMixRotatedCassette.cc:711
g4SimHits_cfi.Material
Material
Definition: g4SimHits_cfi.py:651
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
HGCalProperty::waferOrient
int32_t waferOrient(const int32_t property)
Definition: HGCalProperty.cc:21
HGCalGeomTools::shiftXY
std::pair< double, double > shiftXY(int waferPosition, double waferSize) const
Definition: HGCalGeomTools.cc:234
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
to_string
static std::string to_string(const XMLCh *ch)
Definition: TotemDAQMappingESSourceXML.cc:353
HGCalMixRotatedCassette::materialTop_
std::vector< std::string > materialTop_
Definition: DDHGCalMixRotatedCassette.cc:709
HGCalProperty::waferCassette
int32_t waferCassette(const int32_t property)
Definition: HGCalProperty.cc:25
HGCalTypes::WaferFull
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:35
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:23
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
HGCalMixRotatedCassette::waferFull_
std::vector< std::string > waferFull_
Definition: DDHGCalMixRotatedCassette.cc:697
HGCalProperty::waferThick
int32_t waferThick(const int32_t property)
Definition: HGCalProperty.cc:13
cms::Volume
dd4hep::Volume Volume
Definition: DDFilteredView.h:47
HGCalMixRotatedCassette::passivePart_
std::vector< std::string > passivePart_
Definition: DDHGCalMixRotatedCassette.cc:700
cuy.ii
ii
Definition: cuy.py:589
HGCalTileIndex::tileUnpack
std::tuple< int32_t, int32_t, int32_t > tileUnpack(int32_t index)
Definition: HGCalTileIndex.cc:59
HGCalTypes::WaferPartHDOffset
static constexpr int32_t WaferPartHDOffset
Definition: HGCalTypes.h:58
nano_mu_digi_cff.layer
layer
Definition: nano_mu_digi_cff.py:28
part
part
Definition: HCALResponse.h:20
HGCalProperty::waferPartial
int32_t waferPartial(const int32_t property)
Definition: HGCalProperty.cc:17
HGCalWaferIndex::waferV
int32_t waferV(const int32_t index)
Definition: HGCalWaferIndex.cc:32
HGCalTypes::packTypeUV
static int32_t packTypeUV(int type, int u, int v)
Definition: HGCalTypes.cc:3
HGCalTypes::WaferHDTop
static constexpr int32_t WaferHDTop
Definition: HGCalTypes.h:51
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:128
HGCalCassette::cassetteType
static int cassetteType(int det, int zside, int cassette)
Definition: HGCalCassette.cc:73
HGCalMixRotatedCassette::passiveFull_
std::vector< std::string > passiveFull_
Definition: DDHGCalMixRotatedCassette.cc:699
Matriplex::atan2
MPlex< T, D1, D2, N > atan2(const MPlex< T, D1, D2, N > &y, const MPlex< T, D1, D2, N > &x)
Definition: Matriplex.h:648
angle
T angle(T x1, T y1, T z1, T x2, T y2, T z2)
Definition: angle.h:11
HGCalTypes::layerFrontBack
static constexpr int32_t layerFrontBack(int32_t layerOrient)
Definition: HGCalTypes.h:137

Member Data Documentation

◆ absorbMode_

int HGCalMixRotatedCassette::absorbMode_

Definition at line 684 of file DDHGCalMixRotatedCassette.cc.

◆ alpha_

double HGCalMixRotatedCassette::alpha_

Definition at line 727 of file DDHGCalMixRotatedCassette.cc.

◆ cassette_

HGCalCassette HGCalMixRotatedCassette::cassette_

Definition at line 671 of file DDHGCalMixRotatedCassette.cc.

◆ cassettes_

int HGCalMixRotatedCassette::cassettes_

Definition at line 690 of file DDHGCalMixRotatedCassette.cc.

◆ cassetteShift_

std::vector<double> HGCalMixRotatedCassette::cassetteShift_

Definition at line 720 of file DDHGCalMixRotatedCassette.cc.

◆ copies_

std::unordered_set<int> HGCalMixRotatedCassette::copies_

Definition at line 726 of file DDHGCalMixRotatedCassette.cc.

◆ copyNumber_

std::vector<int> HGCalMixRotatedCassette::copyNumber_

Definition at line 704 of file DDHGCalMixRotatedCassette.cc.

◆ copyNumberCoverTop_

int HGCalMixRotatedCassette::copyNumberCoverTop_

Definition at line 715 of file DDHGCalMixRotatedCassette.cc.

◆ copyNumberTop_

std::vector<int> HGCalMixRotatedCassette::copyNumberTop_

Definition at line 713 of file DDHGCalMixRotatedCassette.cc.

◆ cosAlpha_

double HGCalMixRotatedCassette::cosAlpha_

Definition at line 727 of file DDHGCalMixRotatedCassette.cc.

◆ coverTypeTop_

int HGCalMixRotatedCassette::coverTypeTop_

Definition at line 714 of file DDHGCalMixRotatedCassette.cc.

◆ facingTypes_

int HGCalMixRotatedCassette::facingTypes_

Definition at line 677 of file DDHGCalMixRotatedCassette.cc.

◆ firstLayer_

int HGCalMixRotatedCassette::firstLayer_

Definition at line 683 of file DDHGCalMixRotatedCassette.cc.

◆ forFireworks_

int HGCalMixRotatedCassette::forFireworks_

Definition at line 682 of file DDHGCalMixRotatedCassette.cc.

◆ geomTools_

HGCalGeomTools HGCalMixRotatedCassette::geomTools_

Definition at line 670 of file DDHGCalMixRotatedCassette.cc.

◆ layerOrient_

std::vector<int> HGCalMixRotatedCassette::layerOrient_

Definition at line 716 of file DDHGCalMixRotatedCassette.cc.

◆ layers_

std::vector<int> HGCalMixRotatedCassette::layers_

Definition at line 705 of file DDHGCalMixRotatedCassette.cc.

◆ layerSense_

std::vector<int> HGCalMixRotatedCassette::layerSense_

Definition at line 708 of file DDHGCalMixRotatedCassette.cc.

◆ layerThick_

std::vector<double> HGCalMixRotatedCassette::layerThick_

Definition at line 706 of file DDHGCalMixRotatedCassette.cc.

◆ layerThickTop_

std::vector<double> HGCalMixRotatedCassette::layerThickTop_

Definition at line 711 of file DDHGCalMixRotatedCassette.cc.

◆ layerType_

std::vector<int> HGCalMixRotatedCassette::layerType_

Definition at line 707 of file DDHGCalMixRotatedCassette.cc.

◆ layerTypeTop_

std::vector<int> HGCalMixRotatedCassette::layerTypeTop_

Definition at line 712 of file DDHGCalMixRotatedCassette.cc.

◆ materials_

std::vector<std::string> HGCalMixRotatedCassette::materials_

Definition at line 701 of file DDHGCalMixRotatedCassette.cc.

◆ materialTop_

std::vector<std::string> HGCalMixRotatedCassette::materialTop_

Definition at line 709 of file DDHGCalMixRotatedCassette.cc.

◆ names_

std::vector<std::string> HGCalMixRotatedCassette::names_

Definition at line 702 of file DDHGCalMixRotatedCassette.cc.

◆ namesTop_

std::vector<std::string> HGCalMixRotatedCassette::namesTop_

Definition at line 710 of file DDHGCalMixRotatedCassette.cc.

◆ orientationTypes_

int HGCalMixRotatedCassette::orientationTypes_

Definition at line 678 of file DDHGCalMixRotatedCassette.cc.

◆ partialTypes_

int HGCalMixRotatedCassette::partialTypes_

Definition at line 679 of file DDHGCalMixRotatedCassette.cc.

◆ passiveFull_

std::vector<std::string> HGCalMixRotatedCassette::passiveFull_

Definition at line 699 of file DDHGCalMixRotatedCassette.cc.

◆ passivePart_

std::vector<std::string> HGCalMixRotatedCassette::passivePart_

Definition at line 700 of file DDHGCalMixRotatedCassette.cc.

◆ passiveTypes_

int HGCalMixRotatedCassette::passiveTypes_

Definition at line 676 of file DDHGCalMixRotatedCassette.cc.

◆ phiBinsScint_

int HGCalMixRotatedCassette::phiBinsScint_

Definition at line 681 of file DDHGCalMixRotatedCassette.cc.

◆ placeOffset_

int HGCalMixRotatedCassette::placeOffset_

Definition at line 680 of file DDHGCalMixRotatedCassette.cc.

◆ rMaxFront_

std::vector<double> HGCalMixRotatedCassette::rMaxFront_

Definition at line 696 of file DDHGCalMixRotatedCassette.cc.

◆ rMinFront_

std::vector<double> HGCalMixRotatedCassette::rMinFront_

Definition at line 693 of file DDHGCalMixRotatedCassette.cc.

◆ sectors_

int HGCalMixRotatedCassette::sectors_

Definition at line 689 of file DDHGCalMixRotatedCassette.cc.

◆ sensitiveMode_

int HGCalMixRotatedCassette::sensitiveMode_

Definition at line 685 of file DDHGCalMixRotatedCassette.cc.

◆ slopeB_

std::vector<double> HGCalMixRotatedCassette::slopeB_

Definition at line 691 of file DDHGCalMixRotatedCassette.cc.

◆ slopeT_

std::vector<double> HGCalMixRotatedCassette::slopeT_

Definition at line 694 of file DDHGCalMixRotatedCassette.cc.

◆ thick_

std::vector<double> HGCalMixRotatedCassette::thick_

Definition at line 703 of file DDHGCalMixRotatedCassette.cc.

◆ tileIndex_

std::vector<int> HGCalMixRotatedCassette::tileIndex_

Definition at line 723 of file DDHGCalMixRotatedCassette.cc.

◆ tileLayerStart_

std::vector<int> HGCalMixRotatedCassette::tileLayerStart_

Definition at line 725 of file DDHGCalMixRotatedCassette.cc.

◆ tilePhis_

std::vector<int> HGCalMixRotatedCassette::tilePhis_

Definition at line 724 of file DDHGCalMixRotatedCassette.cc.

◆ tileRMax_

std::vector<double> HGCalMixRotatedCassette::tileRMax_

Definition at line 722 of file DDHGCalMixRotatedCassette.cc.

◆ tileRMin_

std::vector<double> HGCalMixRotatedCassette::tileRMin_

Definition at line 721 of file DDHGCalMixRotatedCassette.cc.

◆ tol2_

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

Definition at line 673 of file DDHGCalMixRotatedCassette.cc.

◆ waferFull_

std::vector<std::string> HGCalMixRotatedCassette::waferFull_

Definition at line 697 of file DDHGCalMixRotatedCassette.cc.

◆ waferIndex_

std::vector<int> HGCalMixRotatedCassette::waferIndex_

Definition at line 717 of file DDHGCalMixRotatedCassette.cc.

◆ waferLayerStart_

std::vector<int> HGCalMixRotatedCassette::waferLayerStart_

Definition at line 719 of file DDHGCalMixRotatedCassette.cc.

◆ waferPart_

std::vector<std::string> HGCalMixRotatedCassette::waferPart_

Definition at line 698 of file DDHGCalMixRotatedCassette.cc.

◆ waferProperty_

std::vector<int> HGCalMixRotatedCassette::waferProperty_

Definition at line 718 of file DDHGCalMixRotatedCassette.cc.

◆ waferSepar_

double HGCalMixRotatedCassette::waferSepar_

Definition at line 688 of file DDHGCalMixRotatedCassette.cc.

◆ waferSize_

double HGCalMixRotatedCassette::waferSize_

Definition at line 687 of file DDHGCalMixRotatedCassette.cc.

◆ waferTypes_

int HGCalMixRotatedCassette::waferTypes_

Definition at line 675 of file DDHGCalMixRotatedCassette.cc.

◆ zFrontB_

std::vector<double> HGCalMixRotatedCassette::zFrontB_

Definition at line 692 of file DDHGCalMixRotatedCassette.cc.

◆ zFrontT_

std::vector<double> HGCalMixRotatedCassette::zFrontT_

Definition at line 695 of file DDHGCalMixRotatedCassette.cc.

◆ zMinBlock_

double HGCalMixRotatedCassette::zMinBlock_

Definition at line 686 of file DDHGCalMixRotatedCassette.cc.

Generated for CMSSW Reference Manual by doxygen 1.8.14