#include <DDPixFwdBlades.h>

Inheritance diagram for DDPixFwdBlades:

Public Member Functions
	DDPixFwdBlades ()

void	execute (DDCompactView &cpv) override

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

	~DDPixFwdBlades () override

Private Member Functions
void	computeNippleParameters (double endcap)

int	issueCopyNumber ()

Private Attributes
double	ancorRadius

double	bladeAngle

double	bladeZShift

std::string	childName

std::string	childRotationName

std::vector< double >	childTranslationVector

std::map< std::string, int >	copyNumbers

double	endcap

std::string	flagSelector

std::string	flagString

std::string	idNameSpace

double	jX

double	jY

double	jZ

double	kX

double	kY

double	kZ

int	nBlades

CLHEP::HepRotation *	nippleRotationZMinus

CLHEP::HepRotation *	nippleRotationZPlus

double	nippleTranslationX

double	nippleTranslationY

double	nippleTranslationZ

double	zPlane

Detailed Description

Definition at line 64 of file DDPixFwdBlades.h.

Constructor & Destructor Documentation

DDPixFwdBlades::DDPixFwdBlades ( )

Definition at line 30 of file DDPixFwdBlades.cc.

30 {}

DDPixFwdBlades::~DDPixFwdBlades ( )

override

Definition at line 31 of file DDPixFwdBlades.cc.

31 {}

Member Function Documentation

void DDPixFwdBlades::computeNippleParameters ( double endcap )

private

Definition at line 218 of file DDPixFwdBlades.cc.

