EIC code displayed by LXR master/​acts/​Core/​src/​Navigation/​NavigationStream.cpp	Source navigation Diff markup Identifier search General search
	Version: master test Revert   Change

File indexing completed on 2026-10-28 08:19:31

 // 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 "Acts/Navigation/NavigationStream.hpp"
 
 #include "Acts/Propagator/NavigationTarget.hpp"
 #include "Acts/Surfaces/BoundaryTolerance.hpp"
 #include "Acts/Surfaces/Surface.hpp"
 #include "Acts/Utilities/Enumerate.hpp"
 #include "Acts/Utilities/StringHelpers.hpp"
 
 #include <algorithm>
 
 namespace Acts {
 
 bool NavigationStream::initialize(const GeometryContext& gctx,
                                   const QueryPoint& queryPoint,
                                   const Logger& logger,
                                   const double onSurfaceTolerance,
                                   const bool candidatesAreUnique) {
   // Position and direction from the query point
   const Vector3& position = queryPoint.position;
   const Vector3& direction = queryPoint.direction;
 
   // De-duplicate by surface pointer. The first entry wins, so its tolerance is
   // the one used.
   if (!candidatesAreUnique) {
     ACTS_VERBOSE("De-duplicate the candidates:" << m_candidates);
     std::size_t writeIdx = 0;
     for (std::size_t readIdx = 0; readIdx < m_candidates.size(); ++readIdx) {
       const Surface* surface = &m_candidates[readIdx].surface();
       bool alreadySeen = false;
       for (std::size_t k = 0; k < writeIdx; ++k) {
         if (&m_candidates[k].surface() == surface) {
           alreadySeen = true;
           break;
         }
       }
       if (!alreadySeen) {
         if (writeIdx != readIdx) {
           m_candidates[writeIdx] = m_candidates[readIdx];
         }
         ++writeIdx;
       }
     }
     m_candidates.erase(m_candidates.begin() + writeIdx, m_candidates.end());
   }
 
   // Collect additional candidates for the second valid intersection. Reuse the
   // member scratch buffer to avoid a heap allocation on every call.
   std::vector<NavigationTarget>& additionalCandidates = m_additionalCandidates;
   additionalCandidates.clear();
   for (auto& candidate : m_candidates) {
     // Get the surface from the object intersection
     const Surface& surface = candidate.surface();
     // Intersect the surface
     auto multiIntersection =
         surface.intersect(gctx, position, direction,
                           candidate.boundaryTolerance(), onSurfaceTolerance);
 
     bool firstValid = multiIntersection.at(0).isValid();
     bool secondValid = multiIntersection.at(1).isValid();
     if (firstValid && !secondValid) {
       if (multiIntersection.at(0).pathLength() < -onSurfaceTolerance) {
         continue;
       }
       candidate.intersection() = multiIntersection.at(0);
       candidate.intersectionIndex() = 0;
     } else if (!firstValid && secondValid) {
       if (multiIntersection.at(1).pathLength() < -onSurfaceTolerance) {
         continue;
       }
       candidate.intersection() = multiIntersection.at(1);
       candidate.intersectionIndex() = 1;
     } else {
       // Split them into valid intersections, keep track of potentially
       // additional candidates
       bool originalCandidateUpdated = false;
       for (auto [intersectionIndex, intersection] :
            enumerate(multiIntersection)) {
         // Skip negative solutions, respecting the on surface tolerance
         if (intersection.pathLength() < -onSurfaceTolerance) {
           continue;
         }
         // Valid solution is either on surface or updates the distance
         if (intersection.isValid()) {
           if (!originalCandidateUpdated) {
             candidate.intersection() = intersection;
             candidate.intersectionIndex() = intersectionIndex;
             originalCandidateUpdated = true;
           } else {
             NavigationTarget additionalCandidate = candidate;
             additionalCandidate.intersection() = intersection;
             additionalCandidate.intersectionIndex() = intersectionIndex;
             additionalCandidates.emplace_back(additionalCandidate);
           }
         }
       }
     }
   }
 
   // Drop the candidates without a valid intersection before sorting,
   // would sort last and be cut off anyway...
   std::erase_if(m_candidates, [](const NavigationTarget& a) {
     return !a.intersection().isValid();
   });
 
   // Append the multi intersection candidates (all valid)
   m_candidates.insert(m_candidates.end(), additionalCandidates.begin(),
                       additionalCandidates.end());
 
   // Sort the candidates by path length
   std::ranges::sort(m_candidates, NavigationTarget::pathLengthOrder);
   ACTS_VERBOSE("Sorted candidates:" << m_candidates);
 
   // Remove duplicate intersections, preserving distinct crossings of the same
   // surface. In particular, an on-surface entry must not discard a later exit.
   auto nonUniqueRange = std::ranges::unique(
       m_candidates.begin(), m_candidates.end(),
       [](const NavigationTarget& a, const NavigationTarget& b) {
         return &a.surface() == &b.surface() &&
                a.intersectionIndex() == b.intersectionIndex();
       });
   m_candidates.erase(nonUniqueRange.begin(), nonUniqueRange.end());
 
   m_currentIndex = 0;
   return isValid();
 }
 
 bool NavigationStream::update(const GeometryContext& gctx,
                               const QueryPoint& queryPoint,
                               const Logger& logger, double onSurfaceTolerance) {
   ACTS_VERBOSE("Update from position " << toString(queryPoint.position)
                                        << " and direction "
                                        << toString(queryPoint.direction));
   // Loop over the (currently valid) candidates and update
   for (; m_currentIndex < m_candidates.size(); ++m_currentIndex) {
     // Get the candidate, and resolve the tuple
     NavigationTarget& candidate = currentCandidate();
     // Get the surface from the object intersection
     const Surface& surface = candidate.surface();
     // (re-)Intersect the surface
     auto multiIntersection =
         surface.intersect(gctx, queryPoint.position, queryPoint.direction,
                           candidate.boundaryTolerance(), onSurfaceTolerance);
     // Split them into valid intersections
     for (auto [intersectionIndex, intersection] :
          enumerate(multiIntersection)) {
       // Skip wrong index solution
       if (intersectionIndex != candidate.intersectionIndex()) {
         continue;
       }
       // Valid solution is either on surface or updates the distance
       if (intersection.isValid()) {
         candidate.intersection() = intersection;
         ACTS_VERBOSE("Updated candidate " << candidate);
         return true;
       }
     }
   }
   // No candidate was reachable
   return false;
 }
 
 void NavigationStream::reset() {
   m_candidates.clear();
   m_currentIndex = 0;
 }
 
 void NavigationStream::addSurfaceCandidate(
     const Surface& surface, const BoundaryTolerance& bTolerance) {
   m_candidates.emplace_back(Intersection3D::Invalid(), 0, surface, bTolerance);
 }
 
 void NavigationStream::addSurfaceCandidates(
     std::span<const Surface*> surfaces, const BoundaryTolerance& bTolerance) {
   m_candidates.reserve(m_candidates.size() + surfaces.size());
   std::ranges::for_each(surfaces, [&](const Surface* surface) {
     m_candidates.emplace_back(Intersection3D::Invalid(), 0, *surface,
                               bTolerance);
   });
 }
 
 void NavigationStream::addPortalCandidate(const Portal& portal) {
   m_candidates.emplace_back(Intersection3D::Invalid(), 0, portal,
                             BoundaryTolerance::None());
 }
 
 AppendOnlyNavigationStream::AppendOnlyNavigationStream(NavigationStream& stream)
     : m_stream{&stream} {}
 
 void AppendOnlyNavigationStream::addPortalCandidate(const Portal& portal) {
   m_stream->addPortalCandidate(portal);
 }
 
 void AppendOnlyNavigationStream::addSurfaceCandidate(
     const Surface& surface, const BoundaryTolerance& bTolerance) {
   m_stream->addSurfaceCandidate(surface, bTolerance);
 }
 
 }  // namespace Acts

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