def make_map_background_irf_with_symmetry(fpi, symmetry="constant"):
axis = MapAxis.from_edges([0.1, 1, 10], name="energy", unit="TeV", interp="log")
return make_map_background_irf(
pointing=fpi,
ontime="42 s",
bkg=bkg_3d_custom(symmetry),
geom=WcsGeom.create(npix=(3, 3), binsz=4, axes=[axis], skydir=fpi.radec),
)
def test_map_background_2d(bkg_2d):
axis = MapAxis.from_edges([0.1, 1, 10], name="energy", unit="TeV", interp="log")
skydir = SkyCoord("0d", "0d", frame="galactic")
geom = WcsGeom.create(npix=(3, 3), binsz=4, axes=[axis], skydir=skydir)
bkg = make_map_background_irf(
pointing=skydir,
ontime="42 s",
bkg=bkg_2d,
geom=geom,
)
assert_allclose(bkg.data[:, 1, 1], [1.807479, 0.183212], rtol=1e-5)
def test_make_map_background_irf(bkg_3d, pars, fixed_pointing_info):
m = make_map_background_irf(
pointing=fixed_pointing_info,
ontime="42 s",
bkg=bkg_3d,
geom=geom(
map_type=pars["map_type"],
ebounds=pars["ebounds"],
skydir=fixed_pointing_info.radec,
),
oversampling=10,
)
assert m.data.shape == pars["shape"]
assert m.unit == ""
assert_allclose(m.data.sum(), pars["sum"], rtol=1e-5)
