Interaction.hpp

/*
 * (c) Copyright 2020 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 <PROPOSAL/PROPOSAL.h>

#include <corsika/framework/core/ParticleProperties.hpp>
#include <corsika/framework/process/InteractionProcess.hpp>
#include <corsika/framework/process/ProcessReturn.hpp>
#include <corsika/framework/random/RNGManager.hpp>
#include <corsika/framework/random/UniformRealDistribution.hpp>

#include <corsika/modules/proposal/ProposalProcessBase.hpp>

namespace corsika::proposal {

  class Interaction : public InteractionProcess<Interaction>, ProposalProcessBase {

    enum { eSECONDARIES, eINTERACTION };
    using calculator_t = std::tuple<std::unique_ptr<PROPOSAL::SecondariesCalculator>,
                                    std::unique_ptr<PROPOSAL::Interaction>>;

    std::unordered_map<calc_key_t, calculator_t, hash>
        calc; 

    void buildCalculator(Code, NuclearComposition const&) final;

  public:
    template <typename TEnvironment>
    Interaction(TEnvironment const& env);

    template <typename TSecondaryView>
    ProcessReturn doInteraction(TSecondaryView&);

    template <typename TParticle>
    GrammageType getInteractionLength(TParticle const& p);
  };
} // namespace corsika::proposal

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