Maintenance will be performed on git.ligo.org, chat.ligo.org, containers.ligo.org, and docs.ligo.org on Tuesday 7th July 2020 starting at approximately 10am PDT and lasting for around 15 minutes. There will be a short period of downtime towards the end of the maintenance window. Please direct any comments, questions, or concerns to uwm-help@cgca.uwm.edu.

Commit ab4c4d56 authored by Leo Pound Singer's avatar Leo Pound Singer

Revert "Fix error in volume element for ligo-skymap-stats"

The original Python code was correct, but the comment was wrong.

This reverts commit b436790d.
parent 11f41489
Pipeline #38606 failed with stages
in 28 minutes and 5 seconds
......@@ -14,9 +14,6 @@ Changelog
now natively output multi-resolution HEALPix files, although they still
natively output flat, fixed-resolution HEALPix files.
- Fix an error in the volume element that resulted in ``ligo-skymap-stats``
underestimating credible volumes by a factor of 2.
- Add the ``ligo-skymap-flatten`` tool to convert multi-resolution HEALPix
files to flat, fixed-resolution, implicitly indexed HEALPix files.
......
......@@ -228,9 +228,8 @@ def find_injection_moc(sky_map, true_ra=None, true_dec=None, true_dist=None,
# 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))
# with radii (r - d_r / 2) and (r + d_r / 2).
dV = (np.square(r) + np.square(d_r) / 12) * d_r * dA.reshape(-1, 1)
# Calculate probability within each voxel.
dP = probdensity.reshape(-1, 1) * dV * np.exp(
......
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