45 vector<double>
cbl::cosmology::Cosmology::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,
const std::string output_root,
const double Delta,
const std::string interpType,
const int norm,
const double k_min,
const double k_max,
const double prec,
const std::string method_SS,
const std::string input_file,
const bool is_parameter_file)
47 double fact = (m_unit) ? 1 : m_hh;
48 double D_N = DN(redshift);
50 vector<double> MASS(
Mass.size(), 0), MF(
Mass.size(), 0);
51 for (
size_t i=0; i<
Mass.size(); i++) {
52 MASS[i] =
Mass[i]*fact;
53 MF[i] = m_MF_generator(MASS[i],
Sigma[i], Dln_Sigma[i], redshift, D_N, model_MF, Delta)*pow(fact, 4.);
57 for (
size_t i=0; i<
Mass.size(); i++)
58 MF[i] *= MF_correction(MASS[i], redshift, method_SS, store_output, output_root, interpType, norm, k_min, k_max, prec, input_file, is_parameter_file);
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)
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
double Sigma(const std::vector< double > vect)
the standard deviation of a std::vector