Inheritance diagram for DDPixBarLayerUpgradeAlgo:

Public Member Functions
	DDPixBarLayerUpgradeAlgo ()

void	execute (DDCompactView &cpv) override

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

	~DDPixBarLayerUpgradeAlgo () override

Private Attributes
double	cool1Offset

double	cool2Offset

double	coolDist

double	coolDz

string	coolMat

string	coolMatHalf

double	coolRadius

double	coolThick

string	genMat

string	idNameSpace

string	ladder

double	ladderOffset

double	ladderThick

double	ladderWidth

double	layerDz

int	number

int	outerFirst

double	phiFineTune

double	rInnerFineTune

double	rOuterFineTune

string	tubeMat

string	tubeMatHalf

Detailed Description

Definition at line 26 of file DDPixBarLayerUpgradeAlgo.cc.

Constructor & Destructor Documentation

◆ DDPixBarLayerUpgradeAlgo()

DDPixBarLayerUpgradeAlgo::DDPixBarLayerUpgradeAlgo ( )

Definition at line 65 of file DDPixBarLayerUpgradeAlgo.cc.

                                                    {
   LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo info: Creating an instance";
 }

References LogDebug.

◆ ~DDPixBarLayerUpgradeAlgo()

DDPixBarLayerUpgradeAlgo::~DDPixBarLayerUpgradeAlgo ( )

override

Definition at line 69 of file DDPixBarLayerUpgradeAlgo.cc.

69 {}

Member Function Documentation

◆ execute()

void DDPixBarLayerUpgradeAlgo::execute ( DDCompactView & cpv )

override

Definition at line 113 of file DDPixBarLayerUpgradeAlgo.cc.

                                                          {
   DDName mother = parent().name();
   const string& idName = mother.name();
  
   double dphi = CLHEP::twopi / number;
   double x2 = coolDist * sin(0.5 * dphi);
   double rtmi = coolDist * cos(0.5 * dphi) - (coolRadius + ladderThick) + rInnerFineTune;
   double rmxh = coolDist * cos(0.5 * dphi) + (coolRadius + ladderThick + ladderOffset) + rOuterFineTune;
   double rtmx = sqrt(rmxh * rmxh + ladderWidth * ladderWidth / 4);
   DDSolid solid = DDSolidFactory::tubs(DDName(idName, idNameSpace), 0.5 * layerDz, rtmi, rtmx, 0, CLHEP::twopi);
   LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << DDName(idName, idNameSpace) << " Tubs made of "
                         << genMat << " from 0 to " << CLHEP::twopi / CLHEP::deg << " with Rin " << rtmi << " Rout "
                         << rtmx << " ZHalf " << 0.5 * layerDz;
   DDName matname(DDSplit(genMat).first, DDSplit(genMat).second);
   DDMaterial matter(matname);
   DDLogicalPart layer(solid.ddname(), matter, solid);
  
   // Full Tubes
   string name = idName + "CoolTube";
   solid = DDSolidFactory::tubs(DDName(name, idNameSpace), 0.5 * coolDz, 0, coolRadius, 0, CLHEP::twopi);
   LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << solid.name() << " Tubs made of " << tubeMat
                         << " from 0 to " << CLHEP::twopi / CLHEP::deg << " with Rout " << coolRadius << " ZHalf "
                         << 0.5 * coolDz;
   matter = DDMaterial(DDName(DDSplit(tubeMat).first, DDSplit(tubeMat).second));
   DDLogicalPart coolTube(solid.ddname(), matter, solid);
  
   // Half Tubes
   name = idName + "CoolTubeHalf";
   solid = DDSolidFactory::tubs(DDName(name, idNameSpace), 0.5 * coolDz, 0, coolRadius, 0, CLHEP::pi);
   LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << solid.name() << " Tubs made of " << tubeMatHalf
                         << " from 0 to " << CLHEP::twopi / CLHEP::deg << " with Rout " << coolRadius << " ZHalf "
                         << 0.5 * coolDz;
   matter = DDMaterial(DDName(DDSplit(tubeMatHalf).first, DDSplit(tubeMatHalf).second));
   DDLogicalPart coolTubeHalf(solid.ddname(), matter, solid);
  
   // Full Coolant
   name = idName + "Coolant";
   solid = DDSolidFactory::tubs(DDName(name, idNameSpace), 0.5 * coolDz, 0, coolRadius - coolThick, 0, CLHEP::twopi);
   LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << solid.name() << " Tubs made of " << tubeMat
                         << " from 0 to " << CLHEP::twopi / CLHEP::deg << " with Rout " << coolRadius - coolThick
                         << " ZHalf " << 0.5 * coolDz;
   matter = DDMaterial(DDName(DDSplit(coolMat).first, DDSplit(coolMat).second));
   DDLogicalPart cool(solid.ddname(), matter, solid);
   cpv.position(cool, coolTube, 1, DDTranslation(0.0, 0.0, 0.0), DDRotation());
   LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << cool.name() << " number 1 positioned in "
                         << coolTube.name() << " at (0,0,0) with no rotation";
  
   // Half Coolant
   name = idName + "CoolantHalf";
   solid = DDSolidFactory::tubs(DDName(name, idNameSpace), 0.5 * coolDz, 0, coolRadius - coolThick, 0, CLHEP::pi);
   LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << solid.name() << " Tubs made of " << tubeMatHalf
                         << " from 0 to " << CLHEP::twopi / CLHEP::deg << " with Rout " << coolRadius - coolThick
                         << " ZHalf " << 0.5 * coolDz;
   matter = DDMaterial(DDName(DDSplit(coolMatHalf).first, DDSplit(coolMatHalf).second));
   DDLogicalPart coolHalf(solid.ddname(), matter, solid);
   cpv.position(coolHalf, coolTubeHalf, 1, DDTranslation(0.0, 0.0, 0.0), DDRotation());
   LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << cool.name() << " number 1 positioned in "
                         << coolTube.name() << " at (0,0,0) with no rotation";
  
   DDName ladderFull(DDSplit(ladder).first, DDSplit(ladder).second);
   int copy = 1, iup = (-1) * outerFirst;
   int copyoffset = number + 2;
   for (int i = 1; i < number + 1; i++) {
     double phi = i * dphi + 90 * CLHEP::deg - 0.5 * dphi + phiFineTune;  //to start with the interface ladder
     double phix, phiy, rrr, rrroffset;
     string rots;
     DDTranslation tran;
     DDRotation rot;
     iup = -iup;
     double dr;
     if ((i == 1) || (i == number / 2 + 1)) {
       dr = coolRadius + 0.5 * ladderThick + ladderOffset;  //interface ladder offset
     } else {
       dr = coolRadius + 0.5 * ladderThick;
     }
     if (i % 2 == 1) {
       rrr = coolDist * cos(0.5 * dphi) + iup * dr + rOuterFineTune;
     } else {
       rrr = coolDist * cos(0.5 * dphi) + iup * dr + rInnerFineTune;
     }
     tran = DDTranslation(rrr * cos(phi), rrr * sin(phi), 0);
     rots = idName + to_string(copy);
     if (iup > 0)
       phix = phi - 90 * CLHEP::deg;
     else
       phix = phi + 90 * CLHEP::deg;
     phiy = phix + 90. * CLHEP::deg;
     LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: Creating a new "
                           << "rotation: " << rots << "\t90., " << phix / CLHEP::deg << ", 90.," << phiy / CLHEP::deg
                           << ", 0, 0";
     rot = DDrot(DDName(rots, idNameSpace), 90 * CLHEP::deg, phix, 90 * CLHEP::deg, phiy, 0., 0.);
     cpv.position(ladderFull, layer, copy, tran, rot);
     LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << ladderFull << " number " << copy << " positioned in "
                           << layer.name() << " at " << tran << " with " << rot;
     copy++;
     rrr = coolDist * cos(0.5 * dphi) + coolRadius / 2.;
     rots = idName + to_string(i + 100);
     phix = phi + 90. * CLHEP::deg;
     if (iup < 0)
       phix += dphi;
     phiy = phix + 90. * CLHEP::deg;
     LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: Creating a new "
                           << "rotation: " << rots << "\t90., " << phix / CLHEP::deg << ", 90.," << phiy / CLHEP::deg
                           << ", 0, 0";
     tran = DDTranslation(rrr * cos(phi) - x2 * sin(phi), rrr * sin(phi) + x2 * cos(phi), 0);
     rot = DDrot(DDName(rots, idNameSpace), 90 * CLHEP::deg, phix, 90 * CLHEP::deg, phiy, 0., 0.);
     cpv.position(coolTubeHalf, layer, i + 1, tran, rot);
     if ((i == 1) || (i == number / 2 + 1)) {
       rrroffset = coolDist * cos(0.5 * dphi) + iup * ladderOffset + rOuterFineTune;
       tran = DDTranslation(
           rrroffset * cos(phi) - cool1Offset * sin(phi), rrroffset * sin(phi) + cool1Offset * cos(phi), 0);
       cpv.position(coolTube, layer, copyoffset, tran, DDRotation());
       copyoffset++;
       tran = DDTranslation(
           rrroffset * cos(phi) - cool2Offset * sin(phi), rrroffset * sin(phi) + cool2Offset * cos(phi), 0);
       cpv.position(coolTube, layer, copyoffset, tran, DDRotation());
       copyoffset++;
     }
     LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo test: " << coolTube.name() << " number " << i + 1
                           << " positioned in " << layer.name() << " at " << tran << " with " << rot;
   }
 }

