Commit 11f41489 authored by Leo Pound Singer's avatar Leo Pound Singer

Not a frustum [ci skip]

parent bb620b63
......@@ -226,8 +226,9 @@ def find_injection_moc(sky_map, true_ra=None, true_dec=None, true_dist=None,
d_r = max_r / n_r
r = d_r * np.arange(1, n_r)
# Calculate volume of frustum-shaped voxels with distance centered on r
# and extending from (r - d_r) to (r + d_r).
# Calculate volume of each voxel, defined as the region within the
# HEALPix pixel and contained within the two centric spherical shells
# with radii (r - d_r) and (r + d_r).
dV = (2 / 3 * (3 * np.square(r) + np.square(d_r))
* d_r * dA.reshape(-1, 1))
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment