Inheritance diagram for DDHCalEndcapAlgo:

Public Member Functions
	DDHCalEndcapAlgo ()

int	equipModule (unsigned int i) const

void	execute (DDCompactView &cpv) override

const std::string &	getAbsMat () const

double	getAngBot () const

double	getAngGap () const

double	getAngTop () const

double	getDrEnd () const

double	getDzShift () const

double	getDzStep () const

int	getEndcaps () const

const std::string &	getGenMat () const

double	getHeboxDepth () const

int	getLayer (unsigned int i, unsigned int j) const

int	getLayerN (unsigned int i) const

const std::string &	getLayerName (unsigned int i) const

int	getLayers () const

double	getLayerT (unsigned int i) const

int	getLayerType (unsigned int i) const

const std::string &	getModMat (unsigned int i) const

const std::string &	getModName (unsigned int i) const

int	getModType (unsigned int i) const

int	getModules () const

int	getNsectors () const

int	getNsectortot () const

int	getPhi () const

const std::string &	getPhiName (unsigned int i) const

const std::string &	getPlastMat () const

double	getRinBlock1 (unsigned i) const

double	getRinBlock2 (unsigned i) const

double	getRinDip () const

double	getRinKink () const

const std::string &	getRotation () const

const std::string &	getRotMat () const

double	getRout () const

double	getRoutBlock1 (unsigned i) const

double	getRoutBlock2 (unsigned i) const

double	getRoutDip () const

const std::string &	getScintMat () const

double	getScintT (unsigned int i) const

int	getSectionModule (unsigned i) const

double	getSlope () const

double	getThick (unsigned int i) const

double	getTolAbs () const

double	getTolPos () const

double	getTrim (unsigned int i, unsigned int j) const

double	getZ0Beam () const

double	getZ1Beam () const

double	getZEnd () const

double	getZFront () const

double	getZiBody () const

double	getZiDip () const

double	getZiKink () const

double	getZiL0Body () const

double	getZiL0Nose () const

double	getZiNose () const

double	getZmaxBlock (unsigned i) const

double	getZminBlock (unsigned i) const

double	getZShift () const

double	getZShiftHac2 () const

void	initialize (const DDNumericArguments &nArgs, const DDVectorArguments &vArgs, const DDMapArguments &mArgs, const DDStringArguments &sArgs, const DDStringVectorArguments &vsArgs) override

	~DDHCalEndcapAlgo () override

Protected Member Functions
void	constructGeneralVolume (DDCompactView &cpv)

void	constructInsideModule (const DDLogicalPart &module, int mod, DDCompactView &cpv)

void	constructInsideModule0 (const DDLogicalPart &module, int mod, DDCompactView &cpv)

void	constructInsideSector (const DDLogicalPart &sector, DDCompactView &cpv)

void	constructScintLayer (const DDLogicalPart &glog, double pDz, double yh, double bl, double tl, double alp, const std::string &name, int id, DDCompactView &cpv)

void	parameterLayer (int iphi, double rinF, double routF, double rinB, double routB, double zi, double zo, double &yh1, double &bl1, double &tl1, double &yh2, double &bl2, double &tl2, double &alp, double &theta, double &phi, double &xpos, double &ypos, double &zcpv)

void	parameterLayer0 (int mod, int layer, int iphi, double &yh, double &bl, double &tl, double &alp, double &xpos, double &ypos, double &zcpv)

Private Attributes
std::string	absMat

double	angBot

double	angGap

double	angTop

double	drEnd

double	dzShift

double	dzStep

std::vector< int >	eModule

std::string	genMaterial

double	heboxDepth

std::string	idName

std::string	idNameSpace

int	idOffset

std::vector< int >	layerN

std::vector< int >	layerN0

std::vector< int >	layerN1

std::vector< int >	layerN2

std::vector< int >	layerN3

std::vector< int >	layerN4

std::vector< int >	layerN5

std::vector< std::string >	layerName

int	layers

std::vector< double >	layerT

std::vector< int >	layerType

std::vector< std::string >	modMat

std::vector< std::string >	modName

std::vector< int >	modType

int	modules

int	nEndcap

int	nsectors

int	nsectortot

std::vector< std::string >	phiName

int	phiSections

std::string	plastMat

double	riDip

double	riKink

std::vector< double >	rinBlock1

std::vector< double >	rinBlock2

double	roDip

std::string	rotation

std::string	rotHalf

std::string	rotmat

std::string	rotns

double	rout

std::vector< double >	routBlock1

std::vector< double >	routBlock2

std::string	scintMat

std::vector< double >	scintT

std::vector< int >	sectionModule

double	slope

std::vector< double >	thick

double	tolAbs

double	tolPos

std::vector< double >	trimLeft

std::vector< double >	trimRight

double	z0Beam

double	z1Beam

double	zEnd

double	zFront

double	ziBody

double	ziDip

double	ziKink

double	ziL0Body

double	ziL0Nose

double	ziNose

std::vector< double >	zmaxBlock

std::vector< double >	zminBlock

double	zShift

double	zShiftHac2

Detailed Description

Definition at line 29 of file DDHCalEndcapAlgo.cc.

Constructor & Destructor Documentation

◆ DDHCalEndcapAlgo()

DDHCalEndcapAlgo::DDHCalEndcapAlgo ( )

Definition at line 229 of file DDHCalEndcapAlgo.cc.

     : modMat(0),
       modType(0),
       sectionModule(0),
       layerN(0),
       layerN0(0),
       layerN1(0),
       layerN2(0),
       layerN3(0),
       layerN4(0),
       layerN5(0),
       thick(0),
       trimLeft(0),
       trimRight(0),
       zminBlock(0),
       zmaxBlock(0),
       rinBlock1(0),
       routBlock1(0),
       rinBlock2(0),
       routBlock2(0),
       layerType(0),
       layerT(0),
       scintT(0) {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Creating an instance";
 #endif
 }

◆ ~DDHCalEndcapAlgo()

DDHCalEndcapAlgo::~DDHCalEndcapAlgo ( )

override

Definition at line 257 of file DDHCalEndcapAlgo.cc.

257 {}

Member Function Documentation

◆ constructGeneralVolume()

void DDHCalEndcapAlgo::constructGeneralVolume ( DDCompactView & cpv )

protected

!!!!!!!!!!!!!!!!Should be zero. And removed as soon as

Definition at line 555 of file DDHCalEndcapAlgo.cc.

References simBeamSpotPI::alpha, constructInsideSector(), angle_units::operators::convertRadToDeg(), DDrot(), DDSplit(), equipModule(), dqmdumpme::first, formatAsDegreesInInteger(), getAbsMat(), getAngBot(), getAngTop(), getDrEnd(), getDzShift(), getEndcaps(), getGenMat(), getNsectors(), getNsectortot(), getRinDip(), getRinKink(), getRotation(), getRout(), getRoutDip(), getSlope(), getZ1Beam(), getZEnd(), getZFront(), getZiBody(), getZiDip(), getZiKink(), getZiL0Body(), getZShift(), getZShiftHac2(), mps_fire::i, idName, idNameSpace, cuy::ii, Skims_PA_cff::name, DDBase< N, C >::name(), class-composition::parent, phi, DDSolidFactory::polyhedra(), DDCompactView::position(), makeMuonMisalignmentScenario::rot, rotation, rotHalf, rotns, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, and funct::tan().

