DDHGCalEEAlgo Class Reference

#include <DDHGCalEEAlgo.h>

Inheritance diagram for DDHGCalEEAlgo:

Classes
struct	HGCalEEPar

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

Protected Member Functions
void	constructLayers (const DDLogicalPart &, DDCompactView &cpv)
HGCalEEPar	parameterLayer (double rinF, double routF, double rinB, double routB, double zi, double zo)
double	rMax (double z)

Private Attributes
std::vector< int >	copyNumber
std::vector< int >	heightType
std::string	idName
std::string	idNameSpace
std::vector< int >	layerType
std::vector< std::string >	materials
std::vector< std::string >	names
std::vector< double >	rMaxFront
std::string	rotstr
int	sectors
double	slopeB
std::vector< double >	slopeT
std::vector< double >	thick
std::vector< double >	thickBlock
std::vector< double >	zFront
double	zMinBlock

Detailed Description

Definition at line 10 of file DDHGCalEEAlgo.h.

Constructor & Destructor Documentation

DDHGCalEEAlgo::DDHGCalEEAlgo

(

)

Definition at line 20 of file DDHGCalEEAlgo.cc.

                             {
  edm::LogInfo("HGCalGeom") << "DDHGCalEEAlgo info: Creating an instance";
}

DDHGCalEEAlgo::~DDHGCalEEAlgo ( )

override

Definition at line 24 of file DDHGCalEEAlgo.cc.

{}

Member Function Documentation

void DDHGCalEEAlgo::constructLayers ( const DDLogicalPart &  module, DDCompactView &  cpv  )

protected

Definition at line 89 of file DDHGCalEEAlgo.cc.

References DDHGCalEEAlgo::HGCalEEPar::alp, DDHGCalEEAlgo::HGCalEEPar::bl1, DDHGCalEEAlgo::HGCalEEPar::bl2, popcon2dropbox::copy(), copyNumber, DDBase< N, C >::ddname(), DDSplit(), plotBeamSpotDB::first, heightType, mps_fire::i, idNameSpace, cuy::ii, layerType, materials, dataset::name, DDBase< N, C >::name(), names, parameterLayer(), DDHGCalEEAlgo::HGCalEEPar::phi, DDCompactView::position(), diffTwoXMLs::r1, rMax(), makeMuonMisalignmentScenario::rot, rotstr, edm::second(), slopeB, AlCaHLTBitMon_QueryRunRegistry::string, DDHGCalEEAlgo::HGCalEEPar::theta, thick, thickBlock, DDHGCalEEAlgo::HGCalEEPar::tl1, DDHGCalEEAlgo::HGCalEEPar::tl2, DDSolidFactory::trap(), DDHGCalEEAlgo::HGCalEEPar::xpos, DDHGCalEEAlgo::HGCalEEPar::yh1, DDHGCalEEAlgo::HGCalEEPar::yh2, DDHGCalEEAlgo::HGCalEEPar::ypos, zMinBlock, DDHGCalEEAlgo::HGCalEEPar::zpos, and geometryCSVtoXML::zz.

Referenced by execute(), and DDHGCalEEAlgo::HGCalEEPar::HGCalEEPar().

                                                                                   {
  
  edm::LogInfo("HGCalGeom") << "DDHGCalEEAlgo test: \t\tInside Layers";

  //Pointers to the Rotation Matrices and to the Materials
  DDRotation rot(DDName(DDSplit(rotstr).first, DDSplit(rotstr).second));

  double  zi(zMinBlock), zz(zMinBlock);
  for (unsigned int i=0; i<layerType.size(); i++) {
    int     ii    = layerType[i];
    int     copy  = copyNumber[ii];
    ++copyNumber[ii];
    double layer_thick = thickBlock[i];

    if (heightType[i] == 0) zz = zi;
    double  zlayer = zz + layer_thick;
    double  zo     = zi + thick[ii];
    double  rinF   = zi * slopeB;
    double  rinB   = zlayer * slopeB;
    double  routF  = (heightType[i] == 0) ? rMax(zi) : rMax(zz);
    double  routB  = rMax(zo);

    std::string name = "HGCal"+names[ii]+std::to_string(copy);
    edm::LogInfo("HGCalGeom") << "DDHGCalEEAlgo test: Layer " << i << ":" 
                  << ii << " Front " << zi << ", " << rinF << ", " 
                  << routF << " Back " << zo << ", " << rinB 
                  << ", " << routB << " superlayer thickness " 
                  << layer_thick;
    DDHGCalEEAlgo::HGCalEEPar parm = parameterLayer(rinF, routF, rinB, 
                            routB, zi, zo);
    DDSolid solid = DDSolidFactory::trap(DDName(name, idNameSpace), 
                     0.5*thick[ii], parm.theta,
                     parm.phi, parm.yh1, parm.bl1, 
                     parm.tl1, parm.alp, parm.yh2,
                     parm.bl2, parm.tl2, parm.alp);

    DDName matName(DDSplit(materials[ii]).first, 
           DDSplit(materials[ii]).second);
    DDMaterial matter(matName);
    DDLogicalPart glog = DDLogicalPart(solid.ddname(), matter, solid);
    edm::LogInfo("HGCalGeom") << "DDHGCalEEAlgo test: " 
                  << solid.name() << " Trap made of " << matName
                  << " of dimensions " << 0.5*thick[ii] << ", "
                  << 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;
    DDTranslation r1(parm.xpos, parm.ypos, parm.zpos);
    cpv.position(glog, module, copy, r1, rot);
    edm::LogInfo("HGCalGeom") << "DDHGCalEEAlgo test: " << glog.name()
                  << " number " << copy << " positioned in " 
                  << module.name() << " at " << r1 << " with " 
                  << rot;
    zi = zo;
  }   // End of loop on layers
}

