def test_pwl_pivot_energy():
pwl = PowerLawSpectralModel(amplitude="5.35510540e-11 cm-2 s-1 TeV-1")
pwl.covariance = [
[0.0318377 ** 2, 6.56889442e-14, 0],
[6.56889442e-14, 0, 0],
[0, 0, 0],
]
assert_quantity_allclose(pwl.pivot_energy, 3.3540034240210987 * u.TeV)
