forked from celeritas-project/celeritas
-
Notifications
You must be signed in to change notification settings - Fork 0
/
PhotoelectricProcess.cc
71 lines (64 loc) · 2.37 KB
/
PhotoelectricProcess.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
//----------------------------------*-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/PhotoelectricProcess.cc
//---------------------------------------------------------------------------//
#include "PhotoelectricProcess.hh"
#include <utility>
#include "corecel/Assert.hh"
#include "corecel/cont/Range.hh"
#include "celeritas/em/model/LivermorePEModel.hh"
#include "celeritas/io/ImportProcess.hh"
#include "celeritas/phys/PDGNumber.hh"
namespace celeritas
{
//---------------------------------------------------------------------------//
/*!
* Construct from host data.
*/
PhotoelectricProcess::PhotoelectricProcess(SPConstParticles particles,
SPConstMaterials materials,
SPConstImported process_data,
ReadData load_data)
: particles_(std::move(particles))
, materials_(std::move(materials))
, imported_(process_data,
particles_,
ImportProcessClass::photoelectric,
{pdg::gamma()})
, load_pe_(std::move(load_data))
{
CELER_EXPECT(particles_);
CELER_EXPECT(materials_);
CELER_EXPECT(load_pe_);
}
//---------------------------------------------------------------------------//
/*!
* Construct the models associated with this process.
*/
auto PhotoelectricProcess::build_models(ActionIdIter start_id) const -> VecModel
{
return {std::make_shared<LivermorePEModel>(
*start_id++, *particles_, *materials_, load_pe_)};
}
//---------------------------------------------------------------------------//
/*!
* Get the interaction cross sections for the given energy range.
*/
auto PhotoelectricProcess::step_limits(Applicability applic) const
-> StepLimitBuilders
{
return imported_.step_limits(std::move(applic));
}
//---------------------------------------------------------------------------//
/*!
* Name of the process.
*/
std::string_view PhotoelectricProcess::label() const
{
return "Photoelectric effect";
}
//---------------------------------------------------------------------------//
} // namespace celeritas