Remove GCTAEventBin::emin() and emax() methods (no interface change in bugfix release)

Jürgen Knödlseder
1 parent c4911769
Showing 5 changed files with 17 additions and 58 deletions
inst/cta/include/GCTAEventBin.hpp
inst/cta/pyext/GCTAEventBin.i
inst/cta/src/GCTAEventBin.cpp
inst/cta/src/GCTAModelIrfBackground.cpp
inst/cta/test/test_CTA.cpp
@@ -80,8 +80,6 @@ public:
     const int&     ipix(void) const;
     const int&     ieng(void) const;
     const double&  solidangle(void) const;
-    GEnergy        emin(void) const;
-    GEnergy        emax(void) const;
     const GEnergy& ewidth(void) const;
     const double&  ontime(void) const;
     const double&  weight(void) const;
@@ -187,18 +185,4 @@ void GCTAEventBin::ieng(const int&amp; ieng)
     return;
 }
  
-
-/***********************************************************************//**
- * @brief Return maximum energy of event bin
- *
- * @return Maximum energy of event bin.
- *
- * Returns the maximum energy of event bin.
- ***************************************************************************/
-inline
-GEnergy GCTAEventBin::emax(void) const
-{
-    return (emin() + ewidth());
-}
-
 #endif /* GCTAEVENTBIN_HPP */
@@ -61,8 +61,6 @@ public:
     const int&     ipix(void) const;
     const int&     ieng(void) const;
     const double&  solidangle(void) const;
-    GEnergy        emin(void) const;
-    GEnergy        emax(void) const;
     const GEnergy& ewidth(void) const;
     const double&  ontime(void) const;
     const double&  weight(void) const;
@@ -342,40 +342,6 @@ const double&amp; GCTAEventBin::solidangle(void) const
  
  
 /***********************************************************************//**
- * @brief Return minimum energy of event bin
- *
- * @return Minimum energy of event bin
- *
- * Returns minimum energy of event bin, computed using
- *
- * \f[
- *    E_{\rm min} = \frac{-\Delta E + sqrt{\Delta E^2 + 4 E^2}}{2}
- * \f]
- *
- * where
- * \f$\Delta E\f$ is the energy bin width, returned by ewidth() and
- * \f$E\f$ is the energy, returned by energy().
- ***************************************************************************/
-GEnergy GCTAEventBin::emin(void) const
-{
-    // Get energy and energy width in MeV
-    double energy = this->energy().MeV();
-    double ewidth = this->ewidth().MeV();
-
-    // Compute minimum energy in MeV
-    double arg   = ewidth * ewidth + 4.0 * energy * energy;
-    double e_min = 0.5 * (-ewidth + std::sqrt(arg));
-
-    // Set minimum energy
-    GEnergy emin;
-    emin.MeV(e_min);
-
-    // Return minimum energy
-    return emin;
-}
-
-
-/***********************************************************************//**
  * @brief Return energy width of event bin
  *
  * @return Energy width of event bin
@@ -360,11 +360,24 @@ double GCTAModelIrfBackground::eval(const GEvent&amp;       event,
         spat        = bgd(logE, inst_dir.detx(), inst_dir.dety());
     }
     else {
+        // Get pointer to event bin
         const GCTAEventBin* bin = static_cast<const GCTAEventBin*>(&event);
-        GEnergy emin  = bin->emin();
-        GEnergy emax  = emin + bin->ewidth();
-        double  norm  = 1.0 / bin->ewidth().MeV();
-        spat          = bgd.rate_ebin(inst_dir, emin, emax) * norm;
+
+        // Get energy and energy width in MeV
+        double energy = bin->energy().MeV();
+        double ewidth = bin->ewidth().MeV();
+
+        // Compute minimum energy in MeV
+        double arg   = ewidth * ewidth + 4.0 * energy * energy;
+        double e_min = 0.5 * (-ewidth + std::sqrt(arg));
+
+        // Set minimum and maximumenergy
+        GEnergy emin;
+        emin.MeV(e_min);
+        GEnergy emax = emin + bin->ewidth();
+
+        // Compute spatial component
+        spat = bgd.rate_ebin(inst_dir, emin, emax) / ewidth;
     }
     #else
     double logE = event.energy().log10TeV();
@@ -2229,8 +2229,6 @@ void TestGCTAObservation::test_event_bin(void)
     test_value(bin.ieng(), 7, "Check ieng() for filled bin");
     test_value(bin.solidangle(), 3.14, 1.0e-10,
                "Check solidangle() for filled bin");
-    test_value(bin.emin().TeV(), 1.996503, "Check emin() for filled bin");
-    test_value(bin.emax().TeV(), 2.003503,  "Check emax() for filled bin");
     test_value(bin.ewidth().GeV(), 7.0, 1.0e-10,
                "Check ewidth() for filled bin");
     test_value(bin.ontime(), 101.0, 1.0e-10, "Check ontime() for filled bin");