forked from celeritas-project/celeritas
/
EPlusAnnihilationProcess.cc
72 lines (63 loc) · 2.29 KB
/
EPlusAnnihilationProcess.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
69
70
71
72
//----------------------------------*-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/EPlusAnnihilationProcess.cc
//---------------------------------------------------------------------------//
#include "EPlusAnnihilationProcess.hh"
#include <memory>
#include <type_traits>
#include <utility>
#include "corecel/cont/Range.hh"
#include "celeritas/em/model/EPlusGGModel.hh"
#include "celeritas/grid/ValueGridBuilder.hh"
#include "celeritas/grid/ValueGridType.hh"
#include "celeritas/phys/Model.hh"
#include "celeritas/phys/PDGNumber.hh"
namespace celeritas
{
//---------------------------------------------------------------------------//
/*!
* Construct from host data.
*/
EPlusAnnihilationProcess::EPlusAnnihilationProcess(SPConstParticles particles,
Options options)
: particles_(std::move(particles))
, positron_id_(particles_->find(pdg::positron()))
, options_(options)
{
CELER_EXPECT(particles_);
CELER_ENSURE(positron_id_);
}
//---------------------------------------------------------------------------//
/*!
* Construct the models associated with this process.
*/
auto EPlusAnnihilationProcess::build_models(ActionIdIter start_id) const
-> VecModel
{
return {std::make_shared<EPlusGGModel>(*start_id++, *particles_)};
}
//---------------------------------------------------------------------------//
/*!
* Get the interaction cross sections for the given energy range.
*/
auto EPlusAnnihilationProcess::step_limits(Applicability range) const
-> StepLimitBuilders
{
CELER_EXPECT(range.particle == positron_id_);
StepLimitBuilders builders;
builders[ValueGridType::macro_xs] = std::make_unique<ValueGridOTFBuilder>();
return builders;
}
//---------------------------------------------------------------------------//
/*!
* Name of the process.
*/
std::string EPlusAnnihilationProcess::label() const
{
return "Positron annihiliation";
}
//---------------------------------------------------------------------------//
} // namespace celeritas