Combining two CMB lensing likelihoods

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Federico Bianchini
Posts: 11
Joined: June 25 2013
Affiliation: University of Melbourne
Contact:

Combining two CMB lensing likelihoods

Hi Antony and all,

I'm trying to combine two CMB lensing likelihoods, the one from Planck 2018 and another one based on the Gaussian likelihood (including linear corrections) in CMBlikes.f90.

To this end, I simply point to the two likelihoods in my CosmoMC parameter files, i.e.

Code: Select all

DEFAULT(CosmoMC/batch3/lensing.ini)
DEFAULT(CosmoMC/batch3/my_very_own_lensing.ini)


However, the second dataset is not included when I run the sampler (in fact, it doesn't appear in the .likelihoods file produced during the run).

I'm assuming the cmb_dataset[lensing] instruction is not read twice but I couldn't exactly pinpoint the source of the problem.

Any feedback would be greatly appreciated!

Cheers,
~federico

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

Re: Combining two CMB lensing likelihoods

I think you need to use something like

Code: Select all

cmb_dataset[my_lensing] = ...
(the thing in brackets is just a name tag for the likelihood; though note you also need to add other relevant xxx[lensing] input parameters to xxx[my_lensing] accordingly)

Federico Bianchini
Posts: 11
Joined: June 25 2013
Affiliation: University of Melbourne
Contact:

Re: Combining two CMB lensing likelihoods

Thanks, that worked like a charm!

~federico