#include <DDTIBRadCableAlgo_MTCC.h>

Inheritance diagram for DDTIBRadCableAlgo_MTCC:

Public Member Functions
	DDTIBRadCableAlgo_MTCC ()

void	execute (DDCompactView &cpv) override

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

	~DDTIBRadCableAlgo_MTCC () override

Private Attributes
std::vector< std::string >	cableMat

double	cableT

double	cylinderT

double	deltaR

std::string	idNameSpace

std::vector< double >	layers

std::vector< double >	layRin

double	rMax

double	rMin

std::vector< std::string >	strucMat

double	supportDR

std::string	supportMat

double	supportT

Detailed Description

Definition at line 10 of file DDTIBRadCableAlgo_MTCC.h.

Constructor & Destructor Documentation

◆ DDTIBRadCableAlgo_MTCC()

DDTIBRadCableAlgo_MTCC::DDTIBRadCableAlgo_MTCC ( )

Definition at line 19 of file DDTIBRadCableAlgo_MTCC.cc.

                                                : layRin(0), cableMat(0), strucMat(0), layers(0) {
   LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC info: Creating an instance";
 }

References LogDebug.

◆ ~DDTIBRadCableAlgo_MTCC()

DDTIBRadCableAlgo_MTCC::~DDTIBRadCableAlgo_MTCC ( )

override

Definition at line 23 of file DDTIBRadCableAlgo_MTCC.cc.

23 {}

Member Function Documentation

◆ execute()

void DDTIBRadCableAlgo_MTCC::execute ( DDCompactView & cpv )

override

Definition at line 70 of file DDTIBRadCableAlgo_MTCC.cc.

                                                        {
   LogDebug("TIBGeom") << "==>> Constructing DDTIBRadCableAlgo_MTCC...";
   DDName diskName = parent().name();
  
   DDSolid solid;
   std::string name;
   double rin, rout;
  
   // Loop over sub disks
   DDName suppName(DDSplit(supportMat).first, DDSplit(supportMat).second);
   DDMaterial suppMatter(suppName);
   double diskDz = 0.5 * (supportT + cableT * layRin.size());
   double dz = 0.5 * supportT;
  
   for (unsigned int i = 0; i < layRin.size(); i++) {
     // fill only layers in layers list
     bool empty = true;
     for (double layer : layers) {
       if (i + 1 == (unsigned int)layer) {
         empty = false;
       }
     }
  
     if (!empty) {
       //Support disks
       name = "TIBSupportSideDisk" + std::to_string(i);
       rin = layRin[i] + 0.5 * (deltaR - cylinderT) - supportDR;
       rout = layRin[i] + 0.5 * (deltaR + cylinderT) + supportDR;
       solid = DDSolidFactory::tubs(DDName(name, idNameSpace), dz, rin, rout, 0, CLHEP::twopi);
       LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC test: " << DDName(name, idNameSpace) << " Tubs made of "
                           << supportMat << " from 0 to " << CLHEP::twopi / CLHEP::deg << " with Rin " << rin << " Rout "
                           << rout << " ZHalf " << dz;
       DDLogicalPart suppLogic(DDName(name, idNameSpace), suppMatter, solid);
  
       DDTranslation r1(0, 0, (dz - diskDz));
       cpv.position(DDName(name, idNameSpace), diskName, i + 1, r1, DDRotation());
       LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC test: " << DDName(name, idNameSpace) << " number " << i + 1
                           << " positioned in " << diskName << " at " << r1 << " with no rotation";
  
       //Open Structure
       name = "TIBOpenZone" + std::to_string(i);
       rout = rin;
       if (i == 0)
         rin = rMin;
       else
         rin = layRin[i - 1] + 0.5 * (deltaR + cylinderT) + supportDR;
       solid = DDSolidFactory::tubs(DDName(name, idNameSpace), dz, rin, rout, 0, CLHEP::twopi);
       LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC test: " << DDName(name, idNameSpace) << " Tubs made of "
                           << strucMat[i] << " from 0 to " << CLHEP::twopi / CLHEP::deg << " with Rin " << rin
                           << " Rout " << rout << " ZHalf " << dz;
       DDName strucName(DDSplit(strucMat[i]).first, DDSplit(strucMat[i]).second);
       DDMaterial strucMatter(strucName);
       DDLogicalPart strucLogic(DDName(name, idNameSpace), strucMatter, solid);
  
       DDTranslation r2(0, 0, (dz - diskDz));
       cpv.position(DDName(name, idNameSpace), diskName, i + 1, r2, DDRotation());
       LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC test: " << DDName(name, idNameSpace) << " number " << i + 1
                           << " positioned in " << diskName << " at " << r2 << " with no rotation";
  
       //Now the radial cable
       name = "TIBRadCable" + std::to_string(i);
       double rv = layRin[i] + 0.5 * deltaR;
       std::vector<double> pgonZ;
       pgonZ.emplace_back(-0.5 * cableT);
       pgonZ.emplace_back(cableT * (rv / rMax - 0.5));
       pgonZ.emplace_back(0.5 * cableT);
       std::vector<double> pgonRmin;
       pgonRmin.emplace_back(rv);
       pgonRmin.emplace_back(rv);
       pgonRmin.emplace_back(rv);
       std::vector<double> pgonRmax;
       pgonRmax.emplace_back(rMax);
       pgonRmax.emplace_back(rMax);
       pgonRmax.emplace_back(rv);
       solid = DDSolidFactory::polycone(DDName(name, idNameSpace), 0, CLHEP::twopi, pgonZ, pgonRmin, pgonRmax);
       LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC test: " << DDName(name, idNameSpace) << " Polycone made of "
                           << cableMat[i] << " from 0 to " << CLHEP::twopi / CLHEP::deg << " and with " << pgonZ.size()
                           << " sections";
       for (unsigned int ii = 0; ii < pgonZ.size(); ii++)
         LogDebug("TIBGeom") << "\t"
                             << "\tZ = " << pgonZ[ii] << "\tRmin = " << pgonRmin[ii] << "\tRmax = " << pgonRmax[ii];
       DDName cableName(DDSplit(cableMat[i]).first, DDSplit(cableMat[i]).second);
       LogDebug("TIBGeom") << " material cableName " << i << " " << cableName;
       DDMaterial cableMatter(cableName);
       DDLogicalPart cableLogic(DDName(name, idNameSpace), cableMatter, solid);
  
       DDTranslation r3(0, 0, (diskDz - (i + 0.5) * cableT));
       cpv.position(DDName(name, idNameSpace), diskName, i + 1, r3, DDRotation());
       LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC test: " << DDName(name, idNameSpace) << " number " << i + 1
                           << " positioned in " << diskName << " at " << r3 << " with no rotation";
     }
     //
  
     // fill only layers in layers list
     empty = true;
     for (double layer : layers) {
       if (i + 1 == (unsigned int)layer) {
         empty = false;
       }
     }
  
     if (!empty) {
       //Now the last open zone
       unsigned int i = layRin.size();
       rin = layRin[i - 1] + 0.5 * (deltaR + cylinderT) + supportDR;
       rout = rMax;
       name = "TIBOpenZone" + std::to_string(i);
       solid = DDSolidFactory::tubs(DDName(name, idNameSpace), dz, rin, rout, 0, CLHEP::twopi);
       LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC test: " << DDName(name, idNameSpace) << " Tubs made of "
                           << strucMat[i] << " from 0 to " << CLHEP::twopi / CLHEP::deg << " with Rin " << rin
                           << " Rout " << rout << " ZHalf " << dz;
       DDName strucName(DDSplit(strucMat[i]).first, DDSplit(strucMat[i]).second);
       DDMaterial strucMatter(strucName);
       DDLogicalPart strucLogic(DDName(name, idNameSpace), strucMatter, solid);
  
       DDTranslation r2(0, 0, (dz - diskDz));
       cpv.position(DDName(name, idNameSpace), diskName, i + 1, r2, DDRotation());
       LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC test: " << DDName(name, idNameSpace) << " number " << i + 1
                           << " positioned in " << diskName << " at " << r2 << " with no rotation";
     }
   }
   LogDebug("TIBGeom") << "<<== End of DDTIBRadCableAlgo_MTCC construction ...";
 }

