54 const int nparameters = 3;
56 vector<statistics::ParameterType> parameterType = {statistics::ParameterType::_Base_, statistics::ParameterType::_Base_, statistics::ParameterType::_Base_};
58 vector<string> parameterName = {
"alpha",
"f*sigma8",
"b*sigma8"};
63 vector<statistics::PriorDistribution> priors = {alpha_prior, fsigma8_prior, bsigma8_prior};
66 m_data_model->poly_order = 0;
81 coutCBL <<
"Setting up the fiducial model for the projected correlation function of the dark matter" << endl;
83 const vector<double> rad =
logarithmic_bin_vector(m_data_model->step, max(m_data_model->r_min, 1.e-4), min(m_data_model->r_max, 100.));
84 vector<double> wpDM(m_data_model->step);
86 for (
size_t i=0; i<(size_t)m_data_model->step; i++)
87 wpDM[i] = m_data_model->cosmology->wp_DM(rad[i], m_data_model->method_Pk, m_data_model->NL, m_data_model->redshift, m_data_model->pi_max,
true, m_data_model->output_root, m_data_model->norm, m_data_model->r_min, m_data_model->r_max, m_data_model->k_min, m_data_model->k_max, m_data_model->aa, m_data_model->GSL, m_data_model->prec, m_data_model->file_par);
89 if(!m_data_model->store_output)
90 m_data_model->cosmology->remove_output_Pk_tables(m_data_model->method_Pk, m_data_model->NL, m_data_model->redshift, m_data_model->output_root);
92 m_data_model->func_xi = make_shared<glob::FuncGrid>(
glob::FuncGrid(rad, wpDM,
"Spline"));
105 set_fiducial_sigma();
108 const int nparameters = 5;
110 vector<statistics::ParameterType> parameterType = {statistics::ParameterType::_Base_};
112 vector<string> parameterName(nparameters);
113 parameterName[0] =
"Mmin";
114 parameterName[1] =
"sigmalgM";
115 parameterName[2] =
"M0";
116 parameterName[3] =
"M1";
117 parameterName[4] =
"alpha";
119 vector<statistics::PriorDistribution> priors = {Mmin_prior, sigmalgM_prior, M0_prior, M1_prior, alpha_prior};
125 m_model = make_shared<statistics::Model1D>(
statistics::Model1D(&
wp_HOD, nparameters, parameterType, parameterName, m_data_HOD));
#define coutCBL
CBL print message.
The class Modelling_TwoPointCorrelation_projected.
void set_fiducial_wpDM()
set the fiducial model for the dark matter projected correlation function
void set_model_HOD(const statistics::PriorDistribution Mmin_prior={}, const statistics::PriorDistribution sigmalgM_prior={}, const statistics::PriorDistribution M0_prior={}, const statistics::PriorDistribution M1_prior={}, const statistics::PriorDistribution alpha_prior={})
set the HOD parameters used to model the full shape of the projected two-point correlation function
void set_model_linearBias(const statistics::PriorDistribution bsigma8_prior)
set the model to fit the projected two-point correlation function assuming a linear bias
The class PriorDistribution.
std::vector< double > xi0_linear(const std::vector< double > rad, const std::shared_ptr< void > inputs, std::vector< double > ¶meter)
model for the monopole of the two-point correlation function
std::vector< double > wp_HOD(const std::vector< double > rp, const std::shared_ptr< void > inputs, std::vector< double > ¶meter)
HOD model of the projected two-point correlation function.
The global namespace of the CosmoBolognaLib
std::vector< T > logarithmic_bin_vector(const size_t nn, const T min, const T max)
fill a std::vector with logarithmically spaced values