CMS 3D CMS Logo
List of all members | Public Member Functions | Protected Member Functions | Private Attributes
DDHGCalModuleAlgo Class Reference

#include <DDHGCalModuleAlgo.h>

Inheritance diagram for DDHGCalModuleAlgo:

Public Member Functions

 DDHGCalModuleAlgo ()
 
void execute (DDCompactView &cpv) override
 
void initialize (const DDNumericArguments &nArgs, const DDVectorArguments &vArgs, const DDMapArguments &mArgs, const DDStringArguments &sArgs, const DDStringVectorArguments &vsArgs) override
 
 ~DDHGCalModuleAlgo () override
 

Protected Member Functions

void constructLayers (const DDLogicalPart &, DDCompactView &cpv)
 
void positionSensitive (DDLogicalPart &glog, double rin, double rout, DDCompactView &cpv)
 
double rMax (double z)
 

Private Attributes

std::unordered_set< int > copies
 
std::vector< int > copyNumber
 
std::string idName
 
std::string idNameSpace
 
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
 
double rMaxFine
 
std::vector< double > rMaxFront
 
int sectors
 
std::vector< double > slopeB
 
std::vector< double > slopeT
 
std::vector< double > thick
 
std::vector< std::string > wafer
 
double waferGap
 
double waferW
 
std::vector< double > zFront
 
double zMinBlock
 

Detailed Description

Definition at line 11 of file DDHGCalModuleAlgo.h.

Constructor & Destructor Documentation

DDHGCalModuleAlgo::DDHGCalModuleAlgo ( )

Definition at line 24 of file DDHGCalModuleAlgo.cc.

24  {
25 #ifdef EDM_ML_DEBUG
26  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: Creating an instance";
27 #endif
28 }
DDHGCalModuleAlgo::~DDHGCalModuleAlgo ( )
override

Definition at line 30 of file DDHGCalModuleAlgo.cc.

30 {}

Member Function Documentation

void DDHGCalModuleAlgo::constructLayers ( const DDLogicalPart &  module,
DDCompactView &  cpv 
)
protected

Definition at line 128 of file DDHGCalModuleAlgo.cc.

References alpha, geant_units::operators::convertRadToDeg(), popcon2dropbox::copy(), funct::cos(), DDBase< N, C >::ddname(), DDSplit(), plotBeamSpotDB::first, mps_fire::i, cuy::ii, gen::k, LayerTriplets::layers(), dataset::name, DDBase< N, C >::name(), names, geant_units::piRadians(), DDSolidFactory::polyhedra(), DDCompactView::position(), diffTwoXMLs::r1, hfnoseDigiStudy_cfi::rMax, makeMuonMisalignmentScenario::rot, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, DDSolidFactory::tubs(), and geometryCSVtoXML::zz.

