def test_non_circular_regions(region, nreg):
pointing = SkyCoord(0.0, 0.0, unit="deg")
finder = ReflectedRegionsFinder(center=pointing,
region=region,
min_distance_input="0 deg")
finder.run()
regions = finder.reflected_regions
assert len(regions) == nreg
def test_bad_on_region(exclusion_mask, on_region):
pointing = SkyCoord(83.63, 22.01, unit="deg", frame="icrs")
finder = ReflectedRegionsFinder(
center=pointing,
region=on_region,
exclusion_mask=exclusion_mask,
min_distance_input="0 deg",
)
finder.run()
regions = finder.reflected_regions
assert len(regions) == 0
# try plotting
with mpl_plot_check():
finder.plot()
def test_find_reflected_regions(exclusion_mask, on_region, pointing_pos, nreg1,
reg3_ra, nreg2, nreg3):
pointing = pointing_pos
finder = ReflectedRegionsFinder(
center=pointing,
region=on_region,
exclusion_mask=exclusion_mask,
min_distance_input="0 deg",
)
finder.run()
regions = finder.reflected_regions
assert len(regions) == nreg1
assert_quantity_allclose(regions[3].center.icrs.ra, reg3_ra, rtol=1e-2)
# Test without exclusion
finder.exclusion_mask = None
finder.run()
regions = finder.reflected_regions
assert len(regions) == nreg2
# Test with too small exclusion
small_mask = exclusion_mask.cutout(pointing, Angle("0.1 deg"))
finder.exclusion_mask = small_mask
finder.run()
regions = finder.reflected_regions
assert len(regions) == nreg3
# Test with maximum number of regions
finder.max_region_number = 5
finder.run()
regions = finder.reflected_regions
assert len(regions) == 5
# Test with an other type of region
on_ellipse_annulus = EllipseAnnulusSkyRegion(
center=on_region.center.galactic,
inner_width=0.1 * u.deg,
outer_width=0.2 * u.deg,
inner_height=0.3 * u.deg,
outer_height=0.6 * u.deg,
angle=130 * u.deg,
)
finder.region = on_ellipse_annulus
finder.reference_map = None
finder.run()
regions = finder.reflected_regions
assert len(regions) == 5
binsz=0.01,
skydir=(83.6, 23.0))
coords = exclusion_mask.geom.get_coord().skycoord
mask = (Angle("23d") < coords.dec) & (coords.dec < Angle("24d"))
exclusion_mask.data = np.invert(mask)
pos = SkyCoord(83.633, 22.014, unit="deg")
radius = Angle(0.3, "deg")
on_region = CircleSkyRegion(pos, radius)
center = SkyCoord(83.633, 24, unit="deg")
# One can impose a minimal distance between ON region and first reflected regions
finder = ReflectedRegionsFinder(
region=on_region,
center=center,
exclusion_mask=exclusion_mask,
min_distance_input="0.2 rad",
)
finder.run()
fig1 = plt.figure(1)
finder.plot(fig=fig1)
# One can impose a minimal distance between two adjacent regions
finder = ReflectedRegionsFinder(
region=on_region,
center=center,
exclusion_mask=exclusion_mask,
min_distance_input="0.2 rad",
min_distance="0.1 rad",
)