DDTIDAxialCableAlgo Class Reference

Inheritance diagram for DDTIDAxialCableAlgo:

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

Private Attributes
vector< double >	angles
string	childName
string	idNameSpace
string	matIn
string	matOut
double	rMax
double	rMin
double	rTop
double	thick
double	width
double	zBend
double	zEnd
vector< double >	zposRing
vector< double >	zposWheel

Detailed Description

Definition at line 23 of file DDTIDAxialCableAlgo.cc.

Constructor & Destructor Documentation

DDTIDAxialCableAlgo()

DDTIDAxialCableAlgo::DDTIDAxialCableAlgo

(

)

Definition at line 55 of file DDTIDAxialCableAlgo.cc.

References LogDebug.

{ LogDebug("TIDGeom") << "DDTIDAxialCableAlgo info: Creating an instance"; }

~DDTIDAxialCableAlgo()

DDTIDAxialCableAlgo::~DDTIDAxialCableAlgo ( )

override

Definition at line 57 of file DDTIDAxialCableAlgo.cc.

{}

Member Function Documentation

execute()

void DDTIDAxialCableAlgo::execute ( DDCompactView &  cpv )

override

Definition at line 100 of file DDTIDAxialCableAlgo.cc.

References particleFlowDisplacedVertex_cfi::angles, DDrot(), DDSplit(), dqmdumpme::first, mps_fire::i, createfilelist::int, dqmdumpme::k, LogDebug, Skims_PA_cff::name, class-composition::parent, DDSolidFactory::polycone(), DDCompactView::position(), photonAnalyzer_cfi::rMax, photonAnalyzer_cfi::rMin, idealTransformation::rotation, edm::second(), theta(), to_string(), DDSolidFactory::tubs(), and ApeEstimator_cff::width.

                                                    {
  DDName mother = parent().name();
  vector<DDName> logs;
  double thk = thick / zposRing.size();
  double r = rMin;
  double thktot = 0;
  double z;

  //Cables between the wheels
  for (int k = 0; k < (int)(zposWheel.size()); k++) {
    vector<double> pconZ, pconRmin, pconRmax;
    for (int i = 0; i < (int)(zposRing.size()); i++) {
      thktot += thk;
      z = zposWheel[k] + zposRing[i] - 0.5 * thk;
      if (i != 0) {
        pconZ.emplace_back(z);
        pconRmin.emplace_back(r);
        pconRmax.emplace_back(rMax);
      }
      r = rMin;
      pconZ.emplace_back(z);
      pconRmin.emplace_back(r);
      pconRmax.emplace_back(rMax);
      z += thk;
      pconZ.emplace_back(z);
      pconRmin.emplace_back(r);
      pconRmax.emplace_back(rMax);
      r = rMax - thktot;
      pconZ.emplace_back(z);
      pconRmin.emplace_back(r);
      pconRmax.emplace_back(rMax);
    }
    if (k >= ((int)(zposWheel.size()) - 1))
      z = zBend;
    else
      z = zposWheel[k + 1] + zposRing[0] - 0.5 * thk;
    pconZ.emplace_back(z);
    pconRmin.emplace_back(r);
    pconRmax.emplace_back(rMax);

    string name = childName + to_string(k);
    DDSolid solid = DDSolidFactory::polycone(DDName(name, idNameSpace), -0.5 * width, width, pconZ, pconRmin, pconRmax);

    LogDebug("TIDGeom") << "DDTIDAxialCableAlgo test: " << DDName(name, idNameSpace) << " Polycone made of " << matIn
                        << " from " << -0.5 * width / CLHEP::deg << " to " << 0.5 * width / CLHEP::deg << " and with "
                        << pconZ.size() << " sections ";
    for (int i = 0; i < (int)(pconZ.size()); i++)
      LogDebug("TIDGeom") << "\t[" << i << "]\tZ = " << pconZ[i] << "\tRmin = " << pconRmin[i]
                          << "\tRmax = " << pconRmax[i];

    DDName mat(DDSplit(matIn).first, DDSplit(matIn).second);
    DDMaterial matter(mat);
    DDLogicalPart genlogic(DDName(name, idNameSpace), matter, solid);
    logs.emplace_back(DDName(name, idNameSpace));
  }

  //Cable in the vertical part
  vector<double> pconZ, pconRmin, pconRmax;
  r = thktot * rMax / rTop;
  z = zBend - thktot;
  LogDebug("TIDGeom") << "DDTIDAxialCableAlgo test: Thk " << thk << " Total " << thktot << " rMax " << rMax << " rTop "
                      << rTop << " dR " << r << " z " << z;
  pconZ.emplace_back(z);
  pconRmin.emplace_back(rMax);
  pconRmax.emplace_back(rMax);
  z = zBend - r;
  pconZ.emplace_back(z);
  pconRmin.emplace_back(rMax);
  pconRmax.emplace_back(rTop);
  pconZ.emplace_back(zBend);
  pconRmin.emplace_back(rMax);
  pconRmax.emplace_back(rTop);

  string name = childName + to_string(zposWheel.size());
  DDSolid solid = DDSolidFactory::polycone(DDName(name, idNameSpace), -0.5 * width, width, pconZ, pconRmin, pconRmax);

  LogDebug("TIDGeom") << "DDTIDAxialCableAlgo test: " << DDName(name, idNameSpace) << " Polycone made of " << matIn
                      << " from " << -0.5 * width / CLHEP::deg << " to " << 0.5 * width / CLHEP::deg << " and with "
                      << pconZ.size() << " sections ";
  for (int i = 0; i < (int)(pconZ.size()); i++)
    LogDebug("TIDGeom") << "\t[" << i << "]\tZ = " << pconZ[i] << "\tRmin = " << pconRmin[i]
                        << "\tRmax = " << pconRmax[i];

  DDMaterial matter(DDName(DDSplit(matIn).first, DDSplit(matIn).second));
  DDLogicalPart genlogic(DDName(name, idNameSpace), matter, solid);
  logs.emplace_back(DDName(name, idNameSpace));

  //Cable in the outer part
  name = childName + to_string(zposWheel.size() + 1);
  r = rTop - r;
  solid = DDSolidFactory::tubs(DDName(name, idNameSpace), 0.5 * (zEnd - zBend), r, rTop, -0.5 * width, width);
  LogDebug("TIDGeom") << "DDTIDAxialCableAlgo test: " << DDName(name, idNameSpace) << " Tubs made of " << matOut
                      << " from " << -0.5 * width / CLHEP::deg << " to " << 0.5 * width / CLHEP::deg << " with Rin "
                      << r << " Rout " << rTop << " ZHalf " << 0.5 * (zEnd - zBend);
  matter = DDMaterial(DDName(DDSplit(matOut).first, DDSplit(matOut).second));
  genlogic = DDLogicalPart(DDName(name, idNameSpace), matter, solid);
  logs.emplace_back(DDName(name, idNameSpace));

  //Position the cables
  double theta = 90. * CLHEP::deg;
  for (int i = 0; i < (int)(angles.size()); i++) {
    double phix = angles[i];
    double phiy = phix + 90. * CLHEP::deg;
    double phideg = phix / CLHEP::deg;

    DDRotation rotation;
    if (phideg != 0) {
      string rotstr = childName + to_string(phideg * 10.);
      rotation = DDRotation(DDName(rotstr, idNameSpace));
      if (!rotation) {
        LogDebug("TIDGeom") << "DDTIDAxialCableAlgo test: Creating a new "
                            << "rotation: " << rotstr << " " << theta / CLHEP::deg << ", " << phix / CLHEP::deg << ", "
                            << theta / CLHEP::deg << ", " << phiy / CLHEP::deg << ", 0, 0";
        rotation = DDrot(DDName(rotstr, idNameSpace), theta, phix, theta, phiy, 0., 0.);
      }
    }

    for (int k = 0; k < (int)(logs.size()); k++) {
      DDTranslation tran(0, 0, 0);
      if (k == ((int)(logs.size()) - 1))
        tran = DDTranslation(0, 0, 0.5 * (zEnd + zBend));
      cpv.position(logs[k], mother, i + 1, tran, rotation);
      LogDebug("TIDGeom") << "DDTIDAxialCableAlgo test " << logs[k] << " number " << i + 1 << " positioned in "
                          << mother << " at " << tran << " with " << rotation;
    }
  }
}

