#include <DDHCalEndcapAlgo.h>

Inheritance diagram for DDHCalEndcapAlgo:

Public Member Functions
	DDHCalEndcapAlgo ()

int	equipModule (unsigned int i) const

void	execute (DDCompactView &cpv) override

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

std::string	getGenMat () const

double	getHeboxDepth () const

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

int	getLayerN (unsigned int i) const

std::string	getLayerName (unsigned int i) const

int	getLayers () const

double	getLayerT (unsigned int i) const

int	getLayerType (unsigned int i) const

std::string	getModMat (unsigned int i) 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

std::string	getPhiName (unsigned int i) const

std::string	getPlastMat () const

double	getRinBlock1 (unsigned i) const

double	getRinBlock2 (unsigned i) const

double	getRinDip () const

double	getRinKink () const

std::string	getRotation () const

std::string	getRotMat () const

double	getRout () const

double	getRoutBlock1 (unsigned i) const

double	getRoutBlock2 (unsigned i) const

double	getRoutDip () 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 10 of file DDHCalEndcapAlgo.h.

Constructor & Destructor Documentation

DDHCalEndcapAlgo::DDHCalEndcapAlgo ( )

Definition at line 22 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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo info: Creating an instance";
 #endif
 }

DDHCalEndcapAlgo::~DDHCalEndcapAlgo ( )

override

Definition at line 33 of file DDHCalEndcapAlgo.cc.

33 {}

Member Function Documentation

void DDHCalEndcapAlgo::constructGeneralVolume ( DDCompactView & cpv )

protected

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

Definition at line 354 of file DDHCalEndcapAlgo.cc.

References alpha, constructInsideSector(), DDrot(), DDSplit(), equipModule(), plotBeamSpotDB::first, 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, dataset::name, DDBase< N, C >::name(), class-composition::parent, phi, pi, DDSolidFactory::polyhedra(), DDCompactView::position(), makeMuonMisalignmentScenario::rot, rotation, rotHalf, rotns, edm::second(), and funct::tan().