References cableMat, cableT, cylinderT, DDSplit(), deltaR, PVValHelper::dz, relativeConstraints::empty, dqmdumpme::first, mps_fire::i, idNameSpace, cuy::ii, layers, layRin, LogDebug, Skims_PA_cff::name, class-composition::parent, DDSolidFactory::polycone(), DDCompactView::position(), diffTwoXMLs::r1, diffTwoXMLs::r2, rMax, rMin, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, strucMat, supportDR, supportMat, supportT, and DDSolidFactory::tubs().

◆ initialize()

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

override

Definition at line 25 of file DDTIBRadCableAlgo_MTCC.cc.

                                                                                {
   idNameSpace = DDCurrentNamespace::ns();
   unsigned int i;
   DDName parentName = parent().name();
   LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC debug: Parent " << parentName << " NameSpace " << idNameSpace;
  
   rMin = nArgs["RMin"];
   rMax = nArgs["RMax"];
   layRin = vArgs["RadiusLo"];
   deltaR = nArgs["DeltaR"];
   LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC debug: Disk Rmin " << rMin << "\tRMax " << rMax
                       << "\tSeparation of layers " << deltaR << " with " << layRin.size() << " layers at R =";
   for (i = 0; i < layRin.size(); i++)
     LogDebug("TIBGeom") << " " << i << " : " << layRin[i];
  
   cylinderT = nArgs["CylinderThick"];
   supportT = nArgs["SupportThick"];
   supportDR = nArgs["SupportDR"];
   supportMat = sArgs["SupportMaterial"];
   LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC debug: SupportCylinder "
                       << "Thickness " << cylinderT << "\tSupportDisk "
                       << "Thickness " << supportT << "\tExtra width along R " << supportDR
                       << "\tMaterial: " << supportMat;
  
   cableT = nArgs["CableThick"];
   cableMat = vsArgs["CableMaterial"];
   LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC debug: Cable Thickness " << cableT << " with materials: ";
   for (i = 0; i < cableMat.size(); i++)
     LogDebug("TIBGeom") << " " << i << " : " << cableMat[i];
  
   strucMat = vsArgs["StructureMaterial"];
   LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC debug: " << strucMat.size() << " materials for open structure:";
   for (i = 0; i < strucMat.size(); i++)
     LogDebug("TIBGeom") << " " << i << " : " << strucMat[i];
  
   layers = vArgs["Layers"];
   LogDebug("TIBGeom") << "DDTIBRadCableAlgo_MTCC debug: " << layers.size() << " layers:";
   for (i = 0; i < layers.size(); i++)
     LogDebug("TIBGeom") << " " << i << " : " << layers[i];
 }