def METHOD_NAME():
    cat_index, finder_id, ntot, mvir = stat.HaloStatFile(ts2.filename).read("npart","Mvir")
    assert all(finder_id==[1,49,52,58,94,121,127,148,163])
    assert all(cat_index==[1,49,52,58,94,121,127,148,163])
    assert all(ntot==[3273314, 27631, 24654, 22366, 12915, 9831, 9498, 8200, 7256])
    npt.assert_almost_equal(mvir,   [  1.12282280e+12 ,  1.19939950e+10 ,  1.19538740e+10  , 3.07825010e+10,
                                       1.76325820e+10 ,  1.33353700e+10  , 1.28836660e+10 ,  1.11229900e+10,
                                       9.84248220e+09], decimal=5)