GlobalErrorBase< T, ErrorWeightType > Class Template Reference

#include <GlobalErrorBase.h>

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

Public Member Functions
T	ctt () const
T	ctx () const
T	cty () const
T	ctz () const
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 (T c11, T c21, T c22, T c31, T c32, T c33, T c41, T c42, T c43, T c44)
	GlobalErrorBase (const AlgebraicSymMatrix33 &err)
	GlobalErrorBase (const AlgebraicSymMatrix44 &err)
const AlgebraicSymMatrix33	matrix () const
const AlgebraicSymMatrix44 &	matrix4D () 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
AlgebraicSymMatrix44	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;
    theCartesianError(3,0) = 0.;
    theCartesianError(3,1) = 0.;
    theCartesianError(3,2) = 0.;
    theCartesianError(3,3) = 0.;
  }

template<class T, class ErrorWeightType>

GlobalErrorBase< T, ErrorWeightType >::GlobalErrorBase ( T  c11, T  c21, T  c22, T  c31, T  c32, T  c33, T  c41, T  c42, T  c43, T  c44  )

inline

Constructor. The symmetric matrix stored as a lower triangular matrix (4D)

Definition at line 65 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;
    theCartesianError(3,0)=c41;
    theCartesianError(3,1)=c42;
    theCartesianError(3,2)=c43;
    theCartesianError(3,3)=c44;
  }

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 81 of file GlobalErrorBase.h.

                                                   { 
     theCartesianError(0,0) = err(0,0);
     theCartesianError(1,0) = err(1,0);
     theCartesianError(1,1) = err(1,1);
     theCartesianError(2,0) = err(2,0);
     theCartesianError(2,1) = err(2,1);
     theCartesianError(2,2) = err(2,2);
     theCartesianError(3,0) = 0.;
     theCartesianError(3,1) = 0.;
     theCartesianError(3,2) = 0.;
     theCartesianError(3,3) = 0.;
  }

template<class T, class ErrorWeightType>

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

inline

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

Definition at line 97 of file GlobalErrorBase.h.

                                                   : 
  theCartesianError(err) { }

template<class T, class ErrorWeightType>

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

inline

Definition at line 101 of file GlobalErrorBase.h.

{}

Member Function Documentation

template<class T, class ErrorWeightType>

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

inline

Definition at line 139 of file GlobalErrorBase.h.

                {
    return theCartesianError(3,3);
  }

template<class T, class ErrorWeightType>

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

inline

Definition at line 127 of file GlobalErrorBase.h.

                {
    return theCartesianError(3,0);
  }

template<class T, class ErrorWeightType>

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

inline

Definition at line 131 of file GlobalErrorBase.h.

                {
    return theCartesianError(3,1);
  }

template<class T, class ErrorWeightType>

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

inline

Definition at line 135 of file GlobalErrorBase.h.

                {
    return theCartesianError(3,2);
  }

template<class T, class ErrorWeightType>

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

inline

Definition at line 103 of file GlobalErrorBase.h.

Referenced by AlignableModifier::addAlignmentPositionErrorLocal(), AlignmentPositionError::AlignmentPositionError(), CkfDebugger::analyseCompatibleMeasurements(), TestHits::analyze(), TestSmoothHits::analyze(), TestTrackHits::analyze(), RecoVertex::convertError(), VertexFitterResult::fill(), SeedForPhotonConversion1Leg::initialError(), SeedForPhotonConversionFromQuadruplets::initialError(), GlobalErrorBase< double, ErrorMatrixTag >::phierr(), PrimaryVertexProducer::produce(), GlobalErrorBase< double, ErrorMatrixTag >::rerr(), KFBasedPixelFitter::run(), CandKinematicVertexFitter::set(), ErrorFrameTransformer::transform(), and PrimaryVertexProducerAlgorithm::vertices().

                {
    return theCartesianError(0,0);
  }

template<class T, class ErrorWeightType>

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

inline

Definition at line 107 of file GlobalErrorBase.h.

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

                {
    return theCartesianError(1,0);
  }

template<class T, class ErrorWeightType>

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

inline

Definition at line 111 of file GlobalErrorBase.h.

Referenced by AlignableModifier::addAlignmentPositionErrorLocal(), AlignmentPositionError::AlignmentPositionError(), CkfDebugger::analyseCompatibleMeasurements(), TestHits::analyze(), TestSmoothHits::analyze(), TestTrackHits::analyze(), RecoVertex::convertError(), VertexFitterResult::fill(), GlobalErrorBase< double, ErrorMatrixTag >::phierr(), PrimaryVertexProducer::produce(), GlobalErrorBase< double, ErrorMatrixTag >::rerr(), CandKinematicVertexFitter::set(), ErrorFrameTransformer::transform(), and PrimaryVertexProducerAlgorithm::vertices().

                {
    return theCartesianError(1,1);
  }

template<class T, class ErrorWeightType>

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

inline

Definition at line 115 of file GlobalErrorBase.h.

Referenced by AlignableModifier::addAlignmentPositionErrorLocal(), AlignmentPositionError::AlignmentPositionError(), RecoVertex::convertError(), CandKinematicVertexFitter::set(), and ErrorFrameTransformer::transform().

                {
    return theCartesianError(2,0);
  }

template<class T, class ErrorWeightType>

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

inline

