Skip to content

Commit a07ab02

Browse files
committed
fix ptsrc masking failure sometimes
1 parent f8ce829 commit a07ab02

File tree

1 file changed

+3
-2
lines changed

1 file changed

+3
-2
lines changed

src/masking.jl

Lines changed: 3 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -11,12 +11,13 @@ function make_mask(
1111
deg2npix(x) = round(Int, x/θpix*60)
1212
arcmin2npix(x) = round(Int, x/θpix)
1313

14-
ptsrc = .!bleed(sim_ptsrcs(Nside, num_ptsrcs), arcmin2npix(ptsrc_radius_arcmin))
14+
ptsrc = num_ptsrcs==0 ? 1 : .!bleed(sim_ptsrcs(Nside, num_ptsrcs), arcmin2npix(ptsrc_radius_arcmin))
1515
boundary = boundarymask(Nside, deg2npix(edge_padding_deg))
1616
mask_array = if apodization_deg in (false, 0)
1717
boundary .& ptsrc
1818
else
19-
cos_apod(boundary, deg2npix(apodization_deg), deg2npix(edge_rounding_deg)) .* cos_apod(ptsrc, arcmin2npix(ptsrc_radius_arcmin));
19+
apod_ptsrc = num_ptsrcs==0 ? 1 : cos_apod(ptsrc, arcmin2npix(ptsrc_radius_arcmin));
20+
cos_apod(boundary, deg2npix(apodization_deg), deg2npix(edge_rounding_deg)) .* apod_ptsrc
2021
end
2122

2223
FlatMap(Float32.(mask_array), θpix=θpix)

0 commit comments

Comments
 (0)