129  {
130 #ifdef EDM_ML_DEBUG
131  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo test: \t\tInside Layers";
132 #endif
133  double zi(zMinBlock);
134  int laymin(0);
135  const double tol(0.01);
136  for (unsigned int i = 0; i < layers.size(); i++) {
137  double zo = zi + layerThick[i];
138  double routF = rMax(zi);
139  int laymax = laymin + layers[i];
140  double zz = zi;
141  double thickTot(0);
142  for (int ly = laymin; ly < laymax; ++ly) {
143  int ii = layerType[ly];
144  int copy = copyNumber[ii];
145  double rinB = (layerSense[ly] == 0) ? (zo * slopeB[0]) : (zo * slopeB[1]);
146  zz += (0.5 * thick[ii]);
147  thickTot += thick[ii];
148 
149  std::string name = "HGCal" + names[ii] + std::to_string(copy);
150 #ifdef EDM_ML_DEBUG
151  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo test: Layer " << ly
152  << ":" << ii << " Front " << zi << ", "
153  << routF << " Back " << zo << ", " << rinB
154  << " superlayer thickness " << layerThick[i];
155 #endif
156  DDName matName(DDSplit(materials[ii]).first,
157  DDSplit(materials[ii]).second);
158  DDMaterial matter(matName);
159  DDLogicalPart glog;
160  if (layerSense[ly] == 0) {
161  double alpha = geant_units::piRadians / sectors;
162  double rmax = routF * cos(alpha) - tol;
163  std::vector<double> pgonZ, pgonRin, pgonRout;
164  pgonZ.emplace_back(-0.5 * thick[ii]);
165  pgonZ.emplace_back(0.5 * thick[ii]);
166  pgonRin.emplace_back(rinB);
167  pgonRin.emplace_back(rinB);
168  pgonRout.emplace_back(rmax);
169  pgonRout.emplace_back(rmax);
170  DDSolid solid = DDSolidFactory::polyhedra(DDName(name, idNameSpace),
171  sectors, -alpha,
172  2*geant_units::piRadians,
173  pgonZ, pgonRin, pgonRout);
174  glog = DDLogicalPart(solid.ddname(), matter, solid);
175 #ifdef EDM_ML_DEBUG
176  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << solid.name()
177  << " polyhedra of " << sectors
178  << " sectors covering "
179  << convertRadToDeg(-alpha) << ":"
180  << (360.0+convertRadToDeg(-alpha))
181  << " with " << pgonZ.size() << " sections";
182  for (unsigned int k = 0; k < pgonZ.size(); ++k)
183  edm::LogVerbatim("HGCalGeom") << "[" << k << "] z " << pgonZ[k] << " R "
184  << pgonRin[k] << ":" << pgonRout[k];
185 #endif
186  } else {
187  DDSolid solid =
188  DDSolidFactory::tubs(DDName(name, idNameSpace), 0.5 * thick[ii],
189  rinB, routF, 0.0, 2*geant_units::piRadians);
190  glog = DDLogicalPart(solid.ddname(), matter, solid);
191 #ifdef EDM_ML_DEBUG
192  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << solid.name()
193  << " Tubs made of " << matName
194  << " of dimensions " << rinB << ", "
195  << routF << ", " << 0.5 * thick[ii]
196  << ", 0.0, 360.0";
197 #endif
198  positionSensitive(glog, rinB, routF, cpv);
199  }
200  DDTranslation r1(0, 0, zz);
201  DDRotation rot;
202  cpv.position(glog, module, copy, r1, rot);
203  ++copyNumber[ii];
204 #ifdef EDM_ML_DEBUG
205  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo test: " << glog.name() << " number "
206  << copy << " positioned in " << module.name() << " at " << r1
207  << " with " << rot << std::endl;
208 #endif
209  zz += (0.5 * thick[ii]);
210  } // End of loop over layers in a block
211  zi = zo;
212  laymin = laymax;
213  if (fabs(thickTot - layerThick[i]) < 0.00001) {
214  } else if (thickTot > layerThick[i]) {
215  edm::LogError("HGCalGeom")
216  << "Thickness of the partition " << layerThick[i]
217  << " is smaller than thickness " << thickTot
218  << " of all its components **** ERROR ****\n";
219  } else if (thickTot < layerThick[i]) {
220  edm::LogWarning("HGCalGeom")
221  << "Thickness of the partition " << layerThick[i]
222  << " does not match with " << thickTot << " of the components\n";
223  }
224  } // End of loop over blocks
225 }
DDHGCalModuleAlgo::copyNumber
std::vector< int > copyNumber
Definition: DDHGCalModuleAlgo.h:34
DDHGCalModuleAlgo::thick
std::vector< double > thick
Definition: DDHGCalModuleAlgo.h:33
alpha
float alpha
Definition: AMPTWrapper.h:95
mps_fire.i
i
Definition: mps_fire.py:338
DDBase::name
const N & name() const
Definition: DDBase.h:74
DDHGCalModuleAlgo::positionSensitive
void positionSensitive(DDLogicalPart &glog, double rin, double rout, DDCompactView &cpv)
Definition: DDHGCalModuleAlgo.cc:245
popcon2dropbox.copy
def copy(args, dbName)
Definition: popcon2dropbox.py:122
DDMaterial
DDMaterial is used to define and access material information.
Definition: DDMaterial.h:43
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
DDHGCalModuleAlgo::names
std::vector< std::string > names
Definition: DDHGCalModuleAlgo.h:32
DDName
DDName is used to identify DDD entities uniquely.
Definition: DDName.h:15
geant_units::operators::convertRadToDeg
constexpr NumType convertRadToDeg(NumType radians)
Definition: GeantUnits.h:98
geant_units::piRadians
constexpr long double piRadians(M_PI)
DDSolid
A DDSolid represents the shape of a part.
Definition: DDSolid.h:39
DDTranslation
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DDTranslation
Definition: DDTranslation.h:7
DDRotation
Represents a uniquely identifyable rotation matrix.
Definition: DDTransform.h:68
DDHGCalModuleAlgo::materials
std::vector< std::string > materials
Definition: DDHGCalModuleAlgo.h:31
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:215
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
DDLogicalPart
A DDLogicalPart aggregates information concerning material, solid and sensitveness ...
Definition: DDLogicalPart.h:93
DDSolidFactory::tubs
static DDSolid tubs(const DDName &name, double zhalf, double rIn, double rOut, double startPhi, double deltaPhi)
Definition: DDSolid.cc:865
cuy.ii
ii
Definition: cuy.py:590
gen::k
int k[5][pyjets_maxn]
Definition: Cascade2Hadronizer.cc:79
DDHGCalModuleAlgo::slopeB
std::vector< double > slopeB
Definition: DDHGCalModuleAlgo.h:44
DDHGCalModuleAlgo::idNameSpace
std::string idNameSpace
Definition: DDHGCalModuleAlgo.h:49
DDHGCalModuleAlgo::layerType
std::vector< int > layerType
Definition: DDHGCalModuleAlgo.h:37
DDHGCalModuleAlgo::layers
std::vector< int > layers
Definition: DDHGCalModuleAlgo.h:35
DDCompactView::position
void position(const DDLogicalPart &self, const DDLogicalPart &parent, const std::string &copyno, const DDTranslation &trans, const DDRotation &rot, const DDDivision *div=0)
Definition: DDCompactView.cc:84
DDSplit
std::pair< std::string, std::string > DDSplit(const std::string &n)
split into (name,namespace), separator = &#39;:&#39;
Definition: DDSplit.cc:3
DDSolidFactory::polyhedra
static DDSolid polyhedra(const DDName &name, int sides, double startPhi, double deltaPhi, const std::vector< double > &z, const std::vector< double > &rmin, const std::vector< double > &rmax)
Creates a polyhedra (refere to Geant3 or Geant4 documentation)
Definition: DDSolid.cc:730
DDHGCalModuleAlgo::layerThick
std::vector< double > layerThick
Definition: DDHGCalModuleAlgo.h:36
DDHGCalModuleAlgo::rMax
double rMax(double z)
Definition: DDHGCalModuleAlgo.cc:227
DDHGCalModuleAlgo::layerSense
std::vector< int > layerSense
Definition: DDHGCalModuleAlgo.h:38
DDBase::ddname
const N & ddname() const
Definition: DDBase.h:76
void DDHGCalModuleAlgo::execute ( DDCompactView &  cpv)
override

