KVFHelper< N > Class Template Reference

#include <KVFHelper.h>

Public Types
typedef std::pair< bool, double >	BDpair
typedef VertexTrack< N >::RefCountedLinearizedTrackState	RefCountedLinearizedTrackState
typedef VertexTrack< N >::RefCountedRefittedTrackState	RefCountedRefittedTrackState
typedef CachingVertex< N >::RefCountedVertexTrack	RefCountedVertexTrack

Public Member Functions
BDpair	trackParameterChi2 (const RefCountedLinearizedTrackState linTrack, const RefCountedRefittedTrackState refittedTrackState) const
BDpair	trackParameterChi2 (const RefCountedVertexTrack track) const
double	vertexChi2 (const CachingVertex< N > &vertexA, const CachingVertex< N > &vertexB) const
double	vertexChi2 (const VertexState &vertexA, const VertexState &VertexB) const

Detailed Description

template<unsigned int N>
class KVFHelper< N >

Helper class to calculate chi2 increments for the Kalman Filter

Definition at line 11 of file KVFHelper.h.

Member Typedef Documentation

BDpair

template<unsigned int N>

typedef std::pair<bool, double> KVFHelper< N >::BDpair

Definition at line 16 of file KVFHelper.h.

RefCountedLinearizedTrackState

template<unsigned int N>

typedef VertexTrack<N>::RefCountedLinearizedTrackState KVFHelper< N >::RefCountedLinearizedTrackState

Definition at line 14 of file KVFHelper.h.

RefCountedRefittedTrackState

template<unsigned int N>

typedef VertexTrack<N>::RefCountedRefittedTrackState KVFHelper< N >::RefCountedRefittedTrackState

Definition at line 15 of file KVFHelper.h.

RefCountedVertexTrack

template<unsigned int N>

typedef CachingVertex<N>::RefCountedVertexTrack KVFHelper< N >::RefCountedVertexTrack

Definition at line 13 of file KVFHelper.h.

Member Function Documentation

trackParameterChi2() [1/2]

template<unsigned int N>

KVFHelper< N >::BDpair KVFHelper< N >::trackParameterChi2	(	const RefCountedLinearizedTrackState	linTrack,
		const RefCountedRefittedTrackState	refittedTrackState
	)		const

Methode which calculates the chi**2 between the prior and the fitted track parameters.

Parameters

linTrack	The track as linearized
refittedTrackState	The refitted track

Definition at line 37 of file KVFHelper.cc.

                                                                                                                {
  typedef ROOT::Math::SMatrix<double, N, N, ROOT::Math::MatRepSym<double, N> > AlgebraicSymMatrixNN;
  typedef ROOT::Math::SVector<double, N> AlgebraicVectorN;
 
  AlgebraicVectorN parameterResiduals =
      linTrack->predictedStateParameters() - linTrack->refittedParamFromEquation(refittedTrackState);
  linTrack->checkParameters(parameterResiduals);
  int error;
  float lChi2 = ROOT::Math::Similarity(parameterResiduals, linTrack->predictedStateWeight(error));
  if (error != 0)
    return BDpair(false, -1.);
  return BDpair(true, lChi2);
}

trackParameterChi2() [2/2]

template<unsigned int N>

KVFHelper< N >::BDpair KVFHelper< N >::trackParameterChi2

(

const RefCountedVertexTrack

track

)

const

Methode which calculates the chi**2 between the prior and the fitted track parameters. The track must have been smoothed before calling this method. If not, an Exception will be thrown.

Definition at line 32 of file KVFHelper.cc.

                                                                                                    {
  return trackParameterChi2(track->linearizedTrack(), track->refittedState());
}

vertexChi2() [1/2]

template<unsigned int N>

double KVFHelper< N >::vertexChi2	(	const CachingVertex< N > &	vertexA,
		const CachingVertex< N > &	vertexB
	)		const

Definition at line 5 of file KVFHelper.cc.

                                                                                                      {
  return vertexChi2(vertexA.vertexState(), vertexB.vertexState());
}

vertexChi2() [2/2]

template<unsigned int N>

double KVFHelper< N >::vertexChi2	(	const VertexState &	vertexA,
		const VertexState &	VertexB
	)		const

Methode which calculates the chi**2-increment due to the vertices E.g. between the prior and the fitted vertex. The covariance matrix used is that of the first vertex (vertexA). This method will not take into account multiple states, so in case one of the VertexStates is a multi-state vertex, only the mean will be used.

Parameters

vertexA	The prior vertex state
VertexB	The fitted vertex state

Definition at line 10 of file KVFHelper.cc.

                                                                                            {
  // std::cout <<"Start\n";
  GlobalPoint inPosition = vertexA.position();
  GlobalPoint fnPosition = vertexB.position();
  //   std::cout << inPosition<< fnPosition<<std::endl;
 
  AlgebraicVector3 oldVertexPositionV;
  oldVertexPositionV(0) = inPosition.x();
  oldVertexPositionV(1) = inPosition.y();
  oldVertexPositionV(2) = inPosition.z();
 
  AlgebraicVector3 newVertexPositionV;
  newVertexPositionV(0) = fnPosition.x();
  newVertexPositionV(1) = fnPosition.y();
  newVertexPositionV(2) = fnPosition.z();
 
  AlgebraicVector3 positionResidual = newVertexPositionV - oldVertexPositionV;
 
  return ROOT::Math::Similarity(positionResidual, vertexA.weight().matrix());
}