em/process/MuIonizationProcess.hh

0001 //----------------------------------*-C++-*----------------------------------//
0002 // Copyright 2024 UT-Battelle, LLC, and other Celeritas developers.
0003 // See the top-level COPYRIGHT file for details.
0004 // SPDX-License-Identifier: (Apache-2.0 OR MIT)
0005 //---------------------------------------------------------------------------//
0006 //! \file celeritas/em/process/MuIonizationProcess.hh
0007 //---------------------------------------------------------------------------//
0008 #pragma once
0009
0010 #include <memory>
0011
0012 #include "celeritas/phys/Applicability.hh"
0013 #include "celeritas/phys/ImportedProcessAdapter.hh"
0014 #include "celeritas/phys/ParticleParams.hh"
0015 #include "celeritas/phys/Process.hh"
0016
0017 namespace celeritas
0018 {
0019 //---------------------------------------------------------------------------//
0020 /*!
0021  * Ionization process for muons.
0022  */
0023 class MuIonizationProcess : public Process
0024 {
0025   public:
0026     //!@{
0027     //! \name Type aliases
0028     using Energy = units::MevEnergy;
0029     using SPConstParticles = std::shared_ptr<ParticleParams const>;
0030     using SPConstImported = std::shared_ptr<ImportedProcesses const>;
0031     //!@}
0032
0033     // Options for electron and positron ionization
0034     struct Options
0035     {
0036         //! Maximum energy for the Bragg and ICRU73QO models
0037         Energy bragg_icru73qo_upper_limit{0.2};  //!< 200 keV
0038         //! Maximum energy for the Bethe-Bloch model
0039         Energy bethe_bloch_upper_limit{1e3};  //!< 1 GeV
0040         //! Use integral method for sampling discrete interaction length
0041         bool use_integral_xs{true};
0042     };
0043
0044   public:
0045     // Construct with imported data
0046     MuIonizationProcess(SPConstParticles particles,
0047                         SPConstImported process_data,
0048                         Options options);
0049
0050     // Construct the models associated with this process
0051     VecModel build_models(ActionIdIter start_id) const final;
0052
0053     // Get the interaction cross sections for the given energy range
0054     StepLimitBuilders step_limits(Applicability applicability) const final;
0055
0056     //! Whether to use the integral method to sample interaction length
0057     bool use_integral_xs() const final { return options_.use_integral_xs; }
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