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
HGCalEEFileAlgo Struct Reference

Public Member Functions

 HGCalEEFileAlgo ()
 
 HGCalEEFileAlgo (cms::DDParsingContext &ctxt, xml_h e)
 
void positionSensitive (cms::DDParsingContext &ctxt, xml_h e, const dd4hep::Volume &glog, double rin, double rout, double zpos, int layertype, int layer)
 

Public Attributes

int absorbMode_
 
double alpha_
 
std::unordered_set< int > copies_
 
std::vector< int > copyNumber_
 
double cosAlpha_
 
int firstLayer_
 
HGCalGeomTools geomTools_
 
std::vector< int > layerCenter_
 
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_
 
std::vector< double > rMaxFront_
 
std::vector< double > rMinFront_
 
int sectors_
 
int sensitiveMode_
 
std::vector< double > slopeB_
 
std::vector< double > slopeT_
 
std::vector< double > thick_
 
std::vector< int > waferIndex_
 
std::vector< int > waferProperty_
 
std::vector< std::string > wafers_
 
double waferSepar_
 
double waferSize_
 
std::vector< double > zFrontB_
 
std::vector< double > zFrontT_
 
double zMinBlock_
 

Detailed Description

Definition at line 27 of file DDHGCalEEFileAlgo.cc.

Constructor & Destructor Documentation

HGCalEEFileAlgo::HGCalEEFileAlgo ( )
inline

Definition at line 28 of file DDHGCalEEFileAlgo.cc.

References Exception.

28 { throw cms::Exception("HGCalGeom") << "Wrong initialization to HGCalEEFileAlgo"; }
HGCalEEFileAlgo::HGCalEEFileAlgo ( cms::DDParsingContext &  ctxt,
xml_h  e 
)
inline

Definition at line 29 of file DDHGCalEEFileAlgo.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(), mergeVDriftHistosByStation::name, cms::DDNamespace::name(), cms::DDAlgoArguments::parentName(), HGCalGeometryMode::Polyhedra, PixelTestBeamValidation_cfi::Position, cms::DDNamespace::prepend(), diffTwoXMLs::r1, HGCalGeomTools::radius(), AlCaHLTBitMon_QueryRunRegistry::string, cms::DDAlgoArguments::value(), cms::DDNamespace::volume(), HGCalWaferIndex::waferLayer(), HGCalProperty::waferOrient(), HGCalProperty::waferPartial(), HGCalProperty::waferThick(), HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

