Software. Likelihood code for the combined analysis of BOSS DR12 and Planck with COSMOMC
Here we provide the likelihood used to get the results of the combination of Planck 2015 and BOSSdr12 data for the paper:
The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: towards a computationally efficient analysis without informative priors. (Pellejero-Ibañez et al. 2017)
The main advantage of this approximated likelihood is that it is about 60 times faster than running full likelihood of Planck + LSS BOSS ones and that it contains correlations between CMB and LSS clustering probes. In the following plots we can see that this likelihood reproduces the projections of the parameter space.
We refer to the paper for further information.
Download
The code can be downloaded here, as tar file.
How to compile the likelihood
We describe here how to make the code work in COSMOMC (by Anthony Lewis).
- You will need to have a precompiled version of COSMOMC in your machine (the code works on the November 2016 version but can also be used in previous versions until those of 2014).
- Get into your cosmomc folder. Copy all the *.f90 files and the Makefile in the source/ subfolder.
- Copy the params_CMB_DoubProbe.paramnames into the paramnames/ folder and replace it with params_CMB.paramnames (command line: “mv params_CMB_DoubProbe.paramnames params_CMB.paramnames ”)
- Do “make” on the command line (in case it doesn't work, do it in the source/ subfolder).
How to run the likelihood
The COSMOMC code is now using a modify version of the MGS likelihood computed in the source/bao.f90 file so, in order to run the code, switch on the BAO.ini likelihood in the batch2/.
- We provide a .ini file (double_probe.ini) that will be running the likelihood provided. This will run it with the covariance matrix given in table 6 of the paper (file name: Imp_Samp_Double_probe_2z_nosystematcs_planck15_owmnu_map_new_params.covmat) and with the data given in table 5.
- Note that if you want to switch off the computation of all the likelihoods corresponding to BAO and keep only the one corresponding to this work (Planck + BOSSdr12) please, comment out the following lines at batch2/BAO.ini:
- #DEFAULT(BAODR11.ini)
- #bao_dataset[6DF] = data/sdss_6DF_bao.dataset
- In case you want to add back MGS likelihood you only need to open source/bao_DoubProbe.f90 and uncomment the chi2 computation of MGS beginning at: and change the line
- !Uncoment to recover MGS likelihood
- To run it just type the command:
- ./cosmomc double_probe.ini
For additional questions, please contact Marcos Pellejero (mpi[AT]iac.es).