An element model implementing the deformations used for the NSW sTGC as-built parameters. More...

#include <ElementModelSTGC.h>

Inheritance diagram for NswAsBuilt::ElementModelSTGC:

Collaboration diagram for NswAsBuilt::ElementModelSTGC:

Public Types
enum	parameter_t { X =0 , Y =1 , Z =2 , THX =3 , THY =4 , THZ =5 , ROT =6 , OFF =7 , SCL =8 , NPAR =9 }
using	ipar_t = unsigned int
using	VectorSet = Eigen::Matrix<double, 3, Eigen::Dynamic, Eigen::ColMajor\|Eigen::AutoAlign, 3, 5>
using	VectorSetRef = Eigen::Ref<VectorSet>

Public Member Functions
	ElementModelSTGC (double lenX, double lenY, Amg::Vector3D defo0)
	ElementModelSTGC ()=delete
virtual int	nParameters () const override
	The number of parameters used.
virtual void	transform (const ParameterVector &parvec, VectorSetRef local) const override
	Transform a set of vectors expressed in local frame, stored in a matrix.
virtual void	cacheTransform (ParameterVector &parvec) const override
	Cache the rigid component of this deformation model.
virtual ipar_t	getParameterIndex (const std::string &parname) const override
virtual std::string	getParameterName (ipar_t ipar) const override

Private Member Functions
void	applyDeformation (const ParameterVector &parvec, Eigen::Ref< Amg::Vector3D > local) const
void	applyDeformation2 (const ParameterVector &parvec, VectorSetRef local) const
Amg::Vector3D	stgcScale (double scl, const Amg::Vector3D &d0) const
Amg::Vector3D	stgcNonPar (double npar, const Amg::Vector3D &d0) const

Static Private Member Functions
static Amg::Vector3D	stgcOffset (double off)
static Amg::Vector3D	stgcRotation (double rot, const Amg::Vector3D &d0)

Private Attributes
double	m_lenY {0.}
double	m_lenX {0.}
Amg::Vector3D	m_defo0 {Amg::Vector3D::Zero()}

Detailed Description

An element model implementing the deformations used for the NSW sTGC as-built parameters.

The first 6 parameters represent the rigid component of the transform (a 3D translation and 3 Euler angles to be applied in the order ZYX)

The following parameters are providing deformations, in accordance with the: 4-parameter alingment model implemented for the sTGC quadruplets:

dy = off + rot*x + scl*y + npar*x*y

For more information check https://cds.cern.ch/record/2688394?ln=en

Definition at line 26 of file ElementModelSTGC.h.

Member Typedef Documentation

◆ ipar_t

using NswAsBuilt::ElementModel::ipar_t = unsigned int

inherited

Definition at line 35 of file ElementModel.h.

◆ VectorSet

using NswAsBuilt::ElementModel::VectorSet = Eigen::Matrix<double, 3, Eigen::Dynamic, Eigen::ColMajor|Eigen::AutoAlign, 3, 5>

inherited

Definition at line 36 of file ElementModel.h.

◆ VectorSetRef

using NswAsBuilt::ElementModel::VectorSetRef = Eigen::Ref<VectorSet>

inherited

Definition at line 37 of file ElementModel.h.

Member Enumeration Documentation

◆ parameter_t

enum NswAsBuilt::ElementModelSTGC::parameter_t

Enumerator
X
Y
Z
THX
THY
THZ
ROT
OFF
SCL
NPAR

Definition at line 28 of file ElementModelSTGC.h.

                       {
        X=0,
        Y=1,
        Z=2,
        THX=3,
        THY=4,
        THZ=5,
        ROT=6,
        OFF=7,
        SCL=8,
        NPAR=9
      };

Constructor & Destructor Documentation

◆ ElementModelSTGC() [1/2]

ElementModelSTGC::ElementModelSTGC	(	double	lenX,
		double	lenY,
		Amg::Vector3D	defo0 )

Definition at line 11 of file ElementModelSTGC.cxx.

  : m_lenY(lenY),m_lenX(lenX),
  m_defo0(std::move(defo0))
{
}

◆ ElementModelSTGC() [2/2]

NswAsBuilt::ElementModelSTGC::ElementModelSTGC ( )

delete

Member Function Documentation

◆ applyDeformation()

void ElementModelSTGC::applyDeformation	(	const ParameterVector &	parvec,
		Eigen::Ref< Amg::Vector3D >	local ) const

private

Definition at line 110 of file ElementModelSTGC.cxx.

                                                                                                        {
  // Applies the deformation to the set of local vectors given as argument
  // This old-style implementation is the reference implementation
  Amg::Vector3D d0 = local - m_defo0;
  local = local
    + stgcOffset(parvec[OFF])
    + stgcRotation(parvec[ROT], d0)
    + stgcScale(parvec[SCL], d0)
    + stgcNonPar(parvec[NPAR], d0)
    ;
}

◆ applyDeformation2()

void ElementModelSTGC::applyDeformation2	(	const ParameterVector &	parvec,
		VectorSetRef	local ) const

private

Definition at line 122 of file ElementModelSTGC.cxx.

                                                                                                {
  // Applies the deformation to the set of local vectors given as argument
  // This implementation uses Eigen-style algebra (does the same as the method applyDeformation above, but benefits from Eigen's optimizations)
 
  // Temporaries allocated on the stack
  // d0 = local - defo0
  VectorSet d0 = local.colwise() - m_defo0;
 
  double off = parvec[OFF];
  double rot = parvec[ROT];
  double scl = parvec[SCL]/m_lenY;
  double npar = parvec[NPAR]/(m_lenX*m_lenY);
 
  // OFF:
  local.array().colwise() += Eigen::Array3d{0. ,off, 0.};
 
  // ROT
  local.topRows<2>().array() += d0.topRows<2>().array().colwise().reverse().colwise() * Eigen::Array2d{0., rot};
 
  // SCL:
  local.array() += d0.array().colwise() * Eigen::Array3d{0., scl, 0.};
 
  // NPAR:
  local.topRows<2>().array() +=  (d0.topRows<2>().array().colwise().reverse() * d0.topRows<2>().array()).colwise() * Eigen::Array2d{0., npar};
}

◆ cacheTransform()

void ElementModelSTGC::cacheTransform ( ParameterVector & parvec ) const

overridevirtual

Cache the rigid component of this deformation model.

Cache the rigid component of this deformation model Add rotation and offset obtained from the sTGC fit.

Implements NswAsBuilt::ElementModel.

Definition at line 45 of file ElementModelSTGC.cxx.

                                                                   {
  const auto& pars = parvec.parameters;
  parvec.transformCache
    = Eigen::Translation3d(pars[X], pars[Y], pars[Z])
    * Eigen::AngleAxisd(pars[THZ], Eigen::Vector3d::UnitZ())
    * Eigen::AngleAxisd(pars[THY], Eigen::Vector3d::UnitY())
    * Eigen::AngleAxisd(pars[THX], Eigen::Vector3d::UnitX());
  parvec.transformCacheValid = true;
}

◆ getParameterIndex()

ElementModelSTGC::ipar_t ElementModelSTGC::getParameterIndex ( const std::string & parname ) const

overridevirtual

Implements NswAsBuilt::ElementModel.

Definition at line 58 of file ElementModelSTGC.cxx.

                                                                                         {
  if (parname == "x") return parameter_t::X;
  if (parname == "y") return parameter_t::Y;
  if (parname == "z") return parameter_t::Z;
  if (parname == "thx") return parameter_t::THX;
  if (parname == "thy") return parameter_t::THY;
  if (parname == "thz") return parameter_t::THZ;
  if (parname == "rot") return parameter_t::ROT;
  if (parname == "off") return parameter_t::OFF;
  if (parname == "scl") return parameter_t::SCL;
  if (parname == "npar") return parameter_t::NPAR;
  throw std::runtime_error("Invalid parameter name "+parname);
}

◆ getParameterName()

std::string ElementModelSTGC::getParameterName ( ipar_t ipar ) const

overridevirtual

Implements NswAsBuilt::ElementModel.

Definition at line 72 of file ElementModelSTGC.cxx.

                                                              {
  switch (ipar) {
    case X: return "x";
    case Y: return "y";
    case Z: return "z";
    case THX: return "thx";
    case THY: return "thy";
    case THZ: return "thz";
    case ROT: return "rot";
    case OFF: return "off";
    case SCL: return "scl";
    case NPAR: return "npar";
    default: throw std::runtime_error("Invalid parameter");
  }
}