29  {
30  cms::DDNamespace ns(ctxt, e, true);
31  cms::DDAlgoArguments args(ctxt, e);
32 
33 #ifdef EDM_ML_DEBUG
34  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: Creating an instance";
35 #endif
36  static constexpr double tol1 = 0.01 * dd4hep::mm;
37  static constexpr double tol2 = 0.00001 * dd4hep::mm;
38 
39  dd4hep::Volume mother = ns.volume(args.parentName());
40  wafers_ = args.value<std::vector<std::string>>("WaferNames");
41 #ifdef EDM_ML_DEBUG
42  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: " << wafers_.size() << " wafers";
43  for (unsigned int i = 0; i < wafers_.size(); ++i)
44  edm::LogVerbatim("HGCalGeom") << "Wafer[" << i << "] " << wafers_[i];
45 #endif
46  materials_ = args.value<std::vector<std::string>>("MaterialNames");
47  names_ = args.value<std::vector<std::string>>("VolumeNames");
48  thick_ = args.value<std::vector<double>>("Thickness");
49  copyNumber_.resize(materials_.size(), 1);
50 #ifdef EDM_ML_DEBUG
51  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: " << materials_.size() << " types of volumes";
52  for (unsigned int i = 0; i < names_.size(); ++i)
53  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << names_[i] << " of thickness "
54  << cms::convert2mm(thick_[i]) << " filled with " << materials_[i]
55  << " first copy number " << copyNumber_[i];
56 #endif
57  layers_ = args.value<std::vector<int>>("Layers");
58  layerThick_ = args.value<std::vector<double>>("LayerThick");
59 #ifdef EDM_ML_DEBUG
60  edm::LogVerbatim("HGCalGeom") << "There are " << layers_.size() << " blocks";
61  for (unsigned int i = 0; i < layers_.size(); ++i)
62  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] of thickness " << cms::convert2mm(layerThick_[i])
63  << " with " << layers_[i] << " layers";
64 #endif
65  layerType_ = args.value<std::vector<int>>("LayerType");
66  layerSense_ = args.value<std::vector<int>>("LayerSense");
67  firstLayer_ = args.value<int>("FirstLayer");
68  absorbMode_ = args.value<int>("AbsorberMode");
69  sensitiveMode_ = args.value<int>("SensitiveMode");
70 #ifdef EDM_ML_DEBUG
71  edm::LogVerbatim("HGCalGeom") << "First Layer " << firstLayer_ << " and "
72  << "Absober:Sensitive mode " << absorbMode_ << ":" << sensitiveMode_;
73 #endif
74  layerCenter_ = args.value<std::vector<int>>("LayerCenter");
75 #ifdef EDM_ML_DEBUG
76  for (unsigned int i = 0; i < layerCenter_.size(); ++i)
77  edm::LogVerbatim("HGCalGeom") << "LayerCenter [" << i << "] " << layerCenter_[i];
78 #endif
79  if (firstLayer_ > 0) {
80  for (unsigned int i = 0; i < layerType_.size(); ++i) {
81  if (layerSense_[i] > 0) {
82  int ii = layerType_[i];
83  copyNumber_[ii] = firstLayer_;
84 #ifdef EDM_ML_DEBUG
85  edm::LogVerbatim("HGCalGeom") << "First copy number for layer type " << i << ":" << ii << " with "
86  << materials_[ii] << " changed to " << copyNumber_[ii];
87 #endif
88  break;
89  }
90  }
91  } else {
92  firstLayer_ = 1;
93  }
94 #ifdef EDM_ML_DEBUG
95  edm::LogVerbatim("HGCalGeom") << "There are " << layerType_.size() << " layers";
96  for (unsigned int i = 0; i < layerType_.size(); ++i)
97  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerType_[i] << " sensitive class "
98  << layerSense_[i];
99 #endif
100  zMinBlock_ = args.value<double>("zMinBlock");
101  waferSize_ = args.value<double>("waferSize");
102  waferSepar_ = args.value<double>("SensorSeparation");
103  sectors_ = args.value<int>("Sectors");
104  alpha_ = (1._pi) / sectors_;
105  cosAlpha_ = cos(alpha_);
106 #ifdef EDM_ML_DEBUG
107  edm::LogVerbatim("HGCalGeom") << "zStart " << cms::convert2mm(zMinBlock_) << " wafer width "
108  << cms::convert2mm(waferSize_) << " separations " << cms::convert2mm(waferSepar_)
109  << " sectors " << sectors_ << ":" << convertRadToDeg(alpha_) << ":" << cosAlpha_;
110 #endif
111  waferIndex_ = args.value<std::vector<int>>("WaferIndex");
112  waferProperty_ = args.value<std::vector<int>>("WaferProperties");
113 #ifdef EDM_ML_DEBUG
114  edm::LogVerbatim("HGCalGeom") << "waferProperties with " << waferIndex_.size() << " entries";
115  for (unsigned int k = 0; k < waferIndex_.size(); ++k)
116  edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << waferIndex_[k] << " ("
117  << HGCalWaferIndex::waferLayer(waferIndex_[k]) << ", "
118  << HGCalWaferIndex::waferU(waferIndex_[k]) << ", "
119  << HGCalWaferIndex::waferV(waferIndex_[k]) << ") : ("
120  << HGCalProperty::waferThick(waferProperty_[k]) << ":"
121  << HGCalProperty::waferPartial(waferProperty_[k]) << ":"
122  << HGCalProperty::waferOrient(waferProperty_[k]) << ")";
123 #endif
124  slopeB_ = args.value<std::vector<double>>("SlopeBottom");
125  zFrontB_ = args.value<std::vector<double>>("ZFrontBottom");
126  rMinFront_ = args.value<std::vector<double>>("RMinFront");
127  slopeT_ = args.value<std::vector<double>>("SlopeTop");
128  zFrontT_ = args.value<std::vector<double>>("ZFrontTop");
129  rMaxFront_ = args.value<std::vector<double>>("RMaxFront");
130 #ifdef EDM_ML_DEBUG
131  for (unsigned int i = 0; i < slopeB_.size(); ++i)
132  edm::LogVerbatim("HGCalGeom") << "Bottom Block [" << i << "] Zmin " << cms::convert2mm(zFrontB_[i]) << " Rmin "
133  << cms::convert2mm(rMinFront_[i]) << " Slope " << slopeB_[i];
134  for (unsigned int i = 0; i < slopeT_.size(); ++i)
135  edm::LogVerbatim("HGCalGeom") << "Top Block [" << i << "] Zmin " << cms::convert2mm(zFrontT_[i]) << " Rmax "
136  << cms::convert2mm(rMaxFront_[i]) << " Slope " << slopeT_[i];
137  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: NameSpace " << ns.name();
138 #endif
139 
140 #ifdef EDM_ML_DEBUG
141  edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalEEFileAlgo...";
142  copies_.clear();
143 #endif
144 
145  double zi(zMinBlock_);
146  int laymin(0);
147  for (unsigned int i = 0; i < layers_.size(); ++i) {
148  double zo = zi + layerThick_[i];
149  double routF = HGCalGeomTools::radius(zi, zFrontT_, rMaxFront_, slopeT_);
150  int laymax = laymin + layers_[i];
151  double zz = zi;
152  double thickTot(0);
153  for (int ly = laymin; ly < laymax; ++ly) {
154  int ii = layerType_[ly];
155  int copy = copyNumber_[ii];
156  double hthick = 0.5 * thick_[ii];
157  double rinB = HGCalGeomTools::radius(zo - tol1, zFrontB_, rMinFront_, slopeB_);
158  zz += hthick;
159  thickTot += thick_[ii];
160 
161  std::string name = names_[ii] + std::to_string(copy);
162 #ifdef EDM_ML_DEBUG
163  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: Layer " << ly << ":" << ii << " Front "
164  << cms::convert2mm(zi) << ", " << cms::convert2mm(routF) << " Back "
165  << cms::convert2mm(zo) << ", " << cms::convert2mm(rinB)
166  << " superlayer thickness " << cms::convert2mm(layerThick_[i]);
167 #endif
168 
169  dd4hep::Material matter = ns.material(materials_[ii]);
170  dd4hep::Volume glog;
171 
172  if (layerSense_[ly] < 1) {
173  std::vector<double> pgonZ, pgonRin, pgonRout;
174  double rmax = routF * cosAlpha_ - tol1;
175  HGCalGeomTools::radius(zz - hthick,
176  zz + hthick,
177  zFrontB_,
178  rMinFront_,
179  slopeB_,
180  zFrontT_,
181  rMaxFront_,
182  slopeT_,
183  -layerSense_[ly],
184  pgonZ,
185  pgonRin,
186  pgonRout);
187  for (unsigned int isec = 0; isec < pgonZ.size(); ++isec) {
188  pgonZ[isec] -= zz;
189  if (layerSense_[ly] == 0 || absorbMode_ == 0)
190  pgonRout[isec] = rmax;
191  else
192  pgonRout[isec] = pgonRout[isec] * cosAlpha_ - tol1;
193  }
194  dd4hep::Solid solid = dd4hep::Polyhedra(sectors_, -alpha_, 2._pi, pgonZ, pgonRin, pgonRout);
195  ns.addSolidNS(ns.prepend(name), solid);
196  glog = dd4hep::Volume(solid.name(), solid, matter);
197  ns.addVolumeNS(glog);
198 #ifdef EDM_ML_DEBUG
199  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: " << solid.name() << " polyhedra of " << sectors_
200  << " sectors covering " << convertRadToDeg(-alpha_) << ":"
201  << convertRadToDeg(-alpha_ + 2._pi) << " with " << pgonZ.size()
202  << " sections and filled with " << matter.name();
203  for (unsigned int k = 0; k < pgonZ.size(); ++k)
204  edm::LogVerbatim("HGCalGeom") << "[" << k << "] z " << cms::convert2mm(pgonZ[k]) << " R "
205  << cms::convert2mm(pgonRin[k]) << ":" << cms::convert2mm(pgonRout[k]);
206 #endif
207  } else {
208  double rins =
209  (sensitiveMode_ < 1) ? rinB : HGCalGeomTools::radius(zz + hthick - tol1, zFrontB_, rMinFront_, slopeB_);
210  double routs =
211  (sensitiveMode_ < 1) ? routF : HGCalGeomTools::radius(zz - hthick, zFrontT_, rMaxFront_, slopeT_);
212  dd4hep::Solid solid = dd4hep::Tube(rins, routs, hthick, 0.0, 2._pi);
213  ns.addSolidNS(ns.prepend(name), solid);
214  glog = dd4hep::Volume(solid.name(), solid, matter);
215  ns.addVolumeNS(glog);
216 
217 #ifdef EDM_ML_DEBUG
218  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: " << solid.name() << " Tubs made of " << matter.name()
219  << " of dimensions " << cms::convert2mm(rinB) << ":" << cms::convert2mm(rins)
220  << ", " << cms::convert2mm(routF) << ":" << cms::convert2mm(routs) << ", "
221  << cms::convert2mm(hthick) << ", 0.0, 360.0 and position " << glog.name()
222  << " number " << copy << ":" << layerCenter_[copy - firstLayer_];
223 #endif
224  positionSensitive(ctxt, e, glog, rins, routs, zz, layerSense_[ly], (copy - firstLayer_));
225  }
226 
227  dd4hep::Position r1(0, 0, zz);
228  mother.placeVolume(glog, copy, r1);
229  ++copyNumber_[ii];
230 #ifdef EDM_ML_DEBUG
231  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: " << glog.name() << " number " << copy << " positioned in "
232  << mother.name() << " at (0,0," << cms::convert2mm(zz) << ") with no rotation";
233 #endif
234  zz += hthick;
235  } // End of loop over layers in a block
236  zi = zo;
237  laymin = laymax;
238  // Make consistency check of all the partitions of the block
239  if (std::abs(thickTot - layerThick_[i]) >= tol2) {
240  if (thickTot > layerThick_[i]) {
241  edm::LogError("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
242  << " is smaller than " << cms::convert2mm(thickTot)
243  << ": thickness of all its components **** ERROR ****";
244  } else {
245  edm::LogWarning("HGCalGeom") << "Thickness of the partition " << cms::convert2mm(layerThick_[i])
246  << " does not match with " << cms::convert2mm(thickTot) << " of the components";
247  }
248  }
249  } // End of loop over blocks
250 
251 #ifdef EDM_ML_DEBUG
252  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: " << copies_.size() << " different wafer copy numbers";
253  int k(0);
254  for (std::unordered_set<int>::const_iterator itr = copies_.begin(); itr != copies_.end(); ++itr, ++k) {
255  edm::LogVerbatim("HGCalGeom") << "Copy [" << k << "] : " << (*itr);
256  }
257  copies_.clear();
258  edm::LogVerbatim("HGCalGeom") << "<<== End of DDHGCalEEFileAlgo construction...";
259 #endif
260  }
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
HGCalEEFileAlgo::positionSensitive
void positionSensitive(cms::DDParsingContext &ctxt, xml_h e, const dd4hep::Volume &glog, double rin, double rout, double zpos, int layertype, int layer)
Definition: DDHGCalEEFileAlgo.cc:262
HGCalWaferIndex::waferU
int32_t waferU(const int32_t index)
Definition: HGCalWaferIndex.cc:27
HGCalEEFileAlgo::thick_
std::vector< double > thick_
Definition: DDHGCalEEFileAlgo.cc:363
HGCalEEFileAlgo::layerThick_
std::vector< double > layerThick_
Definition: DDHGCalEEFileAlgo.cc:366
HGCalWaferIndex::waferLayer
int32_t waferLayer(const int32_t index)
Definition: HGCalWaferIndex.cc:23
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
HGCalEEFileAlgo::names_
std::vector< std::string > names_
Definition: DDHGCalEEFileAlgo.cc:362
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
HGCalEEFileAlgo::waferProperty_
std::vector< int > waferProperty_
Definition: DDHGCalEEFileAlgo.cc:375
HGCalProperty::waferOrient
int32_t waferOrient(const int32_t property)
Definition: HGCalProperty.cc:17
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HGCalEEFileAlgo::slopeT_
std::vector< double > slopeT_
Definition: DDHGCalEEFileAlgo.cc:382
HGCalEEFileAlgo::layerCenter_
std::vector< int > layerCenter_
Definition: DDHGCalEEFileAlgo.cc:369
HGCalEEFileAlgo::layerSense_
std::vector< int > layerSense_
Definition: DDHGCalEEFileAlgo.cc:368
HGCalEEFileAlgo::copies_
std::unordered_set< int > copies_
Definition: DDHGCalEEFileAlgo.cc:385
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
HGCalProperty::waferThick
int32_t waferThick(const int32_t property)
Definition: HGCalProperty.cc:9
cms::Volume
dd4hep::Volume Volume
Definition: DDFilteredView.h:47
HGCalEEFileAlgo::copyNumber_
std::vector< int > copyNumber_
Definition: DDHGCalEEFileAlgo.cc:364
HGCalEEFileAlgo::slopeB_
std::vector< double > slopeB_
Definition: DDHGCalEEFileAlgo.cc:379
HGCalProperty::waferPartial
int32_t waferPartial(const int32_t property)
Definition: HGCalProperty.cc:13
HGCalEEFileAlgo::zFrontT_
std::vector< double > zFrontT_
Definition: DDHGCalEEFileAlgo.cc:383
HGCalEEFileAlgo::layerType_
std::vector< int > layerType_
Definition: DDHGCalEEFileAlgo.cc:367
HGCalEEFileAlgo::waferIndex_
std::vector< int > waferIndex_
Definition: DDHGCalEEFileAlgo.cc:374
HGCalEEFileAlgo::rMinFront_
std::vector< double > rMinFront_
Definition: DDHGCalEEFileAlgo.cc:381
HGCalWaferIndex::waferV
int32_t waferV(const int32_t index)
Definition: HGCalWaferIndex.cc:32
HGCalEEFileAlgo::rMaxFront_
std::vector< double > rMaxFront_
Definition: DDHGCalEEFileAlgo.cc:384
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
HGCalEEFileAlgo::zFrontB_
std::vector< double > zFrontB_
Definition: DDHGCalEEFileAlgo.cc:380
HGCalEEFileAlgo::materials_
std::vector< std::string > materials_
Definition: DDHGCalEEFileAlgo.cc:361
HGCalEEFileAlgo::layers_
std::vector< int > layers_
Definition: DDHGCalEEFileAlgo.cc:365
HGCalEEFileAlgo::wafers_
std::vector< std::string > wafers_
Definition: DDHGCalEEFileAlgo.cc:360

