Source code for iact_estimator.tests.test_spectral

# Test cases automatically generated by Pynguin (https://www.pynguin.eu).
# Please check them before you use them.

import astropy.units as u
import numpy as np




def test_crab_nebula_spectrum():
    from ..spectral import crab_nebula_spectrum

    spectrum = crab_nebula_spectrum()

    assert spectrum.amplitude.unit == u.Unit("TeV^-1 s^-1 cm^-2")
    assert spectrum.reference.unit == u.TeV

    assert np.allclose(
        [
            spectrum.alpha.value,
            spectrum.beta.value,
            spectrum.amplitude.value,
            spectrum.reference.value,
        ],
        [2.51, 0.21 / np.log(10), 3.39e-11, 1.0],
    )