Add a StereoDispCombiner

docs/changes/2731.feature.rst

@@ -0,0 +1,9 @@
+Add a ``StereoDispCombiner`` implementing DISP-based stereo direction reconstruction.
+
+Adapted from the implementation in magic-cta-pipe and EventDisplay, it combines
+mono DISP reconstruction by:
+- get combinations of 2 telescope for every array event
+- calculating the minimum distance for the 4 possible SIGN pairs for each combination
+- DISP sign score can be optionally used to weight the distances
+- calculating the weighted average for the minimum distance SIGN pair (fov lon/lat) per combination
+- calculating the weighted average of all combinations for an array event with the summed weights

src/ctapipe/image/statistics.py

@@ -104,7 +104,7 @@ def descriptive_statistics(
     )
 
 
-@njit
+@njit(cache=not CTAPIPE_DISABLE_NUMBA_CACHE)
 def n_largest(n, array):
     """return the n largest values of an array"""
     return nlargest(n, array)

src/ctapipe/reco/__init__.py

@@ -17,7 +17,7 @@
     EnergyRegressor,
     ParticleClassifier,
 )
-from .stereo_combination import StereoCombiner, StereoMeanCombiner
+from .stereo_combination import StereoCombiner, StereoDispCombiner, StereoMeanCombiner
 
 __all__ = [
     "Reconstructor",
@@ -32,5 +32,6 @@
     "DispReconstructor",
     "StereoCombiner",
     "StereoMeanCombiner",
+    "StereoDispCombiner",
     "CrossValidator",
 ]