Member Function Documentation

void HGCalEEFileAlgo::positionSensitive ( cms::DDParsingContext &  ctxt,
xml_h  e,
const dd4hep::Volume &  glog,
double  rin,
double  rout,
double  zpos,
int  layertype,
int  layer 
)
inline

Definition at line 262 of file DDHGCalEEFileAlgo.cc.

References funct::abs(), cms::convert2mm(), filterCSVwithJSON::copy, PVValHelper::dy, HGCalWaferType::getType(), gpuVertexFinder::iv, HGCalParameters::k_CornerSize, phase1PixelTopology::layer, N, nin, gpuPixelDoublets::ntot, HGCalTypes::packTypeUV(), PixelTestBeamValidation_cfi::Position, dttmaxenums::R, alignCSCRings::r, mathSSE::sqrt(), findQualityFiles::v, cms::DDNamespace::volume(), HGCalGeomTools::waferCorner(), and HGCalWaferIndex::waferIndex().

269  {
270  cms::DDNamespace ns(ctxt, e, true);
271  static const double sqrt3 = std::sqrt(3.0);
272  int layercenter = layerCenter_[layer];
273  double r = 0.5 * (waferSize_ + waferSepar_);
274  double R = 2.0 * r / sqrt3;
275  double dy = 0.75 * R;
276  int N = (int)(0.5 * rout / r) + 2;
277  const auto& xyoff = geomTools_.shiftXY(layercenter, (waferSize_ + waferSepar_));
278 #ifdef EDM_ML_DEBUG
279  int ium(0), ivm(0), iumAll(0), ivmAll(0), kount(0), ntot(0), nin(0);
280  std::vector<int> ntype(6, 0);
281  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: " << glog.name() << " rin:rout " << cms::convert2mm(rin) << ":"
282  << cms::convert2mm(rout) << " zpos " << cms::convert2mm(zpos) << " N " << N
283  << " for maximum u, v; r " << cms::convert2mm(r) << " R " << cms::convert2mm(R)
284  << " dy " << cms::convert2mm(dy) << " Shift " << cms::convert2mm(xyoff.first) << ":"
285  << cms::convert2mm(xyoff.second) << " WaferSize "
286  << cms::convert2mm((waferSize_ + waferSepar_));
287 #endif
288  for (int u = -N; u <= N; ++u) {
289  for (int v = -N; v <= N; ++v) {
290 #ifdef EDM_ML_DEBUG
291  int iu = std::abs(u);
292  int iv = std::abs(v);
293 #endif
294  int nr = 2 * v;
295  int nc = -2 * u + v;
296  double xpos = xyoff.first + nc * r;
297  double ypos = xyoff.second + nr * dy;
298  const auto& corner = HGCalGeomTools::waferCorner(xpos, ypos, r, R, rin, rout, false);
299 #ifdef EDM_ML_DEBUG
300  ++ntot;
301  if (((corner.first <= 0) && std::abs(u) < 5 && std::abs(v) < 5) || (std::abs(u) < 2 && std::abs(v) < 2)) {
302  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: " << glog.name() << " R " << cms::convert2mm(rin) << ":"
303  << cms::convert2mm(rout) << "\n Z " << cms::convert2mm(zpos) << " LayerType "
304  << layertype << " u " << u << " v " << v << " with " << corner.first
305  << " corners";
306  }
307 #endif
308  int indx = HGCalWaferIndex::waferIndex((layer + firstLayer_), u, v, false);
309  int type = HGCalWaferType::getType(indx, waferIndex_, waferProperty_);
310  if (corner.first > 0 && type >= 0) {
311  int copy = HGCalTypes::packTypeUV(type, u, v);
312 #ifdef EDM_ML_DEBUG
313  edm::LogVerbatim("HGCalGeom") << " DDHGCalEEFileAlgo: " << wafers_[type] << " number " << copy << " type "
314  << type << " layer:u:v:indx " << (layer + firstLayer_) << ":" << u << ":" << v
315  << ":" << indx;
316  if (iu > ium)
317  ium = iu;
318  if (iv > ivm)
319  ivm = iv;
320  kount++;
321  if (copies_.count(copy) == 0)
322  copies_.insert(copy);
323 #endif
324  if (corner.first == (int)(HGCalParameters::k_CornerSize)) {
325 #ifdef EDM_ML_DEBUG
326  if (iu > iumAll)
327  iumAll = iu;
328  if (iv > ivmAll)
329  ivmAll = iv;
330  ++nin;
331 #endif
332  dd4hep::Position tran(xpos, ypos, 0.0);
333  if (layertype > 1)
334  type += 3;
335  glog.placeVolume(ns.volume(wafers_[type]), copy, tran);
336 #ifdef EDM_ML_DEBUG
337  ++ntype[type];
338  edm::LogVerbatim("HGCalGeom")
339  << " DDHGCalEEFileAlgo: " << wafers_[type] << " number " << copy << " type " << layertype << ":" << type
340  << " positioned in " << glog.name() << " at (" << cms::convert2mm(xpos) << ", " << cms::convert2mm(ypos)
341  << ", 0) with no rotation";
342 #endif
343  }
344  }
345  }
346  }
347 
348 #ifdef EDM_ML_DEBUG
349  edm::LogVerbatim("HGCalGeom") << "DDHGCalEEFileAlgo: Maximum # of u " << ium << ":" << iumAll << " # of v " << ivm
350  << ":" << ivmAll << " and " << nin << ":" << kount << ":" << ntot << " wafers ("
351  << ntype[0] << ":" << ntype[1] << ":" << ntype[2] << ":" << ntype[3] << ":"
352  << ntype[4] << ":" << ntype[5] << ") for " << glog.name() << " R "
353  << cms::convert2mm(rin) << ":" << cms::convert2mm(rout);
354 #endif
355  }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
HGCalEEFileAlgo::geomTools_
HGCalGeomTools geomTools_
Definition: DDHGCalEEFileAlgo.cc:358
gpuVertexFinder::iv
int32_t *__restrict__ iv
Definition: gpuClusterTracksDBSCAN.h:42
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
HGCalEEFileAlgo::waferProperty_
std::vector< int > waferProperty_
Definition: DDHGCalEEFileAlgo.cc:375
HGCalParameters::k_CornerSize
static constexpr uint32_t k_CornerSize
Definition: HGCalParameters.h:44
phase1PixelTopology::layer
constexpr std::array< uint8_t, layerIndexSize > layer
Definition: phase1PixelTopology.h:110
nin
int nin
Definition: CascadeWrapper.h:114
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
HGCalWaferType::getType
int getType(double xpos, double ypos, double zpos)
Definition: HGCalWaferType.cc:36
HGCalEEFileAlgo::layerCenter_
std::vector< int > layerCenter_
Definition: DDHGCalEEFileAlgo.cc:369
HGCalEEFileAlgo::copies_
std::unordered_set< int > copies_
Definition: DDHGCalEEFileAlgo.cc:385
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
HGCalWaferIndex::waferIndex
int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
Definition: HGCalWaferIndex.cc:4
N
#define N
Definition: blowfish.cc:9
HGCalEEFileAlgo::waferIndex_
std::vector< int > waferIndex_
Definition: DDHGCalEEFileAlgo.cc:374
HGCalTypes::packTypeUV
static int32_t packTypeUV(int type, int u, int v)
Definition: HGCalTypes.cc:3
gpuPixelDoublets::ntot
__shared__ uint32_t ntot
Definition: gpuPixelDoubletsAlgos.h:67
HGCalEEFileAlgo::wafers_
std::vector< std::string > wafers_
Definition: DDHGCalEEFileAlgo.cc:360

