forked from celeritas-project/celeritas
-
Notifications
You must be signed in to change notification settings - Fork 0
/
GammaConversionProcess.cc
68 lines (61 loc) · 2.26 KB
/
GammaConversionProcess.cc
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
//----------------------------------*-C++-*----------------------------------//
// Copyright 2020-2024 UT-Battelle, LLC, and other Celeritas developers.
// See the top-level COPYRIGHT file for details.
// SPDX-License-Identifier: (Apache-2.0 OR MIT)
//---------------------------------------------------------------------------//
//! \file celeritas/em/process/GammaConversionProcess.cc
//---------------------------------------------------------------------------//
#include "GammaConversionProcess.hh"
#include <utility>
#include "corecel/Assert.hh"
#include "corecel/cont/Range.hh"
#include "celeritas/em/model/BetheHeitlerModel.hh"
#include "celeritas/io/ImportProcess.hh"
#include "celeritas/phys/PDGNumber.hh"
namespace celeritas
{
//---------------------------------------------------------------------------//
/*!
* Construct from particles and imported Geant data.
*/
GammaConversionProcess::GammaConversionProcess(SPConstParticles particles,
SPConstImported process_data,
Options options)
: particles_(std::move(particles))
, imported_(process_data,
particles_,
ImportProcessClass::conversion,
{pdg::gamma()})
, options_(options)
{
CELER_EXPECT(particles_);
}
//---------------------------------------------------------------------------//
/*!
* Construct the models associated with this process.
*/
auto GammaConversionProcess::build_models(ActionIdIter start_id) const
-> VecModel
{
return {std::make_shared<BetheHeitlerModel>(
*start_id++, *particles_, imported_.processes(), options_.enable_lpm)};
}
//---------------------------------------------------------------------------//
/*!
* Get the interaction cross sections for the given energy range.
*/
auto GammaConversionProcess::step_limits(Applicability applic) const
-> StepLimitBuilders
{
return imported_.step_limits(std::move(applic));
}
//---------------------------------------------------------------------------//
/*!
* Name of the process.
*/
std::string_view GammaConversionProcess::label() const
{
return "Photon annihiliation";
}
//---------------------------------------------------------------------------//
} // namespace celeritas