Referenced by execute().

                                                                 {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: General volume...";
 #endif
 
   bool proto = true;
   for (int i = 0; i < 3; i++)
     if (equipModule(i) > 0)
       proto = false;
 
   DDRotation rot;
   if (DDSplit(getRotation()).first == "NULL")
     rot = DDRotation();
   else
     rot = DDRotation(DDName(DDSplit(getRotation()).first, DDSplit(getRotation()).second));
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << " First " << DDSplit(getRotation()).first << " Second "
                                << DDSplit(getRotation()).second << " Rotation " << rot;
 #endif
 
   DDTranslation r0(0, 0, getZShift());
   double alpha = (1._pi) / getNsectors();
   double dphi = getNsectortot() * (2._pi) / getNsectors();
 
   //vertical walls are allowed in SolidPolyhedra
   double delz = 0;
 
   std::vector<double> pgonZ, pgonRmin, pgonRmax;
   if (proto) {
     double zf = getZiBody() + getZShiftHac2();
     pgonZ.emplace_back(zf - getDzShift());
     pgonRmin.emplace_back(zf * tan(getAngBot()));
     pgonRmax.emplace_back((zf - getZ1Beam()) * getSlope());
   } else {
     pgonZ.emplace_back(getZFront() - getDzShift());
     pgonRmin.emplace_back(getZFront() * tan(getAngTop()));
     pgonRmax.emplace_back((getZFront() - getZ1Beam()) * getSlope());
     pgonZ.emplace_back(getZiL0Body() - getDzShift());
     pgonRmin.emplace_back(getZiL0Body() * tan(getAngTop()));
     pgonRmax.emplace_back((getZiL0Body() - getZ1Beam()) * getSlope());
     pgonZ.emplace_back(getZiL0Body() - getDzShift());
     pgonRmin.emplace_back(getZiL0Body() * tan(getAngBot()));
     pgonRmax.emplace_back((getZiL0Body() - getZ1Beam()) * getSlope());
   }
   pgonZ.emplace_back(getZiKink() - getDzShift());
   pgonRmin.emplace_back(getRinKink());
   pgonRmax.emplace_back(getRout());
   pgonZ.emplace_back(getZiDip() - getDzShift());
   pgonRmin.emplace_back(getRinDip());
   pgonRmax.emplace_back(getRout());
   pgonZ.emplace_back(getZiDip() - getDzShift() + delz);
   pgonRmin.emplace_back(getRinDip());
   pgonRmax.emplace_back(getRoutDip());
   pgonZ.emplace_back(getZEnd() - getDzShift());
   pgonRmin.emplace_back(getZEnd() * tan(getAngBot()));
   pgonRmax.emplace_back(getRoutDip());
   pgonZ.emplace_back(getZEnd());
   pgonRmin.emplace_back(getZEnd() * tan(getAngBot()));
   pgonRmax.emplace_back(getRoutDip());
 
   std::string name("Null");
   DDSolid solid;
   solid =
       DDSolidFactory::polyhedra(DDName(idName, idNameSpace), getNsectortot(), -alpha, dphi, pgonZ, pgonRmin, pgonRmax);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << DDName(idName, idNameSpace) << " Polyhedra made of "
                                << getGenMat() << " with " << getNsectortot() << " sectors from "
                                << convertRadToDeg(-alpha) << " to " << convertRadToDeg(-alpha + dphi) << " and with "
                                << pgonZ.size() << " sections";
   for (unsigned int i = 0; i < pgonZ.size(); i++)
     edm::LogVerbatim("HCalGeom") << "\t\tZ = " << pgonZ[i] << "\tRmin = " << pgonRmin[i] << "\tRmax = " << pgonRmax[i];
 #endif
 
   DDName matname(DDSplit(getGenMat()).first, DDSplit(getGenMat()).second);
   DDMaterial matter(matname);
   DDLogicalPart genlogic(DDName(idName, idNameSpace), matter, solid);
 
   DDName parentName = parent().name();
   cpv.position(DDName(idName, idNameSpace), parentName, 1, r0, rot);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << DDName(idName, idNameSpace) << " number 1 positioned in "
                                << parentName << " at " << r0 << " with " << rot;
 #endif
 
   if (getEndcaps() != 1) {
     rot = DDRotation(DDName(rotHalf, rotns));
     cpv.position(DDName(idName, idNameSpace), parentName, 2, r0, rot);
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << DDName(idName, idNameSpace) << " number 2 "
                                  << "positioned in " << parentName << " at " << r0 << " with " << rot;
 #endif
   }
 
   //Forward half
   name = idName + "Front";
   std::vector<double> pgonZMod, pgonRminMod, pgonRmaxMod;
   for (unsigned int i = 0; i < (pgonZ.size() - 1); i++) {
     pgonZMod.emplace_back(pgonZ[i] + getDzShift());
     pgonRminMod.emplace_back(pgonRmin[i]);
     pgonRmaxMod.emplace_back(pgonRmax[i]);
   }
   solid = DDSolidFactory::polyhedra(
       DDName(name, idNameSpace), getNsectortot(), -alpha, dphi, pgonZMod, pgonRminMod, pgonRmaxMod);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << DDName(name, idNameSpace) << " Polyhedra made of "
                                << getGenMat() << " with " << getNsectortot() << " sectors from "
                                << convertRadToDeg(-alpha) << " to " << convertRadToDeg(-alpha + dphi) << " and with "
                                << pgonZMod.size() << " sections ";
   for (unsigned int i = 0; i < pgonZMod.size(); i++)
     edm::LogVerbatim("HCalGeom") << "\t\tZ = " << pgonZMod[i] << "\tRmin = " << pgonRminMod[i]
                                  << "\tRmax = " << pgonRmaxMod[i];
 #endif
 
   DDLogicalPart genlogich(DDName(name, idNameSpace), matter, solid);
 
   cpv.position(genlogich, genlogic, 1, DDTranslation(0.0, 0.0, -getDzShift()), DDRotation());
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << genlogich.name() << " number 1 positioned in "
                                << genlogic.name() << " at (0,0," << -getDzShift() << ") with no rotation";
 #endif
 
   //Construct sector (from -alpha to +alpha)
   name = idName + "Module";
   solid =
       DDSolidFactory::polyhedra(DDName(name, idNameSpace), 1, -alpha, 2 * alpha, pgonZMod, pgonRminMod, pgonRmaxMod);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << DDName(name, idNameSpace) << " Polyhedra made of "
                                << getGenMat() << " with 1 sector from " << convertRadToDeg(-alpha) << " to "
                                << convertRadToDeg(alpha) << " and with " << pgonZMod.size() << " sections";
   for (unsigned int i = 0; i < pgonZMod.size(); i++)
     edm::LogVerbatim("HCalGeom") << "\t\tZ = " << pgonZMod[i] << "\tRmin = " << pgonRminMod[i]
                                  << "\tRmax = " << pgonRmaxMod[i];
 #endif
 
   DDLogicalPart seclogic(DDName(name, idNameSpace), matter, solid);
 
   for (int ii = 0; ii < getNsectortot(); ii++) {
     double phi = ii * 2 * alpha;
     DDRotation rotation;
     std::string rotstr("NULL");
     if (phi != 0) {
       rotstr = "R" + formatAsDegreesInInteger(phi);
       rotation = DDRotation(DDName(rotstr, rotns));
       if (!rotation) {
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Creating a new rotation " << rotstr << "\t 90,"
                                      << convertRadToDeg(phi) << ", 90," << convertRadToDeg(phi + 90._deg) << ", 0, 0";
 #endif
         rotation = DDrot(DDName(rotstr, rotns), 90._deg, phi, 90._deg, (90._deg + phi), 0, 0);
       }  //if !rotation
     }    //if phi!=0
 
     cpv.position(seclogic, genlogich, ii + 1, DDTranslation(0.0, 0.0, 0.0), rotation);
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << seclogic.name() << " number " << ii + 1 << " positioned in "
                                  << genlogich.name() << " at (0,0,0) with " << rotation;
 #endif
   }
 
   //Construct the things inside the sector
   constructInsideSector(seclogic, cpv);
 
   //Backward half
   name = idName + "Back";
   std::vector<double> pgonZBack, pgonRminBack, pgonRmaxBack;
   pgonZBack.emplace_back(getZEnd() - getDzShift());
   pgonRminBack.emplace_back(pgonZBack[0] * tan(getAngBot()) + getDrEnd());
   pgonRmaxBack.emplace_back(getRoutDip());
   pgonZBack.emplace_back(getZEnd());
   pgonRminBack.emplace_back(pgonZBack[1] * tan(getAngBot()) + getDrEnd());
   pgonRmaxBack.emplace_back(getRoutDip());
   solid = DDSolidFactory::polyhedra(
       DDName(name, idNameSpace), getNsectortot(), -alpha, dphi, pgonZBack, pgonRminBack, pgonRmaxBack);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << DDName(name, idNameSpace) << " Polyhedra made of "
                                << getAbsMat() << " with " << getNsectortot() << " sectors from "
                                << convertRadToDeg(-alpha) << " to " << convertRadToDeg(-alpha + dphi) << " and with "
                                << pgonZBack.size() << " sections";
   for (unsigned int i = 0; i < pgonZBack.size(); i++)
     edm::LogVerbatim("HCalGeom") << "\t\tZ = " << pgonZBack[i] << "\tRmin = " << pgonRminBack[i]
                                  << "\tRmax = " << pgonRmaxBack[i];
 #endif
 
   DDName absMatname(DDSplit(getAbsMat()).first, DDSplit(getAbsMat()).second);
   DDMaterial absMatter(absMatname);
   DDLogicalPart glog(DDName(name, idNameSpace), absMatter, solid);
 
   cpv.position(glog, genlogic, 1, DDTranslation(0.0, 0.0, 0.0), DDRotation());
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << glog.name() << " number 1 positioned in " << genlogic.name()
                                << " at (0,0,0) with no rotation";
 #endif
 }

◆ constructInsideModule()

void DDHCalEndcapAlgo::constructInsideModule	(	const DDLogicalPart &	module,
		int	mod,
		DDCompactView &	cpv
	)

protected

Definition at line 1055 of file DDHCalEndcapAlgo.cc.

References simBeamSpotPI::alpha, constructScintLayer(), angle_units::operators::convertRadToDeg(), DDBase< N, C >::ddname(), DDSplit(), dqmdumpme::first, getAngBot(), getAngTop(), getDzStep(), getGenMat(), getLayer(), getLayerN(), getLayerName(), getLayerT(), getLayerType(), getModName(), getNsectors(), getPhi(), getPhiName(), getPlastMat(), getRotMat(), getRoutBlock2(), getScintT(), getSlope(), getThick(), getTrim(), getZ1Beam(), getZminBlock(), mps_fire::i, idNameSpace, idOffset, hcalRecHitTable_cff::iphi, nano_mu_digi_cff::layer, mod(), Skims_PA_cff::name, DDBase< N, C >::name(), parameterLayer(), phi, DDCompactView::position(), diffTwoXMLs::r1, diffTwoXMLs::r2, makeMuonMisalignmentScenario::rot, rotns, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), theta(), and DDSolidFactory::trap().