Definition at line 119 of file GlobalErrorBase.h.

Referenced by AlignableModifier::addAlignmentPositionErrorLocal(), AlignmentPositionError::AlignmentPositionError(), RecoVertex::convertError(), CandKinematicVertexFitter::set(), and ErrorFrameTransformer::transform().

                {
    return theCartesianError(2,1);
  }

template<class T, class ErrorWeightType>

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

inline

Definition at line 123 of file GlobalErrorBase.h.

Referenced by AlignableModifier::addAlignmentPositionErrorLocal(), AlignmentPositionError::AlignmentPositionError(), CkfDebugger::analyseCompatibleMeasurements(), TestHits::analyze(), TestSmoothHits::analyze(), TestTrackHits::analyze(), BarrelDetLayer::compatible(), RecoVertex::convertError(), VertexFitterResult::fill(), TValidTrackingRecHit::globalState(), BaseTrackerRecHit::globalState(), CAHitTripletGenerator::hitNtuplets(), CAHitQuadrupletGenerator::hitNtuplets(), SeedForPhotonConversion1Leg::initialError(), SeedForPhotonConversionFromQuadruplets::initialError(), PrimaryVertexProducer::produce(), PixelFitterByConformalMappingAndLine::run(), KFBasedPixelFitter::run(), CandKinematicVertexFitter::set(), ErrorFrameTransformer::transform(), PrimaryVertexProducerAlgorithm::vertices(), and SimpleNavigableLayer::wellInside().

                {
    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 147 of file GlobalErrorBase.h.

Referenced by BPHCompositeVertexSelect::accept(), BPHFittedVertexSelect::accept(), GsfVertexWeightCalculator::calculate(), VertexDistance3D::compatibility(), VertexDistanceXY::compatibility(), BasicSingleVertexState::computeError(), BasicSingleVertexState::computeWeight(), VertexDistance3D::distance(), VertexDistanceXY::distance(), BasicSingleVertexState::error(), HIPTwoBodyDecayAnalyzer::fitDimuonVertex(), CocoaDBMgr::GetAlignInfoErrorFromOptO(), reco::GhostTrackPrediction::init(), reco::GhostTrackState::lambdaError(), KinematicVertex::operator reco::Vertex(), TemplatedSecondaryVertexProducer< IPTI, VTX >::SVBuilder::operator()(), KinematicVertex::operator==(), KalmanVertexUpdator< N >::positionUpdate(), Onia2MuMuPAT::produce(), PrimaryVertexProducer::produce(), KalmanVertexTrackUpdator< N >::trackRefit(), trackVertexCompat(), vertexAtState(), KVFHelper< N >::vertexChi2(), reco::GhostTrackVertexFinder::vertexCompat(), PrimaryVertexProducerAlgorithm::vertices(), vtxErrorLong(), and BasicSingleVertexState::weight().

                                            {
    AlgebraicSymMatrix33 out;
    out(0,0) = theCartesianError(0,0);
    out(1,0) = theCartesianError(1,0);
    out(1,1) = theCartesianError(1,1);
    out(2,0) = theCartesianError(2,0);
    out(2,1) = theCartesianError(2,1);
    out(2,2) = theCartesianError(2,2);
    return out;
  }  

template<class T, class ErrorWeightType>

const AlgebraicSymMatrix44& GlobalErrorBase< T, ErrorWeightType >::matrix4D ( ) const

inline

Access method to the matrix, /return The SymMatrix 4x4

Definition at line 162 of file GlobalErrorBase.h.

Referenced by BasicSingleVertexState::computeError(), BasicSingleVertexState::computePosition(), BasicSingleVertexState::computeWeight(), BasicSingleVertexState::computeWeightTimesPos(), TransientVertex::operator reco::Vertex(), PrimaryVertexProducer::produce(), BasicSingleVertexState::timeError(), and svhelper::updateVertexTime().

                                               {
    return theCartesianError;
  }

template<class T, class ErrorWeightType>

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

inline

Definition at line 188 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 191 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 177 of file GlobalErrorBase.h.

Referenced by TValidTrackingRecHit::globalState(), BaseTrackerRecHit::globalState(), and GlobalTrackQualityProducer::kink().

                                            {
    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 166 of file GlobalErrorBase.h.

Referenced by BPHMonitor::analyze(), BPHMonitor::case11_selection(), LeptonSkimming::filter(), TValidTrackingRecHit::globalState(), BaseTrackerRecHit::globalState(), CAHitTripletGenerator::hitNtuplets(), CAHitQuadrupletGenerator::hitNtuplets(), HLTmmkkFilter::hltFilter(), HLTmmkFilter::hltFilter(), and PixelFitterByConformalMappingAndLine::run().

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

AlgebraicSymMatrix44 GlobalErrorBase< T, ErrorWeightType >::theCartesianError

private

Definition at line 197 of file GlobalErrorBase.h.

Referenced by GlobalErrorBase< double, ErrorMatrixTag >::ctt(), GlobalErrorBase< double, ErrorMatrixTag >::ctx(), GlobalErrorBase< double, ErrorMatrixTag >::cty(), GlobalErrorBase< double, ErrorMatrixTag >::ctz(), 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 >::matrix4D(), GlobalErrorBase< double, ErrorMatrixTag >::operator+(), and GlobalErrorBase< double, ErrorMatrixTag >::operator-().