xdistributions.h

Go to the documentation of this file.

//
// ePDF
//
 
#pragma once
 
#include "ePDF/ndistributions.h"
 
#include <yaml-cpp/yaml.h>
#include <complex>
#include <vector>
#include <memory>
#include <gsl/gsl_integration.h>
#include <functional>
 
namespace ePDF
{
  struct evol_params
  {
    double                          x;
    double                          Q;
    int                             id;
    std::shared_ptr<NDistributions> Ndist;
    evol_params operator = (evol_params const& p)
    {
      x     = p.x;
      Q     = p.Q;
      id    = p.id;
      Ndist = p.Ndist;
      return *this;
    }
  };
 
  class xDistributions
  {
  public:
    xDistributions(YAML::Node const& config);
 
    ~xDistributions();
 
    std::vector<double> Evolve(double const& x, double const& Q);
 
    std::vector<double> integrand(double const& y) const;
    std::vector<double> talbot(double const& y) const;
    std::vector<double> straight(double const& y) const;
 
    std::vector<double> trapezoid(double const& a, double const& b) const;
    std::vector<double> gauss(double const& a, double const& b) const;
    std::vector<double> gaussGSL(double const& a, double const& b);
 
    std::vector<std::complex<double>> Moments(std::complex<double> const& N, double const& Q) const;
 
    void SetParameters(evol_params const& p) { _p = p; };
 
  private:
    std::string               const _contour;
    std::string               const _integrator;
    double                    const _eps;
    evol_params                     _p;
    gsl_integration_workspace      *_w;
  };
}