def METHOD_NAME(self, psf):
    energy = 1 * u.TeV
    theta = 0 * u.deg
    containment = [0.68, 0.8, 0.9]
    desired = psf.containment_radius(
        energy_true=energy, offset=theta, fraction=containment
    )
    assert_allclose(desired, [0.14775, 0.18675, 0.25075] * u.deg, rtol=1e-3)