CMS 3D CMS Logo

List of all members | Public Member Functions | Public Attributes
HGCalSiliconModule Struct Reference

Public Member Functions

 HGCalSiliconModule ()
 
 HGCalSiliconModule (cms::DDParsingContext &ctxt, xml_h e)
 
void positionSensitive (cms::DDParsingContext &ctxt, xml_h e, const dd4hep::Volume &glog, int layer)
 

Public Attributes

int absorbMode_
 
double alpha_
 
std::unordered_set< int > copies_
 
std::vector< int > copyNumber_
 
double cosAlpha_
 
int facingTypes_
 
int firstLayer_
 
HGCalGeomTools geomTools_
 
std::vector< int > layers_
 
std::vector< int > layerSense_
 
std::vector< double > layerThick_
 
std::vector< int > layerType_
 
std::vector< int > layerTypes_
 
std::vector< std::string > materials_
 
std::vector< std::string > names_
 
int orientationTypes_
 
int partialTypes_
 
std::vector< double > rMaxFront_
 
std::vector< double > rMinFront_
 
std::string rotstr_
 
int sectors_
 
int sensitiveMode_
 
std::vector< double > slopeB_
 
std::vector< double > slopeT_
 
std::vector< double > thick_
 
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_
 

Detailed Description

Definition at line 28 of file DDHGCalSiliconModule.cc.

Constructor & Destructor Documentation

◆ HGCalSiliconModule() [1/2]

HGCalSiliconModule::HGCalSiliconModule ( )
inline

Definition at line 29 of file DDHGCalSiliconModule.cc.

References Exception.

29 { throw cms::Exception("HGCalGeom") << "Wrong initialization to HGCalSiliconModule"; }

◆ HGCalSiliconModule() [2/2]

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

Definition at line 30 of file DDHGCalSiliconModule.cc.

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

