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)