Add a StereoDispCombiner

docs/changes/2731.feature.rst

@@ -0,0 +1,11 @@
+Added a new ``StereoDispCombiner`` implementing DISP-based stereo direction
+reconstruction.
+The algorithm is adapted and generalized from the implementations in
+*magic-cta-pipe* and *EventDisplay* and resolves the head–tail ambiguity by
+evaluating all DISP sign combinations and combining multiple telescope
+combinations via weighted means.
+
+The reconstruction supports optional pre-selection of telescopes using
+``n_best_tels`` and rejection of nearly parallel multiplicity-2 events via
+``min_ang_diff``.
+Only geometry reconstruction is supported.

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",
 ]