change method of multiple scattering
This fixes the underlying issue of #603 (closed).
We now sample a space moliere angle from PROPOSAL (where the interface is much clearer now). We then rotate the initial particle direction by this moliere angle. Afterwards, we introduce a random azimuith between 0 and 2PI with a second rotation.