ModelFunction_TwoPointCorrelation.h

Go to the documentation of this file.

/********************************************************************
 *  Copyright (C) 2016 by Federico Marulli and Alfonso Veropalumbo  *
 *  federico.marulli3@unibo.it                                      *
 *                                                                  *
 *  This program is free software; you can redistribute it and/or   *
 *  modify it under the terms of the GNU General Public License as  *
 *  published by the Free Software Foundation; either version 2 of  *
 *  the License, or (at your option) any later version.             *
 *                                                                  *
 *  This program is distributed in the hope that it will be useful, *
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of  *
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the   *
 *  GNU General Public License for more details.                    *
 *                                                                  *
 *  You should have received a copy of the GNU General Public       *
 *  License along with this program; if not, write to the Free      *
 *  Software Foundation, Inc.,                                      *
 *  59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.       *
 ********************************************************************/
 
#ifndef __MODFUNCTWOP__
#define __MODFUNCTWOP__
 
#include "Cosmology.h"
#include "Modelling_MassObservableRelation.h"
 
 
// ============================================================================
 
 
namespace cbl {
 
  namespace modelling {
 
    namespace twopt {
 
      struct STR_data_model {
 
    int poly_order;
 
    std::shared_ptr<cosmology::Cosmology> cosmology;
 
    std::shared_ptr<cosmology::Cosmology> test_cosmology;
    
    std::shared_ptr<modelling::massobsrel::Modelling_MassObservableRelation> scaling_relation;
    
    double z_abs_err = -1;
    
    double proxy_rel_err = -1;
 
    double redshift;
 
    std::string method_Pk;
 
    std::string output_dir;
 
    bool store_output;
 
    std::string output_root;
 
    bool NL;
 
    double sigmaNL;
 
    double sigmaNL_perp;
 
    double sigmaNL_par;
 
    int norm;
 
    double k_min;
 
    double k_max;
 
    double aa;
 
    bool GSL;
 
    double prec;
 
    std::string file_par;
 
    std::shared_ptr<glob::FuncGrid> func_xi;
 
    std::shared_ptr<glob::FuncGrid> func_xi_;
 
    std::shared_ptr<glob::FuncGrid> func_xi__;
 
    double pi_max;
 
    double r_min;
 
    double r_max;
 
    int step;
 
    double linear_growth_rate_z;
 
    double sigma8_z;
 
    double var;
 
    std::vector<cosmology::CosmologicalParameter> Cpar;
 
    int FV;
 
    int bias_nl;
 
    double bA;
 
    int xiType;
 
    double k_star;
 
    int xiNL;
 
    double v_min;
 
    double v_max;
 
    int step_v;
 
    int xi_real_index;
 
    std::vector<bool> use_pole;
 
    std::vector<int> dataset_order;
 
    int nmultipoles;
 
    int nWedges;
 
    std::vector<std::vector<double>> mu_integral_limits;
 
    std::vector<double> kk;
 
    std::vector<double> rr;
 
    std::shared_ptr<glob::FuncGrid> func_Pk;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_NW;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_nonlin;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_DeltaDelta;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_DeltaTheta;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_ThetaTheta;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_A11;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_A12;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_A22;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_A23;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_A33;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_B12;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_B13;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_B14;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_B22;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_B23;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_B24;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_B33;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_B34;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_B44;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_b2d;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_b2v;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_b22;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_bs2d;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_bs2v;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_b2s2;
 
    std::shared_ptr<glob::FuncGrid> func_Pk_bs22;
 
    std::shared_ptr<glob::FuncGrid> func_sigma32Pklin;
 
    std::shared_ptr<glob::FuncGrid> func_Pk1loop;
 
    std::vector<std::shared_ptr<glob::FuncGrid>> funcs_pk;
 
    std::vector<std::shared_ptr<glob::FuncGrid>> func_multipoles;
 
    std::vector<std::shared_ptr<glob::FuncGrid>> func_wedges;
 
    std::string Pk_mu_model;
 
    double bias;
 
    double DVfid;
 
    double DAfid;
 
    double HHfid;
 
    double Delta;
    
    double Delta_input;
 
    bool isDelta_critical;
 
    std::shared_ptr<cbl::random::NormalRandomNumbers> gau_ran;
 
    std::shared_ptr<cbl::data::Data> cluster_mass_proxy;
 
