From a8dbe52e26b350fce7e88b6f7c90ee8beb4fa054 Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Mon, 28 Mar 2016 17:16:33 -0400
Subject: new changes for noise traces

---
 delta_2photon_noise.xmds | 10 +++++-----
 pp.m                     | 13 +++++++++++++
 2 files changed, 18 insertions(+), 5 deletions(-)

diff --git a/delta_2photon_noise.xmds b/delta_2photon_noise.xmds
index 7723082..2cb3a52 100644
--- a/delta_2photon_noise.xmds
+++ b/delta_2photon_noise.xmds
@@ -21,7 +21,7 @@
 			const complex gbc=.001;
 			const complex Split = .1;
 			const double tset=5;
-			const double alpha = .0;
+			const double alpha = 0.001;
 
 			complex Gab, Gca, Gcb, GCB;
 			complex Ep,Ed,Epc,Edc;
@@ -41,7 +41,7 @@
 		</transverse_dimensions>
 	</geometry>
 
-	<driver name="multi-path" paths="2" />
+	<driver name="multi-path" paths="16" />
 
 	<noise_vector name="drivingNoise" kind="wiener" type="real" method="dsfmt" >
     		<components>Eta</components>
@@ -80,10 +80,10 @@
 				Epc = conj(Ep);
 				Edc = conj(Ed);
 
-				Gab=g+i*(delta_d+delta_2photon+noise);
+				Gab=g+i*(delta_d+delta_2photon+1*noise);
 				Gca=g-i*(delta_d);
-				Gcb=gbc+i*(Split + delta_2photon+noise);
-				GCB=gbc+i*(-Split + delta_2photon+noise);
+				Gcb=gbc+i*(Split + delta_2photon+0*noise);
+				GCB=gbc+i*(-Split + delta_2photon+0*noise);
 
 
 				rba=conj(rab);
diff --git a/pp.m b/pp.m
index 150cb01..f59f19c 100644
--- a/pp.m
+++ b/pp.m
@@ -1,8 +1,21 @@
 delta_2photon_noise
 
+figure(1)
 plot(delta_2photon_1,  mean_mor_signal_1(end,:) )
 hold all
 plot(delta_2photon_1,  mean_eit_signal_1(end,:) )
 
 legend( {'mor', 'eit'} )
 
+figure(2)
+plot(delta_2photon_1,  stderr_mor_signal_1(end,:) )
+hold all
+plot(delta_2photon_1,  stderr_eit_signal_1(end,:) )
+
+legend( {'mor', 'eit'} )
+
+figure(3)
+plot(delta_2photon_1, stderr_eit_signal_1(end,:)./stderr_mor_signal_1(end,:) )
+hold all
+
+legend( {'noise eit/mor'} )
-- 
cgit v1.2.3