CombinedKinematicConstraint.h

Go to the documentation of this file.

#ifndef CombinedKinematicConstraint_H
#define CombinedKinematicConstraint_H

#include "RecoVertex/KinematicFitPrimitives/interface/MultiTrackKinematicConstraint.h"
#include "RecoVertex/KinematicFitPrimitives/interface/KinematicState.h"
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "RecoVertex/VertexPrimitives/interface/VertexException.h"

class CombinedKinematicConstraint : public MultiTrackKinematicConstraint {
public:
  CombinedKinematicConstraint(const std::vector<MultiTrackKinematicConstraint*>& constraintVector)
      : constraints(constraintVector) {
    if (constraints.empty())
      throw VertexException("CombinedKinematicConstraint::<1 constraints passed.");
  }

  AlgebraicVector value(const std::vector<KinematicState>& states, const GlobalPoint& point) const override;

  AlgebraicMatrix parametersDerivative(const std::vector<KinematicState>& states,
                                       const GlobalPoint& point) const override;

  AlgebraicMatrix positionDerivative(const std::vector<KinematicState>& states,
                                     const GlobalPoint& point) const override;

  int numberOfEquations() const override;

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

private:
  std::vector<MultiTrackKinematicConstraint*> constraints;
};
#endif