def METHOD_NAME(self):
    """Test calculation of hmax and interpolation of Xmax tables"""
    image = np.array([1, 1, 1])
    pixel_x = np.array([1, 1, 1]) * u.deg
    pixel_y = np.array([1, 1, 1]) * u.deg
    array_pointing = SkyCoord(alt=0 * u.deg, az=0 * u.deg, frame=self.horizon_frame)
    self.impact_reco.set_event_properties(
        {1: image},
        {1: image},
        {1: pixel_x},
        {1: pixel_y},
        {1: "DUMMY"},
        {1: 0 * u.m},
        {1: 0 * u.m},
        array_direction=array_pointing,
        hillas={1: self.h1},
    )
    shower_max = self.impact_reco.get_shower_max(0, 0, 0, 100, 0)
    assert_allclose(shower_max, 484.2442217190515, rtol=0.01)