#include <KalmanVertexTrackCompatibilityEstimator.h>

Inheritance diagram for KalmanVertexTrackCompatibilityEstimator< N >:

Public Types
typedef std::pair< bool, double >	BDpair

typedef VertexTrack< N >::RefCountedLinearizedTrackState	RefCountedLinearizedTrackState

typedef VertexTrack< N >::RefCountedRefittedTrackState	RefCountedRefittedTrackState

typedef CachingVertex< N >::RefCountedVertexTrack	RefCountedVertexTrack

Public Types inherited from VertexTrackCompatibilityEstimator< N >
typedef std::pair< bool, double >	BDpair

typedef VertexTrack< N >::RefCountedLinearizedTrackState	RefCountedLinearizedTrackState

typedef CachingVertex< N >::RefCountedVertexTrack	RefCountedVertexTrack

Public Member Functions
virtual KalmanVertexTrackCompatibilityEstimator< N > *	clone () const

virtual BDpair	estimate (const CachingVertex< N > &vrt, const RefCountedVertexTrack track, unsigned int hint=UINT_MAX) const

virtual BDpair	estimate (const CachingVertex< N > &v, const RefCountedLinearizedTrackState track, unsigned int hint=UINT_MAX) const

virtual BDpair	estimate (const reco::Vertex &vertex, const reco::TransientTrack &track) const

	KalmanVertexTrackCompatibilityEstimator ()

virtual	~KalmanVertexTrackCompatibilityEstimator ()

Public Member Functions inherited from VertexTrackCompatibilityEstimator< N >
	VertexTrackCompatibilityEstimator ()

virtual	~VertexTrackCompatibilityEstimator ()

Private Member Functions
BDpair	estimateDifference (const CachingVertex< N > &more, const CachingVertex< N > &less, const RefCountedVertexTrack track) const

BDpair	estimateFittedTrack (const CachingVertex< N > &v, const RefCountedVertexTrack track) const

BDpair	estimateNFittedTrack (const CachingVertex< N > &v, const RefCountedVertexTrack track) const

Private Attributes
KVFHelper< N >	helper

KalmanVertexTrackUpdator< N >	trackUpdator

KalmanVertexUpdator< N >	updator

VertexTrackFactory< N >	vTrackFactory

Detailed Description

template<unsigned int N>
class KalmanVertexTrackCompatibilityEstimator< N >

Calculates the compatiblity of a track with respect to a vertex using the Kalman filter algorithms. The compatibility is computed from the squared standardized residuals between the track and the vertex. If track and vertex errors are Gaussian and correct, this quantity is distributed as chi**2(ndf=2)). Can be used to identify outlying tracks.

Definition at line 25 of file KalmanVertexTrackCompatibilityEstimator.h.

Member Typedef Documentation

template<unsigned int N>

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

Definition at line 33 of file KalmanVertexTrackCompatibilityEstimator.h.

template<unsigned int N>

typedef VertexTrack<N>::RefCountedLinearizedTrackState KalmanVertexTrackCompatibilityEstimator< N >::RefCountedLinearizedTrackState

Definition at line 31 of file KalmanVertexTrackCompatibilityEstimator.h.

template<unsigned int N>

typedef VertexTrack<N>::RefCountedRefittedTrackState KalmanVertexTrackCompatibilityEstimator< N >::RefCountedRefittedTrackState

Definition at line 32 of file KalmanVertexTrackCompatibilityEstimator.h.

template<unsigned int N>

typedef CachingVertex<N>::RefCountedVertexTrack KalmanVertexTrackCompatibilityEstimator< N >::RefCountedVertexTrack

Definition at line 30 of file KalmanVertexTrackCompatibilityEstimator.h.

Constructor & Destructor Documentation

template<unsigned int N>

KalmanVertexTrackCompatibilityEstimator< N >::KalmanVertexTrackCompatibilityEstimator ( )

inline

Definition at line 35 of file KalmanVertexTrackCompatibilityEstimator.h.

Referenced by KalmanVertexTrackCompatibilityEstimator< N >::clone().

35 {}

template<unsigned int N>

virtual KalmanVertexTrackCompatibilityEstimator< N >::~KalmanVertexTrackCompatibilityEstimator ( )

inlinevirtual

Definition at line 37 of file KalmanVertexTrackCompatibilityEstimator.h.

