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 > layerOrient_
 
std::vector< int > layers_
 
std::vector< int > layerSense_
 
std::vector< double > layerThick_
 
std::vector< int > layerType_
 
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, 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, dqmdumpme::k, HGCalTypes::layerType(), cms::DDNamespace::material(), g4SimHits_cfi::Material, SiStripPI::min, Skims_PA_cff::name, cms::DDNamespace::name(), HGCalGeometryMode::Polyhedra, PixelTestBeamValidation_cfi::Position, cms::DDNamespace::prepend(), HGCalGeomTools::radius(), makeMuonMisalignmentScenario::rot, cms::DDNamespace::rotation(), cms::rotation_utils::rotName(), AlCaHLTBitMon_QueryRunRegistry::string, to_string(), cms::DDNamespace::volume(), HGCalTypes::WaferCenterR, HGCalWaferIndex::waferLayer(), HGCalProperty::waferOrient(), HGCalProperty::waferPartial(), HGCalProperty::waferThick(), HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), and geometryCSVtoXML::zz.

30  {
31  cms::DDNamespace ns(ctxt, e, true);
32  cms::DDAlgoArguments args(ctxt, e);
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  layerOrient_ = args.value<std::vector<int>>("LayerTypes");
114  for (unsigned int k = 0; k < layerOrient_.size(); ++k)
115  layerOrient_[k] = HGCalTypes::layerType(layerOrient_[k]);
116 #ifdef EDM_ML_DEBUG
117  for (unsigned int i = 0; i < layerOrient_.size(); ++i)
118  edm::LogVerbatim("HGCalGeom") << "LayerTypes [" << i << "] " << layerOrient_[i];
119 #endif
120  if (firstLayer_ > 0) {
121  for (unsigned int i = 0; i < layerType_.size(); ++i) {
122  if (layerSense_[i] > 0) {
123  int ii = layerType_[i];
124  copyNumber_[ii] = (layerSense_[i] == 1) ? firstLayer_ : (firstLayer_ + 1);
125 #ifdef EDM_ML_DEBUG
126  edm::LogVerbatim("HGCalGeom") << "First copy number for layer type " << i << ":" << ii << " with "
127  << materials_[ii] << " changed to " << copyNumber_[ii];
128 #endif
129  }
130  }
131  } else {
132  firstLayer_ = 1;
133  }
134 #ifdef EDM_ML_DEBUG
135  edm::LogVerbatim("HGCalGeom") << "There are " << layerType_.size() << " layers";
136  for (unsigned int i = 0; i < layerType_.size(); ++i)
137  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerType_[i] << " sensitive class "
138  << layerSense_[i];
139 #endif
140  slopeB_ = args.value<std::vector<double>>("SlopeBottom");
141  zFrontB_ = args.value<std::vector<double>>("ZFrontBottom");
142  rMinFront_ = args.value<std::vector<double>>("RMinFront");
143  slopeT_ = args.value<std::vector<double>>("SlopeTop");
144  zFrontT_ = args.value<std::vector<double>>("ZFrontTop");
145  rMaxFront_ = args.value<std::vector<double>>("RMaxFront");
146 #ifdef EDM_ML_DEBUG
147  for (unsigned int i = 0; i < slopeB_.size(); ++i)
148  edm::LogVerbatim("HGCalGeom") << "Bottom Block [" << i << "] Zmin " << cms::convert2mm(zFrontB_[i]) << " Rmin "
149  << cms::convert2mm(rMinFront_[i]) << " Slope " << slopeB_[i];
150  for (unsigned int i = 0; i < slopeT_.size(); ++i)
151  edm::LogVerbatim("HGCalGeom") << "Top Block [" << i << "] Zmin " << cms::convert2mm(zFrontT_[i]) << " Rmax "
152  << cms::convert2mm(rMaxFront_[i]) << " Slope " << slopeT_[i];
153 #endif
154  waferIndex_ = args.value<std::vector<int>>("WaferIndex");
155  waferProperty_ = args.value<std::vector<int>>("WaferProperties");
156  waferLayerStart_ = args.value<std::vector<int>>("WaferLayerStart");
157 #ifdef EDM_ML_DEBUG
158  edm::LogVerbatim("HGCalGeom") << "waferProperties with " << waferIndex_.size() << " entries in "
159  << waferLayerStart_.size() << " layers";
160  for (unsigned int k = 0; k < waferLayerStart_.size(); ++k)
161  edm::LogVerbatim("HGCalGeom") << "LayerStart[" << k << "] " << waferLayerStart_[k];
162  for (unsigned int k = 0; k < waferIndex_.size(); ++k)
163  edm::LogVerbatim("HGCalGeom") << "Wafer[" << k << "] " << waferIndex_[k] << " ("
164  << HGCalWaferIndex::waferLayer(waferIndex_[k]) << ", "
165  << HGCalWaferIndex::waferU(waferIndex_[k]) << ", "
166  << HGCalWaferIndex::waferV(waferIndex_[k]) << ") : ("
167  << HGCalProperty::waferThick(waferProperty_[k]) << ":"
168  << HGCalProperty::waferPartial(waferProperty_[k]) << ":"
169  << HGCalProperty::waferOrient(waferProperty_[k]) << ")";
170  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: NameSpace " << ns.name();
171 #endif
172 
173 #ifdef EDM_ML_DEBUG
174  edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalSiliconModule...";
175  copies_.clear();
176 #endif
177 
178  double zi(zMinBlock_);
179  int laymin(0);
180  for (unsigned int i = 0; i < layers_.size(); ++i) {
181  double zo = zi + layerThick_[i];
182  double routF = HGCalGeomTools::radius(zi, zFrontT_, rMaxFront_, slopeT_);
183  int laymax = laymin + layers_[i];
184  double zz = zi;
185  double thickTot(0);
186  for (int ly = laymin; ly < laymax; ++ly) {
187  int ii = layerType_[ly];
188  int copy = copyNumber_[ii];
189  double hthick = 0.5 * thick_[ii];
190  double rinB = HGCalGeomTools::radius(zo - tol1, zFrontB_, rMinFront_, slopeB_);
191  zz += hthick;
192  thickTot += thick_[ii];
193 
194  std::string name = names_[ii] + std::to_string(copy);
195 #ifdef EDM_ML_DEBUG
196  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: Layer " << ly << ":" << ii << " Front "
197  << cms::convert2mm(zi) << ", " << cms::convert2mm(routF) << " Back "
198  << cms::convert2mm(zo) << ", " << cms::convert2mm(rinB)
199  << " superlayer thickness " << cms::convert2mm(layerThick_[i]);
200 #endif
201 
202  dd4hep::Material matter = ns.material(materials_[ii]);
203  dd4hep::Volume glog;
204 
205  if (layerSense_[ly] < 1) {
206  std::vector<double> pgonZ, pgonRin, pgonRout;
207  double rmax = routF * cosAlpha_ - tol1;
208  HGCalGeomTools::radius(zz - hthick,
209  zz + hthick,
210  zFrontB_,
211  rMinFront_,
212  slopeB_,
213  zFrontT_,
214  rMaxFront_,
215  slopeT_,
216  -layerSense_[ly],
217  pgonZ,
218  pgonRin,
219  pgonRout);
220  for (unsigned int isec = 0; isec < pgonZ.size(); ++isec) {
221  pgonZ[isec] -= zz;
222  if (layerSense_[ly] == 0 || absorbMode_ == 0)
223  pgonRout[isec] = rmax;
224  else
225  pgonRout[isec] = pgonRout[isec] * cosAlpha_ - tol1;
226  }
227  dd4hep::Solid solid = dd4hep::Polyhedra(sectors_, -alpha_, 2._pi, pgonZ, pgonRin, pgonRout);
228  ns.addSolidNS(ns.prepend(name), solid);
229  glog = dd4hep::Volume(solid.name(), solid, matter);
230  ns.addVolumeNS(glog);
231 #ifdef EDM_ML_DEBUG
232  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << solid.name() << " polyhedra of " << sectors_
233  << " sectors covering " << convertRadToDeg(-alpha_) << ":"
234  << convertRadToDeg(-alpha_ + 2._pi) << " with " << pgonZ.size()
235  << " sections and filled with " << matter.name();
236  for (unsigned int k = 0; k < pgonZ.size(); ++k)
237  edm::LogVerbatim("HGCalGeom") << "[" << k << "] z " << cms::convert2mm(pgonZ[k]) << " R "
238  << cms::convert2mm(pgonRin[k]) << ":" << cms::convert2mm(pgonRout[k]);
239 #endif
240  } else {
241  double rins =
242  (sensitiveMode_ < 1) ? rinB : HGCalGeomTools::radius(zz + hthick - tol1, zFrontB_, rMinFront_, slopeB_);
243  double routs =
244  (sensitiveMode_ < 1) ? routF : HGCalGeomTools::radius(zz - hthick, zFrontT_, rMaxFront_, slopeT_);
245  dd4hep::Solid solid = dd4hep::Tube(rins, routs, hthick, 0.0, 2._pi);
246  ns.addSolidNS(ns.prepend(name), solid);
247  glog = dd4hep::Volume(solid.name(), solid, matter);
248  ns.addVolumeNS(glog);
249 
250 #ifdef EDM_ML_DEBUG
251  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << solid.name() << " Tubs made of " << matter.name()
252  << " of dimensions " << cms::convert2mm(rinB) << ":" << cms::convert2mm(rins)
253  << ", " << cms::convert2mm(routF) << ":" << cms::convert2mm(routs) << ", "
254  << cms::convert2mm(hthick) << ", 0.0, 360.0 and position " << glog.name()
255  << " number " << copy << ":" << layerOrient_[copy - firstLayer_];
256 #endif
257  positionSensitive(ctxt, e, glog, (copy - firstLayer_));
258  }
259 
260  dd4hep::Position r1(0, 0, zz);
261  dd4hep::Rotation3D rot;
262 #ifdef EDM_ML_DEBUG
263  std::string rotName("Null");
264 #endif
265  if ((layerSense_[ly] > 0) && (layerOrient_[copy - firstLayer_] == HGCalTypes::WaferCenterR)) {
266  rot = ns.rotation(rotstr_);
267 #ifdef EDM_ML_DEBUG
268  rotName = rotstr_;
269 #endif
270  }
271  mother.placeVolume(glog, copy, dd4hep::Transform3D(rot, r1));
272  int inc = ((layerSense_[ly] > 0) && (facingTypes_ > 1)) ? 2 : 1;
273  copyNumber_[ii] = copy + inc;
274 #ifdef EDM_ML_DEBUG
275  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << glog.name() << " number " << copy
276  << " positioned in " << mother.name() << " at (0,0," << cms::convert2mm(zz)
277  << ") with " << rotName << " rotation";
278 #endif
279  zz += hthick;
280  } // End of loop over layers in a block
281  zi = zo;
282  laymin = laymax;
283  // Make consistency check of all the partitions of the block
284  if (std::abs(thickTot - layerThick_[i]) >= tol2) {
285  if (thickTot > layerThick_[i]) {
286  edm::LogError("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
287  << " is smaller than " << cms::convert2mm(thickTot)
288  << ": thickness of all its components **** ERROR ****";
289  } else {
290  edm::LogWarning("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
291  << " does not match with " << cms::convert2mm(thickTot) << " of the components";
292  }
293  }
294  } // End of loop over blocks
295 
296 #ifdef EDM_ML_DEBUG
297  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << copies_.size() << " different wafer copy numbers";
298  int k(0);
299  for (std::unordered_set<int>::const_iterator itr = copies_.begin(); itr != copies_.end(); ++itr, ++k) {
300  edm::LogVerbatim("HGCalGeom") << "Copy [" << k << "] : " << (*itr);
301  }
302  copies_.clear();
303  edm::LogVerbatim("HGCalGeom") << "<<== End of DDHGCalSiliconModule construction...";
304 #endif
305  }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:134
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
HGCalSiliconModule::copies_
std::unordered_set< int > copies_
Definition: DDHGCalSiliconModule.cc:422
HGCalSiliconModule::slopeT_
std::vector< double > slopeT_
Definition: DDHGCalSiliconModule.cc:415
HGCalSiliconModule::slopeB_
std::vector< double > slopeB_
Definition: DDHGCalSiliconModule.cc:412
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
HGCalSiliconModule::copyNumber_
std::vector< int > copyNumber_
Definition: DDHGCalSiliconModule.cc:407
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
testProducerWithPsetDescEmpty_cfi.i2
i2
Definition: testProducerWithPsetDescEmpty_cfi.py:46
HGCalSiliconModule::layerOrient_
std::vector< int > layerOrient_
Definition: DDHGCalSiliconModule.cc:418
HGCalSiliconModule::layerSense_
std::vector< int > layerSense_
Definition: DDHGCalSiliconModule.cc:411
PixelTestBeamValidation_cfi.Position
Position
Definition: PixelTestBeamValidation_cfi.py:83
g4SimHits_cfi.Material
Material
Definition: g4SimHits_cfi.py:642
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
HGCalSiliconModule::zFrontT_
std::vector< double > zFrontT_
Definition: DDHGCalSiliconModule.cc:416
cms::rotation_utils::rotName
static const std::string & rotName(const T &rot, const cms::DDParsingContext &context)
Definition: DDCoreToDDXMLOutput.cc:81
HGCalSiliconModule::zFrontB_
std::vector< double > zFrontB_
Definition: DDHGCalSiliconModule.cc:413
HGCalSiliconModule::thick_
std::vector< double > thick_
Definition: DDHGCalSiliconModule.cc:406
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
HGCalSiliconModule::names_
std::vector< std::string > names_
Definition: DDHGCalSiliconModule.cc:405
HGCalProperty::waferThick
int32_t waferThick(const int32_t property)
Definition: HGCalProperty.cc:13
HGCalSiliconModule::waferProperty_
std::vector< int > waferProperty_
Definition: DDHGCalSiliconModule.cc:420
cms::Volume
dd4hep::Volume Volume
Definition: DDFilteredView.h:47
cuy.ii
ii
Definition: cuy.py:589
HGCalSiliconModule::waferPart_
std::vector< std::string > waferPart_
Definition: DDHGCalSiliconModule.cc:403
HGCalSiliconModule::rMaxFront_
std::vector< double > rMaxFront_
Definition: DDHGCalSiliconModule.cc:417
HGCalTypes::WaferCenterR
static constexpr int32_t WaferCenterR
Definition: HGCalTypes.h:27
HGCalSiliconModule::waferLayerStart_
std::vector< int > waferLayerStart_
Definition: DDHGCalSiliconModule.cc:421
HGCalSiliconModule::materials_
std::vector< std::string > materials_
Definition: DDHGCalSiliconModule.cc:404
HGCalSiliconModule::layerType_
std::vector< int > layerType_
Definition: DDHGCalSiliconModule.cc:410
testProducerWithPsetDescEmpty_cfi.i1
i1
Definition: testProducerWithPsetDescEmpty_cfi.py:45
HGCalSiliconModule::positionSensitive
void positionSensitive(cms::DDParsingContext &ctxt, xml_h e, const dd4hep::Volume &glog, int layer)
Definition: DDHGCalSiliconModule.cc:307
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
HGCalSiliconModule::layers_
std::vector< int > layers_
Definition: DDHGCalSiliconModule.cc:408
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:128
HGCalTypes::layerType
static int32_t layerType(int type)
Definition: HGCalTypes.cc:42
HGCalSiliconModule::waferFull_
std::vector< std::string > waferFull_
Definition: DDHGCalSiliconModule.cc:402
HGCalSiliconModule::rMinFront_
std::vector< double > rMinFront_
Definition: DDHGCalSiliconModule.cc:414
l1ctLayer1_patternWriters_cff.isec
isec
Definition: l1ctLayer1_patternWriters_cff.py:231
HGCalSiliconModule::waferIndex_
std::vector< int > waferIndex_
Definition: DDHGCalSiliconModule.cc:419
HGCalSiliconModule::layerThick_
std::vector< double > layerThick_
Definition: DDHGCalSiliconModule.cc:409
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17

Member Function Documentation

◆ positionSensitive()

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

Definition at line 307 of file DDHGCalSiliconModule.cc.

References funct::abs(), cms::convert2mm(), filterCSVwithJSON::copy, PVValHelper::dy, MillePedeFileConverter_cfg::e, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, dqmdumpme::k, EgHLTOffHistBins_cfi::nr, HGCalTypes::packTypeUV(), PixelTestBeamValidation_cfi::Position, dttmaxenums::R, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, findQualityFiles::v, cms::DDNamespace::volume(), HGCalTypes::WaferCenterB, HGCalTypes::WaferFull, HGCalWaferIndex::waferLayer(), HGCalProperty::waferOrient(), HGCalProperty::waferPartial(), HGCalProperty::waferThick(), HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

307  {
308  cms::DDNamespace ns(ctxt, e, true);
309  static const double sqrt3 = std::sqrt(3.0);
310  int layercenter = layerOrient_[layer];
311  int layertype = (layerOrient_[layer] == HGCalTypes::WaferCenterB) ? 1 : 0;
312  int firstWafer = waferLayerStart_[layer];
313  int lastWafer = ((layer + 1 < static_cast<int>(waferLayerStart_.size())) ? waferLayerStart_[layer + 1]
314  : static_cast<int>(waferIndex_.size()));
315  double r = 0.5 * (waferSize_ + waferSepar_);
316  double R = 2.0 * r / sqrt3;
317  double dy = 0.75 * R;
318  const auto& xyoff = geomTools_.shiftXY(layercenter, (waferSize_ + waferSepar_));
319 #ifdef EDM_ML_DEBUG
320  int ium(0), ivm(0), kount(0);
321  std::vector<int> ntype(3, 0);
322  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: " << glog.name() << " r " << cms::convert2mm(r) << " R "
323  << cms::convert2mm(R) << " dy " << cms::convert2mm(dy) << " Shift "
324  << cms::convert2mm(xyoff.first) << ":" << cms::convert2mm(xyoff.second)
325  << " WaferSize " << cms::convert2mm((waferSize_ + waferSepar_)) << " index "
326  << firstWafer << ":" << (lastWafer - 1);
327 #endif
328  for (int k = firstWafer; k < lastWafer; ++k) {
329  int u = HGCalWaferIndex::waferU(waferIndex_[k]);
330  int v = HGCalWaferIndex::waferV(waferIndex_[k]);
331 #ifdef EDM_ML_DEBUG
332  int iu = std::abs(u);
333  int iv = std::abs(v);
334 #endif
335  int nr = 2 * v;
336  int nc = -2 * u + v;
337  double xpos = xyoff.first + nc * r;
338  double ypos = xyoff.second + nr * dy;
339  int type = HGCalProperty::waferThick(waferProperty_[k]);
340  int part = HGCalProperty::waferPartial(waferProperty_[k]);
341  int orien = HGCalProperty::waferOrient(waferProperty_[k]);
342  std::string wafer;
343  int i(999);
344  if (part == HGCalTypes::WaferFull) {
345  i = layertype * waferTypes_ + type;
346  wafer = waferFull_[i];
347  } else {
348  i = (part - 1) * waferTypes_ * facingTypes_ * orientationTypes_ + layertype * waferTypes_ * orientationTypes_ +
349  type * orientationTypes_ + orien;
350 #ifdef EDM_ML_DEBUG
351  edm::LogVerbatim("HGCalGeom") << " layertype:type:part:orien:ind " << layertype << ":" << type << ":" << part
352  << ":" << orien << ":" << i << ":" << waferPart_.size();
353 #endif
354  wafer = waferPart_[i];
355  }
356  int copy = HGCalTypes::packTypeUV(type, u, v);
357 #ifdef EDM_ML_DEBUG
358  edm::LogVerbatim("HGCalGeom") << " DDHGCalSiliconModule: Layer" << HGCalWaferIndex::waferLayer(waferIndex_[k])
359  << " Wafer " << wafer << " number " << copy << " type:part:orien:ind " << type
360  << ":" << part << ":" << orien << ":" << i << " layer:u:v:indx "
361  << (layer + firstLayer_) << ":" << u << ":" << v;
362  if (iu > ium)
363  ium = iu;
364  if (iv > ivm)
365  ivm = iv;
366  kount++;
367  if (copies_.count(copy) == 0)
368  copies_.insert(copy);
369 #endif
370  dd4hep::Position tran(xpos, ypos, 0.0);
371  glog.placeVolume(ns.volume(wafer), copy, tran);
372 #ifdef EDM_ML_DEBUG
373  ++ntype[type];
374  edm::LogVerbatim("HGCalGeom") << " DDHGCalSiliconModule: " << wafer << " number " << copy << " type " << layertype
375  << ":" << type << " positioned in " << glog.name() << " at ("
376  << cms::convert2mm(xpos) << "," << cms::convert2mm(ypos) << ",0) with no rotation";
377 #endif
378  }
379 
380 #ifdef EDM_ML_DEBUG
381  edm::LogVerbatim("HGCalGeom") << "DDHGCalSiliconModule: Maximum # of u " << ium << " # of v " << ivm << " and "
382  << kount << " wafers (" << ntype[0] << ":" << ntype[1] << ":" << ntype[2] << ") for "
383  << glog.name();
384 #endif
385  }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:134
EgHLTOffHistBins_cfi.nr
nr
Definition: EgHLTOffHistBins_cfi.py:4
mps_fire.i
i
Definition: mps_fire.py:429
HGCalSiliconModule::copies_
std::unordered_set< int > copies_
Definition: DDHGCalSiliconModule.cc:422
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
cms::convert2mm
constexpr NumType convert2mm(NumType length)
Definition: DDutils.h:7
HGCalSiliconModule::layerOrient_
std::vector< int > layerOrient_
Definition: DDHGCalSiliconModule.cc:418
PixelTestBeamValidation_cfi.Position
Position
Definition: PixelTestBeamValidation_cfi.py:83
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
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
HGCalTypes::WaferFull
static constexpr int32_t WaferFull
Definition: HGCalTypes.h:35
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv
int32_t *__restrict__ iv
Definition: clusterTracksByDensity.h:58
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
HGCalProperty::waferThick
int32_t waferThick(const int32_t property)
Definition: HGCalProperty.cc:13
HGCalSiliconModule::waferProperty_
std::vector< int > waferProperty_
Definition: DDHGCalSiliconModule.cc:420
HGCalSiliconModule::waferPart_
std::vector< std::string > waferPart_
Definition: DDHGCalSiliconModule.cc:403
nano_mu_digi_cff.layer
layer
Definition: nano_mu_digi_cff.py:27
HGCalSiliconModule::waferLayerStart_
std::vector< int > waferLayerStart_
Definition: DDHGCalSiliconModule.cc:421
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
HGCalSiliconModule::waferFull_
std::vector< std::string > waferFull_
Definition: DDHGCalSiliconModule.cc:402
HGCalTypes::WaferCenterB
static constexpr int32_t WaferCenterB
Definition: HGCalTypes.h:26
HGCalSiliconModule::waferIndex_
std::vector< int > waferIndex_
Definition: DDHGCalSiliconModule.cc:419

Member Data Documentation

◆ absorbMode_

int HGCalSiliconModule::absorbMode_

Definition at line 395 of file DDHGCalSiliconModule.cc.

◆ alpha_

double HGCalSiliconModule::alpha_

Definition at line 423 of file DDHGCalSiliconModule.cc.

◆ copies_

std::unordered_set<int> HGCalSiliconModule::copies_

Definition at line 422 of file DDHGCalSiliconModule.cc.

◆ copyNumber_

std::vector<int> HGCalSiliconModule::copyNumber_

Definition at line 407 of file DDHGCalSiliconModule.cc.

◆ cosAlpha_

double HGCalSiliconModule::cosAlpha_

Definition at line 423 of file DDHGCalSiliconModule.cc.

◆ facingTypes_

int HGCalSiliconModule::facingTypes_

Definition at line 391 of file DDHGCalSiliconModule.cc.

◆ firstLayer_

int HGCalSiliconModule::firstLayer_

Definition at line 394 of file DDHGCalSiliconModule.cc.

◆ geomTools_

HGCalGeomTools HGCalSiliconModule::geomTools_

Definition at line 388 of file DDHGCalSiliconModule.cc.

◆ layerOrient_

std::vector<int> HGCalSiliconModule::layerOrient_

Definition at line 418 of file DDHGCalSiliconModule.cc.

◆ layers_

std::vector<int> HGCalSiliconModule::layers_

Definition at line 408 of file DDHGCalSiliconModule.cc.

◆ layerSense_

std::vector<int> HGCalSiliconModule::layerSense_

Definition at line 411 of file DDHGCalSiliconModule.cc.

◆ layerThick_

std::vector<double> HGCalSiliconModule::layerThick_

Definition at line 409 of file DDHGCalSiliconModule.cc.

◆ layerType_

std::vector<int> HGCalSiliconModule::layerType_

Definition at line 410 of file DDHGCalSiliconModule.cc.

◆ materials_

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

Definition at line 404 of file DDHGCalSiliconModule.cc.

◆ names_

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

Definition at line 405 of file DDHGCalSiliconModule.cc.

◆ orientationTypes_

int HGCalSiliconModule::orientationTypes_

Definition at line 393 of file DDHGCalSiliconModule.cc.

◆ partialTypes_

int HGCalSiliconModule::partialTypes_

Definition at line 392 of file DDHGCalSiliconModule.cc.

◆ rMaxFront_

std::vector<double> HGCalSiliconModule::rMaxFront_

Definition at line 417 of file DDHGCalSiliconModule.cc.

◆ rMinFront_

std::vector<double> HGCalSiliconModule::rMinFront_

Definition at line 414 of file DDHGCalSiliconModule.cc.

◆ rotstr_

std::string HGCalSiliconModule::rotstr_

Definition at line 401 of file DDHGCalSiliconModule.cc.

◆ sectors_

int HGCalSiliconModule::sectors_

Definition at line 400 of file DDHGCalSiliconModule.cc.

◆ sensitiveMode_

int HGCalSiliconModule::sensitiveMode_

Definition at line 396 of file DDHGCalSiliconModule.cc.

◆ slopeB_

std::vector<double> HGCalSiliconModule::slopeB_

Definition at line 412 of file DDHGCalSiliconModule.cc.

◆ slopeT_

std::vector<double> HGCalSiliconModule::slopeT_

Definition at line 415 of file DDHGCalSiliconModule.cc.

◆ thick_

std::vector<double> HGCalSiliconModule::thick_

Definition at line 406 of file DDHGCalSiliconModule.cc.

◆ waferFull_

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

Definition at line 402 of file DDHGCalSiliconModule.cc.

◆ waferIndex_

std::vector<int> HGCalSiliconModule::waferIndex_

Definition at line 419 of file DDHGCalSiliconModule.cc.

◆ waferLayerStart_

std::vector<int> HGCalSiliconModule::waferLayerStart_

Definition at line 421 of file DDHGCalSiliconModule.cc.

◆ waferPart_

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

Definition at line 403 of file DDHGCalSiliconModule.cc.

◆ waferProperty_

std::vector<int> HGCalSiliconModule::waferProperty_

Definition at line 420 of file DDHGCalSiliconModule.cc.

◆ waferSepar_

double HGCalSiliconModule::waferSepar_

Definition at line 399 of file DDHGCalSiliconModule.cc.

◆ waferSize_

double HGCalSiliconModule::waferSize_

Definition at line 398 of file DDHGCalSiliconModule.cc.

◆ waferTypes_

int HGCalSiliconModule::waferTypes_

Definition at line 390 of file DDHGCalSiliconModule.cc.

◆ zFrontB_

std::vector<double> HGCalSiliconModule::zFrontB_

Definition at line 413 of file DDHGCalSiliconModule.cc.

◆ zFrontT_

std::vector<double> HGCalSiliconModule::zFrontT_

Definition at line 416 of file DDHGCalSiliconModule.cc.

◆ zMinBlock_

double HGCalSiliconModule::zMinBlock_

Definition at line 397 of file DDHGCalSiliconModule.cc.

Generated for CMSSW Reference Manual by doxygen 1.8.14