30  {
31  cms::DDNamespace ns(ctxt, e, true);
33 
34 #ifdef EDM_ML_DEBUG
35  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: Creating an instance";
36 #endif
37  static constexpr double tol1 = 0.01 * dd4hep::mm;
38  static constexpr double tol2 = 0.00001 * dd4hep::mm;
39 
40  dd4hep::Volume mother = ns.volume(args.parentName());
41  waferTypes_ = args.value<int>("WaferTypes");
42  facingTypes_ = args.value<int>("FacingTypes");
43  partialTypes_ = args.value<int>("PartialTypes");
44  orientationTypes_ = args.value<int>("OrientationTypes");
45 #ifdef EDM_ML_DEBUG
46  edm::LogVerbatim("HGCalGeom") << "Number of types of wafers: " << waferTypes_ << " facings: " << facingTypes_
47  << " partials: " << partialTypes_ << " Orientations: " << orientationTypes_;
48 #endif
49  firstLayer_ = args.value<int>("FirstLayer");
50  absorbMode_ = args.value<int>("AbsorberMode");
51  sensitiveMode_ = args.value<int>("SensitiveMode");
52 #ifdef EDM_ML_DEBUG
53  edm::LogVerbatim("HGCalGeom") << "First Layer " << firstLayer_ << " and "
54  << "Absober:Sensitive mode " << absorbMode_ << ":" << sensitiveMode_;
55 #endif
56  zMinBlock_ = args.value<double>("zMinBlock");
57  waferSize_ = args.value<double>("waferSize");
58  waferSepar_ = args.value<double>("SensorSeparation");
59  sectors_ = args.value<int>("Sectors");
60  alpha_ = (1._pi) / sectors_;
61  cosAlpha_ = cos(alpha_);
62  rotstr_ = args.value<std::string>("LayerRotation");
63 #ifdef EDM_ML_DEBUG
64  edm::LogVerbatim("HGCalGeom") << "zStart " << cms::convert2mm(zMinBlock_) << " wafer width "
65  << cms::convert2mm(waferSize_) << " separations " << cms::convert2mm(waferSepar_)
66  << " sectors " << sectors_ << ":" << convertRadToDeg(alpha_) << ":" << cosAlpha_
67  << " rotation matrix " << rotstr_;
68 #endif
69  waferFull_ = args.value<std::vector<std::string>>("WaferNamesFull");
70  waferPart_ = args.value<std::vector<std::string>>("WaferNamesPartial");
71 #ifdef EDM_ML_DEBUG
72  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << waferFull_.size() << " full and " << waferPart_.size()
73  << " partial modules\nDDHGCalSiliconModule:Full Modules:";
74  unsigned int i1max = static_cast<unsigned int>(waferFull_.size());
75  for (unsigned int i1 = 0; i1 < i1max; i1 += 2) {
76  std::ostringstream st1;
77  unsigned int i2 = std::min((i1 + 2), i1max);
78  for (unsigned int i = i1; i < i2; ++i)
79  st1 << " [" << i << "] " << waferFull_[i];
80  edm::LogVerbatim("HGCalGeom") << st1.str() << std::endl;
81  }
82  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: Partial Modules:";
83  i1max = static_cast<unsigned int>(waferPart_.size());
84  for (unsigned int i1 = 0; i1 < i1max; i1 += 2) {
85  std::ostringstream st1;
86  unsigned int i2 = std::min((i1 + 2), i1max);
87  for (unsigned int i = i1; i < i2; ++i)
88  st1 << " [" << i << "] " << waferPart_[i];
89  edm::LogVerbatim("HGCalGeom") << st1.str() << std::endl;
90  }
91 #endif
92  materials_ = args.value<std::vector<std::string>>("MaterialNames");
93  names_ = args.value<std::vector<std::string>>("VolumeNames");
94  thick_ = args.value<std::vector<double>>("Thickness");
95  copyNumber_.resize(materials_.size(), 1);
96 #ifdef EDM_ML_DEBUG
97  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << materials_.size() << " types of volumes";
98  for (unsigned int i = 0; i < names_.size(); ++i)
99  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << names_[i] << " of thickness "
100  << cms::convert2mm(thick_[i]) << " filled with " << materials_[i]
101  << " first copy number " << copyNumber_[i];
102 #endif
103  layers_ = args.value<std::vector<int>>("Layers");
104  layerThick_ = args.value<std::vector<double>>("LayerThick");
105 #ifdef EDM_ML_DEBUG
106  edm::LogVerbatim("HGCalGeom") << "There are " << layers_.size() << " blocks";
107  for (unsigned int i = 0; i < layers_.size(); ++i)
108  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] of thickness " << cms::convert2mm(layerThick_[i])
109  << " with " << layers_[i] << " layers";
110 #endif
111  layerType_ = args.value<std::vector<int>>("LayerType");
112  layerSense_ = args.value<std::vector<int>>("LayerSense");
113  layerTypes_ = args.value<std::vector<int>>("LayerTypes");
114 #ifdef EDM_ML_DEBUG
115  for (unsigned int i = 0; i < layerTypes_.size(); ++i)
116  edm::LogVerbatim("HGCalGeom") << "LayerTypes [" << i << "] " << layerTypes_[i];
117 #endif
118  if (firstLayer_ > 0) {
119  for (unsigned int i = 0; i < layerType_.size(); ++i) {
120  if (layerSense_[i] > 0) {
121  int ii = layerType_[i];
122  copyNumber_[ii] = (layerSense_[i] == 1) ? firstLayer_ : (firstLayer_ + 1);
123 #ifdef EDM_ML_DEBUG
124  edm::LogVerbatim("HGCalGeom") << "First copy number for layer type " << i << ":" << ii << " with "
125  << materials_[ii] << " changed to " << copyNumber_[ii];
126 #endif
127  }
128  }
129  } else {
130  firstLayer_ = 1;
131  }
132 #ifdef EDM_ML_DEBUG
133  edm::LogVerbatim("HGCalGeom") << "There are " << layerType_.size() << " layers";
134  for (unsigned int i = 0; i < layerType_.size(); ++i)
135  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerType_[i] << " sensitive class "
136  << layerSense_[i];
137 #endif
138  slopeB_ = args.value<std::vector<double>>("SlopeBottom");
139  zFrontB_ = args.value<std::vector<double>>("ZFrontBottom");
140  rMinFront_ = args.value<std::vector<double>>("RMinFront");
141  slopeT_ = args.value<std::vector<double>>("SlopeTop");
142  zFrontT_ = args.value<std::vector<double>>("ZFrontTop");
143  rMaxFront_ = args.value<std::vector<double>>("RMaxFront");
144 #ifdef EDM_ML_DEBUG
145  for (unsigned int i = 0; i < slopeB_.size(); ++i)
146  edm::LogVerbatim("HGCalGeom") << "Bottom Block [" << i << "] Zmin " << cms::convert2mm(zFrontB_[i]) << " Rmin "
147  << cms::convert2mm(rMinFront_[i]) << " Slope " << slopeB_[i];
148  for (unsigned int i = 0; i < slopeT_.size(); ++i)
149  edm::LogVerbatim("HGCalGeom") << "Top Block [" << i << "] Zmin " << cms::convert2mm(zFrontT_[i]) << " Rmax "
150  << cms::convert2mm(rMaxFront_[i]) << " Slope " << slopeT_[i];
151 #endif
152  waferIndex_ = args.value<std::vector<int>>("WaferIndex");
153  waferProperty_ = args.value<std::vector<int>>("WaferProperties");
154  waferLayerStart_ = args.value<std::vector<int>>("WaferLayerStart");
155 #ifdef EDM_ML_DEBUG
156  edm::LogVerbatim("HGCalGeom") << "waferProperties with " << waferIndex_.size() << " entries in "
157  << waferLayerStart_.size() << " layers";
158  for (unsigned int k = 0; k < waferLayerStart_.size(); ++k)
159  edm::LogVerbatim("HGCalGeom") << "LayerStart[" << k << "] " << waferLayerStart_[k];
160  for (unsigned int k = 0; k < waferIndex_.size(); ++k)
161  edm::LogVerbatim("HGCalGeom") << "Wafer[" << k << "] " << waferIndex_[k] << " ("
164  << HGCalWaferIndex::waferV(waferIndex_[k]) << ") : ("
168  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: NameSpace " << ns.name();
169 #endif
170 
171 #ifdef EDM_ML_DEBUG
172  edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalSiliconModule...";
173  copies_.clear();
174 #endif
175 
176  double zi(zMinBlock_);
177  int laymin(0);
178  for (unsigned int i = 0; i < layers_.size(); ++i) {
179  double zo = zi + layerThick_[i];
180  double routF = HGCalGeomTools::radius(zi, zFrontT_, rMaxFront_, slopeT_);
181  int laymax = laymin + layers_[i];
182  double zz = zi;
183  double thickTot(0);
184  for (int ly = laymin; ly < laymax; ++ly) {
185  int ii = layerType_[ly];
186  int copy = copyNumber_[ii];
187  double hthick = 0.5 * thick_[ii];
188  double rinB = HGCalGeomTools::radius(zo - tol1, zFrontB_, rMinFront_, slopeB_);
189  zz += hthick;
190  thickTot += thick_[ii];
191 
193 #ifdef EDM_ML_DEBUG
194  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: Layer " << ly << ":" << ii << " Front "
195  << cms::convert2mm(zi) << ", " << cms::convert2mm(routF) << " Back "
196  << cms::convert2mm(zo) << ", " << cms::convert2mm(rinB)
197  << " superlayer thickness " << cms::convert2mm(layerThick_[i]);
198 #endif
199 
200  dd4hep::Material matter = ns.material(materials_[ii]);
201  dd4hep::Volume glog;
202 
203  if (layerSense_[ly] < 1) {
204  std::vector<double> pgonZ, pgonRin, pgonRout;
205  double rmax = routF * cosAlpha_ - tol1;
206  HGCalGeomTools::radius(zz - hthick,
207  zz + hthick,
208  zFrontB_,
209  rMinFront_,
210  slopeB_,
211  zFrontT_,
212  rMaxFront_,
213  slopeT_,
214  -layerSense_[ly],
215  pgonZ,
216  pgonRin,
217  pgonRout);
218  for (unsigned int isec = 0; isec < pgonZ.size(); ++isec) {
219  pgonZ[isec] -= zz;
220  if (layerSense_[ly] == 0 || absorbMode_ == 0)
221  pgonRout[isec] = rmax;
222  else
223  pgonRout[isec] = pgonRout[isec] * cosAlpha_ - tol1;
224  }
225  dd4hep::Solid solid = dd4hep::Polyhedra(sectors_, -alpha_, 2._pi, pgonZ, pgonRin, pgonRout);
226  ns.addSolidNS(ns.prepend(name), solid);
227  glog = dd4hep::Volume(solid.name(), solid, matter);
228  ns.addVolumeNS(glog);
229 #ifdef EDM_ML_DEBUG
230  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << solid.name() << " polyhedra of " << sectors_
231  << " sectors covering " << convertRadToDeg(-alpha_) << ":"
232  << convertRadToDeg(-alpha_ + 2._pi) << " with " << pgonZ.size()
233  << " sections and filled with " << matter.name();
234  for (unsigned int k = 0; k < pgonZ.size(); ++k)
235  edm::LogVerbatim("HGCalGeom") << "[" << k << "] z " << cms::convert2mm(pgonZ[k]) << " R "
236  << cms::convert2mm(pgonRin[k]) << ":" << cms::convert2mm(pgonRout[k]);
237 #endif
238  } else {
239  double rins =
240  (sensitiveMode_ < 1) ? rinB : HGCalGeomTools::radius(zz + hthick - tol1, zFrontB_, rMinFront_, slopeB_);
241  double routs =
242  (sensitiveMode_ < 1) ? routF : HGCalGeomTools::radius(zz - hthick, zFrontT_, rMaxFront_, slopeT_);
243  dd4hep::Solid solid = dd4hep::Tube(rins, routs, hthick, 0.0, 2._pi);
244  ns.addSolidNS(ns.prepend(name), solid);
245  glog = dd4hep::Volume(solid.name(), solid, matter);
246  ns.addVolumeNS(glog);
247 
248 #ifdef EDM_ML_DEBUG
249  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << solid.name() << " Tubs made of " << matter.name()
250  << " of dimensions " << cms::convert2mm(rinB) << ":" << cms::convert2mm(rins)
251  << ", " << cms::convert2mm(routF) << ":" << cms::convert2mm(routs) << ", "
252  << cms::convert2mm(hthick) << ", 0.0, 360.0 and position " << glog.name()
253  << " number " << copy << ":" << layerTypes_[copy - firstLayer_];
254 #endif
255  positionSensitive(ctxt, e, glog, (copy - firstLayer_));
256  }
257 
258  dd4hep::Position r1(0, 0, zz);
259  dd4hep::Rotation3D rot;
260 #ifdef EDM_ML_DEBUG
261  std::string rotName("Null");
262 #endif
264  rot = ns.rotation(rotstr_);
265 #ifdef EDM_ML_DEBUG
266  rotName = rotstr_;
267 #endif
268  }
269  mother.placeVolume(glog, copy, dd4hep::Transform3D(rot, r1));
270  int inc = ((layerSense_[ly] > 0) && (facingTypes_ > 1)) ? 2 : 1;
271  copyNumber_[ii] = copy + inc;
272 #ifdef EDM_ML_DEBUG
273  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << glog.name() << " number " << copy
274  << " positioned in " << mother.name() << " at (0,0," << cms::convert2mm(zz)
275  << ") with " << rotName << " rotation";
276 #endif
277  zz += hthick;
278  } // End of loop over layers in a block
279  zi = zo;
280  laymin = laymax;
281  // Make consistency check of all the partitions of the block
282  if (std::abs(thickTot - layerThick_[i]) >= tol2) {
283  if (thickTot > layerThick_[i]) {
284  edm::LogError("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
285  << " is smaller than " << cms::convert2mm(thickTot)
286  << ": thickness of all its components **** ERROR ****";
287  } else {
288  edm::LogWarning("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
289  << " does not match with " << cms::convert2mm(thickTot) << " of the components";
290  }
291  }
292  } // End of loop over blocks
293 
294 #ifdef EDM_ML_DEBUG
295  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << copies_.size() << " different wafer copy numbers";
296  int k(0);
297  for (std::unordered_set<int>::const_iterator itr = copies_.begin(); itr != copies_.end(); ++itr, ++k) {
298  edm::LogVerbatim("HGCalGeom") << "Copy [" << k << "] : " << (*itr);
299  }
300  copies_.clear();
301  edm::LogVerbatim("HGCalGeom") << "<<== End of DDHGCalSiliconModule construction...";
302 #endif
303  }
Log< level::Info, true > LogVerbatim
static constexpr int32_t WaferCenteredRotated
Definition: HGCalTypes.h:71
static void radius(double zf, double zb, std::vector< double > const &zFront1, std::vector< double > const &rFront1, std::vector< double > const &slope1, std::vector< double > const &zFront2, std::vector< double > const &rFront2, std::vector< double > const &slope2, int flag, std::vector< double > &zz, std::vector< double > &rin, std::vector< double > &rout)
std::unordered_set< int > copies_
std::vector< double > slopeT_
std::vector< double > slopeB_
int32_t waferU(const int32_t index)
int32_t waferLayer(const int32_t index)
std::vector< int > copyNumber_
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
constexpr NumType convert2mm(NumType length)
Definition: DDutils.h:7
std::string to_string(const V &value)
Definition: OMSAccess.h:71
std::vector< int > layerSense_
Log< level::Error, false > LogError
int32_t waferOrient(const int32_t property)
std::vector< double > zFrontT_
static const std::string & rotName(const T &rot, const cms::DDParsingContext &context)
std::vector< int > layerTypes_
std::vector< double > zFrontB_
std::vector< double > thick_
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
std::vector< std::string > names_
int32_t waferThick(const int32_t property)
std::vector< int > waferProperty_
dd4hep::Volume Volume
ii
Definition: cuy.py:589
std::vector< std::string > waferPart_
std::vector< double > rMaxFront_
std::vector< int > waferLayerStart_
std::vector< std::string > materials_
std::vector< int > layerType_
void positionSensitive(cms::DDParsingContext &ctxt, xml_h e, const dd4hep::Volume &glog, int layer)
int32_t waferPartial(const int32_t property)
int32_t waferV(const int32_t index)
std::vector< int > layers_
Log< level::Warning, false > LogWarning
std::vector< std::string > waferFull_
std::vector< double > rMinFront_
std::vector< int > waferIndex_
std::vector< double > layerThick_

