def METHOD_NAME(hdv2_stack):
    # Test energy bands
    ebands = hdv2_stack.energy_bands([0.1, 60] * u.MeV, 0.1 * u.MeV, dx=1 * u.um)
    # Expected energy bands, in MeV (only in active layers)
    expected = np.array([[3.5, 3.8], [4.6, 4.9], [5.6, 5.7], [6.4, 6.5], [7.1, 7.2]])
    assert np.allclose(ebands.to(u.MeV).value[:5, :], expected, atol=0.15)