InteractionProcess.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/process/BaseProcess.hpp>
#include <corsika/framework/core/PhysicalUnits.hpp>

#include <corsika/detail/framework/process/InteractionProcess.hpp> // for extra traits, method/interface checking

namespace corsika {

  template <typename TDerived>
  class InteractionProcess : public BaseProcess<TDerived> {

  public:
    using BaseProcess<TDerived>::ref;

    template <typename TParticle>
    InverseGrammageType getInverseInteractionLength(TParticle const& particle) {

      // interface checking on TProcess1
      static_assert(
          has_method_getInteractionLength_v<TDerived, GrammageType, TParticle const&>,
          "TDerived has no method with correct signature \"GrammageType "
          "getInteractionLength(TParticle const&)\" required for "
          "InteractionProcess<TDerived>. ");

      return 1. / ref().getInteractionLength(particle);
    }
  };

  template <typename TProcess>
  struct is_interaction_process<
      TProcess, std::enable_if_t<
                    std::is_base_of_v<InteractionProcess<typename std::decay_t<TProcess>>,
                                      typename std::decay_t<TProcess>>>>
      : std::true_type {};

} // namespace corsika