diff options
| -rw-r--r-- | delta_2photon_noise.xmds | 10 |
1 files 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> |