GlobalErrorBase< T, ErrorWeightType > Class Template Reference

#include <GlobalErrorBase.h>

Classes
class	NullMatrix
	Tag to request a null error matrix. More...

Public Member Functions
T	cxx () const
T	cyx () const
T	cyy () const
T	czx () const
T	czy () const
T	czz () const
	GlobalErrorBase ()
	GlobalErrorBase (const NullMatrix &)
	GlobalErrorBase (T c11, T c21, T c22, T c31, T c32, T c33)
	GlobalErrorBase (const AlgebraicSymMatrix33 &err)
const AlgebraicSymMatrix33 &	matrix () const
const AlgebraicSymMatrix33 &	matrix_new () const
GlobalErrorBase	operator+ (const GlobalErrorBase &err) const
GlobalErrorBase	operator- (const GlobalErrorBase &err) const
T	phierr (const GlobalPoint &aPoint) const
T	rerr (const GlobalPoint &aPoint) const
	~GlobalErrorBase ()

Private Attributes
AlgebraicSymMatrix33	theCartesianError

Detailed Description

template<class T, class ErrorWeightType>
class GlobalErrorBase< T, ErrorWeightType >

Templated class representing a symmetric 3*3 matrix describing, according to the ErrorWeightType tag, a (cartesian) covariance matrix or the weight matrix (the inverse of the covariance matrix).

• To have a covariance matrix, the ErrorMatrixTag has to be used, and a typedef is available as GlobalError
• To have a weight matrix, the WeightMatrixTag has to be used, and a typedef is available as Globalweight

The typedefs should be used in the code.

Definition at line 26 of file GlobalErrorBase.h.

Constructor & Destructor Documentation

template<class T, class ErrorWeightType>

GlobalErrorBase< T, ErrorWeightType >::GlobalErrorBase ( )

inline

Default constructor, creating a null 3*3 matrix (all values are 0)

Definition at line 36 of file GlobalErrorBase.h.

Referenced by GlobalErrorBase< double, ErrorMatrixTag >::operator+(), and GlobalErrorBase< double, ErrorMatrixTag >::operator-().

{}

template<class T, class ErrorWeightType>

GlobalErrorBase< T, ErrorWeightType >::GlobalErrorBase ( const NullMatrix &  )

inline

Obsolete Constructor that allocates a null GlobalErrorBase (it does not create the error matrix at all)

Definition at line 41 of file GlobalErrorBase.h.

{}

template<class T, class ErrorWeightType>

GlobalErrorBase< T, ErrorWeightType >::GlobalErrorBase ( T  c11, T  c21, T  c22, T  c31, T  c32, T  c33  )

inline

Constructor. The symmetric matrix stored as a lower triangular matrix

Definition at line 48 of file GlobalErrorBase.h.

                                                            {
    theCartesianError(0,0)=c11;
    theCartesianError(1,0)=c21;
    theCartesianError(1,1)=c22;
    theCartesianError(2,0)=c31;
    theCartesianError(2,1)=c32;
    theCartesianError(2,2)=c33;
  }

template<class T, class ErrorWeightType>

GlobalErrorBase< T, ErrorWeightType >::GlobalErrorBase ( const AlgebraicSymMatrix33 &  err )

inline

Constructor from SymMatrix. The original matrix has to be a 3*3 matrix.

Definition at line 60 of file GlobalErrorBase.h.

                                                      : 
      theCartesianError(err) { }

template<class T, class ErrorWeightType>

GlobalErrorBase< T, ErrorWeightType >::~GlobalErrorBase ( )

inline

Definition at line 63 of file GlobalErrorBase.h.

{}

Member Function Documentation

template<class T, class ErrorWeightType>

T GlobalErrorBase< T, ErrorWeightType >::cxx ( ) const

inline

Definition at line 65 of file GlobalErrorBase.h.

Referenced by AlignmentPositionError::AlignmentPositionError(), RecoVertex::convertError(), GlobalErrorBase< double, ErrorMatrixTag >::phierr(), GlobalErrorBase< double, ErrorMatrixTag >::rerr(), and ErrorFrameTransformer::transform().

                {
    return theCartesianError(0,0);
  }

template<class T, class ErrorWeightType>

T GlobalErrorBase< T, ErrorWeightType >::cyx ( ) const

inline

Definition at line 69 of file GlobalErrorBase.h.

Referenced by AlignmentPositionError::AlignmentPositionError(), RecoVertex::convertError(), GlobalErrorBase< double, ErrorMatrixTag >::phierr(), GlobalErrorBase< double, ErrorMatrixTag >::rerr(), and ErrorFrameTransformer::transform().

                {
    return theCartesianError(1,0);
  }

template<class T, class ErrorWeightType>

T GlobalErrorBase< T, ErrorWeightType >::cyy ( ) const

inline

Definition at line 73 of file GlobalErrorBase.h.

Referenced by AlignmentPositionError::AlignmentPositionError(), RecoVertex::convertError(), GlobalErrorBase< double, ErrorMatrixTag >::phierr(), GlobalErrorBase< double, ErrorMatrixTag >::rerr(), and ErrorFrameTransformer::transform().

                {
    return theCartesianError(1,1);
  }

template<class T, class ErrorWeightType>

T GlobalErrorBase< T, ErrorWeightType >::czx ( ) const

inline

Definition at line 77 of file GlobalErrorBase.h.

