CMS 3D CMS Logo
DDHGCalEEAlgo.cc
Go to the documentation of this file.
1 /*
2  * DD4hep_HGCalEEAlgo.cc
3  *
4  * Created on: 27-Aug-2019
5  * Author: rsehgal
6  *
7  * DD4hep code for, HGCalEEAlgo developed by Sunanda Banerjee
8  */
9 
10 #include <cmath>
11 #include <memory>
12 #include <string>
13 #include <unordered_set>
14 #include <vector>
15 
16 #include "DD4hep/DetFactoryHelper.h"
17 #include "DataFormats/Math/interface/angle_units.h"
18 #include "DetectorDescription/DDCMS/interface/DDPlugins.h"
19 #include "DetectorDescription/DDCMS/interface/DDutils.h"
20 #include "FWCore/MessageLogger/interface/MessageLogger.h"
21 #include "Geometry/HGCalCommonData/interface/HGCalGeomTools.h"
22 #include "Geometry/HGCalCommonData/interface/HGCalParameters.h"
23 #include "Geometry/HGCalCommonData/interface/HGCalTypes.h"
24 #include "Geometry/HGCalCommonData/interface/HGCalWaferType.h"
25 
26 //#define EDM_ML_DEBUG
27 using namespace angle_units::operators;
28 
29 struct HGCalEEAlgo {
30  HGCalGeomTools geomTools_;
31  std::unique_ptr<HGCalWaferType> waferType_;
32  dd4hep::Volume mother_;
33 
34  std::vector<std::string> wafers_; // Wafers
35  std::vector<std::string> materials_; // Materials
36  std::vector<std::string> names_; // Names
37  std::vector<double> thick_; // Thickness of the material
38  std::vector<int> copyNumber_; // Initial copy numbers
39  std::vector<int> layers_; // Number of layers in a section
40  std::vector<double> layerThick_; // Thickness of each section
41  std::vector<int> layerType_; // Type of the layer
42  std::vector<int> layerSense_; // Content of a layer (sensitive?)
43  std::vector<int> layerCenter_; // Centering of the wafers
44  int firstLayer_; // Copy # of the first sensitive layer
45  int absorbMode_; // Absorber mode
46  int sensitiveMode_; // Sensitive mode
47  double zMinBlock_; // Starting z-value of the block
48  std::vector<double> rad100to200_; // Parameters for 120-200mum trans.
49  std::vector<double> rad200to300_; // Parameters for 200-300mum trans.
50  double zMinRadPar_; // Minimum z for radius parametriz.
51  int choiceType_; // Type of parametrization to be used
52  int nCutRadPar_; // Cut off threshold for corners
53  double fracAreaMin_; // Minimum fractional conatined area
54  double waferSize_; // Width of the wafer
55  double waferSepar_; // Sensor separation
56  int sectors_; // Sectors
57  std::vector<double> slopeB_; // Slope at the lower R
58  std::vector<double> zFrontB_; // Starting Z values for the slopes
59  std::vector<double> rMinFront_; // Corresponding rMin's
60  std::vector<double> slopeT_; // Slopes at the larger R
61  std::vector<double> zFrontT_; // Starting Z values for the slopes
62  std::vector<double> rMaxFront_; // Corresponding rMax's
63  std::unordered_set<int> copies_; // List of copy #'s
64  double alpha_, cosAlpha_;
65 
66  HGCalEEAlgo() = delete;
67 
68  HGCalEEAlgo(cms::DDParsingContext& ctxt, xml_h e) {
69  cms::DDNamespace ns(ctxt, e, true);
70  cms::DDAlgoArguments args(ctxt, e);
71 
72  mother_ = ns.volume(args.parentName());
73  wafers_ = args.value<std::vector<std::string>>("WaferNames");
74 #ifdef EDM_ML_DEBUG
75  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: " << wafers_.size() << " wafers";
76  for (unsigned int i = 0; i < wafers_.size(); ++i)
77  edm::LogVerbatim("HGCalGeom") << "Wafer[" << i << "] " << wafers_[i];
78 #endif
79 
80  materials_ = args.value<std::vector<std::string>>("MaterialNames");
81  names_ = args.value<std::vector<std::string>>("VolumeNames");
82  thick_ = args.value<std::vector<double>>("Thickness");
83  copyNumber_.resize(materials_.size(), 1);
84 #ifdef EDM_ML_DEBUG
85  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: " << materials_.size() << " types of volumes";
86  for (unsigned int i = 0; i < names_.size(); ++i)
87  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << names_[i] << " of thickness "
88  << cms::convert2mm(thick_[i]) << " filled with " << materials_[i]
89  << " first copy number " << copyNumber_[i];
90 #endif
91 
92  layers_ = args.value<std::vector<int>>("Layers");
93  layerThick_ = args.value<std::vector<double>>("LayerThick");
94 #ifdef EDM_ML_DEBUG
95  edm::LogVerbatim("HGCalGeom") << "There are " << layers_.size() << " blocks";
96  for (unsigned int i = 0; i < layers_.size(); ++i)
97  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] of thickness " << cms::convert2mm(layerThick_[i])
98  << " with " << layers_[i] << " layers";
99 #endif
100 
101  layerType_ = args.value<std::vector<int>>("LayerType");
102  layerSense_ = args.value<std::vector<int>>("LayerSense");
103  firstLayer_ = args.value<int>("FirstLayer");
104  absorbMode_ = args.value<int>("AbsorberMode");
105  sensitiveMode_ = args.value<int>("SensitiveMode");
106 #ifdef EDM_ML_DEBUG
107  edm::LogVerbatim("HGCalGeom") << "First Layer " << firstLayer_ << " and "
108  << "Absober:Sensitive mode " << absorbMode_ << ":" << sensitiveMode_;
109 #endif
110  layerCenter_ = args.value<std::vector<int>>("LayerCenter");
111 #ifdef EDM_ML_DEBUG
112  for (unsigned int i = 0; i < layerCenter_.size(); ++i)
113  edm::LogVerbatim("HGCalGeom") << "LayerCenter [" << i << "] " << layerCenter_[i];
114 #endif
115  if (firstLayer_ > 0) {
116  for (unsigned int i = 0; i < layerType_.size(); ++i) {
117  if (layerSense_[i] > 0) {
118  int ii = layerType_[i];
119  copyNumber_[ii] = firstLayer_;
120 #ifdef EDM_ML_DEBUG
121  edm::LogVerbatim("HGCalGeom") << "First copy number for layer type " << i << ":" << ii << " with "
122  << materials_[ii] << " changed to " << copyNumber_[ii];
123 #endif
124  break;
125  }
126  }
127  } else {
128  firstLayer_ = 1;
129  }
130 #ifdef EDM_ML_DEBUG
131  edm::LogVerbatim("HGCalGeom") << "There are " << layerType_.size() << " layers";
132  for (unsigned int i = 0; i < layerType_.size(); ++i)
133  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerType_[i] << " sensitive class "
134  << layerSense_[i];
135 #endif
136  zMinBlock_ = args.value<double>("zMinBlock");
137 
138  rad100to200_ = args.value<std::vector<double>>("rad100to200");
139  rad200to300_ = args.value<std::vector<double>>("rad200to300");
140  zMinRadPar_ = args.value<double>("zMinForRadPar");
141  choiceType_ = args.value<int>("choiceType");
142  nCutRadPar_ = args.value<int>("nCornerCut");
143  fracAreaMin_ = args.value<double>("fracAreaMin");
144  waferSize_ = args.value<double>("waferSize");
145  waferSepar_ = args.value<double>("SensorSeparation");
146  sectors_ = args.value<int>("Sectors");
147  alpha_ = (1._pi) / sectors_;
148  cosAlpha_ = cos(alpha_);
149 #ifdef EDM_ML_DEBUG
150  edm::LogVerbatim("HGCalGeom") << "zStart " << cms::convert2mm(zMinBlock_)
151  << " radius for wafer type separation uses " << rad100to200_.size()
152  << " parameters; zmin " << cms::convert2mm(zMinRadPar_) << " cutoff " << choiceType_
153  << ":" << nCutRadPar_ << ":" << fracAreaMin_ << " wafer width "
154  << cms::convert2mm(waferSize_) << " separations " << cms::convert2mm(waferSepar_)
155  << " sectors " << sectors_ << ":" << convertRadToDeg(alpha_) << ":" << cosAlpha_;
156  for (unsigned int k = 0; k < rad100to200_.size(); ++k)
157  edm::LogVerbatim("HGCalGeom") << "[" << k << "] 100-200 " << rad100to200_[k] << " 200-300 " << rad200to300_[k];
158 #endif
159 
160  slopeB_ = args.value<std::vector<double>>("SlopeBottom");
161  zFrontB_ = args.value<std::vector<double>>("ZFrontBottom");
162  rMinFront_ = args.value<std::vector<double>>("RMinFront");
163  slopeT_ = args.value<std::vector<double>>("SlopeTop");
164  zFrontT_ = args.value<std::vector<double>>("ZFrontTop");
165  rMaxFront_ = args.value<std::vector<double>>("RMaxFront");
166 #ifdef EDM_ML_DEBUG
167  for (unsigned int i = 0; i < slopeB_.size(); ++i)
168  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] Zmin " << cms::convert2mm(zFrontB_[i]) << " Rmin "
169  << cms::convert2mm(rMinFront_[i]) << " Slope " << slopeB_[i];
170  for (unsigned int i = 0; i < slopeT_.size(); ++i)
171  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] Zmin " << cms::convert2mm(zFrontT_[i]) << " Rmax "
172  << cms::convert2mm(rMaxFront_[i]) << " Slope " << slopeT_[i];
173 #endif
174 
175 #ifdef EDM_ML_DEBUG
176  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: NameSpace " << ns.name();
177 #endif
178 
179  waferType_ = std::make_unique<HGCalWaferType>(rad100to200_,
180  rad200to300_,
181  cms::convert2mm((waferSize_ + waferSepar_)),
182  cms::convert2mm(zMinRadPar_),
183  choiceType_,
184  nCutRadPar_,
185  fracAreaMin_);
186 
187  ConstructAlgo(ctxt, e);
188  }
189 
190  void ConstructAlgo(cms::DDParsingContext& ctxt, xml_h e) {
191 #ifdef EDM_ML_DEBUG
192  edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalEEAlgo...";
193  copies_.clear();
194 #endif
195  dd4hep::Volume par;
196  ConstructLayers(par, ctxt, e);
197 #ifdef EDM_ML_DEBUG
198  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: " << copies_.size() << " different wafer copy numbers";
199  int k(0);
200  for (std::unordered_set<int>::const_iterator itr = copies_.begin(); itr != copies_.end(); ++itr, ++k) {
201  edm::LogVerbatim("HGCalGeom") << "Copy [" << k << "] : " << (*itr);
202  }
203  copies_.clear();
204  edm::LogVerbatim("HGCalGeom") << "<<== End of DDHGCalEEAlgo construction...";
205 #endif
206  }
207 
208  void ConstructLayers(const dd4hep::Volume module, cms::DDParsingContext& ctxt, xml_h e) {
209  static constexpr double tol1 = 0.01 * dd4hep::mm;
210  static constexpr double tol2 = 0.00001 * dd4hep::mm;
211  cms::DDNamespace ns(ctxt, e, true);
212 
213 #ifdef EDM_ML_DEBUG
214  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: \t\tInside Layers";
215 #endif
216 
217  double zi(zMinBlock_);
218  int laymin(0);
219  for (unsigned int i = 0; i < layers_.size(); i++) {
220  double zo = zi + layerThick_[i];
221  double routF = HGCalGeomTools::radius(zi, zFrontT_, rMaxFront_, slopeT_);
222  int laymax = laymin + layers_[i];
223  double zz = zi;
224  double thickTot(0);
225  for (int ly = laymin; ly < laymax; ++ly) {
226  int ii = layerType_[ly];
227  int copy = copyNumber_[ii];
228  double hthick = 0.5 * thick_[ii];
229  double rinB = HGCalGeomTools::radius(zo - tol1, zFrontB_, rMinFront_, slopeB_);
230  zz += hthick;
231  thickTot += thick_[ii];
232 
233  std::string name = ns.prepend(names_[ii]) + std::to_string(copy);
234 #ifdef EDM_ML_DEBUG
235  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: Layer " << ly << ":" << ii << " Front " << cms::convert2mm(zi)
236  << ", " << cms::convert2mm(routF) << " Back " << cms::convert2mm(zo) << ", "
237  << cms::convert2mm(rinB) << " superlayer thickness "
238  << cms::convert2mm(layerThick_[i]);
239 #endif
240 
241  std::string matName = materials_[ii];
242  dd4hep::Material matter = ns.material(matName);
243  dd4hep::Volume glog;
244  if (layerSense_[ly] < 1) {
245  std::vector<double> pgonZ, pgonRin, pgonRout;
246  if (layerSense_[ly] == 0 || absorbMode_ == 0) {
247  double rmax = routF * cosAlpha_ - tol1;
248  pgonZ.emplace_back(-hthick);
249  pgonZ.emplace_back(hthick);
250  pgonRin.emplace_back(rinB);
251  pgonRin.emplace_back(rinB);
252  pgonRout.emplace_back(rmax);
253  pgonRout.emplace_back(rmax);
254  } else {
255  HGCalGeomTools::radius(zz - hthick,
256  zz + hthick,
257  zFrontB_,
258  rMinFront_,
259  slopeB_,
260  zFrontT_,
261  rMaxFront_,
262  slopeT_,
263  -layerSense_[ly],
264  pgonZ,
265  pgonRin,
266  pgonRout);
267 #ifdef EDM_ML_DEBUG
268  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: z " << cms::convert2mm((zz - hthick)) << ":"
269  << cms::convert2mm((zz + hthick)) << " with " << pgonZ.size() << " palnes";
270  for (unsigned int isec = 0; isec < pgonZ.size(); ++isec)
271  edm::LogVerbatim("HGCalGeom") << "[" << isec << "] z " << cms::convert2mm(pgonZ[isec]) << " R "
272  << cms::convert2mm(pgonRin[isec]) << ":" << cms::convert2mm(pgonRout[isec]);
273 #endif
274  for (unsigned int isec = 0; isec < pgonZ.size(); ++isec) {
275  pgonZ[isec] -= zz;
276  pgonRout[isec] = pgonRout[isec] * cosAlpha_ - tol1;
277  }
278  }
279 
280  dd4hep::Solid solid = dd4hep::Polyhedra(sectors_, -alpha_, 2._pi, pgonZ, pgonRin, pgonRout);
281  ns.addSolidNS(ns.prepend(name), solid);
282  glog = dd4hep::Volume(solid.name(), solid, matter);
283  ns.addVolumeNS(glog);
284 
285 #ifdef EDM_ML_DEBUG
286  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: " << solid.name() << " polyhedra of " << sectors_
287  << " sectors covering " << convertRadToDeg(-alpha_) << ":"
288  << convertRadToDeg(-alpha_ + 2._pi) << " with " << pgonZ.size()
289  << " sections and filled with " << matName;
290 
291  for (unsigned int k = 0; k < pgonZ.size(); ++k)
292  edm::LogVerbatim("HGCalGeom") << "[" << k << "] z " << cms::convert2mm(pgonZ[k]) << " R "
293  << cms::convert2mm(pgonRin[k]) << ":" << cms::convert2mm(pgonRout[k]);
294 #endif
295  } else {
296  double rins =
297  (sensitiveMode_ < 1) ? rinB : HGCalGeomTools::radius(zz + hthick - tol1, zFrontB_, rMinFront_, slopeB_);
298  double routs =
299  (sensitiveMode_ < 1) ? routF : HGCalGeomTools::radius(zz - hthick, zFrontT_, rMaxFront_, slopeT_);
300  dd4hep::Solid solid = dd4hep::Tube(rins, routs, hthick, 0.0, 2._pi);
301  ns.addSolidNS(ns.prepend(name), solid);
302  glog = dd4hep::Volume(solid.name(), solid, matter);
303  ns.addVolumeNS(glog);
304 
305 #ifdef EDM_ML_DEBUG
306  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: " << solid.name() << " Tubs made of " << matter.name()
307  << " of dimensions " << cms::convert2mm(rinB) << ":" << cms::convert2mm(rins)
308  << ", " << cms::convert2mm(routF) << ":" << cms::convert2mm(routs) << ", "
309  << cms::convert2mm(hthick) << ", 0.0, 360.0 and position " << glog.name()
310  << " number " << copy << ":" << layerCenter_[copy - firstLayer_];
311 #endif
312  PositionSensitive(
313  ctxt, e, glog, rins, routs, zz, layerSense_[ly], layerCenter_[copy - firstLayer_]); //, cpv);
314  }
315 
316  dd4hep::Position r1(0, 0, zz);
317  mother_.placeVolume(glog, copy, r1);
318  ++copyNumber_[ii];
319 
320 #ifdef EDM_ML_DEBUG
321  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: " << glog.name() << " number " << copy << " positioned in "
322  << module.name() << " at (0,0," << cms::convert2mm(zz) << ") with no rotation";
323 #endif
324  zz += hthick;
325  } // End of loop over layers in a block
326  zi = zo;
327  laymin = laymax;
328  if (std::abs(thickTot - layerThick_[i]) >= tol2) {
329  if (thickTot > layerThick_[i]) {
330  edm::LogError("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
331  << " is smaller than " << cms::convert2mm(thickTot)
332  << ": thickness of all its components **** ERROR ****";
333  } else {
334  edm::LogWarning("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
335  << " does not match with " << cms::convert2mm(thickTot) << " of the components";
336  }
337  }
338 
339  } // End of loop over layers in a block
340  }
341 
342  void PositionSensitive(cms::DDParsingContext& ctxt,
343  xml_h e,
344  const dd4hep::Volume& glog,
345  double rin,
346  double rout,
347  double zpos,
348  int layertype,
349  int layercenter) {
350  cms::DDNamespace ns(ctxt, e, true);
351  static const double sqrt3 = std::sqrt(3.0);
352  double r = 0.5 * (waferSize_ + waferSepar_);
353  double R = 2.0 * r / sqrt3;
354  double dy = 0.75 * R;
355  int N = (int)(0.5 * rout / r) + 2;
356  const auto& xyoff = geomTools_.shiftXY(layercenter, (waferSize_ + waferSepar_));
357 #ifdef EDM_ML_DEBUG
358  int ium(0), ivm(0), iumAll(0), ivmAll(0), kount(0), ntot(0), nin(0);
359  std::vector<int> ntype(6, 0);
360  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: " << glog.name() << " rout " << cms::convert2mm(rout) << " N " << N
361  << " for maximum u, v; r " << cms::convert2mm(r) << " R " << cms::convert2mm(R)
362  << " dy " << cms::convert2mm(dy) << " Shift " << cms::convert2mm(xyoff.first) << ":"
363  << cms::convert2mm(xyoff.second) << " WaferSize "
364  << cms::convert2mm((waferSize_ + waferSepar_));
365 #endif
366 
367  for (int u = -N; u <= N; ++u) {
368  for (int v = -N; v <= N; ++v) {
369  int nr = 2 * v;
370  int nc = -2 * u + v;
371  double xpos = xyoff.first + nc * r;
372  double ypos = xyoff.second + nr * dy;
373  const auto& corner = HGCalGeomTools::waferCorner(xpos, ypos, r, R, rin, rout, false);
374 #ifdef EDM_ML_DEBUG
375  int iu = std::abs(u);
376  int iv = std::abs(v);
377  ++ntot;
378  if (((corner.first <= 0) && std::abs(u) < 5 && std::abs(v) < 5) || (std::abs(u) < 2 && std::abs(v) < 2)) {
379  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEAlgo: " << glog.name() << " R " << cms::convert2mm(rin) << ":"
380  << cms::convert2mm(rout) << "\n Z " << cms::convert2mm(zpos) << " LayerType "
381  << layertype << " u " << u << " v " << v << " with " << corner.first
382  << " corners";
383  }
384 #endif
385  if (corner.first > 0) {
386  int type = waferType_->getType(cms::convert2mm(xpos), cms::convert2mm(ypos), cms::convert2mm(zpos));
387  int copy = HGCalTypes::packTypeUV(type, u, v);
388 #ifdef EDM_ML_DEBUG
389  if (iu > ium)
390  ium = iu;
391  if (iv > ivm)
392  ivm = iv;
393  kount++;
394  if (copies_.count(copy) == 0)
395  copies_.insert(copy);
396 #endif
397  if (corner.first == (int)(HGCalParameters::k_CornerSize)) {
398 #ifdef EDM_ML_DEBUG
399  if (iu > iumAll)
400  iumAll = iu;
401  if (iv > ivmAll)
402  ivmAll = iv;
403  ++nin;
404 #endif
405 
406  dd4hep::Position tran(xpos, ypos, 0.0);
407  if (layertype > 1)
408  type += 3;
409  glog.placeVolume(ns.volume(wafers_[type]), copy, tran);
410 #ifdef EDM_ML_DEBUG
411  ++ntype[type];
412  edm::LogVerbatim("HGCalGeom")
413  << " DDHGCalEEAlgo: " << wafers_[type] << " number " << copy << " positioned in " << glog.name()
414  << " at (" << cms::convert2mm(xpos) << ", " << cms::convert2mm(ypos) << ",0) with no rotation";
415 #endif
416  }
417  }
418  }
419  }
420 
421 #ifdef EDM_ML_DEBUG
422  edm::LogVerbatim("HGCalGeom") << " DDHGCalEEAlgo: Maximum # of u " << ium << ":" << iumAll << " # of v " << ivm
423  << ":" << ivmAll << " and " << nin << ":" << kount << ":" << ntot << " wafers ("
424  << ntype[0] << ":" << ntype[1] << ":" << ntype[2] << ":" << ntype[3] << ":"
425  << ntype[4] << ":" << ntype[5] << ") for " << glog.name() << " R "
426  << cms::convert2mm(rin) << ":" << cms::convert2mm(rout);
427 #endif
428  }
429 };
430 
431 static long algorithm(dd4hep::Detector& /* description */, cms::DDParsingContext& ctxt, xml_h e) {
432  HGCalEEAlgo eealgo(ctxt, e);
433  return cms::s_executed;
434 }
435 
436 DECLARE_DDCMS_DETELEMENT(DDCMS_hgcal_DDHGCalEEAlgo, algorithm)
HGCalEEAlgo::geomTools_
HGCalGeomTools geomTools_
Definition: DDHGCalEEAlgo.cc:30
HGCalEEAlgo::choiceType_
int choiceType_
Definition: DDHGCalEEAlgo.cc:51
HGCalEEAlgo::sectors_
int sectors_
Definition: DDHGCalEEAlgo.cc:56
writedatasetfile.args
args
Definition: writedatasetfile.py:18
HGCalGeomTools::waferCorner
static std::pair< int32_t, int32_t > waferCorner(double xpos, double ypos, double r, double R, double rMin, double rMax, bool oldBug=false)
Definition: HGCalGeomTools.cc:287
mps_fire.i
i
Definition: mps_fire.py:428
geometryCSVtoXML.zz
zz
Definition: geometryCSVtoXML.py:19
HGCalEEAlgo::thick_
std::vector< double > thick_
Definition: DDHGCalEEAlgo.cc:37
g4SimHits_cfi.Material
Material
Definition: g4SimHits_cfi.py:594
MessageLogger.h
DECLARE_DDCMS_DETELEMENT
#define DECLARE_DDCMS_DETELEMENT(name, func)
Definition: DDPlugins.h:25
cms::DDNamespace::material
dd4hep::Material material(const std::string &name) const
Definition: DDNamespace.cc:154
filterCSVwithJSON.copy
copy
Definition: filterCSVwithJSON.py:36
HGCalEEAlgo::layerSense_
std::vector< int > layerSense_
Definition: DDHGCalEEAlgo.cc:42
HGCalEEAlgo::zFrontT_
std::vector< double > zFrontT_
Definition: DDHGCalEEAlgo.cc:61
cms::DDParsingContext
Definition: DDParsingContext.h:14
HGCalEEAlgo::rad200to300_
std::vector< double > rad200to300_
Definition: DDHGCalEEAlgo.cc:49
HGCalEEAlgo::waferSize_
double waferSize_
Definition: DDHGCalEEAlgo.cc:54
angle_units::operators::convertRadToDeg
constexpr NumType convertRadToDeg(NumType radians)
Definition: angle_units.h:21
cms::DDNamespace
Definition: DDNamespace.h:22
gpuVertexFinder::iv
int32_t *__restrict__ iv
Definition: gpuClusterTracksDBSCAN.h:42
HGCalParameters::k_CornerSize
static constexpr uint32_t k_CornerSize
Definition: HGCalParameters.h:44
nin
int nin
Definition: CascadeWrapper.h:114
HGCalEEAlgo::sensitiveMode_
int sensitiveMode_
Definition: DDHGCalEEAlgo.cc:46
HGCalEEAlgo::layerThick_
std::vector< double > layerThick_
Definition: DDHGCalEEAlgo.cc:40
gpuPixelDoublets::ntot
__shared__ uint32_t ntot
Definition: gpuPixelDoubletsAlgos.h:67
HGCalGeomTools::shiftXY
std::pair< double, double > shiftXY(int waferPosition, double waferSize) const
Definition: HGCalGeomTools.cc:234
findQualityFiles.v
v
Definition: findQualityFiles.py:179
HGCalEEAlgo::slopeT_
std::vector< double > slopeT_
Definition: DDHGCalEEAlgo.cc:60
HGCalEEAlgo::zMinBlock_
double zMinBlock_
Definition: DDHGCalEEAlgo.cc:47
HGCalEEAlgo
Definition: DDHGCalEEAlgo.cc:29
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
HGCalEEAlgo::rMinFront_
std::vector< double > rMinFront_
Definition: DDHGCalEEAlgo.cc:59
HGCalEEAlgo::cosAlpha_
double cosAlpha_
Definition: DDHGCalEEAlgo.cc:64
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
HGCalEEAlgo::HGCalEEAlgo
HGCalEEAlgo(cms::DDParsingContext &ctxt, xml_h e)
Definition: DDHGCalEEAlgo.cc:68
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
cms::DDNamespace::prepend
std::string prepend(const std::string &) const
Definition: DDNamespace.cc:98
HGCalWaferType.h
N
#define N
Definition: blowfish.cc:9
HGCalEEAlgo::wafers_
std::vector< std::string > wafers_
Definition: DDHGCalEEAlgo.cc:34
HGCalEEAlgo::rad100to200_
std::vector< double > rad100to200_
Definition: DDHGCalEEAlgo.cc:48
PixelTestBeamValidation_cfi.Position
Position
Definition: PixelTestBeamValidation_cfi.py:83
dqmdumpme.k
k
Definition: dqmdumpme.py:60
HGCalTypes::packTypeUV
static int32_t packTypeUV(int type, int u, int v)
Definition: HGCalTypes.cc:3
HGCalEEAlgo::copies_
std::unordered_set< int > copies_
Definition: DDHGCalEEAlgo.cc:63
cms::Volume
dd4hep::Volume Volume
Definition: DDFilteredView.h:47
HGCalEEAlgo::waferSepar_
double waferSepar_
Definition: DDHGCalEEAlgo.cc:55
EgHLTOffHistBins_cfi.nr
nr
Definition: EgHLTOffHistBins_cfi.py:4
HGCalEEAlgo::zFrontB_
std::vector< double > zFrontB_
Definition: DDHGCalEEAlgo.cc:58
type
type
Definition: SiPixelVCal_PayloadInspector.cc:39
HGCalEEAlgo::layerType_
std::vector< int > layerType_
Definition: DDHGCalEEAlgo.cc:41
angle_units::operators
Definition: angle_units.h:11
gainCalibHelper::gainCalibPI::type
type
Definition: SiPixelGainCalibHelper.h:40
DDPlugins.h
createfilelist.int
int
Definition: createfilelist.py:10
cms::convert2mm
constexpr NumType convert2mm(NumType length)
Definition: DDutils.h:7
HGCalEEAlgo::copyNumber_
std::vector< int > copyNumber_
Definition: DDHGCalEEAlgo.cc:38
HGCalGeomTools
Definition: HGCalGeomTools.h:8
cms::DDNamespace::addSolidNS
dd4hep::Solid addSolidNS(const std::string &name, dd4hep::Solid solid) const
Definition: DDNamespace.cc:280
cms::DDAlgoArguments
Definition: DDAlgoArguments.h:28
HGCalEEAlgo::materials_
std::vector< std::string > materials_
Definition: DDHGCalEEAlgo.cc:35
PVValHelper::dy
Definition: PVValidationHelpers.h:50
HGCalEEAlgo::firstLayer_
int firstLayer_
Definition: DDHGCalEEAlgo.cc:44
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
HGCalEEAlgo::names_
std::vector< std::string > names_
Definition: DDHGCalEEAlgo.cc:36
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
algorithm
static long algorithm(dd4hep::Detector &, cms::DDParsingContext &ctxt, xml_h e)
Definition: DDHGCalEEAlgo.cc:431
cms::DDNamespace::addVolumeNS
dd4hep::Volume addVolumeNS(dd4hep::Volume vol) const
Definition: DDNamespace.cc:190
cms::s_executed
static constexpr long s_executed
Definition: DDAlgoArguments.h:16
HGCalEEAlgo::nCutRadPar_
int nCutRadPar_
Definition: DDHGCalEEAlgo.cc:52
alignCSCRings.r
r
Definition: alignCSCRings.py:93
callgraph.module
module
Definition: callgraph.py:61
align::Detector
Definition: StructureType.h:92
HGCalGeomTools.h
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
HGCalEEAlgo::slopeB_
std::vector< double > slopeB_
Definition: DDHGCalEEAlgo.cc:57
HGCalEEAlgo::ConstructAlgo
void ConstructAlgo(cms::DDParsingContext &ctxt, xml_h e)
Definition: DDHGCalEEAlgo.cc:190
HGCalEEAlgo::waferType_
std::unique_ptr< HGCalWaferType > waferType_
Definition: DDHGCalEEAlgo.cc:31
HGCalEEAlgo::absorbMode_
int absorbMode_
Definition: DDHGCalEEAlgo.cc:45
diffTwoXMLs.r1
r1
Definition: diffTwoXMLs.py:53
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
HGCalGeometryMode::Polyhedra
Definition: HGCalGeometryMode.h:37
HGCalEEAlgo::rMaxFront_
std::vector< double > rMaxFront_
Definition: DDHGCalEEAlgo.cc:62
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
TCMET_cfi.corner
corner
Definition: TCMET_cfi.py:34
HGCalEEAlgo::layerCenter_
std::vector< int > layerCenter_
Definition: DDHGCalEEAlgo.cc:43
HGCalEEAlgo::ConstructLayers
void ConstructLayers(const dd4hep::Volume module, cms::DDParsingContext &ctxt, xml_h e)
Definition: DDHGCalEEAlgo.cc:208
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
angle_units.h
HGCalTypes.h
HGCalEEAlgo::zMinRadPar_
double zMinRadPar_
Definition: DDHGCalEEAlgo.cc:50
HGCalEEAlgo::mother_
dd4hep::Volume mother_
Definition: DDHGCalEEAlgo.cc:32
edm::Log
Definition: MessageLogger.h:70
cms::DDNamespace::name
std::string_view name() const
Definition: DDNamespace.h:79
dttmaxenums::R
Definition: DTTMax.h:29
HGCalParameters.h
HGCalEEAlgo::fracAreaMin_
double fracAreaMin_
Definition: DDHGCalEEAlgo.cc:53
cuy.ii
ii
Definition: cuy.py:589
cms::DDNamespace::volume
dd4hep::Volume volume(const std::string &name, bool exc=true) const
Definition: DDNamespace.cc:264
HGCalEEAlgo::layers_
std::vector< int > layers_
Definition: DDHGCalEEAlgo.cc:39
HGCalEEAlgo::PositionSensitive
void PositionSensitive(cms::DDParsingContext &ctxt, xml_h e, const dd4hep::Volume &glog, double rin, double rout, double zpos, int layertype, int layercenter)
Definition: DDHGCalEEAlgo.cc:342
DDutils.h
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
Generated for CMSSW Reference Manual by doxygen 1.8.16