evolutionsetup.h

Go to the documentation of this file.

//
// APFEL++ 2017
//
// Author: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include "apfel/lhtoypdfs.h"
 
#include <string>
#include <cmath>
#include <vector>
#include <functional>
#include <map>
 
namespace apfel
{
  struct EvolutionSetup
  {
    enum FlavourScheme: int {VFNS, FFNS};
 
    enum EvolutionTheory: int {QCD, QCD_QED};
 
    enum Virtuality: int {SPACE, TIME};
 
    enum EvolPolarisation: int {UNP, POL, TRANS};
 
    enum CouplingEvolution: int {exact, expanded};
 
    enum PDFEvolution: int {exactmu, exactalpha, expandalpha, truncated};
 
    enum MassRenScheme: int {POLE, MSBAR};
 
    struct GridParameters
    {
      int    nx;                  
      double xmin;                
      int    id;                  
    };
 
    std::string                 name;                
    double                      Q0;                  
    std::vector<GridParameters> GridParameters;      
    int                         nQg;                 
    double                      Qmin;                
    double                      Qmax;                
    int                         InterDegreeQ;        
    double                      Lambda;              
    int                         PerturbativeOrder;   
    FlavourScheme               FlavourScheme;       
    int                         Nf_FF;               
    EvolutionTheory             Theory;              
    Virtuality                  Virtuality;          
    EvolPolarisation            EvolPolarisation;    
    double                      AlphaQCDRef;         
    double                      QQCDRef;             
    double                      AlphaQEDRef;         
    double                      QQEDRef;             
    CouplingEvolution           CouplingEvolution;   
    PDFEvolution                PDFEvolution;        
    double                      xi;                  
    std::vector<double>         Thresholds;          
    std::vector<double>         Masses;              
    MassRenScheme               MassRenScheme;       
    double                      TauMass;             
    double                      GaussAccuracy;       
    std::vector<std::function<std::map<int, double>(double const&, double const&)>> InSet; 
 
    // *INDENT-OFF*
    EvolutionSetup():
      name("default"),
      Q0(sqrt(2)),
      GridParameters{{100, 1e-5, 3}, {100, 1e-1, 3}, {100, 6e-1, 3}, {80, 8.5e-1, 5}},
      nQg(50), Qmin(1), Qmax(1000), InterDegreeQ(3), Lambda(0.25),
      PerturbativeOrder(2),
      FlavourScheme(VFNS),
      Nf_FF(3),
      Theory(QCD),
      Virtuality(SPACE),
      EvolPolarisation(UNP),
      AlphaQCDRef(0.35), QQCDRef(sqrt(2)),
      AlphaQEDRef(7.496252e-3), QQEDRef(1.777),
      CouplingEvolution(exact), PDFEvolution(exactmu),
      xi(1),
      Thresholds{0, 0, 0, sqrt(2), 4.5, 175},
      Masses{Thresholds},
      MassRenScheme(POLE),
      TauMass(1.777),
      GaussAccuracy(1e-5),
      InSet({LHToyPDFs})
    {
    }
    // *INDENT-ON*
  };
}