Referenced by constructInsideSector().

                                                                                                      {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: \t\tInside module ..." << mod;
 #endif
 
   //Pointers to the Rotation Matrices and to the Materials
   std::string rotstr = getRotMat();
   DDRotation rot(DDName(rotstr, rotns));
   DDName matName(DDSplit(getGenMat()).first, DDSplit(getGenMat()).second);
   DDMaterial matter(matName);
   DDName plasName(DDSplit(getPlastMat()).first, DDSplit(getPlastMat()).second);
   DDMaterial matplastic(plasName);
 
   double alpha = (1._pi) / getNsectors();
   double zi = getZminBlock(mod);
 
   for (int i = 0; i < getLayerN(mod); i++) {
     std::string name;
     DDSolid solid;
     DDLogicalPart glog, plog;
     int layer = getLayer(mod, i);
     double zo = zi + 0.5 * getDzStep();
 
     for (int iphi = 0; iphi < getPhi(); iphi++) {
       double ziAir = zo - getThick(mod);
       double rinF, rinB;
       if (layer == 1) {
         rinF = ziAir * tan(getAngTop());
         rinB = zo * tan(getAngTop());
       } else {
         rinF = ziAir * tan(getAngBot());
         rinB = zo * tan(getAngBot());
       }
       double routF = (ziAir - getZ1Beam()) * getSlope();
       double routB = (zo - getZ1Beam()) * getSlope();
       if (routF > getRoutBlock2(mod))
         routF = getRoutBlock2(mod);
       if (routB > getRoutBlock2(mod))
         routB = getRoutBlock2(mod);
 
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Layer " << i << " Phi " << iphi << " Front " << ziAir << ", "
                                    << rinF << ", " << routF << " Back " << zo << ", " << rinB << ", " << routB;
 #endif
 
       double yh1, bl1, tl1, yh2, bl2, tl2, theta, phi, alp;
       double xpos, ypos, zpos;
       parameterLayer(
           iphi, rinF, routF, rinB, routB, ziAir, zo, yh1, bl1, tl1, yh2, bl2, tl2, alp, theta, phi, xpos, ypos, zpos);
 
       name = module.name().name() + getLayerName(layer) + getPhiName(iphi) + "Air";
       solid = DDSolidFactory::trap(
           DDName(name, idNameSpace), 0.5 * getThick(mod), theta, phi, yh1, bl1, tl1, alp, yh2, bl2, tl2, alp);
 
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << solid.name() << " Trap made of " << getGenMat()
                                    << " of dimensions " << 0.5 * getThick(mod) << ", " << convertRadToDeg(theta) << ", "
                                    << convertRadToDeg(phi) << ", " << yh1 << ", " << bl1 << ", " << tl1 << ", "
                                    << convertRadToDeg(alp) << ", " << yh2 << ", " << bl2 << ", " << tl2 << ", "
                                    << convertRadToDeg(alp);
 #endif
 
       glog = DDLogicalPart(solid.ddname(), matter, solid);
 
       DDTranslation r1(xpos, ypos, zpos);
       cpv.position(glog, module, layer + 1, r1, rot);
 
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << glog.name() << " number " << layer + 1
                                    << " positioned in " << module.name() << " at " << r1 << " with " << rot;
 #endif
 
       //Now the plastic with scintillators
       double yh = 0.5 * (routF - rinB) - getTrim(mod, iphi);
       double bl = 0.5 * rinB * tan(alpha) - getTrim(mod, iphi);
       double tl = 0.5 * routF * tan(alpha) - getTrim(mod, iphi);
       name = module.name().name() + getLayerName(layer) + getPhiName(iphi);
       solid = DDSolidFactory::trap(
           DDName(name, idNameSpace), 0.5 * getLayerT(layer), 0, 0, yh, bl, tl, alp, yh, bl, tl, alp);
 
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << solid.name() << " Trap made of " << getPlastMat()
                                    << " of dimensions " << 0.5 * getLayerT(layer) << ", 0, 0, " << yh << ", " << bl
                                    << ", " << tl << ", " << convertRadToDeg(alp) << ", " << yh << ", " << bl << ", "
                                    << tl << ", " << convertRadToDeg(alp);
 #endif
 
       plog = DDLogicalPart(solid.ddname(), matplastic, solid);
 
       ypos = 0.5 * (routF + rinB) - xpos;
       DDTranslation r2(0., ypos, 0.);
       cpv.position(plog, glog, idOffset + layer + 1, r2, DDRotation());
 
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << plog.name() << " number " << idOffset + layer + 1
                                    << " positioned in " << glog.name() << " at " << r2 << " with no rotation";
 #endif
 
       //Constructin the scintillators inside
       int copyNo = layer * 10 + getLayerType(layer);
       name = getModName(mod) + getLayerName(layer) + getPhiName(iphi);
       constructScintLayer(plog, getScintT(layer), yh, bl, tl, alp, name, copyNo, cpv);
       zo += 0.5 * getDzStep();
     }  // End of loop over phi indices
     zi = zo - 0.5 * getDzStep();
   }  // End of loop on layers
 }

◆ constructInsideModule0()

void DDHCalEndcapAlgo::constructInsideModule0	(	const DDLogicalPart &	module,
		int	mod,
		DDCompactView &	cpv
	)

protected

Definition at line 948 of file DDHCalEndcapAlgo.cc.

References constructScintLayer(), angle_units::operators::convertRadToDeg(), DDBase< N, C >::ddname(), DDSplit(), fact, dqmdumpme::first, getAbsMat(), getAngBot(), getAngTop(), getDzStep(), getLayer(), getLayerName(), getLayerT(), getLayerType(), getModName(), getPhi(), getPhiName(), getPlastMat(), getRotMat(), getScintT(), getSlope(), getThick(), getTolAbs(), getZ1Beam(), getZminBlock(), idNameSpace, idOffset, hcalRecHitTable_cff::iphi, nano_mu_digi_cff::layer, mod(), Skims_PA_cff::name, DDBase< N, C >::name(), parameterLayer(), parameterLayer0(), phi, DDCompactView::position(), diffTwoXMLs::r1, diffTwoXMLs::r2, makeMuonMisalignmentScenario::rot, rotns, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), theta(), and DDSolidFactory::trap().

Referenced by constructInsideSector().

                                                                                                       {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: \t\tInside module0 ..." << mod;
 #endif
 
   //Pointers to the Rotation Matrices and to the Materials
   std::string rotstr = getRotMat();
   DDRotation rot(DDName(rotstr, rotns));
   DDName matName(DDSplit(getAbsMat()).first, DDSplit(getAbsMat()).second);
   DDMaterial matabsorbr(matName);
   DDName plasName(DDSplit(getPlastMat()).first, DDSplit(getPlastMat()).second);
   DDMaterial matplastic(plasName);
 
   int layer = getLayer(mod, 0);
   int layer0 = getLayer(mod, 1);
   std::string name;
   double xpos, ypos, zpos;
   DDSolid solid;
   DDLogicalPart glog, plog;
   for (int iphi = 0; iphi < getPhi(); iphi++) {
     double yh, bl, tl, alp;
     parameterLayer0(mod, layer, iphi, yh, bl, tl, alp, xpos, ypos, zpos);
     name = module.name().name() + getLayerName(layer) + getPhiName(iphi);
     solid =
         DDSolidFactory::trap(DDName(name, idNameSpace), 0.5 * getLayerT(layer), 0, 0, yh, bl, tl, alp, yh, bl, tl, alp);
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << solid.name() << " Trap made of " << getPlastMat()
                                  << " of dimensions " << 0.5 * getLayerT(layer) << ", 0, 0, " << yh << ", " << bl
                                  << ", " << tl << ", " << convertRadToDeg(alp) << ", " << yh << ", " << bl << ", " << tl
                                  << ", " << convertRadToDeg(alp);
 #endif
 
     glog = DDLogicalPart(solid.ddname(), matplastic, solid);
 
     DDTranslation r1(xpos, ypos, zpos);
     cpv.position(glog, module, idOffset + layer + 1, r1, rot);
 
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << glog.name() << " number " << idOffset + layer + 1
                                  << " positioned in " << module.name() << " at " << r1 << " with " << rot;
 #endif
 
     //Now construct the layer of scintillator inside this
     int copyNo = layer0 * 10 + getLayerType(layer);
     name = getModName(mod) + getLayerName(layer) + getPhiName(iphi);
     constructScintLayer(glog, getScintT(layer), yh, bl, tl, alp, name, copyNo, cpv);
   }
 
   //Now the absorber layer
   double zi = getZminBlock(mod) + getLayerT(layer);
   double zo = zi + 0.5 * getDzStep();
   double rinF, routF, rinB, routB;
   if (mod == 0) {
     rinF = zi * tan(getAngTop());
     routF = (zi - getZ1Beam()) * getSlope();
     rinB = zo * tan(getAngTop());
     routB = (zo - getZ1Beam()) * getSlope();
   } else {
     rinF = zi * tan(getAngBot());
     routF = zi * tan(getAngTop());
     rinB = zo * tan(getAngBot());
     routB = zo * tan(getAngTop());
   }
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Module " << mod << " Front " << zi << ", " << rinF << ", " << routF
                                << " Back " << zo << ", " << rinB << ", " << routB;
 #endif
 
   double yh1, bl1, tl1, yh2, bl2, tl2, theta, phi, alp;
   parameterLayer(0, rinF, routF, rinB, routB, zi, zo, yh1, bl1, tl1, yh2, bl2, tl2, alp, theta, phi, xpos, ypos, zpos);
   double fact = getTolAbs();
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Trim " << fact << " Param " << yh1 << ", " << bl1 << ", " << tl1
                                << ", " << yh2 << ", " << bl2 << ", " << tl2;
 #endif
 
   bl1 -= fact;
   tl1 -= fact;
   bl2 -= fact;
   tl2 -= fact;
 
   name = module.name().name() + "Absorber";
   solid = DDSolidFactory::trap(
       DDName(name, idNameSpace), 0.5 * getThick(mod), theta, phi, yh1, bl1, tl1, alp, yh2, bl2, tl2, alp);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << solid.name() << " Trap made of " << getAbsMat()
                                << " of dimensions " << 0.5 * getThick(mod) << ", " << convertRadToDeg(theta) << ", "
                                << convertRadToDeg(phi) << ", " << yh1 << ", " << bl1 << ", " << tl1 << ", "
                                << convertRadToDeg(alp) << ", " << yh2 << ", " << bl2 << ", " << tl2 << ", "
                                << convertRadToDeg(alp);
 #endif
 
   glog = DDLogicalPart(solid.ddname(), matabsorbr, solid);
 
   DDTranslation r2(xpos, ypos, zpos);
   cpv.position(glog, module, 1, r2, rot);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << glog.name() << " number 1 positioned in " << module.name()
                                << " at " << r2 << " with " << rot;
 #endif
 }

◆ constructInsideSector()

void DDHCalEndcapAlgo::constructInsideSector	(	const DDLogicalPart &	sector,
		DDCompactView &	cpv
	)

protected

!!!!!!!!!!!!!!!!Should be zero. And removed as soon as

Definition at line 758 of file DDHCalEndcapAlgo.cc.

References simBeamSpotPI::alpha, constructInsideModule(), constructInsideModule0(), angle_units::operators::convertRadToDeg(), DDSplit(), equipModule(), dqmdumpme::first, getAngBot(), getModMat(), getModName(), getModType(), getModules(), getNsectors(), getRinBlock1(), getRinBlock2(), getRinDip(), getRinKink(), getRout(), getRoutBlock1(), getRoutBlock2(), getRoutDip(), getSectionModule(), getSlope(), getZ1Beam(), getZiDip(), getZiKink(), getZmaxBlock(), getZminBlock(), getZShiftHac2(), mps_fire::i, idName, idNameSpace, dqmdumpme::k, Skims_PA_cff::name, DDBase< N, C >::name(), DDSolidFactory::polyhedra(), DDCompactView::position(), edm::second(), nano_mu_digi_cff::sector, AlCaHLTBitMon_QueryRunRegistry::string, and funct::tan().

