SurfaceSortingAlgorithm.cpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2020 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 "ActsExamples/TrackFitting/SurfaceSortingAlgorithm.hpp"

#include "ActsExamples/EventData/ProtoTrack.hpp"
#include "ActsExamples/EventData/SimHit.hpp"
#include "ActsFatras/EventData/Hit.hpp"

#include <cstddef>
#include <memory>
#include <stdexcept>
#include <utility>
#include <vector>

namespace ActsExamples {
struct AlgorithmContext;
}  // namespace ActsExamples

ActsExamples::SurfaceSortingAlgorithm::SurfaceSortingAlgorithm(
    Config cfg, Acts::Logging::Level level)
    : ActsExamples::IAlgorithm("SurfaceSortingAlgorithm", level),
      m_cfg(std::move(cfg)) {
  if (m_cfg.inputProtoTracks.empty()) {
    throw std::invalid_argument("Missing input proto track collection");
  }
  if (m_cfg.inputSimHits.empty()) {
    throw std::invalid_argument("Missing input simulated hits collection");
  }
  if (m_cfg.inputMeasurementSimHitsMap.empty()) {
    throw std::invalid_argument("Missing input measurement sim hits map");
  }
  if (m_cfg.outputProtoTracks.empty()) {
    throw std::invalid_argument("Missing output proto track collection");
  }

  m_inputProtoTracks.initialize(m_cfg.inputProtoTracks);
  m_inputSimHits.initialize(m_cfg.inputSimHits);
  m_inputMeasurementSimHitsMap.initialize(m_cfg.inputMeasurementSimHitsMap);
  m_outputProtoTracks.initialize(m_cfg.outputProtoTracks);
}

ActsExamples::ProcessCode ActsExamples::SurfaceSortingAlgorithm::execute(
    const ActsExamples::AlgorithmContext& ctx) const {
  const auto& protoTracks = m_inputProtoTracks(ctx);
  const auto& simHits = m_inputSimHits(ctx);
  const auto& simHitsMap = m_inputMeasurementSimHitsMap(ctx);

  ProtoTrackContainer sortedTracks;
  sortedTracks.reserve(protoTracks.size());
  TrackHitList trackHitList;

  for (std::size_t itrack = 0; itrack < protoTracks.size(); ++itrack) {
    const auto& protoTrack = protoTracks[itrack];

    ProtoTrack sortedProtoTrack;
    sortedProtoTrack.reserve(protoTrack.size());
    trackHitList.clear();

    if (protoTrack.empty()) {
      continue;
    }

    for (const auto hit : protoTrack) {
      const auto simHitIndex = simHitsMap.find(hit)->second;
      auto simHit = simHits.nth(simHitIndex);
      auto simHitTime = simHit->time();
      trackHitList.insert(std::make_pair(simHitTime, hit));
    }

    for (auto const& [time, hit] : trackHitList) {
      sortedProtoTrack.emplace_back(hit);
    }

    sortedTracks.emplace_back(std::move(sortedProtoTrack));
  }

  m_outputProtoTracks(ctx, std::move(sortedTracks));

  return ActsExamples::ProcessCode::SUCCESS;
}