def METHOD_NAME(simulation_factory, snapshot_factory, special_pair_cls,
                  params, r_cut, force, energy):
    snapshot = snapshot_factory()
    sim = simulation_factory(snapshot)
    potential = special_pair_cls()
    potential.params['A-A'] = params
    potential.r_cut['A-A'] = r_cut
    pickling_check(potential)
    sim.operations.integrator = hoomd.md.Integrator(dt=0.005,
                                                    forces=[potential])
    sim.run(0)
    pickling_check(potential)