Member Function Documentation

◆ positionSensitive()

void HGCalSiliconModule::positionSensitive ( cms::DDParsingContext ctxt,
xml_h  e,
const dd4hep::Volume &  glog,
int  layer 
)
inline

Definition at line 305 of file DDHGCalSiliconModule.cc.

References funct::abs(), cms::convert2mm(), filterCSVwithJSON::copy, HGCalTypes::CornerCenteredLambda, HGCalTypes::CornerCenteredY, HGCalTypes::CornerCenterYm, HGCalTypes::CornerCenterYp, PVValHelper::dy, MillePedeFileConverter_cfg::e, mps_fire::i, gpuVertexFinder::iv, dqmdumpme::k, phase1PixelTopology::layer, EgHLTOffHistBins_cfi::nr, HGCalTypes::packTypeUV(), PixelTestBeamValidation_cfi::Position, dttmaxenums::R, alignCSCRings::r, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, findQualityFiles::v, cms::DDNamespace::volume(), HGCalTypes::WaferCenter, HGCalTypes::WaferCenteredBack, HGCalTypes::WaferFull, HGCalWaferIndex::waferLayer(), HGCalProperty::waferOrient(), HGCalProperty::waferPartial(), HGCalProperty::waferThick(), HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

305  {
306  cms::DDNamespace ns(ctxt, e, true);
307  static const double sqrt3 = std::sqrt(3.0);
308  int layercenter = (layerTypes_[layer] == HGCalTypes::CornerCenteredLambda)
312  int layertype = (layerTypes_[layer] == HGCalTypes::WaferCenteredBack) ? 1 : 0;
313  int firstWafer = waferLayerStart_[layer];
314  int lastWafer = ((layer + 1 < static_cast<int>(waferLayerStart_.size())) ? waferLayerStart_[layer + 1]
315  : static_cast<int>(waferIndex_.size()));
316  double r = 0.5 * (waferSize_ + waferSepar_);
317  double R = 2.0 * r / sqrt3;
318  double dy = 0.75 * R;
319  const auto& xyoff = geomTools_.shiftXY(layercenter, (waferSize_ + waferSepar_));
320 #ifdef EDM_ML_DEBUG
321  int ium(0), ivm(0), kount(0);
322  std::vector<int> ntype(3, 0);
323  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << glog.name() << " r " << cms::convert2mm(r) << " R "
324  << cms::convert2mm(R) << " dy " << cms::convert2mm(dy) << " Shift "
325  << cms::convert2mm(xyoff.first) << ":" << cms::convert2mm(xyoff.second)
326  << " WaferSize " << cms::convert2mm((waferSize_ + waferSepar_)) << " index "
327  << firstWafer << ":" << (lastWafer - 1);
328 #endif
329  for (int k = firstWafer; k < lastWafer; ++k) {
332 #ifdef EDM_ML_DEBUG
333  int iu = std::abs(u);
334  int iv = std::abs(v);
335 #endif
336  int nr = 2 * v;
337  int nc = -2 * u + v;
338  double xpos = xyoff.first + nc * r;
339  double ypos = xyoff.second + nr * dy;
343  std::string wafer;
344  int i(999);
345  if (part == HGCalTypes::WaferFull) {
346  i = layertype * waferTypes_ + type;
347  wafer = waferFull_[i];
348  } else {
350  type * orientationTypes_ + orien;
351 #ifdef EDM_ML_DEBUG
352  edm::LogVerbatim("HGCalGeom") << " layertype:type:part:orien:ind " << layertype << ":" << type << ":" << part
353  << ":" << orien << ":" << i << ":" << waferPart_.size();
354 #endif
355  wafer = waferPart_[i];
356  }
357  int copy = HGCalTypes::packTypeUV(type, u, v);
358 #ifdef EDM_ML_DEBUG
359  edm::LogVerbatim("HGCalGeom") << " DDHGCalSiliconModule: Layer" << HGCalWaferIndex::waferLayer(waferIndex_[k])
360  << " Wafer " << wafer << " number " << copy << " type:part:orien:ind " << type
361  << ":" << part << ":" << orien << ":" << i << " layer:u:v:indx "
362  << (layer + firstLayer_) << ":" << u << ":" << v;
363  if (iu > ium)
364  ium = iu;
365  if (iv > ivm)
366  ivm = iv;
367  kount++;
368  if (copies_.count(copy) == 0)
369  copies_.insert(copy);
370 #endif
371  dd4hep::Position tran(xpos, ypos, 0.0);
372  glog.placeVolume(ns.volume(wafer), copy, tran);
373 #ifdef EDM_ML_DEBUG
374  ++ntype[type];
375  edm::LogVerbatim("HGCalGeom") << " DDHGCalSiliconModule: " << wafer << " number " << copy << " type " << layertype
376  << ":" << type << " positioned in " << glog.name() << " at ("
377  << cms::convert2mm(xpos) << "," << cms::convert2mm(ypos) << ",0) with no rotation";
378 #endif
379  }
380 
381 #ifdef EDM_ML_DEBUG
382  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: Maximum # of u " << ium << " # of v " << ivm << " and "
383  << kount << " wafers (" << ntype[0] << ":" << ntype[1] << ":" << ntype[2] << ") for "
384  << glog.name();
385 #endif
386  }
Log< level::Info, true > LogVerbatim
static constexpr int32_t WaferCenter
Definition: HGCalTypes.h:20
std::unordered_set< int > copies_
int32_t *__restrict__ iv
int32_t waferU(const int32_t index)
int32_t waferLayer(const int32_t index)
static constexpr int32_t CornerCenterYm
Definition: HGCalTypes.h:22
constexpr NumType convert2mm(NumType length)
Definition: DDutils.h:7
int32_t waferOrient(const int32_t property)
std::pair< double, double > shiftXY(int waferPosition, double waferSize) const
constexpr std::array< uint8_t, layerIndexSize > layer
std::vector< int > layerTypes_
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:34
T sqrt(T t)
Definition: SSEVec.h:19
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
int32_t waferThick(const int32_t property)
std::vector< int > waferProperty_
static constexpr int32_t CornerCenterYp
Definition: HGCalTypes.h:21
std::vector< std::string > waferPart_
static constexpr int32_t WaferCenteredBack
Definition: HGCalTypes.h:68
std::vector< int > waferLayerStart_
part
Definition: HCALResponse.h:20
static constexpr int32_t CornerCenteredLambda
Definition: HGCalTypes.h:70
int32_t waferPartial(const int32_t property)
static constexpr int32_t CornerCenteredY
Definition: HGCalTypes.h:69
int32_t waferV(const int32_t index)
static int32_t packTypeUV(int type, int u, int v)
Definition: HGCalTypes.cc:3
std::vector< std::string > waferFull_
std::vector< int > waferIndex_

