EIC code displayed by LXR master/​acts/​Tests/​UnitTests/​Core/​TrackFitting/​Gx2fTests.cpp	Source navigation Diff markup Identifier search General search
	Version: master test Revert   Change

File indexing completed on 2026-10-13 08:37:10

 // This file is part of the ACTS project.
 //
 // Copyright (C) 2016 CERN for the benefit of the ACTS project
 //
 // This Source Code Form is subject to the terms of the Mozilla Public
 // License, v. 2.0. If a copy of the MPL was not distributed with this
 // file, You can obtain one at https://mozilla.org/MPL/2.0/.
 
 #include <boost/test/unit_test.hpp>
 
 #include "Acts/Definitions/Algebra.hpp"
 #include "Acts/Definitions/Units.hpp"
 #include "Acts/EventData/Types.hpp"
 #include "Acts/EventData/VectorMultiTrajectory.hpp"
 #include "Acts/EventData/VectorTrackContainer.hpp"
 #include "Acts/EventData/detail/TestSourceLink.hpp"
 #include "Acts/Geometry/CuboidVolumeBuilder.hpp"
 #include "Acts/Geometry/TrackingGeometry.hpp"
 #include "Acts/Geometry/TrackingGeometryBuilder.hpp"
 #include "Acts/Material/HomogeneousSurfaceMaterial.hpp"
 #include "Acts/Material/MaterialSlab.hpp"
 #include "Acts/Propagator/EigenStepper.hpp"
 #include "Acts/Propagator/Navigator.hpp"
 #include "Acts/Propagator/Propagator.hpp"
 #include "Acts/Propagator/StraightLineStepper.hpp"
 #include "Acts/Surfaces/PlaneSurface.hpp"
 #include "Acts/Surfaces/RectangleBounds.hpp"
 #include "Acts/TrackFitting/GlobalChiSquareFitter.hpp"
 #include "Acts/Utilities/Logger.hpp"
 #include "Acts/Visualization/EventDataView3D.hpp"
 #include "Acts/Visualization/GeometryView3D.hpp"
 #include "Acts/Visualization/ObjVisualization3D.hpp"
 #include "ActsTests/CommonHelpers/DetectorElementStub.hpp"
 #include "ActsTests/CommonHelpers/MeasurementsCreator.hpp"
 #include "ActsTests/CommonHelpers/PredefinedMaterials.hpp"
 
 #include <numbers>
 #include <vector>
 
 #include "FitterTestsCommon.hpp"
 
 using namespace Acts;
 using namespace Acts::Experimental;
 using namespace Acts::UnitLiterals;
 using namespace Acts::detail::Test;
 
 Logging::Level logLevel = Logging::VERBOSE;
 const auto gx2fLogger = getDefaultLogger("Gx2f", logLevel);
 
 namespace ActsTests {
 
 /// @brief Helper function to visualise measurements in a 3D environment.
 ///
 /// This function iterates through the provided measurements and visualises each
 /// one using the specified 3D visualisation helper. The visualisation takes
 /// into account the surface transformations and localisation errors.
 ///
 /// @param helper The 3D visualisation helper used to draw the measurements.
 /// @param measurements A collection of measurements to be visualised, containing source links with parameters and covariance information.
 /// @param geometry A shared pointer to the constant tracking geometry used to find surfaces associated with measurements.
 /// @param geoCtx The geometry context used for transformations and accessing geometry-related information.
 /// @param locErrorScale Scaling factor for localisation errors. Default value is 1.0.
 /// @param viewConfig Configuration settings for the visualisation view. Default value is s_viewMeasurement.
 static void drawMeasurements(
     IVisualization3D& helper, const Measurements& measurements,
     const std::shared_ptr<const TrackingGeometry>& geometry,
     const GeometryContext& geoCtx, double locErrorScale = 1.,
     const ViewConfig& viewConfig = s_viewMeasurement) {
   std::cout << "\n*** Draw measurements ***\n" << std::endl;
 
   for (auto& singleMeasurement : measurements.sourceLinks) {
     auto cov = singleMeasurement.covariance;
     auto lposition = singleMeasurement.parameters;
 
     auto surf = geometry->findSurface(singleMeasurement.m_geometryId);
     auto transf = surf->localToGlobalTransform(geoCtx);
 
     EventDataView3D::drawMeasurement(helper, lposition, cov, transf,
                                      locErrorScale, viewConfig);
   }
 }
 
 //// Construct initial track parameters.
 BoundTrackParameters makeParameters(
     const double x = 0.0_m, const double y = 0.0_m, const double z = 0.0_m,
     const double w = 42_ns, const double phi = 0_degree,
     const double theta = 90_degree, const double p = 2_GeV,
     const double q = 1_e) {
   // create covariance matrix from reasonable standard deviations
   BoundVector stddev;
   stddev[eBoundLoc0] = 100_um;
   stddev[eBoundLoc1] = 100_um;
   stddev[eBoundTime] = 25_ns;
   stddev[eBoundPhi] = 2_degree;
   stddev[eBoundTheta] = 2_degree;
   stddev[eBoundQOverP] = 1 / 100_GeV;
   const BoundMatrix cov = stddev.cwiseProduct(stddev).asDiagonal();
   // define a track in the transverse plane along x
   const Vector4 mPos4(x, y, z, w);
   return BoundTrackParameters::createCurvilinear(mPos4, phi, theta, q / p, cov,
                                                  ParticleHypothesis::pion());
 }
 
 static std::vector<SourceLink> prepareSourceLinks(
     const std::vector<TestSourceLink>& sourceLinks) {
   std::vector<SourceLink> result;
   std::transform(sourceLinks.begin(), sourceLinks.end(),
                  std::back_inserter(result),
                  [](const auto& sl) { return SourceLink{sl}; });
   return result;
 }
 
 /// @brief Create a simple telescope detector
 ///
 /// @param geoCtx
 /// @param nSurfaces Number of surfaces
 /// @param surfaceIndexWithMaterial A set of index of the material surfaces
 std::shared_ptr<const TrackingGeometry> makeToyDetector(
     const GeometryContext& geoCtx, const std::size_t nSurfaces = 5,
     const std::set<std::size_t>& surfaceIndexWithMaterial = {}) {
   if (nSurfaces < 1) {
     throw std::invalid_argument("At least 1 surfaces needs to be created.");
   }
 
   // Define the dimensions of the square surfaces
   const double halfSizeSurface = 1_m;
 
   // Rotation of the surfaces around the y-axis
   const double rotationAngle = std::numbers::pi / 2.;
   const Vector3 xPos(std::cos(rotationAngle), 0., std::sin(rotationAngle));
   const Vector3 yPos(0., 1., 0.);
   const Vector3 zPos(-std::sin(rotationAngle), 0., std::cos(rotationAngle));
 
   // Construct builder
   CuboidVolumeBuilder cvb;
 
   // Create configurations for surfaces
   std::vector<CuboidVolumeBuilder::SurfaceConfig> surfaceConfig;
   for (std::size_t surfPos = 1; surfPos <= nSurfaces; surfPos++) {
     // Position of the surfaces
     CuboidVolumeBuilder::SurfaceConfig cfg;
     cfg.position = {surfPos * UnitConstants::m, 0., 0.};
 
     // Rotation of the surfaces
     cfg.rotation.col(0) = xPos;
     cfg.rotation.col(1) = yPos;
     cfg.rotation.col(2) = zPos;
 
     // Boundaries of the surfaces (shape)
     cfg.rBounds = std::make_shared<const RectangleBounds>(
         RectangleBounds(halfSizeSurface, halfSizeSurface));
 
     // Add material only for selected surfaces
     if (surfaceIndexWithMaterial.contains(surfPos)) {
       // Material of the surfaces
       MaterialSlab matProp(makeSilicon(), 5_mm);
       cfg.surMat = std::make_shared<HomogeneousSurfaceMaterial>(matProp);
     }
 
     // Thickness of the detector element
     cfg.thickness = 1_um;
 
     cfg.detElementConstructor =
         [](const Transform3& trans,
            const std::shared_ptr<const RectangleBounds>& bounds,
            double thickness) {
           return new DetectorElementStub(trans, bounds, thickness);
         };
     surfaceConfig.push_back(cfg);
   }
 
   // Build layer configurations
   std::vector<CuboidVolumeBuilder::LayerConfig> layerConfig;
   for (auto& sCfg : surfaceConfig) {
     CuboidVolumeBuilder::LayerConfig cfg;
     cfg.surfaceCfg = {sCfg};
     cfg.active = true;
     cfg.envelopeX = {-0.1_mm, 0.1_mm};
     cfg.envelopeY = {-0.1_mm, 0.1_mm};
     cfg.envelopeZ = {-0.1_mm, 0.1_mm};
     layerConfig.push_back(cfg);
   }
 
   // Inner Volume - Build volume configuration
   CuboidVolumeBuilder::VolumeConfig volumeConfig;
   volumeConfig.length = {(nSurfaces + 1) * 1_m, 2 * halfSizeSurface,
                          2 * halfSizeSurface};
   volumeConfig.position = {volumeConfig.length.x() / 2, 0., 0.};
   volumeConfig.layerCfg = layerConfig;
   volumeConfig.name = "TestVolume";
 
   // Outer volume - Build TrackingGeometry configuration
   CuboidVolumeBuilder::Config config;
   config.length = {(nSurfaces + 1) * 1_m, 2 * halfSizeSurface,
                    2 * halfSizeSurface};
   config.position = {volumeConfig.length.x() / 2, 0., 0.};
   config.volumeCfg = {volumeConfig};
 
   cvb.setConfig(config);
 
   TrackingGeometryBuilder::Config tgbCfg;
 
   tgbCfg.trackingVolumeBuilders.push_back(
       [=](const auto& context, const auto& inner, const auto&) {
         return cvb.trackingVolume(context, inner, nullptr);
       });
 
   TrackingGeometryBuilder tgb(tgbCfg);
 
   std::unique_ptr<const TrackingGeometry> detector =
       tgb.trackingGeometry(geoCtx);
   return detector;
 }
 
 struct Detector {
   // geometry
   std::shared_ptr<const TrackingGeometry> geometry;
 };
 
 BOOST_AUTO_TEST_SUITE(TrackFittingSuite)
 
 ACTS_LOCAL_LOGGER(getDefaultLogger("Gx2fTests", logLevel))
 
 // Context objects
 const auto geoCtx = GeometryContext::dangerouslyDefaultConstruct();
 const MagneticFieldContext magCtx;
 const CalibrationContext calCtx;
 
 // Measurement resolutions
 const MeasurementResolution resPixel = {MeasurementType::eLoc01,
                                         {25_um, 50_um}};
 const MeasurementResolution resStrip0 = {MeasurementType::eLoc0, {25_um}};
 const MeasurementResolution resStrip1 = {MeasurementType::eLoc1, {50_um}};
 const MeasurementResolutionMap resMapAllPixel = {
     {GeometryIdentifier().withVolume(0), resPixel}};
 
 // This test checks if the call to the fitter works and no errors occur in the
 // framework, without fitting and updating any parameters
 BOOST_AUTO_TEST_CASE(NoFit) {
   ACTS_INFO("*** Test: NoFit -- Start");
 
   std::default_random_engine rng(42);
 
   ACTS_DEBUG("Create the detector");
   const std::size_t nSurfaces = 5;
   Detector detector;
   detector.geometry = makeToyDetector(geoCtx, nSurfaces);
 
   ACTS_DEBUG("Set the start parameters for measurement creation and fit");
   const auto parametersMeasurements = makeParameters();
   const auto startParametersFit = makeParameters(
       7_mm, 11_mm, 15_mm, 42_ns, 10_degree, 80_degree, 1_GeV, 1_e);
 
   ACTS_DEBUG("Create the measurements");
   using SimPropagator = Propagator<StraightLineStepper, Navigator>;
   const SimPropagator simPropagator = makeStraightPropagator(detector.geometry);
   const auto measurements =
       createMeasurements(simPropagator, geoCtx, magCtx, parametersMeasurements,
                          resMapAllPixel, rng);
   const auto sourceLinks = prepareSourceLinks(measurements.sourceLinks);
   ACTS_VERBOSE("sourceLinks.size() = " << sourceLinks.size());
 
   ACTS_DEBUG("Set up the fitter");
   // Reuse the SimPropagator, since we will not actually use it
   using Gx2Fitter = Gx2Fitter<SimPropagator, VectorMultiTrajectory>;
   const Gx2Fitter fitter(simPropagator, gx2fLogger->clone());
 
   const Surface* rSurface = &parametersMeasurements.referenceSurface();
 
   Gx2FitterExtensions<VectorMultiTrajectory> extensions;
   extensions.calibrator
       .connect<&testSourceLinkCalibratorStrict<VectorMultiTrajectory>>();
   TestSourceLink::SurfaceAccessor surfaceAccessor{*detector.geometry};
   extensions.surfaceAccessor
       .connect<&TestSourceLink::SurfaceAccessor::operator()>(&surfaceAccessor);
 
   Gx2FitterOptions gx2fOptions(geoCtx, magCtx, calCtx, extensions,
                                PropagatorPlainOptions(geoCtx, magCtx), rSurface,
                                false, false, FreeToBoundCorrection(false), 0,
                                0);
 
   TrackContainer tracks{VectorTrackContainer{}, VectorMultiTrajectory{}};
 
   ACTS_DEBUG("Fit the track");
   ACTS_VERBOSE("startParameter unsmeared:\n" << parametersMeasurements);
   ACTS_VERBOSE("startParameter fit:\n" << startParametersFit);
   const auto res = fitter.fit(sourceLinks.begin(), sourceLinks.end(),
                               startParametersFit, gx2fOptions, tracks);
 
   BOOST_REQUIRE(res.ok());
 
   const auto& track = *res;
   BOOST_CHECK_EQUAL(track.tipIndex(), kTrackIndexInvalid);
   BOOST_CHECK(track.hasReferenceSurface());
 
   // Track quantities
   BOOST_CHECK_EQUAL(track.chi2(), 0.);
   BOOST_CHECK_EQUAL(track.nDoF(), 0u);
   BOOST_CHECK_EQUAL(track.nHoles(), 0u);
   BOOST_CHECK_EQUAL(track.nMeasurements(), 0u);
   BOOST_CHECK_EQUAL(track.nSharedHits(), 0u);
   BOOST_CHECK_EQUAL(track.nOutliers(), 0u);
 
   // Parameters
   BOOST_CHECK_EQUAL(track.parameters(), startParametersFit.parameters());
   BOOST_CHECK_EQUAL(track.covariance(), BoundMatrix::Identity());
 
   // Convergence
   BOOST_CHECK_EQUAL(
       (track.template component<
           std::uint32_t, hashString(Gx2fConstants::gx2fnUpdateColumn)>()),
       0);
 
   ACTS_INFO("*** Test: NoFit -- Finish");
 }
 
 BOOST_AUTO_TEST_CASE(Fit5Iterations) {
   ACTS_INFO("*** Test: Fit5Iterations -- Start");
 
   std::default_random_engine rng(42);
 
   ACTS_DEBUG("Create the detector");
   const std::size_t nSurfaces = 5;
   Detector detector;
   detector.geometry = makeToyDetector(geoCtx, nSurfaces);
 
   ACTS_DEBUG("Set the start parameters for measurement creation and fit");
   const auto parametersMeasurements = makeParameters();
   const auto startParametersFit = makeParameters(
       7_mm, 11_mm, 15_mm, 42_ns, 10_degree, 80_degree, 1_GeV, 1_e);
 
   ACTS_DEBUG("Create the measurements");
   using SimPropagator = Propagator<StraightLineStepper, Navigator>;
   const SimPropagator simPropagator = makeStraightPropagator(detector.geometry);
   const auto measurements =
       createMeasurements(simPropagator, geoCtx, magCtx, parametersMeasurements,
                          resMapAllPixel, rng);
   const auto sourceLinks = prepareSourceLinks(measurements.sourceLinks);
   ACTS_VERBOSE("sourceLinks.size() = " << sourceLinks.size());
 
   BOOST_REQUIRE_EQUAL(sourceLinks.size(), nSurfaces);
 
   ACTS_DEBUG("Set up the fitter");
   const Surface* rSurface = &parametersMeasurements.referenceSurface();
 
   using RecoStepper = EigenStepper<>;
   const auto recoPropagator =
       makeConstantFieldPropagator<RecoStepper>(detector.geometry, 0_T);
 
   using RecoPropagator = decltype(recoPropagator);
   using Gx2Fitter = Gx2Fitter<RecoPropagator, VectorMultiTrajectory>;
   const Gx2Fitter fitter(recoPropagator, gx2fLogger->clone());
 
   Gx2FitterExtensions<VectorMultiTrajectory> extensions;
   extensions.calibrator
       .connect<&testSourceLinkCalibratorStrict<VectorMultiTrajectory>>();
   TestSourceLink::SurfaceAccessor surfaceAccessor{*detector.geometry};
   extensions.surfaceAccessor
       .connect<&TestSourceLink::SurfaceAccessor::operator()>(&surfaceAccessor);
 
   const Gx2FitterOptions gx2fOptions(geoCtx, magCtx, calCtx, extensions,
                                      PropagatorPlainOptions(geoCtx, magCtx),
                                      rSurface, false, false,
                                      FreeToBoundCorrection(false), 5, 0);
 
   TrackContainer tracks{VectorTrackContainer{}, VectorMultiTrajectory{}};
 
   ACTS_DEBUG("Fit the track");
   ACTS_VERBOSE("startParameter unsmeared:\n" << parametersMeasurements);
   ACTS_VERBOSE("startParameter fit:\n" << startParametersFit);
   const auto res = fitter.fit(sourceLinks.begin(), sourceLinks.end(),
                               startParametersFit, gx2fOptions, tracks);
 
   BOOST_REQUIRE(res.ok());
 
   const auto& track = *res;
 
   BOOST_CHECK_EQUAL(track.tipIndex(), nSurfaces - 1);
   BOOST_CHECK(track.hasReferenceSurface());
 
   // Track quantities
   CHECK_CLOSE_ABS(track.chi2(), 8., 2.);
   BOOST_CHECK_EQUAL(track.nDoF(), nSurfaces * 2);
   BOOST_CHECK_EQUAL(track.nHoles(), 0u);
   BOOST_CHECK_EQUAL(track.nMeasurements(), nSurfaces);
   BOOST_CHECK_EQUAL(track.nSharedHits(), 0u);
   BOOST_CHECK_EQUAL(track.nOutliers(), 0u);
 
   // Parameters
   // We need quite coarse checks here, since on different builds
   // the created measurements differ in the randomness
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc0], -11., 7e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc1], -15., 6e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundPhi], 1e-5, 1e3);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTheta], std::numbers::pi / 2,
                     1e-3);
   BOOST_CHECK_EQUAL(track.parameters()[eBoundQOverP], 1);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTime],
                     startParametersFit.parameters()[eBoundTime], 1e-6);
   BOOST_CHECK_CLOSE(track.covariance().determinant(), 1e-27, 4e0);
 
   // Convergence
   BOOST_CHECK_EQUAL(
       (track.template component<
           std::uint32_t, hashString(Gx2fConstants::gx2fnUpdateColumn)>()),
       4);
 
   ACTS_INFO("*** Test: Fit5Iterations -- Finish");
 }
 
 // This test checks that the fitter reaches a measurement on a free surface,
 // which is not part of the tracking geometry
 BOOST_AUTO_TEST_CASE(FreeMeasurementSurface) {
   ACTS_INFO("*** Test: FreeMeasurementSurface -- Start");
 
   std::default_random_engine rng(42);
 
   ACTS_DEBUG("Create the detector");
   const std::size_t nSurfaces = 5;
   Detector detector;
   detector.geometry = makeToyDetector(geoCtx, nSurfaces);
 
   ACTS_DEBUG("Set the start parameters for measurement creation and fit");
   const auto parametersMeasurements = makeParameters();
   const auto startParametersFit = makeParameters(
       7_mm, 11_mm, 15_mm, 42_ns, 10_degree, 80_degree, 1_GeV, 1_e);
 
   ACTS_DEBUG("Create the measurements");
   using SimPropagator = Propagator<StraightLineStepper, Navigator>;
   const SimPropagator simPropagator = makeStraightPropagator(detector.geometry);
   const auto measurements =
       createMeasurements(simPropagator, geoCtx, magCtx, parametersMeasurements,
                          resMapAllPixel, rng);
   auto sourceLinks = prepareSourceLinks(measurements.sourceLinks);
   BOOST_REQUIRE_EQUAL(sourceLinks.size(), nSurfaces);
 
   ACTS_DEBUG("Add a measurement on a free surface between two layers");
   const Surface* firstSurface = detector.geometry->findSurface(
       measurements.sourceLinks.front().m_geometryId);
   BOOST_REQUIRE(firstSurface != nullptr);
   Transform3 freeTransform = Transform3::Identity();
   freeTransform.linear() =
       firstSurface->localToGlobalTransform(geoCtx).linear();
   freeTransform.translation() = Vector3(2.5_m, 0., 0.);
   const auto freeSurface = Surface::makeShared<PlaneSurface>(
       freeTransform, std::make_shared<RectangleBounds>(1_m, 1_m));
   BOOST_REQUIRE_EQUAL(freeSurface->geometryId(), GeometryIdentifier{});
   sourceLinks.emplace_back(
       TestSourceLink(eBoundLoc0, eBoundLoc1, Vector2::Zero(),
                      Vector2(25_um * 25_um, 50_um * 50_um).asDiagonal()));
 
   // Resolve the source link without a geometry id to the free surface
   struct FreeSurfaceAccessor {
     const TrackingGeometry* geometry{};
     const Surface* freeSurface{};
 
     const Surface* operator()(const SourceLink& sourceLink) const {
       const GeometryIdentifier geoId =
           sourceLink.get<TestSourceLink>().m_geometryId;
       return geoId == GeometryIdentifier{} ? freeSurface
                                            : geometry->findSurface(geoId);
     }
   };
   FreeSurfaceAccessor surfaceAccessor{detector.geometry.get(),
                                       freeSurface.get()};
 
   ACTS_DEBUG("Set up the fitter");
   const Surface* rSurface = &parametersMeasurements.referenceSurface();
 
   using RecoStepper = EigenStepper<>;
   const auto recoPropagator =
       makeConstantFieldPropagator<RecoStepper>(detector.geometry, 0_T);
 
   using RecoPropagator = decltype(recoPropagator);
   using Gx2Fitter = Gx2Fitter<RecoPropagator, VectorMultiTrajectory>;
   const Gx2Fitter fitter(recoPropagator, gx2fLogger->clone());
 
   Gx2FitterExtensions<VectorMultiTrajectory> extensions;
   extensions.calibrator
       .connect<&testSourceLinkCalibratorStrict<VectorMultiTrajectory>>();
   extensions.surfaceAccessor.connect<&FreeSurfaceAccessor::operator()>(
       &surfaceAccessor);
 
   const Gx2FitterOptions gx2fOptions(geoCtx, magCtx, calCtx, extensions,
                                      PropagatorPlainOptions(geoCtx, magCtx),
                                      rSurface, false, false,
                                      FreeToBoundCorrection(false), 5, 0);
 
   TrackContainer tracks{VectorTrackContainer{}, VectorMultiTrajectory{}};
 
   ACTS_DEBUG("Fit the track");
   const auto res = fitter.fit(sourceLinks.begin(), sourceLinks.end(),
                               startParametersFit, gx2fOptions, tracks);
 
   BOOST_REQUIRE(res.ok());
 
   const auto& track = *res;
   BOOST_CHECK_EQUAL(track.nMeasurements(), nSurfaces + 1);
   BOOST_CHECK_EQUAL(track.nDoF(), (nSurfaces + 1) * 2);
   BOOST_CHECK_EQUAL(track.nHoles(), 0u);
   BOOST_CHECK_EQUAL(track.nOutliers(), 0u);
 
   // The free surface carries a measurement state
   bool freeSurfaceMeasured = false;
   for (const auto trackState : track.trackStatesReversed()) {
     if (&trackState.referenceSurface() == freeSurface.get()) {
       freeSurfaceMeasured = trackState.typeFlags().isMeasurement();
     }
   }
   BOOST_CHECK(freeSurfaceMeasured);
 
   ACTS_INFO("*** Test: FreeMeasurementSurface -- Finish");
 }
 
 BOOST_AUTO_TEST_CASE(MixedDetector) {
   ACTS_INFO("*** Test: MixedDetector -- Start");
 
   std::default_random_engine rng(42);
 
   ACTS_DEBUG("Create the detector");
   const std::size_t nSurfaces = 7;
   Detector detector;
   detector.geometry = makeToyDetector(geoCtx, nSurfaces);
 
   ACTS_DEBUG("Set the start parameters for measurement creation and fit");
   const auto parametersMeasurements = makeParameters();
   const auto startParametersFit = makeParameters(
       7_mm, 11_mm, 15_mm, 42_ns, 10_degree, 80_degree, 1_GeV, 1_e);
 
   ACTS_DEBUG("Create the measurements");
   const MeasurementResolutionMap resMap = {
       {GeometryIdentifier().withVolume(2).withLayer(2), resPixel},
       {GeometryIdentifier().withVolume(2).withLayer(4), resStrip0},
       {GeometryIdentifier().withVolume(2).withLayer(6), resStrip1},
       {GeometryIdentifier().withVolume(2).withLayer(8), resPixel},
       {GeometryIdentifier().withVolume(2).withLayer(10), resStrip0},
       {GeometryIdentifier().withVolume(2).withLayer(12), resStrip1},
       {GeometryIdentifier().withVolume(2).withLayer(14), resPixel},
   };
 
   using SimPropagator = Propagator<StraightLineStepper, Navigator>;
   const SimPropagator simPropagator = makeStraightPropagator(detector.geometry);
   const auto measurements = createMeasurements(
       simPropagator, geoCtx, magCtx, parametersMeasurements, resMap, rng);
   const auto sourceLinks = prepareSourceLinks(measurements.sourceLinks);
   ACTS_VERBOSE("sourceLinks.size() = " << sourceLinks.size());
 
   BOOST_REQUIRE_EQUAL(sourceLinks.size(), nSurfaces);
 
   ACTS_DEBUG("Set up the fitter");
   const Surface* rSurface = &parametersMeasurements.referenceSurface();
 
   using RecoStepper = EigenStepper<>;
   const auto recoPropagator =
       makeConstantFieldPropagator<RecoStepper>(detector.geometry, 0_T);
 
   using RecoPropagator = decltype(recoPropagator);
   using Gx2Fitter = Gx2Fitter<RecoPropagator, VectorMultiTrajectory>;
   const Gx2Fitter fitter(recoPropagator, gx2fLogger->clone());
 
   Gx2FitterExtensions<VectorMultiTrajectory> extensions;
   extensions.calibrator
       .connect<&testSourceLinkCalibratorStrict<VectorMultiTrajectory>>();
   TestSourceLink::SurfaceAccessor surfaceAccessor{*detector.geometry};
   extensions.surfaceAccessor
       .connect<&TestSourceLink::SurfaceAccessor::operator()>(&surfaceAccessor);
 
   const Gx2FitterOptions gx2fOptions(geoCtx, magCtx, calCtx, extensions,
                                      PropagatorPlainOptions(geoCtx, magCtx),
                                      rSurface, false, false,
                                      FreeToBoundCorrection(false), 5, 0);
 
   TrackContainer tracks{VectorTrackContainer{}, VectorMultiTrajectory{}};
 
   ACTS_DEBUG("Fit the track");
   ACTS_VERBOSE("startParameter unsmeared:\n" << parametersMeasurements);
   ACTS_VERBOSE("startParameter fit:\n" << startParametersFit);
   const auto res = fitter.fit(sourceLinks.begin(), sourceLinks.end(),
                               startParametersFit, gx2fOptions, tracks);
 
   BOOST_REQUIRE(res.ok());
 
   const auto& track = *res;
   BOOST_CHECK_EQUAL(track.tipIndex(), nSurfaces - 1);
   BOOST_CHECK(track.hasReferenceSurface());
 
   // Track quantities
   CHECK_CLOSE_ABS(track.chi2(), 8.5, 4.);
   BOOST_CHECK_EQUAL(track.nDoF(), 10u);
   BOOST_CHECK_EQUAL(track.nHoles(), 0u);
   BOOST_CHECK_EQUAL(track.nMeasurements(), nSurfaces);
   BOOST_CHECK_EQUAL(track.nSharedHits(), 0u);
   BOOST_CHECK_EQUAL(track.nOutliers(), 0u);
 
   // Parameters
   // We need quite coarse checks here, since on different builds
   // the created measurements differ in the randomness
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc0], -11., 7e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc1], -15., 6e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundPhi], 1e-5, 1e3);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTheta], std::numbers::pi / 2,
                     1e-3);
   BOOST_CHECK_EQUAL(track.parameters()[eBoundQOverP], 1);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTime],
                     startParametersFit.parameters()[eBoundTime], 1e-6);
   BOOST_CHECK_CLOSE(track.covariance().determinant(), 2e-28, 1e0);
 
   // Convergence
   BOOST_CHECK_EQUAL(
       (track.template component<
           std::uint32_t, hashString(Gx2fConstants::gx2fnUpdateColumn)>()),
       4);
 
   ACTS_INFO("*** Test: MixedDetector -- Finish");
 }
 
 // This test checks if we can fit QOverP, when a magnetic field is introduced
 BOOST_AUTO_TEST_CASE(FitWithBfield) {
   ACTS_INFO("*** Test: FitWithBfield -- Start");
 
   std::default_random_engine rng(42);
 
   ACTS_DEBUG("Create the detector");
   const std::size_t nSurfaces = 5;
   Detector detector;
   detector.geometry = makeToyDetector(geoCtx, nSurfaces);
 
   ACTS_DEBUG("Set the start parameters for measurement creation and fit");
   const auto parametersMeasurements = makeParameters();
   const auto startParametersFit = makeParameters(
       7_mm, 11_mm, 15_mm, 42_ns, 10_degree, 80_degree, 1_GeV, 1_e);
 
   ACTS_DEBUG("Create the measurements");
   using SimStepper = EigenStepper<>;
   const auto simPropagator =
       makeConstantFieldPropagator<SimStepper>(detector.geometry, 0.3_T);
 
   const auto measurements =
       createMeasurements(simPropagator, geoCtx, magCtx, parametersMeasurements,
                          resMapAllPixel, rng);
 
   const auto sourceLinks = prepareSourceLinks(measurements.sourceLinks);
   ACTS_VERBOSE("sourceLinks.size() = " << sourceLinks.size());
 
   BOOST_REQUIRE_EQUAL(sourceLinks.size(), nSurfaces);
 
   ACTS_DEBUG("Set up the fitter");
   const Surface* rSurface = &parametersMeasurements.referenceSurface();
 
   // Reuse the SimPropagator, since it already uses the EigenStepper<>
   using SimPropagator = decltype(simPropagator);
   using Gx2Fitter = Gx2Fitter<SimPropagator, VectorMultiTrajectory>;
   const Gx2Fitter fitter(simPropagator, gx2fLogger->clone());
 
   Gx2FitterExtensions<VectorMultiTrajectory> extensions;
   extensions.calibrator
       .connect<&testSourceLinkCalibratorStrict<VectorMultiTrajectory>>();
   TestSourceLink::SurfaceAccessor surfaceAccessor{*detector.geometry};
   extensions.surfaceAccessor
       .connect<&TestSourceLink::SurfaceAccessor::operator()>(&surfaceAccessor);
 
   const Gx2FitterOptions gx2fOptions(geoCtx, magCtx, calCtx, extensions,
                                      PropagatorPlainOptions(geoCtx, magCtx),
                                      rSurface, false, false,
                                      FreeToBoundCorrection(false), 5, 0);
 
   TrackContainer tracks{VectorTrackContainer{}, VectorMultiTrajectory{}};
 
   ACTS_DEBUG("Fit the track");
   ACTS_VERBOSE("startParameter unsmeared:\n" << parametersMeasurements);
   ACTS_VERBOSE("startParameter fit:\n" << startParametersFit);
   const auto res = fitter.fit(sourceLinks.begin(), sourceLinks.end(),
                               startParametersFit, gx2fOptions, tracks);
 
   BOOST_REQUIRE(res.ok());
 
   const auto& track = *res;
   BOOST_CHECK_EQUAL(track.tipIndex(), nSurfaces - 1);
   BOOST_CHECK(track.hasReferenceSurface());
 
   // Track quantities
   CHECK_CLOSE_ABS(track.chi2(), 7.5, 1.5);
   BOOST_CHECK_EQUAL(track.nDoF(), nSurfaces * 2);
   BOOST_CHECK_EQUAL(track.nHoles(), 0u);
   BOOST_CHECK_EQUAL(track.nMeasurements(), nSurfaces);
   BOOST_CHECK_EQUAL(track.nSharedHits(), 0u);
   BOOST_CHECK_EQUAL(track.nOutliers(), 0u);
 
   // Parameters
   // We need quite coarse checks here, since on different builds
   // the created measurements differ in the randomness
   // TODO investigate further the reference values for eBoundPhi and det(cov)
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc0], -11., 8e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc1], -15., 6e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundPhi], 1e-4, 1e3);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTheta], std::numbers::pi / 2,
                     1e-3);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundQOverP], 0.5, 2e-1);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTime],
                     startParametersFit.parameters()[eBoundTime], 1e-6);
   BOOST_CHECK_CLOSE(track.covariance().determinant(), 8e-35, 4e0);
 
   // Convergence
   BOOST_CHECK_EQUAL(
       (track.template component<
           std::uint32_t, hashString(Gx2fConstants::gx2fnUpdateColumn)>()),
       4);
 
   ACTS_INFO("*** Test: FitWithBfield -- Finish");
 }
 
 BOOST_AUTO_TEST_CASE(relChi2changeCutOff) {
   ACTS_INFO("*** Test: relChi2changeCutOff -- Start");
 
   std::default_random_engine rng(42);
 
   ACTS_DEBUG("Create the detector");
   const std::size_t nSurfaces = 5;
   Detector detector;
   detector.geometry = makeToyDetector(geoCtx, nSurfaces);
 
   ACTS_DEBUG("Set the start parameters for measurement creation and fit");
   const auto parametersMeasurements = makeParameters();
   const auto startParametersFit = makeParameters(
       7_mm, 11_mm, 15_mm, 42_ns, 10_degree, 80_degree, 1_GeV, 1_e);
 
   ACTS_DEBUG("Create the measurements");
   // simulation propagator
   using SimPropagator = Propagator<StraightLineStepper, Navigator>;
   const SimPropagator simPropagator = makeStraightPropagator(detector.geometry);
   const auto measurements =
       createMeasurements(simPropagator, geoCtx, magCtx, parametersMeasurements,
                          resMapAllPixel, rng);
   const auto sourceLinks = prepareSourceLinks(measurements.sourceLinks);
   ACTS_VERBOSE("sourceLinks.size() = " << sourceLinks.size());
 
   BOOST_REQUIRE_EQUAL(sourceLinks.size(), nSurfaces);
 
   ACTS_DEBUG("Set up the fitter");
   const Surface* rSurface = &parametersMeasurements.referenceSurface();
 
   using RecoStepper = EigenStepper<>;
   const auto recoPropagator =
       makeConstantFieldPropagator<RecoStepper>(detector.geometry, 0_T);
 
   using RecoPropagator = decltype(recoPropagator);
   using Gx2Fitter = Gx2Fitter<RecoPropagator, VectorMultiTrajectory>;
   const Gx2Fitter fitter(recoPropagator, gx2fLogger->clone());
 
   Gx2FitterExtensions<VectorMultiTrajectory> extensions;
   extensions.calibrator
       .connect<&testSourceLinkCalibratorStrict<VectorMultiTrajectory>>();
   TestSourceLink::SurfaceAccessor surfaceAccessor{*detector.geometry};
   extensions.surfaceAccessor
       .connect<&TestSourceLink::SurfaceAccessor::operator()>(&surfaceAccessor);
 
   const Gx2FitterOptions gx2fOptions(geoCtx, magCtx, calCtx, extensions,
                                      PropagatorPlainOptions(geoCtx, magCtx),
                                      rSurface, false, false,
                                      FreeToBoundCorrection(false), 500, 1e-5);
 
   TrackContainer tracks{VectorTrackContainer{}, VectorMultiTrajectory{}};
 
   ACTS_DEBUG("Fit the track");
   ACTS_VERBOSE("startParameter unsmeared:\n" << parametersMeasurements);
   ACTS_VERBOSE("startParameter fit:\n" << startParametersFit);
   const auto res = fitter.fit(sourceLinks.begin(), sourceLinks.end(),
                               startParametersFit, gx2fOptions, tracks);
 
   BOOST_REQUIRE(res.ok());
 
   const auto& track = *res;
   BOOST_CHECK_EQUAL(track.tipIndex(), nSurfaces - 1);
   BOOST_CHECK(track.hasReferenceSurface());
 
   // Track quantities
   CHECK_CLOSE_ABS(track.chi2(), 8., 2.);
   BOOST_CHECK_EQUAL(track.nDoF(), nSurfaces * 2);
   BOOST_CHECK_EQUAL(track.nHoles(), 0u);
   BOOST_CHECK_EQUAL(track.nMeasurements(), nSurfaces);
   BOOST_CHECK_EQUAL(track.nSharedHits(), 0u);
   BOOST_CHECK_EQUAL(track.nOutliers(), 0u);
 
   // Parameters
   // We need quite coarse checks here, since on different builds
   // the created measurements differ in the randomness
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc0], -11., 7e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc1], -15., 6e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundPhi], 1e-5, 1e3);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTheta], std::numbers::pi / 2,
                     1e-3);
   BOOST_CHECK_EQUAL(track.parameters()[eBoundQOverP], 1);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTime],
                     startParametersFit.parameters()[eBoundTime], 1e-6);
   BOOST_CHECK_CLOSE(track.covariance().determinant(), 1e-27, 4e0);
 
   // Convergence
   BOOST_CHECK_LT(
       (track.template component<
           std::uint32_t, hashString(Gx2fConstants::gx2fnUpdateColumn)>()),
       10);
 
   ACTS_INFO("*** Test: relChi2changeCutOff -- Finish");
 }
 
 BOOST_AUTO_TEST_CASE(DidNotConverge) {
   ACTS_INFO("*** Test: DidNotConverge -- Start");
 
   std::default_random_engine rng(42);
 
   ACTS_DEBUG("Create the detector");
   const std::size_t nSurfaces = 5;
   Detector detector;
   detector.geometry = makeToyDetector(geoCtx, nSurfaces);
 
   ACTS_DEBUG("Set the start parameters for measurement creation and fit");
   const auto parametersMeasurements = makeParameters();
   const auto startParametersFit = makeParameters(
       7_mm, 11_mm, 15_mm, 42_ns, 10_degree, 80_degree, 1_GeV, 1_e);
 
   ACTS_DEBUG("Create the measurements");
   // simulation propagator
   using SimPropagator = Propagator<StraightLineStepper, Navigator>;
   const SimPropagator simPropagator = makeStraightPropagator(detector.geometry);
   const auto measurements =
       createMeasurements(simPropagator, geoCtx, magCtx, parametersMeasurements,
                          resMapAllPixel, rng);
   const auto sourceLinks = prepareSourceLinks(measurements.sourceLinks);
   ACTS_VERBOSE("sourceLinks.size() = " << sourceLinks.size());
 
   BOOST_REQUIRE_EQUAL(sourceLinks.size(), nSurfaces);
 
   ACTS_DEBUG("Set up the fitter");
   const Surface* rSurface = &parametersMeasurements.referenceSurface();
 
   using RecoStepper = EigenStepper<>;
   const auto recoPropagator =
       makeConstantFieldPropagator<RecoStepper>(detector.geometry, 0_T);
 
   using RecoPropagator = decltype(recoPropagator);
   using Gx2Fitter = Gx2Fitter<RecoPropagator, VectorMultiTrajectory>;
   const Gx2Fitter fitter(recoPropagator, gx2fLogger->clone());
 
   Gx2FitterExtensions<VectorMultiTrajectory> extensions;
   extensions.calibrator
       .connect<&testSourceLinkCalibratorStrict<VectorMultiTrajectory>>();
   TestSourceLink::SurfaceAccessor surfaceAccessor{*detector.geometry};
   extensions.surfaceAccessor
       .connect<&TestSourceLink::SurfaceAccessor::operator()>(&surfaceAccessor);
 
   // The relChi2changeCutOff = 0 prevents to stop the fitter after convergence,
   // therefore all updates will be done (even if the result does not change).
   // Since we didn't break due to convergence, we reach nUpdatesMax and
   // therefore fail the fit.
   const Gx2FitterOptions gx2fOptions(geoCtx, magCtx, calCtx, extensions,
                                      PropagatorPlainOptions(geoCtx, magCtx),
                                      rSurface, false, false,
                                      FreeToBoundCorrection(false), 6, 0);
 
   TrackContainer tracks{VectorTrackContainer{}, VectorMultiTrajectory{}};
 
   ACTS_DEBUG("Fit the track");
   ACTS_VERBOSE("startParameter unsmeared:\n" << parametersMeasurements);
   ACTS_VERBOSE("startParameter fit:\n" << startParametersFit);
   const auto res = fitter.fit(sourceLinks.begin(), sourceLinks.end(),
                               startParametersFit, gx2fOptions, tracks);
 
   BOOST_REQUIRE(!res.ok());
   BOOST_CHECK_EQUAL(res.error(), GlobalChiSquareFitterError::DidNotConverge);
 
   ACTS_INFO("*** Test: DidNotConverge -- Finish");
 }
 
 BOOST_AUTO_TEST_CASE(NotEnoughMeasurements) {
   ACTS_INFO("*** Test: NotEnoughMeasurements -- Start");
 
   std::default_random_engine rng(42);
 
   ACTS_DEBUG("Create the detector");
   const std::size_t nSurfaces = 2;
   Detector detector;
   detector.geometry = makeToyDetector(geoCtx, nSurfaces);
 
   ACTS_DEBUG("Set the start parameters for measurement creation and fit");
   const auto parametersMeasurements = makeParameters();
   const auto startParametersFit = makeParameters(
       7_mm, 11_mm, 15_mm, 42_ns, 10_degree, 80_degree, 1_GeV, 1_e);
 
   ACTS_DEBUG("Create the measurements");
   // simulation propagator
   using SimPropagator = Propagator<StraightLineStepper, Navigator>;
   const SimPropagator simPropagator = makeStraightPropagator(detector.geometry);
   const auto measurements =
       createMeasurements(simPropagator, geoCtx, magCtx, parametersMeasurements,
                          resMapAllPixel, rng);
   const auto sourceLinks = prepareSourceLinks(measurements.sourceLinks);
   ACTS_VERBOSE("sourceLinks.size() = " << sourceLinks.size());
 
   BOOST_REQUIRE_EQUAL(sourceLinks.size(), nSurfaces);
 
   ACTS_DEBUG("Set up the fitter");
   const Surface* rSurface = &parametersMeasurements.referenceSurface();
 
   using RecoStepper = EigenStepper<>;
   const auto recoPropagator =
       makeConstantFieldPropagator<RecoStepper>(detector.geometry, 0_T);
 
   using RecoPropagator = decltype(recoPropagator);
   using Gx2Fitter = Gx2Fitter<RecoPropagator, VectorMultiTrajectory>;
   const Gx2Fitter fitter(recoPropagator, gx2fLogger->clone());
 
   Gx2FitterExtensions<VectorMultiTrajectory> extensions;
   extensions.calibrator
       .connect<&testSourceLinkCalibratorStrict<VectorMultiTrajectory>>();
   TestSourceLink::SurfaceAccessor surfaceAccessor{*detector.geometry};
   extensions.surfaceAccessor
       .connect<&TestSourceLink::SurfaceAccessor::operator()>(&surfaceAccessor);
 
   const Gx2FitterOptions gx2fOptions(geoCtx, magCtx, calCtx, extensions,
                                      PropagatorPlainOptions(geoCtx, magCtx),
                                      rSurface, false, false,
                                      FreeToBoundCorrection(false), 6, 0);
 
   TrackContainer tracks{VectorTrackContainer{}, VectorMultiTrajectory{}};
 
   ACTS_DEBUG("Fit the track");
   ACTS_VERBOSE("startParameter unsmeared:\n" << parametersMeasurements);
   ACTS_VERBOSE("startParameter fit:\n" << startParametersFit);
   const auto res = fitter.fit(sourceLinks.begin(), sourceLinks.end(),
                               startParametersFit, gx2fOptions, tracks);
 
   BOOST_REQUIRE(!res.ok());
   BOOST_CHECK_EQUAL(res.error(),
                     GlobalChiSquareFitterError::NotEnoughMeasurements);
 
   ACTS_INFO("*** Test: NotEnoughMeasurements -- Finish");
 }
 
 BOOST_AUTO_TEST_CASE(FindHoles) {
   ACTS_INFO("*** Test: FindHoles -- Start");
 
   std::default_random_engine rng(42);
 
   ACTS_DEBUG("Create the detector");
   //  const std::size_t nSurfaces = 7;
   const std::size_t nSurfaces = 8;
   Detector detector;
   detector.geometry = makeToyDetector(geoCtx, nSurfaces);
 
   ACTS_DEBUG("Set the start parameters for measurement creation and fit");
   const auto parametersMeasurements = makeParameters();
   const auto startParametersFit = makeParameters(
       7_mm, 11_mm, 15_mm, 42_ns, 10_degree, 80_degree, 1_GeV, 1_e);
 
   ACTS_DEBUG("Create the measurements");
   using SimPropagator = Propagator<StraightLineStepper, Navigator>;
   const SimPropagator simPropagator = makeStraightPropagator(detector.geometry);
   const auto measurements =
       createMeasurements(simPropagator, geoCtx, magCtx, parametersMeasurements,
                          resMapAllPixel, rng);
   auto sourceLinks = prepareSourceLinks(measurements.sourceLinks);
   ACTS_VERBOSE("sourceLinks.size() [before] = " << sourceLinks.size());
 
   // We remove the first measurement in the list. This does not create a hole.
   sourceLinks.erase(std::next(sourceLinks.begin(), 0));
   ACTS_VERBOSE(
       "sourceLinks.size() [after first erase] = " << sourceLinks.size());
 
   // We remove the last measurement in the list. This does not create a hole.
   sourceLinks.pop_back();
   ACTS_VERBOSE("sourceLinks.size() [after pop] = " << sourceLinks.size());
 
   // We remove the second to last measurement in the list. This effectively
   // creates a hole on that surface.
   const std::size_t indexHole = sourceLinks.size() - 2;
   ACTS_VERBOSE("Remove measurement " << indexHole);
   sourceLinks.erase(std::next(sourceLinks.begin(), indexHole));
   ACTS_VERBOSE("sourceLinks.size() [after second-to-last erase]= "
                << sourceLinks.size());
 
   // We removed 3 measurements
   //  const std::size_t nMeasurements = nSurfaces - 2;
   const std::size_t nMeasurements = nSurfaces - 3;
   BOOST_REQUIRE_EQUAL(sourceLinks.size(), nMeasurements);
 
   ACTS_DEBUG("Set up the fitter");
   const Surface* rSurface = &parametersMeasurements.referenceSurface();
 
   using RecoStepper = EigenStepper<>;
   const auto recoPropagator =
       makeConstantFieldPropagator<RecoStepper>(detector.geometry, 0_T);
 
   using RecoPropagator = decltype(recoPropagator);
   using Gx2Fitter = Gx2Fitter<RecoPropagator, VectorMultiTrajectory>;
   const Gx2Fitter fitter(recoPropagator, gx2fLogger->clone());
 
   Gx2FitterExtensions<VectorMultiTrajectory> extensions;
   extensions.calibrator
       .connect<&testSourceLinkCalibratorStrict<VectorMultiTrajectory>>();
   TestSourceLink::SurfaceAccessor surfaceAccessor{*detector.geometry};
   extensions.surfaceAccessor
       .connect<&TestSourceLink::SurfaceAccessor::operator()>(&surfaceAccessor);
 
   const Gx2FitterOptions gx2fOptions(geoCtx, magCtx, calCtx, extensions,
                                      PropagatorPlainOptions(geoCtx, magCtx),
                                      rSurface, false, false,
                                      FreeToBoundCorrection(false), 20, 1e-5);
 
   TrackContainer tracks{VectorTrackContainer{}, VectorMultiTrajectory{}};
 
   ACTS_DEBUG("Fit the track");
   ACTS_VERBOSE("startParameter unsmeared:\n" << parametersMeasurements);
   ACTS_VERBOSE("startParameter fit:\n" << startParametersFit);
   const auto res = fitter.fit(sourceLinks.begin(), sourceLinks.end(),
                               startParametersFit, gx2fOptions, tracks);
 
   BOOST_REQUIRE(res.ok());
 
   const auto& track = *res;
 
   // -1, because the index starts at 0
   // -2, because the first and the last surface are not part of the track
   BOOST_CHECK_EQUAL(track.tipIndex(), nSurfaces - 1 - 2);
   BOOST_CHECK(track.hasReferenceSurface());
 
   // Track quantities
   CHECK_CLOSE_ABS(track.chi2(), 6.5, 2.);
   BOOST_CHECK_EQUAL(track.nDoF(), 10u);
   BOOST_CHECK_EQUAL(track.nHoles(), 1u);
   BOOST_CHECK_EQUAL(track.nMeasurements(), nMeasurements);
   BOOST_CHECK_EQUAL(track.nSharedHits(), 0u);
   BOOST_CHECK_EQUAL(track.nOutliers(), 0u);
 
   // Parameters
   // We need quite coarse checks here, since on different builds
   // the created measurements differ in the randomness
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc0], -11., 7e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc1], -15., 6e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundPhi], 1e-5, 1e3);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTheta], std::numbers::pi / 2,
                     1e-3);
   BOOST_CHECK_EQUAL(track.parameters()[eBoundQOverP], 1);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTime],
                     startParametersFit.parameters()[eBoundTime], 1e-6);
   BOOST_CHECK_CLOSE(track.covariance().determinant(), 4.7e-28, 2e0);
 
   ACTS_INFO("*** Test: FindHoles -- Finish");
 }
 
 BOOST_AUTO_TEST_CASE(Material) {
   ACTS_INFO("*** Test: Material -- Start");
 
   std::default_random_engine rng(42);
 
   ACTS_DEBUG("Create the detector");
   const std::size_t nSurfaces = 7;
   const std::set<std::size_t> surfaceIndexWithMaterial = {4};
   Detector detector;
   detector.geometry =
       makeToyDetector(geoCtx, nSurfaces, surfaceIndexWithMaterial);
 
   ACTS_DEBUG("Set the start parameters for measurement creation and fit");
   const auto parametersMeasurements = makeParameters();
   const auto startParametersFit = makeParameters(
       7_mm, 11_mm, 15_mm, 42_ns, 10_degree, 80_degree, 1_GeV, 1_e);
 
   ACTS_DEBUG("Create the measurements");
   using SimPropagator = Propagator<StraightLineStepper, Navigator>;
   const SimPropagator simPropagator = makeStraightPropagator(detector.geometry);
   auto measurements =
       createMeasurements(simPropagator, geoCtx, magCtx, parametersMeasurements,
                          resMapAllPixel, rng);
 
   const Vector2 scatterOffset = {100_mm, 100_mm};
   const std::size_t indexMaterialSurface = 3;
   for (std::size_t iMeas = indexMaterialSurface; iMeas < nSurfaces; iMeas++) {
     // This only works, because our detector is evenly spaced
     const std::size_t offsetFactor = iMeas - indexMaterialSurface;
 
     auto& sl = measurements.sourceLinks[iMeas];
     sl.parameters[0] += scatterOffset[0] * offsetFactor;
     sl.parameters[1] += scatterOffset[1] * offsetFactor;
   }
 
   const auto sourceLinks = prepareSourceLinks(measurements.sourceLinks);
   ACTS_VERBOSE("sourceLinks.size() = " << sourceLinks.size());
 
   BOOST_REQUIRE_EQUAL(sourceLinks.size(), nSurfaces);
 
   ACTS_DEBUG("Set up the fitter");
   const Surface* rSurface = &parametersMeasurements.referenceSurface();
 
   using RecoStepper = EigenStepper<>;
   const auto recoPropagator =
       makeConstantFieldPropagator<RecoStepper>(detector.geometry, 0_T);
 
   using RecoPropagator = decltype(recoPropagator);
   using Gx2Fitter = Gx2Fitter<RecoPropagator, VectorMultiTrajectory>;
   const Gx2Fitter fitter(recoPropagator, gx2fLogger->clone());
 
   Gx2FitterExtensions<VectorMultiTrajectory> extensions;
   extensions.calibrator
       .connect<&testSourceLinkCalibratorStrict<VectorMultiTrajectory>>();
   TestSourceLink::SurfaceAccessor surfaceAccessor{*detector.geometry};
   extensions.surfaceAccessor
       .connect<&TestSourceLink::SurfaceAccessor::operator()>(&surfaceAccessor);
 
   const Gx2FitterOptions gx2fOptions(geoCtx, magCtx, calCtx, extensions,
                                      PropagatorPlainOptions(geoCtx, magCtx),
                                      rSurface, true, false,
                                      FreeToBoundCorrection(false), 5, 0);
 
   TrackContainer tracks{VectorTrackContainer{}, VectorMultiTrajectory{}};
 
   ACTS_DEBUG("Fit the track");
   ACTS_VERBOSE("startParameter unsmeared:\n" << parametersMeasurements);
   ACTS_VERBOSE("startParameter fit:\n" << startParametersFit);
   const auto res = fitter.fit(sourceLinks.begin(), sourceLinks.end(),
                               startParametersFit, gx2fOptions, tracks);
 
   // Helper to visualise the detector
   {
     std::cout << "\n*** Create .obj of Detector ***\n" << std::endl;
     // Only need for obj
     ObjVisualization3D obj;
 
     bool triangulate = true;
     ViewConfig viewSensitive = {.color = {0, 180, 240}};
     viewSensitive.triangulate = triangulate;
     ViewConfig viewPassive = {.color = {240, 280, 0}};
     viewPassive.triangulate = triangulate;
     ViewConfig viewVolume = {.color = {220, 220, 0}};
     viewVolume.triangulate = triangulate;
     ViewConfig viewContainer = {.color = {220, 220, 0}};
     viewContainer.triangulate = triangulate;
     ViewConfig viewGrid = {.color = {220, 0, 0}};
     viewGrid.quarterSegments = 8;
     viewGrid.offset = 3.;
     viewGrid.triangulate = triangulate;
 
     std::string tag = "gx2f_toydet";
 
     const TrackingVolume& tgVolume =
         *(detector.geometry->highestTrackingVolume());
 
     GeometryView3D::drawTrackingVolume(obj, tgVolume, geoCtx, viewContainer,
                                        viewVolume, viewPassive, viewSensitive,
                                        viewGrid, true, tag);
   }
   // Helper to visualise the measurements
   {
     std::cout << "\n*** Create .obj of measurements ***\n" << std::endl;
     ObjVisualization3D obj;
 
     double localErrorScale = 10000000.;
     ViewConfig mcolor{.color = {255, 145, 48}};
     mcolor.offset = 2;
     //  mcolor.visible = true;
 
     drawMeasurements(obj, measurements, detector.geometry, geoCtx,
                      localErrorScale, mcolor);
 
     obj.write("meas");
   }
 
   BOOST_REQUIRE(res.ok());
 
   const auto& track = *res;
 
   BOOST_CHECK_EQUAL(track.tipIndex(), nSurfaces - 1);
   BOOST_CHECK(track.hasReferenceSurface());
 
   // TODO Add material handling to the gx2f, to pass the 6 commented tests
   // Track quantities
   //  CHECK_CLOSE_ABS(track.chi2(), 8., 2.);
   BOOST_CHECK_EQUAL(track.nDoF(), nSurfaces * 2);
   BOOST_CHECK_EQUAL(track.nHoles(), 0u);
   BOOST_CHECK_EQUAL(track.nMeasurements(), nSurfaces);
   BOOST_CHECK_EQUAL(track.nSharedHits(), 0u);
   BOOST_CHECK_EQUAL(track.nOutliers(), 0u);
 
   // Parameters
   // We need quite coarse checks here, since on different builds
   // the created measurements differ in the randomness
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc0], -11., 26e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundLoc1], -15., 15e0);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundPhi], 1e-5, 1.1e3);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTheta], std::numbers::pi / 2,
                     2e-2);
   BOOST_CHECK_EQUAL(track.parameters()[eBoundQOverP], 1);
   BOOST_CHECK_CLOSE(track.parameters()[eBoundTime],
                     startParametersFit.parameters()[eBoundTime], 1e-6);
   BOOST_CHECK_CLOSE(track.covariance().determinant(), 3.5e-27, 1e1);
 
   // Convergence
   BOOST_CHECK_EQUAL(
       (track.template component<
           std::uint32_t, hashString(Gx2fConstants::gx2fnUpdateColumn)>()),
       4);
 
   ACTS_INFO("*** Test: Material -- Finish");
 }
 BOOST_AUTO_TEST_SUITE_END()
 
 }  // namespace ActsTests

This page was automatically generated by the 2.3.7 LXR engine. The LXR team