gtmdbuilder.h

Go to the documentation of this file.

//
// APFEL++ 2017
//
// Author: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include "apfel/grid.h"
#include "apfel/operator.h"
#include "apfel/set.h"
#include "apfel/dglapbuilder.h"
#include "apfel/tabulateobject.h"
#include "apfel/constants.h"
 
namespace apfel
{
  struct GtmdObjects
  {
    double                                    Threshold;
    double                                    xi;
    std::map<int, double>                     Beta;
    std::map<int, double>                     GammaFq;
    std::map<int, double>                     GammaFg;
    std::map<int, double>                     GammaK;
    std::map<int, std::vector<double>>        KCS;
    std::map<int, std::vector<Set<Operator>>> MatchingFunctions;
  };
  struct GtmdObjects {…};
 
 
  std::map<int, GtmdObjects> InitializeGtmdObjects(Grid                const& g,
                                                   std::vector<double> const& Thresholds,
                                                   double              const& xi,
                                                   double              const& IntEps = 1e-5);
 
 
  std::function<Set<Distribution>(double const&, double const&, double const&)> BuildGtmds(std::map<int, GtmdObjects>                      const& GtmdObj,
                                                                                           std::function<Set<Distribution>(double const&)> const& CollGPDs,
                                                                                           std::function<double(double const&)>            const& Alphas,
                                                                                           int                                             const& PerturbativeOrder,
                                                                                           double                                          const& Ci = 1,
                                                                                           double                                          const& IntEps = 1e-7);
 
  std::function<Set<Distribution>(double const&)> MatchGtmds(std::map<int, GtmdObjects>                      const& GtmdObj,
                                                             std::function<Set<Distribution>(double const&)> const& CollGPDs,
                                                             std::function<double(double const&)>            const& Alphas,
                                                             int                                             const& PerturbativeOrder,
                                                             double                                          const& Ci = 1);
 
  std::function<Set<Operator>(double const&)> MatchingFunctions(std::map<int, GtmdObjects>           const& GtmdObj,
                                                                std::function<double(double const&)> const& Alphas,
                                                                int                                  const& PerturbativeOrder,
                                                                double                               const& Ci = 1);
 
  std::function<std::vector<double>(double const&, double const&, double const&, double const&)> EvolutionFactors(std::map<int, GtmdObjects>           const& GtmdObj,
                                                                                                                  std::function<double(double const&)> const& Alphas,
                                                                                                                  int                                  const& PerturbativeOrder,
                                                                                                                  double                               const& Ci = 1,
                                                                                                                  double                               const& IntEps = 1e-7);
 
  std::function<double(double const&, double const&, double const&, double const&)> QuarkEvolutionFactor(std::map<int, GtmdObjects>           const& GtmdObj,
                                                                                                         std::function<double(double const&)> const& Alphas,
                                                                                                         int                                  const& PerturbativeOrder,
                                                                                                         double                               const& Ci = 1,
                                                                                                         double                               const& IntEps = 1e-7);
 
  std::function<double(double const&, double const&, double const&, double const&)> GluonEvolutionFactor(std::map<int, GtmdObjects>           const& GtmdObj,
                                                                                                         std::function<double(double const&)> const& Alphas,
                                                                                                         int                                  const& PerturbativeOrder,
                                                                                                         double                               const& Ci = 1,
                                                                                                         double                               const& IntEps = 1e-7);
}