EnergyDependentMultiGaussPSF¶
-
class
gammapy.irf.
EnergyDependentMultiGaussPSF
(energy_axis_true, offset_axis, sigmas, norms, energy_thresh_lo='0.1 TeV', energy_thresh_hi='100 TeV')[source]¶ Bases:
object
Triple Gauss analytical PSF depending on energy and theta.
To evaluate the PSF call the
to_energy_dependent_table_psf
orpsf_at_energy_and_theta
methods.- Parameters
- energy_axis_true
MapAxis
True energy axis
- offset_axis
MapAxis
Offset axis.
- sigmaslist of ‘numpy.ndarray’
Triple Gauss sigma parameters, where every entry is a two dimensional ‘numpy.ndarray’ containing the sigma value for every given energy and theta.
- normslist of ‘numpy.ndarray’
Triple Gauss norm parameters, where every entry is a two dimensional ‘numpy.ndarray’ containing the norm value for every given energy and theta. Norm corresponds to the value of the Gaussian at theta = 0.
- energy_thresh_lo
Quantity
Lower save energy threshold of the psf.
- energy_thresh_hi
Quantity
Upper save energy threshold of the psf.
- energy_axis_true
Examples
Plot R68 of the PSF vs. theta and energy:
import matplotlib.pyplot as plt from gammapy.irf import EnergyDependentMultiGaussPSF filename = '$GAMMAPY_DATA/cta-1dc/caldb/data/cta/1dc/bcf/South_z20_50h/irf_file.fits' psf = EnergyDependentMultiGaussPSF.read(filename, hdu='POINT SPREAD FUNCTION') psf.plot_containment(0.68) plt.show()
Attributes Summary
Methods Summary
containment_radius
(energy, theta[, fraction])Compute containment for all energy and theta values
from_table_hdu
(hdu)Create
EnergyDependentMultiGaussPSF
from HDU list.info
([fractions, energies, thetas])Print PSF summary info.
peek
([figsize])Quick-look summary plots.
plot_containment
([fraction, ax, add_cbar])Plot containment image with energy and theta axes.
plot_containment_vs_energy
([fractions, …])Plot containment fraction as a function of energy.
psf_at_energy_and_theta
(energy, theta)Get
MultiGauss2D
model for given energy and theta.read
(filename[, hdu])Create
EnergyDependentMultiGaussPSF
from FITS file.to_energy_dependent_table_psf
([theta, rad, …])Convert triple Gaussian PSF ot table PSF.
Convert psf table data to FITS hdu list.
to_psf3d
([rad])Create a PSF3D from an analytical PSF.
write
(filename, *args, **kwargs)Write PSF to FITS file.
Attributes Documentation
-
energy_axis_true
¶
-
offset_axis
¶
-
tag
= 'psf_3gauss'¶
Methods Documentation
-
containment_radius
(energy, theta, fraction=0.68)[source]¶ Compute containment for all energy and theta values
-
classmethod
from_table_hdu
(hdu)[source]¶ Create
EnergyDependentMultiGaussPSF
from HDU list.- Parameters
- hdu
BinTableHDU
HDU
- hdu
-
info
(fractions=[0.68, 0.95], energies=<Quantity [ 1., 10.] TeV>, thetas=<Quantity [0.] deg>)[source]¶ Print PSF summary info.
The containment radius for given fraction, energies and thetas is computed and printed on the command line.
- Parameters
- fractionslist
Containment fraction to compute containment radius for.
- energies
Quantity
Energies to compute containment radius for.
- thetas
Quantity
Thetas to compute containment radius for.
- Returns
- ssstring
Formatted string containing the summary info.
-
plot_containment
(fraction=0.68, ax=None, add_cbar=True, **kwargs)[source]¶ Plot containment image with energy and theta axes.
- Parameters
- fractionfloat
Containment fraction between 0 and 1.
- add_cbarbool
Add a colorbar
-
plot_containment_vs_energy
(fractions=[0.68, 0.95], thetas=<Angle [0., 1.] deg>, ax=None, **kwargs)[source]¶ Plot containment fraction as a function of energy.
-
psf_at_energy_and_theta
(energy, theta)[source]¶ Get
MultiGauss2D
model for given energy and theta.No interpolation is used.
- Parameters
- energy
Quantity
Energy at which a PSF is requested.
- theta
Angle
Offset angle at which a PSF is requested.
- energy
- Returns
- psf
MultiGauss2D
Multigauss PSF object.
- psf
-
classmethod
read
(filename, hdu='PSF_2D_GAUSS')[source]¶ Create
EnergyDependentMultiGaussPSF
from FITS file.- Parameters
- filenamestr
File name
-
to_energy_dependent_table_psf
(theta=None, rad=None, exposure=None)[source]¶ Convert triple Gaussian PSF ot table PSF.
- Parameters
- Returns
- tabe_psf
EnergyDependentTablePSF
Instance of
EnergyDependentTablePSF
.
- tabe_psf
-
to_hdulist
()[source]¶ Convert psf table data to FITS hdu list.
- Returns
- hdu_list
HDUList
PSF in HDU list format.
- hdu_list