Member Data Documentation

◆ absorbMode_

int HGCalSiliconModule::absorbMode_

Definition at line 396 of file DDHGCalSiliconModule.cc.

◆ alpha_

double HGCalSiliconModule::alpha_

Definition at line 424 of file DDHGCalSiliconModule.cc.

◆ copies_

std::unordered_set<int> HGCalSiliconModule::copies_

Definition at line 423 of file DDHGCalSiliconModule.cc.

◆ copyNumber_

std::vector<int> HGCalSiliconModule::copyNumber_

Definition at line 408 of file DDHGCalSiliconModule.cc.

◆ cosAlpha_

double HGCalSiliconModule::cosAlpha_

Definition at line 424 of file DDHGCalSiliconModule.cc.

◆ facingTypes_

int HGCalSiliconModule::facingTypes_

Definition at line 392 of file DDHGCalSiliconModule.cc.

◆ firstLayer_

int HGCalSiliconModule::firstLayer_

Definition at line 395 of file DDHGCalSiliconModule.cc.

◆ geomTools_

HGCalGeomTools HGCalSiliconModule::geomTools_

Definition at line 389 of file DDHGCalSiliconModule.cc.

◆ layers_

std::vector<int> HGCalSiliconModule::layers_

Definition at line 409 of file DDHGCalSiliconModule.cc.

