From 78668acddb2a2fbda01711c18b3eb552ca715bdd Mon Sep 17 00:00:00 2001
From: "Eugeniy E. Mikhailov" <evgmik@gmail.com>
Date: Mon, 12 Oct 2020 20:37:17 -0400
Subject: redone phase diagram

---
 compass_lin_extrema_vs_theta_output_results.m | 48 +++++++++++++++++++++++----
 1 file changed, 41 insertions(+), 7 deletions(-)

diff --git a/compass_lin_extrema_vs_theta_output_results.m b/compass_lin_extrema_vs_theta_output_results.m
index cfb101b..eb151b8 100644
--- a/compass_lin_extrema_vs_theta_output_results.m
+++ b/compass_lin_extrema_vs_theta_output_results.m
@@ -33,7 +33,9 @@ line_colors= [ ...
 	];
 
 figure(1);
+if (exist('keep','var') && keep)
 clf();
+end
 hold off;
 labels={};
 for i=1:N_detunings-1
@@ -58,21 +60,39 @@ print('compass_lin_sidebands_vs_theta.png')
 hold off;
 
 figure(2)
+if (exist('keep','var') && keep)
 clf();
+end
+%   resonance   -3   -2   -1   0   1   2   3
+%   index        1    2    3   4   5   6   7
 j=4; k= 3;
-j=2; k= 7;
+j=6; k= 7;
+j=4; k= 7;
 % plotting parametric line of a sideband amplitude (j) vs another one (k)
 z = zeros(1,N_angles);
 col = thetas;  % This is the color, vary with x in this case.
-surface([zoom_factor*(transmission_matrix(j,:));zoom_factor*(transmission_matrix(j,:))], ...
-  [zoom_factor*(transmission_matrix(k,:));zoom_factor*(transmission_matrix(k,:))], ...
+peak_1 = (transmission_matrix(j,:));
+peak_2 = (transmission_matrix(k,:));
+peak_3 = (transmission_matrix(6,:));
+x_combo = peak_2./(peak_1+peak_2);
+y_combo = peak_1./(peak_1+peak_2);
+eps=1e-5;
+ind= abs(peak_1+peak_2)<eps;
+x_combo(ind)=1;
+y_combo(ind)=0;
+%x_combo = peak_1;
+%y_combo = peak_2;
+surface([x_combo;x_combo], ...
+  [y_combo;y_combo], ...
   [z;z],[col;col],...
   'facecol','no',...
   'edgecol','interp',...
   'linew',2);
 %plot(zoom_factor*(transmission_matrix(j,:)), zoom_factor*(transmission_matrix(k,:)))
-xlabel( labels{j});
-ylabel( labels{k});
+%xlabel( labels{j});
+%ylabel( labels{k});
+xlabel('Sideband ratio: a_3/(a_0+a_3)');
+ylabel('Sideband ratio: a_0/(a_0+a_3)');
 %clabel('\theta');
 colorbar;
 set(gca, 'ZTick', [0,pi/4, pi/2, 3*pi/4, pi])
@@ -80,8 +100,22 @@ set(gca, 'ZTickLabel', {'0', '\pi/4', '\pi/2', '3\pi/4', '\pi'})
 
 %colorbar('Ticks',[0, pi/2, pi], 'TicksLabels', ['0', '\pi/2', '\pi'])
 %ctickslabel(['0', '\pi/2', '\pi'])
-title('One sideband amplitude vs another for different angles theta');
-print('compass_lin_sidebands_combo_vs_theta.png')
+title('Phase diagram for different angles theta');
+print('compass_lin_sidebands_phase_diagram_vs_theta.png')
+
+figure(3)
+if (exist('keep','var') && keep)
+clf();
+end
+
+plot(thetas, (x_combo))
+set(gca, 'XTick', [0,pi/4, pi/2, 3*pi/4, pi])
+set(gca, 'XTickLabel', {'0', '\pi/4', '\pi/2', '3\pi/4', '\pi'})
+title("Sideband ratio: a_3/(a_0+a_3)");
+xlabel('Angle \theta, between B-field and light propagation direction');
+ylabel('Ratio');
+ylim([0:1]);
+print('compass_lin_sidebands_ratio_vs_theta.png')
 
 
 % vim: ts=2:sw=2:fdm=indent
-- 
cgit v1.2.3