3.3. The datasets
Observations available for comparison with model predictions are
registered in the
EPIC/cosmology/observational_data/available_observables.ini
. They
are separated in sections (Hz
, H0
, SNeIa
, BAO
and
CMB
), which contain the name of the CosmologicalSetup
’s module for the
observable theoretical calculation (predicting_function_name
).
Each dataset goes below one of the sections, with the folder or text file
containing the data indicated.
The path is relative to the EPIC/cosmology/observational_data/
folder.
If the folder name is the same as the observable
label it can be ommited. Besides, the Dataset
subclasses are defined
at the beginning of the .ini
file. Each of these classes has its own
methods for initialization and likelihood evaluation.
We choose the datasets by passing a dictionary with observables as keys and
datasets as values to the function choose_from_datasets
from
observations
in EPIC.cosmology
:
The WiggleZ
dataset is available but is commented out because it is
correlated with the SDSS datasets and thus not supposed to be used together
with them.
Refer to the papers published by the authors of the observations for details.
Other incompatible combinations are defined in the
conflicting_dataset_pairs.txt
file in
EPIC/cosmology/observational_data/
.
The code will check for these conflicts prior to proceeding with an analysis.
The returned flattened dictionary of Dataset objects will later be passed to a
MCMC analysis.
Now I describe briefly the datasets made available by the community.
Type Ia supernovae
Two types of analyses can be made with the JLA catalogue.
One can either use the full likelihood (JLA_full
) or a simplified version
based on 30 redshift bins (JLA_simplified
).
Here we are using the binned data consisting of distance modulus estimates at
31 points (defining 30 bins of redshift).
If you want to use the full dataset (which makes the analysis much slower since
it involves three more nuisance parameters and requires the program to invert a
740 by 740 matrix at every iteration for the calculation of the JLA
likelihood), you need to download the covariance matrix data
(covmat_v6.tgz
) from
http://supernovae.in2p3.fr/sdss_snls_jla/ReadMe.html. The covmat
folder must be extracted to the jla_likelihood_v6
folder.
This is not included in EPIC because the data files are too big.
Either way, the data location is the same data folder jla_likelihood_v6
.
Note that the binned dataset introduces one nuisance parameter M
,
representing an overall shift in the absolute magnitudes, and the full dataset
introduces four nuisance parameters related to the light-curve parametrization.
See Betoule et al. (2014) for more details.
CMB distance priors
Constraining models with temperature or polarization anisotropy amplitudes is
not currently implemented.
However, you can include the CMB distance priors from Planck2015
or the updated priors from Planck2018.
The datasets consist of an acoustic scale \(l_A\), a shift parameter
\(R\) and the physical density of baryons \(\Omega_{b0}h^2\).
You can choose between the data for \(\Lambda\text{CDM}\) and
\(w\text{CDM}\) with either Planck2015_distances_LCDM
,
Planck2015_distances_wCDM
, Planck2018_distances_LCDM
, or
Planck2018_distances_wCDM
;
Planck2015_distances_LCDM+Omega_k
Planck2018_distances_LCDM+Omega_k
are also available for when curvature is
supported.
BAO data
Measurements of the baryon acoustic scales from the Six Degree Field Galaxy
Survey (6dF) combined with the most recent data releases of Sloan Digital Sky
Survey (SDSS-MGS), the LOWZ and CMASS galaxy samples of the
Baryon Oscillation Spectroscopic Survey (BOSS-LOWZ and BOSS-CMASS),
data from
the Quasar-Lyman \(\alpha\) cross-correlation,
the distribution of the Lyman \(\alpha\) forest in BOSS
and from the
WiggleZ Dark Energy Survey are available in the BAO
folder,
as well as the latest consensus from the completed SDSS-III BOSS survey.
The observable is based on the value of
the characteristic ratio \(r_s(z_d)/D_V(z)\) between
the sound horizon \(r_s\) at decoupling time (\(z_d\)) and the
effective BAO distance \(D_V\), or some variation of this.
The respective references are given in the headers of the data files.
\(H(z)\) data
These are the cosmic chronometer data.
30 measurements of the Hubble expansion rate \(H(z)\) at redshifts between
0 and 2.
The values of redshift, \(H\) and the uncertainties are given in the file
Hz/Hz_Moresco_et_al_2016.txt
.
\(H_0\) data
The \(2.4\%\) precision local measure of \(H_0\) is present in
H0/local_H0_Riess_et_al_2016.txt
.
Footnotes