pnames = rf.load_param_names(root + "-fisher-full-0.dat")
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
]
excl = [2, 6, 7, 8, 14]
excl += [i for i in range(len(pnames)) if "pk" in pnames[i]]
F, lbls = rf.combined_fisher_matrix(F_list,
expand=zfns,
names=pnames,
exclude=excl)
# Add Planck prior
#Fpl = euclid.add_detf_planck_prior(F, lbls, info=False)
#Fpl = euclid.add_planck_prior(F, lbls, info=False)
if USE_DETF_PLANCK_PRIOR:
# DETF Planck prior
print "*** Using DETF Planck prior ***"
l2 = ['n_s', 'w0', 'wa', 'omega_b', 'omegak', 'omegaDE', 'h']
F_detf = euclid.detf_to_rf("DETF_PLANCK_FISHER.txt", cosmo)
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']