    std::vector<double> cluster_mass_proxy_error;
 
    std::string model_MF;
 
    std::string model_bias;
 
    std::string meanType;
 
    std::shared_ptr<glob::FuncGrid> func_sigma;
 
    std::function<double(const double)> cosmopar_bias_interp_1D;
 
    std::function<double(const double, const double)> cosmopar_bias_interp_2D;
 
    std::shared_ptr<glob::FuncGrid> interp_SelectionFunction_cut;
 
    std::shared_ptr<glob::FuncGrid2D> interp_SelectionFunction;
 
    double z_min;
 
    double z_max;
 
    int mass_step;
 
    double Mass_min;
 
    double Mass_max;
 
    std::vector<double> mass;
 
    cosmology::Cosmology cosmology_mass;
 
    std::vector<double> redshift_source;
 
    STR_data_model () = default;
      };
 
      std::vector<double> true_k_mu_AP (const double kk, const double mu, const double alpha_perp, const double alpha_par);
 
      double Pk_l (const double kk, const int l, const std::string model, const std::vector<double> parameter, const std::vector<std::shared_ptr<glob::FuncGrid>> pk_interp, const double prec=1.e-5, const double alpha_perp = 1., const double alpha_par = 1.);
 
      std::vector<double> Pk_l (const std::vector<double> kk, const int l, const std::string model, const std::vector<double> parameter, const std::vector<std::shared_ptr<glob::FuncGrid>> pk_interp, const double prec=1.e-5, const double alpha_perp = 1., const double alpha_par = 1.);
 
       cbl::glob::FuncGrid Xil_interp (const std::vector<double> kk, const int l, const std::string model, const std::vector<double> parameter, const std::vector<std::shared_ptr<glob::FuncGrid>> pk_interp, const double prec=1.e-5, const double alpha_perp = 1., const double alpha_par = 1.);
 
 
      std::vector<std::vector<double>> Xi_l (const std::vector<double> rr, const int nmultipoles, const std::string model, const std::vector<double> parameter, const std::vector<std::shared_ptr<glob::FuncGrid>> pk_interp, const double prec=1.e-5, const double alpha_perp = 1., const double alpha_par = 1.);
 
      std::vector<double> Xi_l (const std::vector<double> rr, const std::vector<int> dataset_order, const std::vector<bool> use_pole, const std::string model, const std::vector<double> parameter, const std::vector<std::shared_ptr<glob::FuncGrid>> pk_interp, const double prec=1.e-5, const double alpha_perp = 1., const double alpha_par = 1.);
 
      std::vector<std::vector<double>> Xi_rppi (const std::vector<double> rp, const std::vector<double> pi, const std::string model, const std::vector<double> parameter, const std::vector<std::shared_ptr<glob::FuncGrid>> pk_interp, const double prec=1.e-5, const double alpha_perp = 1., const double alpha_par = 1.);
 
      double Xi_polar (const double rad_fid, const double mu_fid, const double alpha_perpendicular, const double alpha_parallel, const std::vector<std::shared_ptr<cbl::glob::FuncGrid>> xi_multipoles);
 
      std::vector<double> wp_from_Xi_rppi (const std::vector<double> rp, const double pimax, const std::string model, const std::vector<double> parameter, const std::vector<std::shared_ptr<glob::FuncGrid>> pk_interp, const double prec=1.e-5, const double alpha_perp = 1., const double alpha_par = 1.);
 
      std::vector<std::vector<double>> damped_Pk_terms (const std::vector<double> kk, const double linear_growth_rate, const double SigmaS, const std::shared_ptr<cbl::glob::FuncGrid> PkDM);
 
      std::vector<double> damped_Xi (const std::vector<double> ss, const double bias, const double linear_growth_rate, const double SigmaS, const std::vector<double> kk, const std::shared_ptr<cbl::glob::FuncGrid> PkDM);
      
      double Pkmu (const double kk, const double mu, const std::string model, const std::vector<double> parameter, const std::vector<std::shared_ptr<glob::FuncGrid>> pk_interp, const double alpha_perp = 1., const double alpha_par = 1.);
      
      double Pkmu_DeWiggled (const double kk, const double mu, const double sigmaNL_perp, const double sigmaNL_par, const double linear_growth_rate, const double bias, const double SigmaS, const std::shared_ptr<cbl::glob::FuncGrid> Pk, const std::shared_ptr<cbl::glob::FuncGrid> Pk_NW);
 
      double Pkmu_ModeCoupling (const double kk, const double mu, const double linear_growth_rate, const double bias, const double sigmav, const double AMC, const std::shared_ptr<cbl::glob::FuncGrid> PkLin, const std::shared_ptr<cbl::glob::FuncGrid> PkMC);
 
      double Pkmu_dispersion (const double kk, const double mu, const std::string DFoG, const double linear_growth_rate, const double bias, const double sigmav, const std::shared_ptr<cbl::glob::FuncGrid> Pklin);
 
      double Pkmu_Scoccimarro (const double kk, const double mu, const std::string DFoG, const double linear_growth_rate, const double bias, const double sigmav, const std::shared_ptr<cbl::glob::FuncGrid> Pk_DeltaDelta, const std::shared_ptr<cbl::glob::FuncGrid> Pk_DeltaTheta, const std::shared_ptr<cbl::glob::FuncGrid> Pk_ThetaTheta);
      
      double Pkmu_Scoccimarro_fitPezzotta (const double kk, const double mu, const std::string DFoG, const double linear_growth_rate, const double bias, const double sigmav, const double kd, const double kt, const std::shared_ptr<cbl::glob::FuncGrid> Pklin, const std::shared_ptr<cbl::glob::FuncGrid> Pknonlin);
 
      double Pkmu_Scoccimarro_fitBel (const double kk, const double mu, const std::string DFoG, const double linear_growth_rate, const double bias, const double sigmav, const double kd, const double bb, const double a1, const double a2, const double a3, const std::shared_ptr<cbl::glob::FuncGrid> Pklin, const std::shared_ptr<cbl::glob::FuncGrid> Pknonlin);
 
      double Pkmu_TNS (const double kk, const double mu, const std::string DFoG, const double linear_growth_rate, const double bias, const double sigmav, const std::shared_ptr<cbl::glob::FuncGrid> Pk_DeltaDelta, const std::shared_ptr<cbl::glob::FuncGrid> Pk_DeltaTheta, const std::shared_ptr<cbl::glob::FuncGrid> Pk_ThetaTheta, const std::shared_ptr<cbl::glob::FuncGrid> Pk_A11, const std::shared_ptr<cbl::glob::FuncGrid> Pk_A12, const std::shared_ptr<cbl::glob::FuncGrid> Pk_A22, const std::shared_ptr<cbl::glob::FuncGrid> Pk_A23, const std::shared_ptr<cbl::glob::FuncGrid> Pk_A33, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B12, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B13, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B14, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B22, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B23, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B24, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B33, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B34, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B44);
 
      double Pkmu_eTNS (const double kk, const double mu, const std::string DFoG, const double linear_growth_rate, const double bias, const double bias2, const double sigmav, const double Ncorr, const std::shared_ptr<cbl::glob::FuncGrid> Pk_DeltaDelta, const std::shared_ptr<cbl::glob::FuncGrid> Pk_DeltaTheta, const std::shared_ptr<cbl::glob::FuncGrid> Pk_ThetaTheta, const std::shared_ptr<cbl::glob::FuncGrid> Pk_A11, const std::shared_ptr<cbl::glob::FuncGrid> Pk_A12, const std::shared_ptr<cbl::glob::FuncGrid> Pk_A22, const std::shared_ptr<cbl::glob::FuncGrid> Pk_A23, const std::shared_ptr<cbl::glob::FuncGrid> Pk_A33, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B12, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B13, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B14, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B22, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B23, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B24, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B33, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B34, const std::shared_ptr<cbl::glob::FuncGrid> Pk_B44, const std::shared_ptr<cbl::glob::FuncGrid> Pk_b2d, const std::shared_ptr<cbl::glob::FuncGrid> Pk_b2v, const std::shared_ptr<cbl::glob::FuncGrid> Pk_b22, const std::shared_ptr<cbl::glob::FuncGrid> Pk_bs2d, const std::shared_ptr<cbl::glob::FuncGrid> Pk_bs2v, const std::shared_ptr<cbl::glob::FuncGrid> Pk_b2s2, const std::shared_ptr<cbl::glob::FuncGrid> Pk_bs22, const std::shared_ptr<cbl::glob::FuncGrid> sigma32Pklin);
 
    }
  }
}
 
#endif