CMS 3D CMS Logo

List of all members | Public Member Functions | Protected Member Functions | Private Attributes | Static Private Attributes
DDHGCalMixRotatedCassette Class Reference
Inheritance diagram for DDHGCalMixRotatedCassette:

Public Member Functions

 DDHGCalMixRotatedCassette ()
 
void execute (DDCompactView &cpv) override
 
void initialize (const DDNumericArguments &nArgs, const DDVectorArguments &vArgs, const DDMapArguments &mArgs, const DDStringArguments &sArgs, const DDStringVectorArguments &vsArgs) override
 

Protected Member Functions

void constructLayers (const DDLogicalPart &, DDCompactView &cpv)
 
void positionMix (const DDLogicalPart &glog, const std::string &nameM, int copyM, double thick, const DDMaterial &matter, int absType, DDCompactView &cpv)
 

Private 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_
 
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::string nameSpace_
 
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 Private Attributes

static constexpr double tol1_ = 0.01
 
static constexpr double tol2_ = 0.00001
 

Detailed Description

Definition at line 37 of file DDHGCalMixRotatedCassette.cc.

Constructor & Destructor Documentation

◆ DDHGCalMixRotatedCassette()

DDHGCalMixRotatedCassette::DDHGCalMixRotatedCassette ( )

Definition at line 120 of file DDHGCalMixRotatedCassette.cc.

120  {
121 #ifdef EDM_ML_DEBUG
122  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Creating an instance";
123 #endif
124 }
Log< level::Info, true > LogVerbatim

Member Function Documentation

◆ constructLayers()

void DDHGCalMixRotatedCassette::constructLayers ( const DDLogicalPart module,
DDCompactView cpv 
)
protected

Definition at line 366 of file DDHGCalMixRotatedCassette.cc.

References funct::abs(), angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, DDBase< N, C >::ddname(), DDSplit(), dqmdumpme::first, mps_fire::i, cuy::ii, l1ctLayer1_patternWriters_cff::isec, dqmdumpme::k, SiStripPI::min, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, Skims_PA_cff::name, DDBase< N, C >::name(), DDSolidFactory::polyhedra(), DDCompactView::position(), HGCalGeomTools::radius(), makeMuonMisalignmentScenario::rot, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, to_string(), DDSolidFactory::tubs(), and geometryCSVtoXML::zz.

366  {
367  double zi(zMinBlock_);
368  int laymin(0);
369  for (unsigned int i = 0; i < layers_.size(); i++) {
370  double zo = zi + layerThick_[i];
371  double routF = HGCalGeomTools::radius(zi, zFrontT_, rMaxFront_, slopeT_);
372  int laymax = laymin + layers_[i];
373  double zz = zi;
374  double thickTot(0);
375  for (int ly = laymin; ly < laymax; ++ly) {
376  int ii = layerType_[ly];
377  int copy = copyNumber_[ii];
378  double hthick = 0.5 * thick_[ii];
379  double rinB = HGCalGeomTools::radius(zo, zFrontB_, rMinFront_, slopeB_);
380  zz += hthick;
381  thickTot += thick_[ii];
382 
384 #ifdef EDM_ML_DEBUG
385  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Layer " << ly << ":" << ii << " Front " << zi << ", "
386  << routF << " Back " << zo << ", " << rinB << " superlayer thickness "
387  << layerThick_[i];
388 #endif
390  DDMaterial matter(matName);
391  DDLogicalPart glog;
392  if (layerSense_[ly] == 0) {
393  std::vector<double> pgonZ, pgonRin, pgonRout;
394  double rmax =
396  HGCalGeomTools::radius(zz - hthick,
397  zz + hthick,
398  zFrontB_,
399  rMinFront_,
400  slopeB_,
401  zFrontT_,
402  rMaxFront_,
403  slopeT_,
404  -layerSense_[ly],
405  pgonZ,
406  pgonRin,
407  pgonRout);
408  for (unsigned int isec = 0; isec < pgonZ.size(); ++isec) {
409  pgonZ[isec] -= zz;
410  if (layerSense_[ly] == 0 || absorbMode_ == 0)
411  pgonRout[isec] = rmax;
412  else
413  pgonRout[isec] = pgonRout[isec] * cosAlpha_ - tol1_;
414  }
415  DDSolid solid =
416  DDSolidFactory::polyhedra(DDName(name, nameSpace_), sectors_, -alpha_, 2._pi, pgonZ, pgonRin, pgonRout);
417  glog = DDLogicalPart(solid.ddname(), matter, solid);
418 #ifdef EDM_ML_DEBUG
419  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << solid.name() << " polyhedra of " << sectors_
420  << " sectors covering " << convertRadToDeg(-alpha_) << ":"
421  << convertRadToDeg(-alpha_ + 2._pi) << " with " << pgonZ.size() << " sections";
422  for (unsigned int k = 0; k < pgonZ.size(); ++k)
423  edm::LogVerbatim("HGCalGeom") << "[" << k << "] z " << pgonZ[k] << " R " << pgonRin[k] << ":" << pgonRout[k];
424 #endif
425  } else {
426  int mode = (layerSense_[ly] > 0) ? sensitiveMode_ : absorbMode_;
427  double rins = (mode < 1) ? rinB : HGCalGeomTools::radius(zz + hthick, zFrontB_, rMinFront_, slopeB_);
428  double routs = (mode < 1) ? routF : HGCalGeomTools::radius(zz - hthick, zFrontT_, rMaxFront_, slopeT_);
429  DDSolid solid = DDSolidFactory::tubs(DDName(name, nameSpace_), hthick, rins, routs, 0.0, 2._pi);
430  glog = DDLogicalPart(solid.ddname(), matter, solid);
431 #ifdef EDM_ML_DEBUG
432  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << solid.name() << " Tubs made of " << matName
433  << " of dimensions " << rinB << ":" << rins << ", " << routF << ":" << routs
434  << ", " << hthick << ", 0.0, 360.0 and positioned in: " << glog.name()
435  << " number " << copy;
436 #endif
437  positionMix(glog, name, copy, thick_[ii], matter, -layerSense_[ly], cpv);
438  }
439  DDTranslation r1(0, 0, zz);
440  DDRotation rot;
441  cpv.position(glog, module, copy, r1, rot);
442  ++copyNumber_[ii];
443 #ifdef EDM_ML_DEBUG
444  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << glog.name() << " number " << copy
445  << " positioned in " << module.name() << " at " << r1 << " with no rotation";
446 #endif
447  zz += hthick;
448  } // End of loop over layers in a block
449  zi = zo;
450  laymin = laymax;
451  // Make consistency check of all the partitions of the block
452  if (std::abs(thickTot - layerThick_[i]) >= tol2_) {
453  if (thickTot > layerThick_[i]) {
454  edm::LogError("HGCalGeom") << "Thickness of the partition " << layerThick_[i] << " is smaller than " << thickTot
455  << ": thickness of all its components **** ERROR ****";
456  } else {
457  edm::LogWarning("HGCalGeom") << "Thickness of the partition " << layerThick_[i] << " does not match with "
458  << thickTot << " of the components";
459  }
460  }
461  } // End of loop over blocks
462 }
Log< level::Info, true > LogVerbatim
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)
void position(const DDLogicalPart &self, const DDLogicalPart &parent, const std::string &copyno, const DDTranslation &trans, const DDRotation &rot, const DDDivision *div=nullptr)
DDMaterial is used to define and access material information.
Definition: DDMaterial.h:45
void positionMix(const DDLogicalPart &glog, const std::string &nameM, int copyM, double thick, const DDMaterial &matter, int absType, DDCompactView &cpv)
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
DDName is used to identify DDD entities uniquely.
Definition: DDName.h:17
std::vector< std::string > names_
Log< level::Error, false > LogError
A DDSolid represents the shape of a part.
Definition: DDSolid.h:39
static std::string to_string(const XMLCh *ch)
Represents a uniquely identifyable rotation matrix.
Definition: DDTransform.h:57
U second(std::pair< T, U > const &p)
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
A DDLogicalPart aggregates information concerning material, solid and sensitveness ...
Definition: DDLogicalPart.h:93
static DDSolid tubs(const DDName &name, double zhalf, double rIn, double rOut, double startPhi, double deltaPhi)
Definition: DDSolid.cc:667
ii
Definition: cuy.py:589
const N & name() const
Definition: DDBase.h:59
const N & ddname() const
Definition: DDBase.h:61
Log< level::Warning, false > LogWarning
std::vector< std::string > materials_
std::pair< std::string, std::string > DDSplit(const std::string &n)
split into (name,namespace), separator = &#39;:&#39;
Definition: DDSplit.cc:3
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DDTranslation
Definition: DDTranslation.h:7
static DDSolid polyhedra(const DDName &name, int sides, double startPhi, double deltaPhi, const std::vector< double > &z, const std::vector< double > &rmin, const std::vector< double > &rmax)
Creates a polyhedra (refere to Geant3 or Geant4 documentation)
Definition: DDSolid.cc:565

