CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
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::HGCalSiliconModule ( )
inline

Definition at line 29 of file DDHGCalSiliconModule.cc.

References Exception.

29 { throw cms::Exception("HGCalGeom") << "Wrong initialization to HGCalSiliconModule"; }
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(), mps_fire::i, cuy::ii, isotrackApplyRegressor::k, cms::DDNamespace::material(), min(), mergeVDriftHistosByStation::name, cms::DDNamespace::name(), cms::DDAlgoArguments::parentName(), HGCalGeometryMode::Polyhedra, PixelTestBeamValidation_cfi::Position, cms::DDNamespace::prepend(), diffTwoXMLs::r1, HGCalGeomTools::radius(), makeMuonMisalignmentScenario::rot, cms::DDNamespace::rotation(), cms::rotation_utils::rotName(), AlCaHLTBitMon_QueryRunRegistry::string, cms::DDAlgoArguments::value(), cms::DDNamespace::volume(), HGCalTypes::WaferCenteredRotated, HGCalWaferIndex::waferLayer(), HGCalProperty::waferOrient(), HGCalProperty::waferPartial(), HGCalProperty::waferThick(), HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

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  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] << " ("
162  << HGCalWaferIndex::waferLayer(waferIndex_[k]) << ", "
163  << HGCalWaferIndex::waferU(waferIndex_[k]) << ", "
164  << HGCalWaferIndex::waferV(waferIndex_[k]) << ") : ("
165  << HGCalProperty::waferThick(waferProperty_[k]) << ":"
166  << HGCalProperty::waferPartial(waferProperty_[k]) << ":"
167  << HGCalProperty::waferOrient(waferProperty_[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 
192  std::string name = names_[ii] + std::to_string(copy);
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
263  if ((layerSense_[ly] > 0) && (layerTypes_[copy - firstLayer_] == HGCalTypes::WaferCenteredRotated)) {
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  }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
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:428
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
HGCalSiliconModule::layerSense_
std::vector< int > layerSense_
Definition: DDHGCalSiliconModule.cc:411
PixelTestBeamValidation_cfi.Position
tuple Position
Definition: PixelTestBeamValidation_cfi.py:83
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
cuy.ii
int ii
Definition: cuy.py:589
HGCalProperty::waferOrient
int32_t waferOrient(const int32_t property)
Definition: HGCalProperty.cc:17
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
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::layerTypes_
std::vector< int > layerTypes_
Definition: DDHGCalSiliconModule.cc:418
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:9
min
T min(T a, T b)
Definition: MathUtil.h:58
HGCalSiliconModule::waferProperty_
std::vector< int > waferProperty_
Definition: DDHGCalSiliconModule.cc:420
cms::Volume
dd4hep::Volume Volume
Definition: DDFilteredView.h:47
HGCalSiliconModule::waferPart_
std::vector< std::string > waferPart_
Definition: DDHGCalSiliconModule.cc:403
HGCalSiliconModule::rMaxFront_
std::vector< double > rMaxFront_
Definition: DDHGCalSiliconModule.cc:417
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
HGCalSiliconModule::positionSensitive
void positionSensitive(cms::DDParsingContext &ctxt, xml_h e, const dd4hep::Volume &glog, int layer)
Definition: DDHGCalSiliconModule.cc:305
HGCalProperty::waferPartial
int32_t waferPartial(const int32_t property)
Definition: HGCalProperty.cc:13
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:122
HGCalSiliconModule::waferFull_
std::vector< std::string > waferFull_
Definition: DDHGCalSiliconModule.cc:402
HGCalSiliconModule::rMinFront_
std::vector< double > rMinFront_
Definition: DDHGCalSiliconModule.cc:414
HGCalSiliconModule::waferIndex_
std::vector< int > waferIndex_
Definition: DDHGCalSiliconModule.cc:419
HGCalSiliconModule::layerThick_
std::vector< double > layerThick_
Definition: DDHGCalSiliconModule.cc:409

Member Function Documentation

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, PVValHelper::dy, mps_fire::i, gpuVertexFinder::iv, isotrackApplyRegressor::k, phase1PixelTopology::layer, HGCalTypes::packTypeUV(), PixelTestBeamValidation_cfi::Position, dttmaxenums::R, alignCSCRings::r, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, findQualityFiles::v, cms::DDNamespace::volume(), 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)
309  ? 1
310  : ((layerTypes_[layer] == HGCalTypes::CornerCenteredY) ? 2 : 0);
311  int layertype = (layerTypes_[layer] == HGCalTypes::WaferCenteredBack) ? 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:128
mps_fire.i
i
Definition: mps_fire.py:428
HGCalSiliconModule::copies_
std::unordered_set< int > copies_
Definition: DDHGCalSiliconModule.cc:422
gpuVertexFinder::iv
int32_t *__restrict__ iv
Definition: gpuClusterTracksDBSCAN.h:42
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
PixelTestBeamValidation_cfi.Position
tuple Position
Definition: PixelTestBeamValidation_cfi.py:83
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
HGCalProperty::waferOrient
int32_t waferOrient(const int32_t property)
Definition: HGCalProperty.cc:17
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
phase1PixelTopology::layer
constexpr std::array< uint8_t, layerIndexSize > layer
Definition: phase1PixelTopology.h:110
HGCalSiliconModule::layerTypes_
std::vector< int > layerTypes_
Definition: DDHGCalSiliconModule.cc:418
HGCalGeomTools::shiftXY
std::pair< double, double > shiftXY(int waferPosition, double waferSize) const
Definition: HGCalGeomTools.cc:234
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:9
HGCalSiliconModule::waferProperty_
std::vector< int > waferProperty_
Definition: DDHGCalSiliconModule.cc:420
HGCalSiliconModule::waferPart_
std::vector< std::string > waferPart_
Definition: DDHGCalSiliconModule.cc:403
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:13
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
HGCalSiliconModule::waferIndex_
std::vector< int > waferIndex_
Definition: DDHGCalSiliconModule.cc:419

Member Data Documentation

int HGCalSiliconModule::absorbMode_

Definition at line 395 of file DDHGCalSiliconModule.cc.

double HGCalSiliconModule::alpha_

Definition at line 423 of file DDHGCalSiliconModule.cc.

std::unordered_set<int> HGCalSiliconModule::copies_

Definition at line 422 of file DDHGCalSiliconModule.cc.

std::vector<int> HGCalSiliconModule::copyNumber_

Definition at line 407 of file DDHGCalSiliconModule.cc.

double HGCalSiliconModule::cosAlpha_

Definition at line 423 of file DDHGCalSiliconModule.cc.

int HGCalSiliconModule::facingTypes_

Definition at line 391 of file DDHGCalSiliconModule.cc.

int HGCalSiliconModule::firstLayer_

Definition at line 394 of file DDHGCalSiliconModule.cc.

HGCalGeomTools HGCalSiliconModule::geomTools_

Definition at line 388 of file DDHGCalSiliconModule.cc.

std::vector<int> HGCalSiliconModule::layers_

Definition at line 408 of file DDHGCalSiliconModule.cc.

std::vector<int> HGCalSiliconModule::layerSense_

Definition at line 411 of file DDHGCalSiliconModule.cc.

std::vector<double> HGCalSiliconModule::layerThick_

Definition at line 409 of file DDHGCalSiliconModule.cc.

std::vector<int> HGCalSiliconModule::layerType_

Definition at line 410 of file DDHGCalSiliconModule.cc.

std::vector<int> HGCalSiliconModule::layerTypes_

Definition at line 418 of file DDHGCalSiliconModule.cc.

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

Definition at line 404 of file DDHGCalSiliconModule.cc.

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

Definition at line 405 of file DDHGCalSiliconModule.cc.

int HGCalSiliconModule::orientationTypes_

Definition at line 393 of file DDHGCalSiliconModule.cc.

int HGCalSiliconModule::partialTypes_

Definition at line 392 of file DDHGCalSiliconModule.cc.

std::vector<double> HGCalSiliconModule::rMaxFront_

Definition at line 417 of file DDHGCalSiliconModule.cc.

std::vector<double> HGCalSiliconModule::rMinFront_

Definition at line 414 of file DDHGCalSiliconModule.cc.

std::string HGCalSiliconModule::rotstr_

Definition at line 401 of file DDHGCalSiliconModule.cc.

int HGCalSiliconModule::sectors_

Definition at line 400 of file DDHGCalSiliconModule.cc.

int HGCalSiliconModule::sensitiveMode_

Definition at line 396 of file DDHGCalSiliconModule.cc.

std::vector<double> HGCalSiliconModule::slopeB_

Definition at line 412 of file DDHGCalSiliconModule.cc.

std::vector<double> HGCalSiliconModule::slopeT_

Definition at line 415 of file DDHGCalSiliconModule.cc.

std::vector<double> HGCalSiliconModule::thick_

Definition at line 406 of file DDHGCalSiliconModule.cc.

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

Definition at line 402 of file DDHGCalSiliconModule.cc.

std::vector<int> HGCalSiliconModule::waferIndex_

Definition at line 419 of file DDHGCalSiliconModule.cc.

std::vector<int> HGCalSiliconModule::waferLayerStart_

Definition at line 421 of file DDHGCalSiliconModule.cc.

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

Definition at line 403 of file DDHGCalSiliconModule.cc.

std::vector<int> HGCalSiliconModule::waferProperty_

Definition at line 420 of file DDHGCalSiliconModule.cc.

double HGCalSiliconModule::waferSepar_

Definition at line 399 of file DDHGCalSiliconModule.cc.

double HGCalSiliconModule::waferSize_

Definition at line 398 of file DDHGCalSiliconModule.cc.

int HGCalSiliconModule::waferTypes_

Definition at line 390 of file DDHGCalSiliconModule.cc.

std::vector<double> HGCalSiliconModule::zFrontB_

Definition at line 413 of file DDHGCalSiliconModule.cc.

std::vector<double> HGCalSiliconModule::zFrontT_

Definition at line 416 of file DDHGCalSiliconModule.cc.

double HGCalSiliconModule::zMinBlock_

Definition at line 397 of file DDHGCalSiliconModule.cc.

Generated for CMSSW Reference Manual by doxygen 1.8.5