Bad Contours with CMB and BAO likelihoods
Posted: March 19 2019
Hello everyone,
I'm new at using cosmomc and I have some problems with it. I'm currently using the 2015 version. My goal is to produce the plot of w - H0 for CMB and BAO likelihoods using a certain w(α) model and parametrization.
More specifically the test.ini file I'm using is,
And the params_CMB_defaults file has the following form (since I wanted to vary some fixed parameters such as w0),
I also made changes to the files equations_ppf and halofit_ppf in camb, setting the new w and density equations in both of them. However the results were less than anticipated, and had the form,
I'm open to any and all suggestions as to where my mistake is.
EDIT: Could it be because the default sample number is not enough for 2 likelihoods?
Thank you in advance
I'm new at using cosmomc and I have some problems with it. I'm currently using the 2015 version. My goal is to produce the plot of w - H0 for CMB and BAO likelihoods using a certain w(α) model and parametrization.
More specifically the test.ini file I'm using is,
Code: Select all
#general settings
#Bicep-Keck-Planck, varying cosmological parameters
#DEFAULT(batch2/BKPlanck.ini)
#Planck 2015, default just include native likelihoods (others require clik)
DEFAULT(batch2/plik_dx11dr2_HM_v18_TT.ini)
DEFAULT(batch2/lowTEB.ini)
#DEFAULT(batch2/lowl.ini)
#DEFAULT(batch2/lensing.ini)
#Other Likelihoods
DEFAULT(batch2/BAO.ini)
#DEFAULT(batch2/WiggleZ_MPK.ini)
#DEFAULT(batch2/MPK.ini)
#DEFAULT(batch2/WL.ini)
#general settings
DEFAULT(batch2/common.ini)
#e.g. to vary r in addition to standard 6:
#(for r>0 also need compute_tensors=T)
#compute_tensors = T
#param[r] = 0.03 0 2 0.04 0.04
#high for new runs
MPI_Max_R_ProposeUpdate = 30
propose_matrix= planck_covmats/base_BAO_TT_lowTEB_plik.covmat
Code: Select all
#param[omegabh2] = 0.0221
param[omegabh2] = 0.023 0.013 0.033 0.0001 0.0001
#param[omegach2] = 0.12
param[omegach2] = 0.12 0.001 0.99 0.001 0.0005
#param[theta] = 1.0411
param[theta] = 1.0411 0.5 10 0.0004 0.0002
#param[tau] = 0.09
param[tau] = 0.09 0.01 0.8 0.01 0.005
num_massive_neutrinos=1
param[mnu] = 0.06
param[meffsterile] = 0
param[omegak] = 0
param[w] = -1 -2 0 0.001 0.001
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
param[ns] = 0.96 0.7 1.3 0.004 0.002
#log[10^10 A_s]
#param[logA] = 3.1
param[logA] = 3.1 1.7 5 0.001 0.001
I'm open to any and all suggestions as to where my mistake is.
EDIT: Could it be because the default sample number is not enough for 2 likelihoods?
Thank you in advance