SpinFunctions

class ompy.SpinFunctions(Ex, J, model, pars)[source]

Bases: object

Calculates spin distributions, spin cuts (…)

Parameters:
  • Ex (Union[float, Sequence]) – Excitation energy

  • J (Union[float, Sequence]) – Spin

  • model (str) – Modelname for the the spincut

  • pars (Dict[str, Any]) – Additional parameters necessary for the spin cut model

Methods Summary

distribution()

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 .

gDisc_and_sigmaSn(sigma2_disc, sigma2_Sn[, Ex])

Linear interpolation of the spin-cut between a spin cut "from the discrete levels" and a set value at the nutron separation energy. Essentially the same as gDisc_and_EB05 but with spin-cut at Sn explicitly set. Reference: Guttormsen et al., PRC 96, 024313 (2017) R. Capote et al., Nucl. Data Sheets 110, 3107-3214 (2009).

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, see sec.

gEB09_emp(mass, Pa_prime[, Ex])

Von Egidy & B PRC80,054310, Eq.(16) FG+CT

gconst(sigma[, Ex])

Constant spin-cutoff parameter

get_sigma2()

Get the square of the spin cut for a specified model

Methods Documentation

distribution()[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 .. rubric:: References

Guttormsen et al., PRC 96, 024313 (2017) R. Capote et al., Nucl. Data Sheets 110, 3107-3214 (2009)

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]

gDisc_and_sigmaSn(sigma2_disc, sigma2_Sn, Ex=None)[source]

Linear interpolation of the spin-cut between a spin cut “from the discrete levels” and a set value at the nutron separation energy. Essentially the same as gDisc_and_EB05 but with spin-cut at Sn explicitly set. Reference:

Guttormsen et al., PRC 96, 024313 (2017) R. Capote et al., Nucl. Data Sheets 110, 3107-3214 (2009)

Parameters:
  • sigma2_disc (Tuple[float, float]) – [float, float] [Energy, sigma2] for the discretes

  • sigma2_Sn (Tuple[float, float]) – [float, float] [Energy, sigma2] at Sn.

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

Parameters:
  • mass (int) – The mass number of the residual nucleus

  • NLDa (float) – Level density parameter

  • Eshift (float) – Energy shift

  • Ex (float or Sequence, optional) – Excitation energy. Defaults to self.Ex

Returns:

Squared spincut

Return type:

Union[float, Sequence]

gEB09_CT(mass, Ex=None)[source]

The constant temperature (CT) formula - Von Egidy & B PRC80,054310, see sec. IV, p7 refering to ref. below - original ref: Von Egidy et al., NPA 481 (1988) 189, Eq. (3)

Parameters:

mass (int) – Excitation energy

Returns:

Squared spincut

Return type:

Union[float, Sequence]

gEB09_emp(mass, Pa_prime, Ex=None)[source]

Von Egidy & B PRC80,054310, Eq.(16) FG+CT

Parameters:
  • mass (int) – Excitation energy

  • Pa_prime (float) – Deuteron pairing energy

  • Ex (float or Sequence, optional) – Excitation energy. Defaults to self.Ex

Returns:

Squared spincut

Return type:

Union[float, Sequence]

gconst(sigma, Ex=None)[source]

Constant spin-cutoff parameter

Parameters:

sigma (int) – Spin cut-off parameter

Returns:

Squared spincut

Return type:

Union[float, Sequence]

get_sigma2()[source]

Get the square of the spin cut for a specified model