From 84063fe7494b6d6a51d9e1db133c04fd18b3ebaa Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Wed, 11 May 2016 14:47:24 -0400
Subject: added ground levels split as executable argument

---
 .../Nlevels_with_MOR.xmds                               | 17 ++++++++++-------
 1 file changed, 10 insertions(+), 7 deletions(-)

diff --git a/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds b/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds
index c58cb11..8c11b5a 100644
--- a/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds
+++ b/MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds
@@ -73,7 +73,6 @@
 				//const double G4=3.0 *(2*M_PI*1e6);
 
 				const double gbc = (2*M_PI)*1e3; // 1/s untits
-				const complex Split = 0;
 				const complex noise = 0;
 				
 				complex Gab, GAB, Gca, GCA, Gcb, GCB;
@@ -102,6 +101,10 @@
 			<!-- Fields detuning in [1/s] -->
       <argument name="delta_one_ph0"  type="real" default_value="0.0" />
       <argument name="delta_two_ph0"  type="real" default_value="0.0" />
+
+			<!-- Ground levels differential split (Mimics magnetic field) [1/s] -->
+      <argument name="Split0"  type="real" default_value="0.0" />
+
 			<!--Pulse duration/width [s] -->
 			<argument name="Pwidth"  type="real" default_value="0.1e-6" />
 			<!--  Atom and cell properties -->
@@ -304,14 +307,14 @@
 					double delta_two_ph = delta_two_ph0;
 					delta_two_ph += abs(Gamma_super)*1*((t)/(t_span/2)); // two photon sweep to mimic EIT shapes
 
-					Gab=Gamma_super+i*(delta_one_ph + delta_two_ph + DopplerShiftPr); // probe right detuning
-					GAB=Gamma_super+i*(delta_one_ph + delta_two_ph + DopplerShiftPr); // probe left detuning
+					Gab=Gamma_super+i*( Split0 + delta_one_ph + delta_two_ph + DopplerShiftPr); // probe right detuning
+					GAB=Gamma_super+i*(-Split0 + delta_one_ph + delta_two_ph + DopplerShiftPr); // probe left detuning
 
-					Gca=Gamma_super-i*(delta_one_ph                + DopplerShiftDr); // drive left detuning
-					GCA=Gamma_super-i*(delta_one_ph                + DopplerShiftDr); // drive right detuning
+					Gca=Gamma_super-i*(-Split0 + delta_one_ph                + DopplerShiftDr); // drive left detuning
+					GCA=Gamma_super-i*( Split0 + delta_one_ph                + DopplerShiftDr); // drive right detuning
 
-					Gcb=gbc+i*( Split + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) );
-					GCB=gbc+i*(-Split + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) );
+					Gcb=gbc+i*( 2*Split0 + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) );
+					GCB=gbc+i*(-2*Split0 + delta_two_ph + (DopplerShiftPr - DopplerShiftDr) );
 
 					complex rba=conj(rab);
 					complex rac=conj(rca);
-- 
cgit v1.2.3