◆ execute()

void DDHGCalMixRotatedCassette::execute ( DDCompactView cpv)
override

Definition at line 349 of file DDHGCalMixRotatedCassette.cc.

References DDHGCalGeom::constructLayers(), dqmdumpme::k, and class-composition::parent.

349  {
350 #ifdef EDM_ML_DEBUG
351  edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalMixRotatedCassette...";
352  copies_.clear();
353 #endif
354  constructLayers(parent(), cpv);
355 #ifdef EDM_ML_DEBUG
356  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << copies_.size() << " different wafer copy numbers";
357  int k(0);
358  for (std::unordered_set<int>::const_iterator itr = copies_.begin(); itr != copies_.end(); ++itr, ++k) {
359  edm::LogVerbatim("HGCalGeom") << "Copy [" << k << "] : " << (*itr);
360  }
361  copies_.clear();
362  edm::LogVerbatim("HGCalGeom") << "<<== End of DDHGCalMixRotatedCassette construction...";
363 #endif
364 }
Log< level::Info, true > LogVerbatim
std::unordered_set< int > copies_
void constructLayers(const DDLogicalPart &, DDCompactView &cpv)

◆ initialize()

void DDHGCalMixRotatedCassette::initialize ( const DDNumericArguments nArgs,
const DDVectorArguments vArgs,
const DDMapArguments mArgs,
const DDStringArguments sArgs,
const DDStringVectorArguments vsArgs 
)
override

Definition at line 126 of file DDHGCalMixRotatedCassette.cc.

References angle_units::operators::convertRadToDeg(), funct::cos(), dbl_to_int(), mps_fire::i, testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, cuy::ii, createfilelist::int, dqmiolumiharvest::j, dqmdumpme::k, HGCalTypes::layerType(), SiStripPI::min, DDCurrentNamespace::ns(), HGCalTileIndex::tileUnpack(), HGCalWaferIndex::waferLayer(), HGCalProperty::waferOrient(), HGCalProperty::waferPartial(), HGCalProperty::waferThick(), HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