References ancorRadius, bladeAngle, bladeZShift, funct::cos(), DDrot(), hiPixelPairStep_cff::deltaPhi, jX, jY, jZ, kX, kY, kZ, LogDebug, nBlades, nippleRotationZMinus, nippleRotationZPlus, nippleTranslationX, nippleTranslationY, nippleTranslationZ, funct::sin(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by execute().

                                                           {
     
   double effBladeAngle = endcap * bladeAngle;
   
   CLHEP::Hep3Vector jC; // Point J in the "cover" blade frame
   CLHEP::Hep3Vector kB; // Point K in the "body" blade frame
   std::string rotNameNippleToCover;
   std::string rotNameCoverToNipple;
   std::string rotNameNippleToBody;
 
   if (endcap > 0.) {
     jC = CLHEP::Hep3Vector(jX, jY + ancorRadius, jZ);
     kB = CLHEP::Hep3Vector(kX, kY + ancorRadius, kZ);
     rotNameNippleToCover = "NippleToCoverZPlus";
     rotNameCoverToNipple = "CoverToNippleZPlus";
     rotNameNippleToBody = "NippleToBodyZPlus";
   } else {
     jC = CLHEP::Hep3Vector(-jX, jY + ancorRadius, jZ);
     kB = CLHEP::Hep3Vector(-kX, kY + ancorRadius, kZ);
     rotNameNippleToCover = "NippleToCoverZMinus";
     rotNameCoverToNipple = "CoverToNippleZMinus";
     rotNameNippleToBody = "NippleToBodyZMinus";
   }
         
   // Z-shift from "cover" to "body" blade frame:
   
   CLHEP::Hep3Vector tCB(bladeZShift*sin(effBladeAngle), 0., bladeZShift*cos(effBladeAngle));
   
   // Rotation from "cover" blade frame into "body" blade frame :
   
   double deltaPhi = endcap*(360./nBlades)*CLHEP::deg;
   CLHEP::HepRotation rCB(CLHEP::Hep3Vector(1.*sin(effBladeAngle), 0., 1.*cos(effBladeAngle)), deltaPhi);
   
   // Transform vector k into "cover" blade frame :
   
   CLHEP::Hep3Vector kC = rCB * (kB + tCB);
   
   // Vector JK in the "cover" blade frame:
   
   CLHEP::Hep3Vector jkC = kC - jC;
   double jkLength = jkC.mag();
   DDConstant JK(DDName("JK", "pixfwdNipple"), std::make_unique<double>( jkLength ));
   LogDebug("PixelGeom") << "+++++++++++++++ DDPixFwdBlades: " << "JK Length " <<  jkLength * CLHEP::mm;
   
   // Position of the center of a nipple in "cover" blade frame :
   
   CLHEP::Hep3Vector nippleTranslation((kC+jC)/2. - CLHEP::Hep3Vector(0., ancorRadius, 0.));
   if (endcap > 0) {
     nippleTranslationX = nippleTranslation.x();
     nippleTranslationY = nippleTranslation.y();
     nippleTranslationZ = nippleTranslation.z();
   }
   LogDebug("PixelGeom") << "Child translation : " << nippleTranslation;
   
   // Rotations from nipple frame to "cover" blade frame and back :
   
   CLHEP::Hep3Vector vZ(0.,0.,1.);
   CLHEP::Hep3Vector axis = vZ.cross(jkC);
   double angleCover = vZ.angle(jkC);
   LogDebug("PixelGeom") << " Angle to Cover: " << angleCover;
   CLHEP::HepRotation* rpCN = new CLHEP::HepRotation(axis, angleCover);
   if (endcap > 0.) {
     nippleRotationZPlus = rpCN;
   } else {
     nippleRotationZMinus = rpCN;
   }
   //( endcap > 0. ? nippleRotationZPlus : nippleRotationZMinus ) = rpCN;
 
   DDrot(DDName(rotNameCoverToNipple, "pixfwdNipple"),
     std::make_unique<DDRotationMatrix>( rpCN->xx(), rpCN->xy(), rpCN->xz(),
                         rpCN->yx(), rpCN->yy(), rpCN->yz(),
                         rpCN->zx(), rpCN->zy(), rpCN->zz()));
   CLHEP::HepRotation rpNC(axis, -angleCover);
 
   DDrot(DDName(rotNameNippleToCover, "pixfwdNipple"),
     std::make_unique<DDRotationMatrix>( rpNC.xx(), rpNC.xy(), rpNC.xz(),
                         rpNC.yx(), rpNC.yy(), rpNC.yz(),
                         rpNC.zx(), rpNC.zy(), rpNC.zz()));
   
   // Rotation from nipple frame to "body" blade frame :
   
   CLHEP::HepRotation rpNB( rpNC * rCB );
 
   DDrot(DDName(rotNameNippleToBody, "pixfwdNipple"),
     std::make_unique<DDRotationMatrix>( rpNB.xx(), rpNB.xy(), rpNB.xz(),
                         rpNB.yx(), rpNB.yy(), rpNB.yz(),
                         rpNB.zx(), rpNB.zy(), rpNB.zz()));
   double angleBody = vZ.angle(rpNB * vZ);
   LogDebug("PixelGeom") << " Angle to body : " << angleBody;  
 }

void DDPixFwdBlades::execute ( DDCompactView & cpv )

override

Definition at line 107 of file DDPixFwdBlades.cc.

References ancorRadius, bladeAngle, bladeZShift, class-composition::child, childName, childRotationName, childTranslationVector, computeNippleParameters(), popcon2dropbox::copy(), DDrot(), DDSplit(), hiPixelPairStep_cff::deltaPhi, endcap, plotBeamSpotDB::first, flagSelector, flagString, idNameSpace, issueCopyNumber(), DDName::name(), nBlades, nippleRotationZMinus, nippleRotationZPlus, nippleTranslationX, nippleTranslationY, nippleTranslationZ, class-composition::parent, phi, DDCompactView::position(), idealTransformation::rotation, DDRotation::rotation(), edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, x, y, z, and zPlane.

                                                {
 
   // -- Compute Nipple parameters if not already computed :
   
   if (!nippleRotationZPlus) {
     computeNippleParameters(1.);   // Z Plus endcap
     computeNippleParameters(-1.);  // Z Minus endcap
   }
   if (childName.empty()) return;
   
   // -- Signed versions of blade angle and z-shift :
   
   double effBladeAngle = - endcap * bladeAngle;
   double effBladeZShift = endcap * bladeZShift;
   
   // -- Names of mother and child volumes :
 
   DDName mother = parent().name();
   DDName child(DDSplit(childName).first, DDSplit(childName).second);
   
   // -- Get translation and rotation from "blade frame" to "child frame", if any :
   
   CLHEP::HepRotation childRotMatrix = CLHEP::HepRotation();
   if (!childRotationName.empty()) {
     DDRotation childRotation = DDRotation(DDName(DDSplit(childRotationName).first, DDSplit(childRotationName).second));
     // due to conversion to ROOT::Math::Rotation3D -- Michael Case
     DD3Vector x, y, z;
     childRotation.rotation().GetComponents(x, y, z); // these are the orthonormal columns.
     CLHEP::HepRep3x3 tr(x.X(), y.X(), z.X(), x.Y(), y.Y(), z.Y(), x.Z(), y.Z(), z.Z());
     childRotMatrix = CLHEP::HepRotation(tr);
   } else if (childName == "pixfwdNipple:PixelForwardNippleZPlus") {
     childRotMatrix = *nippleRotationZPlus;
   } else if (childName == "pixfwdNipple:PixelForwardNippleZMinus") {
     childRotMatrix = *nippleRotationZMinus;
   }
   
   CLHEP::Hep3Vector childTranslation;
   if (childName == "pixfwdNipple:PixelForwardNippleZPlus") {
     childTranslation = CLHEP::Hep3Vector(nippleTranslationX, nippleTranslationY, nippleTranslationZ);
   } else if (childName == "pixfwdNipple:PixelForwardNippleZMinus") {
     childTranslation = CLHEP::Hep3Vector(-nippleTranslationX, nippleTranslationY, nippleTranslationZ);
   } else {
     childTranslation = CLHEP::Hep3Vector(childTranslationVector[0],childTranslationVector[1],childTranslationVector[2]);
   }
   
   // Create a matrix for rotation around blade axis (to "blade frame") :
   
   CLHEP::HepRotation bladeRotMatrix(CLHEP::Hep3Vector(0.,1.,0.),effBladeAngle);
   
   // Cycle over Phi positions, placing copies of the child volume :
 
   double deltaPhi = (360./nBlades)*CLHEP::deg;
   int nQuarter = nBlades/4;
   double zShiftMax = effBladeZShift*((nQuarter-1)/2.);
 
   for (int iBlade=0; iBlade < nBlades; iBlade++) {
     
     // check if this blade position should be skipped :
     
     if (flagString[iBlade] != flagSelector[0]) continue;
     int copy = issueCopyNumber();
     
     // calculate Phi and Z shift for this blade :
 
     double phi = (iBlade + 0.5) * deltaPhi - 90.*CLHEP::deg;
     int iQuarter = iBlade % nQuarter;
     double zShift = - zShiftMax + iQuarter * effBladeZShift;
     
     // compute rotation matrix from mother to blade frame :
     
     CLHEP::HepRotation rotMatrix(CLHEP::Hep3Vector(0.,0.,1.), phi);
     rotMatrix *= bladeRotMatrix;
 
     // convert translation vector from blade frame to mother frame, and add Z shift :
     
     CLHEP::Hep3Vector translation = rotMatrix(childTranslation + CLHEP::Hep3Vector(0., ancorRadius, 0.));
     translation += CLHEP::Hep3Vector(0., 0., zShift + zPlane);
     
     // create DDRotation for placing the child if not already existent :
 
     DDRotation rotation;   
     std::string rotstr = mother.name() + DDSplit(childName).first + std::to_string(copy);
     rotation = DDRotation(DDName(rotstr, idNameSpace));
 
     if (!rotation) {
       rotMatrix *= childRotMatrix;
       rotation = DDrot(DDName(rotstr, idNameSpace),
                std::make_unique<DDRotationMatrix>(rotMatrix.xx(), rotMatrix.xy(), rotMatrix.xz(),
                               rotMatrix.yx(), rotMatrix.yy(), rotMatrix.yz(),
                               rotMatrix.zx(), rotMatrix.zy(), rotMatrix.zz()));
     }
     // position the child :
 
     DDTranslation ddtran(translation.x(), translation.y(), translation.z());
     cpv.position(child, mother, copy, ddtran, rotation);
   }
 
   // End of cycle over Phi positions
 
 }

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

override

Definition at line 35 of file DDPixFwdBlades.cc.

References ancorRadius, bladeAngle, bladeZShift, childName, childRotationName, childTranslationVector, copyNumbers, endcap, flagSelector, flagString, idNameSpace, jX, jY, jZ, kX, kY, kZ, nBlades, nippleRotationZMinus, nippleRotationZPlus, nippleTranslationX, nippleTranslationY, nippleTranslationZ, DDCurrentNamespace::ns(), and zPlane.

                                                   {
                     
   if ( nArgs.find("Endcap") != nArgs.end() ) {
     endcap = nArgs["Endcap"];
   } else {
     endcap = 1.;
   }
 
   if ( sArgs.find("FlagString") != sArgs.end() ) {
     flagString = sArgs["FlagString"];
     flagSelector = sArgs["FlagSelector"];
   } else {
     flagString = "YYYYYYYYYYYYYYYYYYYYYYYY";
     flagSelector = "Y";
   }
 
   if ( sArgs.find("Child") != sArgs.end() ) {
     childName   = sArgs["Child"];
   } else {
     childName   = "";
   }
 
   if ( vArgs.find("ChildTranslation") != vArgs.end() ) {
     childTranslationVector = vArgs["ChildTranslation"];
   } else {
     childTranslationVector = std::vector<double>(3, 0.);
   }
 
   if ( sArgs.find("ChildRotation") != sArgs.end() ) {
     childRotationName = sArgs["ChildRotation"];
   } else {
     childRotationName = "";
   }
 
   idNameSpace = DDCurrentNamespace::ns();
 
   // -- Input geometry parameters :  -----------------------------------------------------
 
   nBlades = 24;                     // Number of blades
   bladeAngle = 20.*CLHEP::deg;   // Angle of blade rotation around its axis
   zPlane = 0.;                   // Common shift in Z for all blades (with respect to disk center plane)
   bladeZShift = 6.*CLHEP::mm;    // Shift in Z between the axes of two adjacent blades
   
   ancorRadius = 54.631*CLHEP::mm;// Distance from beam line to ancor point defining center of "blade frame"
   
   // Coordinates of Nipple ancor points J and K in "blade frame" :
 
   jX = -16.25*CLHEP::mm;
   jY = 96.50*CLHEP::mm;
   jZ = 1.25*CLHEP::mm;
   kX = 16.25*CLHEP::mm;
   kY = 96.50*CLHEP::mm;
   kZ = -1.25*CLHEP::mm;
   
   // -- Static initialization :  -----------------------------------------------------------
 
   nippleRotationZPlus = nullptr;
   nippleRotationZMinus = nullptr;
   nippleTranslationX = 0.;
   nippleTranslationY = 0.;
   nippleTranslationZ = 0.;
 
   copyNumbers.clear();
 
 }

int DDPixFwdBlades::issueCopyNumber ( )

private

Definition at line 210 of file DDPixFwdBlades.cc.

References childName, and copyNumbers.

Referenced by execute().

                                     {
   if (copyNumbers.count(childName) == 0) copyNumbers[childName] = 0;
   return ++copyNumbers[childName];
 }

Member Data Documentation

double DDPixFwdBlades::ancorRadius

private

Definition at line 94 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), execute(), and initialize().

