From 09ed61a30016fb4ac3733b8fcd3f532bd4838db5 Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Wed, 23 Mar 2016 18:50:22 -0400
Subject: better expression forl beamsplitted signals

---
 delta_2photon_noise.xmds | 10 ++++++----
 1 file changed, 6 insertions(+), 4 deletions(-)

diff --git a/delta_2photon_noise.xmds b/delta_2photon_noise.xmds
index 70e9150..f04e1db 100644
--- a/delta_2photon_noise.xmds
+++ b/delta_2photon_noise.xmds
@@ -19,7 +19,7 @@
 			<![CDATA[
 			const complex g=10;
 			const complex gbc=.001;
-			const complex Split = 1;
+			const complex Split = 0;
 			const double tset=5;
 			const double alpha = 0;
 
@@ -120,7 +120,7 @@
 			<group>
 				<sampling basis="delta_2photon" initial_sample="yes">
 					<dependencies>density_matrix </dependencies>
-					<moments>eit_signal mor_signal rab_rlOut rab_imOut rca_rlOut rca_imOut raaOut rbbOut rccOut noiseOut</moments>
+					<moments>eit_signal mor_signal rab_rlOut rab_imOut rca_rlOut rca_imOut raaOut rbbOut rccOut noiseOut PDl PDr</moments>
 					<![CDATA[
 					rab_rlOut = rab.Re();
 					rab_imOut = rab.Im();
@@ -130,8 +130,10 @@
 					rbbOut = rbb.Re();
 					rccOut = rcc.Re();
 					noiseOut = noise.Re();
-					eit_signal = abs(rab + rAB) + abs(rca + rCA);
-					mor_signal = abs(rab + rAB) - abs(rca + rCA);
+					PDl = norm( 1/sqrt(2)*(rab + rAB) + i/sqrt(2)*(rca + rCA) );
+					PDr = norm( 1/sqrt(2)*(rab + rAB) - i/sqrt(2)*(rca + rCA) );
+					eit_signal = PDl + PDr;
+					mor_signal = PDl - PDr;
 					]]>
 				</sampling>
 			</group>
-- 
cgit v1.2.3