Page 1 of 1

Combining two CMB lensing likelihoods

Posted: May 21 2019
by Federico Bianchini
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

Re: Combining two CMB lensing likelihoods

Posted: May 21 2019
by Antony Lewis
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)

Re: Combining two CMB lensing likelihoods

Posted: May 30 2019
by Federico Bianchini
Thanks, that worked like a charm!

~federico