def make_edisp_map(edisp,
pointing,
geom,
exposure_map=None,
use_region_center=True):
"""Make a edisp map for a single observation
Expected axes : migra and true energy in this specific order
The name of the migra MapAxis is expected to be 'migra'
Parameters
----------
edisp : `~gammapy.irf.EnergyDispersion2D`
the 2D Energy Dispersion IRF
pointing : `~astropy.coordinates.SkyCoord`
the pointing direction
geom : `~gammapy.maps.Geom`
the map geom to be used. It provides the target geometry.
migra and true energy axes should be given in this specific order.
exposure_map : `~gammapy.maps.Map`, optional
the associated exposure map.
default is None
use_region_center: Bool
If geom is a RegionGeom, whether to just
consider the values at the region center
or the instead the average over the whole region
Returns
-------
edispmap : `~gammapy.irf.EDispMap`
the resulting EDisp map
"""
# Compute separations with pointing position
if not use_region_center:
coords, weights = geom.get_wcs_coord_and_weights()
else:
coords, weights = geom.get_coord(sparse=True), None
offset = coords.skycoord.separation(pointing)
# Compute EDisp values
data = edisp.evaluate(
offset=offset,
energy_true=coords["energy_true"],
migra=coords["migra"],
)
if not use_region_center:
data = np.average(data, axis=2, weights=weights)
# Create Map and fill relevant entries
edisp_map = Map.from_geom(geom, data=data.to_value(""), unit="")
edisp_map.normalize(axis_name="migra")
return EDispMap(edisp_map, exposure_map)
def make_edisp_map(edisp, pointing, geom, exposure_map=None):
"""Make a edisp map for a single observation
Expected axes : migra and true energy in this specific order
The name of the migra MapAxis is expected to be 'migra'
Parameters
----------
edisp : `~gammapy.irf.EnergyDispersion2D`
the 2D Energy Dispersion IRF
pointing : `~astropy.coordinates.SkyCoord`
the pointing direction
geom : `~gammapy.maps.Geom`
the map geom to be used. It provides the target geometry.
migra and true energy axes should be given in this specific order.
exposure_map : `~gammapy.maps.Map`, optional
the associated exposure map.
default is None
Returns
-------
edispmap : `~gammapy.cube.EDispMap`
the resulting EDisp map
"""
energy_axis = geom.get_axis_by_name("energy_true")
energy = energy_axis.center
migra_axis = geom.get_axis_by_name("migra")
migra = migra_axis.center
# Compute separations with pointing position
offset = geom.separation(pointing)
# Compute EDisp values
edisp_values = edisp.data.evaluate(
offset=offset,
energy_true=energy[:, np.newaxis, np.newaxis, np.newaxis],
migra=migra[:, np.newaxis, np.newaxis],
)
# Create Map and fill relevant entries
data = edisp_values.to_value("")
edispmap = Map.from_geom(geom, data=data, unit="")
return EDispMap(edispmap, exposure_map)