def METHOD_NAME(lt_ctx):
    data = _mk_random(size=(16, 16, 16, 16), dtype='>u2')
    dataset = MemoryDataSet(data=data, tileshape=(8, 16, 16), num_partitions=32)
    expected = data.sum(axis=(0, 1))
    analysis = lt_ctx.create_sum_analysis(dataset=dataset)
    results = lt_ctx.run(analysis)
    assert results['intensity'].raw_data.shape == (16, 16)
    assert np.allclose(results['intensity'].raw_data, expected)
    assert np.allclose(results['intensity_lin'].raw_data, expected)