Referenced by constructGeneralVolume().

                                                                                             {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Modules (" << getModules() << ") ...";
 #endif
 
   double alpha = (1._pi) / getNsectors();
 
   for (int i = 0; i < getModules(); i++) {
     std::string name = idName + getModName(i);
     DDName matname(DDSplit(getModMat(i)).first, DDSplit(getModMat(i)).second);
     DDMaterial matter(matname);
 
     if (equipModule(i) > 0) {
       int nsec = getSectionModule(i);
 
       //vertical walls are allowed in SolidPolyhedra
       double deltaz = 0;
 
       std::vector<double> pgonZ, pgonRmin, pgonRmax;
       if (nsec == 3) {
         double zf = getZminBlock(i) + getZShiftHac2();
         pgonZ.emplace_back(zf);
         pgonRmin.emplace_back(zf * tan(getAngBot()));
         pgonRmax.emplace_back((zf - getZ1Beam()) * getSlope());
         pgonZ.emplace_back(getZiKink());
         pgonRmin.emplace_back(getRinKink());
         pgonRmax.emplace_back(getRout());
       } else {
         pgonZ.emplace_back(getZminBlock(i));
         pgonRmin.emplace_back(getRinBlock1(i));
         pgonRmax.emplace_back(getRoutBlock1(i));
       }
       if (nsec == 4) {
         pgonZ.emplace_back(getZiDip());
         pgonRmin.emplace_back(getRinDip());
         pgonRmax.emplace_back(getRout());
         pgonZ.emplace_back(pgonZ[1] + deltaz);
         pgonRmin.emplace_back(pgonRmin[1]);
         pgonRmax.emplace_back(getRoutDip());
       }
       pgonZ.emplace_back(getZmaxBlock(i));
       pgonRmin.emplace_back(getRinBlock2(i));
       pgonRmax.emplace_back(getRoutBlock2(i));
 
       //Solid & volume
       DDSolid solid;
       solid = DDSolidFactory::polyhedra(DDName(name, idNameSpace), 1, -alpha, 2 * alpha, pgonZ, pgonRmin, pgonRmax);
 
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << DDName(name, idNameSpace) << " Polyhedra made of "
                                    << getModMat(i) << " with 1 sector from " << convertRadToDeg(-alpha) << " to "
                                    << convertRadToDeg(alpha) << " and with " << nsec << " sections";
       for (unsigned int k = 0; k < pgonZ.size(); k++)
         edm::LogVerbatim("HCalGeom") << "\t\tZ = " << pgonZ[k] << "\tRmin = " << pgonRmin[k]
                                      << "\tRmax = " << pgonRmax[k];
 #endif
 
       DDLogicalPart glog(DDName(name, idNameSpace), matter, solid);
 
       cpv.position(glog, sector, i + 1, DDTranslation(0.0, 0.0, 0.0), DDRotation());
 
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << glog.name() << " number " << i + 1 << " positioned in "
                                    << sector.name() << " at (0,0,0) with no rotation";
 #endif
 
       if (getModType(i) == 0)
         constructInsideModule0(glog, i, cpv);
       else
         constructInsideModule(glog, i, cpv);
     }
   }
 }

◆ constructScintLayer()

void DDHCalEndcapAlgo::constructScintLayer	(	const DDLogicalPart &	glog,
		double	pDz,
		double	yh,
		double	bl,
		double	tl,
		double	alp,
		const std::string &	name,
		int	id,
		DDCompactView &	cpv
	)

protected

Definition at line 1164 of file DDHCalEndcapAlgo.cc.

References angle_units::operators::convertRadToDeg(), DDBase< N, C >::ddname(), DDSplit(), hgcalTestNeighbor_cfi::detector, PVValHelper::dz, dqmdumpme::first, getScintMat(), idName, idNameSpace, Skims_PA_cff::name, DDCompactView::position(), edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, and DDSolidFactory::trap().

Referenced by constructInsideModule(), and constructInsideModule0().

                                                                {
   DDName matname(DDSplit(getScintMat()).first, DDSplit(getScintMat()).second);
   DDMaterial matter(matname);
   std::string name = idName + "Scintillator" + nm;
 
   DDSolid solid = DDSolidFactory::trap(DDName(name, idNameSpace), 0.5 * dz, 0, 0, yh, bl, tl, alp, yh, bl, tl, alp);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << DDName(name, idNameSpace) << " Trap made of " << getScintMat()
                                << " of dimensions " << 0.5 * dz << ", 0, 0, " << yh << ", " << bl << ", " << tl << ", "
                                << convertRadToDeg(alp) << ", " << yh << ", " << bl << ", " << tl << ", "
                                << convertRadToDeg(alp);
 #endif
 
   DDLogicalPart glog(solid.ddname(), matter, solid);
 
   cpv.position(glog, detector, id, DDTranslation(0, 0, 0), DDRotation());
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << glog.name() << " number " << id << " positioned in "
                                << detector.name() << " at (0,0,0) with no rotation";
 #endif
 }

◆ equipModule()

int DDHCalEndcapAlgo::equipModule ( unsigned int i ) const

inline

Definition at line 41 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructGeneralVolume(), and constructInsideSector().

41 { return eModule[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::eModule

std::vector< int > eModule

Definition: DDHCalEndcapAlgo.cc:156

◆ execute()

void DDHCalEndcapAlgo::execute ( DDCompactView & cpv )

override

Definition at line 541 of file DDHCalEndcapAlgo.cc.

References constructGeneralVolume().

                                                  {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "==>> Constructing DDHCalEndcapAlgo...";
 #endif
 
   constructGeneralVolume(cpv);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "<<== End of DDHCalEndcapAlgo construction ...";
 #endif
 }

◆ getAbsMat()

const std::string& DDHCalEndcapAlgo::getAbsMat ( ) const

inline

Definition at line 69 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and constructInsideModule0().

69 { return absMat; }

DDHCalEndcapAlgo::absMat

std::string absMat

Definition: DDHCalEndcapAlgo.cc:188

◆ getAngBot()

double DDHCalEndcapAlgo::getAngBot ( ) const

inline

Definition at line 65 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), constructInsideModule(), constructInsideModule0(), constructInsideSector(), and parameterLayer0().

65 { return angBot; }

DDHCalEndcapAlgo::angBot

double angBot

Definition: DDHCalEndcapAlgo.cc:184

◆ getAngGap()

double DDHCalEndcapAlgo::getAngGap ( ) const

inline

Definition at line 66 of file DDHCalEndcapAlgo.cc.

66 { return angGap; }

DDHCalEndcapAlgo::angGap

double angGap

Definition: DDHCalEndcapAlgo.cc:185

◆ getAngTop()

double DDHCalEndcapAlgo::getAngTop ( ) const

inline

Definition at line 64 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), constructInsideModule(), constructInsideModule0(), and parameterLayer0().

64 { return angTop; }

DDHCalEndcapAlgo::angTop

double angTop

Definition: DDHCalEndcapAlgo.cc:183

◆ getDrEnd()

double DDHCalEndcapAlgo::getDrEnd ( ) const

inline

Definition at line 63 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume().

63 { return drEnd; }

DDHCalEndcapAlgo::drEnd

double drEnd

Definition: DDHCalEndcapAlgo.cc:181

◆ getDzShift()

double DDHCalEndcapAlgo::getDzShift ( ) const

inline

Definition at line 55 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume().

55 { return dzShift; }

DDHCalEndcapAlgo::dzShift

double dzShift

Definition: DDHCalEndcapAlgo.cc:173

◆ getDzStep()

double DDHCalEndcapAlgo::getDzStep ( ) const

inline

Definition at line 54 of file DDHCalEndcapAlgo.cc.

Referenced by constructInsideModule(), and constructInsideModule0().

54 { return dzStep; }

DDHCalEndcapAlgo::dzStep

double dzStep

Definition: DDHCalEndcapAlgo.cc:172

◆ getEndcaps()

int DDHCalEndcapAlgo::getEndcaps ( ) const

inline

Definition at line 40 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume().

40 { return nEndcap; }

DDHCalEndcapAlgo::nEndcap

int nEndcap

Definition: DDHCalEndcapAlgo.cc:155

◆ getGenMat()

const std::string& DDHCalEndcapAlgo::getGenMat ( ) const

inline

Definition at line 36 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and constructInsideModule().

36 { return genMaterial; }

DDHCalEndcapAlgo::genMaterial

std::string genMaterial

Definition: DDHCalEndcapAlgo.cc:152

◆ getHeboxDepth()

double DDHCalEndcapAlgo::getHeboxDepth ( ) const

inline

Definition at line 62 of file DDHCalEndcapAlgo.cc.

62 { return heboxDepth; }

DDHCalEndcapAlgo::heboxDepth

double heboxDepth

Definition: DDHCalEndcapAlgo.cc:180

◆ getLayer()

int DDHCalEndcapAlgo::getLayer	(	unsigned int	i,
		unsigned int	j
	)		const

Definition at line 259 of file DDHCalEndcapAlgo.cc.

References mps_fire::i, dqmiolumiharvest::j, layerN0, layerN1, layerN2, layerN3, layerN4, and layerN5.

Referenced by constructInsideModule(), and constructInsideModule0().

                                                                    {
   switch (i) {
     case 0:
       return layerN0[j];
       break;
 
     case 1:
       return layerN1[j];
       break;
 
     case 2:
       return layerN2[j];
       break;
 
     case 3:
       return layerN3[j];
       break;
 
     case 4:
       return layerN4[j];
       break;
 
     case 5:
       return layerN5[j];
       break;
 
     default:
       return 0;
   }
 }

◆ getLayerN()

int DDHCalEndcapAlgo::getLayerN ( unsigned int i ) const

inline

Definition at line 75 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideModule().

