DDHGCalHEAlgo Class Reference

#include <DDHGCalHEAlgo.h>

Inheritance diagram for DDHGCalHEAlgo:

Classes
struct	HGCalHEPar

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

Protected Member Functions
void	constructLayers (DDLogicalPart, DDCompactView &cpv)
HGCalHEPar	parameterLayer (double rinF, double routF, double rinB, double routB, double zi, double zo)
HGCalHEPar	parameterLayer (int type, 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
double	thickModule
std::vector< int >	type
std::vector< double >	zFront
std::vector< double >	zMinBlock

Detailed Description

Definition at line 10 of file DDHGCalHEAlgo.h.

Constructor & Destructor Documentation

DDHGCalHEAlgo::DDHGCalHEAlgo

(

)

Definition at line 20 of file DDHGCalHEAlgo.cc.

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

DDHGCalHEAlgo::~DDHGCalHEAlgo ( )

virtual

Definition at line 24 of file DDHGCalHEAlgo.cc.

{}

Member Function Documentation

void DDHGCalHEAlgo::constructLayers ( DDLogicalPart  module, DDCompactView &  cpv  )

protected

Definition at line 87 of file DDHGCalHEAlgo.cc.

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

Referenced by execute().

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

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

  double zz(zMinBlock[0]);
  for (unsigned int i=0; i<layerType.size(); i++) {
    int   ii       = layerType[i];
    int   copy     = copyNumber[ii];
    ++copyNumber[ii];
    int   ityp     = type[ii];
    type[ii]       =-ityp;
    double zi      = zMinBlock[ii];
    zMinBlock[ii] += thickModule;
    if (heightType[i] == 0) zz = zi;

    double zlayer  = zz + 2*thickModule;
    if((i % 6)>2) zlayer  = zz + thickModule;

    double zo      = zi + thick[ii];
    double rinF    = zi * slopeB;
    double rinB    = zlayer * slopeB;

    double routF   = (heightType[i] == 0) ? rMax(zi) : rMax(zz);
    if((i % 6)>2) routF   = (heightType[i] == 0) ? rMax(zi-thickModule) : rMax(zz-thickModule);
    
    double routB   = rMax(zo);
    std::string name = "HGCal"+names[ii]+std::to_string(copy);
    edm::LogInfo("HGCalGeom") << "DDHGCalHEAlgo test: Layer " << i << ":" 
                  << ii << ":" << ityp << " Front " << zi << ", " 
                  << rinF << ", " << routF << " Back " << zo 
                  << ", " << rinB << ", " << routB;
    DDHGCalHEAlgo::HGCalHEPar parm = (ityp == 0) ?
      parameterLayer(rinF, routF, rinB, routB, zi, zo) :
      parameterLayer(ityp, 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") << "DDHGCalHEAlgo 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") << "DDHGCalHEAlgo test: " << glog.name()
                  << " number " << copy << " positioned in " 
                  << module.name() << " at " << r1 << " with " 
                  << rot;
  }   // End of loop on layers
}

void DDHGCalHEAlgo::execute

(

DDCompactView &

cpv

)

Definition at line 80 of file DDHGCalHEAlgo.cc.

References constructLayers(), and SpecificationBuilder_cfi::parent().

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

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

Definition at line 26 of file DDHGCalHEAlgo.cc.

References copyNumber, dbl_to_int(), heightType, i, idName, idNameSpace, layerType, materials, names, DDCurrentNamespace::ns(), SpecificationBuilder_cfi::parent(), rMaxFront, rotstr, sectors, slopeB, slopeT, thick, thickModule, zFront, and zMinBlock.

                                                         {

  materials     = vsArgs["MaterialNames"];
  names         = vsArgs["VolumeNames"];
  thick         = vArgs["Thickness"];
  type          = dbl_to_int(vArgs["Type"]);
  copyNumber    = dbl_to_int(vArgs["Offsets"]);
  zMinBlock     = vArgs["ZMinType"];
  rotstr        = sArgs["Rotation"];
  layerType     = dbl_to_int(vArgs["LayerType"]);
  heightType    = dbl_to_int(vArgs["HeightType"]);
  thickModule   = nArgs["ThickModule"];
  edm::LogInfo("HGCalGeom") << "DDHGCalHEAlgo: " << materials.size()
                << " volumes to be put with rotation " << rotstr
                << " in " << layerType.size() << " layers and "
                << "module thickness " << thickModule;
  for (unsigned int i=0; i<names.size(); ++i)
    edm::LogInfo("HGCalGeom") << "Volume [" << i << "] " << names[i]
                  << " filled with " << materials[i] << " of type "
                  << type[i] << " thickness " << thick[i]
                  << " first copy number " << copyNumber[i]
                  << " starting at " << zMinBlock[i];
  for (unsigned int i=0; i<layerType.size(); ++i)
    edm::LogInfo("HGCalGeom") << "Layer [" << i << "] with material type " 
                  << layerType[i] << " height type "
                  << heightType[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") << "DDHGCalHEAlgo: 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") << "DDHGCalHEAlgo: Sectors " << sectors
                << "\tNameSpace:Name " << idNameSpace
                << ":" << idName;

}

DDHGCalHEAlgo::HGCalHEPar DDHGCalHEAlgo::parameterLayer ( double  rinF, double  routF, double  rinB, double  routB, double  zi, double  zo  )

protected

Definition at line 155 of file DDHGCalHEAlgo.cc.

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

Referenced by constructLayers().

                                                      {

  DDHGCalHEAlgo::HGCalHEPar parm;
  //Given rin, rout compute parameters of the trapezoid and 
  //position of the trapezoid for a standrd layer
  double alpha = CLHEP::pi/sectors;
  double rout  = routF;
  edm::LogInfo("HGCalGeom") << "Input: Front " << rinF << " " << routF << " "
                << zi << " Back " << rinB << " " << routB << " "
                << zo << " Alpha " << alpha/CLHEP::deg << " Rout "
                << rout;

  parm.yh2  = parm.yh1  = 0.5 * (rout*cos(alpha) - rinB);
  parm.bl2  = parm.bl1  = rinB  * tan(alpha);
  parm.tl2  = parm.tl1  = rout  * sin(alpha);
  parm.xpos = 0.5*(rout*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;
}

DDHGCalHEAlgo::HGCalHEPar DDHGCalHEAlgo::parameterLayer ( int  type, double  rinF, double  routF, double  rinB, double  routB, double  zi, double  zo  )

protected

Definition at line 185 of file DDHGCalHEAlgo.cc.

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

                                                      {

  DDHGCalHEAlgo::HGCalHEPar parm;
  //Given rin, rout compute parameters of the trapezoid and 
  //position of the trapezoid for a standrd layer
  double alpha = CLHEP::pi/sectors;
  double rout  = routF;
  edm::LogInfo("HGCalGeom") << "Input " << type << " Front " << rinF << " " 
                << routF << " " << zi << " Back " << rinB << " " 
                << routB << " " << zo <<" Alpha " 
                << alpha/CLHEP::deg << " rout " << rout;

  parm.yh2  = parm.yh1  = 0.5 * (rout*cos(alpha) - rinB);
  parm.bl2  = parm.bl1  = 0.5 * rinB  * tan(alpha);
  parm.tl2  = parm.tl1  = 0.5 * rout  * sin(alpha);
  double dx = 0.25* (parm.bl2+parm.tl2-parm.bl1-parm.tl1);
  double dy = 0.0;
  parm.xpos = 0.5*(rout*cos(alpha)+rinB);
  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 (type > 0) {
    parm.ypos = -parm.ypos;
  } else {
    parm.alp  = -parm.alp;
    dx        = -dx;
  }
  double r    = sqrt (dx*dx + dy*dy);
  edm::LogInfo("HGCalGeom") << "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("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 DDHGCalHEAlgo::rMax ( double  z )

protected

Definition at line 230 of file DDHGCalHEAlgo.cc.

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

Referenced by constructLayers().

                                   {

  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> DDHGCalHEAlgo::copyNumber

private

Definition at line 46 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<int> DDHGCalHEAlgo::heightType

private

Definition at line 51 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().

std::string DDHGCalHEAlgo::idName

private

Definition at line 59 of file DDHGCalHEAlgo.h.

Referenced by initialize().

std::string DDHGCalHEAlgo::idNameSpace

private

Definition at line 60 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<int> DDHGCalHEAlgo::layerType

private

Definition at line 50 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<std::string> DDHGCalHEAlgo::materials

private

Definition at line 43 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<std::string> DDHGCalHEAlgo::names

private

Definition at line 44 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<double> DDHGCalHEAlgo::rMaxFront

private

Definition at line 57 of file DDHGCalHEAlgo.h.

Referenced by initialize(), and rMax().

std::string DDHGCalHEAlgo::rotstr

private

Definition at line 45 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().

int DDHGCalHEAlgo::sectors

private

Definition at line 53 of file DDHGCalHEAlgo.h.

Referenced by initialize(), and parameterLayer().

double DDHGCalHEAlgo::slopeB

private

Definition at line 54 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<double> DDHGCalHEAlgo::slopeT

private

Definition at line 55 of file DDHGCalHEAlgo.h.

Referenced by initialize(), and rMax().

std::vector<double> DDHGCalHEAlgo::thick

private

Definition at line 47 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().

double DDHGCalHEAlgo::thickModule

private

Definition at line 52 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().

std::vector<int> DDHGCalHEAlgo::type

private

Definition at line 48 of file DDHGCalHEAlgo.h.

Referenced by cuy.ValElement::__init__(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputCommands(), core.autovars.NTupleVariable::makeBranch(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputCommands(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and core.AutoHandle.AutoHandle::ReallyLoad().

std::vector<double> DDHGCalHEAlgo::zFront

private

Definition at line 56 of file DDHGCalHEAlgo.h.

Referenced by initialize(), and rMax().

std::vector<double> DDHGCalHEAlgo::zMinBlock

private

Definition at line 49 of file DDHGCalHEAlgo.h.

Referenced by constructLayers(), and initialize().