MARLEY (Model of Argon Reaction Low Energy Yields)
v1.2.0
A Monte Carlo event generator for tens-of-MeV neutrino interactions
|
Implements the Weisskopf single-particle estimates of the gamma-ray strength functions. More...
#include <WeisskopfSingleParticleModel.hh>
Public Member Functions | |
WeisskopfSingleParticleModel (int Z, int A, double D0=1.) | |
virtual double | strength_function (TransitionType type, int l, double e_gamma) override |
Returns the gamma-ray strength function (MeV –2 \(\ell\)–1) for the requested gamma energy and multipolarity. More... | |
virtual double | transmission_coefficient (TransitionType type, int l, double e_gamma) override |
Returns the gamma-ray transmission coefficient (dimensionless) for the requested gamma energy and multipolarity. More... | |
![]() | |
GammaStrengthFunctionModel (int Z, int A) | |
Additional Inherited Members | |
![]() | |
enum class | TransitionType { electric , magnetic , unphysical } |
Electromagnetic transitions in nuclei may be classified by their multipolarity (electric vs. magnetic multipole radiation) | |
![]() | |
static void | check_multipolarity (int l) |
Check that l > 0 and throw a marley::Error if it is not. More... | |
![]() | |
int | A_ |
Mass number. | |
int | Z_ |
Atomic number. | |
Implements the Weisskopf single-particle estimates of the gamma-ray strength functions.
Under this model, the electric ( \(f_{\text{E}\ell}\)) and magnetic ( \(f_{\text{M}\ell}\)) gamma-ray strength functions are independent of the gamma energy and are given by
\[ f_{\text{E}\ell} = \frac{2\alpha\Lambda}{\text{D}_0} \left(\frac{R}{\hbar c}\right)^{\!2\ell} \]
and
\[ f_{\text{M}\ell} = 10\left(\frac{\hbar c} {R\,m_\text{p}}\right)^{\!2} f_{\text{E}\ell}, \]
where \(\alpha\) is the fine-structure constant, \(\text{D}_0\) is the level spacing parameter (MeV), \(R = (1.2\text{ fm})A^{1/3}\) is the approximate nuclear radius, \(m_\text{p}\) is the proton mass (MeV), and \(\Lambda\) is a function of the multipolarity \(\ell\) given by
\[ \Lambda \equiv \left(\frac{3}{\ell + 3}\right)^{\!2} \left(\frac{\ell + 1}{\ell\left[(2\ell + 1)!!\right]^2} \right).\]
These estimates are typically only good to an order of magnitude, so using a more sophisticated model, e.g., the StandardLorentzianModel, is strongly recommended.
marley::WeisskopfSingleParticleModel::WeisskopfSingleParticleModel | ( | int | Z, |
int | A, | ||
double | D0 = 1. |
||
) |
Z | Atomic number of the desired nuclide |
A | Mass number of the desired nuclide |
D0 | Level spacing parameter (MeV) |
|
overridevirtual |
Returns the gamma-ray strength function (MeV –2 \(\ell\)–1) for the requested gamma energy and multipolarity.
type | Electric or magnetic transition |
l | Multipolarity of the transition |
e_gamma | Gamma-ray energy (MeV) |
Implements marley::GammaStrengthFunctionModel.
|
overridevirtual |
Returns the gamma-ray transmission coefficient (dimensionless) for the requested gamma energy and multipolarity.
The gamma-ray transmission coefficient and strength function are related via \(\text{T}_{\text{X}\ell}(\text{E}_\gamma) = 2\pi f_{\text{X}\ell}(\text{E}_\gamma)\text{E}_\gamma^{(2\ell + 1)},\) where X is the type of transition (electric or magnetic), \(\ell\) is the multipolarity, \(\text{T}_{\text{X}\ell}\) is the transmission coefficient, \(f_{\text{X}\ell}\) is the strength function, and \(\text{E}_\gamma\) is the gamma-ray energy.
type | Electric or magnetic transition |
l | Multipolarity of the transition |
e_gamma | Gamma-ray energy (MeV) |
Implements marley::GammaStrengthFunctionModel.