GsfVertexWeightCalculator Class Reference

#include <GsfVertexWeightCalculator.h>

Public Types
typedef ReferenceCountingPointer < LinearizedTrackState< 5 > >	RefCountedLinearizedTrackState

Public Member Functions
double	calculate (const VertexState &oldVertex, const RefCountedLinearizedTrackState track, double cov) const

Private Types
typedef LinearizedTrackState < 5 >::AlgebraicMatrixN3	AlgebraicMatrixN3
typedef LinearizedTrackState < 5 >::AlgebraicMatrixNM	AlgebraicMatrixNM
typedef LinearizedTrackState < 5 >::AlgebraicSymMatrixNN	AlgebraicSymMatrixNN
typedef LinearizedTrackState < 5 >::AlgebraicVectorN	AlgebraicVectorN

Detailed Description

Calulates the (non-normalised) weight of a component the new mixture of vertices for the Gaussian Sum vertex filter. (c.f. Th.Speer & R. Fruewirth, Comp.Phys.Comm 174, 935 (2006) )

Definition at line 14 of file GsfVertexWeightCalculator.h.

Member Typedef Documentation

typedef LinearizedTrackState<5>::AlgebraicMatrixN3 GsfVertexWeightCalculator::AlgebraicMatrixN3

private

Definition at line 28 of file GsfVertexWeightCalculator.h.

typedef LinearizedTrackState<5>::AlgebraicMatrixNM GsfVertexWeightCalculator::AlgebraicMatrixNM

private

Definition at line 29 of file GsfVertexWeightCalculator.h.

typedef LinearizedTrackState<5>::AlgebraicSymMatrixNN GsfVertexWeightCalculator::AlgebraicSymMatrixNN

private

Definition at line 30 of file GsfVertexWeightCalculator.h.

typedef LinearizedTrackState<5>::AlgebraicVectorN GsfVertexWeightCalculator::AlgebraicVectorN

private

Definition at line 27 of file GsfVertexWeightCalculator.h.

typedef ReferenceCountingPointer<LinearizedTrackState<5> > GsfVertexWeightCalculator::RefCountedLinearizedTrackState

Definition at line 16 of file GsfVertexWeightCalculator.h.

Member Function Documentation

double GsfVertexWeightCalculator::calculate	(	const VertexState &	oldVertex,
		const RefCountedLinearizedTrackState	track,
		double	cov
	)		const

Method to calculate the weight A negative weight is returned in case of error.

Definition at line 7 of file GsfVertexWeightCalculator.cc.

References a, b, change_name::diff, VertexState::error(), funct::exp(), edm::isNotFinite(), M_PI, GlobalErrorBase< T, ErrorWeightType >::matrix(), VertexState::position(), funct::pow(), mathSSE::sqrt(), histoStyle::weight, VertexState::weightInMixture(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by GsfVertexUpdator::createNewComponent(), and GsfVertexSmoother::createNewComponent().

                                                              {
  double previousWeight = oldVertex.weightInMixture() * track->weightInMixture();
  // Jacobians
  const AlgebraicMatrixN3& a = track->positionJacobian();
  const AlgebraicMatrixNM& b = track->momentumJacobian();

  //track information
  AlgebraicVectorN trackParameters = track->predictedStateParameters();
  AlgebraicSymMatrixNN trackParametersError = track->predictedStateError();
  //vertex information
  AlgebraicSymMatrix33 oldVertexError = oldVertex.error().matrix();
  //Vertex position
  GlobalPoint oldVertexPosition = oldVertex.position();
  AlgebraicVector3 oldVertexCoord;
  oldVertexCoord[0] = oldVertexPosition.x();
  oldVertexCoord[1] = oldVertexPosition.y();
  oldVertexCoord[2] = oldVertexPosition.z();

  // prior momentum information
  AlgebraicVector3 priorMomentum = track->predictedStateMomentumParameters();

  AlgebraicSymMatrix33 priorMomentumCov;
  priorMomentumCov(0, 0) = 1.0E-9;
  priorMomentumCov(1, 1) = 1.0E-6;
  priorMomentumCov(2, 2) = 1.0E-6;
  priorMomentumCov *= cov;

  AlgebraicVectorN diff = trackParameters - track->constantTerm() - a * oldVertexCoord - b * priorMomentum;
  track->checkParameters(diff);
  AlgebraicSymMatrixNN sigmaM =
      trackParametersError + ROOT::Math::Similarity(a, oldVertexError) + ROOT::Math::Similarity(b, priorMomentumCov);

  double sigmaDet;
  sigmaM.Det(sigmaDet);
  int ifail = !sigmaM.Invert();
  if (ifail != 0) {
    edm::LogWarning("GsfVertexWeightCalculator") << "S matrix inversion failed";
    return -1.;
  }

  double chi = ROOT::Math::Similarity(diff, sigmaM);
  ;  //SigmaM is now inverted !!!
  double weight = pow(2. * M_PI, -0.5 * 5) * sqrt(1. / sigmaDet) * exp(-0.5 * chi);

  if (edm::isNotFinite(weight) || sigmaDet <= 0.) {
    edm::LogWarning("GsfVertexWeightCalculator") << "Weight is NaN";
    return -1.;
  }

  return weight * previousWeight;
}