Note

You are not reading the most up to date version of Gammapy documentation.
Access the latest stable version v1.3 or the list of Gammapy releases.

IRF Theory

Modeling the expected number of detected events

To model the expected number of events a gamma-ray source should produce on a detector one has to model its effect using an instrument responce function (IRF). In general, such a function gives the probability to detect a photon emitted from true position ptrue on the sky and true energy Etrue at reconstructed position p and energy E and the effective collection area of the detector at position ptrue on the sky and true energy Etrue.

We can write the expected number of detected events N(p,E):

N(p,E)dpdE=tobsEtruedEtrueptruedptrueR(p,E|ptrue,Etrue)×Φ(ptrue,Etrue)

where:

  • R(p,E|ptrue,Etrue) is the instrument response (unit: m2TeV1)

  • Φ(ptrue,Etrue) is the sky flux model (unit: m2s1TeV1sr1)

  • tobs is the observation time: (unit: s)

The Instrument Response Functions

Most of the time, in high-level gamma-ray data (DL3), we assume that the instrument response can be simplified as the product of three independent functions:

R(p,E|ptrue,Etrue)=Aeff(ptrue,Etrue)×PSF(p|ptrue,Etrue)×Edisp(E|ptrue,Etrue),

where:

  • Aeff(ptrue,Etrue) is the effective collection area of the detector (unit: m2). It is the product of the detector collection area times its detection efficiency at true energy Etrue and position ptrue.

  • PSF(p|ptrue,Etrue) is the point spread function (unit: sr1). It gives the probability of measuring a direction p when the true direction is ptrue and the true energy is Etrue. Gamma-ray instruments consider the probability density of the angular separation between true and reconstructed directions δp=ptruep, i.e. PSF(δp|ptrue,Etrue).

  • Edisp(E|ptrue,Etrue) is the energy dispersion (unit: TeV1). It gives the probability to reconstruct the photon at energy E when the true energy is Etrue and the true position ptrue. Gamma-ray instruments consider the probability density of the migration μ=EEtrue, i.e. Edisp(μ|ptrue,Etrue).

The implicit assumption here is that energy dispersion and PSF are completely independent. This is not totally valid in some situations.

These functions are obtained through Monte-Carlo simulations of gamma-ray showers for different observing conditions, e.g. detector configuration, zenith angle of the pointing position, detector state and different event reconstruction and selection schemes. In the DL3 format, the IRF are distributed for each observing run.

Further details on individuals responses and how they are implemented in gammapy are given in Effective area, Energy Dispersion and Point Spread Function.