Calculations/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m

function ret = accelerationDueToSpontaneousEmissionProcess(this, SaturationParameter, TotalSaturationParameter, Detuning, Linewidth, WaveNumber)
    Vector = [2*rand(1)-1,2*rand(1)-1,2*rand(1)-1];
    Vector = Vector./norm(Vector);

    ScatteringRate = (0.5 * SaturationParameter * Linewidth) / (1 + TotalSaturationParameter + (4 * (Detuning / Linewidth)^2));
    NumberOfScatteringEvents = floor(ScatteringRate * this.TimeStep);

    if NumberOfScatteringEvents > 0
        ret = Vector.*((Helper.PhysicsConstants.PlanckConstantReduced * WaveNumber) / ...
            (Helper.PhysicsConstants.Dy164Mass * this.TimeStep)).* sqrt(NumberOfScatteringEvents);
    else
        ret = zeros(1,3);
    end
end