Member Data Documentation

int HGCalEEFileAlgo::absorbMode_

Definition at line 371 of file DDHGCalEEFileAlgo.cc.

double HGCalEEFileAlgo::alpha_

Definition at line 386 of file DDHGCalEEFileAlgo.cc.

std::unordered_set<int> HGCalEEFileAlgo::copies_

Definition at line 385 of file DDHGCalEEFileAlgo.cc.

std::vector<int> HGCalEEFileAlgo::copyNumber_

Definition at line 364 of file DDHGCalEEFileAlgo.cc.

double HGCalEEFileAlgo::cosAlpha_

Definition at line 386 of file DDHGCalEEFileAlgo.cc.

int HGCalEEFileAlgo::firstLayer_

Definition at line 370 of file DDHGCalEEFileAlgo.cc.

HGCalGeomTools HGCalEEFileAlgo::geomTools_

Definition at line 358 of file DDHGCalEEFileAlgo.cc.

std::vector<int> HGCalEEFileAlgo::layerCenter_

Definition at line 369 of file DDHGCalEEFileAlgo.cc.

std::vector<int> HGCalEEFileAlgo::layers_

Definition at line 365 of file DDHGCalEEFileAlgo.cc.

std::vector<int> HGCalEEFileAlgo::layerSense_

