zfns = [
'b_HI',
]
excl = [
'Tb', 'f', 'aperp', 'apar', 'DA', 'H', 'pk*', 'fs8', 'bs8', 'N_eff'
]
F, lbls = rf.combined_fisher_matrix(F_list,
expand=zfns,
names=pnames,
exclude=excl)
# Add Planck prior
print "*** Using DETF Planck prior ***"
l2 = ['n_s', 'w0', 'wa', 'omega_b', 'omegak', 'omegaDE', 'h', 'sigma8']
F_detf = euclid.detf_to_rf("DETF_PLANCK_FISHER.txt", cosmo, omegab=False)
Fpl, lbls = rf.add_fisher_matrices(F, F_detf, lbls, l2, expand=True)
# Remove extra params
Fpl, lbls = rf.combined_fisher_matrix([
Fpl,
],
expand=[],
names=lbls,
exclude=[
'omegak',
'w0',
'wa',
'gamma',
])
print lbls
F_planck = np.genfromtxt(
"/home/phil/oslo/Dropbox/HI_neutrinos/fisher_code/planck/planck_derived_fisher.dat"
).T
F_planck, lbls_planck = rf.combined_fisher_matrix([
F_planck,
],
expand=[],
names=lbls_planck,
exclude=[
'Neff',
])
# Add Planck prior
Fpl, lbls = rf.add_fisher_matrices(F,
F_planck,
lbls,
lbls_planck,
expand=True)
# Decide whether to fix various parameters
fixed_params = []
if not MARGINALISE_CURVATURE: fixed_params += [
'omegak',
]
if not MARGINALISE_INITIAL_PK: fixed_params += ['n_s', 'sigma8']
if not MARGINALISE_OMEGAB: fixed_params += [
'omega_b',
]
# Fix (w0, wa)
#fixed_params += ['w0', 'wa']
np.genfromtxt(root + "-fisher-full-%d.dat" % i) for i in range(Nbins)
]
# EOS FISHER MATRIX
pnames = rf.load_param_names(root + "-fisher-full-0.dat")
zfns = [
'b_HI',
]
excl = ['Tb', 'f', 'aperp', 'apar', 'DA', 'H', 'N_eff', 'pk*']
F, lbls = rf.combined_fisher_matrix(F_list,
expand=zfns,
names=pnames,
exclude=excl)
# Relabel galaxy bias from low-z survey and sum current survey + low-z
F, lbls = rf.add_fisher_matrices(F_lowz, F, lbls_lowz, lbls, expand=True)
print lbls
# Add Planck prior
print "*** Using DETF Planck prior ***"
l2 = ['n_s', 'w0', 'wa', 'omega_b', 'omegak', 'omegaDE', 'h', 'sigma8']
F_detf = euclid.detf_to_rf("DETF_PLANCK_FISHER.txt", cosmo, omegab=False)
Fpl, lbls = rf.add_fisher_matrices(F, F_detf, lbls, l2, expand=True)
# Add BOSS
Fpl, lbls = rf.add_fisher_matrices(F, F_detf, lbls, l2, expand=True)
# Decide whether to fix various parameters
fixed_params = []
if not MARGINALISE_CURVATURE: fixed_params += [
'omegak',
if 'Euclid' in names[k]:
F1 = F
lbl1 = copy.deepcopy(lbls)
else:
F2 = F
lbl2 = copy.deepcopy(lbls)
# Add Planck prior
if USE_DETF_PLANCK_PRIOR:
# DETF Planck prior
print "*** Using DETF Planck prior ***"
l2 = ['n_s', 'w0', 'wa', 'omega_b', 'omegak', 'omegaDE', 'h', 'sigma8']
F_detf = euclid.detf_to_rf("DETF_PLANCK_FISHER.txt",
cosmo,
omegab=False)
Fpl, lbls = rf.add_fisher_matrices(F, F_detf, lbls, l2, expand=True)
else:
# Euclid Planck prior
print "*** Using Euclid (Mukherjee) Planck prior ***"
l2 = ['n_s', 'w0', 'wa', 'omega_b', 'omegak', 'omegaDE', 'h']
Fe = euclid.planck_prior_full
F_eucl = euclid.euclid_to_rf(Fe, cosmo)
Fpl, lbls = rf.add_fisher_matrices(F, F_eucl, lbls, l2, expand=True)
# Decide whether to fix various parameters
fixed_params = []
if not MARGINALISE_CURVATURE: fixed_params += [
'omegak',
]
if not MARGINALISE_INITIAL_PK: fixed_params += ['n_s', 'sigma8']
if not MARGINALISE_OMEGAB: fixed_params += [
