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)