analyticsolutions.h

Go to the documentation of this file.

//
// ePDF
//
 
#pragma once
 
#include "ePDF/alphaem.h"
#include "ePDF/constants.h"
 
namespace ePDF
{
  const std::vector<bool> orderRLL{true, true, true, false, false, false};
  const std::vector<bool> orderRNLL{true, true, true, true, true, true};
 
  class AnalyticSolutions
  {
  public:
 
    AnalyticSolutions(YAML::Node const& config);
 
    AnalyticSolutions(YAML::Node const& config, int const& orderA, std::vector<bool> const orderR);
 
    void SetPhotonMatching(std::vector<double> const& vec);
 
    void SetNumInt(bool const& numint);
 
    void SetRecHat(bool const& rechatON);
 
    double GetSolution(double const& x, double const& Q, int const& id, std::string const& term);
 
    double Evolve(double const& x, double const& Q, int const& id);
 
    std::vector<double> Evolve(double const& x, double const& Q);
 
    double TestPhoton(double const& x, double const& Q, std::string test);
 
  private:
    AlphaQED _aQED;               
    double const _Qi;             
    double const _nl;             
 
    double const _b0;
    double const _b1;
 
    double const _a0twopi;
    double const _L0;
    double _atwopi;
    double _t;
    double _eta0;
 
    std::vector<bool> _orderR;
    int               _orderA;
 
    double _x0photon;
    double _x1photon;
    double _pphoton;
 
    bool _numint;
 
    bool _rechatON;
 
 
    double RecSeries(std::vector<double> const& series) const;
 
    void Warmup(double const& Q);
 
    double AsySolution(double const& x, int const& id) const;
    double AsyBarSolution(double const& x, int const& id) const;
    double RecSolution(double const& x, int const& id) const;
    double RecBarSolution(double const& x, int const& id) const;
 
    std::vector<double> RecBarNS(double const& z) const;
    std::vector<double> RecBarS(double const& z) const;
    std::vector<double> RecBarG(double const& z) const;
 
    std::vector<double> RecHatNS(double const& z) const;
    std::vector<double> RecHatS(double const& z) const;
    std::vector<double> RecHatG(double const& z) const;
 
    double AsyEleAF(double const& z) const;
    double AsyEleAR(double const& z) const;
    std::vector<double> AsyEleBarAF(double const& z) const;
    std::vector<double> AsyEleBarAR(double const& z) const;
 
    double AsyPhoton(double const& z) const;
 
    double sumRiMi(double const& C1, double const& C2, double const& C3,
                   double const& C4, double const& C5,
                   double const& z, double const& k, double const& M1, double const& M2) const;
 
    std::vector<double> AsyPhotonBar(double const& z) const;
 
    double AsyPhotonSIMPLIFIED(double const& z) const;
  };
}