InteractionModel.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 ofp
 * the license.
 */

#pragma once

#include <corsika/modules/qgsjetII/ParticleConversion.hpp>
#include <qgsjet-II-04.hpp>

#include <corsika/framework/core/ParticleProperties.hpp>
#include <corsika/framework/core/PhysicalUnits.hpp>
#include <corsika/framework/random/RNGManager.hpp>
#include <corsika/framework/utility/COMBoost.hpp>
#include <corsika/framework/utility/CorsikaData.hpp>

#include <boost/filesystem/path.hpp>

#include <random>

namespace corsika::qgsjetII {

  class InteractionModel {

  public:
    InteractionModel(boost::filesystem::path dataPath = corsika_data("QGSJetII"));
    ~InteractionModel();

    bool isValid(Code const beamId, Code const targetId, HEPEnergyType const sqrtS) const;

    CrossSectionType getCrossSection(Code const projectile, Code const target,
                                     FourMomentum const& projectileP4,
                                     FourMomentum const& targetP4) const;

    template <typename TSecondaries>
    void doInteraction(TSecondaries&, Code const projectile, Code const target,
                       FourMomentum const& projectileP4, FourMomentum const& targetP4);

  private:
    int count_ = 0;
    QgsjetIIHadronType alternate_ =
        QgsjetIIHadronType::PiPlusType; // for pi0, rho0 projectiles

    corsika::default_prng_type& rng_ =
        corsika::RNGManager<>::getInstance().getRandomStream("qgsjet");
    std::bernoulli_distribution bernoulli_;
    static size_t constexpr maxMassNumber_ = 208;
    static HEPEnergyType constexpr sqrtSmin_ = 10_GeV;
  };

} // namespace corsika::qgsjetII

#include <corsika/detail/modules/qgsjetII/InteractionModel.inl>