PSFMap

class gammapy.cube.PSFMap(psf_map, exposure_map=None)[source]

Bases: object

Class containing the Map of PSFs and allowing to interact with it.

Parameters:
psf_map : Map

the input PSF Map. Should be a Map with 2 non spatial axes. rad and true energy axes should be given in this specific order.

exposure_map : Map

Associated exposure map. Needs to have a consistent map geometry.

Examples

import numpy as np
from astropy import units as u
from astropy.coordinates import SkyCoord
from gammapy.maps import Map, WcsGeom, MapAxis
from gammapy.irf import EnergyDependentMultiGaussPSF, EffectiveAreaTable2D
from gammapy.cube import make_psf_map, PSFMap, make_map_exposure_true_energy

# Define energy axis. Note that the name is fixed.
energy_axis = MapAxis.from_edges(np.logspace(-1., 1., 4), unit='TeV', name='energy')
# Define rad axis. Again note the axis name
rads = np.linspace(0., 0.5, 100) * u.deg
rad_axis = MapAxis.from_edges(rads, unit='deg', name='theta')

# Define parameters
pointing = SkyCoord(0., 0., unit='deg')
max_offset = 4 * u.deg

# Create WcsGeom
geom = WcsGeom.create(binsz=0.25*u.deg, width=10*u.deg, skydir=pointing, axes=[rad_axis, energy_axis])

# Extract EnergyDependentTablePSF from CTA 1DC IRF
filename = '$GAMMAPY_DATA/cta-1dc/caldb/data/cta/1dc/bcf/South_z20_50h/irf_file.fits'
psf = EnergyDependentMultiGaussPSF.read(filename, hdu='POINT SPREAD FUNCTION')
psf3d = psf.to_psf3d(rads)
aeff2d = EffectiveAreaTable2D.read(filename, hdu='EFFECTIVE AREA')

# Create the exposure map
exposure_geom = geom.to_image().to_cube([energy_axis])
exposure_map = make_map_exposure_true_energy(pointing, "1 h", aeff2d, exposure_geom)

# create the PSFMap for the specified pointing
psf_map = make_psf_map(psf3d, pointing, geom, max_offset, exposure_map)

# Get an EnergyDependentTablePSF at any position in the image
psf_table = psf_map.get_energy_dependent_table_psf(SkyCoord(2., 2.5, unit='deg'))

# Write map to disk
psf_map.write('psf_map.fits')

Methods Summary

containment_radius_map(energy[, fraction]) Containment radius map.
from_hdulist(hdulist[, psf_hdu, …]) Convert to HDUList.
get_energy_dependent_table_psf(position) Get energy-dependent PSF at a given position.
get_psf_kernel(position, geom[, max_radius, …]) Returns a PSF kernel at the given position.
read(filename, **kwargs) Read a psf_map from file and create a PSFMap object
stack(other) Stack PSFMap with another one.
to_hdulist([psf_hdu, psf_hdubands, …]) Convert to HDUList.
write(filename[, overwrite]) Write to fits

Methods Documentation

containment_radius_map(energy, fraction=0.68)[source]

Containment radius map.

Parameters:
energy : Quantity

Scalar energy at which to compute the containment radius

fraction : float

the containment fraction (range: 0 to 1)

Returns:
containment_radius_map : Map

Containment radius map

classmethod from_hdulist(hdulist, psf_hdu='PSFMAP', psf_hdubands='BANDSPSF', exposure_hdu='EXPMAP', exposure_hdubands='BANDSEXP')[source]

Convert to HDUList.

Parameters:
psf_hdu : str

Name or index of the HDU with the psf_map data.

psf_hdubands : str

Name or index of the HDU with the psf_map BANDS table.

exposure_hdu : str

Name or index of the HDU with the exposure_map data.

exposure_hdubands : str

Name or index of the HDU with the exposure_map BANDS table.

get_energy_dependent_table_psf(position)[source]

Get energy-dependent PSF at a given position.

Parameters:
position : SkyCoord

the target position. Should be a single coordinates

Returns:
psf_table : EnergyDependentTablePSF

the table PSF

get_psf_kernel(position, geom, max_radius=None, factor=4)[source]

Returns a PSF kernel at the given position.

The PSF is returned in the form a WcsNDMap defined by the input MapGeom.

Parameters:
position : SkyCoord

the target position. Should be a single coordinate

geom : MapGeom

the target geometry to use

max_radius : Angle

maximum angular size of the kernel map

factor : int

oversampling factor to compute the PSF

Returns:
kernel : PSFKernel

the resulting kernel

classmethod read(filename, **kwargs)[source]

Read a psf_map from file and create a PSFMap object

stack(other)[source]

Stack PSFMap with another one.

The other PSFMap is projected on the current PSFMap geometry.

Parameters:
other : PSFMap

the psfmap to be stacked with this one.

Returns:
new : PSFMap

the stacked psfmap

to_hdulist(psf_hdu='PSFMAP', psf_hdubands='BANDSPSF', exposure_hdu='EXPMAP', exposure_hdubands='BANDSEXP')[source]

Convert to HDUList.

Parameters:
psf_hdu : str

Name or index of the HDU with the psf_map data.

psf_hdubands : str

Name or index of the HDU with the psf_map BANDS table.

exposure_hdu : str

Name or index of the HDU with the exposure_map data.

exposure_hdubands : str

Name or index of the HDU with the exposure_map BANDS table.

Returns:
hdu_list : HDUList
write(filename, overwrite=False, **kwargs)[source]

Write to fits