From b5169b7e979491075e674a995a581cd69a3a6458 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Tue, 10 Jun 2025 18:42:37 +0200 Subject: [PATCH] Minor modifications --- Data-Analyzer/analyzePhaseTransitionSimulation.m | 8 ++++---- Data-Analyzer/matchTFRadii.m | 2 +- Estimations/LinearizeScatteringLengthScan.m | 16 ++++++++-------- 3 files changed, 13 insertions(+), 13 deletions(-) diff --git a/Data-Analyzer/analyzePhaseTransitionSimulation.m b/Data-Analyzer/analyzePhaseTransitionSimulation.m index dc94bb6..e9a3ae9 100644 --- a/Data-Analyzer/analyzePhaseTransitionSimulation.m +++ b/Data-Analyzer/analyzePhaseTransitionSimulation.m @@ -1,6 +1,6 @@ %% Extract Images -baseDir = 'D:/Results - Numerics/Data_Full3D/PhaseTransition/DTS/'; +baseDir = 'D:/Results - Numerics/Data_Full3D/PhaseTransition/DTS/Point2'; JobNumber = 0; runFolder = sprintf('Run_%03d', JobNumber); movieFileName = 'DropletsToStripes.mp4'; % Output file name @@ -10,7 +10,7 @@ reverseOrder = false; % Set this to true to reverse the theta ordering TitleString = 'Change across transition: Droplets To Stripes'; %% -baseDir = 'D:/Results - Numerics/Data_Full3D/PhaseTransition/STD/'; +baseDir = 'D:/Results - Numerics/Data_Full3D/PhaseTransition/STD/Point2'; JobNumber = 0; runFolder = sprintf('Run_%03d', JobNumber); movieFileName = 'StripesToDroplets.mp4'; % Output file name @@ -183,6 +183,7 @@ for k = 1:nimgs y_max = max(y); x_min = min(x); x_max = max(x); + % Display the cropped OD image ax1 = nexttile; @@ -234,8 +235,7 @@ for k = 1:nimgs % Plot the angular distribution nexttile - spectral_contrast(k) = computeSpectralContrast(fft_imgs{k}, 10, 25, Threshold); - [theta_vals, S_theta] = computeNormalizedAngularSpectralDistribution(fft_imgs{k}, 10, 25, N_bins, Threshold, Sigma); + spectral_contrast(k) = computeSpectralContrast(fft_imgs{k}, 10, 35, Threshold); spectral_weight(k) = trapz(theta_vals, S_theta); plot(theta_vals/pi, S_theta,'Linewidth',2); axis square; diff --git a/Data-Analyzer/matchTFRadii.m b/Data-Analyzer/matchTFRadii.m index 8b5d4d2..99be214 100644 --- a/Data-Analyzer/matchTFRadii.m +++ b/Data-Analyzer/matchTFRadii.m @@ -1,7 +1,7 @@ % ---------------------------- % Experimental data % ---------------------------- -PixelSize = 5.86; % microns +PixelSize = 4.6; % microns AtomNumbers = [9.14950197, 8.6845907 , 8.82521245, 8.5899089 , 8.21675841, ... 8.96234044, 8.8636914 , 8.70332154, 8.82930706, 8.91869919, ... diff --git a/Estimations/LinearizeScatteringLengthScan.m b/Estimations/LinearizeScatteringLengthScan.m index 7ef6933..32e6a85 100644 --- a/Estimations/LinearizeScatteringLengthScan.m +++ b/Estimations/LinearizeScatteringLengthScan.m @@ -1,5 +1,5 @@ % --- User setup --- -resIndex = 3; % resonance index (1-based) +resIndex = 1; % resonance index (1-based) duration = 0.5; % total ramp time [s] N_points = 300; % number of time points FR_choice = 1; % resonance data choice (1=new, 0=old) @@ -16,14 +16,14 @@ doPlot = true; % show resonance and ramp region figure; subplot(2,1,1); plot(t, B_t, 'b-', 'LineWidth', 1.5); -ylabel('B (G)'); -title('Magnetic field B(t)'); +ylabel('B (G)', 'Interpreter', 'tex'); +title('Magnetic field B(t)', 'Interpreter', 'tex'); subplot(2,1,2); plot(t, a_t, 'r-', 'LineWidth', 1.5); -ylabel('Scattering length a (a_0)'); -xlabel('Time (s)'); -title('Linear ramp of a(t)'); +ylabel('Scattering length a ($a_0$)', 'Interpreter', 'latex'); +xlabel('Time (s)', 'Interpreter', 'tex'); +title('Linear ramp of a(t)', 'Interpreter', 'tex'); %% Function: makeLinearScatteringRamp @@ -134,8 +134,8 @@ if nargin >= 8 && doPlot hold on; plot(B_ramp, a_ramp, 'r-', 'LineWidth', 2, 'DisplayName', 'Ramp a(t)'); xlabel('B (G)'); - ylabel('a (a_0)'); - title(sprintf('Feshbach Resonance #%d at B0=%.3f G (%s side)', resIndex, B0, side)); + ylabel('a ($a_0$)', 'Interpreter', 'latex'); + title(sprintf('Feshbach Resonance #%d at B0=%.3f G (%s side)', resIndex, B0, side), 'Interpreter', 'tex'); legend('Location','best'); grid on; end