DDHCalEndcapModuleAlgo Class Reference

#include <DDHCalEndcapModuleAlgo.h>

Inheritance diagram for DDHCalEndcapModuleAlgo:

Classes
struct	HcalEndcapPar

Public Member Functions
	DDHCalEndcapModuleAlgo ()
void	execute (DDCompactView &cpv)
void	initialize (const DDNumericArguments &nArgs, const DDVectorArguments &vArgs, const DDMapArguments &mArgs, const DDStringArguments &sArgs, const DDStringVectorArguments &vsArgs)
virtual	~DDHCalEndcapModuleAlgo ()

Protected Member Functions
void	constructInsideModule (DDLogicalPart module, DDCompactView &cpv)
void	constructInsideModule0 (DDLogicalPart module, DDCompactView &cpv)
void	constructScintLayer (DDLogicalPart detector, double dz, DDHCalEndcapModuleAlgo::HcalEndcapPar parm, 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 ( )

virtual

Definition at line 25 of file DDHCalEndcapModuleAlgo.cc.

{}

Member Function Documentation

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

protected

Definition at line 227 of file DDHCalEndcapModuleAlgo.cc.

References DDHCalEndcapModuleAlgo::HcalEndcapPar::alp, alpha, DDHCalEndcapModuleAlgo::HcalEndcapPar::bl1, DDHCalEndcapModuleAlgo::HcalEndcapPar::bl2, constructScintLayer(), DDBase< N, C >::ddname(), DDSplit(), dzStep, plotBeamSpotDB::first, genMaterial, getRout(), getTrim(), mps_fire::i, idName, idNameSpace, idOffset, layerName, layerNumber, layerThick, layerType, modName, modNumber, moduleThick, dataset::name, DDBase< N, C >::name(), parameterLayer(), DDHCalEndcapModuleAlgo::HcalEndcapPar::phi, phiName, pi, plasticMat, DDCompactView::position(), diffTwoXMLs::r1, diffTwoXMLs::r2, makeMuonMisalignmentScenario::rot, rotstr, scintThick, edm::second(), sectors, slopeBot, slopeTopF, AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), DDHCalEndcapModuleAlgo::HcalEndcapPar::theta, DDHCalEndcapModuleAlgo::HcalEndcapPar::tl1, DDHCalEndcapModuleAlgo::HcalEndcapPar::tl2, DDSolidFactory::trap(), DDHCalEndcapModuleAlgo::HcalEndcapPar::xpos, DDHCalEndcapModuleAlgo::HcalEndcapPar::yh1, DDHCalEndcapModuleAlgo::HcalEndcapPar::yh2, DDHCalEndcapModuleAlgo::HcalEndcapPar::ypos, zMinBlock, and DDHCalEndcapModuleAlgo::HcalEndcapPar::zpos.

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 ( DDLogicalPart  module, DDCompactView &  cpv  )

protected

Definition at line 125 of file DDHCalEndcapModuleAlgo.cc.

References absorberMat, DDHCalEndcapModuleAlgo::HcalEndcapPar::alp, DDHCalEndcapModuleAlgo::HcalEndcapPar::bl1, DDHCalEndcapModuleAlgo::HcalEndcapPar::bl2, constructScintLayer(), DDBase< N, C >::ddname(), DDSplit(), dzStep, plotBeamSpotDB::first, idName, idNameSpace, idOffset, layerName, layerNumber, layerThick, layerType, modName, modNumber, moduleThick, dataset::name, DDBase< N, C >::name(), parameterLayer(), parameterLayer0(), DDHCalEndcapModuleAlgo::HcalEndcapPar::phi, phiName, plasticMat, DDCompactView::position(), diffTwoXMLs::r1, diffTwoXMLs::r2, makeMuonMisalignmentScenario::rot, rotstr, scintThick, edm::second(), slopeBot, slopeTop, slopeTopF, AlCaHLTBitMon_QueryRunRegistry::string, DDHCalEndcapModuleAlgo::HcalEndcapPar::theta, DDHCalEndcapModuleAlgo::HcalEndcapPar::tl1, DDHCalEndcapModuleAlgo::HcalEndcapPar::tl2, tolAbs, DDSolidFactory::trap(), DDHCalEndcapModuleAlgo::HcalEndcapPar::xpos, DDHCalEndcapModuleAlgo::HcalEndcapPar::yh1, DDHCalEndcapModuleAlgo::HcalEndcapPar::yh2, DDHCalEndcapModuleAlgo::HcalEndcapPar::ypos, z1Beam, zMinBlock, and DDHCalEndcapModuleAlgo::HcalEndcapPar::zpos.

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 ( DDLogicalPart  detector, double  dz, DDHCalEndcapModuleAlgo::HcalEndcapPar  parm, 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

)

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
	)

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.