130  {
131  waferTypes_ = static_cast<int>(nArgs["WaferTypes"]);
132  passiveTypes_ = static_cast<int>(nArgs["PassiveTypes"]);
133  facingTypes_ = static_cast<int>(nArgs["FacingTypes"]);
134  orientationTypes_ = static_cast<int>(nArgs["OrientationTypes"]);
135  partialTypes_ = static_cast<int>(nArgs["PartialTypes"]);
136  placeOffset_ = static_cast<int>(nArgs["PlaceOffset"]);
137  phiBinsScint_ = static_cast<int>(nArgs["NPhiBinScint"]);
138 #ifdef EDM_ML_DEBUG
139  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette::Number of types of wafers: " << waferTypes_
140  << " passives: " << passiveTypes_ << " facings: " << facingTypes_
141  << " Orientations: " << orientationTypes_ << " PartialTypes: " << partialTypes_
142  << " PlaceOffset: " << placeOffset_ << "; number of cells along phi " << phiBinsScint_;
143 #endif
144  firstLayer_ = (int)(nArgs["FirstLayer"]);
145  absorbMode_ = (int)(nArgs["AbsorberMode"]);
146  sensitiveMode_ = (int)(nArgs["SensitiveMode"]);
147 #ifdef EDM_ML_DEBUG
148  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette::First Layer " << firstLayer_ << " and "
149  << "Absober:Sensitive mode " << absorbMode_ << ":" << sensitiveMode_;
150 #endif
151  zMinBlock_ = nArgs["zMinBlock"];
152  waferSize_ = nArgs["waferSize"];
153  waferSepar_ = nArgs["SensorSeparation"];
154  sectors_ = static_cast<int>(nArgs["Sectors"]);
155  cassettes_ = static_cast<int>(nArgs["Cassettes"]);
156  alpha_ = (1._pi) / sectors_;
157  cosAlpha_ = cos(alpha_);
158 #ifdef EDM_ML_DEBUG
159  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: zStart " << zMinBlock_ << " wafer width " << waferSize_
160  << " separations " << waferSepar_ << " sectors " << sectors_ << ":"
161  << convertRadToDeg(alpha_) << ":" << cosAlpha_ << " with " << cassettes_
162  << " cassettes";
163 #endif
164  slopeB_ = vArgs["SlopeBottom"];
165  zFrontB_ = vArgs["ZFrontBottom"];
166  rMinFront_ = vArgs["RMinFront"];
167  slopeT_ = vArgs["SlopeTop"];
168  zFrontT_ = vArgs["ZFrontTop"];
169  rMaxFront_ = vArgs["RMaxFront"];
170 #ifdef EDM_ML_DEBUG
171  for (unsigned int i = 0; i < slopeB_.size(); ++i)
172  edm::LogVerbatim("HGCalGeom") << "Bottom Block [" << i << "] Zmin " << zFrontB_[i] << " Rmin " << rMinFront_[i]
173  << " Slope " << slopeB_[i];
174  for (unsigned int i = 0; i < slopeT_.size(); ++i)
175  edm::LogVerbatim("HGCalGeom") << "Top Block [" << i << "] Zmin " << zFrontT_[i] << " Rmax " << rMaxFront_[i]
176  << " Slope " << slopeT_[i];
177 #endif
178  waferFull_ = vsArgs["WaferNamesFull"];
179  waferPart_ = vsArgs["WaferNamesPartial"];
180 #ifdef EDM_ML_DEBUG
181  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << waferFull_.size() << " full and "
182  << waferPart_.size() << " partial modules\nDDHGCalMixRotatedCassette:Full Modules:";
183  unsigned int i1max = static_cast<unsigned int>(waferFull_.size());
184  for (unsigned int i1 = 0; i1 < i1max; i1 += 2) {
185  std::ostringstream st1;
186  unsigned int i2 = std::min((i1 + 2), i1max);
187  for (unsigned int i = i1; i < i2; ++i)
188  st1 << " [" << i << "] " << waferFull_[i];
189  edm::LogVerbatim("HGCalGeom") << st1.str();
190  }
191  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Partial Modules:";
192  i1max = static_cast<unsigned int>(waferPart_.size());
193  for (unsigned int i1 = 0; i1 < i1max; i1 += 2) {
194  std::ostringstream st1;
195  unsigned int i2 = std::min((i1 + 2), i1max);
196  for (unsigned int i = i1; i < i2; ++i)
197  st1 << " [" << i << "] " << waferPart_[i];
198  edm::LogVerbatim("HGCalGeom") << st1.str();
199  }
200 #endif
201  passiveFull_ = vsArgs["PassiveNamesFull"];
202  passivePart_ = vsArgs["PassiveNamesPartial"];
203 #ifdef EDM_ML_DEBUG
204  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconRotatedCassette: " << passiveFull_.size() << " full and "
205  << passivePart_.size() << " partial passive modules";
206  i1max = static_cast<unsigned int>(passiveFull_.size());
207  for (unsigned int i1 = 0; i1 < i1max; i1 += 2) {
208  std::ostringstream st1;
209  unsigned int i2 = std::min((i1 + 2), i1max);
210  for (unsigned int i = i1; i < i2; ++i)
211  st1 << " [" << i << "] " << passiveFull_[i];
212  edm::LogVerbatim("HGCalGeom") << st1.str();
213  }
214  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconRotatedCassette: Partial Modules:";
215  i1max = static_cast<unsigned int>(passivePart_.size());
216  for (unsigned int i1 = 0; i1 < i1max; i1 += 2) {
217  std::ostringstream st1;
218  unsigned int i2 = std::min((i1 + 2), i1max);
219  for (unsigned int i = i1; i < i2; ++i)
220  st1 << " [" << i << "] " << passivePart_[i];
221  edm::LogVerbatim("HGCalGeom") << st1.str();
222  }
223 #endif
224  materials_ = vsArgs["MaterialNames"];
225  names_ = vsArgs["VolumeNames"];
226  thick_ = vArgs["Thickness"];
227  copyNumber_.resize(materials_.size(), 1);
228 #ifdef EDM_ML_DEBUG
229  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << materials_.size() << " types of volumes";
230  for (unsigned int i = 0; i < names_.size(); ++i)
231  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << names_[i] << " of thickness " << thick_[i]
232  << " filled with " << materials_[i] << " first copy number " << copyNumber_[i];
233 #endif
234  layers_ = dbl_to_int(vArgs["Layers"]);
235  layerThick_ = vArgs["LayerThick"];
236 #ifdef EDM_ML_DEBUG
237  edm::LogVerbatim("HGCalGeom") << "There are " << layers_.size() << " blocks";
238  for (unsigned int i = 0; i < layers_.size(); ++i)
239  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] of thickness " << layerThick_[i] << " with " << layers_[i]
240  << " layers";
241 #endif
242  layerType_ = dbl_to_int(vArgs["LayerType"]);
243  layerSense_ = dbl_to_int(vArgs["LayerSense"]);
244  layerOrient_ = dbl_to_int(vArgs["LayerTypes"]);
245  for (unsigned int k = 0; k < layerOrient_.size(); ++k)
247 #ifdef EDM_ML_DEBUG
248  for (unsigned int i = 0; i < layerOrient_.size(); ++i)
249  edm::LogVerbatim("HGCalGeom") << "LayerOrient [" << i << "] " << layerOrient_[i];
250 #endif
251  if (firstLayer_ > 0) {
252  for (unsigned int i = 0; i < layerType_.size(); ++i) {
253  if (layerSense_[i] != 0) {
254  int ii = layerType_[i];
256 #ifdef EDM_ML_DEBUG
257  edm::LogVerbatim("HGCalGeom") << "First copy number for layer type " << i << ":" << ii << " with "
258  << materials_[ii] << " changed to " << copyNumber_[ii];
259 #endif
260  }
261  }
262  } else {
263  firstLayer_ = 1;
264  }
265 #ifdef EDM_ML_DEBUG
266  edm::LogVerbatim("HGCalGeom") << "There are " << layerType_.size() << " layers";
267  for (unsigned int i = 0; i < layerType_.size(); ++i)
268  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerType_[i] << " sensitive class "
269  << layerSense_[i];
270 #endif
271  materialTop_ = vsArgs["TopMaterialNames"];
272  namesTop_ = vsArgs["TopVolumeNames"];
273  layerThickTop_ = vArgs["TopLayerThickness"];
274  layerTypeTop_ = dbl_to_int(vArgs["TopLayerType"]);
275  copyNumberTop_.resize(materialTop_.size(), firstLayer_);
276  coverTypeTop_ = static_cast<int>(nArgs["TopCoverLayerType"]);
278 #ifdef EDM_ML_DEBUG
279  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << materialTop_.size()
280  << " types of volumes in the top part; cover Type " << coverTypeTop_
281  << " with initial copy number " << copyNumberCoverTop_;
282  for (unsigned int i = 0; i < materialTop_.size(); ++i)
283  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << namesTop_[i] << " of thickness " << layerThickTop_[i]
284  << " filled with " << materialTop_[i] << " first copy number " << copyNumberTop_[i];
285  edm::LogVerbatim("HGCalGeom") << "There are " << layerTypeTop_.size() << " layers in the top part";
286  for (unsigned int i = 0; i < layerTypeTop_.size(); ++i)
287  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerTypeTop_[i];
288 #endif
289  waferIndex_ = dbl_to_int(vArgs["WaferIndex"]);
290  waferProperty_ = dbl_to_int(vArgs["WaferProperties"]);
291  waferLayerStart_ = dbl_to_int(vArgs["WaferLayerStart"]);
292  cassetteShift_ = vArgs["CassetteShift"];
293 #ifdef EDM_ML_DEBUG
294  edm::LogVerbatim("HGCalGeom") << "waferProperties with " << waferIndex_.size() << " entries in "
295  << waferLayerStart_.size() << " layers";
296  for (unsigned int k = 0; k < waferLayerStart_.size(); ++k)
297  edm::LogVerbatim("HGCalGeom") << "LayerStart[" << k << "] " << waferLayerStart_[k];
298  for (unsigned int k = 0; k < waferIndex_.size(); ++k)
299  edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << waferIndex_[k] << " ("
302  << HGCalWaferIndex::waferV(waferIndex_[k]) << ") : ("
306  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << cassetteShift_.size()
307  << " elements for cassette shifts";
308  unsigned int j1max = cassetteShift_.size();
309  for (unsigned int j1 = 0; j1 < j1max; j1 += 6) {
310  std::ostringstream st1;
311  unsigned int j2 = std::min((j1 + 6), j1max);
312  for (unsigned int j = j1; j < j2; ++j)
313  st1 << " [" << j << "] " << std::setw(9) << cassetteShift_[j];
314  edm::LogVerbatim("HGCalGeom") << st1.str();
315  }
316 #endif
317  tileRMin_ = vArgs["TileRMin"];
318  tileRMax_ = vArgs["TileRMax"];
319  tileIndex_ = dbl_to_int(vArgs["TileLayerRings"]);
320  tilePhis_ = dbl_to_int(vArgs["TilePhiRange"]);
321  tileLayerStart_ = dbl_to_int(vArgs["TileLayerStart"]);
322 #ifdef EDM_ML_DEBUG
323  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette:: with " << tileRMin_.size() << " rings";
324  for (unsigned int k = 0; k < tileRMin_.size(); ++k)
325  edm::LogVerbatim("HGCalGeom") << "Ring[" << k << "] " << tileRMin_[k] << " : " << tileRMax_[k];
326  edm::LogVerbatim("HGCalGeom") << "TileProperties with " << tileIndex_.size() << " entries in "
327  << tileLayerStart_.size() << " layers";
328  for (unsigned int k = 0; k < tileLayerStart_.size(); ++k)
329  edm::LogVerbatim("HGCalGeom") << "LayerStart[" << k << "] " << tileLayerStart_[k];
330  for (unsigned int k = 0; k < tileIndex_.size(); ++k)
331  edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << tileIndex_[k] << " ("
332  << "Layer " << std::get<0>(HGCalTileIndex::tileUnpack(tileIndex_[k])) << " Ring "
333  << std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[k])) << ":"
334  << std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[k])) << ") Phi "
335  << std::get<1>(HGCalTileIndex::tileUnpack(tilePhis_[k])) << ":"
336  << std::get<2>(HGCalTileIndex::tileUnpack(tilePhis_[k]));
337 #endif
339 #ifdef EDM_ML_DEBUG
340  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: NameSpace " << nameSpace_ << ":";
341 #endif
343 }
Log< level::Info, true > LogVerbatim
int32_t waferU(const int32_t index)
int32_t waferLayer(const int32_t index)
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
std::vector< std::string > names_
static std::string & ns()
int32_t waferOrient(const int32_t property)
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
std::vector< std::string > namesTop_
std::vector< std::string > waferFull_
int32_t waferThick(const int32_t property)
void setParameter(int cassette, const std::vector< double > &shifts, bool both=true)
Definition: HGCalCassette.cc:8
ii
Definition: cuy.py:589
std::vector< int > dbl_to_int(const std::vector< double > &vecdbl)
Converts a std::vector of doubles to a std::vector of int.
Definition: DDutils.h:7
std::tuple< int32_t, int32_t, int32_t > tileUnpack(int32_t index)
std::vector< std::string > passivePart_
std::vector< std::string > materialTop_
std::vector< std::string > passiveFull_
int32_t waferPartial(const int32_t property)
int32_t waferV(const int32_t index)
static int32_t layerType(int type)
Definition: HGCalTypes.cc:42
std::vector< std::string > materials_
std::vector< std::string > waferPart_

