Properly handle keep_sidechains flag in case of disulfid bonds

The ResolveCysteins function is used to optimize the disulfid bond
network. The function uses RotamerGroupPtr as input.
To allow fixed cysteins also being part of the optimization,
the rotamer groups get newly constructed containing only one
rotamer with the same SG position as in the input structure.