Commit 7d33bb30579fc9add561188a221eb7d65fe665af
1 parent
2336cb83
Increase elliptical gaussian cut-off from 2.5 to 3.0 as the results for 2.5 stil…
…l seem to be biased in the ellipse size.
Showing
1 changed file
with
5 additions
and
5 deletions
src/model/GModelSpatialEllipticalGauss.cpp
| ... | ... | @@ -37,7 +37,7 @@ |
| 37 | 37 | |
| 38 | 38 | /* __ Constants __________________________________________________________ */ |
| 39 | 39 | namespace { |
| 40 | - const double c_theta_max = 2.5; //!< semiaxis multiplied for theta_max() | |
| 40 | + const double c_theta_max = 3.0; //!< semiaxis multiplied for theta_max() | |
| 41 | 41 | const double c_max_exponent = 0.5 * c_theta_max * c_theta_max; |
| 42 | 42 | const double c_fraction = 1.0 - std::exp(-c_max_exponent); |
| 43 | 43 | } |
| ... | ... | @@ -290,7 +290,7 @@ GModelSpatialEllipticalGauss* GModelSpatialEllipticalGauss::clone(void) const |
| 290 | 290 | * |
| 291 | 291 | * @warning |
| 292 | 292 | * For numerical reasons the elliptical Gaussian will be truncated for |
| 293 | - * \f$\theta\f$ angles that correspond to 2.5 times the effective ellipse | |
| 293 | + * \f$\theta\f$ angles that correspond to 3.0 times the effective ellipse | |
| 294 | 294 | * radius. |
| 295 | 295 | ***************************************************************************/ |
| 296 | 296 | double GModelSpatialEllipticalGauss::eval(const double& theta, |
| ... | ... | @@ -395,7 +395,7 @@ double GModelSpatialEllipticalGauss::eval_gradients(const double& theta, |
| 395 | 395 | * |
| 396 | 396 | * @warning |
| 397 | 397 | * For numerical reasons the elliptical Gaussian will be truncated for |
| 398 | - * \f$\theta\f$ angles that correspond to 2.5 times the effective ellipse | |
| 398 | + * \f$\theta\f$ angles that correspond to 3.0 times the effective ellipse | |
| 399 | 399 | * radius. |
| 400 | 400 | ***************************************************************************/ |
| 401 | 401 | GSkyDir GModelSpatialEllipticalGauss::mc(const GEnergy& energy, |
| ... | ... | @@ -468,10 +468,10 @@ bool GModelSpatialEllipticalGauss::contains(const GSkyDir& dir, |
| 468 | 468 | * |
| 469 | 469 | * @return Returns maximum model radius. |
| 470 | 470 | * |
| 471 | - * Returns the maximum of 2.5 semimajor() and 2.5 semiminor() as | |
| 471 | + * Returns the maximum of 3.0 semimajor() and 3.0 semiminor() as | |
| 472 | 472 | * approximate edge of the Gaussian. This limit is of course arbitrary, but |
| 473 | 473 | * allows to limit the integration region for response computation. The value |
| 474 | - * of 2.5 has been determined by experiment (#1561). | |
| 474 | + * of 3.0 has been determined by experiment (#1561). | |
| 475 | 475 | ***************************************************************************/ |
| 476 | 476 | double GModelSpatialEllipticalGauss::theta_max(void) const |
| 477 | 477 | { | ... | ... |