IRFDOC_ELECTRON_SPECTRUM

ctao_cr_spectra.IRFDOC_ELECTRON_SPECTRUM = PowerLawWithExponentialGaussian(2.385e-09 1 / (TeV s sr cm2) * (E / 1.0 TeV)**-3.43 * (1 + 1.95 * (exp(Gauss(log10(E / 1.0 TeV), -0.101, 0.741)) - 1))

A power law with an additional Gaussian bump.

Beware that the Gaussian is not normalized!

\[\Phi(E, \Phi_0, \gamma, f, \mu, \sigma, E_\text{ref}) = \Phi_0 \left( \frac{E}{E_\text{ref}} \right)^{\gamma} \cdot \left( 1 + f \cdot \left( \exp\left( \operatorname{Gauss}(\log_{10}(E / E_\text{ref}), \mu, \sigma) \right) - 1 \right) \right)\]

Where \(\operatorname{Gauss}\) is the unnormalized Gaussian distribution:

\[\operatorname{Gauss}(x, \mu, \sigma) = \exp\left( -\frac{1}{2} \left(\frac{x - \mu}{\sigma}\right)^2 \right)\]

Attributes:

normalization: astropy.units.Quantity[flux]

\(\Phi_0\),

a: float

\(\alpha\)

b: float

\(\beta\)

e_ref: astropy.units.Quantity[energy]

\(E_\text{ref}\)