I'm trying to verify the statement that the corrections to the Planck lensing likelihood (detailed in C.5 of Planck Lensing 2015 1502.01591) should be zero when the input cosmology is equal to the fiducial cosmology used in the lensing analysis.
To do this, I use the "test likelihood" function in CosmoMC and enable only "lensing.ini". I then set the cosmological parameters to the fiducials specifed in 1502.01591 . I run the likelihood test and it tells me Planck lensing provides a loglike of 5.3.
Next, I go into source/CMBLikes.f90 and comment out the correction by changing, inside GetBinnedMapCls,
Code: Select all
Cls = Cls + (correctionCl - this%FiducialCorrection(:,bin))
Code: Select all
Cls = Cls
I run the CosmoMC test again and it tells me Planck lensing provides a loglike of 4.67 .
So my question is, shouldn't these two loglikes be exactly equal, if not very close?
The input parameters I'm using are:
Code: Select all
param[omegabh2] = 0.0222
param[omegach2] = 0.1203
param[omeganh2] = 0.00064
param[theta] = 1.041281
param[tau] = 0.065
num_massive_neutrinos=1
param[mnu] = 0.06
param[meffsterile] = 0
param[omegak] = 0
param[w] = -1
param[nrun] = 0
param[nrunrun] = 0
param[r] = 0
param[wa] = 0
param[nnu] = 3.046
param[yhe] = 0.24
param[alpha1] = 0
param[deltazrei] = 0.5
param[Alens] = 1
param[Alensf]=-1
param[fdm] = 0
param[ns] = 0.96
#log[10^10 A_s]
param[logA] = 3.04
param[Aphiphi] = 1
param[calPlanck]= 1
Thanks,
Mat