clgg_slash_binned[zbin],
lbmax=i) for i in xrange(1, clgg[zbin].size)
]
# embed()
# exit()
cw = '#E9724C'
cnow = '#255F85'
z_ = np.logspace(-4, 3, 1000)
cosmo_ = Cosmo(Giannantonio15Params)
d0 = cosmo_.D_z_norm(z_)
s8 = cosmo_.sigma_Rz(8. / cosmo_.h)
om = cosmo_.omegam
def D_z_norm(z, omegam, gamma=0.55):
if np.isscalar(z) or (np.size(z) == 1):
def func(x, omegam, gamma):
return f_z(x, omegam=omegam, gamma=gamma) / (1 + x)
return np.exp(-integrate.quad(func, 0, z, args=(
omegam,
gamma,
))[0])
else:
return np.asarray([D_z_norm(tz, omegam, gamma=gamma) for tz in z])
clkg_slash_binned[zbin], clgg_slash_binned[zbin], lbmax=i) for i in xrange(1,clgg[zbin].size)] # embed() # exit() cw = '#E9724C' cnow = '#255F85' z_ = np.logspace(-4,3,1000) cosmo_ = Cosmo(Giannantonio15Params) d0 = cosmo_.D_z_norm(z_) s8 = cosmo_.sigma_Rz(8./cosmo_.h) om = cosmo_.omegam def D_z_norm(z, omegam, gamma=0.55): if np.isscalar(z) or (np.size(z) == 1): def func(x, omegam, gamma): return f_z(x, omegam=omegam, gamma=gamma)/(1+x) return np.exp( -integrate.quad( func, 0, z, args=(omegam,gamma,))[0]) else: return np.asarray([ D_z_norm(tz, omegam, gamma=gamma) for tz in z ]) def f_z(z, omegam, gamma=0.55): # [unitless] return (omegam*(1+z)**3/E_z(z, omegam)**2)**gamma def E_z(z, omegam):
# plt.ylabel(r'$D_G$')#, size=15)
# plt.tight_layout()
# if fsky:
# plt.savefig('plots/D_G_dl'+str(delta_ell)+'_lmin_'+str(lmin)+'_lmax'+str(lmax)+'_KS_min_'+str(K_S_min)+'_KS_max_'+str(K_S_max)+'_nside256_v2_single_weighted_fsky.pdf')
# else:
# plt.savefig('plots/D_G_dl'+str(delta_ell)+'_lmin_'+str(lmin)+'_lmax'+str(lmax)+'_KS_min_'+str(K_S_min)+'_KS_max_'+str(K_S_max)+'_nside256_v2_single_weighted.pdf')
# # plt.show()
# plt.close()
# exit()
z_ = np.logspace(-4, 3, 1000)
cosmo = Cosmo(Giannantonio15Params)
d0 = cosmo.D_z_norm(z_)
s8 = cosmo.sigma_Rz(8. / cosmo.h)
om = cosmo.omegam
gamma0 = 0.25
gammaa = -1.
w = -1.2
dgamma0 = cosmo.D_z_norm(z_, gamma0=gamma0)
dgamma0 = dgamma0 * d0[-1] / dgamma0[-1] #* (cos)/(om*s8)
dgammaa = cosmo.D_z_norm(z_, gammaa=gammaa)
dgammaa = dgammaa * d0[-1] / dgamma0[-1]
params_w = Giannantonio15Params.copy()
params_w['w'] = w
cosmo_w = Cosmo(params_w)