◆ positionMix()

void DDHGCalMixRotatedCassette::positionMix ( const DDLogicalPart glog,
const std::string &  nameM,
int  copyM,
double  thick,
const DDMaterial matter,
int  absType,
DDCompactView cpv 
)
protected

Definition at line 464 of file DDHGCalMixRotatedCassette.cc.

References funct::abs(), angle(), HGCalCassette::cassetteType(), HGCalCell::cellPlacementIndex(), angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, DDBase< N, C >::ddname(), DDSplit(), PVValHelper::dy, dqmdumpme::first, mps_fire::i, cuy::ii, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, dqmdumpme::k, HGCalTypes::layerFrontBack(), Skims_PA_cff::name, DDBase< N, C >::name(), EgHLTOffHistBins_cfi::nr, HGCalTypes::packTypeUV(), DDCompactView::position(), diffTwoXMLs::r2, makeMuonMisalignmentScenario::rot, edm::second(), mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, HGCalTileIndex::tileUnpack(), to_string(), DDSolidFactory::tubs(), findQualityFiles::v, 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().

470  {
471  DDRotation rot;
472 
473  // Make the top part first
474  for (unsigned int ly = 0; ly < layerTypeTop_.size(); ++ly) {
475  int ii = layerTypeTop_[ly];
476  copyNumberTop_[ii] = copyM;
477  }
478  double hthick = 0.5 * thick;
479  double dphi = (2._pi) / phiBinsScint_;
480  double thickTot(0), zpos(-hthick);
481  if (absType < 0) {
482  for (unsigned int ly = 0; ly < layerTypeTop_.size(); ++ly) {
483  int ii = layerTypeTop_[ly];
484  int copy = copyNumberTop_[ii];
485  int layer = copy - firstLayer_;
486  double hthickl = 0.5 * layerThickTop_[ii];
487  thickTot += layerThickTop_[ii];
488  zpos += hthickl;
490  DDMaterial matter1(matName);
491  unsigned int k = 0;
492  int firstTile = tileLayerStart_[layer];
493  int lastTile = ((layer + 1 < static_cast<int>(tileLayerStart_.size())) ? tileLayerStart_[layer + 1]
494  : static_cast<int>(tileIndex_.size()));
495 #ifdef EDM_ML_DEBUG
496  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Layer " << ly << ":" << ii << " Copy " << copy
497  << " Tiles " << firstTile << ":" << lastTile;
498 #endif
499  for (int ti = firstTile; ti < lastTile; ++ti) {
500  double r1 = tileRMin_[std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
501  double r2 = tileRMax_[std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
502  int cassette = std::get<0>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
503  int fimin = std::get<1>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
504  int fimax = std::get<2>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
505  double phi1 = dphi * (fimin - 1);
506  double phi2 = dphi * (fimax - fimin + 1);
507  auto cshift = cassette_.getShift(layer + 1, 1, cassette);
508 #ifdef EDM_ML_DEBUG
509  int cassette0 = HGCalCassette::cassetteType(2, 1, cassette); //
510  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Layer " << copy << ":" << (layer + 1) << " iR "
511  << std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) << ":"
512  << std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) << " R " << r1 << ":"
513  << r2 << " Thick " << (2.0 * hthickl) << " phi " << fimin << ":" << fimax << ":"
514  << convertRadToDeg(phi1) << ":" << convertRadToDeg(phi2) << " cassette "
515  << cassette << ":" << cassette0 << " Shift " << cshift.first << ":"
516  << cshift.second;
517 #endif
519  ++k;
520  DDSolid solid = DDSolidFactory::tubs(DDName(name, nameSpace_), hthickl, r1, r2, phi1, phi2);
521  DDLogicalPart glog1 = DDLogicalPart(solid.ddname(), matter1, solid);
522 #ifdef EDM_ML_DEBUG
523  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << glog1.name() << " Tubs made of " << matName
524  << " of dimensions " << r1 << ", " << r2 << ", " << hthickl << ", "
525  << convertRadToDeg(phi1) << ", " << convertRadToDeg(phi2);
526 #endif
527  DDTranslation tran(-cshift.first, cshift.second, zpos);
528  cpv.position(glog1, glog, copy, tran, rot);
529 #ifdef EDM_ML_DEBUG
530  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Position " << glog1.name() << " number " << copy
531  << " in " << glog.name() << " at " << tran << " with no rotation";
532 #endif
533  }
534  ++copyNumberTop_[ii];
535  zpos += hthickl;
536  }
537  if (std::abs(thickTot - thick) >= tol2_) {
538  if (thickTot > thick) {
539  edm::LogError("HGCalGeom") << "DDHGCalMixRotatedCassette: Thickness of the partition " << thick
540  << " is smaller than " << thickTot
541  << ": thickness of all its components in the top part **** ERROR ****";
542  } else {
543  edm::LogWarning("HGCalGeom") << "DDHGCalMixRotatedCassette: Thickness of the partition " << thick
544  << " does not match with " << thickTot << " of the components in top part";
545  }
546  }
547  } else {
548  int ii = coverTypeTop_;
550  int layer = copy - firstLayer_;
551  double hthickl = 0.5 * layerThickTop_[ii];
552  zpos += hthickl;
554  DDMaterial matter1(matName);
555  unsigned int k = 0;
556  int firstTile = tileLayerStart_[layer];
557  int lastTile = ((layer + 1 < static_cast<int>(tileLayerStart_.size())) ? tileLayerStart_[layer + 1]
558  : static_cast<int>(tileIndex_.size()));
559 #ifdef EDM_ML_DEBUG
560  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Layer " << layer << ":" << ii << " Copy " << copy
561  << " Tiles " << firstTile << ":" << lastTile;
562 #endif
563  for (int ti = firstTile; ti < lastTile; ++ti) {
564  double r1 = tileRMin_[std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
565  double r2 = tileRMax_[std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) - 1];
566  int cassette = std::get<0>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
567  int fimin = std::get<1>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
568  int fimax = std::get<2>(HGCalTileIndex::tileUnpack(tilePhis_[ti]));
569  double phi1 = dphi * (fimin - 1);
570  double phi2 = dphi * (fimax - fimin + 1);
571  auto cshift = cassette_.getShift(layer + 1, 1, cassette);
572 #ifdef EDM_ML_DEBUG
573  int cassette0 = HGCalCassette::cassetteType(2, 1, cassette); //
574  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Layer " << copy << ":" << (layer + 1) << " iR "
575  << std::get<1>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) << ":"
576  << std::get<2>(HGCalTileIndex::tileUnpack(tileIndex_[ti])) << " R " << r1 << ":"
577  << r2 << " Thick " << (2.0 * hthickl) << " phi " << fimin << ":" << fimax << ":"
578  << convertRadToDeg(phi1) << ":" << convertRadToDeg(phi2) << " cassette " << cassette
579  << ":" << cassette0 << " Shift " << cshift.first << ":" << cshift.second;
580 #endif
582  ++k;
583  DDSolid solid = DDSolidFactory::tubs(DDName(name, nameSpace_), hthickl, r1, r2, phi1, phi2);
584  DDLogicalPart glog1 = DDLogicalPart(solid.ddname(), matter1, solid);
585 #ifdef EDM_ML_DEBUG
586  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: " << glog1.name() << " Tubs made of " << matName
587  << " of dimensions " << r1 << ", " << r2 << ", " << hthickl << ", "
588  << convertRadToDeg(phi1) << ", " << convertRadToDeg(phi2);
589 #endif
590  DDTranslation tran(-cshift.first, cshift.second, zpos);
591  cpv.position(glog1, glog, copy, tran, rot);
592 #ifdef EDM_ML_DEBUG
593  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Position " << glog1.name() << " number " << copy
594  << " in " << glog.name() << " at " << tran << " with no rotation";
595 #endif
596  }
598  }
599 
600  // Make the bottom part next
601  int layer = (copyM - firstLayer_);
602  static const double sqrt3 = std::sqrt(3.0);
603  int layercenter = layerOrient_[layer];
605  int firstWafer = waferLayerStart_[layer];
606  int lastWafer = ((layer + 1 < static_cast<int>(waferLayerStart_.size())) ? waferLayerStart_[layer + 1]
607  : static_cast<int>(waferIndex_.size()));
608  double delx = 0.5 * (waferSize_ + waferSepar_);
609  double dely = 2.0 * delx / sqrt3;
610  double dy = 0.75 * dely;
611  const auto& xyoff = geomTools_.shiftXY(layercenter, (waferSize_ + waferSepar_));
612 #ifdef EDM_ML_DEBUG
613  int ium(0), ivm(0), kount(0);
614  std::vector<int> ntype(3, 0);
615  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette::Bottom: " << glog.ddname() << " r " << delx << " R "
616  << dely << " dy " << dy << " Shift " << xyoff.first << ":" << xyoff.second
617  << " WaferSize " << (waferSize_ + waferSepar_) << " index " << firstWafer << ":"
618  << (lastWafer - 1) << " Copy " << copyM << ":" << layer;
619 #endif
620  for (int k = firstWafer; k < lastWafer; ++k) {
623 #ifdef EDM_ML_DEBUG
624  int iu = std::abs(u);
625  int iv = std::abs(v);
626 #endif
627  int nr = 2 * v;
628  int nc = -2 * u + v;
633  int place = HGCalCell::cellPlacementIndex(1, layertype, orien);
634 #ifdef EDM_ML_DEBUG
635  edm::LogVerbatim("HGCalGeom")
636  << "DDHGCalMixRotatedCassette::index:Property:layertype:type:part:orien:cassette:place:offsets:ind " << k << ":"
637  << waferProperty_[k] << ":" << layertype << ":" << type << ":" << part << ":" << orien << ":" << cassette << ":"
638  << place;
639 #endif
640  auto cshift = cassette_.getShift(layer + 1, -1, cassette);
641  double xpos = xyoff.first - cshift.first + nc * delx;
642  double ypos = xyoff.second + cshift.second + nr * dy;
643 #ifdef EDM_ML_DEBUG
644  double xorig = xyoff.first + nc * delx;
645  double yorig = xyoff.second + nr * dy;
646  double angle = std::atan2(yorig, xorig);
647  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette::Wafer: layer " << layer + 1 << " cassette " << cassette
648  << " Shift " << cshift.first << ":" << cshift.second << " Original " << xorig << ":"
649  << yorig << ":" << convertRadToDeg(angle) << " Final " << xpos << ":" << ypos;
650 #endif
651  std::string wafer;
652  int i(999);
653  if (absType < 0) {
654  if (part == HGCalTypes::WaferFull) {
656 #ifdef EDM_ML_DEBUG
657  edm::LogVerbatim("HGCalGeom") << " FullWafer type:place:ind " << type << ":" << place << ":" << i << ":"
658  << waferFull_.size();
659 #endif
660  wafer = waferFull_[i];
661  } else {
662  int partoffset =
664  i = (part - partoffset) * facingTypes_ * orientationTypes_ +
666 #ifdef EDM_ML_DEBUG
667  edm::LogVerbatim("HGCalGeom") << " layertype:type:part:orien:cassette:place:offsets:ind " << layertype << ":"
668  << type << ":" << part << ":" << orien << ":" << cassette << ":" << place << ":"
669  << partoffset << ":" << HGCalTypes::WaferTypeOffset[type] << ":" << i << ":"
670  << waferPart_.size();
671 #endif
672  wafer = waferPart_[i];
673  }
674  } else {
675  type = absType;
676  if (part == HGCalTypes::WaferFull) {
677  i = absType - 1;
678  wafer = passiveFull_[i];
679 #ifdef EDM_ML_DEBUG
680  edm::LogVerbatim("HGCalGeom") << " layertype:abstype:part:orien:cassette:offsets:ind " << layertype << ":"
681  << absType << ":" << part << ":" << orien << ":" << cassette << ":"
682  << ":" << partialTypes_ << ":" << orientationTypes_ << " passive " << i << ":"
683  << wafer;
684 #endif
685  } else {
686  int partoffset = (part >= HGCalTypes::WaferHDTop)
689  i = (part - partoffset) * facingTypes_ * orientationTypes_ +
690  (absType - 1) * facingTypes_ * orientationTypes_ * partialTypes_ + place - placeOffset_;
691 #ifdef EDM_ML_DEBUG
692  edm::LogVerbatim("HGCalGeom") << " layertype:abstype:part:orien:cassette:3Types:offset:ind " << layertype << ":"
693  << absType << ":" << part << ":" << orien << ":" << cassette << ":"
694  << partialTypes_ << ":" << facingTypes_ << ":" << orientationTypes_ << ":"
695  << partoffset << ":" << i << ":" << passivePart_.size();
696 #endif
697  wafer = passivePart_[i];
698  }
699  }
700  int copy = HGCalTypes::packTypeUV(type, u, v);
701 #ifdef EDM_ML_DEBUG
702  edm::LogVerbatim("HGCalGeom") << " DDHGCalMixRotatedCassette: Layer " << HGCalWaferIndex::waferLayer(waferIndex_[k])
703  << " Wafer " << wafer << " number " << copy << " type :part:orien:ind " << type << ":"
704  << part << ":" << orien << ":" << i << " layer:u:v " << (layer + firstLayer_) << ":"
705  << u << ":" << v;
706  if (iu > ium)
707  ium = iu;
708  if (iv > ivm)
709  ivm = iv;
710  kount++;
711  if (copies_.count(copy) == 0)
712  copies_.insert(copy);
713 #endif
714  DDTranslation tran(xpos, ypos, 0.0);
715  DDName name = DDName(DDSplit(wafer).first, DDSplit(wafer).second);
716  cpv.position(name, glog.ddname(), copy, tran, rot);
717 #ifdef EDM_ML_DEBUG
718  ++ntype[type];
719  edm::LogVerbatim("HGCalGeom") << " DDHGCalMixRotatedCassette: " << name << " number " << copy << " type "
720  << layertype << ":" << type << " positioned in " << glog.ddname() << " at " << tran
721  << " with no rotation";
722 #endif
723  }
724 #ifdef EDM_ML_DEBUG
725  edm::LogVerbatim("HGCalGeom") << "DDHGCalMixRotatedCassette: Maximum # of u " << ium << " # of v " << ivm << " and "
726  << kount << " wafers (" << ntype[0] << ":" << ntype[1] << ":" << ntype[2] << ") for "
727  << glog.ddname();
728 #endif
729 }
Log< level::Info, true > LogVerbatim
static int32_t cellPlacementIndex(int32_t iz, int32_t frontBack, int32_t orient)
Definition: HGCalCell.cc:239
static constexpr int32_t WaferPartLDOffset
Definition: HGCalTypes.h:57
std::unordered_set< int > copies_
void position(const DDLogicalPart &self, const DDLogicalPart &parent, const std::string &copyno, const DDTranslation &trans, const DDRotation &rot, const DDDivision *div=nullptr)
static constexpr int32_t WaferTypeOffset[3]
Definition: HGCalTypes.h:61
DDMaterial is used to define and access material information.
Definition: DDMaterial.h:45
int32_t waferU(const int32_t index)
int32_t waferLayer(const int32_t index)
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
std::pair< double, double > getShift(int layer, int zside, int cassette, bool scnt=false) const
DDName is used to identify DDD entities uniquely.
Definition: DDName.h:17
Log< level::Error, false > LogError
int32_t waferOrient(const int32_t property)
std::pair< double, double > shiftXY(int waferPosition, double waferSize) const
A DDSolid represents the shape of a part.
Definition: DDSolid.h:39
static std::string to_string(const XMLCh *ch)
Represents a uniquely identifyable rotation matrix.
Definition: DDTransform.h:57
U second(std::pair< T, U > const &p)
int32_t waferCassette(const int32_t property)
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:35
T sqrt(T t)
Definition: SSEVec.h:23
std::vector< std::string > namesTop_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
A DDLogicalPart aggregates information concerning material, solid and sensitveness ...
Definition: DDLogicalPart.h:93
std::vector< std::string > waferFull_
static DDSolid tubs(const DDName &name, double zhalf, double rIn, double rOut, double startPhi, double deltaPhi)
Definition: DDSolid.cc:667
int32_t waferThick(const int32_t property)
ii
Definition: cuy.py:589
const N & name() const
Definition: DDBase.h:59
std::tuple< int32_t, int32_t, int32_t > tileUnpack(int32_t index)
static constexpr int32_t WaferPartHDOffset
Definition: HGCalTypes.h:58
std::vector< std::string > passivePart_
part
Definition: HCALResponse.h:20
const N & ddname() const
Definition: DDBase.h:61
std::vector< std::string > materialTop_
std::vector< std::string > passiveFull_
int32_t waferPartial(const int32_t property)
int32_t waferV(const int32_t index)
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
std::pair< std::string, std::string > DDSplit(const std::string &n)
split into (name,namespace), separator = &#39;:&#39;
Definition: DDSplit.cc:3
std::vector< std::string > waferPart_
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DDTranslation
Definition: DDTranslation.h:7
static int cassetteType(int det, int zside, int cassette)
T angle(T x1, T y1, T z1, T x2, T y2, T z2)
Definition: angle.h:11
static constexpr int32_t layerFrontBack(int32_t layerOrient)
Definition: HGCalTypes.h:137

