Merge pull request #54 from VERITAS-Observatory/reduce-exclusion-list

fix example; print out max wobble distance

Author: GernotMaier

examples/reflected_region.yml

@@ -19,8 +19,8 @@ datasets:
   stack: false
   geom:
     axes:
-    energy: {min: 0.1 TeV, max: 100 TeV, nbins: 30}
-    energy_true: {min: 0.05 TeV, max: 100 TeV, nbins: 40}
+      energy: {min: 0.1 TeV, max: 100 TeV, nbins: 30}
+      energy_true: {min: 0.05 TeV, max: 100 TeV, nbins: 40}
   exclusion_region:
     on_radius: 0.5 deg
     magnitude_B: 7
@@ -43,7 +43,7 @@ flux_points:
   flux: {min: 1.e-15, max: 1.e-9}
 
 light_curve:
-  energy: {min: 1. TeV, max: 100 TeV}
+  energy: {min: 0.5 TeV, max: 100 TeV}
   time_zone: -7
   n_jobs: 4
   #    time_bin_files:

v2dl5/sky_regions.py

@@ -174,7 +174,8 @@ def _read_bright_star_catalogue(self, exclusion_region_dict=None, max_wobble_dis
             return _exclusion_regions
 
         self._logger.info(
-            "Reading bright star catalogue from %s", exclusion_region_dict["star_file"]
+            f"Reading bright star catalogue from {exclusion_region_dict['star_file']}"
+            f" (max wobble distance: {max_wobble_distance})"
         )
         star_file = files("v2dl5.data").joinpath("data/" + exclusion_region_dict["star_file"])
         hip = fits.open(star_file)
@@ -189,6 +190,7 @@ def _read_bright_star_catalogue(self, exclusion_region_dict=None, max_wobble_dis
         )
 
         catalogue = catalogue[catalogue["angular_separation"] < u.Quantity(max_wobble_distance)]
+        catalogue = catalogue[[row is not None for row in catalogue]]
 
         self._logger.info(
             "Number of stars in the catalogue passing cuts on magnitude and FOV: %d", len(catalogue)