diff --git a/heracles/fields.py b/heracles/fields.py index 76c60c0..1a9eb54 100644 --- a/heracles/fields.py +++ b/heracles/fields.py @@ -375,8 +375,9 @@ async def __call__( val /= wbar # compute bias from variance (per object) - variance = var / wmean**2 - neff = ngal / (4 * np.pi * fsky) + variance = var / w2mean + deff = w2mean / wmean**2 + neff = ngal / (4 * np.pi * fsky) / deff bias = fsky * variance / neff # set metadata of array