GsfVertexTrackCompatibilityEstimator Class Reference

#include <GsfVertexTrackCompatibilityEstimator.h>

Inheritance diagram for GsfVertexTrackCompatibilityEstimator:

Public Types
typedef CachingVertex< 5 >::RefCountedVertexTrack	RefCountedVertexTrack
Public Types inherited from VertexTrackCompatibilityEstimator< 5 >
typedef std::pair< bool, double >	BDpair
typedef VertexTrack< N >::RefCountedLinearizedTrackState	RefCountedLinearizedTrackState
typedef CachingVertex< N >::RefCountedVertexTrack	RefCountedVertexTrack

Public Member Functions
GsfVertexTrackCompatibilityEstimator *	clone () const override
BDpair	estimate (const CachingVertex< 5 > &vrt, const RefCountedVertexTrack track, unsigned int hint=UINT_MAX) const override
BDpair	estimate (const CachingVertex< 5 > &v, const RefCountedLinearizedTrackState track, unsigned int hint=UINT_MAX) const override
BDpair	estimate (const reco::Vertex &vertex, const reco::TransientTrack &track) const override
	GsfVertexTrackCompatibilityEstimator ()
	~GsfVertexTrackCompatibilityEstimator () override
Public Member Functions inherited from VertexTrackCompatibilityEstimator< 5 >
virtual BDpair	estimate (const CachingVertex< N > &v, const RefCountedLinearizedTrackState track, unsigned int hint=UINT_MAX) const=0
virtual BDpair	estimate (const CachingVertex< N > &v, const RefCountedVertexTrack track, unsigned int hint=UINT_MAX) const=0
	VertexTrackCompatibilityEstimator ()
virtual	~VertexTrackCompatibilityEstimator ()

Private Member Functions
BDpair	estimateFittedTrack (const CachingVertex< 5 > &v, const RefCountedVertexTrack track) const
BDpair	estimateNFittedTrack (const CachingVertex< 5 > &v, const RefCountedVertexTrack track) const

Private Attributes
MultiPerigeeLTSFactory	lTrackFactory
GsfVertexUpdator	updator
VertexTrackFactory< 5 >	vTrackFactory

Detailed Description

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 GsfVertexTrackCompatibilityEstimator.h.

Member Typedef Documentation

RefCountedVertexTrack

typedef CachingVertex<5>::RefCountedVertexTrack GsfVertexTrackCompatibilityEstimator::RefCountedVertexTrack

Definition at line 26 of file GsfVertexTrackCompatibilityEstimator.h.

Constructor & Destructor Documentation

GsfVertexTrackCompatibilityEstimator()

GsfVertexTrackCompatibilityEstimator::GsfVertexTrackCompatibilityEstimator ( )

inline

Definition at line 28 of file GsfVertexTrackCompatibilityEstimator.h.

Referenced by clone().

{}

~GsfVertexTrackCompatibilityEstimator()

GsfVertexTrackCompatibilityEstimator::~GsfVertexTrackCompatibilityEstimator ( )

inlineoverride

Definition at line 30 of file GsfVertexTrackCompatibilityEstimator.h.

{}

Member Function Documentation

clone()

GsfVertexTrackCompatibilityEstimator* GsfVertexTrackCompatibilityEstimator::clone ( void  ) const

inlineoverridevirtual

Clone method

Implements VertexTrackCompatibilityEstimator< 5 >.

Definition at line 50 of file GsfVertexTrackCompatibilityEstimator.h.

References GsfVertexTrackCompatibilityEstimator().

                                                               {
    return new GsfVertexTrackCompatibilityEstimator(*this);
  }

estimate() [1/3]

std::pair< bool, double > GsfVertexTrackCompatibilityEstimator::estimate ( const CachingVertex< 5 > &  vrt, const RefCountedVertexTrack  track, unsigned int  hint = UINT_MAX  ) const

override

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.

Definition at line 20 of file GsfVertexTrackCompatibilityEstimator.cc.

