def METHOD_NAME(self):
    s = self.s
    s.data[1, 2, 2, 3] = 1000
    s = self.s.as_lazy()
    mask = np.zeros((3, 3), dtype=bool)
    mask[1, 2] = True
    s_mask = Signal2D(mask)
    s_beam_shift = s.get_direct_beam_position(method="blur", sigma=1)
    s_beam_shift.compute()
    s_beam_shift.make_linear_plane(mask=s_mask)
    s.center_direct_beam(shifts=s_beam_shift)
    s.compute()
    np.testing.assert_almost_equal(s.data[:, :, 8, 8], 10, decimal=5)
    s.data[:, :, 8, 8] = 0
    s.data[1, 2, 3, 2] = 0
    np.testing.assert_almost_equal(s.data, 0.0, decimal=4)