distribution.h

Go to the documentation of this file.

//
// APFEL++ 2017
//
// Author: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include "apfel/lagrangeinterpolator.h"
 
#include <map>
#include <functional>
 
namespace apfel
{
  class Distribution: public LagrangeInterpolator
  {
  public:
    Distribution() = delete;
 
 
    Distribution(Grid const& g);
 
    Distribution(Distribution const& obj);
 
    Distribution(Distribution                     const& obj,
                 std::vector<std::vector<double>> const& distsubgrid,
                 std::vector<double>              const& distjointgrid);
 
    Distribution(Grid                             const& g,
                 std::vector<std::vector<double>> const& distsubgrid,
                 std::vector<double>              const& distjointgrid);
 
    Distribution(Grid                                 const& g,
                 std::function<double(double const&)> const& InDistFunc);
 
    Distribution(Grid                                                const& g,
                 std::function<double(double const&, double const&)> const& InDistFunc,
                 double                                              const& Q);
 
    Distribution(Grid                                             const& g,
                 std::function<double(int const&, double const&)> const& InDistFunc,
                 int                                              const& ipdf);
 
    Distribution(Grid                                                            const& g,
                 std::function<double(int const&, double const&, double const&)> const& InDistFunc,
                 int                                                             const& ipdf,
                 double                                                          const& Q);
 
    void SetJointGrid(int const& ix, double const& x);
 
    void SetSubGrid(int const& ig, int const& ix, double const& x);
 
    void SetJointGrid(std::vector<double> const& jg);
 
    void SetSubGrid(int const& ig, std::vector<double> const& sg);
 
    void SetSubGrids(std::vector<std::vector<double>> const& sgs);
 
    Distribution Derivative() const;
 
    Distribution& operator  = (Distribution const& d);                         
    Distribution& operator *= (double const& s);                               
    Distribution& operator *= (std::function<double(double const&)> const& f); 
    Distribution& operator /= (double const& s);                               
    Distribution& operator *= (Distribution const& d);                         
    Distribution& operator += (Distribution const& d);                         
    Distribution& operator -= (Distribution const& d);                         
  };
 
  Distribution operator * (double const& s, Distribution rhs);                               
  Distribution operator * (Distribution lhs, double const& s);                               
  Distribution operator * (std::function<double(double const&)> const& f, Distribution rhs); 
  Distribution operator * (Distribution lhs, std::function<double(double const&)> const& f); 
  Distribution operator / (Distribution lhs, double const& s);                               
  Distribution operator + (Distribution lhs, Distribution const& rhs);                       
  Distribution operator - (Distribution lhs, Distribution const& rhs);                       
  Distribution operator * (Distribution lhs, Distribution const& rhs);                       
 
 
  std::map<int, Distribution> DistributionMap(Grid                                                               const& g,
                                              std::function<std::map<int, double>(double const&, double const&)> const& InDistFunc,
                                              double                                                             const& Q,
                                              std::vector<int>                                                   const& skip = {});
 
  std::map<int, Distribution> DistributionMap(Grid                                                const& g,
                                              std::function<std::map<int, double>(double const&)> const& InDistFunc,
                                              std::vector<int>                                    const& skip = {});
 
  std::map<int, Distribution> DistributionMap(Grid                                              const& g,
                                              std::function<std::vector<double>(double const&)> const& InDistFunc,
                                              int                                               const& NOutputs = 0);
 
  double Sum(Distribution const& InDist);
 
  double InnerProduct(Distribution const& d1, Distribution const& d2, double const& offset = 0);
}