Member Data Documentation

◆ absorbMode_

int DDHGCalMixRotatedCassette::absorbMode_
private

Definition at line 73 of file DDHGCalMixRotatedCassette.cc.

◆ alpha_

double DDHGCalMixRotatedCassette::alpha_
private

Definition at line 117 of file DDHGCalMixRotatedCassette.cc.

◆ cassette_

HGCalCassette DDHGCalMixRotatedCassette::cassette_
private

Definition at line 60 of file DDHGCalMixRotatedCassette.cc.

◆ cassettes_

int DDHGCalMixRotatedCassette::cassettes_
private

Definition at line 79 of file DDHGCalMixRotatedCassette.cc.

◆ cassetteShift_

std::vector<double> DDHGCalMixRotatedCassette::cassetteShift_
private

Definition at line 109 of file DDHGCalMixRotatedCassette.cc.

◆ copies_

std::unordered_set<int> DDHGCalMixRotatedCassette::copies_
private

Definition at line 116 of file DDHGCalMixRotatedCassette.cc.

◆ copyNumber_

std::vector<int> DDHGCalMixRotatedCassette::copyNumber_
private

Definition at line 93 of file DDHGCalMixRotatedCassette.cc.

◆ copyNumberCoverTop_

int DDHGCalMixRotatedCassette::copyNumberCoverTop_
private

