PortalHelper.cpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2023 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/Detector/detail/PortalHelper.hpp"

#include "Acts/Detector/Portal.hpp"
#include "Acts/Navigation/DetectorVolumeUpdators.hpp"
#include "Acts/Navigation/NavigationDelegates.hpp"
#include "Acts/Surfaces/Surface.hpp"
#include "Acts/Utilities/Helpers.hpp"

#include <array>
#include <stdexcept>
#include <utility>

void Acts::Experimental::detail::PortalHelper::attachDetectorVolumeUpdator(
    Portal& portal, const std::shared_ptr<DetectorVolume>& volume,
    const Direction& direction) {
  // Create a shared link instance & delegate
  auto volumeLinkImpl =
      std::make_unique<const Acts::Experimental::SingleDetectorVolumeImpl>(
          volume.get());
  Acts::Experimental::DetectorVolumeUpdator volumeLink;
  volumeLink.connect<&Acts::Experimental::SingleDetectorVolumeImpl::update>(
      std::move(volumeLinkImpl));
  // Update the volume link and the store
  portal.assignDetectorVolumeUpdator(direction, std::move(volumeLink),
                                     {volume});
}

void Acts::Experimental::detail::PortalHelper::attachDetectorVolumesUpdator(
    const GeometryContext& gctx, Portal& portal,
    const std::vector<std::shared_ptr<DetectorVolume>>& volumes,
    const Direction& direction, const std::vector<ActsScalar>& boundaries,
    const BinningValue& binning) {
  // Check if the boundaries need a transform
  const auto pTransform = portal.surface().transform(gctx);
  // Creating a link to the mother
  auto volumes1D = std::make_unique<const BoundVolumesGrid1Impl>(
      boundaries, binning, unpack_shared_const_vector(volumes),
      pTransform.inverse());
  DetectorVolumeUpdator dVolumeUpdator;
  dVolumeUpdator.connect<&BoundVolumesGrid1Impl::update>(std::move(volumes1D));
  portal.assignDetectorVolumeUpdator(direction, std::move(dVolumeUpdator),
                                     volumes);
}

void Acts::Experimental::detail::PortalHelper::attachDetectorVolumeUpdators(
    const GeometryContext& gctx,
    const std::vector<std::shared_ptr<DetectorVolume>>& volumes,
    std::vector<PortalReplacement>& pReplacements) {
  // Unpack to navigation bare points
  auto cVolumes = unpack_shared_const_vector(volumes);
  // Set to the constructed portals (p), at index (i), in direction (d)
  // using boundaries and binning
  for (auto& [p, i, dir, boundaries, binning] : pReplacements) {
    // Check if the boundaries need a transform
    const auto pTransform = p->surface().transform(gctx);
    // Creating a link to the mother
    auto volumes1D = std::make_unique<const BoundVolumesGrid1Impl>(
        boundaries, binning, cVolumes, pTransform.inverse());
    DetectorVolumeUpdator dVolumeUpdator;
    dVolumeUpdator.connect<&BoundVolumesGrid1Impl::update>(
        std::move(volumes1D));
    p->assignDetectorVolumeUpdator(dir, std::move(dVolumeUpdator), volumes);
  }
}

std::vector<std::shared_ptr<Acts::Experimental::DetectorVolume>>
Acts::Experimental::detail::PortalHelper::attachedDetectorVolumes(
    Portal& portal) noexcept(false) {
  auto& attachedVolumes = portal.attachedDetectorVolumes();
  if (not attachedVolumes[0u].empty() and not attachedVolumes[1u].empty()) {
    throw std::invalid_argument(
        "PortalHelper: trying to get attachedVolumes from already populated "
        "portal.");
  }
  unsigned int iu = attachedVolumes[0u].empty() ? 1u : 0u;
  return attachedVolumes[iu];
}