|
|
|
static constexpr double | tol1_ = 0.01 |
|
static constexpr double | tol2_ = 0.00001 |
|
Definition at line 32 of file DDHGCalEEAlgo.cc.
◆ DDHGCalEEAlgo()
DDHGCalEEAlgo::DDHGCalEEAlgo |
( |
| ) |
|
◆ ~DDHGCalEEAlgo()
DDHGCalEEAlgo::~DDHGCalEEAlgo |
( |
| ) |
|
|
override |
◆ constructLayers()
Definition at line 233 of file DDHGCalEEAlgo.cc.
239 for (
unsigned int i = 0;
i <
layers_.size();
i++) {
245 for (
int ly = laymin; ly < laymax; ++ly) {
255 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEAlgo: Layer " << ly <<
":" <<
ii <<
" Front " << zi <<
", " << routF
256 <<
" Back " << zo <<
", " << rinB <<
" superlayer thickness " <<
layerThick_[
i];
262 std::vector<double> pgonZ, pgonRin, pgonRout;
265 pgonZ.emplace_back(-hthick);
266 pgonZ.emplace_back(hthick);
267 pgonRin.emplace_back(rinB);
268 pgonRin.emplace_back(rinB);
269 pgonRout.emplace_back(rmax);
270 pgonRout.emplace_back(rmax);
285 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEAlgo: z " << (
zz - hthick) <<
":" << (
zz + hthick) <<
" with "
286 << pgonZ.size() <<
" palnes";
287 for (
unsigned int isec = 0; isec < pgonZ.size(); ++isec)
289 <<
"[" << isec <<
"] z " << pgonZ[isec] <<
" R " << pgonRin[isec] <<
":" << pgonRout[isec];
291 for (
unsigned int isec = 0; isec < pgonZ.size(); ++isec) {
303 <<
" sections and filled with " << matName <<
":" << &matter;
304 for (
unsigned int k = 0;
k < pgonZ.size(); ++
k)
305 edm::LogVerbatim(
"HGCalGeom") <<
"[" <<
k <<
"] z " << pgonZ[
k] <<
" R " << pgonRin[
k] <<
":" << pgonRout[
k];
315 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEFileAlgo: " << solid.
name() <<
" Tubs made of " << matName <<
":"
316 << &matter <<
" of dimensions " << rinB <<
":" << rins <<
", " << routF <<
":"
317 << routs <<
", " << hthick <<
", 0.0, 360.0 and position " << glog.
name()
328 << module.
name() <<
" at " <<
r1 <<
" with " <<
rot;
337 <<
": thickness of all its components **** ERROR ****";
340 << 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()
◆ positionSensitive()
void DDHGCalEEAlgo::positionSensitive |
( |
const DDLogicalPart & |
glog, |
|
|
double |
rin, |
|
|
double |
rout, |
|
|
double |
zpos, |
|
|
int |
layertype, |
|
|
int |
layercenter, |
|
|
DDCompactView & |
cpv |
|
) |
| |
|
protected |
Definition at line 346 of file DDHGCalEEAlgo.cc.
353 static const double sqrt3 =
std::sqrt(3.0);
355 double R = 2.0 *
r / sqrt3;
356 double dy = 0.75 *
R;
357 int N = (
int)(0.5 * rout /
r) + 2;
360 int ium(0), ivm(0), iumAll(0), ivmAll(0), kount(0), ntot(0),
nin(0);
361 std::vector<int> ntype(6, 0);
362 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEAlgo: " << glog.
ddname() <<
" rin:rout " << rin <<
":" << rout <<
" zpos "
363 << zpos <<
" N " <<
N <<
" for maximum u, v; r " <<
r <<
" R " <<
R <<
" dy " <<
dy
364 <<
" Shift " << xyoff.first <<
":" << xyoff.second <<
" WaferSize "
367 for (
int u = -
N; u <=
N; ++u) {
368 for (
int v = -
N;
v <=
N; ++
v) {
371 double xpos = xyoff.first + nc *
r;
372 double ypos = xyoff.second +
nr *
dy;
379 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEAlgo: " << glog.
ddname() <<
" R " << rin <<
":" << rout <<
"\n Z "
380 << zpos <<
" LayerType " << layertype <<
" u " << u <<
" v " <<
v <<
" with "
381 <<
corner.first <<
" corners";
420 edm::LogVerbatim(
"HGCalGeom") <<
"DDHGCalEEAlgo: Maximum # of u " << ium <<
":" << iumAll <<
" # of v " << ivm <<
":"
421 << ivmAll <<
" and " <<
nin <<
":" << kount <<
":" << ntot <<
" wafers (" << ntype[0]
422 <<
":" << ntype[1] <<
":" << ntype[2] <<
":" << ntype[3] <<
":" << ntype[4] <<
":"
423 << 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, 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, and HGCalGeomTools::waferCorner().
◆ absorbMode_
int DDHGCalEEAlgo::absorbMode_ |
|
private |
◆ alpha_
double DDHGCalEEAlgo::alpha_ |
|
private |
◆ choiceType_
int DDHGCalEEAlgo::choiceType_ |
|
private |
◆ copies_
std::unordered_set<int> DDHGCalEEAlgo::copies_ |
|
private |
◆ copyNumber_
std::vector<int> DDHGCalEEAlgo::copyNumber_ |
|
private |
◆ cosAlpha_
double DDHGCalEEAlgo::cosAlpha_ |
|
private |
◆ firstLayer_
int DDHGCalEEAlgo::firstLayer_ |
|
private |
◆ fracAreaMin_
double DDHGCalEEAlgo::fracAreaMin_ |
|
private |
◆ geomTools_
◆ layerCenter_
std::vector<int> DDHGCalEEAlgo::layerCenter_ |
|
private |
◆ layers_
std::vector<int> DDHGCalEEAlgo::layers_ |
|
private |
◆ layerSense_
std::vector<int> DDHGCalEEAlgo::layerSense_ |
|
private |
◆ layerThick_
std::vector<double> DDHGCalEEAlgo::layerThick_ |
|
private |
◆ layerType_
std::vector<int> DDHGCalEEAlgo::layerType_ |
|
private |
◆ materials_
std::vector<std::string> DDHGCalEEAlgo::materials_ |
|
private |
◆ names_
std::vector<std::string> DDHGCalEEAlgo::names_ |
|
private |
◆ nameSpace_
std::string DDHGCalEEAlgo::nameSpace_ |
|
private |
◆ nCutRadPar_
int DDHGCalEEAlgo::nCutRadPar_ |
|
private |
◆ rad100to200_
std::vector<double> DDHGCalEEAlgo::rad100to200_ |
|
private |
◆ rad200to300_
std::vector<double> DDHGCalEEAlgo::rad200to300_ |
|
private |
◆ rMaxFront_
std::vector<double> DDHGCalEEAlgo::rMaxFront_ |
|
private |
◆ rMinFront_
std::vector<double> DDHGCalEEAlgo::rMinFront_ |
|
private |
◆ sectors_
int DDHGCalEEAlgo::sectors_ |
|
private |
◆ sensitiveMode_
int DDHGCalEEAlgo::sensitiveMode_ |
|
private |
◆ slopeB_
std::vector<double> DDHGCalEEAlgo::slopeB_ |
|
private |
◆ slopeT_
std::vector<double> DDHGCalEEAlgo::slopeT_ |
|
private |
◆ thick_
std::vector<double> DDHGCalEEAlgo::thick_ |
|
private |
◆ tol1_
constexpr double DDHGCalEEAlgo::tol1_ = 0.01 |
|
staticconstexprprivate |
◆ tol2_
constexpr double DDHGCalEEAlgo::tol2_ = 0.00001 |
|
staticconstexprprivate |
◆ wafers_
std::vector<std::string> DDHGCalEEAlgo::wafers_ |
|
private |
◆ waferSepar_
double DDHGCalEEAlgo::waferSepar_ |
|
private |
◆ waferSize_
double DDHGCalEEAlgo::waferSize_ |
|
private |
◆ waferType_
◆ zFrontB_
std::vector<double> DDHGCalEEAlgo::zFrontB_ |
|
private |
◆ zFrontT_
std::vector<double> DDHGCalEEAlgo::zFrontT_ |
|
private |
◆ zMinBlock_
double DDHGCalEEAlgo::zMinBlock_ |
|
private |
◆ zMinRadPar_
double DDHGCalEEAlgo::zMinRadPar_ |
|
private |
std::vector< double > layerThick_
std::vector< double > slopeT_
std::vector< int > layerCenter_
std::vector< double > rMinFront_
DDName is used to identify DDD entities uniquely.
HGCalGeomTools geomTools_
constexpr NumType convertRadToDeg(NumType radians)
static constexpr uint32_t k_CornerSize
U second(std::pair< T, U > const &p)
std::vector< int > layerType_
void constructLayers(const DDLogicalPart &, DDCompactView &cpv)
std::vector< int > dbl_to_int(const std::vector< double > &vecdbl)
Converts a std::vector of doubles to a std::vector of int.
Log< level::Warning, false > LogWarning
void positionSensitive(const DDLogicalPart &glog, double rin, double rout, double zpos, int layertype, int layercenter, DDCompactView &cpv)
std::vector< std::string > materials_
std::unique_ptr< HGCalWaferType > waferType_
std::vector< int > layers_
DDMaterial is used to define and access material information.
std::vector< double > rad200to300_
std::vector< double > thick_
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)
std::vector< int > copyNumber_
std::vector< int > layerSense_
static DDSolid tubs(const DDName &name, double zhalf, double rIn, double rOut, double startPhi, double deltaPhi)
static constexpr double tol1_
static int32_t packTypeUV(int type, int u, int v)
A DDLogicalPart aggregates information concerning material, solid and sensitveness ....
std::vector< double > rad100to200_
std::vector< double > rMaxFront_
Log< level::Error, false > LogError
std::vector< std::pair< float, float > >::iterator itr
std::vector< double > slopeB_
static constexpr double tol2_
Log< level::Info, true > LogVerbatim
static std::string & ns()
std::vector< double > zFrontT_
std::unordered_set< int > copies_
A DDSolid represents the shape of a part.
Represents a uniquely identifyable rotation matrix.
std::vector< std::string > wafers_
Abs< T >::type abs(const T &t)
std::vector< double > zFrontB_
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)
std::vector< std::string > names_