feat(geo): Gen3 portal links

Core/include/Acts/Geometry/CompositePortalLink.hpp

@@ -0,0 +1,99 @@
+// This file is part of the Acts project.
+//
+// Copyright (C) 2024 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 http://mozilla.org/MPL/2.0/.
+
+#pragma once
+
+#include "Acts/Definitions/Tolerance.hpp"
+#include "Acts/Geometry/PortalLinkBase.hpp"
+
+#include <iosfwd>
+
+#include <boost/container/small_vector.hpp>
+
+namespace Acts {
+
+/// Composite portal links can graft together other portal link instances, for
+/// example grids that could not be merged due to invalid binnings.
+///
+/// +-------+      +-------+
+/// |       |      |       |
+/// |       |      |       |
+/// |       |      |       |
+/// +-------+      |       |
+/// |       |      |       |
+/// |       |  +   +-------+
+/// |       |      |       |
+/// +-------+      |       |
+/// |       |      |       |
+/// |       |      +-------+
+/// |       |      |       |
+/// +-------+      +-------+
+///
+/// During resolution, it will consult each of it's children and return result
+/// of the resolution on the first surface where the lookup position is within
+/// bounds.
+class CompositePortalLink final : public PortalLinkBase {
+ public:
+  /// Construct a composite portal from two arbitrary other portal link. The
+  /// only requirement is that the portal link surfaces are mergeable.
+  /// @param a The first portal link
+  /// @param b The second portal link
+  /// @param direction The binning direction
+  /// @param flatten If true, the composite will flatten any nested composite
+  CompositePortalLink(std::unique_ptr<PortalLinkBase> a,
+                      std::unique_ptr<PortalLinkBase> b, BinningValue direction,
+                      bool flatten = true);
+
+  /// Print the composite portal link
+  /// @param os The output stream
+  void toStream(std::ostream& os) const final;
+
+  /// Resolve the volume for a 2D position
+  /// @note This will transform the position to global coordinates before
+  ///       consulting its children.
+  /// @note @p position is assumed to be on surface
+  /// @param gctx The geometry context
+  /// @param position The 2D position
+  /// @param tolerance The tolerance
+  Result<const TrackingVolume*> resolveVolume(
+      const GeometryContext& gctx, const Vector2& position,
+      double tolerance = s_onSurfaceTolerance) const final;
+
+  /// Resolve the volume for a 3D position
+  /// @note @p position is assumed to be on surface
+  /// @param gctx The geometry context
+  /// @param position The 3D position
+  /// @param tolerance The tolerance
+  Result<const TrackingVolume*> resolveVolume(
+      const GeometryContext& gctx, const Vector3& position,
+      double tolerance = s_onSurfaceTolerance) const final;
+
+  /// Get the depth of the composite tree
+  /// @return The depth
+  std::size_t depth() const;
+
+  /// Get the number of children
+  /// @return The number of children
+  std::size_t size() const;
+
+ private:
+  /// Helper function to construct a merged surface from two portal links along
+  /// a given direction
+  /// @param a The first portal link
+  /// @param b The second portal link
+  /// @param direction The merging direction
+  /// @return The merged surface
+  static std::shared_ptr<RegularSurface> mergedSurface(const PortalLinkBase* a,
+                                                       const PortalLinkBase* b,
+                                                       BinningValue direction);
+
+  boost::container::small_vector<std::unique_ptr<PortalLinkBase>, 4>
+      m_children{};
+};
+
+}  // namespace Acts