Definition at line 104 of file DDHGCalMixRotatedCassette.cc.

◆ copyNumberTop_

std::vector<int> DDHGCalMixRotatedCassette::copyNumberTop_
private

Definition at line 102 of file DDHGCalMixRotatedCassette.cc.

◆ cosAlpha_

double DDHGCalMixRotatedCassette::cosAlpha_
private

Definition at line 117 of file DDHGCalMixRotatedCassette.cc.

◆ coverTypeTop_

int DDHGCalMixRotatedCassette::coverTypeTop_
private

Definition at line 103 of file DDHGCalMixRotatedCassette.cc.

◆ facingTypes_

int DDHGCalMixRotatedCassette::facingTypes_
private

Definition at line 67 of file DDHGCalMixRotatedCassette.cc.

◆ firstLayer_

int DDHGCalMixRotatedCassette::firstLayer_
private

Definition at line 72 of file DDHGCalMixRotatedCassette.cc.

◆ geomTools_

HGCalGeomTools DDHGCalMixRotatedCassette::geomTools_
private

Definition at line 59 of file DDHGCalMixRotatedCassette.cc.

◆ layerOrient_

std::vector<int> DDHGCalMixRotatedCassette::layerOrient_
private

Definition at line 105 of file DDHGCalMixRotatedCassette.cc.

