#include <DDHCalEndcapModuleAlgo.h>

Inheritance diagram for DDHCalEndcapModuleAlgo:

Classes
struct	HcalEndcapPar

Public Member Functions
	DDHCalEndcapModuleAlgo ()

void	execute (DDCompactView &cpv) override

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

	~DDHCalEndcapModuleAlgo () override

Protected Member Functions
void	constructInsideModule (const DDLogicalPart &module, DDCompactView &cpv)

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

void	constructScintLayer (const DDLogicalPart &detector, double dz, DDHCalEndcapModuleAlgo::HcalEndcapPar parm, const std::string &nm, int id, DDCompactView &cpv)

double	getRout (double z) const

double	getTrim (unsigned int j) const

HcalEndcapPar	parameterLayer (unsigned int iphi, double rinF, double routF, double rinB, double routB, double zi, double zo)

HcalEndcapPar	parameterLayer0 (unsigned int iphi)

Private Attributes
std::string	absorberMat

double	dzStep

std::string	genMaterial

std::string	idName

std::string	idNameSpace

int	idOffset

std::vector< std::string >	layerName

std::vector< int >	layerNumber

double	layerThick

int	layerType

std::string	modName

int	modNumber

int	modType

double	moduleThick

std::vector< std::string >	phiName

std::string	plasticMat

double	rMaxBack

double	rMaxFront

std::string	rotstr

std::string	scintMat

double	scintThick

int	sectors

double	slopeBot

double	slopeTop

double	slopeTopF

double	tolAbs

double	trimLeft

double	trimRight

double	z1Beam

double	ziDip

double	zMaxBlock

double	zMinBlock

Detailed Description

Definition at line 10 of file DDHCalEndcapModuleAlgo.h.

Constructor & Destructor Documentation

DDHCalEndcapModuleAlgo::DDHCalEndcapModuleAlgo ( )

Definition at line 21 of file DDHCalEndcapModuleAlgo.cc.

                                                {
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo info: Creating an instance";
 }

DDHCalEndcapModuleAlgo::~DDHCalEndcapModuleAlgo ( )

override

Definition at line 25 of file DDHCalEndcapModuleAlgo.cc.

25 {}

Member Function Documentation

void DDHCalEndcapModuleAlgo::constructInsideModule	(	const DDLogicalPart &	module,
		DDCompactView &	cpv
	)

protected

Definition at line 227 of file DDHCalEndcapModuleAlgo.cc.