void DDHGCalEEAlgo::execute ( DDCompactView &  cpv )

override

Definition at line 82 of file DDHGCalEEAlgo.cc.

References constructLayers(), and class-composition::parent.

Referenced by DDHGCalEEAlgo::HGCalEEPar::HGCalEEPar().

                                              {
  
  edm::LogInfo("HGCalGeom") << "==>> Constructing DDHGCalEEAlgo...";
  constructLayers (parent(), cpv);
  edm::LogInfo("HGCalGeom") << "<<== End of DDHGCalEEAlgo construction ...";
}

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

override

Definition at line 26 of file DDHGCalEEAlgo.cc.

References copyNumber, dbl_to_int(), heightType, mps_fire::i, idName, idNameSpace, createfilelist::int, layerType, materials, names, DDCurrentNamespace::ns(), class-composition::parent, rMaxFront, rotstr, sectors, slopeB, slopeT, thick, thickBlock, zFront, and zMinBlock.

Referenced by DDHGCalEEAlgo::HGCalEEPar::HGCalEEPar().

                                                         {

  materials     = vsArgs["MaterialNames"];
  names         = vsArgs["VolumeNames"];
  thick         = vArgs["Thickness"];
  rotstr        = sArgs["Rotation"];
  layerType     = dbl_to_int(vArgs["LayerType"]);
  heightType    = dbl_to_int(vArgs["HeightType"]);
  thickBlock    = vArgs["LayerThick"];
  zMinBlock     = nArgs["zMinBlock"];
  for (unsigned int i=0; i<materials.size(); ++i) {
    copyNumber.emplace_back(1);
  }
  edm::LogInfo("HGCalGeom") << "DDHGCalEEAlgo: " << materials.size()
                << " volumes to be put with rotation " << rotstr
                << " starting at " << zMinBlock;
  for (unsigned int i=0; i<names.size(); ++i)
    edm::LogInfo("HGCalGeom") << "Volume [" << i << "] " << names[i]
                  << " of thickness " << thick[i]
                  << " filled with " << materials[i]
                  << " first copy number " << copyNumber[i];
  edm::LogInfo("HGCalGeom") << "DDHGCalEEAlgo: " << layerType.size()
                << " layers";
  for (unsigned int i=0; i<layerType.size(); ++i)
    edm::LogInfo("HGCalGeom") << "Layer [" << i << "] with material type " 
                  << layerType[i] << " height type "
                  << heightType[i] << " block thickness "
                  << thickBlock[i];

  sectors       = (int)(nArgs["Sectors"]);
  slopeB        = nArgs["SlopeBottom"];
  slopeT        = vArgs["SlopeTop"];
  zFront        = vArgs["ZFront"];
  rMaxFront     = vArgs["RMaxFront"];
  idName        = parent().name().name();
  idNameSpace   = DDCurrentNamespace::ns();
  edm::LogInfo("HGCalGeom") << "DDHGCalEEAlgo: Bottom slope " << slopeB
                << " " << slopeT.size() << " slopes for top";
  for (unsigned int i=0; i<slopeT.size(); ++i)
    edm::LogInfo("HGCalGeom") << "Block [" << i << "] Zmin " << zFront[i]
                  << " Rmax " << rMaxFront[i] << " Slope " 
                  << slopeT[i];
  edm::LogInfo("HGCalGeom") << "DDHGCalEEAlgo: Sectors " << sectors
                << "\tNameSpace:Name " << idNameSpace
                << ":" << idName;

}

