From 2628e7c8bc51b8c8bc18a16087930c9212264be7 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Thu, 15 Jul 2021 16:52:58 +0200 Subject: [PATCH] Minor changes during testing. --- test_MOTCaptureProcessSimulation.m | 35 ++++++++++++++++++++---------- 1 file changed, 23 insertions(+), 12 deletions(-) diff --git a/test_MOTCaptureProcessSimulation.m b/test_MOTCaptureProcessSimulation.m index 95db29e..24054cf 100644 --- a/test_MOTCaptureProcessSimulation.m +++ b/test_MOTCaptureProcessSimulation.m @@ -24,10 +24,20 @@ MOT2D = Simulator.TwoDimensionalMOT(options{:}); Beams = MOT2D.Beams; %% - Run Simulation -poolobj = gcp('nocreate'); % Check if pool is open -if isempty(poolobj) - parpool; -end +% poolobj = gcp('nocreate'); % Check if pool is open +% if isempty(poolobj) +% parpool; +% end +MOT2D.NumberOfAtoms = 5000; +MOT2D.Sideband = false; +CoolingBeam = Beams{cellfun(@(x) strcmpi(x.Alias, 'Blue'), Beams)}; +CoolingBeam.Power = 0.4; +CoolingBeam.Waist = 13.3e-03; +CoolingBeam.Detuning = -1.67*Helper.PhysicsConstants.BlueLinewidth; +PushBeam = Beams{cellfun(@(x) strcmpi(x.Alias, 'Push'), Beams)}; +PushBeam.Power = 0.025; +PushBeam.Waist = 0.81e-03; +PushBeam.Detuning = 0; [LoadingRate, ~] = MOT2D.runSimulation(Oven); %% - Plot initial distribution % - sampling the position distribution @@ -118,13 +128,13 @@ QuantityOfInterestArray = LoadingRateArray; OptionsStruct = struct; OptionsStruct.RescalingFactorForParameter = 1000; OptionsStruct.XLabelString = 'Cooling Beam Power (mW)'; -OptionsStruct.RescalingFactorForYQuantity = 1e-10; +OptionsStruct.RescalingFactorForYQuantity = 1e-11; OptionsStruct.ErrorsForYQuantity = true; OptionsStruct.ErrorsArray = StandardErrorArray; OptionsStruct.CIForYQuantity = true; OptionsStruct.CIArray = ConfidenceIntervalArray; OptionsStruct.RemoveOutliers = true; -OptionsStruct.YLabelString = 'Loading rate (x 10^{10} atoms/s)'; +OptionsStruct.YLabelString = 'Loading rate (x 10^{11} atoms/s)'; OptionsStruct.TitleString = sprintf('Magnetic Gradient = %.0f (G/cm)', MOT2D.MagneticGradient * 100); options = Helper.convertstruct2cell(OptionsStruct); @@ -135,9 +145,9 @@ clear OptionsStruct %% - Scan parameters: Two-Parameter Scan -MOT2D.NumberOfAtoms = 5000; +MOT2D.NumberOfAtoms = 50; MOT2D.TotalPower = 0.6; -MOT2D.Sideband = true; +MOT2D.Sideband = false; SidebandBeam = Beams{cellfun(@(x) strcmpi(x.Alias, 'BlueSideband'), Beams)}; NumberOfPointsForFirstParam = 10; %iterations of the simulation @@ -145,8 +155,9 @@ NumberOfPointsForSecondParam = 10; % Scan Sideband Detuning and Power Ratio DetuningArray = linspace(-0.5,-10, NumberOfPointsForFirstParam) * Helper.PhysicsConstants.BlueLinewidth; -SidebandPowerArray = linspace(0.1,0.9, NumberOfPointsForSecondParam) * MOT2D.TotalPower; -BluePowerArray = MOT2D.TotalPower - SidebandPowerArray; +% SidebandPowerArray = linspace(0.1,0.9, NumberOfPointsForSecondParam) * MOT2D.TotalPower; +% BluePowerArray = MOT2D.TotalPower - SidebandPowerArray; +BluePowerArray = linspace(0.1,0.9, NumberOfPointsForSecondParam) * MOT2D.TotalPower; LoadingRateArray = zeros(NumberOfPointsForFirstParam, NumberOfPointsForSecondParam); StandardErrorArray = zeros(NumberOfPointsForFirstParam, NumberOfPointsForSecondParam); @@ -177,8 +188,8 @@ OptionsStruct.RescalingFactorForFirstParameter = (Helper.PhysicsConstants.B OptionsStruct.XLabelString = 'Sideband Detuning (\Delta/\Gamma)'; OptionsStruct.RescalingFactorForSecondParameter = 1000; OptionsStruct.YLabelString = 'Sideband Power (mW)'; -OptionsStruct.RescalingFactorForQuantityOfInterest = 1e-10; -OptionsStruct.ZLabelString = 'Loading rate (x 10^{10} atoms/s)'; +OptionsStruct.RescalingFactorForQuantityOfInterest = 1e-11; +OptionsStruct.ZLabelString = 'Loading rate (x 10^{11} atoms/s)'; OptionsStruct.TitleString = sprintf('Magnetic Gradient = %.0f (G/cm)', MOT2D.MagneticGradient * 100); options = Helper.convertstruct2cell(OptionsStruct);