initialize()

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

override

Definition at line 59 of file DDTIDAxialCableAlgo.cc.

References particleFlowDisplacedVertex_cfi::angles, mps_fire::i, createfilelist::int, LogDebug, DDCurrentNamespace::ns(), class-composition::parent, photonAnalyzer_cfi::rMax, photonAnalyzer_cfi::rMin, and ApeEstimator_cff::width.

                                                                     {
  zBend = nArgs["ZBend"];
  zEnd = nArgs["ZEnd"];
  rMin = nArgs["RMin"];
  rMax = nArgs["RMax"];
  rTop = nArgs["RTop"];
  width = nArgs["Width"];
  thick = nArgs["Thick"];
  angles = vArgs["Angles"];
  zposWheel = vArgs["ZPosWheel"];
  zposRing = vArgs["ZPosRing"];

  LogDebug("TIDGeom") << "DDTIDAxialCableAlgo debug: Parameters for creating " << (zposWheel.size() + 2)
                      << " axial cables and position"
                      << "ing " << angles.size() << " copies in Service volume"
                      << "\n                            zBend " << zBend << " zEnd " << zEnd << " rMin " << rMin
                      << " rMax " << rMax << " Cable width " << width / CLHEP::deg << " thickness " << thick
                      << " with Angles";
  for (int i = 0; i < (int)(angles.size()); i++)
    LogDebug("TIDGeom") << "\tangles[" << i << "] = " << angles[i] / CLHEP::deg;
  LogDebug("TIDGeom") << "                          Wheels " << zposWheel.size() << " at Z";
  for (int i = 0; i < (int)(zposWheel.size()); i++)
    LogDebug("TIDGeom") << "\tzposWheel[" << i << "] = " << zposWheel[i];
  LogDebug("TIDGeom") << "                          each with " << zposRing.size() << " Rings at Z";
  for (int i = 0; i < (int)(zposRing.size()); i++)
    LogDebug("TIDGeom") << "\tzposRing[" << i << "] = " << zposRing[i];

  idNameSpace = DDCurrentNamespace::ns();
  childName = sArgs["ChildName"];
  matIn = sArgs["MaterialIn"];
  matOut = sArgs["MaterialOut"];

  DDName parentName = parent().name();
  LogDebug("TIDGeom") << "DDTIDAxialCableAlgo debug: Parent " << parentName << "\tChild " << childName << " NameSpace "
                      << idNameSpace << "\tMaterial " << matIn << " and " << matOut;
}