double DDPixFwdBlades::bladeAngle

private

Definition at line 90 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), execute(), and initialize().

double DDPixFwdBlades::bladeZShift

private

Definition at line 92 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), execute(), and initialize().

std::string DDPixFwdBlades::childName

private

Definition at line 112 of file DDPixFwdBlades.h.

Referenced by execute(), initialize(), and issueCopyNumber().

std::string DDPixFwdBlades::childRotationName

private

Definition at line 115 of file DDPixFwdBlades.h.

Referenced by execute(), and initialize().

std::vector<double> DDPixFwdBlades::childTranslationVector

private

Definition at line 114 of file DDPixFwdBlades.h.

Referenced by execute(), and initialize().

std::map<std::string, int> DDPixFwdBlades::copyNumbers

private

Definition at line 121 of file DDPixFwdBlades.h.

Referenced by initialize(), and issueCopyNumber().

double DDPixFwdBlades::endcap

private

Definition at line 107 of file DDPixFwdBlades.h.

Referenced by execute(), geometryXMLparser.CSCAlignable::index(), and initialize().

std::string DDPixFwdBlades::flagSelector

private

Definition at line 110 of file DDPixFwdBlades.h.

Referenced by execute(), and initialize().

std::string DDPixFwdBlades::flagString

private

Definition at line 109 of file DDPixFwdBlades.h.