Definition at line 368 of file DDHGCalEEFileAlgo.cc.

std::vector<double> HGCalEEFileAlgo::layerThick_

Definition at line 366 of file DDHGCalEEFileAlgo.cc.

std::vector<int> HGCalEEFileAlgo::layerType_

Definition at line 367 of file DDHGCalEEFileAlgo.cc.

std::vector<std::string> HGCalEEFileAlgo::materials_

Definition at line 361 of file DDHGCalEEFileAlgo.cc.

std::vector<std::string> HGCalEEFileAlgo::names_

Definition at line 362 of file DDHGCalEEFileAlgo.cc.

std::vector<double> HGCalEEFileAlgo::rMaxFront_

Definition at line 384 of file DDHGCalEEFileAlgo.cc.

std::vector<double> HGCalEEFileAlgo::rMinFront_

Definition at line 381 of file DDHGCalEEFileAlgo.cc.

int HGCalEEFileAlgo::sectors_

Definition at line 378 of file DDHGCalEEFileAlgo.cc.

int HGCalEEFileAlgo::sensitiveMode_

Definition at line 372 of file DDHGCalEEFileAlgo.cc.

std::vector<double> HGCalEEFileAlgo::slopeB_

Definition at line 379 of file DDHGCalEEFileAlgo.cc.

