HadronicElasticModel.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 <corsika/framework/core/ParticleProperties.hpp>
#include <corsika/framework/random/RNGManager.hpp>

#include <corsika/framework/core/PhysicalConstants.hpp>
#include <corsika/framework/core/PhysicalUnits.hpp>

namespace corsika {

  class HadronicElasticInteraction
      : public InteractionProcess<HadronicElasticInteraction> {
  private:
    using SquaredHEPEnergyType = decltype(HEPEnergyType() * HEPEnergyType());

    using eV2 = decltype(square(electronvolt));
    using inveV2 = decltype(1 / square(electronvolt));

    inveV2 B(eV2 s) const;
    CrossSectionType getCrossSection(SquaredHEPEnergyType s) const;

  public:
    HadronicElasticInteraction( // x & y values taken from DL for pp collisions
        CrossSectionType x = 0.0217 * barn, CrossSectionType y = 0.05608 * barn);

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

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

  private:
    CrossSectionType parX_, parY_;

    static double constexpr gfEpsilon = 0.0808;
    static double constexpr gfEta = 0.4525;
    // Froissart-Martin is not violated up for sqrt s < 10^32 eV with these values [DL].

    default_prng_type& RNG_ =
        RNGManager<>::getInstance().getRandomStream("HadronicElasticModel");
  };

} // namespace corsika

#include <corsika/detail/modules/HadronicElasticModel.inl>