CosmoMC: help with DES likelihood and priors

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
Graeme Addison
Posts: 27
Joined: July 17 2014
Affiliation: Johns Hopkins University

CosmoMC: help with DES likelihood and priors

Post by Graeme Addison » April 17 2019

I was looking at the DES files (batch3/DES.ini, data/DES/*) in the master branch and am having some trouble figuring out whether/how I can reproduce the DES or DES+Planck lensing results shown in the Planck 2018 papers. For example, what do I need to set as priors to reproduce the green, grey, and red contours from Fig 19 of Planck 2018 VI. https://arxiv.org/abs/1807.06209? In footnote 22 on page 27 it says "Here we assume consistent (DES) priors for DES and CMB lensing results". What are these priors? DES did not fit using the CosmoMC params (ombh2, omch2, theta etc.), right?

Sorry if there's some documentation on this and I missed it...

Antony Lewis
Posts: 1485
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: CosmoMC: help with DES likelihood and priors

Post by Antony Lewis » April 17 2019

You want to use

Code: Select all

batch3/DES_astro_priors.ini
This switches the CosmoMC parameterization to the "astro" parameterization which can mimic what DES did.

Graeme Addison
Posts: 27
Joined: July 17 2014
Affiliation: Johns Hopkins University

Re: CosmoMC: help with DES likelihood and priors

Post by Graeme Addison » April 17 2019

Ah, thank you, I did indeed miss that.

Is there anywhere I could find a .covmat to use as a proposal matrix for the astro parameterization with DES? If not, and it's not too much trouble, could you send me one?

Graeme Addison
Posts: 27
Joined: July 17 2014
Affiliation: Johns Hopkins University

Re: CosmoMC: help with DES likelihood and priors

Post by Graeme Addison » April 18 2019

Two more questions (posting here instead of replying to your email, Antony, in case this is useful for other people reading the forum):

1) Does the DES likelihood include the corrections to the covariance matrix described in Appendix C of the DES https://arxiv.org/abs/1708.01530? I guess yes, since just from .likestats for the full 3x2pt shear+galaxy-galaxy+galaxy clustering I get a chi-squared around 500, and in the paper they say updating the covariance takes them from 572 to 497.

2) What's the best way to make a plot showing the DES data points and one or more theory curves (e.g. show two curves with two different theory models)? I see there's a function that seems to write the theory prediction called WL_WriteLikelihoodData in source/wl.f90. Do I just need to set test_output_root and run action=4?

Antony Lewis
Posts: 1485
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: CosmoMC: help with DES likelihood and priors

Post by Antony Lewis » April 19 2019

1. Yes, it should have their final published covariance

2. I think you can do that. You can also use the python version under python/planck/DES.py with something like

Code: Select all

  
 from planck.DES import DES_like
 DESlens = DES_like(DESdataset_file, dataset_params={'used_data_types': 'xip xim'})
 lens_theory = DESlens.get_theory_for_params(param_dictionary)
 fig, axs = plt.subplots(DESlens.nzbins, DESlens.nzbins, figsize=(figsize, figsize), sharex='col', sharey='row')  
 for ax in axs.reshape(-1): ax.axis('off')
 for f1, f2 in DESlens.bin_pairs[0]:
 	ax = axs[f2, f1]
 	ax.axis('on')
        xip = DESlens.data_arrays[0][f1, f2]
        fac = 1e4 * DESlens.theta_bins
        xip = xip * fac
        ax.errorbar(DESlens.theta_bins, xip, fac * DESlens.errors[0][f1, f2], color='C0', fmt='.')
        ax.semilogx(DESlens.theta_bins, fac * lens_theory[0][f1, f2], color='k', ls='--')  
(the full script I used to make the figure in the Planck parameter paper should be added to cosmomc when the +DES chains are finally published with the likelihoods)

Graeme Addison
Posts: 27
Joined: July 17 2014
Affiliation: Johns Hopkins University

Re: CosmoMC: help with DES likelihood and priors

Post by Graeme Addison » April 19 2019

Oh, neat, I haven't been paying attention to the CosmoMC python stuff. Thanks! For the param_dictionary in this line is there a quick way to read a set of the DES nuisance params from DES chain files files (e.g. a .likestats or .minimum file)?

Code: Select all

lens_theory = DESlens.get_theory_for_params(param_dictionary)

Antony Lewis
Posts: 1485
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: CosmoMC: help with DES likelihood and priors

Post by Antony Lewis » April 19 2019

You can use e.g.

Code: Select all

from getdist.types import BestFit
param_dic=BestFit(r'chain_name.minimum', want_fixed=True).getParamDict()

Post Reply