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