Python Cosmo.sigma_Rz Examples

Example #1

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

Example #2

File: TWOmpz_calculate_Dg_planck_chains.py Project: fbianchini/2MPZ_DG

							   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):

Example #3

File: TWOmpz_calculate_Dg.py Project: fbianchini/2MPZ_DG

# 	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)