Added switch to alchemy.py, but will have to separate Lennard-Jones and electrostatics into separate CustomNonbondedForce terms because long-range correction can only be applied to Lennard-Jones.