def GetCrossSpectrumMagLim(cosmo,
a,
z0,
b,
bias=1,
alpha=1.,
lmax=1000,
sigma_zph=0.,
compute_at_z0=False):
DNDZ = dNdzMagLim(a, z0, b, sigma_zph=sigma_zph)
limb = Limber(cosmo, lmin=0, lmax=lmax, compute_at_z0=compute_at_z0)
gals = GalsTomo(cosmo, DNDZ, b=bias, alpha=alpha)
return limb.GetCl(gals, k2=LensCMB(cosmo))
def GetAutoSpectrumMagLim(cosmo,
a,
z0,
b,
bias=1,
alpha=1.,
lmax=1000,
sigma_zph=0.,
compute_at_z0=False):
# DNDZ = dNdzInterpolation(z, dndz, nbins=1, z_min=z[0], z_max=z[-1])
DNDZ = dNdzMagLim(a, z0, b, sigma_zph=sigma_zph, z_min=0, z_max=1)
limb = Limber(cosmo, lmin=0, lmax=lmax, compute_at_z0=compute_at_z0)
gals = GalsTomo(cosmo, DNDZ, b=bias, alpha=alpha)
return limb.GetCl(gals)
