matchedevolution.h

Go to the documentation of this file.

//
// APFEL++ 2017
//
// Author: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include <cmath>
#include <vector>
 
namespace apfel
{
  template<class T>
  class MatchedEvolution
  {
  public:
    MatchedEvolution() = delete;
 
    virtual ~MatchedEvolution() = default;
 
    MatchedEvolution(T                   const& ObjRef,
                     double              const& MuRef,
                     std::vector<double> const& Thresholds,
                     int                 const& nsteps = 10);
 
    virtual T EvolveObject(int const& nf, double const& t0, double const& t1, T const& Obj0) const;
 
    virtual T MatchObject(bool const& Up, int const& nf, T const& Obj) const = 0;
 
    virtual T Derivative(int const& nf, double const& Mu, T const& Obj) const = 0;
 
    T Evaluate(double const& mu) const;
 
 
    T GetObjectRef() const { return _ObjRef; }
 
    double GetMuRef() const { return _MuRef; }
 
    std::vector<double> GetThresholds() const { return _Thresholds; }
 
    int GetNumberOfSteps() const { return _nsteps; }
 
 
    void SetObjectRef(T ObjRef) { _ObjRef = std::move(ObjRef); }
 
    void SetMuRef(double const& MuRef) { _MuRef2 = MuRef * MuRef; _LogMuRef2 = log(_MuRef2); }
 
    void SetNumberOfSteps(int const& nsteps) { _nsteps = nsteps; }
  protected:
    T                   _ObjRef;         //<! Reference value of the object
    double              _MuRef;          //<! Reference scale of the object
    double              _MuRef2;         //<! Squared reference scale of the object
    double              _LogMuRef2;      //<! Log of the squared reference scale of the object
    std::vector<double> _Thresholds;     //<! Values of the thresholds
    int                 _nsteps;         //<! Number of steps of the RK algorithm
    std::vector<double> _Thresholds2;    //<! Squared quark threholds
    std::vector<double> _LogThresholds2; //<! Log of the squared quark threholds
  };
}