From c747db82dd2300cbad3d56831c888f1c48e97586 Mon Sep 17 00:00:00 2001
From: Benjamin Huth <benjamin.huth@cern.ch>
Date: Tue, 10 Dec 2024 13:50:28 +0100
Subject: [PATCH] try external surfaces for GSF with standard navigation

---
 .../Acts/TrackFitting/GaussianSumFitter.hpp        | 14 ++++++++++++--
 1 file changed, 12 insertions(+), 2 deletions(-)

diff --git a/Core/include/Acts/TrackFitting/GaussianSumFitter.hpp b/Core/include/Acts/TrackFitting/GaussianSumFitter.hpp
index 4efcd82e4c4..5ab9b688a79 100644
--- a/Core/include/Acts/TrackFitting/GaussianSumFitter.hpp
+++ b/Core/include/Acts/TrackFitting/GaussianSumFitter.hpp
@@ -150,7 +150,7 @@ struct GaussianSumFitter {
     static_assert(std::is_same_v<Navigator, typename propagator_t::Navigator>);
 
     // Initialize the forward propagation with the DirectNavigator
-    auto fwdPropInitializer = [this](const auto& opts) {
+    auto fwdPropInitializer = [&](const auto& opts) {
       using Actors =
           ActorList<GsfActor, EndOfWorldReached, NavigationBreakAborter>;
       using PropagatorOptions = typename propagator_t::template Options<Actors>;
@@ -159,6 +159,11 @@ struct GaussianSumFitter {
 
       propOptions.setPlainOptions(opts.propagatorPlainOptions);
 
+      for (auto it = begin; it != end; ++it) {
+        propOptions.navigation.insertExternalSurface(
+            options.extensions.surfaceAccessor(SourceLink{*it})->geometryId());
+      }
+
       propOptions.actorList.template get<GsfActor>()
           .m_cfg.bethe_heitler_approx = &m_betheHeitlerApproximation;
 
@@ -166,7 +171,7 @@ struct GaussianSumFitter {
     };
 
     // Initialize the backward propagation with the DirectNavigator
-    auto bwdPropInitializer = [this](const auto& opts) {
+    auto bwdPropInitializer = [&](const auto& opts) {
       using Actors = ActorList<GsfActor, EndOfWorldReached>;
       using PropagatorOptions = typename propagator_t::template Options<Actors>;
 
@@ -174,6 +179,11 @@ struct GaussianSumFitter {
 
       propOptions.setPlainOptions(opts.propagatorPlainOptions);
 
+      for (auto it = begin; it != end; ++it) {
+        propOptions.navigation.insertExternalSurface(
+            options.extensions.surfaceAccessor(SourceLink{*it})->geometryId());
+      }
+
       propOptions.actorList.template get<GsfActor>()
           .m_cfg.bethe_heitler_approx = &m_betheHeitlerApproximation;