parameterisation.h

Go to the documentation of this file.

//
// Authors: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include <vector>
#include <iostream>
#include <functional>
#include <map>
#include <apfel/apfelxx.h>
 
namespace NangaParbat
{
  class Parameterisation
  {
  public:
    Parameterisation(std::string const& name = "", int const& nfuncs = 0, std::vector<double> pars = {}, bool const& anders = false);
 
    virtual ~Parameterisation() {};
 
    virtual void SetParameters(std::vector<double> const& pars) { _pars = pars; };
 
    virtual double Evaluate(double const& x, double const& b, double const& zeta, int const& ifunc) const { return 0; };
    virtual void EvaluateOnGrid() {};
    std::function<double(double const&, double const&, double const&, int const&)> Function() const;
 
    virtual double Derive(double const& x, double const& b, double const& zeta, int const& ifunc, int const& ipar) const { return 0; };
    virtual void DeriveOnGrid() {};
 
    std::function<double(double const&, double const&, double const&, int const&, int const&)> Derivative() const;
 
    virtual std::string LatexFormula() const { return ""; };
 
    virtual std::vector<std::string> GetParameterNames() const { return {}; };
 
    virtual std::string GetDescription() const { return ""; };
 
    std::string                 GetName()              const { return _name; }
    int                         GetNumberOfFunctions() const { return _nfuncs; }
    virtual std::vector<double> GetParameters()        const { return _pars; }
    virtual int                 GetParameterNumber()   const { return _pars.size(); }
    bool                        HasGradient()          const { return _anders; }
 
    // FF_SIDIS
    virtual std::function<std::map<int, double>(double const&, double const&)> LHAPDF_Function()           const { return nullptr; }
    virtual std::function<std::map<int, double>(double const&, double const&)> LHAPDF_Derivative(int)      const { return nullptr; }
    virtual std::function<apfel::Set<apfel::Distribution>(double const&)>      DistributionFunction()      const { return nullptr; }
    virtual std::function<apfel::Set<apfel::Distribution>(double const&)>      DistributionDerivative(int) const { return nullptr; }
 
  protected:
    std::string         _name;   
    int                 _nfuncs; 
    std::vector<double> _pars;   
    bool                _anders; 
  };
}