PointingKinematicConstraint.h

Go to the documentation of this file.

#ifndef PointingKinematicConstraint_H
#define PointingKinematicConstraint_H

#include "RecoVertex/KinematicFitPrimitives/interface/KinematicConstraint.h"


class PointingKinematicConstraint : public KinematicConstraint
{
 public:
 
  PointingKinematicConstraint(const GlobalPoint& ref):refPoint(ref)
  {}

 virtual std::pair<AlgebraicVector, AlgebraicVector> value(const AlgebraicVector& exPoint) const;

 virtual std::pair<AlgebraicMatrix, AlgebraicVector> derivative(const AlgebraicVector& exPoint) const;

 virtual std::pair<AlgebraicMatrix, AlgebraicVector> derivative(const std::vector<RefCountedKinematicParticle> &par) const;

 virtual std::pair<AlgebraicVector, AlgebraicVector> value(const std::vector<RefCountedKinematicParticle> &par) const;

 virtual AlgebraicVector deviations(int nStates) const;

 virtual int numberOfEquations() const;

 virtual PointingKinematicConstraint * clone() const
 {return new PointingKinematicConstraint(*this);}
 
 private:

 std::pair<AlgebraicVector,AlgebraicVector> makeValue(const AlgebraicVector& exPoint)const ; 
 std::pair<AlgebraicMatrix, AlgebraicVector> makeDerivative(const AlgebraicVector& exPoint) const;
 
 GlobalPoint  refPoint;

};

#endif