NswAsBuilt::ElementModelSTGC Class Reference

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. More...
virtual void	transform (const ParameterVector &parvec, VectorSetRef local) const override
	Transform a set of vectors expressed in local frame, stored in a matrix. More...
virtual void	cacheTransform (ParameterVector &parvec) const override
	Cache the rigid component of this deformation model. More...
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.

{ 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.

m_lenX

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

private

Definition at line 82 of file ElementModelSTGC.h.

m_lenY

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

private

Definition at line 81 of file ElementModelSTGC.h.

---

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

• ElementModelSTGC.h
• ElementModelSTGC.cxx