Code: Select all
roots = ['../chains1/test', '../chains2/test']
for r in roots:
samples = getdist.loadMCSamples(r)
p = samples.getParams()
omegab = p.omegabh2/ p.H0**2 * 100**2
samples.addDerived(omegab, name='omegab', label=r'\Omega_b')
samples.updateBaseStatistics()
labels=['chain1', 'chain2']
params = ['H0','ns','omegam','omegab','tau','logA']
g = plots.getSubplotPlotter(width_inch=5)
g.settings.axes_fontsize = 9
g.settings.legend_fontsize = 10
g.plots_1d(roots,params, nx = 3, legend_labels=labels)
g.export()