summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2016-05-13 12:03:33 -0400
committerEugeniy Mikhailov <evgmik@gmail.com>2016-05-13 12:03:33 -0400
commitb9fe2cd7f27cdf07cbd4c9e127f2a458520b42b4 (patch)
treeddd1db4135e94de92af4d291820e02ec1c93f46b
parent25e8ec6e381a1c4896c83ce3d35af7918132b2da (diff)
downloadnoisy_eit_xmds-b9fe2cd7f27cdf07cbd4c9e127f2a458520b42b4.tar.gz
noisy_eit_xmds-b9fe2cd7f27cdf07cbd4c9e127f2a458520b42b4.zip
noise parameters can be send as arguments now
-rw-r--r--MOR_5_levels_with_doppler_and_propagation/Nlevels_with_MOR.xmds11
1 files changed, 8 insertions, 3 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 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);