#include <BowedSurfaceDeformation.h>

Inheritance diagram for BowedSurfaceDeformation:

Public Member Functions
bool	add (const SurfaceDeformation &other) override

	BowedSurfaceDeformation (const std::vector< double > &parameters)

	BowedSurfaceDeformation (double sagittaX, double sagittaXY, double sagittaY)
	constructor from sagittae, i.e. coefficients of Legendre polynomials More...

BowedSurfaceDeformation *	clone () const override

std::vector< double >	parameters () const override
	parameters, i.e. sagittae as given in the constructor More...

Local2DVector	positionCorrection (const Local2DPoint &localPos, const LocalTrackAngles &localAngles, double length, double width) const override

int	type () const override
	specific type, i.e. SurfaceDeformationFactory::kBowedSurface More...

Public Member Functions inherited from SurfaceDeformation
virtual	~SurfaceDeformation ()

Static Public Member Functions
static constexpr unsigned int	maxParameterSize ()
	maximum size of vector that is accepted by constructor from vector More...

static constexpr unsigned int	minParameterSize ()
	minimum size of vector that is accepted by constructor from vector More...

static constexpr unsigned int	parameterSize ()

Static Public Attributes
static constexpr unsigned int	parSize = 3

Private Attributes
double	theSagittaX

double	theSagittaXY

double	theSagittaY

Additional Inherited Members
Public Types inherited from SurfaceDeformation
typedef Topology::Local2DPoint	Local2DPoint

typedef Vector2DBase< double, LocalTag >	Local2DVector

typedef Topology::LocalTrackAngles	LocalTrackAngles

typedef Topology::MathVector2D	MathVector2D

typedef Topology::Vector2D	Vector2D

Detailed Description

BowedSurfaceAlignmentParameters

Class to apply corrections to local positions resulting from a non-planar surface. The bows are parametrised using Legendre polynomials up to second order, excluding 0th and 1st order that are already treated by local w shift and rotations around local u and v axes.

Author: : Gero Flucke date : October 2010

Definition at line 20 of file BowedSurfaceDeformation.h.

Constructor & Destructor Documentation

◆ BowedSurfaceDeformation() [1/2]

BowedSurfaceDeformation::BowedSurfaceDeformation	(	double	sagittaX,
		double	sagittaXY,
		double	sagittaY
	)

inline

constructor from sagittae, i.e. coefficients of Legendre polynomials

Definition at line 23 of file BowedSurfaceDeformation.h.

24 : theSagittaX(sagittaX), theSagittaY(sagittaY), theSagittaXY(sagittaXY) {}

Referenced by clone().

◆ BowedSurfaceDeformation() [2/2]

BowedSurfaceDeformation::BowedSurfaceDeformation ( const std::vector< double > & pars )

constructor from vector of sagittae, parameters.size() must be between minParameterSize() and maxParameterSize()

Author: : Gero Flucke date : October 2010

Definition at line 13 of file BowedSurfaceDeformation.cc.

     : theSagittaX(!pars.empty() ? pars[0] : 0.),
       theSagittaY(pars.size() > 2 ? pars[2] : 0.),
       theSagittaXY(pars.size() > 1 ? pars[1] : 0.) {
   if (pars.size() != minParameterSize()) {
     edm::LogError("BadSetup") << "@SUB=BowedSurfaceDeformation"
                               << "Input vector of wrong size " << pars.size() << " instead of " << minParameterSize()
                               << ", filled up with zeros!";
   }
 }

References minParameterSize().

Member Function Documentation

◆ add()

bool BowedSurfaceDeformation::add ( const SurfaceDeformation & other )

overridevirtual

update information with parameters of 'other', false in case the type or some parameters do not match and thus the information cannot be used (then no changes are done), true if merge was successful

Implements SurfaceDeformation.

Definition at line 71 of file BowedSurfaceDeformation.cc.

                                                                  {
   if (other.type() == this->type()) {
     const std::vector<double> otherParams(other.parameters());
     if (otherParams.size() == parameterSize()) {  // double check!
       theSagittaX += otherParams[0];              // bows can simply be added up
       theSagittaXY += otherParams[1];
       theSagittaY += otherParams[2];
  
       return true;
     }
   }
  
   return false;
 }

References trackingPlots::other, parameterSize(), theSagittaX, theSagittaXY, theSagittaY, and type().

Referenced by counter.Counter::register().

◆ clone()

BowedSurfaceDeformation * BowedSurfaceDeformation::clone ( void ) const

overridevirtual

Implements SurfaceDeformation.

Definition at line 25 of file BowedSurfaceDeformation.cc.

                                                               {
   return new BowedSurfaceDeformation(theSagittaX, theSagittaXY, theSagittaY);
 }

References BowedSurfaceDeformation(), theSagittaX, theSagittaXY, and theSagittaY.

◆ maxParameterSize()

static constexpr unsigned int BowedSurfaceDeformation::maxParameterSize ( )

inlinestaticconstexpr

maximum size of vector that is accepted by constructor from vector

Definition at line 61 of file BowedSurfaceDeformation.h.

61 { return parameterSize(); }

References parameterSize().

Referenced by SurfaceDeformationFactory::create().

◆ minParameterSize()

static constexpr unsigned int BowedSurfaceDeformation::minParameterSize ( )

inlinestaticconstexpr

minimum size of vector that is accepted by constructor from vector

Definition at line 59 of file BowedSurfaceDeformation.h.

59 { return parameterSize(); }

References parameterSize().

Referenced by BowedSurfaceDeformation(), and SurfaceDeformationFactory::create().

◆ parameters()

std::vector< double > BowedSurfaceDeformation::parameters ( void ) const

overridevirtual

parameters, i.e. sagittae as given in the constructor

Implements SurfaceDeformation.

Definition at line 87 of file BowedSurfaceDeformation.cc.

                                                             {
   std::vector<double> result(parameterSize());
   result[0] = theSagittaX;
   result[1] = theSagittaXY;
   result[2] = theSagittaY;
  
   return result;
 }

References parameterSize(), mps_fire::result, theSagittaX, theSagittaXY, and theSagittaY.

◆ parameterSize()

static constexpr unsigned int BowedSurfaceDeformation::parameterSize ( )

inlinestaticconstexpr

Definition at line 56 of file BowedSurfaceDeformation.h.

56 { return parSize; }

References parSize.

Referenced by add(), maxParameterSize(), minParameterSize(), and parameters().

◆ positionCorrection()

SurfaceDeformation::Local2DVector BowedSurfaceDeformation::positionCorrection	(	const Local2DPoint &	localPos,
		const LocalTrackAngles &	localAngles,
		double	length,
		double	width
	)		const

overridevirtual

correction to add to local position depending on

track parameters in local frame (from LocalTrajectoryParameters):
- track position as Local2DPoint(x,y)
- track angles as LocalTrackAngles(dxdz, dydz)
length of surface (local y-coordinate)
width of surface (local x-coordinate)

Implements SurfaceDeformation.

Definition at line 33 of file BowedSurfaceDeformation.cc.

                                                                                                   {
   // different widthes at high/low y could somehow be treated by theRelWidthLowY
   //   if (widthLowY > 0. && widthHighY != widthLowY) {
   //     // TEC would always create a warning...
   //     edm::LogWarning("UnusableData") << "@SUB=BowedSurfaceDeformation::positionCorrection"
   //                    << "Cannot yet deal with different widthes, take "
   //                    << widthHighY << " not " << widthLowY;
   //   }
   //   const double width = widthHighY;
  
   double uRel = (width ? 2. * localPos.x() / width : 0.);    // relative u (-1 .. +1)
   double vRel = (length ? 2. * localPos.y() / length : 0.);  // relative v (-1 .. +1)
   // 'range check':
   const double cutOff = 1.5;
   if (uRel < -cutOff) {
     uRel = -cutOff;
   } else if (uRel > cutOff) {
     uRel = cutOff;
   }
   if (vRel < -cutOff) {
     vRel = -cutOff;
   } else if (vRel > cutOff) {
     vRel = cutOff;
   }
  
   // apply coefficients to Legendre polynomials
   // to get local height relative to 'average'
   const double dw =
       (uRel * uRel - 1. / 3.) * theSagittaX + uRel * vRel * theSagittaXY + (vRel * vRel - 1. / 3.) * theSagittaY;
  
   // positive dxdz/dydz and positive dw mean negative shift in x/y:
   return Local2DVector(-dw * localAngles);
 }