75 { return layerN[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::layerN

std::vector< int > layerN

Definition: DDHCalEndcapAlgo.cc:194

◆ getLayerName()

const std::string& DDHCalEndcapAlgo::getLayerName ( unsigned int i ) const

inline

Definition at line 89 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideModule(), and constructInsideModule0().

89 { return layerName[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::layerName

std::vector< std::string > layerName

Definition: DDHCalEndcapAlgo.cc:214

◆ getLayers()

int DDHCalEndcapAlgo::getLayers ( ) const

inline

Definition at line 88 of file DDHCalEndcapAlgo.cc.

References hgcalTBTopologyTester_cfi::layers.

88 { return layers; }

DDHCalEndcapAlgo::layers

int layers

Definition: DDHCalEndcapAlgo.cc:213

◆ getLayerT()

double DDHCalEndcapAlgo::getLayerT ( unsigned int i ) const

inline

Definition at line 91 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideModule(), constructInsideModule0(), and parameterLayer0().

91 { return layerT[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::layerT

std::vector< double > layerT

Definition: DDHCalEndcapAlgo.cc:216

◆ getLayerType()

int DDHCalEndcapAlgo::getLayerType ( unsigned int i ) const

inline

Definition at line 90 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideModule(), and constructInsideModule0().

90 { return layerType[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::layerType

std::vector< int > layerType

Definition: DDHCalEndcapAlgo.cc:215

◆ getModMat()

const std::string& DDHCalEndcapAlgo::getModMat ( unsigned int i ) const

inline

Definition at line 72 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideSector().

72 { return modMat[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::modMat

std::vector< std::string > modMat

Definition: DDHCalEndcapAlgo.cc:191

◆ getModName()

const std::string& DDHCalEndcapAlgo::getModName ( unsigned int i ) const

inline

Definition at line 71 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideModule(), constructInsideModule0(), and constructInsideSector().

71 { return modName[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::modName

std::vector< std::string > modName

Definition: DDHCalEndcapAlgo.cc:190

◆ getModType()

int DDHCalEndcapAlgo::getModType ( unsigned int i ) const

inline

Definition at line 73 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideSector().

73 { return modType[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::modType

std::vector< int > modType

Definition: DDHCalEndcapAlgo.cc:192

◆ getModules()

int DDHCalEndcapAlgo::getModules ( ) const

inline

Definition at line 70 of file DDHCalEndcapAlgo.cc.

References LogMessageMonitor_cff::modules.

Referenced by constructInsideSector().

70 { return modules; }

DDHCalEndcapAlgo::modules

int modules

Definition: DDHCalEndcapAlgo.cc:189

◆ getNsectors()

int DDHCalEndcapAlgo::getNsectors ( ) const

inline

Definition at line 38 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), constructInsideModule(), constructInsideSector(), parameterLayer(), and parameterLayer0().

38 { return nsectors; }

DDHCalEndcapAlgo::nsectors

int nsectors

Definition: DDHCalEndcapAlgo.cc:153

◆ getNsectortot()

int DDHCalEndcapAlgo::getNsectortot ( ) const

inline

Definition at line 39 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume().

39 { return nsectortot; }

DDHCalEndcapAlgo::nsectortot

int nsectortot

Definition: DDHCalEndcapAlgo.cc:154

◆ getPhi()

int DDHCalEndcapAlgo::getPhi ( ) const

inline

Definition at line 86 of file DDHCalEndcapAlgo.cc.

Referenced by constructInsideModule(), and constructInsideModule0().

86 { return phiSections; }

DDHCalEndcapAlgo::phiSections

int phiSections

Definition: DDHCalEndcapAlgo.cc:211

◆ getPhiName()

const std::string& DDHCalEndcapAlgo::getPhiName ( unsigned int i ) const

inline

Definition at line 87 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideModule(), and constructInsideModule0().

87 { return phiName[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::phiName

std::vector< std::string > phiName

Definition: DDHCalEndcapAlgo.cc:212

◆ getPlastMat()

const std::string& DDHCalEndcapAlgo::getPlastMat ( ) const

inline

Definition at line 93 of file DDHCalEndcapAlgo.cc.

Referenced by constructInsideModule(), and constructInsideModule0().

93 { return plastMat; }

DDHCalEndcapAlgo::plastMat

std::string plastMat

Definition: DDHCalEndcapAlgo.cc:218

◆ getRinBlock1()

double DDHCalEndcapAlgo::getRinBlock1 ( unsigned i ) const

inline

Definition at line 81 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideSector().

81 { return rinBlock1[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::rinBlock1

std::vector< double > rinBlock1

Definition: DDHCalEndcapAlgo.cc:206

◆ getRinBlock2()

double DDHCalEndcapAlgo::getRinBlock2 ( unsigned i ) const

inline

Definition at line 82 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideSector().

82 { return rinBlock2[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::rinBlock2

std::vector< double > rinBlock2

Definition: DDHCalEndcapAlgo.cc:208

◆ getRinDip()

double DDHCalEndcapAlgo::getRinDip ( ) const

inline

Definition at line 60 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and constructInsideSector().

60 { return riDip; }

DDHCalEndcapAlgo::riDip

double riDip

Definition: DDHCalEndcapAlgo.cc:178

◆ getRinKink()

double DDHCalEndcapAlgo::getRinKink ( ) const

inline

Definition at line 59 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and constructInsideSector().

59 { return riKink; }

DDHCalEndcapAlgo::riKink

double riKink

Definition: DDHCalEndcapAlgo.cc:177

◆ getRotation()

const std::string& DDHCalEndcapAlgo::getRotation ( ) const

inline

Definition at line 37 of file DDHCalEndcapAlgo.cc.

References idealTransformation::rotation.

Referenced by constructGeneralVolume().

37 { return rotation; }

DDHCalEndcapAlgo::rotation

std::string rotation

Definition: DDHCalEndcapAlgo.cc:159

◆ getRotMat()

const std::string& DDHCalEndcapAlgo::getRotMat ( ) const

inline

Definition at line 95 of file DDHCalEndcapAlgo.cc.

Referenced by constructInsideModule(), and constructInsideModule0().

95 { return rotmat; }

DDHCalEndcapAlgo::rotmat

std::string rotmat

Definition: DDHCalEndcapAlgo.cc:220

◆ getRout()

double DDHCalEndcapAlgo::getRout ( ) const

inline

Definition at line 58 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and constructInsideSector().

58 { return rout; }

DDHCalEndcapAlgo::rout

double rout

Definition: DDHCalEndcapAlgo.cc:176

◆ getRoutBlock1()

double DDHCalEndcapAlgo::getRoutBlock1 ( unsigned i ) const

inline

Definition at line 83 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideSector().

83 { return routBlock1[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::routBlock1

std::vector< double > routBlock1

Definition: DDHCalEndcapAlgo.cc:207

◆ getRoutBlock2()

double DDHCalEndcapAlgo::getRoutBlock2 ( unsigned i ) const

inline

Definition at line 84 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideModule(), and constructInsideSector().

84 { return routBlock2[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::routBlock2

std::vector< double > routBlock2

Definition: DDHCalEndcapAlgo.cc:209

◆ getRoutDip()

double DDHCalEndcapAlgo::getRoutDip ( ) const

inline

Definition at line 61 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and constructInsideSector().

61 { return roDip; }

DDHCalEndcapAlgo::roDip

double roDip

Definition: DDHCalEndcapAlgo.cc:179

◆ getScintMat()

const std::string& DDHCalEndcapAlgo::getScintMat ( ) const

inline

Definition at line 94 of file DDHCalEndcapAlgo.cc.

Referenced by constructScintLayer().

94 { return scintMat; }

DDHCalEndcapAlgo::scintMat

std::string scintMat

Definition: DDHCalEndcapAlgo.cc:219

◆ getScintT()

double DDHCalEndcapAlgo::getScintT ( unsigned int i ) const

inline

Definition at line 92 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideModule(), and constructInsideModule0().

92 { return scintT[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::scintT

std::vector< double > scintT

Definition: DDHCalEndcapAlgo.cc:217

◆ getSectionModule()

int DDHCalEndcapAlgo::getSectionModule ( unsigned i ) const

inline

Definition at line 74 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideSector().

74 { return sectionModule[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::sectionModule

std::vector< int > sectionModule

Definition: DDHCalEndcapAlgo.cc:193

◆ getSlope()

double DDHCalEndcapAlgo::getSlope ( ) const

inline

Definition at line 67 of file DDHCalEndcapAlgo.cc.

References slope.

Referenced by constructGeneralVolume(), constructInsideModule(), constructInsideModule0(), constructInsideSector(), and parameterLayer0().

67 { return slope; }

DDHCalEndcapAlgo::slope

double slope

Definition: DDHCalEndcapAlgo.cc:186

◆ getThick()

double DDHCalEndcapAlgo::getThick ( unsigned int i ) const

inline

Definition at line 77 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideModule(), and constructInsideModule0().

77 { return thick[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::thick

std::vector< double > thick

Definition: DDHCalEndcapAlgo.cc:201

◆ getTolAbs()

double DDHCalEndcapAlgo::getTolAbs ( ) const

inline

Definition at line 97 of file DDHCalEndcapAlgo.cc.

Referenced by constructInsideModule0().

97 { return tolAbs; }

DDHCalEndcapAlgo::tolAbs

double tolAbs

Definition: DDHCalEndcapAlgo.cc:226

◆ getTolPos()

double DDHCalEndcapAlgo::getTolPos ( ) const

inline

Definition at line 96 of file DDHCalEndcapAlgo.cc.

96 { return tolPos; }

DDHCalEndcapAlgo::tolPos

double tolPos

Definition: DDHCalEndcapAlgo.cc:226

◆ getTrim()

double DDHCalEndcapAlgo::getTrim	(	unsigned int	i,
		unsigned int	j
	)		const

Definition at line 290 of file DDHCalEndcapAlgo.cc.

References mps_fire::i, dqmiolumiharvest::j, trimLeft, and trimRight.

Referenced by constructInsideModule(), and parameterLayer0().

                                                                      {
   if (j == 0)
     return trimLeft[i];
   else
     return trimRight[j];
 }

◆ getZ0Beam()

double DDHCalEndcapAlgo::getZ0Beam ( ) const

inline

Definition at line 51 of file DDHCalEndcapAlgo.cc.

51 { return z0Beam; }

DDHCalEndcapAlgo::z0Beam

double z0Beam

Definition: DDHCalEndcapAlgo.cc:169

◆ getZ1Beam()

double DDHCalEndcapAlgo::getZ1Beam ( ) const

inline

Definition at line 52 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), constructInsideModule(), constructInsideModule0(), constructInsideSector(), and parameterLayer0().

52 { return z1Beam; }

DDHCalEndcapAlgo::z1Beam

double z1Beam

Definition: DDHCalEndcapAlgo.cc:170

◆ getZEnd()

double DDHCalEndcapAlgo::getZEnd ( ) const

inline

Definition at line 45 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume().

45 { return zEnd; }

DDHCalEndcapAlgo::zEnd

double zEnd

Definition: DDHCalEndcapAlgo.cc:163

◆ getZFront()

double DDHCalEndcapAlgo::getZFront ( ) const

inline

Definition at line 44 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume().

44 { return zFront; }

DDHCalEndcapAlgo::zFront

double zFront

Definition: DDHCalEndcapAlgo.cc:162

◆ getZiBody()

double DDHCalEndcapAlgo::getZiBody ( ) const

inline

Definition at line 48 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume().

48 { return ziBody; }

DDHCalEndcapAlgo::ziBody

double ziBody

Definition: DDHCalEndcapAlgo.cc:166

◆ getZiDip()

double DDHCalEndcapAlgo::getZiDip ( ) const

inline

Definition at line 53 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and constructInsideSector().

53 { return ziDip; }

DDHCalEndcapAlgo::ziDip

double ziDip

Definition: DDHCalEndcapAlgo.cc:171

◆ getZiKink()

double DDHCalEndcapAlgo::getZiKink ( ) const

inline

Definition at line 50 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and constructInsideSector().

50 { return ziKink; }

DDHCalEndcapAlgo::ziKink

double ziKink

Definition: DDHCalEndcapAlgo.cc:168

◆ getZiL0Body()

double DDHCalEndcapAlgo::getZiL0Body ( ) const

inline

Definition at line 49 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume().

49 { return ziL0Body; }

DDHCalEndcapAlgo::ziL0Body

double ziL0Body

Definition: DDHCalEndcapAlgo.cc:167

◆ getZiL0Nose()

double DDHCalEndcapAlgo::getZiL0Nose ( ) const

inline

Definition at line 47 of file DDHCalEndcapAlgo.cc.

47 { return ziL0Nose; }

DDHCalEndcapAlgo::ziL0Nose

double ziL0Nose

Definition: DDHCalEndcapAlgo.cc:165

◆ getZiNose()

double DDHCalEndcapAlgo::getZiNose ( ) const

inline

Definition at line 46 of file DDHCalEndcapAlgo.cc.

46 { return ziNose; }

DDHCalEndcapAlgo::ziNose

double ziNose

Definition: DDHCalEndcapAlgo.cc:164

◆ getZmaxBlock()

double DDHCalEndcapAlgo::getZmaxBlock ( unsigned i ) const

inline

Definition at line 80 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideSector(), and parameterLayer0().

80 { return zmaxBlock[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::zmaxBlock

std::vector< double > zmaxBlock

Definition: DDHCalEndcapAlgo.cc:205

◆ getZminBlock()

double DDHCalEndcapAlgo::getZminBlock ( unsigned i ) const

inline

Definition at line 79 of file DDHCalEndcapAlgo.cc.

References mps_fire::i.

Referenced by constructInsideModule(), constructInsideModule0(), constructInsideSector(), and parameterLayer0().

79 { return zminBlock[i]; }

mps_fire.i

i

Definition: mps_fire.py:429

DDHCalEndcapAlgo::zminBlock

std::vector< double > zminBlock

Definition: DDHCalEndcapAlgo.cc:204

◆ getZShift()

double DDHCalEndcapAlgo::getZShift ( ) const

inline

Definition at line 42 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume().

42 { return zShift; }

DDHCalEndcapAlgo::zShift

double zShift

Definition: DDHCalEndcapAlgo.cc:160

◆ getZShiftHac2()

double DDHCalEndcapAlgo::getZShiftHac2 ( ) const

inline

Definition at line 56 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and constructInsideSector().

56 { return zShiftHac2; }

DDHCalEndcapAlgo::zShiftHac2

double zShiftHac2

Definition: DDHCalEndcapAlgo.cc:174

◆ initialize()

void DDHCalEndcapAlgo::initialize	(	const DDNumericArguments &	nArgs,
		const DDVectorArguments &	vArgs,
		const DDMapArguments &	mArgs,
		const DDStringArguments &	sArgs,
		const DDStringVectorArguments &	vsArgs
	)

override

Definition at line 297 of file DDHCalEndcapAlgo.cc.

References absMat, angBot, angGap, angTop, angle_units::operators::convertRadToDeg(), dbl_to_int(), drEnd, dzShift, dzStep, eModule, TrackingDataMCValidation_Standalone_cff::etamin, JetChargeProducer_cfi::exp, HLT_2023v12_cff::gap, genMaterial, heboxDepth, mps_fire::i, idName, idNameSpace, idOffset, createfilelist::int, dqmiolumiharvest::j, layerN, layerN0, layerN1, layerN2, layerN3, layerN4, layerN5, layerName, layers, layerT, layerType, modMat, modName, modType, callgraph::module, modules, nEndcap, DDCurrentNamespace::ns(), nsectors, nsectortot, class-composition::parent, phiName, phiSections, plastMat, diffTwoXMLs::r1, riDip, riKink, rinBlock1, rinBlock2, roDip, rotation, rotHalf, rotmat, rotns, rout, routBlock1, routBlock2, scintMat, scintT, sectionModule, funct::sin(), slope, funct::tan(), thick, tolAbs, tolPos, trimLeft, trimRight, z0Beam, z1Beam, zEnd, zFront, ziBody, ziDip, ziKink, ziL0Body, ziL0Nose, ziNose, zmaxBlock, zminBlock, zShift, and zShiftHac2.

                                                                          {
   genMaterial = sArgs["MaterialName"];
   rotation = sArgs["Rotation"];
   nsectors = int(nArgs["Sector"]);
   nsectortot = int(nArgs["SectorTot"]);
   nEndcap = int(nArgs["Endcap"]);
   rotHalf = sArgs["RotHalf"];
   rotns = sArgs["RotNameSpace"];
   zShift = nArgs["ZShift"];
 
   zFront = nArgs["ZFront"];
   zEnd = nArgs["ZEnd"];
   ziNose = nArgs["ZiNose"];
   ziL0Nose = nArgs["ZiL0Nose"];
   ziBody = nArgs["ZiBody"];
   ziL0Body = nArgs["ZiL0Body"];
   z0Beam = nArgs["Z0Beam"];
   ziDip = nArgs["ZiDip"];
   dzStep = nArgs["DzStep"];
   zShiftHac2 = nArgs["ZShiftHac2"];
   double gap = nArgs["Gap"];
   double z1 = nArgs["Z1"];
   double r1 = nArgs["R1"];
   rout = nArgs["Rout"];
   heboxDepth = nArgs["HEboxDepth"];
   drEnd = nArgs["DrEnd"];
   double etamin = nArgs["Etamin"];
   angBot = nArgs["AngBot"];
   angGap = nArgs["AngGap"];
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: General material " << genMaterial << "\tSectors " << nsectors
                                << ",  " << nsectortot << "\tEndcaps " << nEndcap << "\tRotation matrix for half "
                                << rotns << ":" << rotHalf << "\n\tzFront " << zFront << " zEnd " << zEnd << " ziNose "
                                << ziNose << " ziL0Nose " << ziL0Nose << " ziBody " << ziBody << " ziL0Body " << ziL0Body
                                << " z0Beam " << z0Beam << " ziDip " << ziDip << " dzStep " << dzStep << " Gap " << gap
                                << " z1 " << z1 << "\n\tr1 " << r1 << " rout " << rout << " HeboxDepth " << heboxDepth
                                << " drEnd " << drEnd << "\tetamin " << etamin << " Bottom angle " << angBot
                                << " Gap angle " << angGap << " Z-Shift " << zShift << " " << zShiftHac2;
 #endif
 
   //Derived quantities
   angTop = 2.0 * atan(exp(-etamin));
   slope = tan(angGap);
   z1Beam = z1 - r1 / slope;
   ziKink = z1Beam + rout / slope;
   riKink = ziKink * tan(angBot);
   riDip = ziDip * tan(angBot);
   roDip = rout - heboxDepth;
   dzShift = (z1Beam - z0Beam) - gap / sin(angGap);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: angTop " << convertRadToDeg(angTop) << "\tSlope " << slope
                                << "\tDzShift " << dzShift << "\n\tz1Beam " << z1Beam << "\tziKink" << ziKink
                                << "\triKink " << riKink << "\triDip " << riDip << "\n\troDip " << roDip << "\tRotation "
                                << rotation;
 #endif
 
   //Modules
   absMat = sArgs["AbsMat"];
   modules = int(nArgs["Modules"]);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Number of modules " << modules << " and absorber material "
                                << absMat;
 #endif
 
   modName = vsArgs["ModuleName"];
   modMat = vsArgs["ModuleMat"];
   modType = dbl_to_int(vArgs["ModuleType"]);
   sectionModule = dbl_to_int(vArgs["SectionModule"]);
   thick = vArgs["ModuleThick"];
   trimLeft = vArgs["TrimLeft"];
   trimRight = vArgs["TrimRight"];
   eModule = dbl_to_int(vArgs["EquipModule"]);
   layerN = dbl_to_int(vArgs["LayerN"]);
   layerN0 = dbl_to_int(vArgs["LayerN0"]);
   layerN1 = dbl_to_int(vArgs["LayerN1"]);
   layerN2 = dbl_to_int(vArgs["LayerN2"]);
   layerN3 = dbl_to_int(vArgs["LayerN3"]);
   layerN4 = dbl_to_int(vArgs["LayerN4"]);
   layerN5 = dbl_to_int(vArgs["LayerN5"]);
 
 #ifdef EDM_ML_DEBUG
   for (int i = 0; i < modules; i++) {
     edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: " << modName[i] << " type " << modType[i] << " Sections "
                                  << sectionModule[i] << " thickness of absorber/air " << thick[i] << " trim "
                                  << trimLeft[i] << ", " << trimRight[i] << " equip module " << eModule[i] << " with "
                                  << layerN[i] << " layers";
     if (i == 0) {
       for (int j = 0; j < layerN[i]; j++) {
         edm::LogVerbatim("HCalGeom") << "\t " << layerN0[j] << "/" << layerN0[j + 1];
       }
     } else if (i == 1) {
       for (int j = 0; j < layerN[i]; j++) {
         edm::LogVerbatim("HCalGeom") << "\t " << layerN1[j] << "/" << layerN1[j + 1];
       }
     } else if (i == 2) {
       for (int j = 0; j < layerN[i]; j++) {
         edm::LogVerbatim("HCalGeom") << "\t " << layerN2[j];
       }
     } else if (i == 3) {
       for (int j = 0; j < layerN[i]; j++) {
         edm::LogVerbatim("HCalGeom") << "\t " << layerN3[j];
       }
     } else if (i == 4) {
       for (int j = 0; j < layerN[i]; j++) {
         edm::LogVerbatim("HCalGeom") << "\t " << layerN4[j];
       }
     } else if (i == 5) {
       for (int j = 0; j < layerN[i]; j++) {
         edm::LogVerbatim("HCalGeom") << "\t " << layerN5[j];
       }
     }
   }
 #endif
 
   //Layers
   phiSections = int(nArgs["PhiSections"]);
   phiName = vsArgs["PhiName"];
   layers = int(nArgs["Layers"]);
   layerName = vsArgs["LayerName"];
   layerType = dbl_to_int(vArgs["LayerType"]);
   layerT = vArgs["LayerT"];
   scintT = vArgs["ScintT"];
   scintMat = sArgs["ScintMat"];
   plastMat = sArgs["PlastMat"];
   rotmat = sArgs["RotMat"];
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Phi Sections " << phiSections;
   for (int i = 0; i < phiSections; i++)
     edm::LogVerbatim("HCalGeom") << "\tName[" << i << "] : " << phiName[i];
   edm::LogVerbatim("HCalGeom") << "\tPlastic: " << plastMat << "\tScintillator: " << scintMat << "\tRotation matrix "
                                << rotns << ":" << rotmat << "\n\tNumber of layers " << layers;
   for (int i = 0; i < layers; i++) {
     edm::LogVerbatim("HCalGeom") << "\t" << layerName[i] << "\tType " << layerType[i] << "\tThickness " << layerT[i]
                                  << "\tScint.Thick " << scintT[i];
   }
 #endif
 
   // Derive bounding of the modules
   int module = 0;
   // Layer 0 (Nose)
   if (modules > 0) {
     zminBlock.emplace_back(ziL0Nose);
     zmaxBlock.emplace_back(zminBlock[module] + layerT[0] + 0.5 * dzStep);
     rinBlock1.emplace_back(zminBlock[module] * tan(angTop));
     rinBlock2.emplace_back(zmaxBlock[module] * tan(angTop));
     routBlock1.emplace_back((zminBlock[module] - z1Beam) * slope);
     routBlock2.emplace_back((zmaxBlock[module] - z1Beam) * slope);
     module++;
   }
 
   // Layer 0 (Body)
   if (modules > 1) {
     zminBlock.emplace_back(ziL0Body);
     zmaxBlock.emplace_back(zminBlock[module] + layerT[0] + 0.5 * dzStep);
     rinBlock1.emplace_back(zminBlock[module] * tan(angBot));
     rinBlock2.emplace_back(zmaxBlock[module] * tan(angBot));
     routBlock1.emplace_back(zminBlock[module] * tan(angTop));
     routBlock2.emplace_back(zmaxBlock[module] * tan(angTop));
     module++;
   }
 
   // Hac1
   if (modules > 2) {
     zminBlock.emplace_back(ziNose);
     zmaxBlock.emplace_back(ziBody);
     rinBlock1.emplace_back(zminBlock[module] * tan(angTop));
     rinBlock2.emplace_back(zmaxBlock[module] * tan(angTop));
     routBlock1.emplace_back((zminBlock[module] - z1Beam) * slope);
     routBlock2.emplace_back((zmaxBlock[module] - z1Beam) * slope);
     module++;
   }
 
   // Hac2
   if (modules > 3) {
     zminBlock.emplace_back(ziBody);
     zmaxBlock.emplace_back(zminBlock[module] + layerN[3] * dzStep);
     rinBlock1.emplace_back(zminBlock[module] * tan(angBot));
     rinBlock2.emplace_back(zmaxBlock[module] * tan(angBot));
     routBlock1.emplace_back((zmaxBlock[module - 1] - z1Beam) * slope);
     routBlock2.emplace_back(rout);
     module++;
   }
 
   // Hac3
   if (modules > 4) {
     zminBlock.emplace_back(zmaxBlock[module - 1]);
     zmaxBlock.emplace_back(zminBlock[module] + layerN[4] * dzStep);
     rinBlock1.emplace_back(zminBlock[module] * tan(angBot));
     rinBlock2.emplace_back(zmaxBlock[module] * tan(angBot));
     routBlock1.emplace_back(rout);
     routBlock2.emplace_back(rout);
     module++;
   }
 
   // Hac4
   if (modules > 5) {
     zminBlock.emplace_back(zmaxBlock[module - 1]);
     zmaxBlock.emplace_back(zminBlock[module] + layerN[5] * dzStep);
     rinBlock1.emplace_back(zminBlock[module] * tan(angBot));
     rinBlock2.emplace_back(zmaxBlock[module] * tan(angBot));
     routBlock1.emplace_back(rout);
     routBlock2.emplace_back(roDip);
     module++;
   }
 
 #ifdef EDM_ML_DEBUG
   for (int i = 0; i < module; i++)
     edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Module " << i << "\tZ/Rin/Rout " << zminBlock[i] << ", "
                                  << zmaxBlock[i] << "/ " << rinBlock1[i] << ", " << rinBlock2[i] << "/ "
                                  << routBlock1[i] << ", " << routBlock2[i];
 #endif
 
   idName = sArgs["MotherName"];
   idNameSpace = DDCurrentNamespace::ns();
   idOffset = int(nArgs["IdOffset"]);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Parent " << parent().name() << " idName " << idName
                                << " NameSpace " << idNameSpace << " Offset " << idOffset;
 #endif
 
   tolPos = nArgs["TolPos"];
   tolAbs = nArgs["TolAbs"];
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "DDHCalEndcapAlgo: Tolerances - Positioning " << tolPos << " Absorber " << tolAbs;
 #endif
 }

◆ parameterLayer()

void DDHCalEndcapAlgo::parameterLayer	(	int	iphi,
		double	rinF,
		double	routF,
		double	rinB,
		double	routB,
		double	zi,
		double	zo,
		double &	yh1,
		double &	bl1,
		double &	tl1,
		double &	yh2,
		double &	bl2,
		double &	tl2,
		double &	alp,
		double &	theta,
		double &	phi,
		double &	xpos,
		double &	ypos,
		double &	zcpv
	)

protected

Definition at line 890 of file DDHCalEndcapAlgo.cc.

References simBeamSpotPI::alpha, angle_units::operators::convertRadToDeg(), PVValHelper::dx, PVValHelper::dy, getNsectors(), hcalRecHitTable_cff::iphi, phi, alignCSCRings::r, mathSSE::sqrt(), funct::tan(), and theta().

Referenced by constructInsideModule(), and constructInsideModule0().

                                                     {
   //Given rin, rout compute parameters of the trapezoid and
   //position of the trapezoid for a standrd layer
   double alpha = (1._pi) / getNsectors();
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "Input " << iphi << " Front " << rinF << " " << routF << " " << zi << " Back " << rinB
                                << " " << routB << " " << zo << " Alpha " << convertRadToDeg(alpha);
 #endif
 
   yh1 = 0.5 * (routF - rinB);
   bl1 = 0.5 * rinB * tan(alpha);
   tl1 = 0.5 * routF * tan(alpha);
   yh2 = 0.5 * (routF - rinB);
   bl2 = 0.5 * rinB * tan(alpha);
   tl2 = 0.5 * routF * tan(alpha);
   double dx = 0.25 * (bl2 + tl2 - bl1 - tl1);
   double dy = 0.5 * (rinB + routF - rinB - routF);
   xpos = 0.25 * (rinB + routF + rinB + routF);
   ypos = 0.25 * (bl2 + tl2 + bl1 + tl1);
   zpos = 0.5 * (zi + zo);
   alp = atan(0.5 * tan(alpha));
   //  ypos-= getTolPos();
   if (iphi == 0) {
     ypos = -ypos;
   } else {
     alp = -alp;
     dx = -dx;
   }
   double r = sqrt(dx * dx + dy * dy);
   theta = atan(r / (zo - zi));
   phi = atan2(dy, dx);
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "Output Dimensions " << yh1 << " " << bl1 << " " << tl1 << " " << yh2 << " " << bl2
                                << " " << tl2 << " " << convertRadToDeg(alp) << " " << convertRadToDeg(theta) << " "
                                << convertRadToDeg(phi) << " Position " << xpos << " " << ypos << " " << zpos;
 #endif
 }

◆ parameterLayer0()

void DDHCalEndcapAlgo::parameterLayer0	(	int	mod,
		int	layer,
		int	iphi,
		double &	yh,
		double &	bl,
		double &	tl,
		double &	alp,
		double &	xpos,
		double &	ypos,
		double &	zcpv
	)

protected

Definition at line 833 of file DDHCalEndcapAlgo.cc.

References simBeamSpotPI::alpha, angle_units::operators::convertRadToDeg(), getAngBot(), getAngTop(), getLayerT(), getNsectors(), getSlope(), getTrim(), getZ1Beam(), getZmaxBlock(), getZminBlock(), hcalRecHitTable_cff::iphi, nano_mu_digi_cff::layer, mod(), rout, and funct::tan().

Referenced by constructInsideModule0().

                                                      {
   //Given module and layer number compute parameters of trapezoid
   //and positioning parameters
   double alpha = (1._pi) / getNsectors();
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "Input " << iphi << " " << layer << " " << iphi << " Alpha "
                                << convertRadToDeg(alpha);
 #endif
 
   double zi, zo;
   if (iphi == 0) {
     zi = getZminBlock(mod);
     zo = zi + getLayerT(layer);
   } else {
     zo = getZmaxBlock(mod);
     zi = zo - getLayerT(layer);
   }
   double rin, rout;
   if (mod == 0) {
     rin = zo * tan(getAngTop());
     rout = (zi - getZ1Beam()) * getSlope();
   } else {
     rin = zo * tan(getAngBot());
     rout = zi * tan(getAngTop());
   }
   yh = 0.5 * (rout - rin);
   bl = 0.5 * rin * tan(alpha);
   tl = 0.5 * rout * tan(alpha);
   xpos = 0.5 * (rin + rout);
   ypos = 0.5 * (bl + tl);
   zpos = 0.5 * (zi + zo);
   yh -= getTrim(mod, iphi);
   bl -= getTrim(mod, iphi);
   tl -= getTrim(mod, iphi);
   alp = atan(0.5 * tan(alpha));
   if (iphi == 0) {
     ypos = -ypos;
   } else {
     alp = -alp;
   }
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "Output Dimensions " << yh << " " << bl << " " << tl << " " << convertRadToDeg(alp)
                                << " Position " << xpos << " " << ypos << " " << zpos;
 #endif
 }

Member Data Documentation

◆ absMat

std::string DDHCalEndcapAlgo::absMat

private

Definition at line 188 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ angBot

double DDHCalEndcapAlgo::angBot

private

Definition at line 184 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ angGap

double DDHCalEndcapAlgo::angGap

private

Definition at line 185 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ angTop

double DDHCalEndcapAlgo::angTop

private

Definition at line 183 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ drEnd

double DDHCalEndcapAlgo::drEnd

private

Definition at line 181 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ dzShift

double DDHCalEndcapAlgo::dzShift

private

Definition at line 173 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ dzStep

double DDHCalEndcapAlgo::dzStep

private

Definition at line 172 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ eModule

std::vector<int> DDHCalEndcapAlgo::eModule

private

Definition at line 156 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ genMaterial

std::string DDHCalEndcapAlgo::genMaterial

private

Definition at line 152 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ heboxDepth

double DDHCalEndcapAlgo::heboxDepth

private

Definition at line 180 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ idName

std::string DDHCalEndcapAlgo::idName

private

Definition at line 222 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), constructInsideSector(), constructScintLayer(), and initialize().

◆ idNameSpace

std::string DDHCalEndcapAlgo::idNameSpace

private

Definition at line 223 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), constructInsideModule(), constructInsideModule0(), constructInsideSector(), constructScintLayer(), and initialize().

◆ idOffset

int DDHCalEndcapAlgo::idOffset

private

Definition at line 224 of file DDHCalEndcapAlgo.cc.

Referenced by constructInsideModule(), constructInsideModule0(), and initialize().

◆ layerN

std::vector<int> DDHCalEndcapAlgo::layerN

private

Definition at line 194 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ layerN0

std::vector<int> DDHCalEndcapAlgo::layerN0

private

Definition at line 195 of file DDHCalEndcapAlgo.cc.

Referenced by getLayer(), and initialize().

◆ layerN1

std::vector<int> DDHCalEndcapAlgo::layerN1

private

Definition at line 196 of file DDHCalEndcapAlgo.cc.

Referenced by getLayer(), and initialize().

◆ layerN2

std::vector<int> DDHCalEndcapAlgo::layerN2

private

Definition at line 197 of file DDHCalEndcapAlgo.cc.

Referenced by getLayer(), and initialize().

◆ layerN3

std::vector<int> DDHCalEndcapAlgo::layerN3

private

Definition at line 198 of file DDHCalEndcapAlgo.cc.

Referenced by getLayer(), and initialize().

◆ layerN4

std::vector<int> DDHCalEndcapAlgo::layerN4

private

Definition at line 199 of file DDHCalEndcapAlgo.cc.

Referenced by getLayer(), and initialize().

◆ layerN5

std::vector<int> DDHCalEndcapAlgo::layerN5

private

Definition at line 200 of file DDHCalEndcapAlgo.cc.

Referenced by getLayer(), and initialize().

◆ layerName

std::vector<std::string> DDHCalEndcapAlgo::layerName

private

Definition at line 214 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ layers

int DDHCalEndcapAlgo::layers

private

Definition at line 213 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ layerT

std::vector<double> DDHCalEndcapAlgo::layerT

private

Definition at line 216 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ layerType

std::vector<int> DDHCalEndcapAlgo::layerType

private

Definition at line 215 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ modMat

std::vector<std::string> DDHCalEndcapAlgo::modMat

private

Definition at line 191 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ modName

std::vector<std::string> DDHCalEndcapAlgo::modName

private

Definition at line 190 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ modType

std::vector<int> DDHCalEndcapAlgo::modType

private

Definition at line 192 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ modules

int DDHCalEndcapAlgo::modules

private

Definition at line 189 of file DDHCalEndcapAlgo.cc.

Referenced by initialize(), and postprocessor.PostProcessor::run().

◆ nEndcap

int DDHCalEndcapAlgo::nEndcap

private

Definition at line 155 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ nsectors

int DDHCalEndcapAlgo::nsectors

private

Definition at line 153 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ nsectortot

int DDHCalEndcapAlgo::nsectortot

private

Definition at line 154 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ phiName

std::vector<std::string> DDHCalEndcapAlgo::phiName

private

Definition at line 212 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ phiSections

int DDHCalEndcapAlgo::phiSections

private

Definition at line 211 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ plastMat

std::string DDHCalEndcapAlgo::plastMat

private

Definition at line 218 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ riDip

double DDHCalEndcapAlgo::riDip

private

Definition at line 178 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ riKink

double DDHCalEndcapAlgo::riKink

private

Definition at line 177 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ rinBlock1

std::vector<double> DDHCalEndcapAlgo::rinBlock1

private

Definition at line 206 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ rinBlock2

std::vector<double> DDHCalEndcapAlgo::rinBlock2

private

Definition at line 208 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ roDip

double DDHCalEndcapAlgo::roDip

private

Definition at line 179 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ rotation

std::string DDHCalEndcapAlgo::rotation

private

Definition at line 159 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and initialize().

◆ rotHalf

std::string DDHCalEndcapAlgo::rotHalf

private

Definition at line 157 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), and initialize().

◆ rotmat

std::string DDHCalEndcapAlgo::rotmat

private

Definition at line 220 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ rotns

std::string DDHCalEndcapAlgo::rotns

private

Definition at line 158 of file DDHCalEndcapAlgo.cc.

Referenced by constructGeneralVolume(), constructInsideModule(), constructInsideModule0(), and initialize().

◆ rout

double DDHCalEndcapAlgo::rout

private

Definition at line 176 of file DDHCalEndcapAlgo.cc.

Referenced by initialize(), and parameterLayer0().

◆ routBlock1

std::vector<double> DDHCalEndcapAlgo::routBlock1

private

Definition at line 207 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ routBlock2

std::vector<double> DDHCalEndcapAlgo::routBlock2

private

Definition at line 209 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ scintMat

std::string DDHCalEndcapAlgo::scintMat

private

Definition at line 219 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ scintT

std::vector<double> DDHCalEndcapAlgo::scintT

private

Definition at line 217 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ sectionModule

std::vector<int> DDHCalEndcapAlgo::sectionModule

private

Definition at line 193 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ slope

double DDHCalEndcapAlgo::slope

private

Definition at line 186 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ thick

std::vector<double> DDHCalEndcapAlgo::thick

private

Definition at line 201 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ tolAbs

double DDHCalEndcapAlgo::tolAbs

private

Definition at line 226 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ tolPos

double DDHCalEndcapAlgo::tolPos

private

Definition at line 226 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ trimLeft

std::vector<double> DDHCalEndcapAlgo::trimLeft

private

Definition at line 202 of file DDHCalEndcapAlgo.cc.

Referenced by getTrim(), and initialize().

◆ trimRight

std::vector<double> DDHCalEndcapAlgo::trimRight

private

Definition at line 203 of file DDHCalEndcapAlgo.cc.

Referenced by getTrim(), and initialize().

◆ z0Beam

double DDHCalEndcapAlgo::z0Beam

private

Definition at line 169 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ z1Beam

double DDHCalEndcapAlgo::z1Beam

private

Definition at line 170 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ zEnd

double DDHCalEndcapAlgo::zEnd

private

Definition at line 163 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ zFront

double DDHCalEndcapAlgo::zFront

private

Definition at line 162 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ ziBody

double DDHCalEndcapAlgo::ziBody

private

Definition at line 166 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ ziDip

double DDHCalEndcapAlgo::ziDip

private

Definition at line 171 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ ziKink

double DDHCalEndcapAlgo::ziKink

private

Definition at line 168 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ ziL0Body

double DDHCalEndcapAlgo::ziL0Body

private

Definition at line 167 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ ziL0Nose

double DDHCalEndcapAlgo::ziL0Nose

private

Definition at line 165 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ ziNose

double DDHCalEndcapAlgo::ziNose

private

Definition at line 164 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ zmaxBlock

std::vector<double> DDHCalEndcapAlgo::zmaxBlock

private

Definition at line 205 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ zminBlock

std::vector<double> DDHCalEndcapAlgo::zminBlock

private

Definition at line 204 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ zShift

double DDHCalEndcapAlgo::zShift

private

Definition at line 160 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

◆ zShiftHac2

double DDHCalEndcapAlgo::zShiftHac2

private

Definition at line 174 of file DDHCalEndcapAlgo.cc.

Referenced by initialize().

Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ DDHCalEndcapAlgo()

◆ ~DDHCalEndcapAlgo()

Member Function Documentation

◆ constructGeneralVolume()

◆ constructInsideModule()

◆ constructInsideModule0()

◆ constructInsideSector()

◆ constructScintLayer()

◆ equipModule()

◆ execute()

◆ getAbsMat()

◆ getAngBot()

◆ getAngGap()

◆ getAngTop()

◆ getDrEnd()

◆ getDzShift()

◆ getDzStep()

◆ getEndcaps()

◆ getGenMat()

◆ getHeboxDepth()

◆ getLayer()

◆ getLayerN()

◆ getLayerName()

◆ getLayers()

◆ getLayerT()

◆ getLayerType()

◆ getModMat()

◆ getModName()

◆ getModType()

◆ getModules()

◆ getNsectors()

◆ getNsectortot()

◆ getPhi()

◆ getPhiName()

◆ getPlastMat()

◆ getRinBlock1()

◆ getRinBlock2()

◆ getRinDip()

◆ getRinKink()

◆ getRotation()

◆ getRotMat()

◆ getRout()

◆ getRoutBlock1()

◆ getRoutBlock2()

◆ getRoutDip()

◆ getScintMat()

◆ getScintT()

◆ getSectionModule()

◆ getSlope()

◆ getThick()

◆ getTolAbs()

◆ getTolPos()

◆ getTrim()

◆ getZ0Beam()

◆ getZ1Beam()

◆ getZEnd()

◆ getZFront()

◆ getZiBody()

◆ getZiDip()

◆ getZiKink()

◆ getZiL0Body()

◆ getZiL0Nose()

◆ getZiNose()

◆ getZmaxBlock()

◆ getZminBlock()

◆ getZShift()

◆ getZShiftHac2()

◆ initialize()

◆ parameterLayer()

◆ parameterLayer0()

Member Data Documentation

◆ absMat

◆ angBot

◆ angGap

◆ angTop