◆ layerSense_

std::vector<int> HGCalSiliconModule::layerSense_

Definition at line 412 of file DDHGCalSiliconModule.cc.

◆ layerThick_

std::vector<double> HGCalSiliconModule::layerThick_

Definition at line 410 of file DDHGCalSiliconModule.cc.

◆ layerType_

std::vector<int> HGCalSiliconModule::layerType_

Definition at line 411 of file DDHGCalSiliconModule.cc.

◆ layerTypes_

std::vector<int> HGCalSiliconModule::layerTypes_

Definition at line 419 of file DDHGCalSiliconModule.cc.

◆ materials_

std::vector<std::string> HGCalSiliconModule::materials_

Definition at line 405 of file DDHGCalSiliconModule.cc.

◆ names_

std::vector<std::string> HGCalSiliconModule::names_

Definition at line 406 of file DDHGCalSiliconModule.cc.

◆ orientationTypes_

int HGCalSiliconModule::orientationTypes_

Definition at line 394 of file DDHGCalSiliconModule.cc.

◆ partialTypes_

int HGCalSiliconModule::partialTypes_

Definition at line 393 of file DDHGCalSiliconModule.cc.

◆ rMaxFront_

std::vector<double> HGCalSiliconModule::rMaxFront_

Definition at line 418 of file DDHGCalSiliconModule.cc.

◆ rMinFront_

std::vector<double> HGCalSiliconModule::rMinFront_

Definition at line 415 of file DDHGCalSiliconModule.cc.

◆ rotstr_

std::string HGCalSiliconModule::rotstr_

Definition at line 402 of file DDHGCalSiliconModule.cc.

◆ sectors_

int HGCalSiliconModule::sectors_

Definition at line 401 of file DDHGCalSiliconModule.cc.

◆ sensitiveMode_

int HGCalSiliconModule::sensitiveMode_

Definition at line 397 of file DDHGCalSiliconModule.cc.

◆ slopeB_

std::vector<double> HGCalSiliconModule::slopeB_

Definition at line 413 of file DDHGCalSiliconModule.cc.

◆ slopeT_

std::vector<double> HGCalSiliconModule::slopeT_

Definition at line 416 of file DDHGCalSiliconModule.cc.

◆ thick_

std::vector<double> HGCalSiliconModule::thick_

Definition at line 407 of file DDHGCalSiliconModule.cc.

◆ waferFull_

std::vector<std::string> HGCalSiliconModule::waferFull_

Definition at line 403 of file DDHGCalSiliconModule.cc.

◆ waferIndex_

std::vector<int> HGCalSiliconModule::waferIndex_

Definition at line 420 of file DDHGCalSiliconModule.cc.

