SimplePointingConstraint.h

Go to the documentation of this file.

#ifndef SimplePointingConstraint_H
#define SimplePointingConstraint_H

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

class SimplePointingConstraint : public KinematicConstraint {
public:
  SimplePointingConstraint(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;

  SimplePointingConstraint* clone() const override { return new SimplePointingConstraint(*this); }

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

  GlobalPoint refPoint;
};

#endif