SpinFunctions¶
-
class
ompy.SpinFunctions(Ex, J, model, pars)[source]¶ Bases:
objectCalculates spin distributions, spin cuts (…)
- Parameters
Methods Summary
Get spin distribution
gDisc_and_EB05(mass, NLDa, Eshift, Sn, …)Linear interpolation of the spin-cut between a spin cut “from the discrete levels” and EB05 RReference: Guttormsen et al., 2017, PRC 96, 024313
gEB05(mass, NLDa, Eshift[, Ex])Von Egidy & B PRC72,044311(2005), Eq.
gEB09_CT(mass[, Ex])The constant temperature (CT) formula - Von Egidy & B PRC80,054310, below Eq.
gEB09_emp(mass, Pa_prime[, Ex])Von Egidy & B PRC80,054310, Eq.(16) FG+CT
Get the square of the spin cut for a specified model
Methods Documentation
-
distibution()[source]¶ Get spin distribution
Note: Assuming equal parity
- Returns
- Spin distribution. Shape
depends on input Ex and J and is squeezed if only one of them is an array. If both are arrays: spinDist[Ex,J]
- Return type
spinDist (Tuple[float, np.ndarray])
-
gDisc_and_EB05(mass, NLDa, Eshift, Sn, sigma2_disc, Ex=None)[source]¶ Linear interpolation of the spin-cut between a spin cut “from the discrete levels” and EB05 RReference: Guttormsen et al., 2017, PRC 96, 024313
Note
We set sigma2(E<E_discrete) = sigma2(E_discrete). This is not specified in the article, and may have been done differently before.
- Parameters
mass (int) – The mass number of the residual nucleus
NLDa (float) – Level density parameter
Eshift (float) – Energy shift
Sn (float) – Neutron separation energy
sigma2_disc (Tuple[float, float]) – [float, float] [Energy, sigma2] from the discretes
Ex (float or Sequence, optional) – Excitation energy. Defaults to self.Ex
- Returns
Squared spincut
- Return type
Union[float, Sequence]
-
gEB05(mass, NLDa, Eshift, Ex=None)[source]¶ Von Egidy & B PRC72,044311(2005), Eq. (4) The rigid moment of inertia formula (RMI) FG+CT
-
gEB09_CT(mass, Ex=None)[source]¶ The constant temperature (CT) formula - Von Egidy & B PRC80,054310, below Eq. (8) - original ref: Von Egidy et al., NPA 481 (1988) 189, Eq. (3)
-
distibution()[source] Get spin distribution
Note: Assuming equal parity
- Returns
- Spin distribution. Shape
depends on input Ex and J and is squeezed if only one of them is an array. If both are arrays: spinDist[Ex,J]
- Return type
spinDist (Tuple[float, np.ndarray])
-
gDisc_and_EB05(mass, NLDa, Eshift, Sn, sigma2_disc, Ex=None)[source] Linear interpolation of the spin-cut between a spin cut “from the discrete levels” and EB05 RReference: Guttormsen et al., 2017, PRC 96, 024313
Note
We set sigma2(E<E_discrete) = sigma2(E_discrete). This is not specified in the article, and may have been done differently before.
- Parameters
mass (int) – The mass number of the residual nucleus
NLDa (float) – Level density parameter
Eshift (float) – Energy shift
Sn (float) – Neutron separation energy
sigma2_disc (Tuple[float, float]) – [float, float] [Energy, sigma2] from the discretes
Ex (float or Sequence, optional) – Excitation energy. Defaults to self.Ex
- Returns
Squared spincut
- Return type
Union[float, Sequence]
-
gEB05(mass, NLDa, Eshift, Ex=None)[source] Von Egidy & B PRC72,044311(2005), Eq. (4) The rigid moment of inertia formula (RMI) FG+CT
-
gEB09_CT(mass, Ex=None)[source] The constant temperature (CT) formula - Von Egidy & B PRC80,054310, below Eq. (8) - original ref: Von Egidy et al., NPA 481 (1988) 189, Eq. (3)
-
gEB09_emp(mass, Pa_prime, Ex=None)[source] Von Egidy & B PRC80,054310, Eq.(16) FG+CT
-
get_sigma2()[source] Get the square of the spin cut for a specified model