Determines the distribution of times spent by atom in the MOT interaction volume to be used to compute the collision probability between the incoming hot atoms and the atoms slowed near the center of the MOT.

This commit is contained in:
Karthik 2021-07-11 14:38:37 +02:00
parent b6d56e9e7b
commit 951c2e3935

View File

@ -0,0 +1,20 @@
function T = computeTimeSpentInInteractionRegion(this, r)
% INPUT:
% r : N x 3 array. N is the number of time steps
% OUTPUT
% T : gives the distribution of time spent in the interaction region
% USAGE:
% T = this.computeTimeSpentInInteractionRegion(r)
T = 0;
NumberOfTimeSteps = int64(this.SimulationTime/this.TimeStep);
for n = 1:(NumberOfTimeSteps - 1)
dr = Helper.calculateDistanceFromPointToLine(r(n+1, :), [0 0 0], [1 0 0]);
if dr < this.CoolingBeamRadius
A = 1;
else
A = 0;
end
T = T + A * this.TimeStep;
end
end