Beersheba: no interpolation creates an image of different size

[gonzaponte]

See this test for example

@fixture(scope="session")
def new_grid():
    det_db   = DataSiPM('new', 0)
    det_grid = [ np.arange( det_db[var].min() + bs/2
                          , det_db[var].max() - bs/2 + np.finfo(np.float32).eps, bs)
                 for var, bs in zip("X Y".split(), [1., 1.])]
    return det_grid


@mark.parametrize("interp_method", InterpolationMethod)
def test_deconvolution_input_interpolation_method(data_hdst, new_grid, interp_method):
    """
    Check that it runs with all interpolation methods
    """
    hdst = load_dst(data_hdst, 'RECO', 'Events')
    hdst = hdst[(hdst.event == 3021916) & (hdst.npeak == 0) & (hdst.Q>40)]

    deconvolve = deconvolution_input([10., 10.], new_grid, interp_method)
    output     = deconvolve((hdst.X, hdst.Y), hdst.Q)

    assert output[0].shape == (70, 100)
    assert len(output[1])  == 2
    assert output[1][0].shape == (7000,)
    assert output[1][1].shape == (7000,)