From 960bd02f8b0b69f50f7433bfb9b6f1d6e042e883 Mon Sep 17 00:00:00 2001
From: Karthik Chandrashekara <karthik@physi.uni-heidelberg.de>
Date: Sun, 11 Jul 2021 06:15:02 +0200
Subject: [PATCH] Cosmetic changes like the change of variable names to reflect
 their correct definition and use.

---
 .../accelerationDueToSpontaneousEmissionProcess.m           | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m b/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m
index 2f02c09..f6a3d6b 100644
--- a/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m	
+++ b/MOT Capture Process Simulation/@MOTSimulator/accelerationDueToSpontaneousEmissionProcess.m	
@@ -1,12 +1,12 @@
-function ret = accelerationDueToSpontaneousEmissionProcess(this, SaturationIntensity, TotalSaturationIntensity, Detuning, Linewidth, WaveVector)
+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 * SaturationIntensity * Linewidth) / (1 + TotalSaturationIntensity + (4 * (Detuning/Linewidth)^2));
+    ScatteringRate = (0.5 * SaturationParameter * Linewidth) / (1 + TotalSaturationParameter + (4 * (Detuning / Linewidth)^2));
     NumberOfScatteringEvents = floor(ScatteringRate * this.TimeStep);
 
     if NumberOfScatteringEvents > 0
-        ret = Vector.*((Helper.PhysicsConstants.PlanckConstantReduced * WaveVector) / ...
+        ret = Vector.*((Helper.PhysicsConstants.PlanckConstantReduced * WaveNumber) / ...
             (Helper.PhysicsConstants.Dy164Mass * this.TimeStep)).* sqrt(NumberOfScatteringEvents);
     else
         ret = zeros(1,3);