diff --git a/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m b/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m index 02f0d4c..9cf4a31 100644 --- a/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m +++ b/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m @@ -2,7 +2,7 @@ function ret = accelerationDueToSpontaneousEmissionProcess(this, SaturationParam Vector = [2*rand(1)-1,2*rand(1)-1,2*rand(1)-1]; Vector = Vector./norm(Vector); - ScatteringRate = (0.5 * SaturationParameter * Linewidth) / (1 + TotalSaturationParameter); + ScatteringRate = (0.5 * Linewidth) * (SaturationParameter / (1 + TotalSaturationParameter)); NumberOfScatteringEvents = floor(ScatteringRate * this.TimeStep); if NumberOfScatteringEvents > 0