VertexDistanceXY Class Reference

#include <VertexDistanceXY.h>

Inheritance diagram for VertexDistanceXY:

Public Member Functions
VertexDistanceXY *	clone () const override
virtual float	compatibility (const GlobalPoint &vtx1Position, const GlobalError &vtx1PositionError, const GlobalPoint &vtx2Position, const GlobalError &vtx2PositionError) const=0
virtual float	compatibility (const reco::Vertex &, const reco::Vertex &) const
virtual float	compatibility (const reco::Vertex &, const VertexState &) const
virtual float	compatibility (const VertexState &, const reco::Vertex &) const
virtual float	compatibility (const VertexState &, const VertexState &) const
virtual Measurement1D	distance (const GlobalPoint &vtx1Position, const GlobalError &vtx1PositionError, const GlobalPoint &vtx2Position, const GlobalError &vtx2PositionError) const=0
Measurement1D	distance (const reco::Vertex &, const reco::Vertex &) const
Measurement1D	distance (const reco::Vertex &, const VertexState &) const
Measurement1D	distance (const VertexState &, const reco::Vertex &) const
Measurement1D	distance (const VertexState &, const VertexState &) const
Measurement1D	signedDistance (const reco::Vertex &primVtx, const reco::Vertex &secVtx, const GlobalVector &momentum) const override
	VertexDistanceXY ()
Public Member Functions inherited from VertexDistance
virtual float	compatibility (const reco::Vertex &, const reco::Vertex &) const
virtual float	compatibility (const reco::Vertex &, const VertexState &) const
virtual float	compatibility (const VertexState &, const reco::Vertex &) const
virtual float	compatibility (const VertexState &, const VertexState &) const
Measurement1D	distance (const reco::Vertex &, const reco::Vertex &) const
Measurement1D	distance (const reco::Vertex &, const VertexState &) const
Measurement1D	distance (const VertexState &, const reco::Vertex &) const
Measurement1D	distance (const VertexState &, const VertexState &) const
virtual	~VertexDistance ()

Private Member Functions
float	compatibility (const GlobalPoint &vtx1Position, const GlobalError &vtx1PositionError, const GlobalPoint &vtx2Position, const GlobalError &vtx2PositionError) const override
Measurement1D	distance (const GlobalPoint &vtx1Position, const GlobalError &vtx1PositionError, const GlobalPoint &vtx2Position, const GlobalError &vtx2PositionError) const override

Private Attributes
AlgebraicSymMatrix22	theNullMatrix

Additional Inherited Members

Detailed Description

Computes the distance between two vertices as the chi-squared formed with their positions in the transverse plane.

Definition at line 11 of file VertexDistanceXY.h.

Constructor & Destructor Documentation

VertexDistanceXY()

VertexDistanceXY::VertexDistanceXY ( )

inline

Definition at line 13 of file VertexDistanceXY.h.

{}

Referenced by clone().

Member Function Documentation

clone()

VertexDistanceXY* VertexDistanceXY::clone ( void  ) const

inlineoverridevirtual

Implements VertexDistance.

Definition at line 27 of file VertexDistanceXY.h.

{ return new VertexDistanceXY(*this); }

References VertexDistanceXY().

compatibility() [1/6]

float VertexDistanceXY::compatibility ( const GlobalPoint &  vtx1Position, const GlobalError &  vtx1PositionError, const GlobalPoint &  vtx2Position, const GlobalError &  vtx2PositionError  ) const

overrideprivatevirtual

Implements VertexDistance.

Definition at line 40 of file VertexDistanceXY.cc.

                                                                                  {
  // error matrix of residuals
  AlgebraicSymMatrix33 err1 = vtx1PositionError.matrix();
  AlgebraicSymMatrix33 err2 = vtx2PositionError.matrix();
  AlgebraicSymMatrix22 error;
  error[0][0] = err1[0][0] + err2[0][0];
  error[0][1] = err1[0][1] + err2[0][1];
  error[1][1] = err1[1][1] + err2[1][1];
  if (error == theNullMatrix)
    return FLT_MAX;
 
  // position residuals
  GlobalVector diff = vtx2Position - vtx1Position;
  AlgebraicVector2 vDiff;
  vDiff[0] = diff.x();
  vDiff[1] = diff.y();
 
  // Invert error matrix of residuals
  bool ifail = !error.Invert();
  if (ifail) {
    throw cms::Exception("VertexDistanceXY::matrix inversion problem");
  }
 
  return ROOT::Math::Similarity(error, vDiff);
}

References change_name::diff, relativeConstraints::error, Exception, and GlobalErrorBase< T, ErrorWeightType >::matrix().

compatibility() [2/6]

virtual float VertexDistance::compatibility

compatibility() [3/6]

float VertexDistance::compatibility

Definition at line 51 of file VertexDistance.cc.

                                                                                {
  return compatibility(GlobalPoint(Basic3DVector<float>(vtx1.position())),
                       GlobalError(vtx1.covariance()),
                       GlobalPoint(Basic3DVector<float>(vtx2.position())),
                       GlobalError(vtx2.covariance()));
}