References KalmanVertexTrackCompatibilityEstimator< N >::estimate(), HiIsolationCommonParameters_cff::track, and findQualityFiles::v.

37 {}

Member Function Documentation

template<unsigned int N>

virtual KalmanVertexTrackCompatibilityEstimator<N>* KalmanVertexTrackCompatibilityEstimator< N >::clone ( ) const

inlinevirtual

Clone method

Implements VertexTrackCompatibilityEstimator< N >.

Definition at line 57 of file KalmanVertexTrackCompatibilityEstimator.h.

References KalmanVertexTrackCompatibilityEstimator< N >::estimateDifference(), KalmanVertexTrackCompatibilityEstimator< N >::estimateFittedTrack(), KalmanVertexTrackCompatibilityEstimator< N >::estimateNFittedTrack(), and KalmanVertexTrackCompatibilityEstimator< N >::KalmanVertexTrackCompatibilityEstimator().

   {
     return new KalmanVertexTrackCompatibilityEstimator(* this);
   }

template<unsigned int N>

KalmanVertexTrackCompatibilityEstimator< N >::BDpair KalmanVertexTrackCompatibilityEstimator< N >::estimate	(	const CachingVertex< N > &	vrt,
		const RefCountedVertexTrack	track,
		unsigned int	hint = `UINT_MAX`
	)		const

virtual

Track-to-vertex compatibility. The track weight is taken into account.

Parameters

track	The track for which the chi**2 has to be estimated.
v	The vertex against which the chi**2 has to be estimated.

Returns: The chi**2.

Implements VertexTrackCompatibilityEstimator< N >.

Definition at line 10 of file KalmanVertexTrackCompatibilityEstimator.cc.

References edmIntegrityCheck::d, l1t::tracks, and CachingVertex< N >::tracksRef().

Referenced by MultiVertexFitter::updateWeights(), and KalmanVertexTrackCompatibilityEstimator< N >::~KalmanVertexTrackCompatibilityEstimator().

 {
 //checking if the track passed really belongs to the vertex
 
   const std::vector<RefCountedVertexTrack> &tracks = vertex.tracksRef();
 
   if ( tracks.size()==0)
    return estimateNFittedTrack(vertex,tr);
 
   if (hint<tracks.size() ) {
     VertexTrackEqual<N> d(tr);
     if ( d(tracks[hint]))
       return estimateFittedTrack(vertex,tracks[hint]);
   }
 
   typename std::vector<RefCountedVertexTrack>::const_iterator pos 
     = find_if(tracks.begin(), tracks.end(), VertexTrackEqual<N>(tr));
   if(pos != tracks.end()) {
     return estimateFittedTrack(vertex,*pos);
   } else {
     return estimateNFittedTrack(vertex,tr);
   }
 } 

template<unsigned int N>

KalmanVertexTrackCompatibilityEstimator< N >::BDpair KalmanVertexTrackCompatibilityEstimator< N >::estimate	(	const CachingVertex< N > &	v,
		const RefCountedLinearizedTrackState	track,
		unsigned int	hint = `UINT_MAX`
	)		const

virtual

Methods giving back the compatibility estimation

Implements VertexTrackCompatibilityEstimator< N >.

Definition at line 39 of file KalmanVertexTrackCompatibilityEstimator.cc.

References CachingVertex< N >::vertexState().

 {
   RefCountedVertexTrack vertexTrack = vTrackFactory.vertexTrack(track,
                          vertex.vertexState());
   return estimate(vertex, vertexTrack,hint);
 }

template<unsigned int N>

KalmanVertexTrackCompatibilityEstimator< N >::BDpair KalmanVertexTrackCompatibilityEstimator< N >::estimate	(	const reco::Vertex &	vertex,
		const reco::TransientTrack &	track
	)		const

virtual

Implements VertexTrackCompatibilityEstimator< N >.

Definition at line 51 of file KalmanVertexTrackCompatibilityEstimator.cc.

