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

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

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

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

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

 AlgebraicVector deviations(int nStates) const override;

 int numberOfEquations() const override;

 PointingKinematicConstraint * clone() const override
 {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