Interaction.hpp

/*
 * (c) Copyright 2018 CORSIKA Project, corsika-project@lists.kit.edu
 *
 * This software is distributed under the terms of the GNU General Public
 * Licence version 3 (GPL Version 3). See file LICENSE for a full version of
 * the license.
 */

#pragma once

#include <corsika/framework/core/ParticleProperties.hpp>
#include <corsika/framework/core/PhysicalUnits.hpp>
#include <corsika/framework/random/RNGManager.hpp>
#include <corsika/framework/process/InteractionProcess.hpp>
#include <tuple>

namespace corsika::epos {

  class Interaction : public InteractionProcess<Interaction> {
    std::string data_path_;
    unsigned int count_ = 0;
    bool epos_listing_;

  public:
    Interaction(std::string const& dataPath = "", bool const epos_printout_on = false);
    ~Interaction();

    std::tuple<CrossSectionType, CrossSectionType> calcCrossSectionCoM(
        Code const, int const, int const, Code const, int const, int const,
        HEPEnergyType const) const;

    std::tuple<CrossSectionType, CrossSectionType> readCrossSectionTableLab(
        Code const, int const, int const, Code const, HEPEnergyType const) const;

    std::tuple<CrossSectionType, CrossSectionType> getCrossSectionLab(
        Code const, int const, int const, Code const, int const, int const,
        HEPEnergyType const) const;

    template <typename TParticle>
    GrammageType getInteractionLength(TParticle const&) const;

    template <typename TSecondaries>
    void doInteraction(TSecondaries&);

    bool isValidCoMEnergy(HEPEnergyType const ecm) const {
      return (minEnergyCoM_ <= ecm) && (ecm <= maxEnergyCoM_);
    }

    bool isValidTarget(Code const) const;

    void initialize() const;
    void initializeEventCoM(Code const, int const, int const, Code const, int const,
                            int const, HEPEnergyType const) const;
    void initializeEventLab(Code const, int const, int const, Code const, int const,
                            int const, HEPEnergyType const) const;
    void configureParticles(Code const, int const, int const, Code const, int const,
                            int const) const;
    void setParticlesStable() const;

  private:
    default_prng_type& RNG_ = RNGManager<>::getInstance().getRandomStream("epos");
    std::shared_ptr<spdlog::logger> logger_ = get_logger("corsika_epos_Interaction");
    HEPEnergyType const minEnergyCoM_ = 6 * 1e9 * electronvolt;
    HEPEnergyType const maxEnergyCoM_ = 2.e6 * 1e9 * electronvolt;
    int const maxTargetMassNumber_ = 20;
    int const minNuclearTargetA_ = 4;
  };

} // namespace corsika::epos

#include <corsika/detail/modules/epos/Interaction.inl>