MultipleKinematicConstraint.h

Go to the documentation of this file.

#ifndef MultipleKinematicConstraint_H
#define MultipleKinematicConstraint_H

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


class MultipleKinematicConstraint : public KinematicConstraint
{
public:

  MultipleKinematicConstraint()
  {em = true;}


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

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

  int numberOfEquations() const override;
  
  void addConstraint(KinematicConstraint * newConst) const; 

  AlgebraicVector deviations(int nStates) const override;
  
  bool isEmpty() const
  {return em;}
  
  MultipleKinematicConstraint * clone() const override
  {return new MultipleKinematicConstraint(*this);}
  
private:

  mutable std::vector<KinematicConstraint *>  cts;
  
  mutable bool em;
};
#endif