Referenced by execute(), and DDHCalEndcapModuleAlgo::HcalEndcapPar::HcalEndcapPar().

                                                                                                   {
   
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: \t\tInside module";
 
   //Pointers to the Rotation Matrices and to the Materials
   DDRotation rot(DDName(DDSplit(rotstr).first, DDSplit(rotstr).second));
   DDName matName(DDSplit(genMaterial).first, DDSplit(genMaterial).second);
   DDMaterial matter(matName);
   DDName plasName(DDSplit(plasticMat).first, DDSplit(plasticMat).second);
   DDMaterial matplastic(plasName);
 
   double  alpha = CLHEP::pi/sectors;
   double  zi    = zMinBlock;
 
   for (unsigned int i=0; i<layerName.size(); i++) {
     std::string name;
     DDSolid solid;
     DDLogicalPart glog, plog;
     int     layer  = layerNumber[i];
     double  zo     = zi + 0.5*dzStep;
 
     for (unsigned int iphi=0; iphi<phiName.size(); iphi++) {
       double ziAir = zo - moduleThick;
       double rinF, rinB;
       if (modNumber == 0) {
     rinF  = ziAir * slopeTopF;
     rinB  = zo    * slopeTopF;
       } else {
     rinF  = ziAir * slopeBot;
     rinB  = zo    * slopeBot;
       }
       double routF = getRout(ziAir);
       double routB = getRout(zo);
       edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: Layer " << i 
                    << " Phi "  << iphi << " Front " << ziAir <<", "
                    << rinF << ", " << routF << " Back " << zo 
                    << ", " << rinB << ", " << routB;
       DDHCalEndcapModuleAlgo::HcalEndcapPar parm = parameterLayer(iphi, rinF, 
                                   routF, rinB, 
                                   routB, ziAir,
                                   zo);
       
       name  = idName+modName+layerName[i]+phiName[iphi]+"Air";
       solid = DDSolidFactory::trap(DDName(name, idNameSpace), 
                    0.5*moduleThick, parm.theta, parm.phi, 
                    parm.yh1, parm.bl1, parm.tl1, parm.alp, 
                    parm.yh2, parm.bl2, parm.tl2, parm.alp);
       edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: " 
                    << solid.name() << " Trap made of " << matName
                    << " of dimensions " << 0.5*moduleThick << ", "
                    << parm.theta/CLHEP::deg << ", " 
                    << parm.phi/CLHEP::deg << ", " << parm.yh1 
                    << ", " << parm.bl1 << ", " << parm.tl1 << ", "
                    << parm.alp/CLHEP::deg << ", " << parm.yh2 
                    << ", " << parm.bl2 << ", " << parm.tl2 << ", " 
                    << parm.alp/CLHEP::deg;
       glog = DDLogicalPart(solid.ddname(), matter, solid);
 
       DDTranslation r1(parm.xpos, parm.ypos, parm.zpos);
       cpv.position(glog, module, layer+1, r1, rot);
       edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: " <<glog.name()
                    << " number " << layer+1 << " positioned in " 
                    << module.name() << " at " << r1 << " with " 
                    << rot;
 
       //Now the plastic with scintillators
       parm.yh1 = 0.5 * (routF - rinB) - getTrim(iphi);
       parm.bl1 = 0.5 * rinB  * tan(alpha) - getTrim(iphi);
       parm.tl1 = 0.5 * routF * tan(alpha) - getTrim(iphi);
       name  = idName+modName+layerName[i]+phiName[iphi];
       solid = DDSolidFactory::trap(DDName(name, idNameSpace), 
                    0.5*layerThick, 0, 0, parm.yh1,
                    parm.bl1, parm.tl1, parm.alp, parm.yh1, 
                    parm.bl1, parm.tl1, parm.alp);
       edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: "<<solid.name()
                    << " Trap made of " << plasName
                    << " of dimensions " << 0.5*layerThick
                    << ", 0, 0, " << parm.yh1 << ", " << parm.bl1 
                    << ", " << parm.tl1 <<", "<< parm.alp/CLHEP::deg
                    << ", " << parm.yh1 << ", " << parm.bl1 << ", " 
                    << parm.tl1 << ", " << parm.alp/CLHEP::deg;
       plog = DDLogicalPart(solid.ddname(), matplastic, solid);
 
       double ypos = 0.5*(routF+rinB) - parm.xpos;
       DDTranslation r2(0., ypos, 0.);
       cpv.position(plog, glog, idOffset+layer+1, r2, DDRotation());
       edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: " <<plog.name()
                    << " number " << idOffset+layer+1 
                    << " positioned in " << glog.name() << " at " 
                    << r2 << " with no rotation";
 
       //Constructin the scintillators inside
       int copyNo = layer*10 + layerType;
       name = modName+layerName[i]+phiName[iphi];
       constructScintLayer (plog, scintThick, parm, name, copyNo, cpv);
       zo += 0.5*dzStep;
     } // End of loop over phi indices
     zi = zo - 0.5*dzStep;
   }   // End of loop on layers
 }

void DDHCalEndcapModuleAlgo::constructInsideModule0	(	const DDLogicalPart &	module,
		DDCompactView &	cpv
	)

protected

Definition at line 125 of file DDHCalEndcapModuleAlgo.cc.

