Skip to content

Commit

Permalink
use angular_error_floor as sigma
Browse files Browse the repository at this point in the history
  • Loading branch information
G-Sommani committed Oct 9, 2024
1 parent e8bb97d commit 2b4f5b5
Showing 1 changed file with 7 additions and 6 deletions.
13 changes: 7 additions & 6 deletions skyreader/plot/plot.py
Original file line number Diff line number Diff line change
Expand Up @@ -421,9 +421,9 @@ def create_plot_zoomed(
circular=False,
circular_err50=0.2,
circular_err90=0.7,
angular_error_floor=False,
llh_map=True,
neutrinofloor_sigma=0.2
angular_error_floor=None, # if not None, sigma of the
# gaussian to convolute the map with in deg.
llh_map=True
):
"""Uses healpy to plot a map."""

Expand Down Expand Up @@ -541,11 +541,12 @@ def bounding_box(ra, dec, theta, phi):
min_map = np.nanmin(equatorial_map)
equatorial_map[np.isnan(equatorial_map)] = min_map

if angular_error_floor:
# convolute with a gaussian with 0.2 deg as sigma
if angular_error_floor is not None:
# convolute with a gaussian. angular_error_floor is the
# sigma in deg.
equatorial_map = healpy.smoothing(
equatorial_map,
sigma=np.deg2rad(neutrinofloor_sigma),
sigma=np.deg2rad(angular_error_floor),
)

# normalize map
Expand Down

0 comments on commit 2b4f5b5

Please sign in to comment.