◆ waferLayerStart_

std::vector<int> HGCalSiliconModule::waferLayerStart_

Definition at line 422 of file DDHGCalSiliconModule.cc.

◆ waferPart_

std::vector<std::string> HGCalSiliconModule::waferPart_

Definition at line 404 of file DDHGCalSiliconModule.cc.

◆ waferProperty_

std::vector<int> HGCalSiliconModule::waferProperty_

Definition at line 421 of file DDHGCalSiliconModule.cc.

◆ waferSepar_

double HGCalSiliconModule::waferSepar_

Definition at line 400 of file DDHGCalSiliconModule.cc.

◆ waferSize_

double HGCalSiliconModule::waferSize_

Definition at line 399 of file DDHGCalSiliconModule.cc.

◆ waferTypes_

int HGCalSiliconModule::waferTypes_

Definition at line 391 of file DDHGCalSiliconModule.cc.

◆ zFrontB_

std::vector<double> HGCalSiliconModule::zFrontB_

Definition at line 414 of file DDHGCalSiliconModule.cc.

◆ zFrontT_

std::vector<double> HGCalSiliconModule::zFrontT_

Definition at line 417 of file DDHGCalSiliconModule.cc.

◆ zMinBlock_

double HGCalSiliconModule::zMinBlock_

Definition at line 398 of file DDHGCalSiliconModule.cc.