CosmoBolognaLib
Free Software C++/Python libraries for cosmological calculations
cbl::modelling::cosmo::Modelling_Cosmology Class Reference

The class Modelling_Cosmology. More...

#include <Modelling_Cosmology.h>

Inheritance diagram for cbl::modelling::cosmo::Modelling_Cosmology:
Collaboration diagram for cbl::modelling::cosmo::Modelling_Cosmology:

Public Member Functions

void set_fiducial_cosmology (const cbl::cosmology::Cosmology cosmology)
 set the fiducial cosmology More...
 
std::shared_ptr< cbl::cosmology::Cosmologyfiducial_cosmology ()
 set the fiducial cosmology More...
 
void set_cosmological_parameters (const std::vector< cbl::cosmology::CosmologicalParameter > cosmoPar_name, const std::vector< cbl::statistics::PriorDistribution > cosmoPar_prior, const std::string distance_prior=par::defaultString, const std::vector< std::string > external_dataset={})
 set the cosmological parameters More...
 
Constructors/destructors
 Modelling_Cosmology ()=default
 default constuctor
 
 Modelling_Cosmology (const std::shared_ptr< cbl::data::Data > dataset, const std::vector< std::string > data_type)
 constuctor More...
 
 ~Modelling_Cosmology ()=default
 default destructor
 
- Public Member Functions inherited from cbl::modelling::Modelling
void m_set_posterior (const int seed)
 set the interal variable m_posterior More...
 
 Modelling ()=default
 default constuctor
 
virtual ~Modelling ()=default
 default destructor
 
std::shared_ptr< data::Datadata ()
 return the dataset More...
 
std::shared_ptr< data::Datadata_fit ()
 return the dataset More...
 
std::shared_ptr< statistics::Likelihoodlikelihood ()
 return the likelihood parameters More...
 
std::shared_ptr< statistics::Posteriorposterior ()
 return the posterior parameters More...
 
std::shared_ptr< statistics::ModelParameterslikelihood_parameters ()
 return the likelihood parameters More...
 
std::shared_ptr< statistics::ModelParametersposterior_parameters ()
 return the posterior parameters More...
 
virtual void set_parameter_from_string (const std::string parameter, const double value)
 set the value of a parameter providing its name string More...
 
virtual double get_parameter_from_string (const std::string parameter) const
 get the value of a parameter providing its name string More...
 
std::shared_ptr< statistics::PriorDistributionget_prior (const int i)
 get the internal variable m_parameter_priors More...
 
std::shared_ptr< statistics::Modelget_response_function ()
 return the response function used to compute the super-sample covariance More...
 
void reset_fit_range ()
 reset the fit range More...
 
void set_fit_range (const double xmin, const double xmax)
 set the fit range More...
 
void set_fit_range (const double xmin, const double xmax, const double ymin, const double ymax)
 set the fit range More...
 
void set_data (const std::shared_ptr< data::Data > dataset)
 set the dataset More...
 
void set_likelihood (const statistics::LikelihoodType likelihood_type, const std::vector< size_t > x_index={0, 2}, const int w_index=-1, const double prec=1.e-10, const int Nres=-1)
 set the likelihood function More...
 
void set_likelihood (const cbl::statistics::Likelihood_function log_likelihood_function)
 set the likelihood function, given a user-defined log-likelihood function More...
 
void maximize_likelihood (const std::vector< double > start, const std::vector< std::vector< double >> parameter_limits, const unsigned int max_iter=10000, const double tol=1.e-6, const double epsilon=1.e-3)
 function that maximizes the posterior, finds the best-fit parameters and stores them in the model More...
 
void maximize_posterior (const std::vector< double > start, const unsigned int max_iter=10000, const double tol=1.e-6, const double epsilon=1.e-3, const int seed=666)
 function that maximizes the posterior, finds the best-fit parameters and stores them in the model More...
 
void sample_posterior (const int chain_size, const int nwalkers, const int seed=666, const double aa=2, const bool parallel=true)
 sample the posterior, initializing the chains by drawing from the prior distributions More...
 
void sample_posterior (const int chain_size, const int nwalkers, const double radius, const std::vector< double > start, const unsigned int max_iter=10000, const double tol=1.e-6, const double epsilon=1.e-3, const int seed=666, const double aa=2, const bool parallel=true)
 sample the posterior, initializing the chains in a ball around the posterior best-fit parameters values More...
 
void sample_posterior (const int chain_size, const int nwalkers, std::vector< double > &value, const double radius, const int seed=666, const double aa=2, const bool parallel=true)
 sample the posterior, initializing the chains by drawing from the prior distributions More...
 
void sample_posterior (const int chain_size, const std::vector< std::vector< double >> chain_value, const int seed=666, const double aa=2, const bool parallel=true)
 sample the posterior, initializing the chains with input values More...
 
void sample_posterior (const int chain_size, const int nwalkers, const std::string input_dir, const std::string input_file, const int seed=666, const double aa=2, const bool parallel=true)
 sample the posterior, initializing the chains reading the input values from an input file More...
 
void importance_sampling (const std::string input_dir, const std::string input_file, const int seed=666, const std::vector< size_t > column={}, const int header_lines_to_skip=1, const bool is_FITS_format=false, const bool apply_to_likelihood=false)
 perform importance sampling More...
 
void write_chain (const std::string output_dir, const std::string output_file, const int start=0, const int thin=1, const bool is_FITS_format=false, const int prec=5, const int ww=14)
 write the chains obtained after the MCMC sampling More...
 
