def METHOD_NAME():
    ds = fake_sph_grid_ds()
    @particle_filter("central_gas", requires=["particle_position"], filtered_type="io")
    def _filter(pfilter, data):
        coords = np.abs(data[pfilter.filtered_type, "particle_position"])
        return (coords[:, 0] < 1.6) & (coords[:, 1] < 1.6) & (coords[:, 2] < 1.6)
    ds.add_particle_filter("central_gas")
    plot = ProjectionPlot(ds, "z", ("central_gas", "density"))
    tmpdir = tempfile.mkdtemp()
    curdir = os.getcwd()
    os.chdir(tmpdir)
    plot.save()
    os.chdir(curdir)
    shutil.rmtree(tmpdir)