]
excl = [2, 6, 7, 8, 14] # omega_k free 4,5
excl += [i for i in range(len(pnames)) if "pk" in pnames[i]]
F, lbls = rf.combined_fisher_matrix(F_list[:l],
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,
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)
# Apply 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)
# Add Planck H_0 prior
#if 'H_0' in labels[k]:
# ph = lbls.index('h')
# Fpl[ph, ph] += 1./(0.012)**2.
