From b9fe2cd7f27cdf07cbd4c9e127f2a458520b42b4 Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Fri, 13 May 2016 12:03:33 -0400
Subject: noise parameters can be send as arguments now

---
 .../Nlevels_with_MOR.xmds                                     | 11 ++++++++---
 1 file changed, 8 insertions(+), 3 deletions(-)

(limited to 'MOR_5_levels_with_doppler_and_propagation')

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 ece5454..86aa117 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
@@ -101,6 +101,11 @@
 
       <argument name="EpLo" type="real" default_value="0.05*(2*M_PI*1e6)" /> 
       <argument name="EpRo" type="real" default_value="0.05*(2*M_PI*1e6)" />
+
+			<!-- Light field Noise terms strength -->
+      <argument name="commAmplNoise" type="real" default_value="0" />
+      <argument name="onePhotonPhaseNoise" type="real" default_value="0" />
+      <argument name="twoPhotonPhaseNoise" type="real" default_value="0" />
  
 
 			<!-- Fields detuning in [1/s] -->
@@ -190,9 +195,9 @@
 			double r1 = 2.0*(1.0*rand()/RAND_MAX) - 1.0; // random number -1 .. +1
 			double r2 = 2.0*(1.0*rand()/RAND_MAX) - 1.0; // random number -1 .. +1
 			double r3 = 2.0*(1.0*rand()/RAND_MAX) - 1.0; // random number -1 .. +1
-			double common_ampl_noise_frac = 1+.0*r1;
-			double one_phot_noise_frac = .0*r2;
-			double two_phot_noise_frac = .3*r3;
+			double common_ampl_noise_frac = 1+commAmplNoise*r1;
+			double one_phot_noise_frac = onePhotonPhaseNoise*r2;
+			double two_phot_noise_frac = twoPhotonPhaseNoise*r3;
 			EdL=EdLo * common_ampl_noise_frac * exp(i*one_phot_noise_frac);
 			EdR=EdRo * common_ampl_noise_frac * exp(i*one_phot_noise_frac);
 			EpL=EpLo * common_ampl_noise_frac * exp(i*one_phot_noise_frac) * exp(i*two_phot_noise_frac);
-- 
cgit v1.2.3