bremsstrahlung_thick_target#

sunkit_spex.emission.bremsstrahlung_thick_target(photon_energies, p, eebrk, q, eelow, eehigh, integrator=None)[source]#

Computes the thick-target bremsstrahlung x-ray/gamma-ray spectrum from an isotropic electron distribution function provided in broken_powerlaw_f. The units of the computed flux is photons per second per keV per square centimeter.

The electron flux distribution function is a double power law in electron energy with a low-energy cutoff and a high-energy cutoff.

Parameters:

photon_energies (numpy.array) – Array of photon energies to evaluate flux at
p (float) – Slope below the break energy
eebrk (float) – Break energy
q (float) – Slope above the break energy
eelow (float) – Low energy electron cut off
eehigh (float) – High energy electron cut off
integrator (callable) – A Python function or method to integrate must support vector limits and match signture fun(x, a, b, n, *args, **kwargs)

Returns:

flux The computed bremsstrahlung photon flux at the given photon energies. Array of photon fluxes (in photons s^-1 keV^-1 cm^-2), when multiplied by a0 * 1.0d+35, corresponding to the photon energies in the input array eph. The detector is assumed to be 1 AU from the source. a0 is the total integrated electron flux, in units of 10^35 electrons s^-1.

Return type:

numpy.array

Notes

If you want to plot the derivative of the flux, or the spectral index of the photon spectrum as a function of photon energy, you should set RERR to 1.d-6, because it is more sensitive to RERR than the flux.

Adapted from SSW Brm2_ThickTarget