EnergyDependentTablePSF

class gammapy.irf.EnergyDependentTablePSF(energy_axis_true, rad_axis, exposure=None, psf_value=None, interp_kwargs=None)

Bases: object

Energy-dependent radially-symmetric table PSF (gtpsf format).

TODO: add references and explanations.

Parameters

energy_axis_trueMapAxis

Energy axis

rad_axisMapAxis

Offset angle wrt source position axis

exposureQuantity

Exposure (1-dim)

psf_valueQuantity

PSF (2-dim with axes: psf[energy_index, offset_index]

interp_kwargsdict

Interpolation keyword arguments pass to ScaledRegularGridInterpolator.

Attributes Summary

energy_axis_true
rad_axis

Methods Summary

containment(energy, rad_max)	Compute containment of the PSF.
containment_radius(energy[, fraction])	Containment radius.
evaluate([energy, rad, method])	Evaluate the PSF at a given energy and offset
from_hdulist(hdu_list)	Create EnergyDependentTablePSF from gtpsf format HDU list.
info()	Print basic info
plot_containment_vs_energy([ax, fractions])	Plot containment versus energy.
plot_exposure_vs_energy()	Plot exposure versus energy.
plot_psf_vs_rad([energies, ax])	Plot PSF vs radius.
read(filename)	Create EnergyDependentTablePSF from gtpsf-format FITS file.
stack(psf)	Stack two EnergyDependentTablePSF objects.s
table_psf_at_energy(energy[, method])	Create TablePSF at one given energy.
table_psf_in_energy_range(energy_range[, …])	Average PSF in a given energy band.
to_hdulist()	Convert to FITS HDU list format.
write(filename, *args, **kwargs)	Write to FITS file.

Attributes Documentation

energy_axis_true

rad_axis

Methods Documentation

containment(energy, rad_max)

Compute containment of the PSF.

Parameters

energyQuantity

Energy

rad_maxAngle

Maximum offset angle.

Returns

fractionarray_like

Containment fraction (in range 0 .. 1)

containment_radius(energy, fraction=0.68)

Containment radius.

Parameters

energyQuantity

Energy

fractionfloat

Containment fraction.

Returns

radQuantity

Containment radius in deg

evaluate(energy=None, rad=None, method='linear')

Evaluate the PSF at a given energy and offset

Parameters

energyQuantity

Energy value

radAngle

Offset wrt source position

method{“linear”, “nearest”}

Linear or nearest neighbour interpolation.

Returns

valuesQuantity

Interpolated value

classmethod from_hdulist(hdu_list)

Create EnergyDependentTablePSF from gtpsf format HDU list.

Parameters

hdu_listHDUList

HDU list with THETA and PSF extensions.

info()

Print basic info

plot_containment_vs_energy(ax=None, fractions=[0.68, 0.8, 0.95], **kwargs)

Plot containment versus energy.

plot_exposure_vs_energy()

Plot exposure versus energy.

plot_psf_vs_rad(energies=None, ax=None, **kwargs)

Plot PSF vs radius.

Parameters

energyQuantity

Energies where to plot the PSF.

**kwargsdict

Keyword arguments pass to plot.

classmethod read(filename)

Create EnergyDependentTablePSF from gtpsf-format FITS file.

Parameters

filenamestr

File name

stack(psf)

Stack two EnergyDependentTablePSF objects.s

Parameters

psfEnergyDependentTablePSF

PSF to stack.

Returns

stacked_psfEnergyDependentTablePSF

Stacked PSF.

table_psf_at_energy(energy, method='linear', **kwargs)

Create TablePSF at one given energy.

Parameters

energyQuantity

Energy

method{“linear”, “nearest”}

Linear or nearest neighbour interpolation.

Returns

psfTablePSF

Table PSF

table_psf_in_energy_range(energy_range, spectrum=None, n_bins=11, **kwargs)

Average PSF in a given energy band.

Expected counts in sub energy bands given the given exposure and spectrum are used as weights.

Parameters

energy_rangeQuantity

Energy band

spectrumSpectralModel

Spectral model used for weighting the PSF. Default is a power law with index=2.

n_binsint

Number of energy points in the energy band, used to compute the weigthed PSF.

Returns

psfTablePSF

Table PSF

to_hdulist()

Convert to FITS HDU list format.

Returns

hdu_listHDUList

PSF in HDU list format.

write(filename, *args, **kwargs)

Write to FITS file.

Calls writeto, forwarding all arguments.