Referenced by execute(), and initialize().

std::string DDPixFwdBlades::idNameSpace

private

Definition at line 119 of file DDPixFwdBlades.h.

Referenced by execute(), and initialize().

double DDPixFwdBlades::jX

private

Definition at line 98 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), and initialize().

double DDPixFwdBlades::jY

private

Definition at line 99 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), and initialize().

double DDPixFwdBlades::jZ

private

Definition at line 100 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), and initialize().

double DDPixFwdBlades::kX

private

Definition at line 101 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), and initialize().

double DDPixFwdBlades::kY

private

Definition at line 102 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), and initialize().

double DDPixFwdBlades::kZ

private

Definition at line 103 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), and initialize().

int DDPixFwdBlades::nBlades

private

Definition at line 89 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), execute(), and initialize().

CLHEP::HepRotation* DDPixFwdBlades::nippleRotationZMinus

private

Definition at line 124 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), execute(), and initialize().

CLHEP::HepRotation* DDPixFwdBlades::nippleRotationZPlus

private

Definition at line 123 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), execute(), and initialize().

double DDPixFwdBlades::nippleTranslationX

private

Definition at line 125 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), execute(), and initialize().

double DDPixFwdBlades::nippleTranslationY

private

Definition at line 125 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), execute(), and initialize().

double DDPixFwdBlades::nippleTranslationZ

private

Definition at line 125 of file DDPixFwdBlades.h.

Referenced by computeNippleParameters(), execute(), and initialize().

double DDPixFwdBlades::zPlane

private

Definition at line 91 of file DDPixFwdBlades.h.

Referenced by execute(), and initialize().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation