From 3a6732a4a64ab9244e1413ffb5ec46f6231e1630 Mon Sep 17 00:00:00 2001
From: Karthik Chandrashekara <karthik@physi.uni-heidelberg.de>
Date: Sun, 11 Jul 2021 21:16:07 +0200
Subject: [PATCH] Corrected scattering rate formula and removed detuning as an
 argument to be passed.

---
 .../accelerationDueToSpontaneousEmissionProcess.m             | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

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