Inheritance diagram for DDTrackerIrregularRingAlgo:

Public Member Functions
	DDTrackerIrregularRingAlgo ()

void	execute (DDCompactView &cpv) override

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

	~DDTrackerIrregularRingAlgo () override

Private Attributes
vector< double >	center

string	childName

string	idNameSpace

int	incrCopyNo

bool	isFlipped

bool	isZPlus

int	n

vector< double >	phiAngles

double	radius

vector< double >	radiusValues

double	rangeAngle

double	startAngle

int	startCopyNo

double	tiltAngle

vector< double >	yawAngles

Detailed Description

Definition at line 54 of file DDTrackerIrregularRingAlgo.cc.

Constructor & Destructor Documentation

DDTrackerIrregularRingAlgo::DDTrackerIrregularRingAlgo ( )

Definition at line 87 of file DDTrackerIrregularRingAlgo.cc.

References LogDebug.

                                                        {
   LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo info: Creating an instance";
 }

DDTrackerIrregularRingAlgo::~DDTrackerIrregularRingAlgo ( )

overridedefault

Member Function Documentation

void DDTrackerIrregularRingAlgo::execute ( DDCompactView & cpv )

override

Definition at line 125 of file DDTrackerIrregularRingAlgo.cc.

References submitPVValidationJobs::child, angle_units::operators::convertDegToRad(), angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, funct::cos(), DDrot(), DDSplit(), first, mps_fire::i, LogDebug, DDRotation::matrix(), dqmiodumpmetadata::n, SpecificationBuilder_cfi::parent(), DDCompactView::position(), CosmicsPD_Skims::radius, edm::second(), funct::sin(), theta(), and cond::impl::to_string().

                                                            {
   DDRotation flipRot, tiltRot, phiOwnAxisRot, phiRot, globalRot;                                // Identity
   DDRotationMatrix flipMatrix, tiltMatrix, phiOwnAxisRotMatrix, phiRotMatrix, globalRotMatrix;  // Identity matrix
   string rotstr = "RTrackerRingAlgo";
 
   // flipMatrix calculus
   if (isFlipped) {
     string flipRotstr = rotstr + "Flip";
     flipRot = DDRotation(DDName(flipRotstr, idNameSpace));
     if (!flipRot) {
       LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new rotation: " << flipRotstr
                               << "\t90., 180., "
                               << "90., 90., "
                               << "180., 0.";
       flipRot = DDrot(DDName(flipRotstr, idNameSpace), 90._deg, 180._deg, 90._deg, 90._deg, 180._deg, 0.);
     }
     flipMatrix = flipRot.matrix();
   }
   // tiltMatrix calculus
   if (isZPlus) {
     string tiltRotstr = rotstr + "Tilt" + to_string(convertRadToDeg(tiltAngle)) + "ZPlus";
     tiltRot = DDRotation(DDName(tiltRotstr, idNameSpace));
     if (!tiltRot) {
       LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new rotation: " << tiltRotstr
                               << "\t90., 90., " << convertRadToDeg(tiltAngle) << ", 180., "
                               << 90. - convertRadToDeg(tiltAngle) << ", 0.";
       tiltRot = DDrot(DDName(tiltRotstr, idNameSpace), 90._deg, 90._deg, tiltAngle, 180._deg, 90._deg - tiltAngle, 0.);
     }
     tiltMatrix = tiltRot.matrix();
     if (isFlipped) {
       tiltMatrix *= flipMatrix;
     }
   } else {
     string tiltRotstr = rotstr + "Tilt" + to_string(convertRadToDeg(tiltAngle)) + "ZMinus";
     tiltRot = DDRotation(DDName(tiltRotstr, idNameSpace));
     if (!tiltRot) {
       LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new rotation: " << tiltRotstr
                               << "\t90., 90., " << convertRadToDeg(tiltAngle) << ", 0., "
                               << 90. + convertRadToDeg(tiltAngle) << ", 0.";
       tiltRot = DDrot(DDName(tiltRotstr, idNameSpace), 90._deg, 90._deg, tiltAngle, 0., 90._deg + tiltAngle, 0.);
     }
     tiltMatrix = tiltRot.matrix();
     if (isFlipped) {
       tiltMatrix *= flipMatrix;
     }
   }
 
   // Loops for all phi values
   DDName mother = parent().name();
   DDName child(DDSplit(childName).first, DDSplit(childName).second);
   double theta = 90._deg;
   int copy = startCopyNo;
   //double phi = startAngle;
 
   for (int i = 0; i < n; i++) {
     // phiRotMatrix calculus
     //double phix = phi;
     //double phix_ownaxis = 0._deg;
     double phix = convertDegToRad(phiAngles.at(i));
     double phix_ownaxis = convertDegToRad(yawAngles.at(i));
     radius = radiusValues.at(i);
     double phiy = phix + 90._deg;
     double phiy_ownaxis = phix_ownaxis + 90._deg;
     double phideg = convertRadToDeg(phix);
     double phideg_ownaxis = convertRadToDeg(phix_ownaxis);
     if (phideg_ownaxis != 0) {
       string phiOwnAxisRotstr = rotstr + "PhiOwnAxis" + to_string(phideg_ownaxis * 10.);
       phiOwnAxisRot = DDRotation(DDName(phiOwnAxisRotstr, idNameSpace));
       if (!phiOwnAxisRot) {
         LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new rotation: " << phiOwnAxisRotstr
                                 << "\t90., " << convertRadToDeg(phix_ownaxis) << ", 90.,"
                                 << convertRadToDeg(phiy_ownaxis) << ", 0., 0.";
         phiOwnAxisRot = DDrot(DDName(phiOwnAxisRotstr, idNameSpace), theta, phix_ownaxis, theta, phiy_ownaxis, 0., 0.);
       }
       phiOwnAxisRotMatrix = phiOwnAxisRot.matrix();
     }
     if (phideg != 0) {
       string phiRotstr = rotstr + "Phi" + to_string(phideg * 10.);
       phiRot = DDRotation(DDName(phiRotstr, idNameSpace));
       if (!phiRot) {
         LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new rotation: " << phiRotstr
                                 << "\t90., " << convertRadToDeg(phix) << ", 90.," << convertRadToDeg(phiy)
                                 << ", 0., 0.";
         phiRot = DDrot(DDName(phiRotstr, idNameSpace), theta, phix, theta, phiy, 0., 0.);
       }
       phiRotMatrix = phiRot.matrix();
     }
 
     // globalRot def
     string globalRotstr = rotstr + "Phi" + to_string(phideg * 10.) + "Tilt" + to_string(convertRadToDeg(tiltAngle));
     if (isZPlus) {
       globalRotstr += "ZPlus";
       if (isFlipped) {
         globalRotstr += "Flip";
       }
     } else {
       globalRotstr += "ZMinus";
       if (isFlipped) {
         globalRotstr += "Flip";
       }
     }
     globalRot = DDRotation(DDName(globalRotstr, idNameSpace));
     if (!globalRot) {
       LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new "
                               << "rotation: " << globalRotstr;
       globalRotMatrix = phiOwnAxisRotMatrix * phiRotMatrix * tiltMatrix;
       globalRot = DDrot(DDName(globalRotstr, idNameSpace), make_unique<DDRotationMatrix>(globalRotMatrix));
     }
 
     // translation def
     double xpos = radius * cos(phix) + center[0];
     double ypos = radius * sin(phix) + center[1];
     double zpos = center[2];
     DDTranslation tran(xpos, ypos, zpos);
 
     // Positions child with respect to parent
     cpv.position(child, mother, copy, tran, globalRot);
     LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo test " << child << " number " << copy << " positioned in "
                             << mother << " at " << tran << " with " << globalRot;
 
     copy += incrCopyNo;
   }
 }

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