DDHGCalEEAlgo::HGCalEEPar DDHGCalEEAlgo::parameterLayer ( double  rinF, double  routF, double  rinB, double  routB, double  zi, double  zo  )

protected

Definition at line 151 of file DDHGCalEEAlgo.cc.

References DDHGCalEEAlgo::HGCalEEPar::alp, alpha, DDHGCalEEAlgo::HGCalEEPar::bl1, DDHGCalEEAlgo::HGCalEEPar::bl2, funct::cos(), DDHGCalEEAlgo::HGCalEEPar::phi, pi, sectors, funct::sin(), funct::tan(), DDHGCalEEAlgo::HGCalEEPar::theta, DDHGCalEEAlgo::HGCalEEPar::tl1, DDHGCalEEAlgo::HGCalEEPar::tl2, DDHGCalEEAlgo::HGCalEEPar::xpos, DDHGCalEEAlgo::HGCalEEPar::yh1, DDHGCalEEAlgo::HGCalEEPar::yh2, DDHGCalEEAlgo::HGCalEEPar::ypos, and DDHGCalEEAlgo::HGCalEEPar::zpos.

Referenced by constructLayers(), and DDHGCalEEAlgo::HGCalEEPar::HGCalEEPar().

                                                      {

  DDHGCalEEAlgo::HGCalEEPar 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("HGCalGeom") << "Input: Front " << rinF << " " << routF << " "
                << zi << " Back " << rinB << " " << routB << " "
                << zo << " Alpha " << alpha/CLHEP::deg;

  parm.yh2  = parm.yh1  = 0.5 * (routF*cos(alpha) - rinB);
  parm.bl2  = parm.bl1  = rinB  * tan(alpha);
  parm.tl2  = parm.tl1  = routF * sin(alpha);
  parm.xpos = 0.5*(routF*cos(alpha)+rinB);
  parm.ypos = 0.0;
  parm.zpos = 0.5*(zi+zo);
  parm.alp  = parm.theta  = parm.phi = 0;
  edm::LogInfo("HGCalGeom") << "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;
}

double DDHGCalEEAlgo::rMax ( double  z )

protected

Definition at line 178 of file DDHGCalEEAlgo.cc.

References gen::k, alignCSCRings::r, rMaxFront, slopeT, and zFront.

Referenced by constructLayers(), and DDHGCalEEAlgo::HGCalEEPar::HGCalEEPar().

                                   {

  double r(0);
  unsigned int ik(0);
  for (unsigned int k=0; k<slopeT.size(); ++k) {
    if (z < zFront[k]) break;
    r  = rMaxFront[k] + (z - zFront[k]) * slopeT[k];
    ik = k;
  }
  edm::LogInfo("HGCalGeom") << "rMax : " << z << ":" << ik << ":" << r ;
  return r;
}

Member Data Documentation

std::vector<int> DDHGCalEEAlgo::copyNumber

private

Definition at line 46 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<int> DDHGCalEEAlgo::heightType

private

Definition at line 45 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().

std::string DDHGCalEEAlgo::idName

private

Definition at line 55 of file DDHGCalEEAlgo.h.

Referenced by initialize().

std::string DDHGCalEEAlgo::idNameSpace

private

Definition at line 56 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<int> DDHGCalEEAlgo::layerType

private

Definition at line 44 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<std::string> DDHGCalEEAlgo::materials

private

Definition at line 41 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<std::string> DDHGCalEEAlgo::names

private

Definition at line 42 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<double> DDHGCalEEAlgo::rMaxFront

private

Definition at line 54 of file DDHGCalEEAlgo.h.

Referenced by initialize(), and rMax().

std::string DDHGCalEEAlgo::rotstr

private

Definition at line 43 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().

int DDHGCalEEAlgo::sectors

private

Definition at line 50 of file DDHGCalEEAlgo.h.

Referenced by initialize(), and parameterLayer().

double DDHGCalEEAlgo::slopeB

private

Definition at line 51 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<double> DDHGCalEEAlgo::slopeT

private

Definition at line 52 of file DDHGCalEEAlgo.h.

Referenced by initialize(), and rMax().

std::vector<double> DDHGCalEEAlgo::thick

private

Definition at line 47 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<double> DDHGCalEEAlgo::thickBlock

private

Definition at line 48 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<double> DDHGCalEEAlgo::zFront

private

Definition at line 53 of file DDHGCalEEAlgo.h.

Referenced by initialize(), and rMax().

double DDHGCalEEAlgo::zMinBlock

private

Definition at line 49 of file DDHGCalEEAlgo.h.

Referenced by constructLayers(), and initialize().