std::vector<double> HGCalEEFileAlgo::slopeT_

Definition at line 382 of file DDHGCalEEFileAlgo.cc.

std::vector<double> HGCalEEFileAlgo::thick_

Definition at line 363 of file DDHGCalEEFileAlgo.cc.

std::vector<int> HGCalEEFileAlgo::waferIndex_

Definition at line 374 of file DDHGCalEEFileAlgo.cc.

std::vector<int> HGCalEEFileAlgo::waferProperty_

Definition at line 375 of file DDHGCalEEFileAlgo.cc.

std::vector<std::string> HGCalEEFileAlgo::wafers_

Definition at line 360 of file DDHGCalEEFileAlgo.cc.

double HGCalEEFileAlgo::waferSepar_

Definition at line 377 of file DDHGCalEEFileAlgo.cc.

double HGCalEEFileAlgo::waferSize_

Definition at line 376 of file DDHGCalEEFileAlgo.cc.

std::vector<double> HGCalEEFileAlgo::zFrontB_

Definition at line 380 of file DDHGCalEEFileAlgo.cc.

std::vector<double> HGCalEEFileAlgo::zFrontT_

Definition at line 383 of file DDHGCalEEFileAlgo.cc.

double HGCalEEFileAlgo::zMinBlock_

Definition at line 373 of file DDHGCalEEFileAlgo.cc.

Generated for CMSSW Reference Manual by doxygen 1.8.5