◆ nParameters()

virtual int NswAsBuilt::ElementModelSTGC::nParameters ( ) const

inlineoverridevirtual

The number of parameters used.

Implements NswAsBuilt::ElementModel.

Definition at line 47 of file ElementModelSTGC.h.

47{ return 10; }

◆ stgcNonPar()

Amg::Vector3D ElementModelSTGC::stgcNonPar	(	double	npar,
		const Amg::Vector3D &	d0 ) const

private

Definition at line 104 of file ElementModelSTGC.cxx.

                                                                                 {
// Non-parallelism measured by CMM/Faro at construction sites
  double delta = npar*d0[0]*d0[1]/(m_lenX*m_lenY);
  return Amg::Vector3D(0., delta, 0.);
}

◆ stgcOffset()

Amg::Vector3D ElementModelSTGC::stgcOffset ( double off )

staticprivate

Definition at line 88 of file ElementModelSTGC.cxx.

                                                   {
// Offset extracted from combined fit of X-ray data and by CMM/Faro measurement at construction sites
  return off * Amg::Vector3D::UnitY();
}

◆ stgcRotation()

Amg::Vector3D ElementModelSTGC::stgcRotation	(	double	rot,
		const Amg::Vector3D &	d0 )

staticprivate

Definition at line 93 of file ElementModelSTGC.cxx.

                                                                            {
// Rotation extracted from combined fit of X-ray data and by CMM/Faro measurement at construction sites
  return -rot * d0.cross(Amg::Vector3D::UnitZ());
}

◆ stgcScale()

Amg::Vector3D ElementModelSTGC::stgcScale	(	double	scl,
		const Amg::Vector3D &	d0 ) const

private

Definition at line 98 of file ElementModelSTGC.cxx.

                                                                               {
// Scale measured by CMM/Faro at construction sites
  return d0.array() * Eigen::Array3d{0., scl/m_lenY, 0.};
}

◆ transform()

void ElementModelSTGC::transform	(	const ParameterVector &	parvec,
		VectorSetRef	local ) const

overridevirtual

Transform a set of vectors expressed in local frame, stored in a matrix.

Implements NswAsBuilt::ElementModel.

Definition at line 20 of file ElementModelSTGC.cxx.

                                                                                        {
 
  if (!parvec.transformCacheValid) {
    throw std::runtime_error("Should call Element::cacheTransforms() first");
  }
 
  // Apply the deformation component
  
  // old-style (reference) implementation: not optimized, provided for comparison purposes
  // for (int i=0; i< local.cols(); ++i) {
  //   applyDeformation(parvec, local.col(i));
  // }
 
 
  // Eigen-style implementation: optimized
  applyDeformation2(parvec, local);
 
  // Apply the rigid component
  local = (parvec.transformCache * local).eval(); // Needs eval to avoid aliasing (?)
}

Member Data Documentation

◆ m_defo0

Amg::Vector3D NswAsBuilt::ElementModelSTGC::m_defo0 {Amg::Vector3D::Zero()}

private

Definition at line 83 of file ElementModelSTGC.h.

83{Amg::Vector3D::Zero()};

◆ m_lenX

double NswAsBuilt::ElementModelSTGC::m_lenX {0.}

private

Definition at line 82 of file ElementModelSTGC.h.

82{0.};

◆ m_lenY

double NswAsBuilt::ElementModelSTGC::m_lenY {0.}

private

Definition at line 81 of file ElementModelSTGC.h.

81{0.};

The documentation for this class was generated from the following files:

Public Types

Public Member Functions

Private Member Functions

Static Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ ipar_t

◆ VectorSet

◆ VectorSetRef

Member Enumeration Documentation

◆ parameter_t

Constructor & Destructor Documentation

◆ ElementModelSTGC() [1/2]

◆ ElementModelSTGC() [2/2]

Member Function Documentation

◆ applyDeformation()

◆ applyDeformation2()

◆ cacheTransform()

◆ getParameterIndex()

◆ getParameterName()

◆ nParameters()

◆ stgcNonPar()

◆ stgcOffset()

◆ stgcRotation()

◆ stgcScale()

◆ transform()

Member Data Documentation

◆ m_defo0

◆ m_lenX

◆ m_lenY