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 += [