override

Definition at line 93 of file DDTrackerIrregularRingAlgo.cc.

References angle_units::operators::convertRadToDeg(), LogDebug, dqmiodumpmetadata::n, DDCurrentNamespace::ns(), SpecificationBuilder_cfi::parent(), and CosmicsPD_Skims::radius.

                                                                             {
   n = int(nArgs["N"]);
   startCopyNo = int(nArgs["StartCopyNo"]);
   incrCopyNo = int(nArgs["IncrCopyNo"]);
   rangeAngle = nArgs["RangeAngle"];
   startAngle = nArgs["StartAngle"];
   radius = nArgs["Radius"];
   center = vArgs["Center"];
   yawAngles = vArgs["yawAngleValues"];
   phiAngles = vArgs["phiAngleValues"];
   radiusValues = vArgs["radiusValues"];
   isZPlus = bool(nArgs["IsZPlus"]);
   tiltAngle = nArgs["TiltAngle"];
   isFlipped = bool(nArgs["IsFlipped"]);
 
   LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo debug: Parameters for position"
                           << "ing:: n " << n << " Start, Range " << convertRadToDeg(startAngle) << " "
                           << convertRadToDeg(rangeAngle) << " Radius " << radius << " Centre " << center[0] << ", "
                           << center[1] << ", " << center[2];
 
   idNameSpace = DDCurrentNamespace::ns();
   childName = sArgs["ChildName"];
 
   DDName parentName = parent().name();
   LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo debug: Parent " << parentName << "\tChild " << childName
                           << " NameSpace " << idNameSpace;
 }