bremsstrahlung_thick_target#

sunkit_spex.legacy.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