References reco::Vertex::chi2(), reco::Vertex::covariance(), KalmanVertexTrackCompatibilityEstimator< N >::estimateFittedTrack(), spr::find(), LinearizedTrackStateFactory::linearizedTrackState(), reco::Vertex::position(), reco::TransientTrack::trackBaseRef(), reco::Vertex::tracks_begin(), and reco::Vertex::tracks_end().

 {   
 //   GlobalPoint linP(vertex.position().x(), vertex.position().z(),vertex.position().z());
     GlobalPoint linP(Basic3DVector<float> (vertex.position()));
 
   LinearizedTrackStateFactory lTrackFactory;
   RefCountedLinearizedTrackState linTrack = 
             lTrackFactory.linearizedTrackState(linP, track);
   GlobalError err(vertex.covariance());
   VertexState vState(linP, err);
   RefCountedVertexTrack vertexTrack = vTrackFactory.vertexTrack(linTrack, vState);
 
   std::vector<RefCountedVertexTrack> initialTracks(1, vertexTrack);
   CachingVertex<N> cachingVertex(linP, err, initialTracks,
                 vertex.chi2());
   // FIXME: this should work also for tracks without a persistent ref.
 //   return estimateNFittedTrack(cachingVertex, vertexTrack);
   if (find(vertex.tracks_begin(), vertex.tracks_end(), track.trackBaseRef()) != vertex.tracks_end())
   {
     return estimateFittedTrack(cachingVertex, vertexTrack);
   } else {
     return estimateNFittedTrack(cachingVertex, vertexTrack);
   }
 }

template<unsigned int N>

KalmanVertexTrackCompatibilityEstimator< N >::BDpair KalmanVertexTrackCompatibilityEstimator< N >::estimateDifference	(	const CachingVertex< N > &	more,
		const CachingVertex< N > &	less,
		const RefCountedVertexTrack	track
	)		const

private

Definition at line 116 of file KalmanVertexTrackCompatibilityEstimator.cc.

Referenced by KalmanVertexTrackCompatibilityEstimator< N >::clone(), and KalmanVertexTrackCompatibilityEstimator< N >::estimateNFittedTrack().

 {
   BDpair trackRes = helper.trackParameterChi2(track);
   return BDpair(trackRes.first, trackRes.second + helper.vertexChi2(less, more)) ;
 }

template<unsigned int N>

KalmanVertexTrackCompatibilityEstimator< N >::BDpair KalmanVertexTrackCompatibilityEstimator< N >::estimateFittedTrack	(	const CachingVertex< N > &	v,
		const RefCountedVertexTrack	track
	)		const

private

Definition at line 85 of file KalmanVertexTrackCompatibilityEstimator.cc.

References KalmanVertexTrackCompatibilityEstimator< N >::estimateNFittedTrack().

Referenced by KalmanVertexTrackCompatibilityEstimator< N >::clone(), and KalmanVertexTrackCompatibilityEstimator< N >::estimate().

 {
   //remove track from the vertex using the vertex updator
   // Using the update instead of the remove methode, we can specify a weight which
   // is different than then one which the vertex track has been defined with.
   //CachingVertex rVert = updator.remove(v, track);
   RefCountedVertexTrack newSmoothedTrack = trackUpdator.update(v, track);
 //   std::cout << newSmoothedTrack->smoothedChi2()<<" "<<estimateDifference(v,rVert,newSmoothedTrack)<<std::endl;
 //   return estimateDifference(v,rVert,newSmoothedTrack);
   return BDpair(true, newSmoothedTrack->smoothedChi2());
 }

template<unsigned int N>

KalmanVertexTrackCompatibilityEstimator< N >::BDpair KalmanVertexTrackCompatibilityEstimator< N >::estimateNFittedTrack	(	const CachingVertex< N > &	v,
		const RefCountedVertexTrack	track
	)		const

private

Definition at line 102 of file KalmanVertexTrackCompatibilityEstimator.cc.

References KalmanVertexTrackCompatibilityEstimator< N >::estimateDifference(), CachingVertex< N >::isValid(), CachingVertex< N >::totalChiSquared(), and gsfElectronCkfTrackCandidateMaker_cff::updator.

Referenced by KalmanVertexTrackCompatibilityEstimator< N >::clone(), and KalmanVertexTrackCompatibilityEstimator< N >::estimateFittedTrack().

 {
   // Using the update instead of the add methode, we can specify a weight which
   // is different than then one which the vertex track has been defined with.
   CachingVertex<N> rVert = updator.add(v, track);
   if (!rVert.isValid()) return BDpair(false,-1.);
   return BDpair(true, rVert.totalChiSquared()-v.totalChiSquared());
 }   