MintFcn.h

Go to the documentation of this file.

#ifndef MINTFCN_HH
#define MINTFCN_HH

#include <iomanip>

#include "Minuit2/FCNBase.h"

#include "Mint/IMinimisable.h"

namespace MINT
{
  class MintFcn : public ROOT::Minuit2::FCNBase
  {
   public:
    MintFcn( IMinimisable* theFunction )
      :_theFunction(theFunction)
    {}

    virtual double operator()( const std::vector<double>& par ) const
    {
      for( unsigned int i=0; i<par.size(); ++i )
    _theFunction->getParSet()->getParPtr(i)->setCurrentFitVal(par[i]);
      _theFunction->parametersChanged();

      const double negll = _theFunction->getVal()-_best_min;

      std::cout << std::setprecision(std::numeric_limits<double>::digits10);
      std::cout << "-2logL = " << negll << std::endl;

      return negll;
    }

    void SetBestMin( const double& best_min )
    { _best_min += best_min; }

  private:
    virtual double Up() const
    { return _error_def; }

    IMinimisable* _theFunction;

    const double _error_def = 1.0;

    double _best_min = 0.0;
  };
} //namespace MINT

#endif //MINTFCN_HH