ContinuousProcess.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 <PROPOSAL/PROPOSAL.h>

#include <corsika/media/Environment.hpp>

#include <corsika/framework/core/ParticleProperties.hpp>
#include <corsika/framework/process/ContinuousProcess.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>

#include <unordered_map>

namespace corsika::proposal {

  class ContinuousProcess
      : public corsika::ContinuousProcess<proposal::ContinuousProcess>,
        ProposalProcessBase {

    struct Calculator {
      std::unique_ptr<PROPOSAL::Displacement> disp;
      std::unique_ptr<PROPOSAL::Scattering> scatter;
    };

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

    HEPEnergyType energy_lost_ = 0 * electronvolt;

    void buildCalculator(Code, NuclearComposition const&) final;

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

    template <typename TParticle>
    void scatter(TParticle&, HEPEnergyType const&, GrammageType const&);

    template <typename TParticle, typename TTrack>
    ProcessReturn doContinuous(TParticle&, TTrack const& track, bool const limitFlag);

    template <typename TParticle, typename TTrack>
    LengthType getMaxStepLength(TParticle const&, TTrack const&);

    void showResults() const;
    void reset();
    HEPEnergyType getEnergyLost() const { return energy_lost_; }
  };
} // namespace corsika::proposal

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