Add fit status fields to output container of MuonIntensityFitter

ctapipe/containers.py

@@ -984,6 +984,11 @@ class MuonEfficiencyContainer(Container):
     impact_x = Field(nan * u.m, "impact parameter x position")
     impact_y = Field(nan * u.m, "impact parameter y position")
     optical_efficiency = Field(nan, "optical efficiency muon")
+    is_valid = Field(False, "True if the fit converged successfully")
+    parameters_at_limit = Field(
+        False, "True if any bounded parameter was fitted close to a bound"
+    )
+    likelihood_value = Field(nan, "cost function value at the minimum")
 
 
 class MuonParametersContainer(Container):

ctapipe/image/muon/intensity_fitter.py

@@ -545,4 +545,7 @@ def __call__(self, tel_id, center_x, center_y, radius, image, pedestal, mask=Non
             impact_y=result["impact_parameter"] * np.sin(result["phi"]) * u.m,
             width=u.Quantity(np.rad2deg(result["ring_width"]), u.deg),
             optical_efficiency=result["optical_efficiency_muon"],
+            is_valid=minuit.valid,
+            parameters_at_limit=minuit.fmin.has_parameters_at_limit,
+            likelihood_value=minuit.fval,
         )

ctapipe/image/muon/tests/test_intensity_fit.py

@@ -89,6 +89,9 @@ def test_muon_efficiency_fit(prod5_lst):
     assert u.isclose(result.impact, impact_parameter, rtol=0.05)
     assert u.isclose(result.width, ring_width, rtol=0.05)
     assert u.isclose(result.optical_efficiency, efficiency, rtol=0.05)
+    assert result.is_valid
+    assert not result.parameters_at_limit
+    assert np.isfinite(result.likelihood_value)
 
 
 def test_scts(prod5_sst):

docs/changes/2381.feature.rst

@@ -0,0 +1 @@
+Added new fields to the ``MuonEfficiencyContainer`` - ``is_valid`` to check if fit converged successfully, ``parameters_at_limit`` to check if parameters were fitted close to a bound and ``likelihood_value`` which represents cost function value atthe minimum. These fields were added to the output of the ``MuonIntensityFitter``.