Inheritance diagram for DDPixFwdBladesNew:

Public Member Functions
	DDPixFwdBladesNew ()

void	execute (DDCompactView &cpv) override

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

	~DDPixFwdBladesNew () override=default

Private Member Functions
CLHEP::HepRotation	getRotation ()

CLHEP::Hep3Vector	getTranslation ()

Private Attributes
double	ancorRadius_

double	bladeAngle_

double	bladeZShift_

std::string	childName_

std::string	childRotationName_

std::vector< double >	childTranslationVector_

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_

int	nippleType_

int	startCopy_

double	zPlane_

Detailed Description

Definition at line 30 of file DDPixFwdBladesNew.cc.

Constructor & Destructor Documentation

◆ DDPixFwdBladesNew()

DDPixFwdBladesNew::DDPixFwdBladesNew ( )

inline

Definition at line 32 of file DDPixFwdBladesNew.cc.

32 {}

◆ ~DDPixFwdBladesNew()

DDPixFwdBladesNew::~DDPixFwdBladesNew ( )

overridedefault

Member Function Documentation

◆ execute()

void DDPixFwdBladesNew::execute ( DDCompactView & cpv )

override

Definition at line 104 of file DDPixFwdBladesNew.cc.

References ancorRadius_, bladeAngle_, bladeZShift_, submitPVValidationJobs::child, childName_, childRotationName_, childTranslationVector_, filterCSVwithJSON::copy, DDrot(), DDSplit(), SiPixelRawToDigiRegional_cfi::deltaPhi, endcap_, first, flagSelector_, flagString_, getRotation(), getTranslation(), idNameSpace_, DDName::name(), nBlades_, nippleType_, class-composition::parent, phi, DDCompactView::position(), idealTransformation::rotation, DDRotation::rotation(), edm::second(), startCopy_, AlCaHLTBitMon_QueryRunRegistry::string, cond::impl::to_string(), x, y, z, and zPlane_.

                                                   {
   // -- 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 (nippleType_ == 1) {
     childRotMatrix = getRotation();
   } else 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);
   }
 
   CLHEP::Hep3Vector childTranslation =
       (nippleType_ == 1)
           ? getTranslation()
           : 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.);
   int copy(startCopy_);
 
   for (int iBlade = 0; iBlade < nBlades_; iBlade++) {
     // check if this blade position should be skipped :
 
     if (flagString_[iBlade] != flagSelector_[0])
       continue;
 
     // 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_));
     edm::LogVerbatim("PixelGeom") << "DDPixFwdBlades: Rotation " << rotstr << " : " << rotation;
 
     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);
     edm::LogVerbatim("PixelGeom") << "DDPixFwdBlades::Position " << child << " copy " << copy << " in " << mother
                                   << " with translation " << ddtran << " and rotation " << rotation;
     ++copy;
   }
 
   // End of cycle over Phi positions
 }

◆ getRotation()

CLHEP::HepRotation DDPixFwdBladesNew::getRotation ( )

private

Definition at line 222 of file DDPixFwdBladesNew.cc.

References ancorRadius_, bladeAngle_, bladeZShift_, funct::cos(), SiPixelRawToDigiRegional_cfi::deltaPhi, endcap_, jX_, jY_, jZ_, kX_, kY_, kZ_, nBlades_, and funct::sin().

Referenced by execute().

                                                 {
   double effBladeAngle = endcap_ * bladeAngle_;
 
   CLHEP::Hep3Vector jC =
       CLHEP::Hep3Vector(endcap_ * jX_, jY_ + ancorRadius_, jZ_);  // Point J in the "cover" blade frame
   CLHEP::Hep3Vector kB =
       CLHEP::Hep3Vector(endcap_ * kX_, kY_ + ancorRadius_, kZ_);  // Point K in the "body" blade frame
 
   // 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);
   CLHEP::Hep3Vector jkC = kC - jC;
   edm::LogVerbatim("PixelGeom") << "+++++++++++++++ DDPixFwdBlades: "
                                 << "JK Length " << jkC.mag() * CLHEP::mm;
 
   // 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);
   edm::LogVerbatim("PixelGeom") << " Angle to Cover: " << angleCover;
   CLHEP::HepRotation rpCN(axis, angleCover);
   return rpCN;
 }