◆ layers_

std::vector<int> DDHGCalMixRotatedCassette::layers_
private

Definition at line 94 of file DDHGCalMixRotatedCassette.cc.

◆ layerSense_

std::vector<int> DDHGCalMixRotatedCassette::layerSense_
private

Definition at line 97 of file DDHGCalMixRotatedCassette.cc.

◆ layerThick_

std::vector<double> DDHGCalMixRotatedCassette::layerThick_
private

Definition at line 95 of file DDHGCalMixRotatedCassette.cc.

◆ layerThickTop_

std::vector<double> DDHGCalMixRotatedCassette::layerThickTop_
private

Definition at line 100 of file DDHGCalMixRotatedCassette.cc.

◆ layerType_

std::vector<int> DDHGCalMixRotatedCassette::layerType_
private

Definition at line 96 of file DDHGCalMixRotatedCassette.cc.

◆ layerTypeTop_

std::vector<int> DDHGCalMixRotatedCassette::layerTypeTop_
private

Definition at line 101 of file DDHGCalMixRotatedCassette.cc.

◆ materials_

std::vector<std::string> DDHGCalMixRotatedCassette::materials_
private

Definition at line 90 of file DDHGCalMixRotatedCassette.cc.

◆ materialTop_

std::vector<std::string> DDHGCalMixRotatedCassette::materialTop_
private

Definition at line 98 of file DDHGCalMixRotatedCassette.cc.

◆ names_

std::vector<std::string> DDHGCalMixRotatedCassette::names_
private

Definition at line 91 of file DDHGCalMixRotatedCassette.cc.

◆ nameSpace_

std::string DDHGCalMixRotatedCassette::nameSpace_
private

Definition at line 115 of file DDHGCalMixRotatedCassette.cc.

◆ namesTop_

std::vector<std::string> DDHGCalMixRotatedCassette::namesTop_
private

Definition at line 99 of file DDHGCalMixRotatedCassette.cc.

◆ orientationTypes_

int DDHGCalMixRotatedCassette::orientationTypes_
private

Definition at line 68 of file DDHGCalMixRotatedCassette.cc.

◆ partialTypes_

int DDHGCalMixRotatedCassette::partialTypes_
private

Definition at line 69 of file DDHGCalMixRotatedCassette.cc.

◆ passiveFull_

std::vector<std::string> DDHGCalMixRotatedCassette::passiveFull_
private

Definition at line 88 of file DDHGCalMixRotatedCassette.cc.

◆ passivePart_

std::vector<std::string> DDHGCalMixRotatedCassette::passivePart_
private

Definition at line 89 of file DDHGCalMixRotatedCassette.cc.

◆ passiveTypes_

int DDHGCalMixRotatedCassette::passiveTypes_
private

Definition at line 66 of file DDHGCalMixRotatedCassette.cc.

◆ phiBinsScint_

int DDHGCalMixRotatedCassette::phiBinsScint_
private

Definition at line 71 of file DDHGCalMixRotatedCassette.cc.

◆ placeOffset_

int DDHGCalMixRotatedCassette::placeOffset_
private

Definition at line 70 of file DDHGCalMixRotatedCassette.cc.

◆ rMaxFront_

std::vector<double> DDHGCalMixRotatedCassette::rMaxFront_
private

Definition at line 85 of file DDHGCalMixRotatedCassette.cc.

◆ rMinFront_

std::vector<double> DDHGCalMixRotatedCassette::rMinFront_
private

Definition at line 82 of file DDHGCalMixRotatedCassette.cc.

◆ sectors_

int DDHGCalMixRotatedCassette::sectors_
private

Definition at line 78 of file DDHGCalMixRotatedCassette.cc.

◆ sensitiveMode_

int DDHGCalMixRotatedCassette::sensitiveMode_
private

Definition at line 74 of file DDHGCalMixRotatedCassette.cc.

◆ slopeB_

std::vector<double> DDHGCalMixRotatedCassette::slopeB_
private

Definition at line 80 of file DDHGCalMixRotatedCassette.cc.

◆ slopeT_

std::vector<double> DDHGCalMixRotatedCassette::slopeT_
private

Definition at line 83 of file DDHGCalMixRotatedCassette.cc.

◆ thick_

std::vector<double> DDHGCalMixRotatedCassette::thick_
private

Definition at line 92 of file DDHGCalMixRotatedCassette.cc.

◆ tileIndex_

std::vector<int> DDHGCalMixRotatedCassette::tileIndex_
private

Definition at line 112 of file DDHGCalMixRotatedCassette.cc.

◆ tileLayerStart_

std::vector<int> DDHGCalMixRotatedCassette::tileLayerStart_
private

Definition at line 114 of file DDHGCalMixRotatedCassette.cc.

◆ tilePhis_

std::vector<int> DDHGCalMixRotatedCassette::tilePhis_
private

Definition at line 113 of file DDHGCalMixRotatedCassette.cc.

◆ tileRMax_

std::vector<double> DDHGCalMixRotatedCassette::tileRMax_
private

Definition at line 111 of file DDHGCalMixRotatedCassette.cc.

◆ tileRMin_

std::vector<double> DDHGCalMixRotatedCassette::tileRMin_
private

Definition at line 110 of file DDHGCalMixRotatedCassette.cc.

◆ tol1_

constexpr double DDHGCalMixRotatedCassette::tol1_ = 0.01
staticprivate

Definition at line 62 of file DDHGCalMixRotatedCassette.cc.

◆ tol2_

constexpr double DDHGCalMixRotatedCassette::tol2_ = 0.00001
staticprivate

Definition at line 63 of file DDHGCalMixRotatedCassette.cc.

◆ waferFull_

std::vector<std::string> DDHGCalMixRotatedCassette::waferFull_
private

Definition at line 86 of file DDHGCalMixRotatedCassette.cc.

◆ waferIndex_

std::vector<int> DDHGCalMixRotatedCassette::waferIndex_
private

Definition at line 106 of file DDHGCalMixRotatedCassette.cc.

◆ waferLayerStart_

std::vector<int> DDHGCalMixRotatedCassette::waferLayerStart_
private

Definition at line 108 of file DDHGCalMixRotatedCassette.cc.

◆ waferPart_

std::vector<std::string> DDHGCalMixRotatedCassette::waferPart_
private

Definition at line 87 of file DDHGCalMixRotatedCassette.cc.

◆ waferProperty_

std::vector<int> DDHGCalMixRotatedCassette::waferProperty_
private

Definition at line 107 of file DDHGCalMixRotatedCassette.cc.

◆ waferSepar_

double DDHGCalMixRotatedCassette::waferSepar_
private

Definition at line 77 of file DDHGCalMixRotatedCassette.cc.

◆ waferSize_

double DDHGCalMixRotatedCassette::waferSize_
private

Definition at line 76 of file DDHGCalMixRotatedCassette.cc.

◆ waferTypes_

int DDHGCalMixRotatedCassette::waferTypes_
private

Definition at line 65 of file DDHGCalMixRotatedCassette.cc.

◆ zFrontB_

std::vector<double> DDHGCalMixRotatedCassette::zFrontB_
private

Definition at line 81 of file DDHGCalMixRotatedCassette.cc.

◆ zFrontT_

std::vector<double> DDHGCalMixRotatedCassette::zFrontT_
private

Definition at line 84 of file DDHGCalMixRotatedCassette.cc.

◆ zMinBlock_

double DDHGCalMixRotatedCassette::zMinBlock_
private

Definition at line 75 of file DDHGCalMixRotatedCassette.cc.