In SESAME the direct simulation model is implemented as follows (R represents a random number between 0 and 1):
The polar angle of deflection is calculated equivalently to the expression A.16 in the single scattering model.
The energy loss is determined by a random number and the
following equation: