add option to turn off observability check

skyportal · Nov 24, 2023 · b96d25e · b96d25e
1 parent 235c2d5
commit b96d25e
Show file tree

Hide file tree

Showing 3 changed files with 11 additions and 0 deletions.
diff --git a/gwemopt/args.py b/gwemopt/args.py
@@ -145,4 +145,6 @@ def parse_args(args):
         "--galactic_limit", help="Galactic limit.", default=0.0, type=float
     )
 
+    parser.add_argument("--ignore_observability", action="store_true", default=False)
+
     return parser.parse_args(args=args)
diff --git a/gwemopt/params.py b/gwemopt/params.py
@@ -105,6 +105,10 @@ def params_struct(opts):
     params["exposuretimes"] = np.array(opts.exposuretimes.split(","), dtype=float)
     params["catalogDir"] = Path(opts.catalogDir)
 
+    params["ignore_observability"] = (
+        opts.ignore_observability if hasattr(opts, "ignore_observability") else False
+    )
+
     if params["catalog"] is not None:
         params["galaxy_limit"] = int(opts.galaxy_limit)
 

diff --git a/gwemopt/segments.py b/gwemopt/segments.py
@@ -285,6 +285,11 @@ def get_segments_tiles(params, config_struct, tile_struct):
         ras.append(tile_struct[key]["ra"])
         decs.append(tile_struct[key]["dec"])
 
+    if params["ignore_observability"]:
+        for ii, key in enumerate(keys):
+            tile_struct[key]["segmentlist"] = copy.deepcopy(segmentlist)
+        return tile_struct
+
     # Convert to RA, Dec.
     radecs = astropy.coordinates.SkyCoord(
         ra=np.array(ras) * u.degree, dec=np.array(decs) * u.degree, frame="icrs"