GenericApproachDescriptor.cpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file GenericApproachDescriptor.cpp

// This file is part of the Acts project.
//
// Copyright (C) 2018 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/.

#include "Acts/Geometry/GenericApproachDescriptor.hpp"

#include "Acts/Surfaces/Surface.hpp"
#include "Acts/Utilities/Intersection.hpp"

#include <algorithm>

#include <boost/container/small_vector.hpp>

void Acts::GenericApproachDescriptor::registerLayer(const Layer& lay) {
  // go through the surfaces
  for (auto& sf : m_surfaceCache) {
    auto mutableSf = const_cast<Surface*>(sf);
    mutableSf->associateLayer(lay);
  }
}

Acts::SurfaceIntersection Acts::GenericApproachDescriptor::approachSurface(
    const GeometryContext& gctx, const Vector3& position,
    const Vector3& direction, const BoundaryCheck& bcheck, double pLimit,
    double oLimit, double tolerance) const {
  // almost always 2
  boost::container::small_vector<SurfaceIntersection, 2> sIntersections;
  sIntersections.reserve(m_surfaceCache.size());
  for (const auto& sf : m_surfaceCache) {
    auto sfIntersection = sf->intersect(gctx, position, direction, bcheck);
    for (const auto& intersection : sfIntersection.split()) {
      if (intersection &&
          detail::checkIntersection(intersection, pLimit, oLimit, tolerance)) {
        sIntersections.push_back(intersection);
      }
    }
  }
  if (sIntersections.empty()) {
    return SurfaceIntersection::invalid();
  }
  return *std::min_element(sIntersections.begin(), sIntersections.end(),
                           SurfaceIntersection::forwardOrder);
}

const std::vector<const Acts::Surface*>&
Acts::GenericApproachDescriptor::containedSurfaces() const {
  return m_surfaceCache;
}

std::vector<const Acts::Surface*>&
Acts::GenericApproachDescriptor::containedSurfaces() {
  return m_surfaceCache;
}