FourMomentumKinematicConstraint.h

Go to the documentation of this file.

#ifndef FourMomentumKinematicConstraint_H
#define FourMomentumKinematicConstraint_H

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

class FourMomentumKinematicConstraint : public KinematicConstraint
{

public:

 FourMomentumKinematicConstraint(const AlgebraicVector& momentum,
                                 const AlgebraicVector& deviation);

virtual std::pair<AlgebraicVector,AlgebraicVector> value(const AlgebraicVector& exPoint) const;
 
virtual std::pair<AlgebraicMatrix, AlgebraicVector> derivative(const AlgebraicVector& exPoint) const;


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

virtual std::pair<AlgebraicVector, AlgebraicVector> value(const std::vector<RefCountedKinematicParticle> par) const;

 
virtual int numberOfEquations() const;

virtual AlgebraicVector deviations(int nStates) const;

virtual FourMomentumKinematicConstraint * clone() const
 {return new FourMomentumKinematicConstraint(*this);}

private:

AlgebraicVector mm;
AlgebraicVector dd;

};
#endif