summaryrefslogtreecommitdiff
path: root/compass_lin.m
diff options
context:
space:
mode:
authorEugeniy E. Mikhailov <evgmik@gmail.com>2020-09-21 19:19:27 -0400
committerEugeniy E. Mikhailov <evgmik@gmail.com>2020-09-21 19:19:27 -0400
commitcccae19837f97ab4442ee1fa393c40cae9eadecf (patch)
tree08acd19b52fa3003ba8cc1df4eb6c0a6e9a04d2a /compass_lin.m
parented45e1268bddae21bd88f4ebbb9993767fbeb529 (diff)
downloadmulti_mode_eit-cccae19837f97ab4442ee1fa393c40cae9eadecf.tar.gz
multi_mode_eit-cccae19837f97ab4442ee1fa393c40cae9eadecf.zip
added elipticity tracking for lin polarized compass
Diffstat (limited to 'compass_lin.m')
-rw-r--r--compass_lin.m20
1 files changed, 15 insertions, 5 deletions
diff --git a/compass_lin.m b/compass_lin.m
index 0ebd20c..18316f8 100644
--- a/compass_lin.m
+++ b/compass_lin.m
@@ -70,10 +70,14 @@ phi=pi*2/8;
theta=0/4;
theta=pi/2;
%theta=65/180*pi;
+
%small ellipticity angle psi (0 to pi/2)
% 0 will give linear polarization
-% pi/2 circular polarization
-psi=pi/2;
+% pi/4 circular polarization
+%psi_el=pi/4;
+psi_el=.2*pi/4;
+%psi_el=0*pi/4;
+
@@ -97,11 +101,17 @@ for phi=phis;
wm=wd-(wp-wd);
light_positive_freq=[ wp, wd, wp-wd];
- % we define light as linearly polarized
- % where phi is angle between light polarization and axis x
- [E_field_lab_pos_freq.x, E_field_lab_pos_freq.y] = rotXpolarization(phi, E_field_lab_pos_freq.linear);
+ % we define light as linearly polarized along x
+ E_field_lab_pos_freq.x = E_field_lab_pos_freq.linear;
+ E_field_lab_pos_freq.y = 0;
+ % now we add small elasticity
+ E_field_lab_pos_freq.y=sin(psi_el)*E_field_lab_pos_freq.x*(1i); % order is important
+ E_field_lab_pos_freq.x=cos(psi_el)*E_field_lab_pos_freq.x;
+ % set phi angle between light polarization and axis x
+ [E_field_lab_pos_freq.x, E_field_lab_pos_freq.y] = rotLinPolarization(phi, E_field_lab_pos_freq.x, E_field_lab_pos_freq.y);
E_field_lab_pos_freq.z=E_field_zero;
+
% now we transfor x,y,z, to x',y', and z' with respect to magnetic field az z' axis
E_field_pos_freq=xyz_lin2atomic_axis_polarization(theta, E_field_lab_pos_freq);