Definition at line 113 of file DDHGCalModuleAlgo.cc.

References class-composition::parent.

113  {
114 #ifdef EDM_ML_DEBUG
115  edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalModuleAlgo...";
116 #endif
117  copies.clear();
118  constructLayers(parent(), cpv);
119 #ifdef EDM_ML_DEBUG
120  edm::LogVerbatim("HGCalGeom") << copies.size()<<" different wafer copy numbers";
121 #endif
122  copies.clear();
123 #ifdef EDM_ML_DEBUG
124  edm::LogVerbatim("HGCalGeom") << "<<== End of DDHGCalModuleAlgo construction";
125 #endif
126 }
DDHGCalModuleAlgo::copies
std::unordered_set< int > copies
Definition: DDHGCalModuleAlgo.h:50
DDHGCalModuleAlgo::constructLayers
void constructLayers(const DDLogicalPart &, DDCompactView &cpv)
Definition: DDHGCalModuleAlgo.cc:128
void DDHGCalModuleAlgo::initialize ( const DDNumericArguments &  nArgs,
const DDVectorArguments &  vArgs,
const DDMapArguments &  mArgs,
const DDStringArguments &  sArgs,
const DDStringVectorArguments &  vsArgs 
)
override

Definition at line 32 of file DDHGCalModuleAlgo.cc.

References dbl_to_int(), mps_fire::i, createfilelist::int, LayerTriplets::layers(), names, and DDCurrentNamespace::ns().

