From e80c68ef2a8a73e636748aa2b6e4d0d9dd7727dd Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Sun, 11 Jul 2021 06:00:08 +0200 Subject: [PATCH] Changed sampling of angle such that it is up to NozzleExitDivergence not MOTExitDivergence, made some other cosmetic changes. --- .../+Plotting/plotCaptureVelocityVsAngle.m | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/MOT Capture Process Simulation/+Plotting/plotCaptureVelocityVsAngle.m b/MOT Capture Process Simulation/+Plotting/plotCaptureVelocityVsAngle.m index bfff80e..7352270 100644 --- a/MOT Capture Process Simulation/+Plotting/plotCaptureVelocityVsAngle.m +++ b/MOT Capture Process Simulation/+Plotting/plotCaptureVelocityVsAngle.m @@ -1,6 +1,6 @@ function plotCaptureVelocityVsAngle(obj) - theta = linspace(0, obj.MOTExitDivergence, 40); + theta = linspace(0, obj.NozzleExitDivergence, 1000); CaptureVelocity = zeros(length(theta),3); for i=1:length(theta) @@ -20,16 +20,16 @@ function plotCaptureVelocityVsAngle(obj) screensize = get(0,'ScreenSize'); f_h.Position = [[screensize(3)/3.5 screensize(4)/3.5] 750 600]; - plot(theta, sqrt(CaptureVelocity(:,1).^2+CaptureVelocity(:,2).^2+CaptureVelocity(:,3).^2), 'Linewidth', 1.5) + plot(theta .* 1e+03, sqrt(CaptureVelocity(:,1).^2+CaptureVelocity(:,2).^2+CaptureVelocity(:,3).^2), 'Linewidth', 1.5) - hXLabel = xlabel('\theta (rad)'); + hXLabel = xlabel('\theta (mrad)'); hYLabel = ylabel('Velocity (m/s)'); - hTitle = sgtitle('Capture Velocity of an atomic beam from (0,0,0)'); + hTitle = sgtitle('Capture velocity for different angles of divergence'); set([hXLabel, hYLabel] , ... 'FontSize' , 14 ); set( hTitle , ... - 'FontSize' , 18 ); + 'FontSize' , 14 ); grid on Helper.bringFiguresWithTagInForeground();