diff --git a/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m b/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m index f6a3d6b..02f0d4c 100644 --- a/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m +++ b/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m @@ -1,8 +1,8 @@ -function ret = accelerationDueToSpontaneousEmissionProcess(this, SaturationParameter, TotalSaturationParameter, Detuning, Linewidth, WaveNumber) +function ret = accelerationDueToSpontaneousEmissionProcess(this, SaturationParameter, TotalSaturationParameter, 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)); + ScatteringRate = (0.5 * SaturationParameter * Linewidth) / (1 + TotalSaturationParameter); NumberOfScatteringEvents = floor(ScatteringRate * this.TimeStep); if NumberOfScatteringEvents > 0