def METHOD_NAME(simulation_factory, two_particle_snapshot_factory):
    """Creates simulation with state initialized.
    Note:
        The positions of the particles correspond correctly to the xy plane
        manifold used for tests. Do not change positions unless sure that the
        manifold tests with running simulation will not error.
    """
    def sim_constructor(active_force=None, rd_updater=None):
        sim = simulation_factory(two_particle_snapshot_factory())
        if isinstance(active_force, hoomd.md.force.Active):
            method = hoomd.md.methods.ConstantVolume(hoomd.filter.All())
        else:
            method = hoomd.md.methods.rattle.NVE(hoomd.filter.All(),
                                                 hoomd.md.manifold.Plane(0.1))
        sim.operations.integrator = hoomd.md.Integrator(dt=0.005,
                                                        methods=[method])
        if active_force is not None:
            sim.operations.integrator.forces.append(active_force)
        if rd_updater is not None:
            sim.operations.updaters.append(rd_updater)
        return sim
    return sim_constructor