hazma.decay.charged_pion¶

hazma.decay.charged_pion(eng_gam, eng_pi)¶

Compute dNdE from charged pion decay.

Compute dNdE from decay \(\pi^{\pm} \to \mu^{\pm} + \nu_{\mu} \to e^{\pm} + \nu_{e} + \nu_{\mu} + \gamma\) in the laborartory frame given a gamma ray engergy of eng_gam and muon energy of eng_pi.

Parameters:	eng_gam (double or numpy.ndarray) – Gamma ray energy(ies) in laboratory frame. eng_pi (double) – Charged pion energy in laboratory frame.
Returns:	spec (double np.ndarray) – List of gamma ray spectrum values, dNdE, evaluated at eng_gams given charged pion energy eng_pi.

Examples

Calculate spectrum for single gamma ray energy

>>> from hazma import decay
>>> eng_gam, eng_pi = 200., 1000.
>>> spec = decay.charged_pion(eng_gam, eng_pi)

Calculate spectrum for array of gamma ray energies

>>> from hazma import decay
>>> import numpy as np
>>> eng_gams = np.logspace(0.0, 3.0, num=200, dtype=float)
>>> eng_pi = 1000.
>>> spec = decay.charged_pion(eng_gams, eng_pi)