References filterCSVwithJSON::copy, funct::cos(), DDBase< N, C >::ddname(), DDrot(), DDSplit(), flavorHistoryFilter_cfi::dr, first, mps_fire::i, versionedElectronIDProducer_cfi::idName, PVValHelper::ladder, phase1PixelTopology::layer, LogDebug, Skims_PA_cff::name, DDName::name(), DDBase< N, C >::name(), contentValuesFiles::number, class-composition::parent, pi, DDCompactView::position(), makeMuonMisalignmentScenario::rot, edm::second(), funct::sin(), mathSSE::sqrt(), DDSolidFactory::tubs(), and testProducerWithPsetDescEmpty_cfi::x2.

◆ initialize()

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

override

Definition at line 71 of file DDPixBarLayerUpgradeAlgo.cc.

                                                                                  {
   idNameSpace = DDCurrentNamespace::ns();
   DDName parentName = parent().name();
  
   genMat = sArgs["GeneralMaterial"];
   number = int(nArgs["Ladders"]);
   layerDz = nArgs["LayerDz"];
   coolDz = nArgs["CoolDz"];
   coolThick = nArgs["CoolThick"];
   coolRadius = nArgs["CoolRadius"];
   coolDist = nArgs["CoolDist"];
   cool1Offset = nArgs["Cool1Offset"];
   cool2Offset = nArgs["Cool2Offset"];
   coolMat = sArgs["CoolMaterial"];
   tubeMat = sArgs["CoolTubeMaterial"];
   coolMatHalf = sArgs["CoolMaterialHalf"];
   tubeMatHalf = sArgs["CoolTubeMaterialHalf"];
   phiFineTune = nArgs["PitchFineTune"];
   rOuterFineTune = nArgs["OuterOffsetFineTune"];
   rInnerFineTune = nArgs["InnerOffsetFineTune"];
  
   LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo debug: Parent " << parentName << " NameSpace " << idNameSpace
                         << "\n"
                         << "\tLadders " << number << "\tGeneral Material " << genMat << "\tLength " << layerDz
                         << "\tSpecification of Cooling Pieces:\n"
                         << "\tLength " << coolDz << " Thickness of Shell " << coolThick << " Radial distance "
                         << coolDist << " Materials " << coolMat << ", " << tubeMat;
  
   ladder = sArgs["LadderName"];
   ladderWidth = nArgs["LadderWidth"];
   ladderThick = nArgs["LadderThick"];
   ladderOffset = nArgs["LadderOffset"];
   outerFirst = int(nArgs["OuterFirst"]);
  
   LogDebug("PixelGeom") << "DDPixBarLayerUpgradeAlgo debug: Full Ladder " << ladder << " width/thickness "
                         << ladderWidth << ", " << ladderThick;
 }