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']