Kernel for IRF offest angle integration of the diffuse source model. More...
#include <GCTAResponse_helpers.hpp>
Public Member Functions | |
| cta_irf_diffuse_kern_theta (const GCTAResponseIrf *rsp, const GModelSpatial *model, const GMatrix *rot, const double &theta, const double &phi, const double &zenith, const double &azimuth, const GEnergy &srcEng, const GTime &srcTime, const double &srcLogEng, const GEnergy &obsEng, const double &eta, const int &min_iter, const int &max_iter, const double &resolution) | |
| double | eval (const double &theta) |
| Kernel for IRF offest angle integration of the diffuse source model. | |
 Public Member Functions inherited from GFunction | |
| GFunction (void) | |
| Void constructor. | |
| GFunction (const GFunction &function) | |
| Copy constructor. | |
| virtual | ~GFunction (void) |
| Destructor. | |
| GFunction & | operator= (const GFunction &function) |
| Assignment operator. | |
Protected Attributes | |
| const GCTAResponseIrf * | m_rsp |
| CTA response. | |
| const GModelSpatial * | m_model |
| Spatial model. | |
| const GMatrix * | m_rot |
| Rotation matrix. | |
| double | m_theta |
| Photon offset angle. | |
| double | m_phi |
| Photon azimuth angle. | |
| double | m_zenith |
| Pointing zenith angle. | |
| double | m_azimuth |
| Pointing azimuth angle. | |
| GEnergy | m_srcEng |
| True photon energy. | |
| GTime | m_srcTime |
| True photon arrival time. | |
| double | m_srcLogEng |
| True photon log energy. | |
| GEnergy | m_obsEng |
| Measured event energy. | |
| double | m_sin_eta |
| Sine of angular distance between. | |
| double | m_cos_eta |
| Cosine of angular distance between. | |
| int | m_min_iter |
| int | m_max_iter |
| double | m_resolution |
Additional Inherited Members | |
| Protected Member Functions inherited from GFunction | |
| void | init_members (void) |
| Initialise class members. | |
| void | copy_members (const GFunction &function) |
| Copy class members. | |
| void | free_members (void) |
| Delete class members. | |
Detailed Description
Kernel for IRF offest angle integration of the diffuse source model.
This class implements the integration kernel \(K(\theta)\) for the integration
\[ \int_{0}^{\theta_{\rm max}} K(\theta | E, t) d\theta \]
of diffuse models. The eval() method computes
\[ K(\theta | E, t) = \sin \theta \times PSF(\theta) \int_{0}^{2\pi} S_{\rm p}(\theta, \phi | E, t) \, Aeff(\theta, \phi) \, Edisp(\theta, \phi) d\phi \]
where
- \(S_{\rm p}(\theta, \phi | E, t)\) is the diffuse model,
- \(PSF(\theta)\) is the azimuthally symmetric Point Spread Function,
- \(Aeff(\theta, \phi)\) is the effective area,
- \(Edisp(\theta, \phi)\) is the energy dispersion,
- \(\theta\) is the distance from the PSF centre, and
- \(\phi\) is the azimuth angle.
Definition at line 774 of file GCTAResponse_helpers.hpp.
Constructor & Destructor Documentation
◆ cta_irf_diffuse_kern_theta()
|
inline |
Definition at line 776 of file GCTAResponse_helpers.hpp.
Member Function Documentation
◆ eval()
|
virtual |
Kernel for IRF offest angle integration of the diffuse source model.
- Parameters
-
[in] theta Offset angle with respect to observed photon direction (radians).
Computes
\[ K(\theta | E, t) = \sin \theta \times PSF(\theta) \int_{0}^{2\pi} S_{\rm p}(\theta, \phi | E, t) \, Aeff(\theta, \phi) \, Edisp(\theta, \phi) d\phi \]
The PSF is assumed to be azimuthally symmetric, hence the PSF is computed outside the azimuthal integration.
Note that the integration is only performed for \(\theta>0\). Otherwise zero is returned.
Implements GFunction.
Definition at line 1288 of file GCTAResponse_helpers.cpp.
References GIntegral::fixed_iter(), gammalib::is_infinite(), gammalib::is_notanumber(), gammalib::iter_phi(), m_azimuth, m_cos_eta, m_max_iter, m_min_iter, m_model, m_obsEng, m_phi, m_resolution, m_rot, m_rsp, m_sin_eta, m_srcEng, m_srcLogEng, m_srcTime, m_theta, m_zenith, GIntegral::print(), GCTAResponseIrf::psf(), gammalib::rad2deg, GIntegral::romberg(), and gammalib::twopi.
Member Data Documentation
◆ m_azimuth
|
protected |
Pointing azimuth angle.
Definition at line 815 of file GCTAResponse_helpers.hpp.
Referenced by eval().
◆ m_cos_eta
|
protected |
Cosine of angular distance between.
Definition at line 823 of file GCTAResponse_helpers.hpp.
Referenced by eval().
◆ m_max_iter
|
protected |
Definition at line 827 of file GCTAResponse_helpers.hpp.
Referenced by eval().
◆ m_min_iter
|
protected |
Definition at line 826 of file GCTAResponse_helpers.hpp.
Referenced by eval().
◆ m_model
|
protected |
◆ m_obsEng
|
protected |
Measured event energy.
Definition at line 819 of file GCTAResponse_helpers.hpp.
Referenced by eval().
◆ m_phi
|
protected |
Photon azimuth angle.
Definition at line 813 of file GCTAResponse_helpers.hpp.
Referenced by eval().
◆ m_resolution
|
protected |
Definition at line 828 of file GCTAResponse_helpers.hpp.
Referenced by eval().
◆ m_rot
|
protected |
◆ m_rsp
|
protected |
◆ m_sin_eta
|
protected |
Sine of angular distance between.
Definition at line 820 of file GCTAResponse_helpers.hpp.
Referenced by eval().
◆ m_srcEng
|
protected |
◆ m_srcLogEng
|
protected |
True photon log energy.
Definition at line 818 of file GCTAResponse_helpers.hpp.
Referenced by eval().
◆ m_srcTime
|
protected |
True photon arrival time.
Definition at line 817 of file GCTAResponse_helpers.hpp.
Referenced by eval().
◆ m_theta
|
protected |
◆ m_zenith
|
protected |
Pointing zenith angle.
Definition at line 814 of file GCTAResponse_helpers.hpp.
Referenced by eval().
The documentation for this class was generated from the following files:
- inst/cta/src/GCTAResponse_helpers.hpp
- inst/cta/src/GCTAResponse_helpers.cpp