36  {
37  wafer = vsArgs["WaferName"];
38 #ifdef EDM_ML_DEBUG
39  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << wafer.size()
40  << " wafers";
41  for (unsigned int i = 0; i < wafer.size(); ++i)
42  edm::LogVerbatim("HGCalGeom") << "Wafer[" << i << "] " << wafer[i];
43 #endif
44  materials = vsArgs["MaterialNames"];
45  names = vsArgs["VolumeNames"];
46  thick = vArgs["Thickness"];
47  for (unsigned int i = 0; i < materials.size(); ++i) {
48  copyNumber.emplace_back(1);
49  }
50 #ifdef EDM_ML_DEBUG
51  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << materials.size()
52  << " types of volumes";
53  for (unsigned int i = 0; i < names.size(); ++i)
54  edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << names[i]
55  << " of thickness " << thick[i]
56  << " filled with " << materials[i]
57  << " first copy number " << copyNumber[i];
58 #endif
59  layers = dbl_to_int(vArgs["Layers"]);
60  layerThick = vArgs["LayerThick"];
61 #ifdef EDM_ML_DEBUG
62  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << layers.size()
63  << " blocks";
64  for (unsigned int i = 0; i < layers.size(); ++i)
65  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] of thickness "
66  << layerThick[i] << " with " << layers[i]
67  << " layers";
68 #endif
69  layerType = dbl_to_int(vArgs["LayerType"]);
70  layerSense = dbl_to_int(vArgs["LayerSense"]);
71 #ifdef EDM_ML_DEBUG
72  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << layerType.size()
73  << " layers";
74  for (unsigned int i = 0; i < layerType.size(); ++i)
75  edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type "
76  << layerType[i] << " sensitive class "
77  << layerSense[i];
78 #endif
79  zMinBlock = nArgs["zMinBlock"];
80  rMaxFine = nArgs["rMaxFine"];
81  waferW = nArgs["waferW"];
82  waferGap = nArgs["waferGap"];
83  sectors = (int)(nArgs["Sectors"]);
84 #ifdef EDM_ML_DEBUG
85  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: zStart " << zMinBlock
86  << " rFineCoarse " << rMaxFine << " wafer width "
87  << waferW << " gap among wafers " << waferGap
88  << " sectors " << sectors;
89 #endif
90  slopeB = vArgs["SlopeBottom"];
91  slopeT = vArgs["SlopeTop"];
92  zFront = vArgs["ZFront"];
93  rMaxFront = vArgs["RMaxFront"];
94 #ifdef EDM_ML_DEBUG
95  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: Bottom slopes "
96  << slopeB[0] << ":" << slopeB[1] << " and "
97  << slopeT.size() << " slopes for top";
98  for (unsigned int i = 0; i < slopeT.size(); ++i)
99  edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] Zmin " << zFront[i]
100  << " Rmax " << rMaxFront[i] << " Slope "
101  << slopeT[i];
102 #endif
103  idNameSpace = DDCurrentNamespace::ns();
104 #ifdef EDM_ML_DEBUG
105  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: NameSpace " << idNameSpace;
106 #endif
107 }
DDHGCalModuleAlgo::copyNumber
std::vector< int > copyNumber
Definition: DDHGCalModuleAlgo.h:34
DDHGCalModuleAlgo::thick
std::vector< double > thick
Definition: DDHGCalModuleAlgo.h:33
mps_fire.i
i
Definition: mps_fire.py:338
DDHGCalModuleAlgo::names
std::vector< std::string > names
Definition: DDHGCalModuleAlgo.h:32
DDCurrentNamespace::ns
static std::string & ns()
Definition: DDCurrentNamespace.cc:4
DDHGCalModuleAlgo::materials
std::vector< std::string > materials
Definition: DDHGCalModuleAlgo.h:31
DDHGCalModuleAlgo::rMaxFront
std::vector< double > rMaxFront
Definition: DDHGCalModuleAlgo.h:47
dbl_to_int
std::vector< int > dbl_to_int(const std::vector< double > &vecdbl)
Converts a std::vector of doubles to a std::vector of int.
Definition: DDutils.h:7
DDHGCalModuleAlgo::slopeB
std::vector< double > slopeB
Definition: DDHGCalModuleAlgo.h:44
DDHGCalModuleAlgo::idNameSpace
std::string idNameSpace
Definition: DDHGCalModuleAlgo.h:49
DDHGCalModuleAlgo::layerType
std::vector< int > layerType
Definition: DDHGCalModuleAlgo.h:37
DDHGCalModuleAlgo::layers
std::vector< int > layers
Definition: DDHGCalModuleAlgo.h:35
DDHGCalModuleAlgo::slopeT
std::vector< double > slopeT
Definition: DDHGCalModuleAlgo.h:45
DDHGCalModuleAlgo::wafer
std::vector< std::string > wafer
Definition: DDHGCalModuleAlgo.h:30
DDHGCalModuleAlgo::zFront
std::vector< double > zFront
Definition: DDHGCalModuleAlgo.h:46
DDHGCalModuleAlgo::layerThick
std::vector< double > layerThick
Definition: DDHGCalModuleAlgo.h:36
DDHGCalModuleAlgo::layerSense
std::vector< int > layerSense
Definition: DDHGCalModuleAlgo.h:38
void DDHGCalModuleAlgo::positionSensitive ( DDLogicalPart &  glog,
double  rin,
double  rout,
DDCompactView &  cpv 
)
protected

