92 inline std::vector<std::string>
CosmologicalModelNames () {
return {
"WMAP5",
"WMAP7",
"WMAP9",
"Planck13",
"Planck15",
"Planck15_TT",
"Planck18"}; }
119 inline std::vector<CosmologicalModel>
CosmologicalModelCast (
const std::vector<int> cosmologicalModelIndeces) {
return castFromValues<CosmologicalModel>(cosmologicalModelIndeces); }
209 inline std::vector<std::string>
CosmologicalParameterNames () {
return {
"Omega_matter_LCDM",
"Omega_matter",
"Omega_baryon",
"Omega_baryon_h2",
"Omega_matter",
"massless_neutrinos",
"massive_neutrinos",
"neutrino_mass",
"Omega_DE",
"Omega_radiation",
"H0",
"hh",
"scalar_amp",
"ln_scalar_amp",
"scalar_pivot",
"n_spec",
"w0",
"wa",
"fNL",
"sigma8",
"tau",
"rs"}; }
239 inline std::vector<CosmologicalParameter>
CosmologicalParameterCast (
const std::vector<int> cosmologicalParameterIndeces) {
return castFromValues<CosmologicalParameter>(cosmologicalParameterIndeces); }
440 double m_func_sigma (
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true, std::function<
double(
double)> filter={},
const bool unit1=
false)
const;
500 double m_sigma2R_notNormalised (
const double radius,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const bool unit1=
false)
const;
559 double m_sigma2M_notNormalised (
const double mass,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const bool unit1=
false)
const;
634 double m_MF_generator (
const double Mass,
const double Sigma,
const double Dln_Sigma,
const double redshift,
const std::string model_MF,
const double Delta=200.,
const bool default_delta=
true,
const double delta_t=1.686);
692 double m_MF_generator (
const double Mass,
const double Sigma,
const double Dln_Sigma,
const double redshift,
const double D_N,
const std::string model_MF,
const double Delta=200.,
const bool default_delta=
true,
const double delta_t=1.686);
716 double m_bias_halo_generator (
const double Sigma,
const double redshift,
const std::string author,
const double Delta=200.)
const;
742 double m_bias_halo_generator (
const double Sigma,
const double redshift,
const double D_N,
const std::string author,
const double Delta=200.)
const;
765 double m_elf_dz (
const double phi)
const;
817 double m_asn_dz (
const double ss)
const;
842 double m_serf_dz (
const double yy)
const;
873 void m_Table_Pk_CAMB_MPTbreeze (
const std::string code,
const bool NL, std::vector<double> &lgkk, std::vector<double> &lgPk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_max=100.)
const;
906 void m_Table_Pk_CAMB_MPTbreeze (
const std::string code,
const bool NL, std::vector<std::vector<double>> &lgkk, std::vector<std::vector<double>> &lgPk,
const std::vector<double> redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_max=100.)
const;
933 void m_Table_Pk_CLASS (
const bool NL, std::vector<double> &lgkk, std::vector<double> &lgPk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_max=100.)
const;
961 void m_Table_Pk_CLASS (
const bool NL, std::vector<std::vector<double>> &lgkk, std::vector<std::vector<double>> &lgPk,
const std::vector<double> redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_max=100.)
const;
994 void m_Table_Pk_parameterFile (
const std::string code,
const std::string file_par,
const bool NL, std::vector<double> &lgkk, std::vector<double> &lgPk,
const double redshift,
const std::string output_root=
"test")
const;
1069 Cosmology (
const double Omega_matter=0.27,
const double Omega_baryon=0.046,
const double Omega_neutrinos=0.,
const double massless_neutrinos=3.04,
const int massive_neutrinos=0,
const double Omega_DE=0.73,
const double Omega_radiation=0.,
const double hh=0.7,
const double scalar_amp=2.46e-9,
const double scalar_pivot=0.05,
const double n_spec=0.96,
const double w0=-1.,
const double wa=0.,
const double fNL=0.,
const int type_NG=1,
const double tau=0.09,
const std::string
model=
"LCDM",
const bool unit=
true);
1369 void set_parameters (
const std::vector<CosmologicalParameter> parameter,
const std::vector<double>
value);
1516 void set_hh (
const double hh=0.7,
const bool warn=
true)
1518 if (
m_unit && warn)
WarningMsgCBL(
"if unit=true then H0=100 (by internal definition)",
"set_hh",
"Cosmology.h");
1533 void set_H0 (
const double H0=70.,
const bool warn=
true)
1535 if (
m_unit && warn)
WarningMsgCBL(
"if unit=true then H0=100 (by internal definition)",
"set_H0",
"Cosmology.h");
1655 double OmegaM (
const double redshift=0.)
const;
1664 double OmegaDE (
const double redshift=0.)
const;
1673 double OmegaR (
const double redshift=0.)
const;
1692 double OmegaNu (
const double redshift=0.)
const;
1701 double OmegaK (
const double redshift=0.)
const;
1710 double Omega (
const double redshift=0.)
const;
1752 double rho_crit (
const double redshift,
const bool unit1=
false)
const;
1770 double rho_m (
const double redshift=0.,
const bool unit1=
false,
const bool nu=
false)
const;
1816 double Delta_c (
const double redshift,
const std::string author=
"BryanNorman")
const;
1844 double M_vir (
const double r_vir,
const double redshift,
const std::string author=
"BryanNorman",
const bool unit1=
false)
const;
1871 double r_vir (
const double M_vir,
const double redshift,
const std::string author=
"BryanNorman",
const bool unit1=
false)
const;
1882 double w_CPL (
const double redshift=0.)
const;
1895 double f_DE (
const double redshift=0.)
const;
1902 double EE (
const double redshift=0.)
const;
1910 double EE_inv (
const double redshift=0.)
const
1911 {
return 1./
EE(redshift); }
1920 {
return 1./(1.+redshift)/
EE(redshift); }
1930 double redshift = 1./aa-1.;
1931 return (1.+redshift)/
EE(redshift);
1939 double HH (
const double redshift=0.)
const;
1965 double DN (
const double redshift,
const double redshift_norm=0.,
const double prec=1.e-4)
const;
1987 double DD (
const double redshift)
const;
2005 double gg (
const double redshift=0.)
const;
2014 double sigma8 (
const double redshift)
const;
2028 double cosmic_time (
const double redshift=0.)
const;
2037 double EE2 (
const double redshift=0.)
const;
2044 double qq (
const double redshift=0.)
const;
2051 double Hdot (
const double redshift=0.)
const;
2078 double z_acc ()
const;
2087 double z_eq ()
const;
2096 double z_eq_rad (
const double T_CMB=2.7255)
const;
2108 double sound_speed (
const double redshift,
const double T_CMB=2.7255)
const;
2120 double rs_integrand (
const double redshift,
const double T_CMB=2.7255)
const;
2132 double rs (
const double redshift,
const double T_CMB=2.7255)
const;
2150 double Lum_bol (
const double redshift=0.,
const double flux=1.)
const;
2169 double Redshift (
const double d_c=1.,
const double z1_guess=0.,
const double z2_guess=10.,
const double prec=0.0001)
const;
2204 double Redshift_LCDM (
const double d_c=1.,
const double z1_guess=0.,
const double z2_guess=10.,
const bool go_fast=1,
const double prec=0.0001)
const;
2237 double Redshift (
const double mm,
const double redshift,
const double ff,
const std::string method_SS,
const double wwf,
const bool store_output=
true,
const std::string output_root=
"test")
const;
2246 double Redshift_time (
const double time,
const double z1_guess,
const double z2_guess)
const;
2262 double deltac (
const double redshift)
const;
2289 double deltav_L (
const double deltav_NL,
const double b_eff,
double slope=0.854,
double offset=0.420)
const;
2308 double deltav_NL (
const double deltav=-2.71)
const;
2322 double r_rL (
const double deltav = -2.71)
const;
2340 double D_C (
const double redshift)
const;
2392 double D_C_LCDM (
const double redshift)
const;
2412 void D_C_table (
const std::string file_table,
const double z_min,
const double z_max,
const int step, std::vector<double> &
Redshift, std::vector<double> &dc)
const;
2419 double D_M (
const double redshift)
const;
2426 double D_A (
const double redshift)
const;
2451 double D_A (
const double z1,
const double z2)
const;
2458 double D_L (
const double redshift)
const;
2466 double D_V (
const double redshift)
const;
2473 double F_AP (
const double redshift)
const;
2488 double Distance (
const double redshift,
const std::string distance_type)
const;
2507 double Volume (
const double z1,
const double z2,
const double Area)
const;
2532 double Volume (
const double z1,
const double z2,
const double RA_min,
const double RA_max,
const double Dec_min,
const double Dec_max)
const;
2548 double Volume (
const double zz)
const;
2558 double max_redshift (
const double Volume,
const double Area,
const double z_min)
const;
2571 double dV_dZdOmega (
const double redshift,
const bool angle_rad)
const;
2624 std::string
create_grid_sigmaM (
const std::string method_SS,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true)
const;
2724 double mass_function (
const double Mass,
const double redshift,
const std::string model_MF,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const bool default_delta=
true,
const double delta_t=1.686);
2814 double mass_function_fR (
const double Mass,
const double redshift,
const std::string model_MF,
const double f_R0=0.,
const bool store_output=
true,
const double Delta=200.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const bool default_delta=
true,
const double delta_t=1.686);
2908 double mass_function_fast (
const double Mass,
const double redshift,
const std::string model_MF,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
3007 double mass_function (
const double Mass,
const double Sigma,
const double Dln_Sigma,
const double redshift,
const std::string model_MF,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string method_SS=
"CAMB",
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
3108 double mass_function (
const double Mass,
const double Sigma,
const double Dln_Sigma,
const double redshift,
const double D_N,
const std::string model_MF,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string method_SS=
"CAMB",
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
3208 std::vector<double>
mass_function (
const std::vector<double>
Mass,
const std::vector<double>
Sigma,
const std::vector<double> Dln_Sigma,
const double redshift,
const std::string model_MF,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string method_SS=
"CAMB",
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
3296 double n_haloes (
const double Mass_min,
const double Mass_max,
const double z_min,
const double z_max,
const bool angle_rad,
const std::string model_MF,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200,
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
3407 double n_haloes (
const double Mass_min,
const double Mass_max,
const double Volume,
const double redshift,
const std::string model_MF,
const std::string method_SS,
const int nbin_mass=0,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const bool default_delta=
true,
const double delta_t=1.686);
3504 double n_haloes_selection_function (
const double Mass_min,
const double Mass_max,
const double z_min,
const double z_max,
const bool angle_rad,
const std::string model_MF,
const std::string method_SS,
const std::string selection_function_file,
const std::vector<int> column={},
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200,
const bool isDelta_critical=
false,
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
3587 std::vector<double>
mass_function (
const std::vector<double> mass,
const double z_min,
const double z_max,
const std::string model_MF,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200,
const bool isDelta_critical=
false,
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
3675 std::vector<double>
mass_function_selection_function_vector (
const std::vector<double> mass,
const double z_min,
const double z_max,
const std::string model_MF,
const std::string method_SS,
const std::string selection_function_file,
const std::vector<int> column={},
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200,
const bool isDelta_critical=
false,
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
3764 std::vector<double>
redshift_distribution_haloes (
const double z_min,
const double z_max,
const int step_z,
const double Area_degrees,
const double Mass_min,
const double Mass_max,
const std::string model_MF,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200,
const bool isDelta_critical=
false,
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
3856 std::vector<double>
redshift_distribution_haloes_selection_function (
const std::vector<double> redshift,
const double Area_degrees,
const double Mass_min,
const double Mass_max,
const std::string model_MF,
const std::string method_SS,
const std::string selection_function_file,
const std::vector<int> column={},
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200,
const bool isDelta_critical=
false,
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
3947 double mean_redshift_haloes_selection_function (
const double z_min,
const double z_max,
const double Mass_min,
const double Mass_max,
const std::string model_MF,
const std::string method_SS,
const std::string selection_function_file,
const std::vector<int> column={},
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200,
const bool isDelta_critical=
false,
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
4043 double MhaloMin (
const int n_halo,
const double Area,
const bool angle_rad,
const double z_min,
const double z_max,
const double Mmax,
const double lgM1_guess,
const double lgM2_guess,
const std::string model_MF,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200,
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true)
const;
4074 double converted_mass (
const double mass,
const cosmology::Cosmology cosmology,
const double redshift,
const double redshift_source=-1.)
const;
4100 double pw (
const double ww,
const double ff,
const std::string author)
const;
4134 double pz (
const double m0,
const double z0,
const double frac,
const double redshift,
const std::string model_model,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test")
const;
4150 double cumPw (
const double ww,
const double ff,
const std::string author)
const;
4161 void medianwf (
const double ff,
const std::string model_model, std::vector<double> &
wf)
const;
4190 void medianzf (
const double ff,
const double mass,
const double z0,
const std::string model_model,
const std::string method_SS, std::vector<double> &zf,
const bool store_output=
true,
const std::string output_root=
"test")
const;
4223 double wf (
const double mm,
const double redshift,
const double ff,
const double zf,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test")
const;
4254 double As (
const double sigma8)
const;
4303 std::string
Pk_output_file (
const std::string code,
const bool NL,
const double redshift,
const bool run=0,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_max=100.,
const std::string file_par=
par::defaultString);
4363 void run_CAMB (std::vector<double> &lgkk, std::vector<double> &lgPk,
const bool NL,
const double redshift,
const std::string output_root=
"test",
const std::string output_dir=
par::defaultString,
const double k_max=100.,
const std::string file_par=
par::defaultString)
const;
4398 void Table_PkCodes (
const std::string code,
const bool NL, std::vector<double> &lgkk, std::vector<double> &lgPk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_max=100.,
const std::string file_par=
par::defaultString)
const;
4435 void Table_PkCodes (
const std::string code,
const bool NL, std::vector<std::vector<double>> &lgkk, std::vector<std::vector<double>> &lgPk,
const std::vector<double> redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_max=100.,
const std::string file_par=
par::defaultString)
const;
4470 void Table_XiCodes (
const std::string code,
const bool NL, std::vector<double> &rr, std::vector<double> &xi,
const double redshift,
const bool store_output,
const std::string output_root,
const double k_max, std::string file_par)
const;
4510 void Pk_0 (
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
4565 std::vector<double>
Pk_matter (
const std::vector<double> kk,
const std::string method_Pk,
const bool NL,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString,
const bool unit1=
false);
4620 std::vector<std::vector<double>>
Pk_matter (
const std::vector<double> kk,
const std::string method_Pk,
const bool NL,
const std::vector<double> redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString,
const bool unit1=
false);
4658 std::vector<double>
Pk_matter_Linear (
const std::string method,
const std::vector<double> kk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const bool norm=1,
const double prec=1.e-4);
4704 std::vector<double>
Pk_matter_NoWiggles_gaussian (
const std::vector<double> kk,
const std::vector<double> PkLin,
const std::vector<double> PkApprox,
const double lambda,
const std::string method);
4732 std::vector<double>
Pk_matter_NoWiggles_bspline (
const std::vector<double> kk,
const std::vector<double> PkLin,
const std::vector<double> PkApprox,
const int order,
const int nknots);
4791 std::vector<double>
Pk_matter_NoWiggles (
const std::string method,
const std::vector<double> kk,
const double redshift,
const std::string linear_method=
"CAMB",
const int order=4,
const int nknots=10,
const double lambda=0.25,
const bool store_output=
true,
const std::string output_root=
"test",
const bool norm=1,
const double prec=1.e-4);
4842 std::vector<double>
Pk_matter_DeWiggled (
const std::string linear_method,
const std::string nowiggles_method,
const std::vector<double> kk,
const double redshift,
const double sigma_NL,
const int order=4,
const int nknots=10,
const double lambda=0.25,
const bool store_output=
true,
const std::string output_root=
"test",
const bool norm=1,
const double prec=1.e-4);
4863 void remove_output_Pk_tables (
const std::string code,
const bool NL,
const double redshift,
const std::string output_root=
"test")
const;
4920 double sigma2R (
const double radius,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const bool unit1=
false)
const;
4977 double sigma2M (
const double mass,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const bool unit1=
false)
const;
5034 double dnsigma2R (
const int nd,
const double radius,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const bool unit1=
false)
const;
5091 double dnsigma2M (
const int nd,
const double mass,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const bool unit1=
false)
const;
5159 double c_vir (
const double c200,
const double redshift,
const std::string author=
"BryanNorman")
const;
5228 double xi_matter (
const double rr,
const std::string method_Pk,
const bool NL,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double aa=0.,
const bool GSL=
false,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
5300 double wtheta_DM (
const double theta,
const std::vector<double> zz,
const std::vector<double> phiz,
const std::string interpolationMethod,
const CoordinateUnits coordUnits =
CoordinateUnits::_degrees_,
const bool GSL=
false,
const std::string method_Pk=
"CAMB",
const bool NL=
false,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=1.e-4,
const double k_max=100,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
5352 double wtheta_DM (
const double theta,
const std::vector<double> kk,
const std::vector<double> Pk,
const std::vector<double> zz,
const std::vector<double> nz,
const std::vector<double> phiz,
const std::string interpolationType=
"Spline",
const CoordinateUnits coordUnits =
CoordinateUnits::_degrees_,
const bool GSL=
false,
const double redshift_Pk=0);
5411 std::vector<double>
C_l_DM (
const int lmax,
const std::vector<double> zz,
const std::vector<double> phiz,
const std::string interpolationMethod,
const std::string method_Pk=
"CAMB",
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=1.e-4,
const double k_max=100,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
5461 double xi_matter_DeWiggle (
const double rr,
const double redshift,
const double sigma_NL,
const bool store_output=
true,
const std::string output_root=
"test",
const bool norm=1,
const double k_min=0.001,
const double k_max=100.,
const double aa=1.,
const double prec=1.e-2);
5532 void get_xi (std::vector<double> &rr, std::vector<double> &Xi,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const bool xiType=0,
const double k_star=-1.,
const bool xiNL=0,
const int norm=-1,
const double r_min=0.1,
const double r_max=150.,
const double k_min=0.001,
const double k_max=100.,
const double aa=0.,
const bool GSL=
false,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
5596 void get_barred_xi (std::vector<double> rr, std::vector<double> Xi, std::vector<double> &Xi_, std::vector<double> &Xi__,
const std::string method_Pk,
const double redshift,
const bool xiType=0,
const double k_star=-1.,
const bool xiNL=0,
const int norm=-1,
const double r_min=0.1,
const double r_max=150.,
const double k_min=0.001,
const double k_max=100.,
const double aa=0.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString)
const;
5664 double wp_DM (
const double rp,
const std::string method_Pk,
const bool NL,
const double redshift,
const double pimax,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double r_min=1.e-3,
const double r_max=350.,
const double k_min=0.001,
const double k_max=100.,
const double aa=0.,
const bool GSL=
false,
const double prec=1.e-2,
const std::string file_par=
cbl::par::defaultString);
5700 double k_star (
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_max=100.,
const std::string file_par=
par::defaultString)
const;
5769 double sigmaR_DM (
const double RR,
const int corrType,
const std::string method_Pk,
const double redshift,
const double pimax=40,
const bool store_output=
true,
const std::string output_root=
"test",
const bool NL=1,
const int norm=-1,
const double r_min=1.e-3,
const double r_max=350.,
const double k_min=0.001,
const double k_max=100.,
const double aa=0.,
const bool GSL=
false,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
5815 double sigma8_Pk (
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const bool NL=0,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString)
const;
5887 double bias_halo (
const double Mass,
const double redshift,
const std::string author,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double Delta=200.,
const double kk=-1.,
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
5961 double bias_halo (
const double Mass,
const double Sigma,
const double redshift,
const std::string model_bias,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double Delta=200.,
const double kk=-1.,
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string method_SS=
"CAMB",
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
6038 double bias_halo (
const double Mass,
const double Sigma,
const double redshift,
const double DN,
const std::string model_bias,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double Delta=200.,
const double kk=-1.,
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string method_SS=
"CAMB",
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
6112 std::vector<double>
bias_halo (
const std::vector<double>
Mass,
const std::vector<double>
Sigma,
const double redshift,
const std::string model_bias,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double Delta=200.,
const double kk=-1.,
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string method_SS=
"CAMB",
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
6226 double bias_eff (
const double Mass_min,
const double Mass_max,
const double redshift,
const std::string model_bias,
const std::string model_MF,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
6319 double bias_eff (
const std::vector<double> MM,
const std::vector<double> MF,
const double redshift,
const std::string model_bias,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
6415 std::vector<double>
bias_eff_mass_grid (
const std::vector<double> MM,
const std::vector<double> redshift,
const std::string model_bias,
const std::string method_SS,
const std::string meanType=
"mean_bias",
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta_crit=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
6507 std::vector<double>
bias_eff_mass (
const std::vector<double> MM,
const std::vector<double> redshift,
const std::string model_bias,
const std::string method_SS,
const std::string meanType=
"mean_bias",
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
6603 std::vector<double>
bias_eff_mass (
const std::vector<double> mass,
const std::vector<double> mass_grid,
const std::vector<double> redshift,
const std::string model_bias,
const std::string method_SS,
const std::string meanType=
"mean_bias",
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
6730 void generate_bias_eff_grid_one_cosmopar (std::vector<double> ¶meter, std::vector<double> &
bias_eff,
const std::string dir_output,
const std::string file_bias_eff_grid,
const cbl::cosmology::CosmologicalParameter cosmoPar,
const double min_par,
const double max_par,
const int nbin_par,
const std::vector<double> mass,
const std::vector<double> mass_grid,
const std::vector<double> redshift,
const std::string model_bias,
const std::string method_SS,
const std::string meanType=
"mean_bias",
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta_crit=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const cbl::cosmology::Cosmology cosmology_mass={},
const std::vector<double> redshift_source={});
6861 void generate_bias_eff_grid_one_cosmopar (std::vector<double> ¶meter, std::vector<double> &
bias_eff,
const std::string dir_output,
const std::string file_bias_eff_grid,
const cbl::cosmology::CosmologicalParameter cosmoPar,
const double min_par,
const double max_par,
const int nbin_par,
const double redshift,
const double Mass_min,
const double Mass_max,
const std::string model_bias,
const std::string model_MF,
const std::string method_SS,
const std::string selection_function_file,
const std::vector<int> column={},
const double alpha=1.,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta_crit=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
7003 void generate_bias_eff_grid_two_cosmopars (std::vector<double> ¶meter1, std::vector<double> ¶meter2, std::vector<std::vector<double>> &
bias_eff,
const std::string dir_output,
const std::string file_bias_eff_grid,
const cbl::cosmology::CosmologicalParameter cosmoPar1,
const double min_par1,
const double max_par1,
const int nbin_par1,
const cbl::cosmology::CosmologicalParameter cosmoPar2,
const double min_par2,
const double max_par2,
const int nbin_par2,
const std::vector<double> mass,
const std::vector<double> mass_grid,
const std::vector<double> redshift,
const std::string model_bias,
const std::string method_SS,
const std::string meanType=
"mean_bias",
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true,
const cbl::cosmology::Cosmology cosmology_mass={},
const std::vector<double> redshift_source={});
7129 std::vector<double>
bias_eff_selection_function (
const glob::FuncGrid interp_sigma,
const glob::FuncGrid interp_DnSigma,
const glob::FuncGrid interp_SF,
const double Mass_min,
const double Mass_max,
const std::vector<double> redshift,
const std::string model_bias,
const std::string model_MF,
const std::string method_SS,
const double alpha=1.,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta_crit=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
7255 std::vector<double>
bias_eff_selection_function (
const glob::FuncGrid interp_sigma,
const glob::FuncGrid interp_DnSigma,
const glob::FuncGrid2D interp_SF,
const double Mass_min,
const double Mass_max,
const std::vector<double> redshift,
const std::string model_bias,
const std::string model_MF,
const std::string method_SS,
const double alpha=1.,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta_crit=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
7370 std::vector<double>
bias_eff_selection_function (
const double Mass_min,
const double Mass_max,
const std::vector<double> redshift,
const std::string model_bias,
const std::string model_MF,
const std::string method_SS,
const std::string selection_function_file,
const std::vector<int> column={},
const double alpha=1.,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta_crit=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
7445 double fsigma8 (
const double redshift,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const bool NL=
false,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString)
const;
7454 double beta (
const double redshift,
const double bias)
const;
7464 double error_beta (
const double redshift,
const double bias,
const double err_bias)
const;
7559 double beta (
const double Mass_min,
const double Mass_max,
const double redshift,
const std::string model_bias,
const std::string model_MF,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
7655 double error_beta (
const double Mass_min,
const double Mass_max,
const double redshift,
const std::string model_bias,
const std::string model_MF,
const std::string method_SS,
const double err_bias,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
7735 double beta (
const std::vector<double> MM,
const std::vector<double> MF,
const double redshift,
const std::string model_bias,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
7816 double error_beta (
const std::vector<double> MM,
const std::vector<double> MF,
const double redshift,
const std::string model_bias,
const std::string method_SS,
const double err_bias,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
7916 double error_beta_measured (
const double Volume,
const double density,
const double Mass_min,
const double Mass_max,
const double redshift,
const std::string model_bias,
const std::string model_MF,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
8010 double quadrupole (
const double Mass_min,
const double Mass_max,
const double redshift,
const std::string model_bias,
const std::string model_MF,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
8084 double quadrupole (
const std::vector<double> MM,
const std::vector<double> MF,
const double redshift,
const std::string model_bias,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const double kk=-1.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
8127 double square_bulk_flow (
const double rr,
const double k_int_min,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
8145 double square_bulk_flow_Table (
const double rr,
const double k_int_min,
const std::vector<double> lgkk,
const std::vector<double> lgPk,
const double redshift)
const;
8188 double square_velocity_dispersion (
const double rr,
const double k_int_min,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
8224 double CMN (
const double rr,
const double k_int_min,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_max=100.,
const std::string file_par=
par::defaultString)
const;
8275 double Sn_PT (
const int nn,
const double RR,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true)
const;
8327 double Sigman_PT (
const int nn,
const double RR,
const std::string method_SS,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true)
const;
8413 double xi0_Kaiser (
const double rad,
const double f_sigma8,
const double bias_sigma8,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const bool xiType=0,
const double k_star=-1.,
const bool NL=
false,
const int norm=-1,
const double r_min=0.1,
const double r_max=150.,
const double k_min=0.001,
const double k_max=100.,
const double aa=0.,
const bool GSL=
false,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
8478 std::vector<double>
xi0_Kaiser (
const std::vector<double> rad,
const double bias,
const std::string method_Pk,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const bool NL=
false,
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
8589 double xi2D_dispersionModel (
const double rp,
const double pi,
const double f_sigma8,
const double bias_sigma8,
const double sigmav,
const std::string method_Pk,
const double redshift,
const int FV,
const bool NL, std::vector<double> rr, std::vector<double> &Xi, std::vector<double> &Xi_, std::vector<double> &Xi__,
const bool store_output=
true,
const std::string output_root=
"test",
const int index=-1,
const bool bias_nl=0,
const double bA=-1.,
const bool xiType=0,
const double k_star=-1.,
const bool xiNL=0,
const double v_min=-3000.,
const double v_max=3000.,
const int step_v=500,
const int norm=-1,
const double r_min=0.1,
const double r_max=150.,
const double k_min=0.001,
const double k_max=100.,
const double aa=0.,
const bool GSL=
false,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
8628 double xi_star (
const double rr,
const double redshift,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_star=-1.,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
8671 double xisnl_gnw (
const double rp,
const double pi,
const double f_sigma8,
const double bias_sigma8,
const double bA,
const double redshift, std::vector<double> rr, std::vector<double> Xi, std::vector<double> &Xi_, std::vector<double> &Xi__,
const bool store_output=
true,
const std::string output_root=
"test");
8723 double xis_gBAO (
const double rp,
const double pi,
const double f_sigma8,
const double bias_sigma8,
const double redshift, std::vector<double> rr, std::vector<double> Xi, std::vector<double> &Xi_, std::vector<double> &Xi__,
const bool store_output=
true,
const std::string output_root=
"test",
const double k_star=-1.,
const double x_min=-3000.,
const double x_max=3000.,
const int step_x=500);
8840 double xi2D_CW (
const double rp,
const double pi,
const double beta,
const double bias_lin,
const double bA,
const double sigmav0,
const double cmu,
const double cs1,
const double cs2,
const double redshift, std::vector<double> rr1, std::vector<double> Xi1, std::vector<double> rr2, std::vector<double>
Xi2, std::vector<double> &Xi1_, std::vector<double> &Xi1__, std::vector<double> &Xi2_, std::vector<double> &Xi2__,
const bool store_output=
true,
const std::string output_root=
"test",
const bool BAO=1,
const bool xiType=0,
const double k_star=-1,
const bool xiNL=0,
const double r_min=0.1,
const double r_max=150.,
const double v_min=-3000.,
const double v_max=3000.,
const int step_v=500,
const double k_min=0.001,
const double k_max=100.,
const double x_min=-3000.,
const double x_max=3000.,
const int step_x=500,
const double aa=0.,
const bool GSL=
false,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
8866 double rs (
const std::string method_Pk,
const double T_CMB=
par::TCMB)
const;
8913 double ys (
const double redshift,
const std::string method_Pk,
const double T_CMB=
par::TCMB)
const;
8925 double Az (
const double redshift)
const;
8950 std::vector<double>
linear_point (
const double redshift,
const double rmin=60.,
const double rmax=150.,
const int nbinr=100,
const std::string interpType=
"Spline");
9004 double Am (
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
9050 double potential_spectral_amplitude (
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
9098 double bispectrum (
const std::vector<double> kk,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
9149 double mrk (
const double kk,
const double mass,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString);
9213 double frk (
const double kk,
const double mass,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
9216 double bias_kernel (
const double,
void *);
9218 double frk_test (
const double,
const double,
const std::string,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
9283 double bias_correction (
const double kk,
const double mass,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
9344 double skewness (
const double mass,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
9405 double dskewnessdM (
const double mass,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
9467 double MF_correction (
const double mass,
const double redshift,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
9493 double f_nu (
const double SS,
const double del_v,
const double del_c)
const;
9563 double size_function (
const double RV,
const double redshift,
const std::string
model,
const double b_eff,
double slope=0.854,
double offset=0.420,
const double deltav_NL=-0.795,
const double del_c=1.69,
const std::string method_Pk=
"EisensteinHu",
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true)
const;
9633 std::vector<double>
size_function (
const std::vector<double> RV,
const double redshift,
const std::string
model,
const double b_eff,
double slope=0.854,
double offset=0.420,
const double deltav_NL=-0.795,
const double del_c=1.69,
const std::string method_Pk=
"EisensteinHu",
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true)
const;
9713 std::vector<std::vector<double>>
Nvoids (
const double min_r,
const double max_r,
const int num_bins,
const double mean_z,
const double Volume,
const std::string
model,
const double b_eff,
double slope=0.854,
double offset=0.420,
const double deltav_NL=-0.795,
const double del_c=1.69,
const std::string method_Pk=
"EisensteinHu",
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true)
const;
9797 std::vector<std::vector<double>>
Nvoids (
const double min_r,
const double max_r,
const int num_bins,
const double min_z,
const double max_z,
const double mean_z,
const double Area,
const std::string
model,
const double b_eff,
double slope=0.854,
double offset=0.420,
const double deltav_NL=-0.795,
const double del_c=1.69,
const std::string method_Pk=
"EisensteinHu",
const bool store_output=
true,
const std::string output_root=
"test",
const std::string interpType=
"Linear",
const double k_max=100.,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true)
const;
9885 double size_function (
const double RV,
const double redshift,
const std::string model_mf,
const double del_v,
const std::string model_sf,
const std::string method_Pk=
"EisensteinHu",
const bool store_output=
true,
const std::string output_root=
"test",
const double Delta=200.,
const std::string interpType=
"Linear",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string input_file=
par::defaultString,
const bool is_parameter_file=
true);
9944 std::vector<std::vector<double>>
XiMultipoles (
const int nbins,
const double rMin,
const double rMax,
const std::vector<double> kk,
const std::vector<double> Pk0,
const std::vector<double> Pk2,
const std::vector<double> Pk4,
const int IntegrationMethod=1);
9971 std::vector<std::vector<double>>
XiMonopole_covariance (
const int nbins,
const double rMin,
const double rMax,
const double nn,
const double Volume,
const std::vector<double> kk,
const std::vector<double> Pk0,
const int IntegrationMethod=1);
10004 std::vector<std::vector<double>>
XiMultipoles_covariance (
const int nbins,
const double rMin,
const double rMax,
const double nn,
const double Volume,
const std::vector<double> kk,
const std::vector<double> Pk0,
const std::vector<double> Pk2,
const std::vector<double> Pk4,
const int IntegrationMethod=1);
10032 double F2 (
const double k,
const double q,
const double kq);
10054 double G2 (
const double k,
const double q,
const double kq);
10083 double f_k (
const double k,
const std::shared_ptr<cbl::glob::FuncGrid> PkLin,
const double qmin,
const double qmax,
const double prec=1.e-3);
10112 double g_k (
const double k,
const std::shared_ptr<cbl::glob::FuncGrid> PkLin,
const double qmin,
const double qmax,
const double prec=1.e-3);
10151 double Pk_1loop (
const double kk,
const std::shared_ptr<cbl::glob::FuncGrid> PkLin,
const int corrtype,
const double qmin,
const double qmax,
const double prec=1.e-3);
10185 double Pk_DeltaDelta (
const double kk,
const std::shared_ptr<cbl::glob::FuncGrid> Pk,
const double qmin,
const double qmax,
const double prec=1.e-3);
10250 std::vector<double>
Pk_DeltaDelta (
const std::vector<double> kk,
const double redshift,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString,
const bool unit1=
false);
10284 double Pk_DeltaTheta (
const double kk,
const std::shared_ptr<cbl::glob::FuncGrid> Pk,
const double qmin,
const double qmax,
const double prec=1.e-3);
10349 std::vector<double>
Pk_DeltaTheta (
const std::vector<double> kk,
const double redshift,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString,
const bool unit1=
false);
10382 double Pk_ThetaTheta (
const double kk,
const std::shared_ptr<cbl::glob::FuncGrid> Pk,
const double qmin,
const double qmax,
const double prec=1.e-3);
10447 std::vector<double>
Pk_ThetaTheta (
const std::vector<double> kk,
const double redshift,
const std::string method_Pk,
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2,
const std::string file_par=
par::defaultString,
const bool unit1=
false);
10509 double Pk_DeltaDelta_fitting_function (
const double kk,
const std::string method_Pk,
const double redshift,
const std::string author,
const bool store_output,
const std::string output_root,
const int norm,
double k_min,
double k_max,
const double prec,
const std::string file_par,
const bool unit1);
10578 double Pk_DeltaTheta_fitting_function (
const double kk,
const std::string method_Pk,
const double redshift,
const std::string author,
const bool store_output,
const std::string output_root,
const bool NL,
const int norm,
double k_min,
double k_max,
const double prec,
const std::string file_par,
const bool unit1);
10651 double Pk_ThetaTheta_fitting_function (
const double kk,
const std::string method_Pk,
const double redshift,
const std::string author,
const bool store_output,
const std::string output_root,
const bool NL,
const int norm,
double k_min,
double k_max,
const double prec,
const std::string file_par,
const bool unit1);
10713 double sigma_v (
const double redshift=0.,
const std::string method_Pk=
"CAMB",
const bool store_output=
true,
const std::string output_root=
"test",
const int norm=-1,
const double k_min=0.001,
const double k_max=100.,
const int bin_k=512,
const double prec=1.e-2,
const std::string file_par=
par::defaultString,
const bool unit1=
false);
10775 std::vector<std::vector<double>>
Pk_TNS_AB_multipoles (std::vector<double> kk,
const std::string method,
const double redshift,
const bool store_output,
const std::string output_root,
const int norm,
const double k_min,
const double k_max,
const double prec);
10838 std::vector<std::vector<double>>
Pk_TNS_AB_1loop (std::vector<double> kk,
const double mu,
const std::string method,
const double redshift,
const bool store_output,
const std::string output_root,
const int norm,
const double k_min,
const double k_max,
const double prec);
10898 std::vector<std::vector<double>>
Pk_TNS_AB_terms_1loop (std::vector<double> kk,
const std::string method,
const double redshift,
const bool store_output,
const std::string output_root,
const int norm,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2);
10972 std::vector<std::vector<double>>
Pk_eTNS_terms_1loop (std::vector<double> kk,
const std::string method,
const double redshift,
const bool store_output,
const std::string output_root,
const int norm,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2);
11039 std::vector<std::vector<double>>
Pk_TNS_AB_1loop (std::vector<double> kk,
const double mu,
const double linear_growth_rate,
const double bias,
const std::string method,
const double redshift,
const bool store_output,
const std::string output_root,
const int norm,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2);
11105 std::vector<std::vector<double>>
Pk_TNS_dd_dt_tt (std::vector<double> kk,
const std::string method,
const double redshift,
const bool store_output,
const std::string output_root,
const int norm,
const double k_min=0.001,
const double k_max=100.,
const double prec=1.e-2);
11144 double denominator_Q (
const double r1,
const double r2,
const double theta,
const std::vector<double> rr,
const std::vector<double>
xi_matter)
const;
11183 void integrals_Q_nonLocal (std::vector<double> &
xi_matter, std::vector<double> &Phi,
const std::vector<double> rr,
const std::vector<double> kk,
const std::vector<double>
Pk_matter,
const double prec)
const;
11214 double Gamma_3PCF (
const double r1,
const double r2,
const double theta,
const std::vector<double> xi,
const std::vector<double> dPhi)
const;
11258 double Q_nonLocal (
const double r1,
const double r2,
const double theta, std::vector<double> &rr, std::vector<double> &
xi_matter, std::vector<double> &Phi,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11284 std::vector<double>
Q_nonLocal (
const double r1,
const double r2,
const std::vector<double> theta,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11319 void integrals_zeta_Slepian (std::vector<double> &
xi_matter, std::vector<double> &xi_matter_m1, std::vector<double> &xi_matter_p1, std::vector<double> &xi_matter_2,
const std::vector<double> rr,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11379 double zeta_precyclic_Slepian (
const double r1,
const double r2,
const double mu,
const double b1,
const double b2,
const glob::FuncGrid interp_xi_matter,
const glob::FuncGrid interp_xi_matter_m1,
const glob::FuncGrid interp_xi_matter_p1,
const glob::FuncGrid interp_xi_matter_2)
const;
11443 double zeta_precyclic_Slepian (
const double r1,
const double r2,
const double r3,
const double deltaR,
const double b1,
const double b2,
const glob::FuncGrid interp_xi_matter,
const glob::FuncGrid interp_xi_matter_m1,
const glob::FuncGrid interp_xi_matter_p1,
const glob::FuncGrid interp_xi_matter_2)
const;
11489 std::vector<double>
zeta_expansion_Slepian (
const double r1,
const double r2,
const double b1,
const double b2, std::vector<double> &rr, std::vector<double> &
xi_matter, std::vector<double> &xi_matter_m1, std::vector<double> &xi_matter_p1, std::vector<double> &xi_matter_2,
const int norders=9,
const double prec=1.e-3)
const;
11535 double zeta_DM_Slepian (
const double r1,
const double r2,
const double theta, std::vector<double> &rr, std::vector<double> &
xi_matter, std::vector<double> &xi_matter_m1, std::vector<double> &xi_matter_p1, std::vector<double> &xi_matter_2,
const std::vector<double> kk,
const std::vector<double>
Pk_matter,
const int norders=9,
const double prec=1.e-3)
const;
11584 double Q_DM_Slepian (
const double r1,
const double r2,
const double theta, std::vector<double> &rr, std::vector<double> &
xi_matter, std::vector<double> &xi_matter_m1, std::vector<double> &xi_matter_p1, std::vector<double> &xi_matter_2,
const std::vector<double> kk,
const std::vector<double>
Pk_matter,
const int norders=9,
const double prec=1.e-3)
const;
11652 double zeta_single_BarrigaGatzanaga (
const double r1,
const double r2,
const double theta,
const std::vector<double> xi,
const std::vector<double> dxi,
const std::vector<double> dPhi)
const;
11696 double zeta_DM_BarrigaGatzanaga (
const double r1,
const double r2,
const double theta, std::vector<double> &rr, std::vector<double> &
xi_matter, std::vector<double> &Phi,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11735 double Q_DM_BarrigaGatzanaga (
const double r1,
const double r2,
const double theta, std::vector<double> &rr, std::vector<double> &
xi_matter, std::vector<double> &Phi,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11762 std::vector<double>
zeta_DM (
const double r1,
const double r2,
const std::vector<double> theta,
const std::string
model,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11790 std::vector<double>
Q_DM (
const double r1,
const double r2,
const std::vector<double> theta,
const std::string
model,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11832 std::vector<double>
zeta_halo (
const double r1,
const double r2,
const std::vector<double> theta,
const double b1,
const double b2,
const std::string
model,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11872 std::vector<double>
Q_halo (
const double r1,
const double r2,
const std::vector<double> theta,
const double b1,
const double b2,
const std::string
model,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11916 std::vector<double>
Q_halo (
const double r1,
const double r2,
const std::vector<double> theta,
const double b1,
const double b2,
const double g2,
const std::string
model,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11938 std::vector<double>
zeta_DM_eq (
const std::vector<double> rr,
const std::string
model,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
11961 std::vector<double>
Q_DM_eq (
const std::vector<double> rr,
const std::string
model,
const std::vector<double> kk,
const std::vector<double>
Pk_matter)
const;
12025 double zeta_multipoles_covariance (
const double Volume,
const double nObjects,
const int l,
const int l_prime,
const double r1,
const double r2,
const double r1_prime,
const double r2_prime,
const double deltaR,
const std::vector<double> kk,
const std::vector<double> Pk,
const std::vector<double> rr,
const std::vector<double> Xi,
const double prec=1.e-3);
12078 std::vector<std::vector<double>>
zeta_covariance (
const double Volume,
const double nObjects,
const std::vector<double> theta,
const double r1,
const double r2,
const double deltaR,
const std::vector<double> kk,
const std::vector<double> Pk,
const int norders=10,
const double prec=1.e-3,
const bool method=
false,
const int nExtractions=10000,
const std::vector<double> mean={},
const int seed=543);
12101 void xi_r_n (std::vector<double> &xi_n,
const std::vector<double> rr,
const int nn,
const std::vector<double> kk,
const std::vector<double> Pk);
12126 void xi_r_n_pm (std::vector<double> &xi_n_p, std::vector<double> &xi_n_m,
const std::vector<double> rr,
const int nn,
const std::vector<double> kk,
const std::vector<double> Pk);
12151 void eff_l_l1 (std::vector<std::vector<double>> &eff,
const std::vector<double> rr,
const int l,
const int l1,
const std::vector<double> kk,
const std::vector<double> Pk);
12181 void I_ELL_ell (std::vector<std::vector<double>> &II,
const std::vector<double> rr,
const int ll,
const int LL,
const std::vector<double> kk,
const std::vector<double> Pk);
12208 void k_ell (std::vector<std::vector<double>> &KK,
const std::vector<double> rr,
const int ll,
const std::vector<double> kk,
const std::vector<double> Pk);
12444 double zeta_ell_precyclic (
const double r1,
const double r2,
const int ell,
const double b1,
const double b2,
const double bt,
const double beta, std::vector<std::shared_ptr<glob::FuncGrid>> interp_xi_ell,
const bool use_k, std::shared_ptr<glob::FuncGrid2D> interp_k_ell);
12481 std::vector<double>
zeta_RSD (
const double r1,
const double r2,
const int ntheta,
const double b1,
const double b2,
const double bt,
const double beta,
const std::vector<double> rr,
const std::vector<double> kk,
const std::vector<double> Pk,
const bool include_limits=
false,
const int max_ll=4,
const bool use_k=
false);
12539 std::vector<double>
zeta_RSD (
const double r1,
const double r2,
const int ntheta,
const double b1,
const double b2,
const double bt,
const double redshift,
const std::string method_Pk,
const int step_r,
const int step_k,
const bool store_output=
true,
const std::string output_root=
"test",
const bool force_RealSpace =
false,
const bool include_limits=
false,
const int max_ll=4,
const bool use_k=
false);
12554 double GSL_bias_kernel_wrapper (
double,
void *);
12556 double func_xi_EH_GSL (
double,
void *);
12557 double func_sigma2M_EH_GSL (
double,
void *);
12559 double bias_kernel2 (
double,
void *);
12560 double skewness_kernel (
double *,
size_t,
void *);
12563 struct GSL_f_pars {
12566 std::string method_Pk;
12568 std::string output_root;
12574 std::string file_par;
12575 cosmology::Cosmology *pt_Cosmology;
12580 double Omega_matter;
12581 double Omega_baryon;
12582 double Omega_neutrinos;
12583 double massless_neutrinos;
12584 double massive_neutrinos;
12586 double Omega_radiation;
12589 double scalar_pivot;
12598 std::string method_Pk;
12604 struct STR_sigma2M_EH
12606 double Omega_matter;
12607 double Omega_baryon;
12608 double Omega_neutrinos;
12609 double massless_neutrinos;
12610 double massive_neutrinos;
12612 double Omega_radiation;
12615 double scalar_pivot;
12624 std::string method_Pk;
12633 std::string model_MF;
12634 std::string method_SS;
12636 std::string output_root;
12638 std::string interpType;
12643 std::string input_file;
12644 bool is_parameter_file;
12645 bool default_delta;
12651 double Omega_matter;
12652 double Omega_baryon;
12653 double Omega_neutrinos;
12654 double massless_neutrinos;
12655 double massive_neutrinos;
12657 double Omega_radiation;
12660 double scalar_pivot;
12671 std::string method_Pk;
12673 std::string output_root;
12679 std::string file_par;
class CAMB that wrap CAMB routines
Class functions used by Numerical methods inside the class Cosmology.
double Pk_DeltaDelta(const double kk, const std::shared_ptr< cbl::glob::FuncGrid > Pk, const double qmin, const double qmax, const double prec=1.e-3)
the real-space matter non-linear power spectrum , computed at 1-loop
double Pk_1loop(const double kk, const std::shared_ptr< cbl::glob::FuncGrid > PkLin, const int corrtype, const double qmin, const double qmax, const double prec=1.e-3)
the one-loop power spectrum
void Table_PkCodes(const std::string code, const bool NL, std::vector< double > &lgkk, std::vector< double > &lgPk, const double redshift, const bool store_output=true, const std::string output_root="test", const double k_max=100., const std::string file_par=par::defaultString) const
write or read the table where the dark matter power spectrum is stored
std::vector< double > Pk_matter_DeWiggled(const std::string linear_method, const std::string nowiggles_method, const std::vector< double > kk, const double redshift, const double sigma_NL, const int order=4, const int nknots=10, const double lambda=0.25, const bool store_output=true, const std::string output_root="test", const bool norm=1, const double prec=1.e-4)
the dark matter power spectrum, de-wiggled (see e.g. Anderson et al 2014)
void integrals_zeta_BarrigaGatzanaga(std::vector< double > &xi_matter, std::vector< double > &Phi, const std::vector< double > rr, const std::vector< double > kk, const std::vector< double > Pk_matter) const
integrals used to compute the Barriga & Gatzanaga al. 2002 three-point correlation function model
void I_ELL_ell(std::vector< std::vector< double >> &II, const std::vector< double > rr, const int ll, const int LL, const std::vector< double > kk, const std::vector< double > Pk)
compute the quantity
double rho_crit(const double redshift, const bool unit1=false) const
the critical cosmic density
double EE_inv3(const double aa) const
inverse of the auxiliary function used to compute the Hubble function integrand of the cosmic time
double zeta_ell_1_factor(const double b1, const double beta)
the multiplicative factor for , with local bias
double DN(const double redshift, const double redshift_norm=0., const double prec=1.e-4) const
the normalised amplitude of the growing mode at a given redshift,
double m_hh
: the Hubble parameter,
bool m_unit
false quantities are provided in phyical units; true quantities are provided in cosmological units ...
std::vector< double > linear_point(const double redshift, const double rmin=60., const double rmax=150., const int nbinr=100, const std::string interpType="Spline")
the linear point
double MF_correction(const double mass, const double redshift, const std::string method_Pk, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
correction to the halo mass in non-Gaussian cosmologies
double m_Omega_k
: the density of curvature energy
void integrals_zeta_Slepian(std::vector< double > &xi_matter, std::vector< double > &xi_matter_m1, std::vector< double > &xi_matter_p1, std::vector< double > &xi_matter_2, const std::vector< double > rr, const std::vector< double > kk, const std::vector< double > Pk_matter) const
integrals used to compute the Slepian et al. 2015 three-point correlation function model
void set_Omega_all(const double OmegaB, const double OmegaCDM, const double OmegaNu, const double OmegaR, const double OmegaDE)
set the value of , , , , , and consequently the values of
double F2(const double k, const double q, const double kq)
function used to estimate the non-linear power spectrum
void set_parameters(const std::vector< CosmologicalParameter > parameter, const std::vector< double > value)
set the value of some cosmological paramters
std::vector< double > zeta_RSD(const double r1, const double r2, const int ntheta, const double b1, const double b2, const double bt, const double beta, const std::vector< double > rr, const std::vector< double > kk, const std::vector< double > Pk, const bool include_limits=false, const int max_ll=4, const bool use_k=false)
the
double scalar_amp() const
get the private member Cosmology::m_scalar_amp
double m_D_H
: the Hubble distance
~Cosmology()=default
default destructor
double m_func_sigma(const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true, std::function< double(double)> filter={}, const bool unit1=false) const
function to compute the not-yet-normalised mass variances and their derivatives
void set_w0(const double w0=-1.)
set the value of w0
double c_vir(const double c200, const double redshift, const std::string author="BryanNorman") const
virial halo concentration given
std::vector< double > Pk_matter_NoWiggles(const std::string method, const std::vector< double > kk, const double redshift, const std::string linear_method="CAMB", const int order=4, const int nknots=10, const double lambda=0.25, const bool store_output=true, const std::string output_root="test", const bool norm=1, const double prec=1.e-4)
the dark matter power spectrum without BAO wiggles
void set_Omega_radiation(const double Omega_radiation)
set the private member Cosmology::m_Omega_radiation
double neutrino_mass() const
the total neutrino mass
double square_velocity_dispersion(const double rr, const double k_int_min, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString)
the mean square velocity dispersion
double Omega_radiation() const
get the private member Cosmology::m_Omega_radiation
void generate_bias_eff_grid_one_cosmopar(std::vector< double > ¶meter, std::vector< double > &bias_eff, const std::string dir_output, const std::string file_bias_eff_grid, const cbl::cosmology::CosmologicalParameter cosmoPar, const double min_par, const double max_par, const int nbin_par, const std::vector< double > mass, const std::vector< double > mass_grid, const std::vector< double > redshift, const std::string model_bias, const std::string method_SS, const std::string meanType="mean_bias", const bool store_output=true, const std::string output_root="test", const double Delta_crit=200., const double kk=-1., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true, const cbl::cosmology::Cosmology cosmology_mass={}, const std::vector< double > redshift_source={})
compute the effective bias of dark matter haloes, by averaging the bias of a set of haloes,...
double m_Omega_CDM
: the cold dark matter density at z=0
double square_bulk_flow_Table(const double rr, const double k_int_min, const std::vector< double > lgkk, const std::vector< double > lgPk, const double redshift) const
the mean square bulk flow
double zeta_multipoles_covariance(const double Volume, const double nObjects, const int l, const int l_prime, const double r1, const double r2, const double r1_prime, const double r2_prime, const double deltaR, const std::vector< double > kk, const std::vector< double > Pk, const std::vector< double > rr, const std::vector< double > Xi, const double prec=1.e-3)
the dark matter three-point correlation function multipoles covariance model, by Slepian et al....
double Omega_k() const
get the private member Cosmology::m_Omega_k
double m_scalar_pivot
the scalar pivot k in
void get_barred_xi(std::vector< double > rr, std::vector< double > Xi, std::vector< double > &Xi_, std::vector< double > &Xi__, const std::string method_Pk, const double redshift, const bool xiType=0, const double k_star=-1., const bool xiNL=0, const int norm=-1, const double r_min=0.1, const double r_max=150., const double k_min=0.001, const double k_max=100., const double aa=0., const double prec=1.e-2, const std::string file_par=par::defaultString) const
get the barred dark matter correlation functions
double wtheta_DM(const double theta, const std::vector< double > zz, const std::vector< double > phiz, const std::string interpolationMethod, const CoordinateUnits coordUnits=CoordinateUnits::_degrees_, const bool GSL=false, const std::string method_Pk="CAMB", const bool NL=false, const bool store_output=true, const std::string output_root="test", const int norm=-1, const double k_min=1.e-4, const double k_max=100, const double prec=1.e-2, const std::string file_par=par::defaultString)
the dark matter angular two-point correlation function
double rs_EH(const double T_CMB=par::TCMB) const
the sound horizon at the drag epoch predicted by Eisentein & Hu 1998
double m_fNL
: the non-Gaussian amplitude
double t_H() const
get the private member Cosmology::m_t_H
double unevolved_mass_function(const double mass_accr) const
the unevolved mass function
double xis_gBAO(const double rp, const double pi, const double f_sigma8, const double bias_sigma8, const double redshift, std::vector< double > rr, std::vector< double > Xi, std::vector< double > &Xi_, std::vector< double > &Xi__, const bool store_output=true, const std::string output_root="test", const double k_star=-1., const double x_min=-3000., const double x_max=3000., const int step_x=500)
the function ξg,BAO(s) of the Chuang & Wang 2012 model
double rho_m(const double redshift=0., const bool unit1=false, const bool nu=false) const
the mean cosmic background density
double Mag_Volume_limited(const double z_max=1., const double mag_lim=-20.) const
maximum absolute magnitude to have a volume-limited catalogue
double xi2D_dispersionModel(const double rp, const double pi, const double f_sigma8, const double bias_sigma8, const double sigmav, const std::string method_Pk, const double redshift, const int FV, const bool NL, std::vector< double > rr, std::vector< double > &Xi, std::vector< double > &Xi_, std::vector< double > &Xi__, const bool store_output=true, const std::string output_root="test", const int index=-1, const bool bias_nl=0, const double bA=-1., const bool xiType=0, const double k_star=-1., const bool xiNL=0, const double v_min=-3000., const double v_max=3000., const int step_v=500, const int norm=-1, const double r_min=0.1, const double r_max=150., const double k_min=0.001, const double k_max=100., const double aa=0., const bool GSL=false, const double prec=1.e-2, const std::string file_par=par::defaultString)
2D correlation function, ξ(rp,π), predicted by the dispersion model
double dskewnessdM(const double mass, const std::string method_Pk, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
the derivative of the skewness, ds/dM
void get_xi(std::vector< double > &rr, std::vector< double > &Xi, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const bool xiType=0, const double k_star=-1., const bool xiNL=0, const int norm=-1, const double r_min=0.1, const double r_max=150., const double k_min=0.001, const double k_max=100., const double aa=0., const bool GSL=false, const double prec=1.e-2, const std::string file_par=par::defaultString)
get the dark matter two-point correlation function
double m_wa
: the parameter of the dark energy equation of state (CPL parameterisation)
double deltav_NL(const double deltav=-2.71) const
Non-Linear (under)density contrast.
double zeta_ell_k_factor(const double b1, const double beta)
the multiplicative factor for , with local bias
double Omega(const double redshift=0.) const
the cosmic density at a given redshift
double bispectrum(const std::vector< double > kk, const std::string method_Pk, const bool store_output=true, const std::string output_root="test", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString)
the bispectrum
double z_drag() const
redshift of drag epoch
double m_RhoZero
: the mean density of the Universe at z=0 [Msun*Mpc^-3]
double square_bulk_flow(const double rr, const double k_int_min, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString)
the mean square bulk flow
double m_Omega_radiation
: the radiation density at z=0
double xi0_Kaiser(const double rad, const double f_sigma8, const double bias_sigma8, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const bool xiType=0, const double k_star=-1., const bool NL=false, const int norm=-1, const double r_min=0.1, const double r_max=150., const double k_min=0.001, const double k_max=100., const double aa=0., const bool GSL=false, const double prec=1.e-2, const std::string file_par=par::defaultString)
monopole of the redshift-space two-point correlation function in the Kaiser limit
std::vector< double > redshift_distribution_haloes_selection_function(const std::vector< double > redshift, const double Area_degrees, const double Mass_min, const double Mass_max, const std::string model_MF, const std::string method_SS, const std::string selection_function_file, const std::vector< int > column={}, const bool store_output=true, const std::string output_root="test", const double Delta=200, const bool isDelta_critical=false, const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true)
redshift distribution of dark matter haloes, given a selection function
void set_model(const std::string model="LCDM")
set the cosmologial model used to compute distances
double m_Pk0_CLASS
the normalisation of the power spectrum for CLASS [http://class-code.net/]
double n_spec() const
get the private member Cosmology::m_n_spec
double converted_mass(const double mass, const cosmology::Cosmology cosmology, const double redshift, const double redshift_source=-1.) const
convert a cluster mass estimated in a different cosmology
std::vector< double > bias_eff_selection_function(const glob::FuncGrid interp_sigma, const glob::FuncGrid interp_DnSigma, const glob::FuncGrid interp_SF, const double Mass_min, const double Mass_max, const std::vector< double > redshift, const std::string model_bias, const std::string model_MF, const std::string method_SS, const double alpha=1., const bool store_output=true, const std::string output_root="test", const double Delta_crit=200., const double kk=-1., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
effective bias of dark matter haloes, computed using a given selection function; σ(mass) and dlnσ/dM ...
void set_type_NG(const int type_NG=1)
set the value of the non-Gaussian shape
double Delta_c(const double redshift, const std::string author="BryanNorman") const
the critical overdensity
double m_sigma2M_notNormalised(const double mass, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true, const bool unit1=false) const
the not-yet-normalised mass variance,
double G2(const double k, const double q, const double kq)
function used to estimate the non-linear power spectrum
double m_n_spec
: the primordial spectral index
double zeta_ell_precyclic(const double r1, const double r2, const int ell, const double b1, const double b2, const double bt, const double beta, std::vector< std::shared_ptr< glob::FuncGrid >> interp_xi_ell, const bool use_k, std::shared_ptr< glob::FuncGrid2D > interp_k_ell)
the pre-cyclic
void generate_bias_eff_grid_two_cosmopars(std::vector< double > ¶meter1, std::vector< double > ¶meter2, std::vector< std::vector< double >> &bias_eff, const std::string dir_output, const std::string file_bias_eff_grid, const cbl::cosmology::CosmologicalParameter cosmoPar1, const double min_par1, const double max_par1, const int nbin_par1, const cbl::cosmology::CosmologicalParameter cosmoPar2, const double min_par2, const double max_par2, const int nbin_par2, const std::vector< double > mass, const std::vector< double > mass_grid, const std::vector< double > redshift, const std::string model_bias, const std::string method_SS, const std::string meanType="mean_bias", const bool store_output=true, const std::string output_root="test", const double Delta=200., const double kk=-1., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true, const cbl::cosmology::Cosmology cosmology_mass={}, const std::vector< double > redshift_source={})
effective bias of dark matter haloes, computed by averaging the bias of a set of haloes,...
std::vector< double > Q_DM_eq(const std::vector< double > rr, const std::string model, const std::vector< double > kk, const std::vector< double > Pk_matter) const
the dark matter equilateral reduced three-point correlation function
std::vector< std::vector< double > > XiMonopole_covariance(const int nbins, const double rMin, const double rMax, const double nn, const double Volume, const std::vector< double > kk, const std::vector< double > Pk0, const int IntegrationMethod=1)
the covariance matrix of the first three non-null multipoles of the two-point correlation function
double sigma_v(const double redshift=0., const std::string method_Pk="CAMB", const bool store_output=true, const std::string output_root="test", const int norm=-1, const double k_min=0.001, const double k_max=100., const int bin_k=512, const double prec=1.e-2, const std::string file_par=par::defaultString, const bool unit1=false)
the linear-order one-dimensional pairwise velocity dispersion,
std::vector< std::vector< double > > Pk_TNS_AB_1loop(std::vector< double > kk, const double mu, const std::string method, const double redshift, const bool store_output, const std::string output_root, const int norm, const double k_min, const double k_max, const double prec)
the A and B correction terms for the TNS model at 1-loop from the multipole expansion
void set_sigma8(const double sigma8=-1.)
set the value of σ8
double z_acc() const
redshift at which the Universe begins to accelerate
double potential_spectral_amplitude(const std::string method_Pk, const bool store_output=true, const std::string output_root="test", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString)
the potential spectral amplitude
double gg(const double redshift=0.) const
the linear growth factor at a given redshift,
double MhaloMin(const int n_halo, const double Area, const bool angle_rad, const double z_min, const double z_max, const double Mmax, const double lgM1_guess, const double lgM2_guess, const std::string model_MF, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const double Delta=200, const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true) const
minimum halo mass, given the number of haloes in a given region of sky
double RhoZero() const
get the private member Cosmology::m_RhoZero
std::vector< std::vector< double > > Pk_TNS_dd_dt_tt(std::vector< double > kk, const std::string method, const double redshift, const bool store_output, const std::string output_root, const int norm, const double k_min=0.001, const double k_max=100., const double prec=1.e-2)
the non-linear , , matter power spectra
double Omega_baryon() const
get the private member Cosmology::m_Omega_baryon
double scalar_pivot() const
get the private member Cosmology::m_scalar_pivot
void k_ell(std::vector< std::vector< double >> &KK, const std::vector< double > rr, const int ll, const std::vector< double > kk, const std::vector< double > Pk)
compute the quantity
double D_C_LCDM(const double redshift) const
the comoving line-of-sight distance at a given redshift
std::vector< double > C_l_DM(const int lmax, const std::vector< double > zz, const std::vector< double > phiz, const std::string interpolationMethod, const std::string method_Pk="CAMB", const bool store_output=true, const std::string output_root="test", const int norm=-1, const double k_min=1.e-4, const double k_max=100, const double prec=1.e-2, const std::string file_par=par::defaultString)
the dark matter angular linear power spectrum .
std::vector< double > zeta_halo(const double r1, const double r2, const std::vector< double > theta, const double b1, const double b2, const std::string model, const std::vector< double > kk, const std::vector< double > Pk_matter) const
the local-bias model of the three-point correlation function of dark matter haloes
double wp_DM(const double rp, const std::string method_Pk, const bool NL, const double redshift, const double pimax, const bool store_output=true, const std::string output_root="test", const int norm=-1, const double r_min=1.e-3, const double r_max=350., const double k_min=0.001, const double k_max=100., const double aa=0., const bool GSL=false, const double prec=1.e-2, const std::string file_par=cbl::par::defaultString)
the dark matter projected correlation function
double rs_integrand(const double redshift, const double T_CMB=2.7255) const
the sound horizon integrand
double m_Omega_baryon
: the baryon density at z=0
void set_tau(const double tau=0.09)
set the value of the τ
double w0() const
get the private member Cosmology::m_w0
double tau() const
get the private member Cosmology::m_tau
void set_parameter(const CosmologicalParameter parameter, const double value)
set the value of one cosmological paramter
std::string create_grid_sigmaM(const std::string method_SS, const double redshift, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true) const
auxiliary function to create a grid file with σ(M)
double H0() const
get the private member Cosmology::m_H0
double r_rL(const double deltav=-2.71) const
expansion factor
double w_CPL(const double redshift=0.) const
the DE equation of state in the CPL parameterisation, as a function of redshift
double m_elf_dz(const double phi) const
the incomplete elliptic integral
double size_function(const double RV, const double redshift, const std::string model, const double b_eff, double slope=0.854, double offset=0.420, const double deltav_NL=-0.795, const double del_c=1.69, const std::string method_Pk="EisensteinHu", const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true) const
the void size function
double sigma8() const
get the private member Cosmology::m_sigma8
double Pk0_EH() const
get the private member Cosmology::m_Pk0_EH
double massless_neutrinos() const
get the private member Cosmology::m_massless_neutrinos
bool unit() const
get the private member Cosmology::m_unit
double zeta_DM_BarrigaGatzanaga(const double r1, const double r2, const double theta, std::vector< double > &rr, std::vector< double > &xi_matter, std::vector< double > &Phi, const std::vector< double > kk, const std::vector< double > Pk_matter) const
the dark matter three-point correlation function model by Barriga & Gatzanaga et al....
double rs_CAMB() const
the sound horizon at the drag epoch estimated with CAMB [http://camb.info/], analytical formula by Au...
double fsigma8(const double redshift, const std::string method_Pk, const bool store_output=true, const std::string output_root="test", const bool NL=false, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString) const
f*σ8: the linear growth rate times the dark matter rms mass fluctuation within 8 Mpc/h
int m_type_NG
the non-Gaussian shape (type=1 local, type=2 equilateral, type=3 enfolded, type=4 orthogonal)
void set_OmegaNu(const double Omega_neutrinos=0., const double massless_neutrinos=3.04, const int massive_neutrinos=0)
set the value of Ων
double sound_speed(const double redshift, const double T_CMB=2.7255) const
the sound speed
std::vector< double > Q_halo(const double r1, const double r2, const std::vector< double > theta, const double b1, const double b2, const std::string model, const std::vector< double > kk, const std::vector< double > Pk_matter) const
the local-bias model of the reduced three-point correlation function of dark matter haloes
double m_Pk0_MPTbreeze
the normalisation of the power spectrum for MPTbreeze [http://arxiv.org/abs/1207.1465]
std::vector< double > bias_eff_mass_grid(const std::vector< double > MM, const std::vector< double > redshift, const std::string model_bias, const std::string method_SS, const std::string meanType="mean_bias", const bool store_output=true, const std::string output_root="test", const double Delta_crit=200., const double kk=-1., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
effective bias of dark matter haloes, computed by averaging the bias of a set of haloes,...
double max_redshift(const double Volume, const double Area, const double z_min) const
maximum redshift for a given volume, sky area and minimum redshift
double quadrupole(const double Mass_min, const double Mass_max, const double redshift, const std::string model_bias, const std::string model_MF, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const double Delta=200., const double kk=-1., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
the normalised quadrupole Q
void remove_output_Pk_tables(const std::string code, const bool NL, const double redshift, const std::string output_root="test") const
remove the output generated by the methods CAMB, MPTbreeze or CLASS
void set_OmegaDE(const double Omega_DE=0.73)
set the value of ΩDE
double Volume(const double z1, const double z2, const double Area) const
comoving volume for a given redshift range and sky area
double EE(const double redshift=0.) const
auxiliary function used to compute the Hubble function
double zeta_ell_3_factor(const double b1, const double beta)
the multiplicative factor for , with local bias
double n_haloes(const double Mass_min, const double Mass_max, const double z_min, const double z_max, const bool angle_rad, const std::string model_MF, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const double Delta=200, const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true)
number of dark matter haloes per steradian or square degree, for a given redshift range
double m_bias_halo_generator(const double Sigma, const double redshift, const std::string author, const double Delta=200.) const
auxiliary function to compute the halo bias
void medianwf(const double ff, const std::string model_model, std::vector< double > &wf) const
median formation w
double m_Pk0_EH
the normalisation of the power spectrum for Eisenstein & Hu [http://background.uchicago....
double zeta_DM_Slepian(const double r1, const double r2, const double theta, std::vector< double > &rr, std::vector< double > &xi_matter, std::vector< double > &xi_matter_m1, std::vector< double > &xi_matter_p1, std::vector< double > &xi_matter_2, const std::vector< double > kk, const std::vector< double > Pk_matter, const int norders=9, const double prec=1.e-3) const
the dark matter three-point correlation function model by Slepian et al. 2015
double F_AP(const double redshift) const
F_AP, the ALCOCK-PACZYNSKI distortion parameter.
void Table_XiCodes(const std::string code, const bool NL, std::vector< double > &rr, std::vector< double > &xi, const double redshift, const bool store_output, const std::string output_root, const double k_max, std::string file_par) const
write or read the table where the dark matter two-point correlation function is stored
std::vector< double > zeta_expansion_Slepian(const double r1, const double r2, const double b1, const double b2, std::vector< double > &rr, std::vector< double > &xi_matter, std::vector< double > &xi_matter_m1, std::vector< double > &xi_matter_p1, std::vector< double > &xi_matter_2, const int norders=9, const double prec=1.e-3) const
the terms of the expansion
double Q_nonLocal(const double r1, const double r2, const double theta, std::vector< double > &rr, std::vector< double > &xi_matter, std::vector< double > &Phi, const std::vector< double > kk, const std::vector< double > Pk_matter) const
the non-local contribution to the reduced dark matter three-point correlation function
double Pk0_CAMB() const
get the private member Cosmology::m_Pk0_CAMB
double error_beta(const double redshift, const double bias, const double err_bias) const
the error on the specific growth rate β
double D_V(const double redshift) const
the average distance at a given redshift, used to rescale the correlation function
double bias_correction(const double kk, const double mass, const std::string method_Pk, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
correction to the halo bias in non-Gaussian cosmologies
int type_NG() const
get the private member Cosmology::m_type_NG
double z_eq_rad(const double T_CMB=2.7255) const
redshift of matter-radiation equality
double Pk_ThetaTheta(const double kk, const std::shared_ptr< cbl::glob::FuncGrid > Pk, const double qmin, const double qmax, const double prec=1.e-3)
the real-space matter non-linear power spectrum , computed at 1-loop
double m_sigma8
: the power spectrum normalisation
std::vector< std::vector< double > > Pk_eTNS_terms_1loop(std::vector< double > kk, const std::string method, const double redshift, const bool store_output, const std::string output_root, const int norm, const double k_min=0.001, const double k_max=100., const double prec=1.e-2)
The expanded correction terms for the extended TNS model (eTNS)
double qq(const double redshift=0.) const
the deceleration parameter at a given redshift
double m_H0
: the Hubble constant at z=0 [km/sec/Mpc]
double m_rs
the sound horizon
double m_massless_neutrinos
: the effective number (for QED + non-instantaneous decoupling)
double m_sigma2R_notNormalised(const double radius, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true, const bool unit1=false) const
the not-yet-normalised mass variance,
double deltac(const double redshift) const
spherical collapse density threshold at a given redshift
double OmegaR(const double redshift=0.) const
the radiation density at a given redshift
double hh() const
get the private member Cosmology::m_hh
Cosmology(const double Omega_matter=0.27, const double Omega_baryon=0.046, const double Omega_neutrinos=0., const double massless_neutrinos=3.04, const int massive_neutrinos=0, const double Omega_DE=0.73, const double Omega_radiation=0., const double hh=0.7, const double scalar_amp=2.46e-9, const double scalar_pivot=0.05, const double n_spec=0.96, const double w0=-1., const double wa=0., const double fNL=0., const int type_NG=1, const double tau=0.09, const std::string model="LCDM", const bool unit=true)
constructor
std::vector< std::vector< double > > Pk_TNS_AB_terms_1loop(std::vector< double > kk, const std::string method, const double redshift, const bool store_output, const std::string output_root, const int norm, const double k_min=0.001, const double k_max=100., const double prec=1.e-2)
the expanded A and B correction terms for the TNS model
double pw(const double ww, const double ff, const std::string author) const
differential distribution
double Omega_neutrinos() const
get the private member Cosmology::m_Omega_neutrinos
void set_scalar_pivot(const double scalar_pivot=0.05)
set the value of the scalar pivot
double dV_dZdOmega(const double redshift, const bool angle_rad) const
the derivative of the comoving volume, d2V/(dz*dΩ) at a given redshift
void set_unit(const bool unit=true)
set the value of unit
double deltav_L(const double deltav_NL, const double b_eff, double slope=0.854, double offset=0.420) const
Linear (under)density contrast.
double wa() const
get the private member Cosmology::m_wa
double frk(const double kk, const double mass, const std::string method_Pk, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
auxiliary function to estimate cosmological quantities in non-Gaussian cosmologies
double zeta_ell_0_factor_tidal(const double gamma_t, const double beta)
the multiplicative factor for , with non-local bias
void set_fNL(const double fNL=0.)
set the value of fNL
double cumPw(const double ww, const double ff, const std::string author) const
cumulative distribution
double m_asn_dz(const double ss) const
the inverse sine amplitude of the Jacobian elliptic function
double ys(const double redshift, const std::string method_Pk, const double T_CMB=par::TCMB) const
the fiducial cosmology independent ratio rs/DV, valid choices for method_Pk are: EisensteinHu [http:/...
double skewness(const double mass, const std::string method_Pk, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
the skewness
double zeta_ell_2_factor(const double b1, const double gamma, const double beta)
the multiplicative factor for , with local bias
double dnsigma2M(const int nd, const double mass, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true, const bool unit1=false) const
the first derivative of the mass variance,
double zeta_precyclic_Slepian(const double r1, const double r2, const double mu, const double b1, const double b2, const glob::FuncGrid interp_xi_matter, const glob::FuncGrid interp_xi_matter_m1, const glob::FuncGrid interp_xi_matter_p1, const glob::FuncGrid interp_xi_matter_2) const
the pre-cyclic three-point correlation function as described in Slepian et al. 2015
double Pk0_CLASS() const
get the private member Cosmology::m_Pk0_CLASS
double m_Omega_neutrinos
: the density of massive neutrinos at z=0
void xi_r_n_pm(std::vector< double > &xi_n_p, std::vector< double > &xi_n_m, const std::vector< double > rr, const int nn, const std::vector< double > kk, const std::vector< double > Pk)
compute the power spectrum integral transform
double m_t_H
: the Hubble time
double mrk(const double kk, const double mass, const std::string method_Pk, const bool store_output=true, const std::string output_root="test", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString)
auxiliary function to estimate cosmological quantities in non-Gaussian cosmologies
double HH(const double redshift=0.) const
the Hubble function
double m_Pk0_CAMB
the normalisation of the power spectrum for CAMB [http://camb.info/]
double rs() const
get the sound horizon at recombination
double As(const double sigma8) const
amplitude of the curvature perturbations
double EE2(const double redshift=0.) const
auxiliary function used to compute the deceleration parameter
double m_acn_dz(const double cc) const
the inverse cosine amplitude of the Jacobian elliptic function
double CMN(const double rr, const double k_int_min, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const double k_max=100., const std::string file_par=par::defaultString) const
the Cosmic Mach Number
double m_MF_generator(const double Mass, const double Sigma, const double Dln_Sigma, const double redshift, const std::string model_MF, const double Delta=200., const bool default_delta=true, const double delta_t=1.686)
auxiliary function to compute the mass function
double pz(const double m0, const double z0, const double frac, const double redshift, const std::string model_model, const std::string method_SS, const bool store_output=true, const std::string output_root="test") const
formation probability
double m_w0
: the parameter of the dark energy equation of state (CPL parameterisation)
double beta(const double redshift, const double bias) const
the specific growth rate β
double f_k(const double k, const std::shared_ptr< cbl::glob::FuncGrid > PkLin, const double qmin, const double qmax, const double prec=1.e-3)
function used to estimate the non-linear power spectrum
double D_A(const double redshift) const
the angular diameter distance at a given redshift
double Redshift(const double d_c=1., const double z1_guess=0., const double z2_guess=10., const double prec=0.0001) const
redshift at a given comoving distance
std::string model() const
get the private member Cosmology::m_model
void print_parameters() const
print the values of the private members on the screen
double EE_inv2(const double redshift=0.) const
inverse of the auxiliary function used to compute the Hubble function, integrand of the lookback time
double m_serf_dz(const double yy) const
the inverse truncated series necessary to compute sn-1(s|m) in ASN_DZ
double D_H() const
get the private member Cosmology::m_D_H
double bias_halo(const double Mass, const double redshift, const std::string author, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double Delta=200., const double kk=-1., const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
bias of dark matter haloes
std::vector< double > Pk_matter_NoWiggles_gaussian(const std::vector< double > kk, const std::vector< double > PkLin, const std::vector< double > PkApprox, const double lambda, const std::string method)
the dark matter power spectrum without BAO wiggles.
void medianzf(const double ff, const double mass, const double z0, const std::string model_model, const std::string method_SS, std::vector< double > &zf, const bool store_output=true, const std::string output_root="test") const
median formation z
std::vector< double > Pk_matter(const std::vector< double > kk, const std::string method_Pk, const bool NL, const double redshift, const bool store_output=true, const std::string output_root="test", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString, const bool unit1=false)
the dark matter power spectrum
std::vector< double > AP_corr(const cbl::cosmology::Cosmology cosm_true, const std::vector< double > redshift)
Supplementary function to compute a correction factor to apply to the void size function,...
void set_scalar_amp(const double scalar_amp=2.46e-9)
set the value of As
double xi_matter(const double rr, const std::string method_Pk, const bool NL, const double redshift, const bool store_output=true, const std::string output_root="test", const int norm=-1, const double k_min=0.001, const double k_max=100., const double aa=0., const bool GSL=false, const double prec=1.e-2, const std::string file_par=par::defaultString)
the dark matter two-point correlation function
double xi_matter_DeWiggle(const double rr, const double redshift, const double sigma_NL, const bool store_output=true, const std::string output_root="test", const bool norm=1, const double k_min=0.001, const double k_max=100., const double aa=1., const double prec=1.e-2)
the dark matter two-point correlation function, de-wiggled (see e.g. Anderson et al 2014)
double error_beta_measured(const double Volume, const double density, const double Mass_min, const double Mass_max, const double redshift, const std::string model_bias, const std::string model_MF, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const double Delta=200., const double kk=-1., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
the error on the specific growth rate β from Bianchi et al. 2012
double f_DE(const double redshift=0.) const
auxiliary function used to compute the Hubble function
void set_OmegaB_h2(const double Omega_baryonh2=0.0222)
set the value of Ωb, keeping ΩCDM=ΩM-Ωb
void m_Table_Pk_parameterFile(const std::string code, const std::string file_par, const bool NL, std::vector< double > &lgkk, std::vector< double > &lgPk, const double redshift, const std::string output_root="test") const
write and read the table where the dark matter power spectrum is stored; it is used when a parameter ...
double cosmic_time(const double redshift=0.) const
cosmic time at a given redshift
double Omega_matter() const
get the private member Cosmology::m_Omega_matter
double m_Omega_DE
: the dark energy density at z=0
double lookback_time(const double redshift=0.) const
lookback time at a given redshift
double Redshift_LCDM(const double d_c=1., const double z1_guess=0., const double z2_guess=10., const bool go_fast=1, const double prec=0.0001) const
redshift at a given comoving distance
double zeta_ell_4_factor_tidal(const double gamma_t, const double beta)
the multiplicative factor for , with non-local bias
double mass_function_fast(const double Mass, const double redshift, const std::string model_MF, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const double Delta=200., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
the mass function of dark matter haloes (filaments and sheets) computed quickly using a grid
double Pk_ThetaTheta_fitting_function(const double kk, const std::string method_Pk, const double redshift, const std::string author, const bool store_output, const std::string output_root, const bool NL, const int norm, double k_min, double k_max, const double prec, const std::string file_par, const bool unit1)
the dark matter velocity divergence power spectrum
void set_default()
internal function to set default values
double Redshift_time(const double time, const double z1_guess, const double z2_guess) const
redshift at a given cosmic time
double EE_inv(const double redshift=0.) const
inverse of the auxiliary function used to compute the Hubble function integrand of the comoving dista...
double zeta_ell_2_factor_tidal(const double gamma_t, const double beta)
the multiplicative factor for , with non-local bias
void set_H0(const double H0=70., const bool warn=true)
set the value of H0
std::string m_model
the cosmologial model used to compute distances
double dnsigma2R(const int nd, const double radius, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true, const bool unit1=false) const
the nth-order derivative of the mass variance,
void Pk_0(const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString)
normalisation of the power spectrum
std::vector< double > mass_function_selection_function_vector(const std::vector< double > mass, const double z_min, const double z_max, const std::string model_MF, const std::string method_SS, const std::string selection_function_file, const std::vector< int > column={}, const bool store_output=true, const std::string output_root="test", const double Delta=200, const bool isDelta_critical=false, const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true)
mass function given a selection function
double zeta_single_BarrigaGatzanaga(const double r1, const double r2, const double theta, const std::vector< double > xi, const std::vector< double > dxi, const std::vector< double > dPhi) const
the single term of the dark matter three-point correlation function model by Barriga & Gatzanaga et a...
std::vector< std::vector< double > > XiMultipoles(const int nbins, const double rMin, const double rMax, const std::vector< double > kk, const std::vector< double > Pk0, const std::vector< double > Pk2, const std::vector< double > Pk4, const int IntegrationMethod=1)
the first three non-null multipoles of the two-point correlation function
double Lum_bol(const double redshift=0., const double flux=1.) const
bolometric luminosity
double Gamma_3PCF(const double r1, const double r2, const double theta, const std::vector< double > xi, const std::vector< double > dPhi) const
function to compute non-local contribution to three-point correlation function; specifically,...
double OmegaR_zeq(const double z_eq=3395.) const
the radiation density, as a function of the redshift of radiation-matter equality
double sigmaR_DM(const double RR, const int corrType, const std::string method_Pk, const double redshift, const double pimax=40, const bool store_output=true, const std::string output_root="test", const bool NL=1, const int norm=-1, const double r_min=1.e-3, const double r_max=350., const double k_min=0.001, const double k_max=100., const double aa=0., const bool GSL=false, const double prec=1.e-2, const std::string file_par=par::defaultString)
the dark matter rms mass fluctuation
double Hdot(const double redshift=0.) const
derivative of the Hubble function at a given redshift
void integrals_Q_nonLocal(std::vector< double > &xi_matter, std::vector< double > &Phi, const std::vector< double > rr, const std::vector< double > kk, const std::vector< double > Pk_matter, const double prec) const
integral functions for the three-point correlation model
double mass_function_fR(const double Mass, const double redshift, const std::string model_MF, const double f_R0=0., const bool store_output=true, const double Delta=200., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true, const bool default_delta=true, const double delta_t=1.686)
the mass function of dark matter haloes in f(R) cosmologies (see Hu & Sawicki 2007) computed with the...
std::vector< double > Pk_matter_Linear(const std::string method, const std::vector< double > kk, const double redshift, const bool store_output=true, const std::string output_root="test", const bool norm=1, const double prec=1.e-4)
the dark matter linear power spectrum.
double Q_DM_Slepian(const double r1, const double r2, const double theta, std::vector< double > &rr, std::vector< double > &xi_matter, std::vector< double > &xi_matter_m1, std::vector< double > &xi_matter_p1, std::vector< double > &xi_matter_2, const std::vector< double > kk, const std::vector< double > Pk_matter, const int norders=9, const double prec=1.e-3) const
the dark matter reduced three-point correlation function model by Slepian et al. 2015
double m_mass_function(const double Mass, std::shared_ptr< void > mass_function_params)
auxiliary function to compute the mass function of dark matter haloes (filaments and sheets)
double m_scalar_amp
: the initial scalar amplitude of the power spectrum
std::vector< double > bias_eff_mass(const std::vector< double > MM, const std::vector< double > redshift, const std::string model_bias, const std::string method_SS, const std::string meanType="mean_bias", const bool store_output=true, const std::string output_root="test", const double Delta=200., const double kk=-1., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
effective bias of dark matter haloes, computed by averaging the bias of a set of haloes
double sigma2M(const double mass, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true, const bool unit1=false) const
the mass variance,
int massive_neutrinos() const
get the private member Cosmology::m_massive_neutrinos
double xi_star(const double rr, const double redshift, const bool store_output=true, const std::string output_root="test", const double k_star=-1., const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString)
the function ξ* of the Chuang & Wang 2012 model
double wf(const double mm, const double redshift, const double ff, const double zf, const std::string method_SS, const bool store_output=true, const std::string output_root="test") const
rescaled variable w as in Lacey and Coles 1993
double m_Omega_matter
: the density of baryons, cold dark matter and massive neutrinos (in units of the critical density) a...
double M_vir(const double r_vir, const double redshift, const std::string author="BryanNorman", const bool unit1=false) const
the virial mass, given the virial radius and the redshift
double D_C(const double redshift) const
the comoving line-of-sight distance at a given redshift
double k_star(const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const double k_max=100., const std::string file_par=par::defaultString) const
the k* parameter
void set_wa(const double wa=0.)
set the value of wa
void set_Omega(const double Omega_matter=0.27)
set the value of ΩM, keeping ΩDE=1-ΩM-Ωrad-Ωk
double Am(const std::string method_Pk, const bool store_output=true, const std::string output_root="test", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString)
the amplitude of the matter power spectrum
double Sn_PT(const int nn, const double RR, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true) const
the hierarchical moments Sn
double fNL() const
get the private member Cosmology::m_fNL
void D_C_table(const std::string file_table, const double z_min, const double z_max, const int step, std::vector< double > &Redshift, std::vector< double > &dc) const
create a table of [redshift, comoving line-of-sight distance]
void eff_l_l1(std::vector< std::vector< double >> &eff, const std::vector< double > rr, const int l, const int l1, const std::vector< double > kk, const std::vector< double > Pk)
compute the power spectrum integral transform
double bias_eff(const double Mass_min, const double Mass_max, const double redshift, const std::string model_bias, const std::string model_MF, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const double Delta=200., const double kk=-1., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true)
the effective bias of dark matter haloes, with masses in a given range and at a given mean redshift
void set_n_spec(const double n_spec)
set the value of nspec
double Pk_DeltaTheta_fitting_function(const double kk, const std::string method_Pk, const double redshift, const std::string author, const bool store_output, const std::string output_root, const bool NL, const int norm, double k_min, double k_max, const double prec, const std::string file_par, const bool unit1)
the dark matter cross power spectrum
double mass_function(const double Mass, const double redshift, const std::string model_MF, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const double Delta=200., const std::string interpType="Linear", const int norm=-1, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string input_file=par::defaultString, const bool is_parameter_file=true, const bool default_delta=true, const double delta_t=1.686)
the mass function of dark matter haloes (filaments and sheets)
double denominator_Q(const double r1, const double r2, const double theta, const std::vector< double > rr, const std::vector< double > xi_matter) const
the normalization factor for reduced three-point correlation function
std::vector< std::vector< double > > zeta_covariance(const double Volume, const double nObjects, const std::vector< double > theta, const double r1, const double r2, const double deltaR, const std::vector< double > kk, const std::vector< double > Pk, const int norders=10, const double prec=1.e-3, const bool method=false, const int nExtractions=10000, const std::vector< double > mean={}, const int seed=543)
the dark matter three-point correlation function covariance model
double mean_redshift_haloes_selection_function(const double z_min, const double z_max, const double Mass_min, const double Mass_max, const std::string model_MF, const std::string method_SS, const std::string selection_function_file, const std::vector< int > column={}, const bool store_output=true, const std::string output_root="test", const double Delta=200, const bool isDelta_critical=false, const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true)
the mean redshift of a dark matter haloe sample, given a selection function
double zeta_ell_0_factor(const double b1, const double gamma, const double beta)
the multiplicative factor for , with local bias
std::string Pk_output_file(const std::string code, const bool NL, const double redshift, const bool run=0, const bool store_output=true, const std::string output_root="test", const double k_max=100., const std::string file_par=par::defaultString)
return the path to the power spectrum output
double Pk0_MPTbreeze() const
get the private member Cosmology::m_Pk0_MPTbreeze
void set_RhoZero(const double RhoZero=7.5e10)
set the value of ρ0
double D_L(const double redshift) const
the luminosity distance at a given redshift
void set_rs(const double rs=-1)
set the value of the ;
double sigma2R(const double radius, const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true, const bool unit1=false) const
the mass variance,
double OmegaK(const double redshift=0.) const
the density of curvature energy at a given redshift
double Q_DM_BarrigaGatzanaga(const double r1, const double r2, const double theta, std::vector< double > &rr, std::vector< double > &xi_matter, std::vector< double > &Phi, const std::vector< double > kk, const std::vector< double > Pk_matter) const
the dark matter reduced three-point correlation function model by Barriga & Gatzanaga et al....
double linear_growth_rate(const double redshift, const double prec=1.e-4) const
the linear growth rate at a given redshift,
double xisnl_gnw(const double rp, const double pi, const double f_sigma8, const double bias_sigma8, const double bA, const double redshift, std::vector< double > rr, std::vector< double > Xi, std::vector< double > &Xi_, std::vector< double > &Xi__, const bool store_output=true, const std::string output_root="test")
the function ξg,nw(s) of the Chuang & Wang 2012 model
void set_OmegaM(const double Omega_matter=0.27)
set the value of ΩM
double Sigman_PT(const int nn, const double RR, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true) const
the deprojected hierarchical moments Σn
double Omega_CDM() const
get the private member Cosmology::m_Omega_CDM
double Az(const double redshift) const
the acoustic parameter
void m_Table_Pk_CLASS(const bool NL, std::vector< double > &lgkk, std::vector< double > &lgPk, const double redshift, const bool store_output=true, const std::string output_root="test", const double k_max=100.) const
write and read the table where the dark matter power spectrum computed with CLASS is stored
double f_nu(const double SS, const double del_v, const double del_c) const
(approximation)
double r_vir(const double M_vir, const double redshift, const std::string author="BryanNorman", const bool unit1=false) const
the virial radius, given the virial mass and the redshift
void run_CAMB(const bool NL, const double redshift, const std::string output_root=par::defaultString, const std::string output_dir=par::defaultString, const double k_max=100., const std::string file_par=par::defaultString) const
run CAMB [http://camb.info/]
std::vector< std::vector< double > > Pk_TNS_AB_multipoles(std::vector< double > kk, const std::string method, const double redshift, const bool store_output, const std::string output_root, const int norm, const double k_min, const double k_max, const double prec)
the multipoles of the A and B correction terms for the TNS model
std::vector< double > Q_DM(const double r1, const double r2, const std::vector< double > theta, const std::string model, const std::vector< double > kk, const std::vector< double > Pk_matter) const
the dark matter reduced three-point correlation function
double xi2D_CW(const double rp, const double pi, const double beta, const double bias_lin, const double bA, const double sigmav0, const double cmu, const double cs1, const double cs2, const double redshift, std::vector< double > rr1, std::vector< double > Xi1, std::vector< double > rr2, std::vector< double > Xi2, std::vector< double > &Xi1_, std::vector< double > &Xi1__, std::vector< double > &Xi2_, std::vector< double > &Xi2__, const bool store_output=true, const std::string output_root="test", const bool BAO=1, const bool xiType=0, const double k_star=-1, const bool xiNL=0, const double r_min=0.1, const double r_max=150., const double v_min=-3000., const double v_max=3000., const int step_v=500, const double k_min=0.001, const double k_max=100., const double x_min=-3000., const double x_max=3000., const int step_x=500, const double aa=0., const bool GSL=false, const double prec=1.e-2, const std::string file_par=par::defaultString)
2D correlation function, ξ(rp,π), predicted by the Chuang & Wang model
std::vector< double > redshift_distribution_haloes(const double z_min, const double z_max, const int step_z, const double Area_degrees, const double Mass_min, const double Mass_max, const std::string model_MF, const std::string method_SS, const bool store_output=true, const std::string output_root="test", const double Delta=200, const bool isDelta_critical=false, const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true)
redshift distribution of dark matter haloes
void m_Table_Pk_CAMB_MPTbreeze(const std::string code, const bool NL, std::vector< double > &lgkk, std::vector< double > &lgPk, const double redshift, const bool store_output=true, const std::string output_root="test", const double k_max=100.) const
write and read the table where the dark matter power spectrum, computed with either CAMB or MPTbreeze...
double z_decoupling() const
redshift at wich occurs baryon photon decoupling
double Distance(const double redshift, const std::string distance_type) const
the distance at a given redshift
double OmegaNu(const double redshift=0.) const
the neutrino density at a given redshift
void set_hh(const double hh=0.7, const bool warn=true)
set the value of h
double Omega_DE() const
get the private member Cosmology::m_Omega_DE
int m_massive_neutrinos
the number of degenerate massive neutrino species
double Pk_DeltaTheta(const double kk, const std::shared_ptr< cbl::glob::FuncGrid > Pk, const double qmin, const double qmax, const double prec=1.e-3)
the real-space matter non-linear power spectrum , computed at 1-loop
std::vector< double > zeta_DM_eq(const std::vector< double > rr, const std::string model, const std::vector< double > kk, const std::vector< double > Pk_matter) const
the dark matter equilateral three-point correlation function
double OmegaM(const double redshift=0.) const
the matter density at a given redshift
double n_haloes_selection_function(const double Mass_min, const double Mass_max, const double z_min, const double z_max, const bool angle_rad, const std::string model_MF, const std::string method_SS, const std::string selection_function_file, const std::vector< int > column={}, const bool store_output=true, const std::string output_root="test", const double Delta=200, const bool isDelta_critical=false, const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true)
number of dark matter haloes per steradian or square degree, for a given redshift range and with sele...
double D_M(const double redshift) const
the comoving transverse distance at a given redshift
double zeta_ell_4_factor(const double b1, const double beta)
the multiplicative factor for , with local bias
std::vector< std::vector< double > > Nvoids(const double min_r, const double max_r, const int num_bins, const double mean_z, const double Volume, const std::string model, const double b_eff, double slope=0.854, double offset=0.420, const double deltav_NL=-0.795, const double del_c=1.69, const std::string method_Pk="EisensteinHu", const bool store_output=true, const std::string output_root="test", const std::string interpType="Linear", const double k_max=100., const std::string input_file=par::defaultString, const bool is_parameter_file=true) const
number of voids computed from the void size function model for bins of radii spaced in log scale and ...
double m_tau
: Thomson scattering optical depth due to reionization
void set_OmegaB(const double Omega_baryon=0.046)
set the value of Ωb, keeping ΩCDM=ΩM-Ωb
double sigma8_Pk(const std::string method_Pk, const double redshift, const bool store_output=true, const std::string output_root="test", const bool NL=0, const double k_min=0.001, const double k_max=100., const double prec=1.e-2, const std::string file_par=par::defaultString) const
the dark matter rms mass fluctuation within 8 Mpc/h
double OmegaDE(const double redshift=0.) const
the dark energy density at a given redshift
double DD(const double redshift) const
the amplitude of the growing mode at a given redshift,
double value(const CosmologicalParameter parameter) const
get the private member specified by the enum CosmologicalParameter
std::vector< double > zeta_DM(const double r1, const double r2, const std::vector< double > theta, const std::string model, const std::vector< double > kk, const std::vector< double > Pk_matter) const
the dark matter three-point correlation function
double g_k(const double k, const std::shared_ptr< cbl::glob::FuncGrid > PkLin, const double qmin, const double qmax, const double prec=1.e-3)
function used to estimate the non-linear power spectrum
std::vector< double > Pk_matter_NoWiggles_bspline(const std::vector< double > kk, const std::vector< double > PkLin, const std::vector< double > PkApprox, const int order, const int nknots)
the dark matter power spectrum without BAO wiggles.
double Pk_DeltaDelta_fitting_function(const double kk, const std::string method_Pk, const double redshift, const std::string author, const bool store_output, const std::string output_root, const int norm, double k_min, double k_max, const double prec, const std::string file_par, const bool unit1)
the non-linear dark matter power spectrum using fitting functions given by Bel et....
std::vector< std::vector< double > > XiMultipoles_covariance(const int nbins, const double rMin, const double rMax, const double nn, const double Volume, const std::vector< double > kk, const std::vector< double > Pk0, const std::vector< double > Pk2, const std::vector< double > Pk4, const int IntegrationMethod=1)
the covariance matrix of the first three non-null multipole moments of the two-point correlation func...
void xi_r_n(std::vector< double > &xi_n, const std::vector< double > rr, const int nn, const std::vector< double > kk, const std::vector< double > Pk)
compute the power spectrum integral transform
double z_eq() const
redshift of matter-dark energy equality
double sigma8_interpolated(const double redshift) const
σ8
double concentration_NFW_Duffy(const double Mass, const double redshift, const std::string halo_def="vir") const
the halo concentration-mass relation for NFW prfile and Duffy model
static const std::string defaultString
default std::string value
static const double pi
: the ratio of a circle's circumference to its diameter
static const double alpha
: the fine-structure constant
static const double cc
: the speed of light in vacuum (the value is exact) [km/sec]
static const double TCMB
: the present day CMB temperature [K]
std::vector< std::string > CosmologicalParameterNames()
return a vector containing the CosmologicalParameter names
std::vector< std::string > CosmologicalModelNames()
return a vector containing the CosmologicalModel names
CosmologicalParameter
the cosmological parameters
@ _sigma8_
: the power spectrum normalisation
@ _massive_neutrinos_
the number of degenerate massive neutrino species
@ _hh_
: the Hubble constant at z=0 divided by 100
@ _H0_
: the Hubble constant at z=0 [km/sec/Mpc]
@ _Omega_baryon_h2_
: the baryon density times at z=0
@ _ln_scalar_amp_
: the logarithm of 1e10 times the initial scalar amplitude of the power spectrum
@ _Omega_radiation_
: the radiation density at z=0
@ _massless_neutrinos_
: the effective number (for QED + non-instantaneous decoupling)
@ _Omega_matter_LCDM_
: the density of baryons, cold dark matter and massive neutrinos (in units of the critical density) a...
@ _fNL_
: the non-Gaussian amplitude
@ _neutrino_mass_
the total neutrino mass
@ _scalar_pivot_
the scalar pivot k in
@ _Omega_baryon_
: the baryon density at z=0
@ _scalar_amp_
: the initial scalar amplitude of the power spectrum
@ _wa_
: the parameter of the dark energy equation of state (CPL parameterisation)
@ _Omega_matter_
: the density of baryons, cold dark matter and massive neutrinos (in units of the critical density) a...
@ _w0_
: the parameter of the dark energy equation of state (CPL parameterisation)
@ _tau_
: Thomson scattering optical depth due to reionization
@ _Omega_DE_
: the dark energy density at z=0
@ _Omega_neutrinos_
: the density of massive neutrinos at z=0
@ _n_spec_
: the primordial spectral index
CosmologicalParameter CosmologicalParameterCast(const int cosmologicalParameterIndex)
cast an enum of type CosmologicalParameter from its index
CosmologicalModel
built-in cosmological models
@ _Planck15_TT_
Planck collaboration 2015, paper XIII: Table 4, TT+lowP+lensing.
@ _WMAP7_
Komatsu et al. 2011: Table 1, WMAP Seven-year Mean.
@ _Planck13_
Planck collaboration 2013, paper XVI: Table 3, Planck+WP.
@ _WMAP5_
Komatsu et al. 2009: Table 1, WMAP 5 Year Mean.
@ _Planck18_
Planck collaboration 2018, Paper VI: Table 2, TT,TE,EE+lowE+lensing.
@ _Planck15_
Planck collaboration 2015, paper XIII: Table 4, TT,TE,EE+lowP+lensing.
@ _WMAP9_
Hinshaw et al. 2013: Table 3, WMAP-only Nine-year.
CosmologicalModel CosmologicalModelCast(const int cosmologicalModelIndex)
cast an enum of type CosmologicalModel from its index
std::string CosmologicalParameter_name(const CosmologicalParameter parameter)
name of the cosmological parameter
double bias(const double Mmin, const double sigmalgM, const double M0, const double M1, const double alpha, const std::shared_ptr< void > inputs)
the mean galaxy bias
The global namespace of the CosmoBolognaLib
T Mass(const T RR, const T Rho)
the mass of a sphere of a given radius and density
std::vector< double > Xi2(const std::vector< double > rr, const std::vector< double > kk, const std::vector< double > Pk2, const double k_cut=0.58, const double cut_pow=4, const int IntegrationMethod=1)
function to obtain the two point correlation function quadrupole
double Sigma(const std::vector< double > vect)
the standard deviation of a std::vector
CoordinateUnits
the coordinate units
@ _degrees_
angle in degrees
void WarningMsgCBL(const std::string msg, const std::string functionCBL, const std::string fileCBL)
internal CBL warning message
1.9.1