Planck lensing likelihood corrections (arXiv:1502.01591)
Posted: November 24 2015
Hi,
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,
to
because that's what I thought would mess with the rest of the code (and data) the least.
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:
Note that since 1502.01591 specifies H0 instead of theta, I used CAMB to convert \{others,H0\} to \{others,theta\}. Also I've set param[calPlanck]= 1 .
Thanks,
Mat
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