Member Data Documentation

angles

vector<double> DDTIDAxialCableAlgo::angles

private

Definition at line 45 of file DDTIDAxialCableAlgo.cc.

childName

string DDTIDAxialCableAlgo::childName

private

Definition at line 50 of file DDTIDAxialCableAlgo.cc.

idNameSpace

string DDTIDAxialCableAlgo::idNameSpace

private

Definition at line 49 of file DDTIDAxialCableAlgo.cc.

matIn

string DDTIDAxialCableAlgo::matIn

private

Definition at line 51 of file DDTIDAxialCableAlgo.cc.

matOut

string DDTIDAxialCableAlgo::matOut

private

Definition at line 52 of file DDTIDAxialCableAlgo.cc.

rMax

double DDTIDAxialCableAlgo::rMax

private

Definition at line 41 of file DDTIDAxialCableAlgo.cc.

rMin

double DDTIDAxialCableAlgo::rMin

private

Definition at line 40 of file DDTIDAxialCableAlgo.cc.

rTop

double DDTIDAxialCableAlgo::rTop

private

Definition at line 42 of file DDTIDAxialCableAlgo.cc.

thick

double DDTIDAxialCableAlgo::thick

private

Definition at line 44 of file DDTIDAxialCableAlgo.cc.

width

double DDTIDAxialCableAlgo::width

private

Definition at line 43 of file DDTIDAxialCableAlgo.cc.

Referenced by python.rootplot.utilities.Hist::__init_TGraph(), python.rootplot.root2matplotlib.Hist::bar(), python.rootplot.root2matplotlib.Hist::barh(), and python.rootplot.utilities.Hist::delete_bin().

zBend

double DDTIDAxialCableAlgo::zBend

private

Definition at line 38 of file DDTIDAxialCableAlgo.cc.

zEnd

double DDTIDAxialCableAlgo::zEnd

private

Definition at line 39 of file DDTIDAxialCableAlgo.cc.

zposRing

vector<double> DDTIDAxialCableAlgo::zposRing

private

Definition at line 47 of file DDTIDAxialCableAlgo.cc.

zposWheel

vector<double> DDTIDAxialCableAlgo::zposWheel

private

Definition at line 46 of file DDTIDAxialCableAlgo.cc.