KalmanVertexTrackCompatibilityEstimator< N > Class Template Reference

#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
KalmanVertexTrackCompatibilityEstimator< N > *	clone () const override
BDpair	estimate (const CachingVertex< N > &v, const RefCountedLinearizedTrackState track, unsigned int hint=UINT_MAX) const override
BDpair	estimate (const CachingVertex< N > &vrt, const RefCountedVertexTrack track, unsigned int hint=UINT_MAX) const override
BDpair	estimate (const reco::Vertex &vertex, const reco::TransientTrack &track) const override
	KalmanVertexTrackCompatibilityEstimator ()
	~KalmanVertexTrackCompatibilityEstimator () override
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 24 of file KalmanVertexTrackCompatibilityEstimator.h.

Member Typedef Documentation

BDpair

template<unsigned int N>

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

Definition at line 29 of file KalmanVertexTrackCompatibilityEstimator.h.

RefCountedLinearizedTrackState

template<unsigned int N>

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

Definition at line 27 of file KalmanVertexTrackCompatibilityEstimator.h.

RefCountedRefittedTrackState

template<unsigned int N>

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

Definition at line 28 of file KalmanVertexTrackCompatibilityEstimator.h.

RefCountedVertexTrack

template<unsigned int N>

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

Definition at line 26 of file KalmanVertexTrackCompatibilityEstimator.h.

Constructor & Destructor Documentation

KalmanVertexTrackCompatibilityEstimator()

template<unsigned int N>

KalmanVertexTrackCompatibilityEstimator< N >::KalmanVertexTrackCompatibilityEstimator ( )

inline

Definition at line 31 of file KalmanVertexTrackCompatibilityEstimator.h.

{}

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

~KalmanVertexTrackCompatibilityEstimator()

template<unsigned int N>

KalmanVertexTrackCompatibilityEstimator< N >::~KalmanVertexTrackCompatibilityEstimator ( )

inlineoverride

Definition at line 33 of file KalmanVertexTrackCompatibilityEstimator.h.

{}

Member Function Documentation

clone()

template<unsigned int N>

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

inlineoverridevirtual

Clone method

Implements VertexTrackCompatibilityEstimator< N >.

Definition at line 53 of file KalmanVertexTrackCompatibilityEstimator.h.

                                                                     {
    return new KalmanVertexTrackCompatibilityEstimator(*this);
  }

References KalmanVertexTrackCompatibilityEstimator< N >::KalmanVertexTrackCompatibilityEstimator().

estimate() [1/3]

template<unsigned int N>

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

overridevirtual

Methods giving back the compatibility estimation

Implements VertexTrackCompatibilityEstimator< N >.

Definition at line 33 of file KalmanVertexTrackCompatibilityEstimator.cc.

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

References HLT_FULL_cff::track, and bphysicsOniaDQM_cfi::vertex.

estimate() [2/3]

template<unsigned int N>

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

overridevirtual

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 8 of file KalmanVertexTrackCompatibilityEstimator.cc.

                                                                                             {
  //checking if the track passed really belongs to the vertex
 
  const std::vector<RefCountedVertexTrack> &tracks = vertex.tracksRef();
 
  if (tracks.empty())
    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);
  }
}

References ztail::d, tracks, and bphysicsOniaDQM_cfi::vertex.

Referenced by MultiVertexFitter::updateWeights().

estimate() [3/3]

template<unsigned int N>

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

overridevirtual

Implements VertexTrackCompatibilityEstimator< N >.

Definition at line 40 of file KalmanVertexTrackCompatibilityEstimator.cc.

                                                                     {
  //   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);
  }
}

References submitPVResolutionJobs::err, spr::find(), LinearizedTrackStateFactory::linearizedTrackState(), HLT_FULL_cff::track, and bphysicsOniaDQM_cfi::vertex.

estimateDifference()

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 95 of file KalmanVertexTrackCompatibilityEstimator.cc.

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

References gpuVertexFinder::more, and HLT_FULL_cff::track.

estimateFittedTrack()

template<unsigned int N>

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

private

Definition at line 67 of file KalmanVertexTrackCompatibilityEstimator.cc.

                                                                                                         {
  //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());
}

References HLT_FULL_cff::track, and findQualityFiles::v.

estimateNFittedTrack()

template<unsigned int N>

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

private

Definition at line 83 of file KalmanVertexTrackCompatibilityEstimator.cc.

                                                                                                          {
  // 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());
}

References CachingVertex< N >::isValid(), CachingVertex< N >::totalChiSquared(), HLT_FULL_cff::track, HLT_FULL_cff::updator, and findQualityFiles::v.

Member Data Documentation

helper

template<unsigned int N>

KVFHelper<N> KalmanVertexTrackCompatibilityEstimator< N >::helper

private

Definition at line 66 of file KalmanVertexTrackCompatibilityEstimator.h.

trackUpdator

template<unsigned int N>

KalmanVertexTrackUpdator<N> KalmanVertexTrackCompatibilityEstimator< N >::trackUpdator

private

Definition at line 64 of file KalmanVertexTrackCompatibilityEstimator.h.

updator

template<unsigned int N>

KalmanVertexUpdator<N> KalmanVertexTrackCompatibilityEstimator< N >::updator

private

Definition at line 63 of file KalmanVertexTrackCompatibilityEstimator.h.

vTrackFactory

template<unsigned int N>

VertexTrackFactory<N> KalmanVertexTrackCompatibilityEstimator< N >::vTrackFactory

private

Definition at line 65 of file KalmanVertexTrackCompatibilityEstimator.h.