initialiseevolution.h

Go to the documentation of this file.

//
// APFEL++ 2017
//
// Author: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include "apfel/evolutionsetup.h"
#include "apfel/tabulateobject.h"
#include "apfel/dglapbuilder.h"
#include "apfel/config.h"
 
namespace apfel
{
  struct LHKnotArray
  {
    std::vector<double> xs;     
    std::vector<double> q2s;    
    std::vector<double> xfs;    
  };
 
  class InitialiseEvolution
  {
  public:
    InitialiseEvolution(EvolutionSetup const& setup, bool const& WriteGrid = false, std::string const& GridHeader = "");
 
    bool CheckSetup() const;
 
    void ReportSetup() const;
 
    void InitialiseCouplings();
 
    void InitialiseDglapObject();
 
    void TabulateEvolution(std::function<std::map<int, double>(double const&, double const&)> const& InSet);
 
    void WriteGridInfo();
 
    void WriteGrid();
 
    double Alphas(double const& mu) const { return _as(mu); }
 
    std::map<double, std::map<int, LHKnotArray>> KnotArray() const { return _KnotArray; }
 
    TabulateObject<Set<Distribution>> TabulatedDistributions() const { return *_TabulatedDists; }
 
    EvolutionSetup GetEvolutionSetup() const { return _setup; }
 
  private:
    EvolutionSetup                                           _setup;          
    bool                                                     _WriteGrid;      
    std::string                                              _GridHeader;     
    std::unique_ptr<const Grid>                              _g;              
    std::function<double(double const&)>                     _as;             
    std::map<int, DglapObjects>                              _DglapObj;       
    std::map<double, std::map<int, LHKnotArray>>             _KnotArray;      
    std::unique_ptr<const TabulateObject<Set<Distribution>>> _TabulatedDists; 
  };
}