convolutiontable.h

Go to the documentation of this file.

//
// Author: Valerio Bertone: valerio.bertone@cern.ch
//
 
#pragma once
 
#include <string>
#include <functional>
#include <map>
#include <vector>
#include <valarray>
#include <yaml-cpp/yaml.h>
 
#include <NangaParbat/parameterisation.h>
#include <NangaParbat/cut.h>
 
namespace NangaParbat
{
  class ConvolutionTable
  {
  public:
    ConvolutionTable();
 
    //ConvolutionTable(YAML::Node const& table, double const& qToQmax = 100, std::vector<std::shared_ptr<Cut>> const& cuts = {}, double const& acc = 1e-7);
    ConvolutionTable(YAML::Node const& table, std::vector<double> const& cutParam = {}, std::vector<std::shared_ptr<Cut>> const& cuts = {}, double const& acc = 1e-7);
 
    //ConvolutionTable(std::string const& infile, double const& qToQmax = 100, std::vector<std::shared_ptr<Cut>> const& cuts = {}, double const& acc = 1e-7);
    ConvolutionTable(std::string const& infile, std::vector<double> const& cutParam = {}, std::vector<std::shared_ptr<Cut>> const& cuts = {}, double const& acc = 1e-7);
 
    std::map<double, double> ConvoluteDY(std::function<double(double const&, double const&, double const&)> const& fNP) const;
 
    std::map<double, double> ConvoluteSIDIS(std::function<double(double const&, double const&, double const&)> const& fNP,
                                            std::function<double(double const&, double const&, double const&)> const& DNP) const;
 
    virtual std::vector<double> GetPredictions(std::function<double(double const&, double const&, double const&)> const& fNP1,
                                               std::function<double(double const&, double const&, double const&)> const& fNP2) const;
 
    virtual std::vector<double> GetPredictions(std::function<double(double const&, double const&, double const&, int const&)> const& fNP) const;
 
    virtual  std::vector<double> GetPredictions(std::function<double(double const&, double const&, double const&, int const&)> const& fNP,
                                                std::function<double(double const&, double const&, double const&, int const&)> const& dNP) const;
 
    std::string                      GetName()       const { return _name; }
    int                              GetProcess()    const { return _proc; }
    double                           GetCME()        const { return _Vs; }
    std::vector<std::vector<double>> GetqTBins()     const { return _qTmap; }
    //double                           GetCutqToverQ() const { return _qToQmax; }
    std::vector<double>              GetcutParam()   const { return _cutParam; }
 
    std::valarray<bool> GetCutMask() const { return _cutmask; };
 
  protected:
    std::string                                                           const _name;     
    int                                                                   const _proc;     
    double                                                                const _Vs;       
    bool                                                                  const _IntqT;    
    std::vector<double>                                                   const _qTv;      
    std::vector<std::vector<double>>                                      const _qTmap;    
    std::vector<double>                                                   const _qTfact;   
    double                                                                const _prefact;  
    double                                                                const _prefact2; 
    std::vector<double>                                                   const _zOgata;   
    std::vector<double>                                                   const _Qg;       
    std::vector<double>                                                         _xig;      
    std::vector<double>                                                         _xbg;      
    std::vector<double>                                                         _zg;       
    std::map<double,std::vector<std::vector<double>>>                           _PSRed;    
    std::map<double,std::vector<std::vector<double>>>                           _dPSRed;   
    std::map<double,std::vector<std::vector<std::vector<double>>>>              _WDY;      
    std::map<double,std::vector<std::vector<std::vector<std::vector<double>>>>> _WSIDIS;   
    //double                                                                      _qToQmax;  //!< Maximum value allowed for the ratio qT / Q for DY
    std::vector<double>                                                   const _cutParam; 
    double                                                                      _acc;      
    std::vector<std::shared_ptr<Cut>>                                           _cuts;     
    std::valarray<bool>                                                         _cutmask;  
 
  public:
    virtual void SetInputFFs(std::function<std::map<int, double>(double const &, double const &)> const &InDistFunc) {};
    virtual void SetInputFFs(std::function<apfel::Set<apfel::Distribution>(double const&)> const& InDistFunc) {};
    virtual std::vector<double> GetPredictions(std::function<double(double const&, double const&, double const&)> const&) const { return {}; };
  };
}