Definition at line 245 of file DDHGCalModuleAlgo.cc.

References popcon2dropbox::copy(), TCMET_cfi::corner, DDBase< N, C >::ddname(), DDSplit(), PVValHelper::dx, PVValHelper::dy, createfilelist::int, HGCalParameters::k_CornerSize, dataset::name, create_public_lumi_plots::ncol, DDCompactView::position(), idealTransformation::rotation, findQualityFiles::rr, mathSSE::sqrt(), funct::tan(), and HGCalGeomTools::waferCorner().

246  {
247  double ww = (waferW + waferGap);
248  double dx = 0.5 * ww;
249  double dy = 3.0 * dx * tan(30._deg);
250  double rr = 2.0 * dx * tan(30._deg);
251  int ncol = (int)(2.0 * rout / ww) + 1;
252  int nrow = (int)(rout / (ww * tan(30._deg))) + 1;
253  int incm(0), inrm(0), kount(0);
254 #ifdef EDM_ML_DEBUG
255  edm::LogVerbatim("HGCalGeom") << glog.ddname() << " rout " << rout << " Row "
256  << nrow << " Column " << ncol;
257 #endif
258  for (int nr = -nrow; nr <= nrow; ++nr) {
259  int inr = (nr >= 0) ? nr : -nr;
260  for (int nc = -ncol; nc <= ncol; ++nc) {
261  int inc = (nc >= 0) ? nc : -nc;
262  if (inr % 2 == inc % 2) {
263  double xpos = nc * dx;
264  double ypos = nr * dy;
265  std::pair<int, int> corner =
266  HGCalGeomTools::waferCorner(xpos, ypos, dx, rr, rin, rout, true);
267  if (corner.first == (int)(HGCalParameters::k_CornerSize)) {
268  double rpos = std::sqrt(xpos * xpos + ypos * ypos);
269  DDTranslation tran(xpos, ypos, 0.0);
270  DDRotation rotation;
271  int copy = inr * 100 + inc;
272  if (nc < 0) copy += 10000;
273  if (nr < 0) copy += 100000;
274  DDName name =
275  (rpos < rMaxFine)
276  ? DDName(DDSplit(wafer[0]).first, DDSplit(wafer[0]).second)
277  : DDName(DDSplit(wafer[1]).first, DDSplit(wafer[1]).second);
278  cpv.position(name, glog.ddname(), copy, tran, rotation);
279  if (inc > incm) incm = inc;
280  if (inr > inrm) inrm = inr;
281  kount++;
282  if (copies.count(copy) == 0) copies.insert(copy);
283 #ifdef EDM_ML_DEBUG
284  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << name << " number " << copy
285  << " positioned in " << glog.ddname() << " at " << tran
286  << " with " << rotation << std::endl;
287 #endif
288  }
289  }
290  }
291  }
292 #ifdef EDM_ML_DEBUG
293  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: # of columns " << incm
294  << " # of rows " << inrm << " and " << kount
295  << " wafers for " << glog.ddname();
296 #endif
297 }
DDHGCalModuleAlgo::copies
std::unordered_set< int > copies
Definition: DDHGCalModuleAlgo.h:50
popcon2dropbox.copy
def copy(args, dbName)
Definition: popcon2dropbox.py:122
DDName
DDName is used to identify DDD entities uniquely.
Definition: DDName.h:15
DDTranslation
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DDTranslation
Definition: DDTranslation.h:7
DDRotation
Represents a uniquely identifyable rotation matrix.
Definition: DDTransform.h:68
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:187
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
HGCalParameters::k_CornerSize
static uint32_t k_CornerSize
Definition: HGCalParameters.h:21
DDCompactView::position
void position(const DDLogicalPart &self, const DDLogicalPart &parent, const std::string &copyno, const DDTranslation &trans, const DDRotation &rot, const DDDivision *div=0)
Definition: DDCompactView.cc:84
create_public_lumi_plots.ncol
ncol
Definition: create_public_lumi_plots.py:1823
DDSplit
std::pair< std::string, std::string > DDSplit(const std::string &n)
split into (name,namespace), separator = &#39;:&#39;
Definition: DDSplit.cc:3
DDHGCalModuleAlgo::wafer
std::vector< std::string > wafer
Definition: DDHGCalModuleAlgo.h:30
TCMET_cfi.corner
corner
Definition: TCMET_cfi.py:34
DDBase::ddname
const N & ddname() const
Definition: DDBase.h:76
double DDHGCalModuleAlgo::rMax ( double  z)
protected

Definition at line 227 of file DDHGCalModuleAlgo.cc.

References gen::k, and alignCSCRings::r.

227  {
228  double r(0);
229 #ifdef EDM_ML_DEBUG
230  unsigned int ik(0);
231 #endif
232  for (unsigned int k = 0; k < slopeT.size(); ++k) {
233  if (z < zFront[k]) break;
234  r = rMaxFront[k] + (z - zFront[k]) * slopeT[k];
235 #ifdef EDM_ML_DEBUG
236  ik = k;
237 #endif
238  }
239 #ifdef EDM_ML_DEBUG
240  edm::LogVerbatim("HGCalGeom") << "rMax : " << z << ":" << ik << ":" << r << std::endl;
241 #endif
242  return r;
243 }
DDHGCalModuleAlgo::rMaxFront
std::vector< double > rMaxFront
Definition: DDHGCalModuleAlgo.h:47
gen::k
int k[5][pyjets_maxn]
Definition: Cascade2Hadronizer.cc:79
DDHGCalModuleAlgo::slopeT
std::vector< double > slopeT
Definition: DDHGCalModuleAlgo.h:45
DDHGCalModuleAlgo::zFront
std::vector< double > zFront
Definition: DDHGCalModuleAlgo.h:46

Member Data Documentation

std::unordered_set<int> DDHGCalModuleAlgo::copies
private

Definition at line 50 of file DDHGCalModuleAlgo.h.

std::vector<int> DDHGCalModuleAlgo::copyNumber
private

Definition at line 34 of file DDHGCalModuleAlgo.h.

std::string DDHGCalModuleAlgo::idName
private

Definition at line 48 of file DDHGCalModuleAlgo.h.

std::string DDHGCalModuleAlgo::idNameSpace
private

Definition at line 49 of file DDHGCalModuleAlgo.h.

std::vector<int> DDHGCalModuleAlgo::layers
private

Definition at line 35 of file DDHGCalModuleAlgo.h.

std::vector<int> DDHGCalModuleAlgo::layerSense
private

Definition at line 38 of file DDHGCalModuleAlgo.h.

std::vector<double> DDHGCalModuleAlgo::layerThick
private

Definition at line 36 of file DDHGCalModuleAlgo.h.

std::vector<int> DDHGCalModuleAlgo::layerType
private

Definition at line 37 of file DDHGCalModuleAlgo.h.

std::vector<std::string> DDHGCalModuleAlgo::materials
private

Definition at line 31 of file DDHGCalModuleAlgo.h.

std::vector<std::string> DDHGCalModuleAlgo::names
private

Definition at line 32 of file DDHGCalModuleAlgo.h.

double DDHGCalModuleAlgo::rMaxFine
private

Definition at line 40 of file DDHGCalModuleAlgo.h.

std::vector<double> DDHGCalModuleAlgo::rMaxFront
private

Definition at line 47 of file DDHGCalModuleAlgo.h.

int DDHGCalModuleAlgo::sectors
private

Definition at line 43 of file DDHGCalModuleAlgo.h.

std::vector<double> DDHGCalModuleAlgo::slopeB
private

Definition at line 44 of file DDHGCalModuleAlgo.h.

std::vector<double> DDHGCalModuleAlgo::slopeT
private

Definition at line 45 of file DDHGCalModuleAlgo.h.

std::vector<double> DDHGCalModuleAlgo::thick
private

Definition at line 33 of file DDHGCalModuleAlgo.h.

std::vector<std::string> DDHGCalModuleAlgo::wafer
private

Definition at line 30 of file DDHGCalModuleAlgo.h.

double DDHGCalModuleAlgo::waferGap
private

Definition at line 42 of file DDHGCalModuleAlgo.h.

double DDHGCalModuleAlgo::waferW
private

Definition at line 41 of file DDHGCalModuleAlgo.h.

std::vector<double> DDHGCalModuleAlgo::zFront
private

Definition at line 46 of file DDHGCalModuleAlgo.h.

double DDHGCalModuleAlgo::zMinBlock
private

Definition at line 39 of file DDHGCalModuleAlgo.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11