In [1]:
%load_ext autoreload
%autoreload 2
from __future__ import print_function
import numpy as np
from orphics import io,stats
In [13]:
# Load fisher matrix
F = np.loadtxt("../../pyfisher/data/Feb26_FisherMat_Planck_notau_lens_fsky0.6_lcdm.csv")
# Define parameter order
plist = ["H0","ombh2","omch2","tau","As","ns"]
# Define fiducials (centers for plots)
fdict = dict(zip(plist,[70.,0.01,0.1,0.06,2.2e-9,0.965]))
print(F.shape)
# Initialize plotter
fp = io.FisherPlots()
# add a section with its latex labels
fp.addSection("default",plist,paramLatexList=plist,fidDict=fdict)
# add a fisher matrix to the section
fp.addFisher("default","set1",F)
# start figure
fp.startFig()
# plot a pair of parameters
fp.plotPair("default",["H0","omch2"],["set1"])
In [ ]: