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/).

bin/raven_calc_signif_gracedb

@@ -57,8 +57,18 @@ parser = OptionParser(
             help="Gracedb url"),
         Option("-f", "--far_grb", metavar="1e-6", default=None,
             help="False alarm rate of GRB"),
-        Option("-r", "--em_rate", metavar="1e-3", default=None,
-            help="False alarm rate of GRB")
+        Option("-n", "--em_rate", metavar="1e-3", default=None,
+            help="False alarm rate of GRB"),
+        Option("-m", "--se_moc", dest="se_moc", action="store_true",
+            help="Assume the superevent sky map is multi-ordered (MOC)."),
+        Option("-M", "--ext_moc", dest="ext_moc", action="store_true",
+            help="Assume the external sky map is multi-ordered (MOC)."),
+        Option("-r", "--se_ring", dest="se_ring", action="store_true",
+            help="Assume the superevent map uses RING ordering rather than nested."),
+        Option("-R", "--ext_ring", dest="ext_ring", action="store_true",
+            help="Assume the external sky map uses RING ordering rather than nested."),
+        Option("-c", "--use_radec", dest="use_radec", action="store_true",
+            help="Choose to use RA and dec of external sky map if a single pixel sky map.")
     ])
 opts, args = parser.parse_args()
 print(opts)
@@ -100,7 +110,7 @@ if opts.em_rate:
 else:
     em_rate = None
 
-if se_fitsfile and ext_fitsfile:
+if se_fitsfile and (ext_fitsfile or opts.use_radec):
     incl_sky = True
 else:
     incl_sky = False
@@ -110,5 +120,7 @@ results = search.calc_signif_gracedb(
               se_id, exttrig_id, tl, th, grb_search=grb_search,
               se_fitsfile=se_fitsfile, ext_fitsfile=ext_fitsfile,
               incl_sky=incl_sky, gracedb=gracedb, far_grb=far_grb,
-              em_rate=em_rate)
+              em_rate=em_rate, se_moc=opts.se_moc, ext_moc=opts.ext_moc,
+              se_nested=not opts.se_ring, ext_nested=not opts.ext_ring,
+              use_radec=opts.use_radec)
 print(results)