Modify wrapper for doubling time to return the doubling coordinates a…

…s an astropy.Time instead of an astropy.units.Quanity Signed-off-by: cgalelli <[email protected]>
cgalelli · Dec 4, 2023 · a12b15d · a12b15d
1 parent 20f25bf
commit a12b15d
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/gammapy/estimators/utils.py b/gammapy/estimators/utils.py
@@ -439,10 +439,12 @@ def compute_lightcurve_doublingtime(lightcurve, flux_quantity="flux"):
             energies[1:],
             doubling_dict["doubling"],
             doubling_dict["doubling_err"],
-            doubling_dict["doubling_coord"],
+            lightcurve.geom.axes["time"].reference_time
+            + doubling_dict["doubling_coord"],
             doubling_dict["halving"],
             doubling_dict["halving_err"],
-            doubling_dict["halving_coord"],
+            lightcurve.geom.axes["time"].reference_time
+            + doubling_dict["halving_coord"],
         ],
         names=(
             "min_energy",