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