DDPixFwdBlades Class Reference

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
string	childName
string	childRotationName
vector< double >	childTranslationVector
map< string, int >	copyNumbers
double	endcap
string	flagSelector
string	flagString
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 83 of file DDPixFwdBlades.cc.

Constructor & Destructor Documentation

DDPixFwdBlades()

DDPixFwdBlades::DDPixFwdBlades

(

)

Definition at line 134 of file DDPixFwdBlades.cc.

{}

~DDPixFwdBlades()

DDPixFwdBlades::~DDPixFwdBlades ( )

override

Definition at line 135 of file DDPixFwdBlades.cc.

{}

Member Function Documentation

computeNippleParameters()

void DDPixFwdBlades::computeNippleParameters ( double  endcap )

private

Definition at line 320 of file DDPixFwdBlades.cc.

                                                          {
  double effBladeAngle = endcap * bladeAngle;
 
  CLHEP::Hep3Vector jC;  // Point J in the "cover" blade frame
  CLHEP::Hep3Vector kB;  // Point K in the "body" blade frame
  string rotNameNippleToCover;
  string rotNameCoverToNipple;
  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"), 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"),
      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"),
        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"),
        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;
}

References funct::cos(), DDrot(), SiPixelRawToDigiRegional_cfi::deltaPhi, makeMuonMisalignmentScenario::endcap, LogDebug, and funct::sin().

execute()

void DDPixFwdBlades::execute ( DDCompactView &  cpv )

override

Definition at line 205 of file DDPixFwdBlades.cc.

                                               {
  // -- 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;
    string rotstr = mother.name() + DDSplit(childName).first + to_string(copy);
    rotation = DDRotation(DDName(rotstr, idNameSpace));
 
    if (!rotation) {
      rotMatrix *= childRotMatrix;
      rotation = DDrot(DDName(rotstr, idNameSpace),
                       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
}

References submitPVValidationJobs::child, filterCSVwithJSON::copy, DDrot(), DDSplit(), SiPixelRawToDigiRegional_cfi::deltaPhi, makeMuonMisalignmentScenario::endcap, first, DDName::name(), class-composition::parent, DDCompactView::position(), idealTransformation::rotation, DDRotation::rotation(), and edm::second().

initialize()

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

override

Definition at line 137 of file DDPixFwdBlades.cc.

                                                                {
  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 = 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();
}

References makeMuonMisalignmentScenario::endcap, and DDCurrentNamespace::ns().

issueCopyNumber()

int DDPixFwdBlades::issueCopyNumber ( )

private

Definition at line 312 of file DDPixFwdBlades.cc.

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

Member Data Documentation

ancorRadius

double DDPixFwdBlades::ancorRadius

private

Definition at line 102 of file DDPixFwdBlades.cc.

bladeAngle

double DDPixFwdBlades::bladeAngle

private

Definition at line 98 of file DDPixFwdBlades.cc.

bladeZShift

double DDPixFwdBlades::bladeZShift

private

Definition at line 100 of file DDPixFwdBlades.cc.

childName

string DDPixFwdBlades::childName

private

Definition at line 118 of file DDPixFwdBlades.cc.

childRotationName

string DDPixFwdBlades::childRotationName

private

Definition at line 121 of file DDPixFwdBlades.cc.

childTranslationVector

vector<double> DDPixFwdBlades::childTranslationVector

private

Definition at line 120 of file DDPixFwdBlades.cc.

copyNumbers

map<string, int> DDPixFwdBlades::copyNumbers

private

Definition at line 124 of file DDPixFwdBlades.cc.

endcap

double DDPixFwdBlades::endcap

private

Definition at line 113 of file DDPixFwdBlades.cc.

Referenced by geometryXMLparser.CSCAlignable::index().

flagSelector

string DDPixFwdBlades::flagSelector

private

Definition at line 116 of file DDPixFwdBlades.cc.

flagString

string DDPixFwdBlades::flagString

private

Definition at line 115 of file DDPixFwdBlades.cc.

idNameSpace

string DDPixFwdBlades::idNameSpace

private

Definition at line 122 of file DDPixFwdBlades.cc.

jX

double DDPixFwdBlades::jX

private

Definition at line 106 of file DDPixFwdBlades.cc.

jY

double DDPixFwdBlades::jY

private

Definition at line 107 of file DDPixFwdBlades.cc.

jZ

double DDPixFwdBlades::jZ

private

Definition at line 108 of file DDPixFwdBlades.cc.

kX

double DDPixFwdBlades::kX

private

Definition at line 109 of file DDPixFwdBlades.cc.

kY

double DDPixFwdBlades::kY

private

Definition at line 110 of file DDPixFwdBlades.cc.

kZ

double DDPixFwdBlades::kZ

private

Definition at line 111 of file DDPixFwdBlades.cc.

nBlades

int DDPixFwdBlades::nBlades

private

Definition at line 97 of file DDPixFwdBlades.cc.

nippleRotationZMinus

CLHEP::HepRotation* DDPixFwdBlades::nippleRotationZMinus

private

Definition at line 127 of file DDPixFwdBlades.cc.

nippleRotationZPlus

CLHEP::HepRotation* DDPixFwdBlades::nippleRotationZPlus

private

Definition at line 126 of file DDPixFwdBlades.cc.

nippleTranslationX

double DDPixFwdBlades::nippleTranslationX

private

Definition at line 128 of file DDPixFwdBlades.cc.

nippleTranslationY

double DDPixFwdBlades::nippleTranslationY

private

Definition at line 128 of file DDPixFwdBlades.cc.

nippleTranslationZ

double DDPixFwdBlades::nippleTranslationZ

private

Definition at line 128 of file DDPixFwdBlades.cc.

zPlane

double DDPixFwdBlades::zPlane

private

Definition at line 99 of file DDPixFwdBlades.cc.