References DDHCalEndcapModuleAlgo::HcalEndcapPar::alp, alpha, DDHCalEndcapModuleAlgo::HcalEndcapPar::bl1, DDHCalEndcapModuleAlgo::HcalEndcapPar::bl2, MillePedeFileConverter_cfg::e, DDHCalEndcapModuleAlgo::HcalEndcapPar::phi, pi, alignCSCRings::r, sectors, mathSSE::sqrt(), funct::tan(), DDHCalEndcapModuleAlgo::HcalEndcapPar::theta, DDHCalEndcapModuleAlgo::HcalEndcapPar::tl1, DDHCalEndcapModuleAlgo::HcalEndcapPar::tl2, DDHCalEndcapModuleAlgo::HcalEndcapPar::xpos, DDHCalEndcapModuleAlgo::HcalEndcapPar::yh1, DDHCalEndcapModuleAlgo::HcalEndcapPar::yh2, DDHCalEndcapModuleAlgo::HcalEndcapPar::ypos, and DDHCalEndcapModuleAlgo::HcalEndcapPar::zpos.

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.

References DDHCalEndcapModuleAlgo::HcalEndcapPar::alp, alpha, DDHCalEndcapModuleAlgo::HcalEndcapPar::bl1, DDHCalEndcapModuleAlgo::HcalEndcapPar::bl2, getTrim(), layerThick, modNumber, pi, sectors, slopeBot, slopeTop, slopeTopF, funct::tan(), DDHCalEndcapModuleAlgo::HcalEndcapPar::tl1, DDHCalEndcapModuleAlgo::HcalEndcapPar::tl2, DDHCalEndcapModuleAlgo::HcalEndcapPar::xpos, DDHCalEndcapModuleAlgo::HcalEndcapPar::yh1, DDHCalEndcapModuleAlgo::HcalEndcapPar::yh2, DDHCalEndcapModuleAlgo::HcalEndcapPar::ypos, z1Beam, zMaxBlock, zMinBlock, and DDHCalEndcapModuleAlgo::HcalEndcapPar::zpos.

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

Member Data Documentation

std::string DDHCalEndcapModuleAlgo::absorberMat

private

Definition at line 48 of file DDHCalEndcapModuleAlgo.h.

Referenced by constructInsideModule0(), and initialize().

double DDHCalEndcapModuleAlgo::dzStep

private

Definition at line 57 of file DDHCalEndcapModuleAlgo.h.

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

std::string DDHCalEndcapModuleAlgo::genMaterial

private

Definition at line 47 of file DDHCalEndcapModuleAlgo.h.

Referenced by constructInsideModule(), and initialize().

std::string DDHCalEndcapModuleAlgo::idName

private

Definition at line 76 of file DDHCalEndcapModuleAlgo.h.

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

std::string DDHCalEndcapModuleAlgo::idNameSpace

private

Definition at line 77 of file DDHCalEndcapModuleAlgo.h.

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

int DDHCalEndcapModuleAlgo::idOffset

private

Definition at line 79 of file DDHCalEndcapModuleAlgo.h.

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

std::vector<std::string> DDHCalEndcapModuleAlgo::layerName

private

Definition at line 74 of file DDHCalEndcapModuleAlgo.h.

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

std::vector<int> DDHCalEndcapModuleAlgo::layerNumber

private

Definition at line 72 of file DDHCalEndcapModuleAlgo.h.

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

double DDHCalEndcapModuleAlgo::layerThick

private

Definition at line 59 of file DDHCalEndcapModuleAlgo.h.

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

int DDHCalEndcapModuleAlgo::layerType

private

Definition at line 71 of file DDHCalEndcapModuleAlgo.h.

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

std::string DDHCalEndcapModuleAlgo::modName

private

Definition at line 78 of file DDHCalEndcapModuleAlgo.h.

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

int DDHCalEndcapModuleAlgo::modNumber

private

Definition at line 70 of file DDHCalEndcapModuleAlgo.h.

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

int DDHCalEndcapModuleAlgo::modType

private

Definition at line 69 of file DDHCalEndcapModuleAlgo.h.

Referenced by execute(), and initialize().

double DDHCalEndcapModuleAlgo::moduleThick

private

Definition at line 58 of file DDHCalEndcapModuleAlgo.h.

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

std::vector<std::string> DDHCalEndcapModuleAlgo::phiName

private

Definition at line 73 of file DDHCalEndcapModuleAlgo.h.

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

std::string DDHCalEndcapModuleAlgo::plasticMat

private

Definition at line 49 of file DDHCalEndcapModuleAlgo.h.

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

double DDHCalEndcapModuleAlgo::rMaxBack

private

Definition at line 61 of file DDHCalEndcapModuleAlgo.h.

Referenced by getRout(), and initialize().

double DDHCalEndcapModuleAlgo::rMaxFront

private

Definition at line 62 of file DDHCalEndcapModuleAlgo.h.

Referenced by getRout(), and initialize().

std::string DDHCalEndcapModuleAlgo::rotstr

private

Definition at line 51 of file DDHCalEndcapModuleAlgo.h.

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

std::string DDHCalEndcapModuleAlgo::scintMat

private

Definition at line 50 of file DDHCalEndcapModuleAlgo.h.

Referenced by constructScintLayer(), and initialize().

double DDHCalEndcapModuleAlgo::scintThick

private

Definition at line 60 of file DDHCalEndcapModuleAlgo.h.

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

int DDHCalEndcapModuleAlgo::sectors

private

Definition at line 52 of file DDHCalEndcapModuleAlgo.h.

Referenced by constructInsideModule(), initialize(), parameterLayer(), and parameterLayer0().

double DDHCalEndcapModuleAlgo::slopeBot

private

Definition at line 63 of file DDHCalEndcapModuleAlgo.h.

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

double DDHCalEndcapModuleAlgo::slopeTop

private

Definition at line 64 of file DDHCalEndcapModuleAlgo.h.

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

double DDHCalEndcapModuleAlgo::slopeTopF

private

Definition at line 65 of file DDHCalEndcapModuleAlgo.h.

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

double DDHCalEndcapModuleAlgo::tolAbs

private

Definition at line 68 of file DDHCalEndcapModuleAlgo.h.

Referenced by constructInsideModule0(), and initialize().

double DDHCalEndcapModuleAlgo::trimLeft

private

Definition at line 66 of file DDHCalEndcapModuleAlgo.h.

Referenced by getTrim(), and initialize().

double DDHCalEndcapModuleAlgo::trimRight

private

Definition at line 67 of file DDHCalEndcapModuleAlgo.h.

Referenced by getTrim(), and initialize().

double DDHCalEndcapModuleAlgo::z1Beam

private

Definition at line 55 of file DDHCalEndcapModuleAlgo.h.

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

double DDHCalEndcapModuleAlgo::ziDip

private

Definition at line 56 of file DDHCalEndcapModuleAlgo.h.

Referenced by getRout(), and initialize().

double DDHCalEndcapModuleAlgo::zMaxBlock

private

Definition at line 54 of file DDHCalEndcapModuleAlgo.h.

Referenced by initialize(), and parameterLayer0().

double DDHCalEndcapModuleAlgo::zMinBlock

private

Definition at line 53 of file DDHCalEndcapModuleAlgo.h.

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