void read_chain (const std::string input_dir, const std::string input_file, const int nwalkers, const std::vector< size_t > columns={}, const int skip_header=1, const bool fits=false)
 read the chains More...
 
void show_results (const int start=0, const int thin=1, const int nbins=50, const bool show_mode=false, const int ns=-1)
 show the results of the MCMC sampling on screen More...
 
void write_results (const std::string output_dir, const std::string root_file, const int start=0, const int thin=1, const int nbins=50, const bool fits=false, const bool compute_mode=false, const int ns=-1)
 write the results of the MCMC sampling to file More...
 
virtual void write_model (const std::string output_dir, const std::string output_file, const std::vector< double > xx, const std::vector< double > parameters)
 write the model at xx for given parameters More...
 
virtual void write_model (const std::string output_dir, const std::string output_file, const std::vector< double > xx, const std::vector< double > yy, const std::vector< double > parameters)
 write the model at xx, yy for given parameters More...
 
virtual void write_model_at_bestfit (const std::string output_dir, const std::string output_file, const std::vector< double > xx)
 write the model at xx with best-fit parameters obtained from posterior maximization More...
 
virtual void write_model_at_bestfit (const std::string output_dir, const std::string output_file, const std::vector< double > xx, const std::vector< double > yy)
 write the model at xx, yy with best-fit parameters obtained from likelihood maximization More...
 
virtual void write_model_from_chains (const std::string output_dir, const std::string output_file, const std::vector< double > xx, const int start=0, const int thin=1)
 write the model at xx computing 16th, 50th and 84th percentiles from the chains More...
 
virtual void write_model_from_chains (const std::string output_dir, const std::string output_file, const std::vector< double > xx, const std::vector< double > yy, const int start=0, const int thin=1)
 write the model at xx, yy computing 16th, 50th and 84th percentiles from the chains More...
 
double reduced_chi2 (const std::vector< double > parameter={})
 the reduced \(\chi^2\) More...
 

Protected Attributes

std::shared_ptr< cbl::cosmology::Cosmologym_cosmology
 cosmology
 
std::map< cbl::cosmology::CosmologicalParameter, int > m_map_cosmoPar
 map with cosmological parameters
 
std::vector< std::string > m_data_type
 the cosmological measurements
 
STR_data_model_cosmology m_data_model
 the container of parameters for model computation
 
- Protected Attributes inherited from cbl::modelling::Modelling
std::shared_ptr< data::Datam_data = NULL
 input data to be modelled
 
bool m_fit_range = false
 check if fit range has been set
 
std::shared_ptr< data::Datam_data_fit
 input data restricted to the range used for the fit
 
std::shared_ptr< statistics::Modelm_model = NULL
 input model
 
std::shared_ptr< statistics::Modelm_response_func = NULL
 response function for the computation of the super-sample covariance
 
std::shared_ptr< statistics::Likelihoodm_likelihood = NULL
 likelihood
 
std::vector< std::shared_ptr< statistics::PriorDistribution > > m_parameter_priors
 prior
 
std::shared_ptr< statistics::Posteriorm_posterior = NULL
 posterior
 

Additional Inherited Members

- Protected Member Functions inherited from cbl::modelling::Modelling
void m_set_prior (std::vector< statistics::PriorDistribution > prior_distribution)
 set the internal variable m_parameter_priors More...
 
void m_isSet_response ()
 check if the response function used to compute the super-sample covariance is set
 

Detailed Description

The class Modelling_Cosmology.

"Headers/Modelling_Cosmology.h"

This file defines the interface of the base class Modelling_Cosmology, used for modelling cosmological measurements

Examples
model_cosmology.cpp.

Definition at line 69 of file Modelling_Cosmology.h.

Constructor & Destructor Documentation

◆ Modelling_Cosmology()

cbl::modelling::cosmo::Modelling_Cosmology::Modelling_Cosmology ( const std::shared_ptr< cbl::data::Data dataset,
const std::vector< std::string >  data_type 
)

constuctor

Parameters
datasetthe dataset containing the input measurements to be modelled
data_typethe type of cosmological measurements to be modelled

Definition at line 63 of file Modelling_Cosmology.cpp.

Member Function Documentation

◆ fiducial_cosmology()

std::shared_ptr<cbl::cosmology::Cosmology> cbl::modelling::cosmo::Modelling_Cosmology::fiducial_cosmology ( )
inline

set the fiducial cosmology

Returns
pointer to object of class Cosmology

Definition at line 127 of file Modelling_Cosmology.h.

◆ set_cosmological_parameters()

void cbl::modelling::cosmo::Modelling_Cosmology::set_cosmological_parameters ( const std::vector< cbl::cosmology::CosmologicalParameter cosmoPar_name,
const std::vector< cbl::statistics::PriorDistribution cosmoPar_prior,
const std::string  distance_prior = par::defaultString,
const std::vector< std::string >  external_dataset = {} 
)

set the cosmological parameters

Parameters
cosmoPar_namevector of cosmological parameter names
cosmoPar_priorvector of cosmological parameter priors
distance_priorthe CMB distance prior name
external_datasetexternal dataset file name
Examples
model_cosmology.cpp.

Definition at line 82 of file Modelling_Cosmology.cpp.

◆ set_fiducial_cosmology()

void cbl::modelling::cosmo::Modelling_Cosmology::set_fiducial_cosmology ( const cbl::cosmology::Cosmology  cosmology)

set the fiducial cosmology

Parameters
cosmologythe fiducial cosmology
Examples
model_cosmology.cpp.

Definition at line 73 of file Modelling_Cosmology.cpp.


The documentation for this class was generated from the following files: