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/gracedb_events.py

@@ -26,7 +26,6 @@ __author__ = "Alex Urban <alexander.urban@ligo.org>"
 
 
 # Imports.
-import healpy as hp
 import tempfile
 
 from gracedb_sdk import Client
@@ -42,7 +41,7 @@ class ExtTrig(object):
     def __init__(self, graceid, gracedb=None, event_dict=None,
                  fitsfile=None, is_moc=True, use_radec=False, nested=True):
         self.graceid = graceid
-        self.fits = fitsfile  #c name of fits file
+        self.fits = fitsfile  # name of fits file
         self.neighbor_type = 'S'  # by default, look for superevents
         self.is_moc = is_moc
         self.use_radec = use_radec
@@ -84,7 +83,7 @@ class ExtTrig(object):
                     self.skymap = skymap_table['PROBDENSITY']
                     self.uniq = skymap_table['UNIQ']
                 else:
-                    self.skymap = fits.read_sky_map(skymapfile.name, 
+                    self.skymap = fits.read_sky_map(skymapfile.name,
                                                     moc=self.is_moc,
                                                     nest=self.nested)[0]
                     self.uniq = []