fcnminuit.h

Go to the documentation of this file.

//
// Authors: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include "NangaParbat/chisquare.h"
 
#if MINUIT2_NOT_BUNDLED_WITH_ROOT
#include <Minuit2/FCNBase.h>
#include <Minuit2/FCNGradientBase.h>
#else
#include <root/Minuit2/FCNBase.h>
#include <root/Minuit2/FCNGradientBase.h>
#endif
 
 
namespace NangaParbat
{
  class FcnMinuit: public ROOT::Minuit2::FCNBase
  {
  public:
    FcnMinuit(ChiSquare const& chi2);
 
    double operator()(std::vector<double> const& pars) const;
 
    double Up() const { return 1. / (double) _chi2.GetDataPointNumber(); }
 
    void SetParameters(std::vector<double> const& pars) { _chi2.SetParameters(pars); };
 
  private:
    mutable ChiSquare _chi2; 
  };
 
  class FcnMinuitGrad: public ROOT::Minuit2::FCNGradientBase
  {
  public:
    FcnMinuitGrad(ChiSquare const& chi2);
 
    double operator()(std::vector<double> const& pars) const;
 
    std::vector<double> Gradient(const std::vector<double>& pars) const;
 
    double Up() const { return 1. / (double) _chi2.GetDataPointNumber(); }
 
    void SetParameters(std::vector<double> const& pars) { _chi2.SetParameters(pars); };
 
  private:
    mutable ChiSquare _chi2; 
  };
}