Referenced by execute(), and DDHCalEndcapModuleAlgo::HcalEndcapPar::HcalEndcapPar().

                                                                                                    {
   
   edm::LogInfo("HCalGeom") <<"DDHCalEndcapModuleAlgo test: \t\tInside module0";
 
   //Pointers to the Rotation Matrices and to the Materials
   DDRotation rot(DDName(DDSplit(rotstr).first, DDSplit(rotstr).second));
   DDName matName(DDSplit(absorberMat).first, DDSplit(absorberMat).second);
   DDMaterial matabsorbr(matName);
   DDName plasName(DDSplit(plasticMat).first, DDSplit(plasticMat).second);
   DDMaterial matplastic(plasName);
 
   int           layer  = layerNumber[0];
   int           layer0 = layerNumber[1];
   std::string   name;
   DDSolid       solid;
   DDLogicalPart glog, plog;
   for (unsigned int iphi=0; iphi<phiName.size(); iphi++) {
     DDHCalEndcapModuleAlgo::HcalEndcapPar parm = parameterLayer0(iphi);
     name  = idName+modName+layerName[0]+phiName[iphi];
     solid = DDSolidFactory::trap(DDName(name, idNameSpace), 
                  0.5*layerThick, 0, 0, parm.yh1, parm.bl1, 
                  parm.tl1, parm.alp, parm.yh2, parm.bl1, 
                  parm.tl2, parm.alp);
     edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: " << solid.name()
                  << " Trap made of " << plasName
                  << " of dimensions " << 0.5*layerThick
                  << ", 0, 0, " << parm.yh1 << ", " << parm.bl1 
                  << ", " << parm.tl1 << ", " << parm.alp/CLHEP::deg
                  << ", " << parm.yh2 << ", " << parm.bl2 << ", " 
                  << parm.tl2 << ", " << parm.alp/CLHEP::deg;
     glog = DDLogicalPart(solid.ddname(), matplastic, solid);
 
     DDTranslation r1(parm.xpos, parm.ypos, parm.zpos);
     cpv.position(glog, module, idOffset+layer+1, r1, rot);
     edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: " << glog.name() 
              << " number " << idOffset+layer+1 << " positioned in "
              << module.name() << " at " << r1 << " with " << rot;
     //Now construct the layer of scintillator inside this
     int copyNo = layer0*10 + layerType;
     name = modName+layerName[0]+phiName[iphi];
     constructScintLayer (glog, scintThick, parm, name, copyNo, cpv);
   }
 
   //Now the absorber layer
   double zi = zMinBlock + layerThick;
   double zo = zi + 0.5*dzStep;
   double rinF, routF, rinB, routB;
   if (modNumber == 0) {
     rinF  = zi * slopeTopF;
     routF = (zi - z1Beam) * slopeTop;
     rinB  = zo * slopeTopF;
     routB = (zo - z1Beam) * slopeTop;
   } else if (modNumber > 0) {
     rinF  = zi * slopeBot;
     routF = zi * slopeTopF;
     rinB  = zo * slopeBot;
     routB = zo * slopeTopF;
   } else {
     rinF  = zi * slopeBot;
     routF = (zi - z1Beam) * slopeTop;
     rinB  = zo * slopeBot;
     routB = (zo - z1Beam) * slopeTop;
   }
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: Front "
                << zi << ", " << rinF << ", " << routF << " Back "
                << zo << ", " << rinB << ", " << routB;
   DDHCalEndcapModuleAlgo::HcalEndcapPar parm = parameterLayer(0, rinF, routF, 
                                   rinB, routB, zi, 
                                   zo);
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: Trim " << tolAbs
                << " Param " << parm.yh1 << ", " << parm.bl1 << ", "
                << parm.tl1 << ", " << parm.yh2 << ", " << parm.bl2 
                << ", " << parm.tl2;
   parm.bl1 -= tolAbs;
   parm.tl1 -= tolAbs;
   parm.bl2 -= tolAbs;
   parm.tl2 -= tolAbs;
 
   name  = idName+modName+layerName[0]+"Absorber";
   solid = DDSolidFactory::trap(DDName(name, idNameSpace), 
                    0.5*moduleThick, parm.theta, parm.phi, parm.yh1,
                    parm.bl1, parm.tl1, parm.alp, parm.yh2, 
                    parm.bl2, parm.tl2, parm.alp);
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: " << solid.name() 
                << " Trap made of " << matName << " of dimensions " 
                << 0.5*moduleThick << ", " << parm.theta/CLHEP::deg 
                << ", " << parm.phi/CLHEP::deg << ", " << parm.yh1 
                << ", " << parm.bl1 << ", "  << parm.tl1 << ", " 
                << parm.alp/CLHEP::deg << ", " << parm.yh2 << ", "
                << parm.bl2 << ", " << parm.tl2 << ", " 
                << parm.alp/CLHEP::deg;
   glog = DDLogicalPart(solid.ddname(), matabsorbr, solid);
 
   DDTranslation r2(parm.xpos, parm.ypos, parm.zpos);
   cpv.position(glog, module, 1, r2, rot);
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: " << glog.name() 
                << " number 1 positioned in " << module.name() 
                << " at " << r2 << " with " << rot;
 }