References DiMuonV_cfg::tracks, and bphysicsOniaDQM_cfi::vertex.

                                                                                                {
  //checking if the track passed really belongs to the vertex
  std::vector<RefCountedVertexTrack> tracks = vertex.tracks();
  std::vector<RefCountedVertexTrack>::iterator pos = find_if(tracks.begin(), tracks.end(), VertexTrackEqual<5>(tr));
  if (pos != tracks.end()) {
    return estimateFittedTrack(vertex, *pos);
  } else {
    return estimateNFittedTrack(vertex, tr);
  }
}

estimate() [2/3]

std::pair< bool, double > GsfVertexTrackCompatibilityEstimator::estimate ( const CachingVertex< 5 > &  v, const RefCountedLinearizedTrackState  track, unsigned int  hint = UINT_MAX  ) const

override

Definition at line 33 of file GsfVertexTrackCompatibilityEstimator.cc.

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

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

estimate() [3/3]

std::pair< bool, double > GsfVertexTrackCompatibilityEstimator::estimate ( const reco::Vertex &  vertex, const reco::TransientTrack &  track  ) const

overridevirtual

Implements VertexTrackCompatibilityEstimator< 5 >.

Definition at line 40 of file GsfVertexTrackCompatibilityEstimator.cc.

References submitPVResolutionJobs::err, spr::find(), HLT_2024v14_cff::track, and bphysicsOniaDQM_cfi::vertex.

                                                                                                              {
  //   GlobalPoint linP(vertex.position().x(), vertex.position().z(),vertex.position().z());
  GlobalPoint linP(Basic3DVector<float>(vertex.position()));

  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<5> cachingVertex(linP, err, initialTracks, vertex.chi2());
  // FIXME: this should work also for tracks without a persistent ref.
  //   const TrackTransientTrack* ttt = dynamic_cast<const TrackTransientTrack*>(track.basicTransientTrack());
  //   if ((ttt!=0) &&
  if (find(vertex.tracks_begin(), vertex.tracks_end(), track.trackBaseRef()) != vertex.tracks_end()) {
    return estimateFittedTrack(cachingVertex, vertexTrack);
  } else {
    return estimateNFittedTrack(cachingVertex, vertexTrack);
  }
}

estimateFittedTrack()

std::pair< bool, double > GsfVertexTrackCompatibilityEstimator::estimateFittedTrack ( const CachingVertex< 5 > &  v, const RefCountedVertexTrack  track  ) const

private

Definition at line 65 of file GsfVertexTrackCompatibilityEstimator.cc.

References HLT_2024v14_cff::track.

                                                                        {
  //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.
  if (track->refittedStateAvailable())
    return BDpair(true, track->smoothedChi2());
  throw VertexException("GsfVertexTrackCompatibilityEstimator::vertex has to be smoothed.");

  //   CachingVertex rVert = updator.remove(v, track);
  //   RefCountedVertexTrack newSmoothedTrack = trackUpdator.update(v, track);
  //   return estimateDifference(v,rVert,newSmoothedTrack);
}

estimateNFittedTrack()

std::pair< bool, double > GsfVertexTrackCompatibilityEstimator::estimateNFittedTrack ( const CachingVertex< 5 > &  v, const RefCountedVertexTrack  track  ) const

private

Definition at line 81 of file GsfVertexTrackCompatibilityEstimator.cc.

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

                                                                        {
  // 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<5> rVert = updator.add(v, track);
  if (!rVert.isValid())
    return BDpair(false, -1.);
  return BDpair(true, rVert.totalChiSquared() - v.totalChiSquared());
}

Member Data Documentation

lTrackFactory

MultiPerigeeLTSFactory GsfVertexTrackCompatibilityEstimator::lTrackFactory

private

Definition at line 60 of file GsfVertexTrackCompatibilityEstimator.h.

updator

GsfVertexUpdator GsfVertexTrackCompatibilityEstimator::updator

private

Definition at line 58 of file GsfVertexTrackCompatibilityEstimator.h.

vTrackFactory

VertexTrackFactory<5> GsfVertexTrackCompatibilityEstimator::vTrackFactory

private

Definition at line 61 of file GsfVertexTrackCompatibilityEstimator.h.