em/process/MuIonizationProcess.hh

0001 //------------------------------- -*- C++ -*- -------------------------------//
0002 // Copyright Celeritas contributors: see top-level COPYRIGHT file for details
0003 // SPDX-License-Identifier: (Apache-2.0 OR MIT)
0004 //---------------------------------------------------------------------------//
0005 //! \file celeritas/em/process/MuIonizationProcess.hh
0006 //---------------------------------------------------------------------------//
0007 #pragma once
0008
0009 #include <memory>
0010
0011 #include "celeritas/phys/Applicability.hh"
0012 #include "celeritas/phys/ImportedProcessAdapter.hh"
0013 #include "celeritas/phys/ParticleParams.hh"
0014 #include "celeritas/phys/Process.hh"
0015
0016 namespace celeritas
0017 {
0018 //---------------------------------------------------------------------------//
0019 /*!
0020  * Ionization process for muons.
0021  */
0022 class MuIonizationProcess : public Process
0023 {
0024   public:
0025     //!@{
0026     //! \name Type aliases
0027     using Energy = units::MevEnergy;
0028     using SPConstParticles = std::shared_ptr<ParticleParams const>;
0029     using SPConstImported = std::shared_ptr<ImportedProcesses const>;
0030     //!@}
0031
0032     // Options for electron and positron ionization
0033     struct Options
0034     {
0035         //! Maximum energy for the Bragg and ICRU73QO models
0036         Energy bragg_icru73qo_upper_limit{0.2};  //!< 200 keV
0037         //! Maximum energy for the Bethe-Bloch model
0038         Energy bethe_bloch_upper_limit{1e3};  //!< 1 GeV
0039     };
0040
0041   public:
0042     // Construct with imported data
0043     MuIonizationProcess(SPConstParticles particles,
0044                         SPConstImported process_data,
0045                         Options options);
0046
0047     // Construct the models associated with this process
0048     VecModel build_models(ActionIdIter start_id) const final;
0049
0050     // Get the interaction cross sections for the given energy range
0051     StepLimitBuilders step_limits(Applicability applicability) const final;
0052
0053     //! Whether the integral method can be used to sample interaction length
0054     bool supports_integral_xs() const final { return true; }
0055
0056     //! Whether the process applies when the particle is stopped
0057     bool applies_at_rest() const final { return imported_.applies_at_rest(); }
0058
0059     // Name of the process
0060     std::string_view label() const final;
0061
0062   private:
0063     SPConstParticles particles_;
0064     ImportedProcessAdapter imported_;
0065     Options options_;
0066 };
0067
0068 //---------------------------------------------------------------------------//
0069 }  // namespace celeritas