void DDHCalEndcapModuleAlgo::constructScintLayer	(	const DDLogicalPart &	detector,
		double	dz,
		DDHCalEndcapModuleAlgo::HcalEndcapPar	parm,
		const std::string &	nm,
		int	id,
		DDCompactView &	cpv
	)

protected

Definition at line 431 of file DDHCalEndcapModuleAlgo.cc.

References DDHCalEndcapModuleAlgo::HcalEndcapPar::alp, DDHCalEndcapModuleAlgo::HcalEndcapPar::bl1, DDBase< N, C >::ddname(), DDSplit(), plotBeamSpotDB::first, idName, idNameSpace, dataset::name, DDBase< N, C >::name(), DDCompactView::position(), scintMat, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, DDHCalEndcapModuleAlgo::HcalEndcapPar::tl1, DDSolidFactory::trap(), and DDHCalEndcapModuleAlgo::HcalEndcapPar::yh1.

Referenced by constructInsideModule(), constructInsideModule0(), and DDHCalEndcapModuleAlgo::HcalEndcapPar::HcalEndcapPar().

                                                                           {
 
   DDName matname(DDSplit(scintMat).first, DDSplit(scintMat).second);
   DDMaterial matter(matname);
   std::string name = idName+"Scintillator"+nm;
 
   DDSolid solid = DDSolidFactory::trap(DDName(name, idNameSpace), 0.5*dz, 0, 0,
                        parm.yh1, parm.bl1, parm.tl1, parm.alp, 
                        parm.yh1, parm.bl1, parm.tl1, parm.alp);
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: " << solid.name()
                << " Trap made of " << scintMat << " of dimensions "
                << 0.5*dz << ", 0, 0, " << parm.yh1 << ", "  
                << parm.bl1 << ", " << parm.tl1 << ", " 
                << parm.alp/CLHEP::deg << ", " << parm.yh1 << ", " 
                << parm.bl1 << ", " << parm.tl1 << ", " 
                << parm.alp/CLHEP::deg;
 
   DDLogicalPart glog(solid.ddname(), matter, solid); 
 
   cpv.position(glog, detector, id, DDTranslation(0,0,0), DDRotation());
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo test: " << glog.name() 
                << " number " << id << " positioned in " 
                << detector.name() <<" at (0,0,0) with no rotation";
 
 }

void DDHCalEndcapModuleAlgo::execute ( DDCompactView & cpv )

override

Definition at line 114 of file DDHCalEndcapModuleAlgo.cc.

References constructInsideModule(), constructInsideModule0(), modType, and class-composition::parent.

Referenced by DDHCalEndcapModuleAlgo::HcalEndcapPar::HcalEndcapPar().

                                                        {
   
   edm::LogInfo("HCalGeom") << "==>> Constructing DDHCalEndcapModuleAlgo...";
   if (modType == 0) 
     constructInsideModule0 (parent(), cpv);
   else
     constructInsideModule  (parent(), cpv);
   edm::LogInfo("HCalGeom") << "<<== End of DDHCalEndcapModuleAlgo construction ...";
 }

double DDHCalEndcapModuleAlgo::getRout ( double z ) const

protected

Definition at line 465 of file DDHCalEndcapModuleAlgo.cc.

References modNumber, alignCSCRings::r, rMaxBack, rMaxFront, slopeTop, slopeTopF, z1Beam, and ziDip.

Referenced by constructInsideModule(), and DDHCalEndcapModuleAlgo::HcalEndcapPar::HcalEndcapPar().

                                                      {
   double r = (modNumber >= 0) ? ((z - z1Beam) * slopeTop) : z * slopeTopF;
   if (z > ziDip) {
     if (r > rMaxBack)  r = rMaxBack;
   } else {
     if (r > rMaxFront) r = rMaxFront;
   }
   return r;
 }

double DDHCalEndcapModuleAlgo::getTrim ( unsigned int j ) const

protected

Definition at line 459 of file DDHCalEndcapModuleAlgo.cc.

References trimLeft, and trimRight.

Referenced by constructInsideModule(), DDHCalEndcapModuleAlgo::HcalEndcapPar::HcalEndcapPar(), and parameterLayer0().

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

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

override

Definition at line 27 of file DDHCalEndcapModuleAlgo.cc.

References absorberMat, dbl_to_int(), dzStep, genMaterial, mps_fire::i, idName, idNameSpace, idOffset, createfilelist::int, layerName, layerNumber, layerThick, layerType, modName, modNumber, modType, moduleThick, DDCurrentNamespace::ns(), class-composition::parent, phiName, plasticMat, rMaxBack, rMaxFront, rotstr, scintMat, scintThick, sectors, slopeBot, slopeTop, slopeTopF, tolAbs, trimLeft, trimRight, z1Beam, ziDip, zMaxBlock, and zMinBlock.

Referenced by DDHCalEndcapModuleAlgo::HcalEndcapPar::HcalEndcapPar().

                                                           {
 
   genMaterial   = sArgs["MaterialName"];
   absorberMat   = sArgs["AbsorberMat"];
   plasticMat    = sArgs["PlasticMat"];
   scintMat      = sArgs["ScintMat"];
   rotstr        = sArgs["Rotation"];
   sectors       = int (nArgs["Sectors"]);
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo: General material " 
                << genMaterial << "\tAbsorber " << absorberMat
                << "\tPlastic " << plasticMat << "\tScintillator "
                << scintMat << "\tRotation " << rotstr
                << "\tSectors "  << sectors;
 
   zMinBlock     = nArgs["ZMinBlock"];
   zMaxBlock     = nArgs["ZMaxBlock"];
   z1Beam        = nArgs["Z1Beam"];
   ziDip         = nArgs["ZiDip"];
   dzStep        = nArgs["DzStep"];
   moduleThick   = nArgs["ModuleThick"];
   layerThick    = nArgs["LayerThick"];
   scintThick    = nArgs["ScintThick"];
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo: Zmin " << zMinBlock
                << "\tZmax " << zMaxBlock << "\tZ1Beam " << z1Beam
                << "\tZiDip " << ziDip << "\tDzStep " << dzStep
                << "\tModuleThick " << moduleThick <<"\tLayerThick "
                << layerThick << "\tScintThick " << scintThick;
 
   rMaxFront     = nArgs["RMaxFront"];
   rMaxBack      = nArgs["RMaxBack"];
   trimLeft      = nArgs["TrimLeft"];
   trimRight     = nArgs["TrimRight"];
   tolAbs        = nArgs["TolAbs"];
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo: RMaxFront " << rMaxFront
                <<"\tRmaxBack " << rMaxBack << "\tTrims " <<trimLeft
                << ":" << trimRight << "\tTolAbs " << tolAbs;
 
   slopeBot      = nArgs["SlopeBottom"];
   slopeTop      = nArgs["SlopeTop"];
   slopeTopF     = nArgs["SlopeTopFront"];
   modType       = (int)(nArgs["ModType"]);
   modNumber     = (int)(nArgs["ModNumber"]);
   layerType     = (int)(nArgs["LayerType"]);
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo: slopeBot " << slopeBot
                << "\tslopeTop " << slopeTop << "\tslopeTopF "
                << slopeTopF << "\tmodType " << modType
                << "\tmodNumber " << modNumber << "\tlayerType "
                << layerType;
 
   layerNumber   = dbl_to_int(vArgs["LayerNumber"]);
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo: " << layerNumber.size()
                << " layer Numbers";
   for (unsigned int i=0; i<layerNumber.size(); ++i)
     edm::LogInfo("HCalGeom") << "LayerNumber[" << i << "] = " <<layerNumber[i];
 
   phiName       = vsArgs["PhiName"];
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo: " << phiName.size()
                << " phi sectors";
   for (unsigned int i=0; i<phiName.size(); ++i)
     edm::LogInfo("HCalGeom") << "PhiName[" << i << "] = " << phiName[i];
 
   layerName     = vsArgs["LayerName"];
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo: " << layerName.size()
                << " layers";
   for (unsigned int i=0; i<layerName.size(); ++i)
     edm::LogInfo("HCalGeom") << "LayerName[" << i << "] = " << layerName[i];
 
   idName      = sArgs["MotherName"];
   idNameSpace = DDCurrentNamespace::ns();
   idOffset    = int (nArgs["IdOffset"]); 
   DDName parentName = parent().name(); 
   modName     = sArgs["ModName"];
   edm::LogInfo("HCalGeom") << "DDHCalEndcapModuleAlgo debug: Parent " 
                << parentName << "   " << modName << " idName " 
                << idName << " NameSpace " << idNameSpace
                << " Offset " << idOffset;
 
 }

DDHCalEndcapModuleAlgo::HcalEndcapPar DDHCalEndcapModuleAlgo::parameterLayer	(	unsigned int	iphi,
		double	rinF,
		double	routF,
		double	rinB,
		double	routB,
		double	zi,
		double	zo
	)

protected

Definition at line 383 of file DDHCalEndcapModuleAlgo.cc.

Referenced by constructInsideModule(), constructInsideModule0(), and DDHCalEndcapModuleAlgo::HcalEndcapPar::HcalEndcapPar().

                                              {
 
   DDHCalEndcapModuleAlgo::HcalEndcapPar parm;
   //Given rin, rout compute parameters of the trapezoid and 
   //position of the trapezoid for a standrd layer
   double alpha = CLHEP::pi/sectors;
   edm::LogInfo("HCalGeom") << "Input " << iphi << " Front " << rinF << " " 
                << routF << " " << zi << " Back " << rinB << " " 
                << routB << " " << zo<<" Alpha " <<alpha/CLHEP::deg;
 
   parm.yh1    = 0.5 * (routF - rinB);
   parm.bl1    = 0.5 * rinB  * tan(alpha);
   parm.tl1    = 0.5 * routF * tan(alpha);
   parm.yh2    = 0.5 * (routF - rinB);
   parm.bl2    = 0.5 * rinB  * tan(alpha);
   parm.tl2    = 0.5 * routF * tan(alpha);
   double dx   = 0.25* (parm.bl2+parm.tl2-parm.bl1-parm.tl1);
   double dy   = 0.5 * (rinB+routF-rinB-routF);
   parm.xpos   = 0.25*(rinB+routF+rinB+routF);
   parm.ypos   = 0.25*(parm.bl2+parm.tl2+parm.bl1+parm.tl1);
   parm.zpos   = 0.5*(zi+zo);
   parm.alp    = atan(0.5 * tan(alpha));
   if (iphi == 0) {
     parm.ypos = -parm.ypos;
   } else {
     parm.alp  = -parm.alp;
     dx        = -dx;
   }
   double r    = sqrt (dx*dx + dy*dy);
   edm::LogInfo("HCalGeom") << "dx|dy|r " << dx << ":" << dy << ":" << r;
   if (r > 1.0e-8) {
     parm.theta  = atan (r/(zo-zi));
     parm.phi    = atan2 (dy, dx);
   } else {
     parm.theta  = parm.phi = 0;
   }
   edm::LogInfo("HCalGeom") << "Output Dimensions " << parm.yh1 << " " 
                << parm.bl1 << " " << parm.tl1 << " " << parm.yh2 
                << " " << parm.bl2 << " " << parm.tl2 << " " 
                << parm.alp/CLHEP::deg <<" " <<parm.theta/CLHEP::deg
                << " " << parm.phi/CLHEP::deg << " Position " 
                << parm.xpos << " " << parm.ypos << " " <<parm.zpos;
   return parm;
 }

DDHCalEndcapModuleAlgo::HcalEndcapPar DDHCalEndcapModuleAlgo::parameterLayer0 ( unsigned int iphi )

protected

Definition at line 331 of file DDHCalEndcapModuleAlgo.cc.

Referenced by constructInsideModule0(), and DDHCalEndcapModuleAlgo::HcalEndcapPar::HcalEndcapPar().

                                                          {
 
   DDHCalEndcapModuleAlgo::HcalEndcapPar parm;
   //Given module and layer number compute parameters of trapezoid
   //and positioning parameters
   double alpha = CLHEP::pi/sectors;
   edm::LogInfo("HCalGeom") << "Input " << iphi << " Alpha " <<alpha/CLHEP::deg;
 
   double zi, zo;
   if (iphi == 0) {
     zi = zMinBlock;
     zo = zi + layerThick;
   } else {
     zo = zMaxBlock;
     zi = zo - layerThick;
   }
   double rin, rout;
   if (modNumber == 0) {
     rin  = zo * slopeTopF;
     rout = (zi - z1Beam) * slopeTop;
   } else if (modNumber > 0) {
     rin  = zo * slopeBot;
     rout = zi * slopeTopF;
   } else {
     rin  = zo * slopeBot;
     rout = (zi - z1Beam) * slopeTop;
   }
   edm::LogInfo("HCalGeom") << "ModNumber " << modNumber << " " << zi << " " << zo << " " << slopeTopF << " " << slopeTop << " " << slopeBot << " " << rin << " " << rout << " " << getTrim(iphi);
   double yh   = 0.5 * (rout - rin);
   double bl   = 0.5 * rin * tan (alpha);
   double tl   = 0.5 * rout * tan(alpha);
   parm.xpos = 0.5 * (rin + rout);
   parm.ypos = 0.5 * (bl + tl);
   parm.zpos = 0.5 * (zi + zo);
   parm.yh1  = parm.yh2 = yh - getTrim(iphi);
   parm.bl1  = parm.bl2 = bl - getTrim(iphi);
   parm.tl1  = parm.tl2 = tl - getTrim(iphi);
   parm.alp  = atan(0.5 * tan(alpha));
   if (iphi == 0) {
     parm.ypos = -parm.ypos;
   } else {
     parm.alp  = -parm.alp;
   }
   edm::LogInfo("HCalGeom") << "Output Dimensions " << parm.yh1 << " " 
                << parm.bl1 << " " << parm.tl1 << " " 
                << parm.alp/CLHEP::deg << " Position " << parm.xpos 
                << " " << parm.ypos << " " << parm.zpos;
   return parm;
 }