Sky map overlap integral can now use MOC (multi-ordered) GW sky maps; fixes #19

This MR adds the ability to handle multi-ordered (MOC) sky maps, compared to be only able to use flat sky maps before. In addition, can handle single pixel secondary sky maps (as would be applicable for Swift, since the localization is on the order or less than a single pixel).

This MR enables the ability to handle the following combinations of sky maps:

1. MOC-MOC
2. MOC-flat
3. flat-flat (previous method)
4. MOC-ra/dec
5. flat-ra/dec

Bin scripts have been updated with these new options.

Appropriate unittests were added to verify these methods give identical results.

These changes were also tested and verified with Gracedb Events (e.g. https://gracedb-playground.ligo.org/superevents/S220317ch/view/).

ligo/raven/search.py

@@ -264,7 +264,6 @@ def skymap_overlap_integral(se_skymap, exttrig_skymap=[],
             # Find ra/dec of external skymap
             level, ipix = ah.uniq_to_level_ipix(ext_skymap_uniq)
             nsides = ah.level_to_nside(level)
-            ext_areas = ah.nside_to_pixel_area(nsides)
             ra_ext, dec_ext = \
                 ah.healpix_to_lonlat(ipix, nsides,
                                      order=ext_order)
@@ -274,7 +273,7 @@ def skymap_overlap_integral(se_skymap, exttrig_skymap=[],
             ext_ind, d2d, d3d = c.match_to_catalog_sky(catalog)
 
             skymap_overlap_integral = \
-                np.sum(se_skymap / u.sr * gw_areas *
+                np.sum(se_skymap / u.sr * areas *
                        exttrig_skymap[ext_ind] / u.sr /
                        sky_prior).to(1).value