◆ getTranslation()

CLHEP::Hep3Vector DDPixFwdBladesNew::getTranslation ( )

private

Definition at line 198 of file DDPixFwdBladesNew.cc.

References ancorRadius_, bladeAngle_, bladeZShift_, funct::cos(), SiPixelRawToDigiRegional_cfi::deltaPhi, endcap_, jX_, jY_, jZ_, kX_, kY_, kZ_, nBlades_, and funct::sin().

Referenced by execute().

                                                   {
   double effBladeAngle = endcap_ * bladeAngle_;
 
   CLHEP::Hep3Vector jC =
       CLHEP::Hep3Vector(endcap_ * jX_, jY_ + ancorRadius_, jZ_);  // Point J in the "cover" blade frame
   CLHEP::Hep3Vector kB =
       CLHEP::Hep3Vector(endcap_ * kX_, kY_ + ancorRadius_, kZ_);  // Point K in the "body" blade frame
 
   // 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);
 
   // Position of the center of a nipple in "cover" blade frame :
   CLHEP::Hep3Vector nippleTranslation((kC + jC) / 2. - CLHEP::Hep3Vector(0., ancorRadius_, 0.));
   edm::LogVerbatim("PixelGeom") << "Child translation : " << nippleTranslation;
   return nippleTranslation;
 }

◆ initialize()

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

override

Definition at line 65 of file DDPixFwdBladesNew.cc.

References ancorRadius_, bladeAngle_, bladeZShift_, childName_, childRotationName_, childTranslationVector_, endcap_, flagSelector_, flagString_, idNameSpace_, jX_, jY_, jZ_, kX_, kY_, kZ_, nBlades_, nippleType_, DDCurrentNamespace::ns(), startCopy_, and zPlane_.

                                                                    {
   endcap_ = nArgs["Endcap"];
   nBlades_ = static_cast<int>(nArgs["Blades"]);  // Number of blades
   bladeAngle_ = nArgs["BladeAngle"];             // Angle of blade rotation around its axis
   bladeZShift_ = nArgs["BladeZShift"];           // Shift in Z between the axes of two adjacent blades
   ancorRadius_ = nArgs["AncorRadius"];  // Distance from beam line to ancor point defining center of "blade frame"
   // Coordinates of Nipple ancor points J and K in "blade frame" :
   nippleType_ = static_cast<int>(nArgs["NippleType"]);
   jX_ = nArgs["JX"];
   jY_ = nArgs["JY"];
   jZ_ = nArgs["JZ"];
   kX_ = nArgs["KX"];
   kY_ = nArgs["KY"];
   kZ_ = nArgs["KZ"];
 
   flagString_ = sArgs["FlagString"];
   flagSelector_ = sArgs["FlagSelector"];
   childName_ = sArgs["Child"];
   startCopy_ = static_cast<int>(nArgs["StartCopy"]);
   childTranslationVector_ = vArgs["ChildTranslation"];
   childRotationName_ = sArgs["ChildRotation"];
 
   idNameSpace_ = DDCurrentNamespace::ns();
 
   edm::LogVerbatim("PixelGeom") << "DDPixFwdBladesNew: Initialize with endcap " << endcap_ << " FlagString "
                                 << flagString_ << " FlagSelector " << flagSelector_ << " Child " << childName_
                                 << " ChildTranslation " << childTranslationVector_[0] << ":"
                                 << childTranslationVector_[1] << ":" << childTranslationVector_[2] << " ChildRotation "
                                 << childRotationName_ << " NameSpace " << idNameSpace_ << "\n  nBlades " << nBlades_
                                 << " bladeAngle " << bladeAngle_ << " zPlane " << zPlane_ << " bladeZShift "
                                 << bladeZShift_ << " ancorRadius " << ancorRadius_ << " NippleType " << nippleType_
                                 << " jX|jY|jZ " << jX_ << ":" << jY_ << ":" << jZ_ << " kX|kY|kZ " << kX_ << ":" << kY_
                                 << ":" << kZ_;
 }