compatibility() [4/6]

float VertexDistance::compatibility

Definition at line 37 of file VertexDistance.cc.

                                                                                     {
  return compatibility(GlobalPoint(Basic3DVector<float>(vtx1.position())),
                       GlobalError(vtx1.covariance()),
                       vtx2.position(),
                       vtx2.error());
}

compatibility() [5/6]

float VertexDistance::compatibility

Definition at line 44 of file VertexDistance.cc.

                                                                                     {
  return compatibility(vtx1.position(),
                       vtx1.error(),
                       GlobalPoint(Basic3DVector<float>(vtx2.position())),
                       GlobalError(vtx2.covariance()));
}

compatibility() [6/6]

float VertexDistance::compatibility

Definition at line 33 of file VertexDistance.cc.

                                                                                          {
  return compatibility(vtx1.position(), vtx1.error(), vtx2.position(), vtx2.error());
}

distance() [1/6]

Measurement1D VertexDistanceXY::distance ( const GlobalPoint &  vtx1Position, const GlobalError &  vtx1PositionError, const GlobalPoint &  vtx2Position, const GlobalError &  vtx2PositionError  ) const

overrideprivatevirtual

Implements VertexDistance.

Definition at line 19 of file VertexDistanceXY.cc.

                                                                                     {
  AlgebraicSymMatrix33 error = vtx1PositionError.matrix() + vtx2PositionError.matrix();
  GlobalVector diff = vtx1Position - vtx2Position;
  AlgebraicVector3 vDiff;
  vDiff[0] = diff.x();
  vDiff[1] = diff.y();
  vDiff[2] = 0.;
 
  double dist = sqrt(pow(diff.x(), 2) + pow(diff.y(), 2));
 
  double err2 = ROOT::Math::Similarity(error, vDiff);
  double err = 0;
  if (dist != 0)
    err = sqrt(err2) / dist;
 
  return Measurement1D(dist, err);
}

References change_name::diff, runTheMatrix::err, relativeConstraints::error, GlobalErrorBase< T, ErrorWeightType >::matrix(), funct::pow(), and mathSSE::sqrt().

Referenced by VertexAnalysis::dist2D(), Onia2MuMuPAT::produce(), VertexTableProducer::produce(), TemplatedInclusiveVertexFinder< InputContainer, VTX >::produce(), and btagbtvdeep::vertexDxy().

distance() [2/6]

virtual Measurement1D VertexDistance::distance

distance() [3/6]

Measurement1D VertexDistance::distance

Definition at line 26 of file VertexDistance.cc.

                                                                                   {
  return distance(GlobalPoint(Basic3DVector<float>(vtx1.position())),
                  GlobalError(vtx1.covariance()),
                  GlobalPoint(Basic3DVector<float>(vtx2.position())),
                  GlobalError(vtx2.covariance()));
}

distance() [4/6]

Measurement1D VertexDistance::distance

Definition at line 12 of file VertexDistance.cc.

                                                                                        {
  return distance(GlobalPoint(Basic3DVector<float>(vtx1.position())),
                  GlobalError(vtx1.covariance()),
                  vtx2.position(),
                  vtx2.error());
}

distance() [5/6]

Measurement1D VertexDistance::distance

Definition at line 19 of file VertexDistance.cc.

                                                                                        {
  return distance(vtx1.position(),
                  vtx1.error(),
                  GlobalPoint(Basic3DVector<float>(vtx2.position())),
                  GlobalError(vtx2.covariance()));
}

distance() [6/6]

Measurement1D VertexDistance::distance

Definition at line 8 of file VertexDistance.cc.

                                                                                             {
  return distance(vtx1.position(), vtx1.error(), vtx2.position(), vtx2.error());
}

signedDistance()

Measurement1D VertexDistanceXY::signedDistance ( const reco::Vertex &  primVtx, const reco::Vertex &  secVtx, const GlobalVector &  momentum  ) const

overridevirtual

The signed distance is computed using a vector from the primary to the secondary vertex and a given reference vector. The sign is determined by the scalar product of the x,y component of the vector connecting the vertices and the reference vector: if the scalar product is greater than zero, the sign is +1, else -1

Implements VertexDistance.

Definition at line 9 of file VertexDistanceXY.cc.

                                                                                   {
  Measurement1D unsignedDistance = distance(vtx1, vtx2);
  Basic3DVector<float> diff = Basic3DVector<float>(vtx2.position()) - Basic3DVector<float>(vtx1.position());
  if ((momentum.x() * diff.x() + momentum.y() * diff.y()) < 0)
    return Measurement1D(-1.0 * unsignedDistance.value(), unsignedDistance.error());
  return unsignedDistance;
}

References change_name::diff, HLT_2018_cff::distance, Measurement1D::error(), reco::Vertex::position(), Measurement1D::value(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Member Data Documentation

theNullMatrix

AlgebraicSymMatrix22 VertexDistanceXY::theNullMatrix

private

Definition at line 33 of file VertexDistanceXY.h.