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);