|
|
|
static constexpr double | tol1_ = 0.01 |
|
static constexpr double | tol2_ = 0.00001 |
|
Definition at line 34 of file DDHGCalEEFileAlgo.cc.
◆ DDHGCalEEFileAlgo()
DDHGCalEEFileAlgo::DDHGCalEEFileAlgo |
( |
| ) |
|
◆ constructLayers()
Definition at line 234 of file DDHGCalEEFileAlgo.cc.
240 for (
unsigned int i = 0;
i <
layers_.size();
i++) {
246 for (
int ly = laymin; ly < laymax; ++ly) {
256 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEFileAlgo: Layer " << ly <<
":" <<
ii <<
" Front " << zi <<
", "
257 << routF <<
" Back " << zo <<
", " << rinB <<
" superlayer thickness "
264 std::vector<double> pgonZ, pgonRin, pgonRout;
267 pgonZ.emplace_back(-hthick);
268 pgonZ.emplace_back(hthick);
269 pgonRin.emplace_back(rinB);
270 pgonRin.emplace_back(rinB);
271 pgonRout.emplace_back(rmax);
272 pgonRout.emplace_back(rmax);
287 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEFileAlgo: z " << (
zz - hthick) <<
":" << (
zz + hthick) <<
" with "
288 << pgonZ.size() <<
" palnes";
289 for (
unsigned int isec = 0; isec < pgonZ.size(); ++isec)
291 <<
"[" << isec <<
"] z " << pgonZ[isec] <<
" R " << pgonRin[isec] <<
":" << pgonRout[isec];
293 for (
unsigned int isec = 0; isec < pgonZ.size(); ++isec) {
305 <<
" sections and filled with " << matName <<
":" << &matter;
306 for (
unsigned int k = 0;
k < pgonZ.size(); ++
k)
307 edm::LogVerbatim(
"HGCalGeom") <<
"[" <<
k <<
"] z " << pgonZ[
k] <<
" R " << pgonRin[
k] <<
":" << pgonRout[
k];
317 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEFileAlgo: " << solid.
name() <<
" Tubs made of " << matName <<
":"
318 << &matter <<
" of dimensions " << rinB <<
":" << rins <<
", " << routF <<
":"
319 << routs <<
", " << hthick <<
", 0.0, 360.0 and position " << glog.
name()
330 << module.
name() <<
" at " <<
r1 <<
" with " <<
rot;
340 <<
": thickness of all its components **** ERROR ****";
343 << thickTot <<
" of the components";
References funct::abs(), angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, DDBase< N, C >::ddname(), DDSplit(), dqmdumpme::first, mps_fire::i, cuy::ii, dqmdumpme::k, Skims_PA_cff::name, DDBase< N, C >::name(), DDSolidFactory::polyhedra(), DDCompactView::position(), diffTwoXMLs::r1, HGCalGeomTools::radius(), makeMuonMisalignmentScenario::rot, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, DDSolidFactory::tubs(), and geometryCSVtoXML::zz.
◆ execute()
◆ initialize()
Definition at line 98 of file DDHGCalEEFileAlgo.cc.
103 wafers_ = vsArgs[
"WaferNames"];
106 for (
unsigned int i = 0;
i <
wafers_.size(); ++
i)
110 names_ = vsArgs[
"VolumeNames"];
111 thick_ = vArgs[
"Thickness"];
115 for (
unsigned int i = 0;
i <
names_.size(); ++
i)
123 for (
unsigned int i = 0;
i <
layers_.size(); ++
i)
147 edm::LogVerbatim(
"HGCalGeom") <<
"First copy number for layer type " <<
i <<
":" <<
ii <<
" with "
193 slopeB_ = vArgs[
"SlopeBottom"];
200 for (
unsigned int i = 0;
i <
slopeB_.size(); ++
i)
203 for (
unsigned int i = 0;
i <
slopeT_.size(); ++
i)
References angle_units::operators::convertRadToDeg(), funct::cos(), dbl_to_int(), mps_fire::i, cuy::ii, createfilelist::int, dqmdumpme::k, DDCurrentNamespace::ns(), HGCalWaferIndex::waferLayer(), HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().
◆ positionSensitive()
void DDHGCalEEFileAlgo::positionSensitive |
( |
const DDLogicalPart & |
glog, |
|
|
double |
rin, |
|
|
double |
rout, |
|
|
double |
zpos, |
|
|
int |
layertype, |
|
|
int |
layer, |
|
|
DDCompactView & |
cpv |
|
) |
| |
|
protected |
Definition at line 349 of file DDHGCalEEFileAlgo.cc.
351 static const double sqrt3 =
std::sqrt(3.0);
354 double R = 2.0 *
r / sqrt3;
355 double dy = 0.75 *
R;
356 int N = (
int)(0.5 * rout /
r) + 2;
359 int ium(0), ivm(0), iumAll(0), ivmAll(0), kount(0), ntot(0),
nin(0);
360 std::vector<int> ntype(6, 0);
361 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEFileAlgo: " << glog.
ddname() <<
" rin:rout " << rin <<
":" << rout
362 <<
" zpos " << zpos <<
" N " <<
N <<
" for maximum u, v; r " <<
r <<
" R " <<
R
363 <<
" dy " <<
dy <<
" Shift " << xyoff.first <<
":" << xyoff.second <<
" WaferSize "
366 for (
int u = -
N; u <=
N; ++u) {
367 for (
int v = -
N;
v <=
N; ++
v) {
374 double xpos = xyoff.first + nc *
r;
375 double ypos = xyoff.second +
nr *
dy;
381 <<
"\n Z " << zpos <<
" LayerType " << layertype <<
" u " << u <<
" v " <<
v
382 <<
" with " <<
corner.first <<
" corners";
391 <<
type <<
" layer:u:v:indx " << (layer +
firstLayer_) <<
":" << u <<
":" <<
v
418 <<
":" <<
type <<
" positioned in " << glog.
ddname() <<
" at " << tran
426 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEFileAlgo: Maximum # of u " << ium <<
":" << iumAll <<
" # of v " << ivm
427 <<
":" << ivmAll <<
" and " <<
nin <<
":" << kount <<
":" << ntot <<
" wafers ("
428 << ntype[0] <<
":" << ntype[1] <<
":" << ntype[2] <<
":" << ntype[3] <<
":" << ntype[4]
429 <<
":" << ntype[5] <<
") for " << glog.
ddname() <<
" R " << rin <<
":" << rout;
References funct::abs(), filterCSVwithJSON::copy, distTCMET_cfi::corner, DDBase< N, C >::ddname(), DDSplit(), PVValHelper::dy, dqmdumpme::first, HGCalWaferType::getType(), createfilelist::int, HGCalParameters::k_CornerSize, N, Skims_PA_cff::name, nin, EgHLTOffHistBins_cfi::nr, HGCalTypes::packTypeUV(), DDCompactView::position(), dttmaxenums::R, alignCSCRings::r, idealTransformation::rotation, edm::second(), mathSSE::sqrt(), findQualityFiles::v, HGCalGeomTools::waferCorner(), and HGCalWaferIndex::waferIndex().
◆ absorbMode_
int DDHGCalEEFileAlgo::absorbMode_ |
|
private |
◆ alpha_
double DDHGCalEEFileAlgo::alpha_ |
|
private |
◆ choiceType_
int DDHGCalEEFileAlgo::choiceType_ |
|
private |
◆ copies_
std::unordered_set<int> DDHGCalEEFileAlgo::copies_ |
|
private |
◆ copyNumber_
std::vector<int> DDHGCalEEFileAlgo::copyNumber_ |
|
private |
◆ cosAlpha_
double DDHGCalEEFileAlgo::cosAlpha_ |
|
private |
◆ firstLayer_
int DDHGCalEEFileAlgo::firstLayer_ |
|
private |
◆ fracAreaMin_
double DDHGCalEEFileAlgo::fracAreaMin_ |
|
private |
◆ geomTools_
◆ layerCenter_
std::vector<int> DDHGCalEEFileAlgo::layerCenter_ |
|
private |
◆ layers_
std::vector<int> DDHGCalEEFileAlgo::layers_ |
|
private |
◆ layerSense_
std::vector<int> DDHGCalEEFileAlgo::layerSense_ |
|
private |
◆ layerThick_
std::vector<double> DDHGCalEEFileAlgo::layerThick_ |
|
private |
◆ layerType_
std::vector<int> DDHGCalEEFileAlgo::layerType_ |
|
private |
◆ materials_
std::vector<std::string> DDHGCalEEFileAlgo::materials_ |
|
private |
◆ names_
std::vector<std::string> DDHGCalEEFileAlgo::names_ |
|
private |
◆ nameSpace_
std::string DDHGCalEEFileAlgo::nameSpace_ |
|
private |
◆ nCutRadPar_
int DDHGCalEEFileAlgo::nCutRadPar_ |
|
private |
◆ rad100to200_
std::vector<double> DDHGCalEEFileAlgo::rad100to200_ |
|
private |
◆ rad200to300_
std::vector<double> DDHGCalEEFileAlgo::rad200to300_ |
|
private |
◆ rMaxFront_
std::vector<double> DDHGCalEEFileAlgo::rMaxFront_ |
|
private |
◆ rMinFront_
std::vector<double> DDHGCalEEFileAlgo::rMinFront_ |
|
private |
◆ sectors_
int DDHGCalEEFileAlgo::sectors_ |
|
private |
◆ sensitiveMode_
int DDHGCalEEFileAlgo::sensitiveMode_ |
|
private |
◆ slopeB_
std::vector<double> DDHGCalEEFileAlgo::slopeB_ |
|
private |
◆ slopeT_
std::vector<double> DDHGCalEEFileAlgo::slopeT_ |
|
private |
◆ thick_
std::vector<double> DDHGCalEEFileAlgo::thick_ |
|
private |
◆ tol1_
constexpr double DDHGCalEEFileAlgo::tol1_ = 0.01 |
|
staticconstexprprivate |
◆ tol2_
constexpr double DDHGCalEEFileAlgo::tol2_ = 0.00001 |
|
staticconstexprprivate |
◆ waferIndex_
std::vector<int> DDHGCalEEFileAlgo::waferIndex_ |
|
private |
◆ wafers_
std::vector<std::string> DDHGCalEEFileAlgo::wafers_ |
|
private |
◆ waferSepar_
double DDHGCalEEFileAlgo::waferSepar_ |
|
private |
◆ waferSize_
double DDHGCalEEFileAlgo::waferSize_ |
|
private |
◆ waferTypes_
std::vector<int> DDHGCalEEFileAlgo::waferTypes_ |
|
private |
◆ zFrontB_
std::vector<double> DDHGCalEEFileAlgo::zFrontB_ |
|
private |
◆ zFrontT_
std::vector<double> DDHGCalEEFileAlgo::zFrontT_ |
|
private |
◆ zMinBlock_
double DDHGCalEEFileAlgo::zMinBlock_ |
|
private |
◆ zMinRadPar_
double DDHGCalEEFileAlgo::zMinRadPar_ |
|
private |
std::vector< double > zFrontT_
static int32_t waferV(const int32_t index)
std::vector< int > layerType_
DDName is used to identify DDD entities uniquely.
std::vector< std::string > wafers_
constexpr NumType convertRadToDeg(NumType radians)
void constructLayers(const DDLogicalPart &, DDCompactView &cpv)
static constexpr uint32_t k_CornerSize
std::vector< double > zFrontB_
U second(std::pair< T, U > const &p)
std::vector< int > layerSense_
std::vector< int > dbl_to_int(const std::vector< double > &vecdbl)
Converts a std::vector of doubles to a std::vector of int.
std::unordered_set< int > copies_
Log< level::Warning, false > LogWarning
std::vector< int > waferIndex_
std::vector< double > rad100to200_
std::vector< std::string > materials_
DDMaterial is used to define and access material information.
std::vector< int > copyNumber_
Cos< T >::type cos(const T &t)
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DDTranslation
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)
HGCalGeomTools geomTools_
static DDSolid tubs(const DDName &name, double zhalf, double rIn, double rOut, double startPhi, double deltaPhi)
std::vector< int > layerCenter_
int getType(double xpos, double ypos, double zpos)
std::vector< double > slopeT_
static int32_t waferU(const int32_t index)
static int32_t packTypeUV(int type, int u, int v)
std::vector< int > waferTypes_
std::vector< double > rad200to300_
static constexpr double tol2_
std::vector< double > rMinFront_
void positionSensitive(const DDLogicalPart &glog, double rin, double rout, double zpos, int layertype, int layer, DDCompactView &cpv)
A DDLogicalPart aggregates information concerning material, solid and sensitveness ....
static constexpr double tol1_
std::vector< std::string > names_
std::vector< double > thick_
Log< level::Error, false > LogError
std::vector< std::pair< float, float > >::iterator itr
std::vector< double > rMaxFront_
std::vector< int > layers_
std::vector< double > slopeB_
Log< level::Info, true > LogVerbatim
static std::string & ns()
static int32_t waferIndex(int32_t layer, int32_t waferU, int32_t waferV, bool old=false)
A DDSolid represents the shape of a part.
Represents a uniquely identifyable rotation matrix.
std::vector< double > layerThick_
Abs< T >::type abs(const T &t)
static int32_t waferLayer(const int32_t index)
std::pair< std::string, std::string > DDSplit(const std::string &n)
split into (name,namespace), separator = ':'
void position(const DDLogicalPart &self, const DDLogicalPart &parent, const std::string ©no, const DDTranslation &trans, const DDRotation &rot, const DDDivision *div=nullptr)