SmartPointingConstraint.h

Go to the documentation of this file.

#ifndef SmartPointingConstraint_H
#define SmartPointingConstraint_H

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

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

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

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

  GlobalPoint refPoint;
};

#endif