EIC code displayed by LXR master/​acts/​Core/​include/​Acts/​EventData/​MultiTrajectory.hpp	Source navigation Diff markup Identifier search General search
	Version: master test Revert   Change

File indexing completed on 2026-05-21 07:46:23

 // 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/.
 
 #pragma once
 
 #include "Acts/Definitions/Algebra.hpp"
 #include "Acts/Definitions/TrackParametrization.hpp"
 #include "Acts/EventData/MeasurementHelpers.hpp"
 #include "Acts/EventData/SourceLink.hpp"
 #include "Acts/EventData/TrackStatePropMask.hpp"
 #include "Acts/EventData/TrackStateProxy.hpp"
 #include "Acts/EventData/Types.hpp"
 #include "Acts/Utilities/HashedString.hpp"
 #include "Acts/Utilities/ThrowAssert.hpp"
 
 #include <cstddef>
 #include <iterator>
 #include <memory>
 #include <optional>
 #include <string_view>
 #include <type_traits>
 
 #include <Eigen/Core>
 
 namespace Acts {
 
 // forward declarations
 template <typename derived_t>
 class MultiTrajectory;
 class Surface;
 
 namespace detail_anytstate {
 template <typename trajectory_t, bool read_only>
 class TrackStateHandler;
 }  // namespace detail_anytstate
 
 namespace detail_lt {
 
 /// Helper type that wraps two iterators
 template <bool reverse, typename trajectory_t, std::size_t M, bool ReadOnly>
 class TrackStateRange {
   using ProxyType = TrackStateProxy<trajectory_t, M, ReadOnly>;
   using IndexType = typename ProxyType::IndexType;
   static constexpr IndexType kInvalid = ProxyType::kInvalid;
 
  public:
   /// Iterator that wraps a track state proxy. The nullopt case signifies the
   /// end of the range, i.e. the "past-the-end" iterator
   struct Iterator {
     std::optional<ProxyType> proxy;
 
     using iterator_category = std::forward_iterator_tag;
     using value_type = ProxyType;
     using difference_type = std::ptrdiff_t;
     using pointer = void;
     using reference = void;
 
     Iterator& operator++() {
       if (!proxy) {
         return *this;
       }
       if constexpr (reverse) {
         if (proxy->hasPrevious()) {
           proxy = proxy->trajectory().getTrackState(proxy->previous());
           return *this;
         } else {
           proxy = std::nullopt;
           return *this;
         }
       } else {
         IndexType next =
             proxy->template component<IndexType, hashString("next")>();
         if (next != kInvalid) {
           proxy = proxy->trajectory().getTrackState(next);
           return *this;
         } else {
           proxy = std::nullopt;
           return *this;
         }
       }
     }
 
     Iterator operator++(int) {
       Iterator tmp(*this);
       operator++();
       return tmp;
     }
 
     bool operator==(const Iterator& other) const {
       if (!proxy && !other.proxy) {
         return true;
       }
       if (proxy && other.proxy) {
         return proxy->index() == other.proxy->index();
       }
       return false;
     }
 
     ProxyType operator*() const { return *proxy; }
     ProxyType operator*() { return *proxy; }
   };
 
   explicit TrackStateRange(ProxyType _begin) : m_begin{_begin} {}
   TrackStateRange() : m_begin{std::nullopt} {}
 
   Iterator begin() { return Iterator{m_begin}; }
   Iterator end() { return Iterator{std::nullopt}; }
 
   Iterator cbegin() const { return Iterator{m_begin}; }
   Iterator cend() const { return Iterator{std::nullopt}; }
 
  private:
   Iterator m_begin;
 };
 
 // implement track state visitor concept
 template <typename T, typename TS>
 concept VisitorConcept = requires(T& t, TS& ts) {
   { t(ts) } -> Concepts::same_as_any_of<void, bool>;
 };
 
 }  // namespace detail_lt
 
 /// This namespace contains typedefs and constant values that are used by
 /// other parts of the @c MultiTrajectory implementation. It extracts these
 /// from @c TrackStateTraits using the default maximum measurement dimension.
 /// @deprecated Use aliased types and constants directly
 namespace MultiTrajectoryTraits {
 
 /// Maximum number of measurement dimensions supported by trajectory
 /// @deprecated Use @ref Acts::kMeasurementSizeMax instead.
 [[deprecated("Use Acts::kMeasurementSizeMax instead.")]]
 constexpr unsigned int MeasurementSizeMax = kMeasurementSizeMax;
 
 /// Type alias for trajectory index type
 /// @deprecated Use @ref Acts::TrackIndexType instead.
 using IndexType [[deprecated("Use Acts::TrackIndexType instead.")]] =
     TrackIndexType;
 
 /// Invalid track state index constant
 /// @deprecated Use @ref Acts::kTrackIndexInvalid instead.
 [[deprecated("Use Acts::kTrackIndexInvalid instead.")]]
 constexpr IndexType kInvalid = kTrackIndexInvalid;
 
 }  // namespace MultiTrajectoryTraits
 
 template <typename T>
 struct IsReadOnlyMultiTrajectory;
 
 /// Store a trajectory of track states with multiple components.
 ///
 /// This container supports both simple, sequential trajectories as well
 /// as combinatorial or multi-component trajectories. Each point can store
 /// a parent point such that the trajectory forms a directed, acyclic graph
 /// of sub-trajectories. From a set of endpoints, all possible sub-components
 /// can be easily identified. Some functionality is provided to simplify
 /// iterating over specific sub-components.
 template <typename derived_t>
 class MultiTrajectory {
  public:
   /// Type alias for derived multi-trajectory implementation
   using Derived = derived_t;
 
   /// Flag indicating whether this multi-trajectory is read-only
   static constexpr bool ReadOnly = IsReadOnlyMultiTrajectory<Derived>::value;
 
   // Pull out type alias and re-expose them for ease of use.
   /// Maximum number of measurement dimensions supported by this trajectory
   static constexpr unsigned int MeasurementSizeMax = kMeasurementSizeMax;
 
   friend class TrackStateProxy<Derived, MeasurementSizeMax, true>;
   friend class TrackStateProxy<Derived, MeasurementSizeMax, false>;
   template <bool R>
   friend class AnyTrackStateProxy;
   template <typename T, bool R>
   friend class detail_anytstate::TrackStateHandler;
   template <typename T>
   friend class MultiTrajectory;
 
   /// Alias for the const version of a track state proxy, with the same
   /// backends as this container
   using ConstTrackStateProxy =
       Acts::TrackStateProxy<Derived, MeasurementSizeMax, true>;
 
   /// Alias for the mutable version of a track state proxy, with the same
   /// backends as this container
   using TrackStateProxy =
       Acts::TrackStateProxy<Derived, MeasurementSizeMax, false>;
 
   /// The index type of the track state container
   using IndexType = TrackIndexType;
 
   /// Sentinel value that indicates an invalid index
   static constexpr IndexType kInvalid = kTrackIndexInvalid;
 
  protected:
   MultiTrajectory() = default;  // pseudo abstract base class
 
  private:
   /// Helper to static cast this to the Derived class for CRTP
   constexpr Derived& self() { return static_cast<Derived&>(*this); }
   /// Helper to static cast this to the Derived class for CRTP. Const version.
   constexpr const Derived& self() const {
     return static_cast<const Derived&>(*this);
   }
 
   /// Helper function to check if a component exists IF it is an optional one.
   /// Used in assertions
   bool checkOptional(HashedString key, IndexType istate) const {
     using namespace Acts::HashedStringLiteral;
     switch (key) {
       case "predicted"_hash:
       case "filtered"_hash:
       case "smoothed"_hash:
       case "calibrated"_hash:
       case "jacobian"_hash:
       case "projector"_hash:
         return self().has_impl(key, istate);
       default:
         return true;
     }
   }
 
  public:
   /// @anchor track_state_container_track_access
   /// @name MultiTrajectory track state (proxy) access and manipulation
   ///
   /// These methods allow accessing track states, i.e. adding or retrieving a
   /// track state proxy that points at a specific track state in the container.
   ///
   /// @{
 
   /// Access a read-only point on the trajectory by index.
   /// @note Only available if the MultiTrajectory is not read-only
   /// @param istate The index to access
   /// @return Read only proxy to the stored track state
   ConstTrackStateProxy getTrackState(IndexType istate) const {
     return {*this, istate};
   }
 
   /// Access a writable point on the trajectory by index.
   /// @note Only available if the MultiTrajectory is not read-only
   /// @param istate The index to access
   /// @return Read-write proxy to the stored track state
   TrackStateProxy getTrackState(IndexType istate)
     requires(!ReadOnly)
   {
     return {*this, istate};
   }
 
   /// Add a track state without providing explicit information. Which components
   /// of the track state are initialized/allocated can be controlled via @p mask
   /// @note Only available if the MultiTrajectory is not read-only
   /// @param mask The bitmask that instructs which components to allocate and
   ///       which to leave invalid
   /// @param iprevious index of the previous state, kInvalid if first
   /// @return Index of the newly added track state
   IndexType addTrackState(TrackStatePropMask mask = TrackStatePropMask::All,
                           IndexType iprevious = kInvalid)
     requires(!ReadOnly)
   {
     return self().addTrackState_impl(mask, iprevious);
   }
 
   /// Add a track state to the container and return a track state proxy to it
   /// This effectively calls @c addTrackState and @c getTrackState
   /// @note Only available if the track state container is not read-only
   /// @param mask Mask indicating which track state components to allocate
   /// @param iprevious Index of the previous track state for linking
   /// @return a track state proxy to the newly added track state
   TrackStateProxy makeTrackState(
       TrackStatePropMask mask = TrackStatePropMask::All,
       IndexType iprevious = kInvalid)
     requires(!ReadOnly)
   {
     return getTrackState(addTrackState(mask, iprevious));
   }
 
   /// @}
 
   /// @anchor track_state_container_iteration
   /// @name MultiTrajectory track state iteration
   /// @{
 
   /// Visit all previous states starting at a given endpoint.
   ///
   /// @param iendpoint  index of the last state
   /// @param callable   non-modifying functor to be called with each point
   template <typename F>
   void visitBackwards(IndexType iendpoint, F&& callable) const
     requires detail_lt::VisitorConcept<F, ConstTrackStateProxy>;
 
   /// Apply a function to all previous states starting at a given endpoint.
   ///
   /// @param iendpoint  index of the last state
   /// @param callable   modifying functor to be called with each point
   ///
   /// @warning If the trajectory contains multiple components with common
   ///          points, this can have an impact on the other components.
   /// @note Only available if the MultiTrajectory is not read-only
   template <typename F>
   void applyBackwards(IndexType iendpoint, F&& callable)
     requires(!ReadOnly) && detail_lt::VisitorConcept<F, TrackStateProxy>
   {
     if (iendpoint == kInvalid) {
       throw std::runtime_error(
           "Cannot apply backwards with kInvalid as endpoint");
     }
 
     while (true) {
       auto ts = getTrackState(iendpoint);
       if constexpr (std::is_same_v<std::invoke_result_t<F, TrackStateProxy>,
                                    bool>) {
         bool proceed = callable(ts);
         // this point has no parent and ends the trajectory, or a break was
         // requested
         if (!proceed || !ts.hasPrevious()) {
           break;
         }
       } else {
         callable(ts);
         // this point has no parent and ends the trajectory
         if (!ts.hasPrevious()) {
           break;
         }
       }
       iendpoint = ts.previous();
     }
   }
 
   /// Range for the track states from @p iendpoint to the trajectory start
   /// @param iendpoint Trajectory entry point to start from
   /// @return Iterator pair to iterate over
   /// @note Const version
   auto reverseTrackStateRange(IndexType iendpoint) const {
     using range_t =
         detail_lt::TrackStateRange<true, Derived, MeasurementSizeMax, true>;
     if (iendpoint == kInvalid) {
       return range_t{};
     }
 
     return range_t{getTrackState(iendpoint)};
   }
 
   /// Range for the track states from @p iendpoint to the trajectory start,
   /// i.e from the outside in.
   /// @note Only available if the MultiTrajectory is not read-only
   /// @param iendpoint Trajectory entry point to start from
   /// @return Iterator pair to iterate over
   /// @note Mutable version
   auto reverseTrackStateRange(IndexType iendpoint)
     requires(!ReadOnly)
   {
     using range_t =
         detail_lt::TrackStateRange<true, Derived, MeasurementSizeMax, false>;
     if (iendpoint == kInvalid) {
       return range_t{};
     }
 
     return range_t{getTrackState(iendpoint)};
   }
 
   /// Range for the track states from @p istartpoint to the trajectory end,
   /// i.e from inside out
   /// @param istartpoint Trajectory state index for the innermost track
   ///        state to start from
   /// @return Iterator pair to iterate over
   /// @note Const version
   auto forwardTrackStateRange(IndexType istartpoint) const {
     using range_t =
         detail_lt::TrackStateRange<false, Derived, MeasurementSizeMax, true>;
     if (istartpoint == kInvalid) {
       return range_t{};
     }
 
     return range_t{getTrackState(istartpoint)};
   }
 
   /// Range for the track states from @p istartpoint to the trajectory end,
   /// i.e from inside out
   /// @note Only available if the MultiTrajectory is not read-only
   /// @param istartpoint Trajectory state index for the innermost track
   ///        state to start from
   /// @return Iterator pair to iterate over
   auto forwardTrackStateRange(IndexType istartpoint)
     requires(!ReadOnly)
   {
     using range_t =
         detail_lt::TrackStateRange<false, Derived, MeasurementSizeMax, false>;
     if (istartpoint == kInvalid) {
       return range_t{};
     }
 
     return range_t{getTrackState(istartpoint)};
   }
 
   /// @}
 
   /// @anchor track_state_container_columns
   /// @name MultiTrajectory column management
   /// MultiTrajectory can manage a set of common static columns, and dynamic
   /// columns that can be added at runtime. This set of methods allows you to
   /// manage the dynamic columns.
   /// @{
 
   /// Add a column to the @c MultiTrajectory
   /// @tparam T Type of the column values to add
   /// @param key the name of the column to be added
   /// @note This takes a string argument rather than a hashed string to maintain
   ///       compatibility with backends.
   /// @note Only available if the MultiTrajectory is not read-only
   template <typename T>
   void addColumn(std::string_view key)
     requires(!ReadOnly)
   {
     self().template addColumn_impl<T>(key);
   }
 
   /// Check if a column with a key @p key exists.
   /// @param key Key to check for a column with
   /// @return True if the column exists, false if not.
   bool hasColumn(HashedString key) const { return self().hasColumn_impl(key); }
 
   /// @}
 
   /// Clear the @c MultiTrajectory. Leaves the underlying storage untouched
   /// @note Only available if the MultiTrajectory is not read-only
   void clear()
     requires(!ReadOnly)
   {
     self().clear_impl();
   }
 
   /// Returns the number of track states contained
   /// @return The number of track states
   IndexType size() const { return self().size_impl(); }
 
  protected:
   // These are internal helper functions which the @c TrackStateProxy class talks to
 
   /// Check for component existence of @p key in track satet @p istate
   /// @param key The key for which to check
   /// @param istate The track state index to check
   /// @return True if the component exists, false if not
   bool has(HashedString key, IndexType istate) const {
     return self().has_impl(key, istate);
   }
 
   /// Check for component existence of @p key in track satet @p istate
   /// @tparam key The key for which to check
   /// @param istate The track state index to check
   /// @return True if the component exists, false if not
   template <HashedString key>
   bool has(IndexType istate) const {
     return self().has_impl(key, istate);
   }
 
   /// Get parameters for a track state
   /// @param parIdx The parameter index
   /// @return Parameters vector
   typename TrackStateProxy::Parameters parameters(IndexType parIdx)
     requires(!ReadOnly)
   {
     return self().parameters_impl(parIdx);
   }
 
   /// Get parameters for a track state (const)
   /// @param parIdx The parameter index
   /// @return Const parameters vector
   typename ConstTrackStateProxy::ConstParameters parameters(
       IndexType parIdx) const {
     return self().parameters_impl(parIdx);
   }
 
   /// Get covariance for a track state
   /// @param covIdx The covariance index
   /// @return Covariance matrix
   typename TrackStateProxy::Covariance covariance(IndexType covIdx)
     requires(!ReadOnly)
   {
     return self().covariance_impl(covIdx);
   }
 
   /// Get covariance for a track state (const)
   /// @param covIdx The covariance index
   /// @return Const covariance matrix
   typename ConstTrackStateProxy::ConstCovariance covariance(
       IndexType covIdx) const {
     return self().covariance_impl(covIdx);
   }
 
   /// Retrieve a jacobian proxy instance for a jacobian at a given index
   /// @param istate The track state
   /// @return Mutable proxy
   typename TrackStateProxy::Jacobian jacobian(IndexType istate)
     requires(!ReadOnly)
   {
     return self().jacobian_impl(istate);
   }
 
   /// Retrieve a jacobian proxy instance for a jacobian at a given index
   /// @param istate The track state
   /// @return Const proxy
   typename ConstTrackStateProxy::ConstJacobian jacobian(
       IndexType istate) const {
     return self().jacobian_impl(istate);
   }
 
   /// Retrieve a calibrated measurement proxy instance for a measurement at a
   /// given index
   /// @tparam measdim the measurement dimension
   /// @param istate The track state
   /// @return Mutable proxy
   template <std::size_t measdim>
   typename TrackStateProxy::template Calibrated<measdim> calibrated(
       IndexType istate)
     requires(!ReadOnly)
   {
     return self().template calibrated_impl<measdim>(istate);
   }
 
   /// Retrieve a calibrated measurement proxy instance for a measurement at a
   /// given index
   /// @tparam measdim the measurement dimension
   /// @param istate The track state
   /// @return Const proxy
   template <std::size_t measdim>
   typename ConstTrackStateProxy::template ConstCalibrated<measdim> calibrated(
       IndexType istate) const {
     return self().template calibrated_impl<measdim>(istate);
   }
 
   /// Retrieve a calibrated measurement covariance proxy instance for a
   /// measurement at a given index
   /// @tparam measdim the measurement dimension
   /// @param istate The track state
   /// @return Mutable proxy
   template <std::size_t measdim>
   typename TrackStateProxy::template CalibratedCovariance<measdim>
   calibratedCovariance(IndexType istate)
     requires(!ReadOnly)
   {
     return self().template calibratedCovariance_impl<measdim>(istate);
   }
 
   /// Retrieve a calibrated measurement covariance proxy instance for a
   /// measurement at a given index
   /// @param istate The track state
   /// @return Mutable proxy
   typename TrackStateProxy::EffectiveCalibrated effectiveCalibrated(
       IndexType istate)
     requires(!ReadOnly)
   {
     // This abuses an incorrectly sized vector / matrix to access the
     // data pointer! This works (don't use the matrix as is!), but be
     // careful!
     return typename TrackStateProxy::EffectiveCalibrated{
         calibrated<eBoundSize>(istate).data(), calibratedSize(istate)};
   }
 
   /// Retrieve a calibrated measurement covariance proxy instance for a
   /// measurement at a given index
   /// @param istate The track state
   /// @return Const proxy
   typename ConstTrackStateProxy::EffectiveCalibrated effectiveCalibrated(
       IndexType istate) const {
     // This abuses an incorrectly sized vector / matrix to access the
     // data pointer! This works (don't use the matrix as is!), but be
     // careful!
     return typename ConstTrackStateProxy::EffectiveCalibrated{
         calibrated<eBoundSize>(istate).data(), calibratedSize(istate)};
   }
 
   /// Retrieve a calibrated measurement covariance proxy instance for a
   /// measurement at a given index
   /// @param istate The track state
   /// @return Mutable proxy
   typename TrackStateProxy::EffectiveCalibratedCovariance
   effectiveCalibratedCovariance(IndexType istate)
     requires(!ReadOnly)
   {
     // This abuses an incorrectly sized vector / matrix to access the
     // data pointer! This works (don't use the matrix as is!), but be
     // careful!
     return typename TrackStateProxy::EffectiveCalibratedCovariance{
         calibratedCovariance<eBoundSize>(istate).data(), calibratedSize(istate),
         calibratedSize(istate)};
   }
 
   /// Retrieve a calibrated measurement covariance proxy instance for a
   /// measurement at a given index
   /// @param istate The track state
   /// @return Const proxy
   typename ConstTrackStateProxy::EffectiveCalibratedCovariance
   effectiveCalibratedCovariance(IndexType istate) const {
     // This abuses an incorrectly sized vector / matrix to access the
     // data pointer! This works (don't use the matrix as is!), but be
     // careful!
     return typename ConstTrackStateProxy::EffectiveCalibratedCovariance{
         calibratedCovariance<eBoundSize>(istate).data(), calibratedSize(istate),
         calibratedSize(istate)};
   }
 
   /// Retrieve a calibrated measurement covariance proxy instance for a
   /// measurement at a given index
   /// @param istate The track state
   /// @return Const proxy
   template <std::size_t measdim>
   typename ConstTrackStateProxy::template ConstCalibratedCovariance<measdim>
   calibratedCovariance(IndexType istate) const {
     return self().template calibratedCovariance_impl<measdim>(istate);
   }
 
   /// Get the calibrated measurement size for a track state
   /// @param istate The track state
   /// @return the calibrated size
   IndexType calibratedSize(IndexType istate) const {
     return self().calibratedSize_impl(istate);
   }
 
   /// Share a shareable component from between track state.
   /// @param iself The track state index to share "into"
   /// @param iother The track state index to share from
   /// @param shareSource Which component to share from
   /// @param shareTarget Which component to share as. This doesn't have to be the same
   ///                    as @p shareSource, e.g. predicted can be shared as filtered.
   /// @note Shareable components are predicted, filtered, smoothed, calibrated, jacobian,
   ///       or projector. See @c TrackStatePropMask.
   /// @note The track states both need to be stored in the
   ///       same @c MultiTrajectory instance
   void shareFrom(IndexType iself, IndexType iother,
                  TrackStatePropMask shareSource, TrackStatePropMask shareTarget)
     requires(!ReadOnly)
   {
     self().shareFrom_impl(iself, iother, shareSource, shareTarget);
   }
 
   /// Unset an optional track state component
   /// @param target The component to unset
   /// @param istate The track state index to operate on
   void unset(TrackStatePropMask target, IndexType istate)
     requires(!ReadOnly)
   {
     self().unset_impl(target, istate);
   }
 
   /// Add additional components to an existing track state
   /// @note Only available if the track state container is not read-only
   /// @param istate The track state index to alter
   /// @param mask The bitmask that instructs which components to allocate
   void addTrackStateComponents(IndexType istate, TrackStatePropMask mask)
     requires(!ReadOnly)
   {
     self().addTrackStateComponents_impl(istate, mask);
   }
 
   /// Retrieve a mutable reference to a component
   /// @tparam T The type of the component to access
   /// @tparam key String key for the component to access
   /// @param istate The track state index to operate on
   /// @return Mutable reference to the component given by @p key
   template <typename T, HashedString key>
   T& component(IndexType istate)
     requires(!ReadOnly)
   {
     assert(checkOptional(key, istate));
     return *std::any_cast<T*>(self().component_impl(key, istate));
   }
 
   /// Retrieve a mutable reference to a component
   /// @tparam T The type of the component to access
   /// @param key String key for the component to access
   /// @param istate The track state index to operate on
   /// @return Mutable reference to the component given by @p key
   template <typename T>
   T& component(HashedString key, IndexType istate)
     requires(!ReadOnly)
   {
     assert(checkOptional(key, istate));
     return *std::any_cast<T*>(self().component_impl(key, istate));
   }
 
   /// Retrieve a const reference to a component
   /// @tparam T The type of the component to access
   /// @tparam key String key for the component to access
   /// @param istate The track state index to operate on
   /// @return Const reference to the component given by @p key
   template <typename T, HashedString key>
   const T& component(IndexType istate) const {
     assert(checkOptional(key, istate));
     return *std::any_cast<const T*>(self().component_impl(key, istate));
   }
 
   /// Retrieve a const reference to a component
   /// @tparam T The type of the component to access
   /// @param key String key for the component to access
   /// @param istate The track state index to operate on
   /// @return Const reference to the component given by @p key
   template <typename T>
   const T& component(HashedString key, IndexType istate) const {
     assert(checkOptional(key, istate));
     return *std::any_cast<const T*>(self().component_impl(key, istate));
   }
 
   /// Allocate storage for a calibrated measurement of specified dimension
   /// @param istate The track state to store for
   /// @param measdim the dimension of the measurement to store
   /// @note In case an allocation is already present, no additional allocation
   ///       will be performed, but the existing allocation will be zeroed.
   void allocateCalibrated(IndexType istate, std::size_t measdim) {
     throw_assert(measdim > 0 && measdim <= eBoundSize,
                  "Invalid measurement dimension detected");
 
     visit_measurement(measdim, [this, istate]<std::size_t DIM>(
                                    std::integral_constant<std::size_t, DIM>) {
       self().allocateCalibrated_impl(
           istate, Vector<DIM>{Vector<DIM>::Zero()},
           SquareMatrix<DIM>{SquareMatrix<DIM>::Zero()});
     });
   }
 
   /// Allocate storage for calibrated measurement
   /// @tparam measdim Measurement dimension
   /// @tparam val_t Value type
   /// @tparam cov_t Covariance type
   /// @param istate State index
   /// @param val Measurement values
   /// @param cov Measurement covariance
   template <std::size_t measdim, typename val_t, typename cov_t>
   void allocateCalibrated(IndexType istate, const Eigen::DenseBase<val_t>& val,
                           const Eigen::DenseBase<cov_t>& cov) {
     self().allocateCalibrated_impl(istate, val, cov);
   }
 
   /// Set the uncalibrated source link for a track state
   /// @param istate State index
   /// @param sourceLink Source link to set
   void setUncalibratedSourceLink(IndexType istate, SourceLink&& sourceLink)
     requires(!ReadOnly)
   {
     self().setUncalibratedSourceLink_impl(istate, std::move(sourceLink));
   }
 
   /// Get the uncalibrated source link for a track state
   /// @param istate State index
   /// @return Source link for the specified state
   SourceLink getUncalibratedSourceLink(IndexType istate) const {
     return self().getUncalibratedSourceLink_impl(istate);
   }
 
   /// Get the reference surface for a track state
   /// @param istate State index
   /// @return Pointer to the reference surface
   const Surface* referenceSurface(IndexType istate) const {
     return self().referenceSurface_impl(istate);
   }
 
   /// Set the reference surface for a track state
   /// @param istate State index
   /// @param surface Shared pointer to the reference surface
   void setReferenceSurface(IndexType istate,
                            std::shared_ptr<const Surface> surface)
     requires(!ReadOnly)
   {
     self().setReferenceSurface_impl(istate, std::move(surface));
   }
 
  private:
   template <typename T>
   void copyDynamicFrom(IndexType dstIdx, const T& src, IndexType srcIdx)
     requires(!ReadOnly)
   {
     const auto& dynamicKeys = src.self().dynamicKeys_impl();
     for (const auto key : dynamicKeys) {
       std::any srcPtr = src.self().component_impl(key, srcIdx);
       self().copyDynamicFrom_impl(dstIdx, key, srcPtr);
     }
   }
 };
 
 }  // namespace Acts
 
 #include "Acts/EventData/MultiTrajectory.ipp"

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