Commit 2d58fa4f4e28956d80580471a2c4347ed4cc77e9
1 parent
c948c52d
update documentation
Showing
1 changed file
with
3 additions
and
23 deletions
src/model/GModel.cpp
| @@ -621,9 +621,7 @@ GModelSpectral* GModel::xml_spectral(const GXmlElement& spectral) const | @@ -621,9 +621,7 @@ GModelSpectral* GModel::xml_spectral(const GXmlElement& spectral) const | ||
| 621 | * @parem[in] grad Compute also model gradients (default=false). | 621 | * @parem[in] grad Compute also model gradients (default=false). |
| 622 | * | 622 | * |
| 623 | * This method computes the spatial model component for a given true photon | 623 | * This method computes the spatial model component for a given true photon |
| 624 | - * energy and arrival time. It distinguishes the case of a point source from | ||
| 625 | - * that of a diffuse model and calls the respective methods() of the | ||
| 626 | - * response class to compute the spatial response. | 624 | + * energy and arrival time. |
| 627 | ***************************************************************************/ | 625 | ***************************************************************************/ |
| 628 | double GModel::spatial(const GEvent& event, | 626 | double GModel::spatial(const GEvent& event, |
| 629 | const GEnergy& srcEng, const GTime& srcTime, | 627 | const GEnergy& srcEng, const GTime& srcTime, |
| @@ -639,26 +637,8 @@ double GModel::spatial(const GEvent& event, | @@ -639,26 +637,8 @@ double GModel::spatial(const GEvent& event, | ||
| 639 | GPointing* pnt = obs.pointing(event.time()); | 637 | GPointing* pnt = obs.pointing(event.time()); |
| 640 | GResponse* rsp = obs.response(event.time()); | 638 | GResponse* rsp = obs.response(event.time()); |
| 641 | 639 | ||
| 642 | - // Allocate IRF value | ||
| 643 | - double irf; | ||
| 644 | - | ||
| 645 | - // If the spatial model is a point source then extract the point | ||
| 646 | - // source location from the spatial model and compute the point | ||
| 647 | - // source IRF for that location. | ||
| 648 | - if (m_spatial->isptsource()) { | ||
| 649 | - | ||
| 650 | - // Get point source location | ||
| 651 | - GSkyDir srcDir = ((GModelSpatialPtsrc*)m_spatial)->dir(); | ||
| 652 | - | ||
| 653 | - // Compute IRF | ||
| 654 | - irf = rsp->irf(event.dir(), event.energy(), event.time(), | ||
| 655 | - srcDir, srcEng, srcTime, *pnt); | ||
| 656 | - | ||
| 657 | - } // endif: Model was a point source | ||
| 658 | - | ||
| 659 | - // ... otherwise compute the diffuse instrument response function. | ||
| 660 | - else | ||
| 661 | - irf = rsp->diffrsp(event, *this, srcEng, srcTime, *pnt); | 640 | + // Get IRF value |
| 641 | + double irf = rsp->irf(event, *this, srcEng, srcTime, *pnt); | ||
| 662 | 642 | ||
| 663 | // Compute source model | 643 | // Compute source model |
| 664 | double source = 1.0; | 644 | double source = 1.0; |