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 25 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 33 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 38 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 44 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 61 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 77 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 93 of file GlobalErrorBase.h.

: theCartesianError(err) {}

template<class T, class ErrorWeightType>

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

inline

Definition at line 95 of file GlobalErrorBase.h.

{}

Member Function Documentation

template<class T, class ErrorWeightType>

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

inline

Definition at line 115 of file GlobalErrorBase.h.

{ return theCartesianError(3, 3); }

template<class T, class ErrorWeightType>

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

inline

Definition at line 109 of file GlobalErrorBase.h.

{ return theCartesianError(3, 0); }

template<class T, class ErrorWeightType>

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

inline

Definition at line 111 of file GlobalErrorBase.h.

{ return theCartesianError(3, 1); }

template<class T, class ErrorWeightType>

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

inline

Definition at line 113 of file GlobalErrorBase.h.

{ return theCartesianError(3, 2); }

template<class T, class ErrorWeightType>

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

inline

Definition at line 97 of file GlobalErrorBase.h.

Referenced by AlignableModifier::addAlignmentPositionErrorLocal(), AlignmentPositionError::AlignmentPositionError(), CkfDebugger::analyseCompatibleMeasurements(), TestHits::analyze(), TestSmoothHits::analyze(), TestTrackHits::analyze(), btagbtvdeep::TrackPairInfoBuilder::buildTrackPairInfo(), 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 99 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 101 of file GlobalErrorBase.h.

Referenced by AlignableModifier::addAlignmentPositionErrorLocal(), AlignmentPositionError::AlignmentPositionError(), CkfDebugger::analyseCompatibleMeasurements(), TestHits::analyze(), TestSmoothHits::analyze(), TestTrackHits::analyze(), btagbtvdeep::TrackPairInfoBuilder::buildTrackPairInfo(), 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 103 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 105 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 107 of file GlobalErrorBase.h.

Referenced by AlignableModifier::addAlignmentPositionErrorLocal(), AlignmentPositionError::AlignmentPositionError(), CkfDebugger::analyseCompatibleMeasurements(), TestHits::analyze(), TestSmoothHits::analyze(), TestTrackHits::analyze(), btagbtvdeep::TrackPairInfoBuilder::buildTrackPairInfo(), 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 121 of file GlobalErrorBase.h.

Referenced by BPHCompositeVertexSelect::accept(), BPHFittedVertexSelect::accept(), GsfVertexWeightCalculator::calculate(), VertexDistanceXY::compatibility(), VertexDistance3D::compatibility(), BasicSingleVertexState::computeError(), BasicSingleVertexState::computeWeight(), VertexDistanceXY::distance(), VertexDistance3D::distance(), BasicSingleVertexState::error(), HIPTwoBodyDecayAnalyzer::fitDimuonVertex(), CocoaDBMgr::GetAlignInfoErrorFromOptO(), reco::GhostTrackPrediction::init(), reco::GhostTrackState::lambdaError(), KinematicVertex::operator reco::Vertex(), KalmanTrackToTrackCovCalculator< N >::operator()(), 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 136 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 160 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 163 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 149 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 138 of file GlobalErrorBase.h.

Referenced by BPHMonitor::analyze(), LeptonSkimming::filter(), TValidTrackingRecHit::globalState(), BaseTrackerRecHit::globalState(), CAHitTripletGenerator::hitNtuplets(), CAHitQuadrupletGenerator::hitNtuplets(), HLTmmkFilter::hltFilter(), HLTmmkkFilter::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 168 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-().