NuclearInteractionOptions.hpp

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

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

#pragma once

#include "Acts/Definitions/Units.hpp"
#include "Acts/Utilities/Logger.hpp"
#include "ActsExamples/Utilities/OptionsFwd.hpp"
#include "ActsFatras/Physics/NuclearInteraction/NuclearInteraction.hpp"
#include "ActsFatras/Physics/NuclearInteraction/NuclearInteractionParameters.hpp"

#include <utility>

#include "boost/program_options.hpp"

namespace ActsExamples {
namespace Options {

void addNuclearInteractionOptions(Description& desc);

ActsFatras::detail::MultiParticleNuclearInteractionParametrisation
readParametrisations(const std::string& fileName);

std::string readNuclearInteractionConfig(
    const boost::program_options::variables_map& variables);

template <typename simulator_t>
void setNuclearInteractionParametrisations(
    const std::string& nuclearInteractionParametrisation,
    simulator_t& simulator) {
  if (nuclearInteractionParametrisation.empty()) {
    return;
  }

  auto& chargedNuclearInteraction =
      simulator.charged.interactions
          .template get<ActsFatras::NuclearInteraction>();
  auto& neutralNuclearInteraction =
      simulator.neutral.interactions
          .template get<ActsFatras::NuclearInteraction>();

  const auto mpp = readParametrisations(nuclearInteractionParametrisation);

  chargedNuclearInteraction.multiParticleParameterisation = mpp;
  neutralNuclearInteraction.multiParticleParameterisation = mpp;
}

}  // namespace Options
}  // namespace ActsExamples