Referenced by AlignmentPositionError::AlignmentPositionError(), RecoVertex::convertError(), and ErrorFrameTransformer::transform().

                {
    return theCartesianError(2,0);
  }

template<class T, class ErrorWeightType>

T GlobalErrorBase< T, ErrorWeightType >::czy ( ) const

inline

Definition at line 81 of file GlobalErrorBase.h.

Referenced by AlignmentPositionError::AlignmentPositionError(), RecoVertex::convertError(), and ErrorFrameTransformer::transform().

                {
    return theCartesianError(2,1);
  }

template<class T, class ErrorWeightType>

T GlobalErrorBase< T, ErrorWeightType >::czz ( ) const

inline

Definition at line 85 of file GlobalErrorBase.h.

Referenced by AlignmentPositionError::AlignmentPositionError(), BarrelDetLayer::compatible(), RecoVertex::convertError(), TValidTrackingRecHit::globalState(), PixelFitterByConformalMappingAndLine::run(), and ErrorFrameTransformer::transform().

                {
    return theCartesianError(2,2);
  }

template<class T, class ErrorWeightType>

const AlgebraicSymMatrix33& GlobalErrorBase< T, ErrorWeightType >::matrix ( void  ) const

inline

Access method to the matrix, /return The SymMatrix

Definition at line 93 of file GlobalErrorBase.h.

Referenced by VertexDistance3D::compatibility(), VertexDistanceXY::compatibility(), BasicSingleVertexState::computeError(), BasicSingleVertexState::computeWeight(), VertexDistanceXY::distance(), VertexDistance3D::distance(), and KinematicVertex::operator==().

                                              {
    return theCartesianError;
  }

template<class T, class ErrorWeightType>

const AlgebraicSymMatrix33& GlobalErrorBase< T, ErrorWeightType >::matrix_new ( ) const

inline

Definition at line 96 of file GlobalErrorBase.h.

Referenced by GsfVertexWeightCalculator::calculate(), BasicSingleVertexState::computePosition(), BasicSingleVertexState::computeWeightTimesPos(), reco::GhostTrackState::lambdaError(), KalmanVertexTrackUpdator< N >::trackRefit(), and KVFHelper< N >::vertexChi2().

                                                  {
    return theCartesianError;
  }

template<class T, class ErrorWeightType>

GlobalErrorBase GlobalErrorBase< T, ErrorWeightType >::operator+ ( const GlobalErrorBase< T, ErrorWeightType > &  err ) const

inline

Definition at line 123 of file GlobalErrorBase.h.

                                                               {
    return GlobalErrorBase(theCartesianError + err.theCartesianError);
  }

template<class T, class ErrorWeightType>

GlobalErrorBase GlobalErrorBase< T, ErrorWeightType >::operator- ( const GlobalErrorBase< T, ErrorWeightType > &  err ) const

inline

Definition at line 126 of file GlobalErrorBase.h.

                                                               {
    return GlobalErrorBase(theCartesianError - err.theCartesianError);
  }

template<class T, class ErrorWeightType>

T GlobalErrorBase< T, ErrorWeightType >::phierr ( const GlobalPoint &  aPoint ) const

inline

Definition at line 112 of file GlobalErrorBase.h.

Referenced by TValidTrackingRecHit::globalState().

                                            {
    T r2 = aPoint.perp2();
    T x2 = aPoint.x()*aPoint.x();
    T y2 = aPoint.y()*aPoint.y();
    T xy = aPoint.x()*aPoint.y();
    if (r2 != 0) 
      return std::max<T>(0, (1./(r2*r2))*(y2*cxx() - 2.*xy*cyx() + x2*cyy()));
    else
      return 0;
  }

template<class T, class ErrorWeightType>

T GlobalErrorBase< T, ErrorWeightType >::rerr ( const GlobalPoint &  aPoint ) const

inline

Definition at line 101 of file GlobalErrorBase.h.

Referenced by TValidTrackingRecHit::globalState(), PixelFitterByConformalMappingAndLine::run(), and RZLine::RZLine().

                                          {
    T r2 = aPoint.perp2();
    T x2 = aPoint.x()*aPoint.x();
    T y2 = aPoint.y()*aPoint.y();
    T xy = aPoint.x()*aPoint.y();
    if(r2 != 0) 
      return std::max<T>(0, (1./r2)*(x2*cxx() + 2.*xy*cyx() + y2*cyy()));
    else 
      return 0.5*(cxx() + cyy());  
  }

Member Data Documentation

template<class T, class ErrorWeightType>

AlgebraicSymMatrix33 GlobalErrorBase< T, ErrorWeightType >::theCartesianError

private

Definition at line 132 of file GlobalErrorBase.h.

Referenced by GlobalErrorBase< double, ErrorMatrixTag >::cxx(), GlobalErrorBase< double, ErrorMatrixTag >::cyx(), GlobalErrorBase< double, ErrorMatrixTag >::cyy(), GlobalErrorBase< double, ErrorMatrixTag >::czx(), GlobalErrorBase< double, ErrorMatrixTag >::czy(), GlobalErrorBase< double, ErrorMatrixTag >::czz(), GlobalErrorBase< double, ErrorMatrixTag >::GlobalErrorBase(), GlobalErrorBase< double, ErrorMatrixTag >::matrix(), GlobalErrorBase< double, ErrorMatrixTag >::matrix_new(), GlobalErrorBase< double, ErrorMatrixTag >::operator+(), and GlobalErrorBase< double, ErrorMatrixTag >::operator-().