Referenced by execute(), and getTolAbs().

                                                                 {
   
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: 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::LogInfo("HCalGeom") << " First " << DDSplit(getRotation()).first
                << " Second " << DDSplit(getRotation()).second 
                << " Rotation " << rot;
 #endif
   
   DDTranslation r0(0,0,getZShift());
   double alpha = CLHEP::pi/getNsectors();
   double dphi  = getNsectortot()*CLHEP::twopi/getNsectors();
 
   //vertical walls are allowed in SolidPolyhedra
   double delz = 0;
 
   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()); 
 
   string name("Null");
   DDSolid solid;
   solid = DDSolidFactory::polyhedra(DDName(idName, idNameSpace),
                     getNsectortot(), -alpha, dphi, pgonZ, 
                     pgonRmin, pgonRmax);
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " 
                << DDName(idName, idNameSpace) 
                << " Polyhedra made of " << getGenMat() << " with " 
                << getNsectortot() << " sectors from " 
                << -alpha/CLHEP::deg << " to " 
                << (-alpha+dphi)/CLHEP::deg << " and with " 
                << pgonZ.size() << " sections";
   for (unsigned int i = 0; i <pgonZ.size(); i++) 
     edm::LogInfo("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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " 
                << 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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " 
                  << DDName(idName, idNameSpace) << " number 2 "
                  << "positioned in " << parentName  << " at " << r0
                  << " with " << rot;
 #endif
 
   }
 
   //Forward half
   name  = idName + "Front";
   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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: "
                << DDName(name,idNameSpace) << " Polyhedra made of "
                << getGenMat() << " with " << getNsectortot() 
                << " sectors from " << -alpha/CLHEP::deg << " to " 
                << (-alpha+dphi)/CLHEP::deg << " and with "
                << pgonZMod.size() << " sections ";
   for (unsigned int i = 0; i < pgonZMod.size(); i++) 
     edm::LogInfo("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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << 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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: "
                << DDName(name,idNameSpace)<< " Polyhedra made of "
                << getGenMat() <<" with 1 sector from " 
                << -alpha/CLHEP::deg << " to " << alpha/CLHEP::deg
                << " and with " << pgonZMod.size() << " sections";
   for (unsigned int i = 0; i < pgonZMod.size(); i++) 
     edm::LogInfo("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;
     double phideg = phi/CLHEP::deg;
     
     DDRotation rotation;
     string rotstr("NULL");
     if (phideg != 0) {
       rotstr = "R"; 
       if (phideg < 100) rotstr = "R0"; 
       rotstr = rotstr + std::to_string(phideg);
       rotation = DDRotation(DDName(rotstr, rotns)); 
       if (!rotation) {
 
 #ifdef EDM_ML_DEBUG
     edm::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: Creating a new "
                  << "rotation " << rotstr << "\t" << 90 << "," 
                  << phideg << ","  << 90 << "," << (phideg+90)
                  << ", 0, 0";
 #endif
     
     rotation = DDrot(DDName(rotstr, rotns), 90*CLHEP::deg, 
              phideg*CLHEP::deg, 90*CLHEP::deg, 
              (90+phideg)*CLHEP::deg, 0*CLHEP::deg,  0*CLHEP::deg);
       } //if !rotation
     } //if phideg!=0
   
    cpv.position(seclogic, genlogich, ii+1, DDTranslation(0.0, 0.0, 0.0), rotation);
 
 #ifdef EDM_ML_DEBUG
     edm::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << 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";
   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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: "
                << DDName(name,idNameSpace) << " Polyhedra made of " 
                << getAbsMat() << " with " << getNsectortot() 
                << " sectors from " << -alpha/CLHEP::deg << " to " 
                << (-alpha+dphi)/CLHEP::deg << " and with " 
                << pgonZBack.size()  << " sections";
   for (unsigned int i = 0; i < pgonZBack.size(); i++) 
     edm::LogInfo("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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << glog.name() 
                << " number 1 positioned in "  << genlogic.name() 
                << " at (0,0,0) with no rotation";
 #endif
 }

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

protected

Definition at line 904 of file DDHCalEndcapAlgo.cc.

References alpha, constructScintLayer(), DDBase< N, C >::ddname(), DDSplit(), plotBeamSpotDB::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, mod(), DDName::name(), dataset::name, DDBase< N, C >::name(), parameterLayer(), phi, pi, DDCompactView::position(), diffTwoXMLs::r1, diffTwoXMLs::r2, makeMuonMisalignmentScenario::rot, rotns, edm::second(), funct::tan(), theta(), and DDSolidFactory::trap().

Referenced by constructInsideSector(), and getTolAbs().

                                                                                                      {
   
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: \t\tInside module ..." 
                << mod;
 #endif
   
   //Pointers to the Rotation Matrices and to the Materials
   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 = CLHEP::pi/getNsectors();
   double  zi    = getZminBlock(mod);
 
   for (int i = 0; i < getLayerN(mod); i++) {
     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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: 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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << solid.name() 
                    << " Trap made of " << getGenMat() 
                    << " of dimensions " << 0.5*getThick(mod) << ", "
                    << theta/CLHEP::deg << ", " << phi/CLHEP::deg
                    << ", " << yh1 << ", " << bl1 << ", " << tl1 
                    << ", " << alp/CLHEP::deg << ", " << yh2 << ", "
                    << bl2 << ", " << tl2 << ", " << alp/CLHEP::deg;
 #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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << 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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << solid.name() 
                    << " Trap made of " << getPlastMat() 
                    << " of dimensions " << 0.5*getLayerT(layer) 
                    << ", 0, 0, " << yh << ", " << bl << ", " << tl
                    << ", " << alp/CLHEP::deg << ", " << yh << ", " 
                    << bl << ", " << tl << ", " << alp/CLHEP::deg;
 #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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << 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
 }

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

protected

Definition at line 782 of file DDHCalEndcapAlgo.cc.

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

Referenced by constructInsideSector(), and getTolAbs().

                                                                                                       {
   
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: \t\tInside module0 ..."
                << mod;
 #endif
   
   //Pointers to the Rotation Matrices and to the Materials
   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);
   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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << solid.name() 
                  << " Trap made of " << getPlastMat() 
                  << " of dimensions " << 0.5*getLayerT(layer) 
                  << ", 0, 0, " << yh << ", " << bl << ", " << tl 
                  << ", " << alp/CLHEP::deg << ", " << yh << ", " 
                  << bl << ", " << tl << ", " << alp/CLHEP::deg;
 #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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << 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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: 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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: 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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << solid.name() 
                << " Trap made of " << getAbsMat() 
                << " of dimensions " << 0.5*getThick(mod) << ", "
                << theta/CLHEP::deg << ", " << phi/CLHEP::deg << ", "
                << yh1 << ", " << bl1 << ", " << tl1 << ", " 
                << alp/CLHEP::deg << ", " << yh2 << ", "
                << bl2 << ", " << tl2 << ", " << alp/CLHEP::deg;
 #endif
   
   glog = DDLogicalPart(solid.ddname(), matabsorbr, solid);
 
   DDTranslation r2(xpos, ypos, zpos);
   cpv.position(glog, module, 1, r2, rot);
 
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << glog.name() 
                << " number 1 positioned in " << module.name() 
                << " at " << r2 << " with " << rot;
 #endif
 }

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

protected

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

Definition at line 594 of file DDHCalEndcapAlgo.cc.

References alpha, constructInsideModule(), constructInsideModule0(), DDSplit(), equipModule(), plotBeamSpotDB::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, gen::k, dataset::name, DDBase< N, C >::name(), pi, DDSolidFactory::polyhedra(), DDCompactView::position(), edm::second(), and funct::tan().

Referenced by constructGeneralVolume(), and getTolAbs().

                                                                                             {
   
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: Modules (" << getModules()
                << ") ...";
 #endif
 
   double alpha = CLHEP::pi/getNsectors();
 
   for (int i = 0; i < getModules(); i++) {
     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;
     
       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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " 
                    << DDName(name,idNameSpace) 
                    << " Polyhedra made of " << getModMat(i) 
                    << " with 1 sector from " << -alpha/CLHEP::deg 
                    << " to " << alpha/CLHEP::deg << " and with "
                    << nsec << " sections";
       for (unsigned int k=0; k<pgonZ.size(); k++)
     edm::LogInfo("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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << 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);
     }
   } 
 }

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 1026 of file DDHCalEndcapAlgo.cc.

References DDBase< N, C >::ddname(), DDSplit(), plotBeamSpotDB::first, getScintMat(), idName, idNameSpace, dataset::name, DDBase< N, C >::name(), DDCompactView::position(), edm::second(), and DDSolidFactory::trap().

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

                                                                                  {
 
   DDName matname(DDSplit(getScintMat()).first, DDSplit(getScintMat()).second);
   DDMaterial matter(matname);
   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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " 
                << DDName(name,idNameSpace) << " Trap made of " 
                << getScintMat() <<" of dimensions " << 0.5*dz 
                << ", 0, 0, " << yh << ", "  << bl << ", "  << tl 
                << ", " << alp/CLHEP::deg << ", " << yh << ", " 
                << bl << ", " << tl << ", " << alp/CLHEP::deg;
 #endif
   
   DDLogicalPart glog(solid.ddname(), matter, solid); 
 
   cpv.position(glog, detector, id, DDTranslation(0,0,0), DDRotation());
 
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HCalGeom") << "DDHCalEndcapAlgo test: " << glog.name() 
                << " number " << id << " positioned in " 
                << detector.name() << " at (0,0,0) with no rotation";
 #endif
 }

int DDHCalEndcapAlgo::equipModule ( unsigned int i ) const

inline

Definition at line 22 of file DDHCalEndcapAlgo.h.

References eModule, and mps_fire::i.

Referenced by constructGeneralVolume(), and constructInsideSector().

22 {return eModule[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::eModule

std::vector< int > eModule

Definition: DDHCalEndcapAlgo.h:112

void DDHCalEndcapAlgo::execute ( DDCompactView & cpv )

override

Definition at line 339 of file DDHCalEndcapAlgo.cc.

References constructGeneralVolume().

Referenced by getTolAbs().

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

std::string DDHCalEndcapAlgo::getAbsMat ( ) const

inline

Definition at line 50 of file DDHCalEndcapAlgo.h.

References absMat.

Referenced by constructGeneralVolume(), and constructInsideModule0().

50 {return absMat;}

DDHCalEndcapAlgo::absMat

std::string absMat

Definition: DDHCalEndcapAlgo.h:144

double DDHCalEndcapAlgo::getAngBot ( ) const

inline

Definition at line 46 of file DDHCalEndcapAlgo.h.

References angBot.

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

46 {return angBot;}

DDHCalEndcapAlgo::angBot

double angBot

Definition: DDHCalEndcapAlgo.h:140

double DDHCalEndcapAlgo::getAngGap ( ) const

inline

Definition at line 47 of file DDHCalEndcapAlgo.h.

References angGap.

47 {return angGap;}

DDHCalEndcapAlgo::angGap

double angGap

Definition: DDHCalEndcapAlgo.h:141

double DDHCalEndcapAlgo::getAngTop ( ) const

inline

Definition at line 45 of file DDHCalEndcapAlgo.h.

References angTop.

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

45 {return angTop;}

DDHCalEndcapAlgo::angTop

double angTop

Definition: DDHCalEndcapAlgo.h:139

double DDHCalEndcapAlgo::getDrEnd ( ) const

inline

Definition at line 44 of file DDHCalEndcapAlgo.h.

References drEnd.

Referenced by constructGeneralVolume().

44 {return drEnd;}

DDHCalEndcapAlgo::drEnd

double drEnd

Definition: DDHCalEndcapAlgo.h:137

double DDHCalEndcapAlgo::getDzShift ( ) const

inline

Definition at line 36 of file DDHCalEndcapAlgo.h.

References dzShift.

Referenced by constructGeneralVolume().

36 {return dzShift;}

DDHCalEndcapAlgo::dzShift

double dzShift

Definition: DDHCalEndcapAlgo.h:129

double DDHCalEndcapAlgo::getDzStep ( ) const

inline

Definition at line 35 of file DDHCalEndcapAlgo.h.

References dzStep.

Referenced by constructInsideModule(), and constructInsideModule0().

35 {return dzStep;}

DDHCalEndcapAlgo::dzStep

double dzStep

Definition: DDHCalEndcapAlgo.h:128

int DDHCalEndcapAlgo::getEndcaps ( ) const

inline

Definition at line 21 of file DDHCalEndcapAlgo.h.

References nEndcap.

Referenced by constructGeneralVolume().

21 {return nEndcap;}

DDHCalEndcapAlgo::nEndcap

int nEndcap

Definition: DDHCalEndcapAlgo.h:111

std::string DDHCalEndcapAlgo::getGenMat ( ) const

inline

Definition at line 17 of file DDHCalEndcapAlgo.h.

References genMaterial.

Referenced by constructGeneralVolume(), and constructInsideModule().

17 {return genMaterial;}

DDHCalEndcapAlgo::genMaterial

std::string genMaterial

Definition: DDHCalEndcapAlgo.h:108

double DDHCalEndcapAlgo::getHeboxDepth ( ) const

inline

Definition at line 43 of file DDHCalEndcapAlgo.h.

References heboxDepth.

43 {return heboxDepth;}

DDHCalEndcapAlgo::heboxDepth

double heboxDepth

Definition: DDHCalEndcapAlgo.h:136

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

Definition at line 35 of file DDHCalEndcapAlgo.cc.

References layerN0, layerN1, layerN2, layerN3, layerN4, and layerN5.

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

                                                                    {
 
   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;
   }
 }

int DDHCalEndcapAlgo::getLayerN ( unsigned int i ) const

inline

Definition at line 56 of file DDHCalEndcapAlgo.h.

References getLayer(), mps_fire::i, and layerN.

Referenced by constructInsideModule().

56 {return layerN[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::layerN

std::vector< int > layerN

Definition: DDHCalEndcapAlgo.h:150

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

inline

Definition at line 70 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and layerName.

Referenced by constructInsideModule(), and constructInsideModule0().

70 {return layerName[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::layerName

std::vector< std::string > layerName

Definition: DDHCalEndcapAlgo.h:170

int DDHCalEndcapAlgo::getLayers ( ) const

inline

Definition at line 69 of file DDHCalEndcapAlgo.h.

References layers.

69 {return layers;}

DDHCalEndcapAlgo::layers

int layers

Definition: DDHCalEndcapAlgo.h:169

double DDHCalEndcapAlgo::getLayerT ( unsigned int i ) const

inline

Definition at line 72 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and layerT.

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

72 {return layerT[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::layerT

std::vector< double > layerT

Definition: DDHCalEndcapAlgo.h:172

int DDHCalEndcapAlgo::getLayerType ( unsigned int i ) const

inline

Definition at line 71 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and layerType.

Referenced by constructInsideModule(), and constructInsideModule0().

71 {return layerType[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::layerType

std::vector< int > layerType

Definition: DDHCalEndcapAlgo.h:171

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

inline

Definition at line 53 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and modMat.

Referenced by constructInsideSector().

53 {return modMat[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::modMat

std::vector< std::string > modMat

Definition: DDHCalEndcapAlgo.h:147

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

inline

Definition at line 52 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and modName.

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

52 {return modName[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::modName

std::vector< std::string > modName

Definition: DDHCalEndcapAlgo.h:146

int DDHCalEndcapAlgo::getModType ( unsigned int i ) const

inline

Definition at line 54 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and modType.

Referenced by constructInsideSector().

54 {return modType[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::modType

std::vector< int > modType

Definition: DDHCalEndcapAlgo.h:148

int DDHCalEndcapAlgo::getModules ( ) const

inline

Definition at line 51 of file DDHCalEndcapAlgo.h.

References modules.

Referenced by constructInsideSector().

51 {return modules;}

DDHCalEndcapAlgo::modules

int modules

Definition: DDHCalEndcapAlgo.h:145

int DDHCalEndcapAlgo::getNsectors ( ) const

inline

Definition at line 19 of file DDHCalEndcapAlgo.h.

References nsectors.

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

19 {return nsectors;}

DDHCalEndcapAlgo::nsectors

int nsectors

Definition: DDHCalEndcapAlgo.h:109

int DDHCalEndcapAlgo::getNsectortot ( ) const

inline

Definition at line 20 of file DDHCalEndcapAlgo.h.

References nsectortot.

Referenced by constructGeneralVolume().

20 {return nsectortot;}

DDHCalEndcapAlgo::nsectortot

int nsectortot

Definition: DDHCalEndcapAlgo.h:110

int DDHCalEndcapAlgo::getPhi ( ) const

inline

Definition at line 67 of file DDHCalEndcapAlgo.h.

References phiSections.

Referenced by constructInsideModule(), and constructInsideModule0().

67 {return phiSections;}

DDHCalEndcapAlgo::phiSections

int phiSections

Definition: DDHCalEndcapAlgo.h:167

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

inline

Definition at line 68 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and phiName.

Referenced by constructInsideModule(), and constructInsideModule0().

68 {return phiName[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::phiName

std::vector< std::string > phiName

Definition: DDHCalEndcapAlgo.h:168

std::string DDHCalEndcapAlgo::getPlastMat ( ) const

inline

Definition at line 74 of file DDHCalEndcapAlgo.h.

References plastMat.

Referenced by constructInsideModule(), and constructInsideModule0().

74 {return plastMat;}

DDHCalEndcapAlgo::plastMat

std::string plastMat

Definition: DDHCalEndcapAlgo.h:174

double DDHCalEndcapAlgo::getRinBlock1 ( unsigned i ) const

inline

Definition at line 62 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and rinBlock1.

Referenced by constructInsideSector().

62 {return rinBlock1[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::rinBlock1

std::vector< double > rinBlock1

Definition: DDHCalEndcapAlgo.h:162

double DDHCalEndcapAlgo::getRinBlock2 ( unsigned i ) const

inline

Definition at line 63 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and rinBlock2.

Referenced by constructInsideSector().

63 {return rinBlock2[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::rinBlock2

std::vector< double > rinBlock2

Definition: DDHCalEndcapAlgo.h:164

double DDHCalEndcapAlgo::getRinDip ( ) const

inline

Definition at line 41 of file DDHCalEndcapAlgo.h.

References riDip.

Referenced by constructGeneralVolume(), and constructInsideSector().

41 {return riDip;}

DDHCalEndcapAlgo::riDip

double riDip

Definition: DDHCalEndcapAlgo.h:134

double DDHCalEndcapAlgo::getRinKink ( ) const

inline

Definition at line 40 of file DDHCalEndcapAlgo.h.

References riKink.

Referenced by constructGeneralVolume(), and constructInsideSector().

40 {return riKink;}

DDHCalEndcapAlgo::riKink

double riKink

Definition: DDHCalEndcapAlgo.h:133

std::string DDHCalEndcapAlgo::getRotation ( ) const

inline

Definition at line 18 of file DDHCalEndcapAlgo.h.

References rotation.

Referenced by constructGeneralVolume().

18 {return rotation;}

DDHCalEndcapAlgo::rotation

std::string rotation

Definition: DDHCalEndcapAlgo.h:115

std::string DDHCalEndcapAlgo::getRotMat ( ) const

inline

Definition at line 76 of file DDHCalEndcapAlgo.h.

References rotmat.

Referenced by constructInsideModule(), and constructInsideModule0().

76 {return rotmat;}

DDHCalEndcapAlgo::rotmat

std::string rotmat

Definition: DDHCalEndcapAlgo.h:176

double DDHCalEndcapAlgo::getRout ( ) const

inline

Definition at line 39 of file DDHCalEndcapAlgo.h.

References rout.

Referenced by constructGeneralVolume(), and constructInsideSector().

39 {return rout;}

DDHCalEndcapAlgo::rout

double rout

Definition: DDHCalEndcapAlgo.h:132

double DDHCalEndcapAlgo::getRoutBlock1 ( unsigned i ) const

inline

Definition at line 64 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and routBlock1.

Referenced by constructInsideSector().

64 {return routBlock1[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::routBlock1

std::vector< double > routBlock1

Definition: DDHCalEndcapAlgo.h:163

double DDHCalEndcapAlgo::getRoutBlock2 ( unsigned i ) const

inline

Definition at line 65 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and routBlock2.

Referenced by constructInsideModule(), and constructInsideSector().

65 {return routBlock2[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::routBlock2

std::vector< double > routBlock2

Definition: DDHCalEndcapAlgo.h:165

double DDHCalEndcapAlgo::getRoutDip ( ) const

inline

Definition at line 42 of file DDHCalEndcapAlgo.h.

References roDip.

Referenced by constructGeneralVolume(), and constructInsideSector().

42 {return roDip;}

DDHCalEndcapAlgo::roDip

double roDip

Definition: DDHCalEndcapAlgo.h:135

std::string DDHCalEndcapAlgo::getScintMat ( ) const

inline

Definition at line 75 of file DDHCalEndcapAlgo.h.

References scintMat.

Referenced by constructScintLayer().

75 {return scintMat;}

DDHCalEndcapAlgo::scintMat

std::string scintMat

Definition: DDHCalEndcapAlgo.h:175

double DDHCalEndcapAlgo::getScintT ( unsigned int i ) const

inline

Definition at line 73 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and scintT.

Referenced by constructInsideModule(), and constructInsideModule0().

73 {return scintT[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::scintT

std::vector< double > scintT

Definition: DDHCalEndcapAlgo.h:173

int DDHCalEndcapAlgo::getSectionModule ( unsigned i ) const

inline

Definition at line 55 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and sectionModule.

Referenced by constructInsideSector().

55 {return sectionModule[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::sectionModule

std::vector< int > sectionModule

Definition: DDHCalEndcapAlgo.h:149

double DDHCalEndcapAlgo::getSlope ( ) const

inline

Definition at line 48 of file DDHCalEndcapAlgo.h.

References slope.

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

48 {return slope;}

DDHCalEndcapAlgo::slope

double slope

Definition: DDHCalEndcapAlgo.h:142

double DDHCalEndcapAlgo::getThick ( unsigned int i ) const

inline

Definition at line 58 of file DDHCalEndcapAlgo.h.

References getTrim(), mps_fire::i, and thick.

Referenced by constructInsideModule(), and constructInsideModule0().

58 {return thick[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::thick

std::vector< double > thick

Definition: DDHCalEndcapAlgo.h:157

double DDHCalEndcapAlgo::getTolAbs ( ) const

inline

Definition at line 78 of file DDHCalEndcapAlgo.h.

References constructGeneralVolume(), constructInsideModule(), constructInsideModule0(), constructInsideSector(), constructScintLayer(), execute(), initialize(), mod(), dataset::name, parameterLayer(), parameterLayer0(), phi, AlCaHLTBitMon_QueryRunRegistry::string, theta(), and tolAbs.

Referenced by constructInsideModule0().

78 {return tolAbs;}

DDHCalEndcapAlgo::tolAbs

double tolAbs

Definition: DDHCalEndcapAlgo.h:182

double DDHCalEndcapAlgo::getTolPos ( ) const

inline

Definition at line 77 of file DDHCalEndcapAlgo.h.

References tolPos.

77 {return tolPos;}

DDHCalEndcapAlgo::tolPos

double tolPos

Definition: DDHCalEndcapAlgo.h:182

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

Definition at line 67 of file DDHCalEndcapAlgo.cc.

References mps_fire::i, trimLeft, and trimRight.

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

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

double DDHCalEndcapAlgo::getZ0Beam ( ) const

inline

Definition at line 32 of file DDHCalEndcapAlgo.h.

References z0Beam.

32 {return z0Beam;}

DDHCalEndcapAlgo::z0Beam

double z0Beam

Definition: DDHCalEndcapAlgo.h:125

double DDHCalEndcapAlgo::getZ1Beam ( ) const

inline

Definition at line 33 of file DDHCalEndcapAlgo.h.

References z1Beam.

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

33 {return z1Beam;}

DDHCalEndcapAlgo::z1Beam

double z1Beam

Definition: DDHCalEndcapAlgo.h:126

double DDHCalEndcapAlgo::getZEnd ( ) const

inline

Definition at line 26 of file DDHCalEndcapAlgo.h.

References zEnd.

Referenced by constructGeneralVolume().

26 {return zEnd;}

DDHCalEndcapAlgo::zEnd

double zEnd

Definition: DDHCalEndcapAlgo.h:119

double DDHCalEndcapAlgo::getZFront ( ) const

inline

Definition at line 25 of file DDHCalEndcapAlgo.h.

References zFront.

Referenced by constructGeneralVolume().

25 {return zFront;}

DDHCalEndcapAlgo::zFront

double zFront

Definition: DDHCalEndcapAlgo.h:118

double DDHCalEndcapAlgo::getZiBody ( ) const

inline

Definition at line 29 of file DDHCalEndcapAlgo.h.

References ziBody.

Referenced by constructGeneralVolume().

29 {return ziBody;}

DDHCalEndcapAlgo::ziBody

double ziBody

Definition: DDHCalEndcapAlgo.h:122

double DDHCalEndcapAlgo::getZiDip ( ) const

inline

Definition at line 34 of file DDHCalEndcapAlgo.h.

References ziDip.

Referenced by constructGeneralVolume(), and constructInsideSector().

34 {return ziDip;}

DDHCalEndcapAlgo::ziDip

double ziDip

Definition: DDHCalEndcapAlgo.h:127

double DDHCalEndcapAlgo::getZiKink ( ) const

inline

Definition at line 31 of file DDHCalEndcapAlgo.h.

References ziKink.

Referenced by constructGeneralVolume(), and constructInsideSector().

31 {return ziKink;}

DDHCalEndcapAlgo::ziKink

double ziKink

Definition: DDHCalEndcapAlgo.h:124

double DDHCalEndcapAlgo::getZiL0Body ( ) const

inline

Definition at line 30 of file DDHCalEndcapAlgo.h.

References ziL0Body.

Referenced by constructGeneralVolume().

30 {return ziL0Body;}

DDHCalEndcapAlgo::ziL0Body

double ziL0Body

Definition: DDHCalEndcapAlgo.h:123

double DDHCalEndcapAlgo::getZiL0Nose ( ) const

inline

Definition at line 28 of file DDHCalEndcapAlgo.h.

References ziL0Nose.

28 {return ziL0Nose;}

DDHCalEndcapAlgo::ziL0Nose

double ziL0Nose

Definition: DDHCalEndcapAlgo.h:121

double DDHCalEndcapAlgo::getZiNose ( ) const

inline

Definition at line 27 of file DDHCalEndcapAlgo.h.

References ziNose.

27 {return ziNose;}

DDHCalEndcapAlgo::ziNose

double ziNose

Definition: DDHCalEndcapAlgo.h:120

double DDHCalEndcapAlgo::getZmaxBlock ( unsigned i ) const

inline

Definition at line 61 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and zmaxBlock.

Referenced by constructInsideSector(), and parameterLayer0().

61 {return zmaxBlock[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::zmaxBlock

std::vector< double > zmaxBlock

Definition: DDHCalEndcapAlgo.h:161

double DDHCalEndcapAlgo::getZminBlock ( unsigned i ) const

inline

Definition at line 60 of file DDHCalEndcapAlgo.h.

References mps_fire::i, and zminBlock.

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

60 {return zminBlock[i];}

mps_fire.i

i

Definition: mps_fire.py:269

DDHCalEndcapAlgo::zminBlock

std::vector< double > zminBlock

Definition: DDHCalEndcapAlgo.h:160

double DDHCalEndcapAlgo::getZShift ( ) const

inline

Definition at line 23 of file DDHCalEndcapAlgo.h.

References zShift.

Referenced by constructGeneralVolume().

23 {return zShift;}

DDHCalEndcapAlgo::zShift

double zShift

Definition: DDHCalEndcapAlgo.h:116

double DDHCalEndcapAlgo::getZShiftHac2 ( ) const

inline

Definition at line 37 of file DDHCalEndcapAlgo.h.

References zShiftHac2.

Referenced by constructGeneralVolume(), and constructInsideSector().

37 {return zShiftHac2;}

DDHCalEndcapAlgo::zShiftHac2

double zShiftHac2

Definition: DDHCalEndcapAlgo.h:130

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

override

Definition at line 75 of file DDHCalEndcapAlgo.cc.

References absMat, angBot, angGap, angTop, dbl_to_int(), drEnd, dzShift, dzStep, eModule, JetChargeProducer_cfi::exp, genMaterial, heboxDepth, mps_fire::i, idName, idNameSpace, idOffset, createfilelist::int, layerN, layerN0, layerN1, layerN2, layerN3, layerN4, layerN5, layerName, layers, layerT, layerType, modMat, modName, modType, python.rootplot.argparse::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.

Referenced by getTolAbs().

                                                           {
 
   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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo debug: 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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo debug: angTop " 
                << angTop/CLHEP::deg  <<"\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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo debug: 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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo debug: " << 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::LogInfo("HCalGeom") << "\t " << layerN0[j] << "/" << layerN0[j+1];
       }
     } else if (i == 1) {
       for (int j = 0; j < layerN[i]; j++) {
     edm::LogInfo("HCalGeom") << "\t " << layerN1[j] << "/" << layerN1[j+1];
       }
     } else if (i == 2) {
       for (int j = 0; j < layerN[i]; j++) {
     edm::LogInfo("HCalGeom") << "\t " << layerN2[j];
       }
     } else if (i == 3) {
       for (int j = 0; j < layerN[i]; j++) {
     edm::LogInfo("HCalGeom") << "\t " << layerN3[j];
       }
     } else if (i == 4) {
       for (int j = 0; j < layerN[i]; j++) {
     edm::LogInfo("HCalGeom") << "\t " << layerN4[j];
       }
     } else if (i == 5) {
       for (int j = 0; j < layerN[i]; j++) {
     edm::LogInfo("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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo debug: Phi Sections " 
                << phiSections;
   for (int i = 0; i < phiSections; i++) 
     edm::LogInfo("HCalGeom") << "\tName[" << i << "] : " << phiName[i];
   edm::LogInfo("HCalGeom") << "\tPlastic: " << plastMat << "\tScintillator: "
                << scintMat << "\tRotation matrix " << rotns << ":" 
                << rotmat << "\n\tNumber of layers " << layers;
   for (int i = 0; i < layers; i++) {
     edm::LogInfo("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::LogInfo("HCalGeom") << "DDHCalEndcapAlgo debug: 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
   DDName parentName = parent().name(); 
   edm::LogInfo("HCalGeom") << "DDHCalEndcapAlgo debug: Parent " << parentName 
                << " idName " << idName << " NameSpace " 
                << idNameSpace << " Offset " << idOffset;
 #endif
 
   tolPos      = nArgs["TolPos"];
   tolAbs      = nArgs["TolAbs"];
 
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HCalGeom") <<"DDHCalEndcapAlgo debug: Tolerances - Positioning "
                << tolPos << " Absorber " << tolAbs;
 #endif
 }

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 730 of file DDHCalEndcapAlgo.cc.

References alpha, PVValHelper::dx, PVValHelper::dy, getNsectors(), pi, alignCSCRings::r, mathSSE::sqrt(), and funct::tan().

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

                                     {
 
   //Given rin, rout compute parameters of the trapezoid and 
   //position of the trapezoid for a standrd layer
   double alpha = CLHEP::pi/getNsectors();
 
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HCalGeom") << "Input " << iphi << " Front " << rinF << " " 
                << routF << " " << zi << " Back " << rinB << " " 
                << routB << " " << zo << " Alpha " 
                << alpha/CLHEP::deg;
 #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::LogInfo("HCalGeom") << "Output Dimensions " << yh1 << " " << bl1 << " "
                << tl1 << " " << yh2 << " " << bl2 << " " << tl2
                << " " << alp/CLHEP::deg << " " << theta/CLHEP::deg
                << " " << phi/CLHEP::deg << " Position " << xpos 
                << " " << ypos << " " << zpos;
 #endif
 }

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 677 of file DDHCalEndcapAlgo.cc.

References alpha, getAngBot(), getAngTop(), getLayerT(), getNsectors(), getSlope(), getTrim(), getZ1Beam(), getZmaxBlock(), getZminBlock(), pi, rout, and funct::tan().

Referenced by constructInsideModule0(), and getTolAbs().

                                                      {
 
   //Given module and layer number compute parameters of trapezoid
   //and positioning parameters
   double alpha = CLHEP::pi/getNsectors();
 
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HCalGeom") << "Input " << iphi << " " << layer << " " << iphi
                << " Alpha " << alpha/CLHEP::deg;
 #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::LogInfo("HCalGeom") << "Output Dimensions " << yh << " " << bl << " "
                << tl << " " << alp/CLHEP::deg << " Position " 
                << xpos << " " << ypos << " " << zpos;
 #endif
 }

Member Data Documentation

std::string DDHCalEndcapAlgo::absMat

private

Definition at line 144 of file DDHCalEndcapAlgo.h.

Referenced by getAbsMat(), and initialize().

double DDHCalEndcapAlgo::angBot

private

Definition at line 140 of file DDHCalEndcapAlgo.h.

Referenced by getAngBot(), and initialize().

double DDHCalEndcapAlgo::angGap

private

Definition at line 141 of file DDHCalEndcapAlgo.h.

Referenced by getAngGap(), and initialize().

double DDHCalEndcapAlgo::angTop

private

Definition at line 139 of file DDHCalEndcapAlgo.h.

Referenced by getAngTop(), and initialize().

double DDHCalEndcapAlgo::drEnd

private

Definition at line 137 of file DDHCalEndcapAlgo.h.

Referenced by getDrEnd(), and initialize().

double DDHCalEndcapAlgo::dzShift

private

Definition at line 129 of file DDHCalEndcapAlgo.h.

Referenced by getDzShift(), and initialize().

double DDHCalEndcapAlgo::dzStep

private

Definition at line 128 of file DDHCalEndcapAlgo.h.

Referenced by getDzStep(), and initialize().

std::vector<int> DDHCalEndcapAlgo::eModule

private

Definition at line 112 of file DDHCalEndcapAlgo.h.

Referenced by equipModule(), and initialize().

std::string DDHCalEndcapAlgo::genMaterial

private

Definition at line 108 of file DDHCalEndcapAlgo.h.

Referenced by getGenMat(), and initialize().

double DDHCalEndcapAlgo::heboxDepth

private

Definition at line 136 of file DDHCalEndcapAlgo.h.

Referenced by getHeboxDepth(), and initialize().

std::string DDHCalEndcapAlgo::idName

private

Definition at line 178 of file DDHCalEndcapAlgo.h.

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

std::string DDHCalEndcapAlgo::idNameSpace

private

Definition at line 179 of file DDHCalEndcapAlgo.h.

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

int DDHCalEndcapAlgo::idOffset

private

Definition at line 180 of file DDHCalEndcapAlgo.h.

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

std::vector<int> DDHCalEndcapAlgo::layerN

private

Definition at line 150 of file DDHCalEndcapAlgo.h.

Referenced by getLayerN(), and initialize().

std::vector<int> DDHCalEndcapAlgo::layerN0

private

Definition at line 151 of file DDHCalEndcapAlgo.h.

Referenced by getLayer(), and initialize().

std::vector<int> DDHCalEndcapAlgo::layerN1

private

Definition at line 152 of file DDHCalEndcapAlgo.h.

Referenced by getLayer(), and initialize().

std::vector<int> DDHCalEndcapAlgo::layerN2

private

Definition at line 153 of file DDHCalEndcapAlgo.h.

Referenced by getLayer(), and initialize().

std::vector<int> DDHCalEndcapAlgo::layerN3

private

Definition at line 154 of file DDHCalEndcapAlgo.h.

Referenced by getLayer(), and initialize().

std::vector<int> DDHCalEndcapAlgo::layerN4

private

Definition at line 155 of file DDHCalEndcapAlgo.h.

Referenced by getLayer(), and initialize().

std::vector<int> DDHCalEndcapAlgo::layerN5

private

Definition at line 156 of file DDHCalEndcapAlgo.h.

Referenced by getLayer(), and initialize().

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

private

Definition at line 170 of file DDHCalEndcapAlgo.h.

Referenced by getLayerName(), and initialize().

int DDHCalEndcapAlgo::layers

private

Definition at line 169 of file DDHCalEndcapAlgo.h.

Referenced by getLayers(), and initialize().

std::vector<double> DDHCalEndcapAlgo::layerT

private

Definition at line 172 of file DDHCalEndcapAlgo.h.

Referenced by getLayerT(), and initialize().

std::vector<int> DDHCalEndcapAlgo::layerType

private

Definition at line 171 of file DDHCalEndcapAlgo.h.

Referenced by getLayerType(), and initialize().

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

private

Definition at line 147 of file DDHCalEndcapAlgo.h.

Referenced by getModMat(), and initialize().

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

private

Definition at line 146 of file DDHCalEndcapAlgo.h.

Referenced by getModName(), and initialize().

std::vector<int> DDHCalEndcapAlgo::modType

private

Definition at line 148 of file DDHCalEndcapAlgo.h.

Referenced by getModType(), and initialize().

int DDHCalEndcapAlgo::modules

private

Definition at line 145 of file DDHCalEndcapAlgo.h.

Referenced by getModules(), and initialize().

int DDHCalEndcapAlgo::nEndcap

private

Definition at line 111 of file DDHCalEndcapAlgo.h.

Referenced by getEndcaps(), and initialize().

int DDHCalEndcapAlgo::nsectors

private

Definition at line 109 of file DDHCalEndcapAlgo.h.

Referenced by getNsectors(), and initialize().

int DDHCalEndcapAlgo::nsectortot

private

Definition at line 110 of file DDHCalEndcapAlgo.h.

Referenced by getNsectortot(), and initialize().

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

private

Definition at line 168 of file DDHCalEndcapAlgo.h.

Referenced by getPhiName(), and initialize().

int DDHCalEndcapAlgo::phiSections

private

Definition at line 167 of file DDHCalEndcapAlgo.h.

Referenced by getPhi(), and initialize().

std::string DDHCalEndcapAlgo::plastMat

private

Definition at line 174 of file DDHCalEndcapAlgo.h.

Referenced by getPlastMat(), and initialize().

double DDHCalEndcapAlgo::riDip

private

Definition at line 134 of file DDHCalEndcapAlgo.h.

Referenced by getRinDip(), and initialize().

double DDHCalEndcapAlgo::riKink

private

Definition at line 133 of file DDHCalEndcapAlgo.h.

Referenced by getRinKink(), and initialize().

std::vector<double> DDHCalEndcapAlgo::rinBlock1

private

Definition at line 162 of file DDHCalEndcapAlgo.h.

Referenced by getRinBlock1(), and initialize().

std::vector<double> DDHCalEndcapAlgo::rinBlock2

private

Definition at line 164 of file DDHCalEndcapAlgo.h.

Referenced by getRinBlock2(), and initialize().

double DDHCalEndcapAlgo::roDip

private

Definition at line 135 of file DDHCalEndcapAlgo.h.

Referenced by getRoutDip(), and initialize().

std::string DDHCalEndcapAlgo::rotation

private

Definition at line 115 of file DDHCalEndcapAlgo.h.

Referenced by constructGeneralVolume(), getRotation(), and initialize().

std::string DDHCalEndcapAlgo::rotHalf

private

Definition at line 113 of file DDHCalEndcapAlgo.h.

Referenced by constructGeneralVolume(), and initialize().

std::string DDHCalEndcapAlgo::rotmat

private

Definition at line 176 of file DDHCalEndcapAlgo.h.

Referenced by getRotMat(), and initialize().

std::string DDHCalEndcapAlgo::rotns

private

Definition at line 114 of file DDHCalEndcapAlgo.h.

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

double DDHCalEndcapAlgo::rout

private

Definition at line 132 of file DDHCalEndcapAlgo.h.

Referenced by getRout(), initialize(), and parameterLayer0().

std::vector<double> DDHCalEndcapAlgo::routBlock1

private

Definition at line 163 of file DDHCalEndcapAlgo.h.

Referenced by getRoutBlock1(), and initialize().

std::vector<double> DDHCalEndcapAlgo::routBlock2

private

Definition at line 165 of file DDHCalEndcapAlgo.h.

Referenced by getRoutBlock2(), and initialize().

std::string DDHCalEndcapAlgo::scintMat

private

Definition at line 175 of file DDHCalEndcapAlgo.h.

Referenced by getScintMat(), and initialize().

std::vector<double> DDHCalEndcapAlgo::scintT

private

Definition at line 173 of file DDHCalEndcapAlgo.h.

Referenced by getScintT(), and initialize().

std::vector<int> DDHCalEndcapAlgo::sectionModule

private

Definition at line 149 of file DDHCalEndcapAlgo.h.

Referenced by getSectionModule(), and initialize().

double DDHCalEndcapAlgo::slope

private

Definition at line 142 of file DDHCalEndcapAlgo.h.

Referenced by Vispa.Views.LineDecayView.DecayLine::containsPoint(), getSlope(), initialize(), and Vispa.Views.LineDecayView.DecayLine::transform().

std::vector<double> DDHCalEndcapAlgo::thick

private

Definition at line 157 of file DDHCalEndcapAlgo.h.

Referenced by getThick(), and initialize().

double DDHCalEndcapAlgo::tolAbs

private

Definition at line 182 of file DDHCalEndcapAlgo.h.

Referenced by getTolAbs(), and initialize().

double DDHCalEndcapAlgo::tolPos

private

Definition at line 182 of file DDHCalEndcapAlgo.h.

Referenced by getTolPos(), and initialize().

std::vector<double> DDHCalEndcapAlgo::trimLeft

private

Definition at line 158 of file DDHCalEndcapAlgo.h.

Referenced by getTrim(), and initialize().

std::vector<double> DDHCalEndcapAlgo::trimRight

private

Definition at line 159 of file DDHCalEndcapAlgo.h.

Referenced by getTrim(), and initialize().

double DDHCalEndcapAlgo::z0Beam

private

Definition at line 125 of file DDHCalEndcapAlgo.h.

Referenced by getZ0Beam(), and initialize().

double DDHCalEndcapAlgo::z1Beam

private

Definition at line 126 of file DDHCalEndcapAlgo.h.

Referenced by getZ1Beam(), and initialize().

double DDHCalEndcapAlgo::zEnd

private

Definition at line 119 of file DDHCalEndcapAlgo.h.

Referenced by getZEnd(), and initialize().

double DDHCalEndcapAlgo::zFront

private

Definition at line 118 of file DDHCalEndcapAlgo.h.

Referenced by getZFront(), and initialize().

double DDHCalEndcapAlgo::ziBody

private

Definition at line 122 of file DDHCalEndcapAlgo.h.

Referenced by getZiBody(), and initialize().

double DDHCalEndcapAlgo::ziDip

private

Definition at line 127 of file DDHCalEndcapAlgo.h.

Referenced by getZiDip(), and initialize().

double DDHCalEndcapAlgo::ziKink

private

Definition at line 124 of file DDHCalEndcapAlgo.h.

Referenced by getZiKink(), and initialize().

double DDHCalEndcapAlgo::ziL0Body

private

Definition at line 123 of file DDHCalEndcapAlgo.h.

Referenced by getZiL0Body(), and initialize().

double DDHCalEndcapAlgo::ziL0Nose

private

Definition at line 121 of file DDHCalEndcapAlgo.h.

Referenced by getZiL0Nose(), and initialize().

double DDHCalEndcapAlgo::ziNose

private

Definition at line 120 of file DDHCalEndcapAlgo.h.

Referenced by getZiNose(), and initialize().

std::vector<double> DDHCalEndcapAlgo::zmaxBlock

private

Definition at line 161 of file DDHCalEndcapAlgo.h.

Referenced by getZmaxBlock(), and initialize().

std::vector<double> DDHCalEndcapAlgo::zminBlock

private

Definition at line 160 of file DDHCalEndcapAlgo.h.

Referenced by getZminBlock(), and initialize().

double DDHCalEndcapAlgo::zShift

private

Definition at line 116 of file DDHCalEndcapAlgo.h.

Referenced by getZShift(), and initialize().

double DDHCalEndcapAlgo::zShiftHac2

private

Definition at line 130 of file DDHCalEndcapAlgo.h.

Referenced by getZShiftHac2(), and initialize().

Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation