#!/usr/bin/gnuplot
set terminal postscript eps color enhanced solid font "Helvetica,22" size 12.6cm,8cm
set output "../tx.eps"
#set output "tx_.eps"
set border linewidth 1.5
set border 3

#f1 = '2016-07-05/tx_2016_07_05'
#f1 = '2016-07-18/2016-07-18'
#f1 = '2016-07-19/2016-07-19'
f1 = '2016-07-21_tx'
#f2='2016-07-05/sat_spec_2016_07_05'
#f2='2016-07-17_sat_refl/2016-07-17_sat_refl'
#f2 = '2016-07-22_refl'
set samples 5000
set multiplot 
set bar 2
unset key

##### tx diff plots ####
set xtics nomirror
set ytics nomirror

tx=0.15
gamma0 = 7
x0 = 53

lorentz(x) = 1 - tx / ( 4*( (x0-x) / (gamma0) )**2 + 1 ) 

#tx1=0.15
#gamma1 = 7
#x1 = 53
#dis=0.02
#lorentz_dis(x) = 1 - tx1 / ( 4*( (x1-x) / (gamma1) )**2 + 1 ) + dis * (x1-x) / ( 4*( (x1-x) / (gamma1) )**2 + 1 ) 

#lorentz_dis_beta(x) = (b**2 + (4*b**2*Lambda_p**2 - 4*b**2*Lambda_p*cos(phi))*gamma2**2/ ( 4*(x2-x)**2 + (gamma2)**2 ) - ( 8*b**2*Lambda_p*sin(phi)*(x-x2)*gamma2 ) / ( 4*(x2-x)**2 + gamma2**2 ) ) / b**2
lorentz_dis_beta(x) = (1 + (4*Lambda_p**2 - 4*Lambda_p*cos(phi))*gamma2**2/ ( 4*(x2-x)**2 + (gamma2)**2 ) - ( 8*Lambda_p*sin(phi)*(x-x2)*gamma2 ) / ( 4*(x2-x)**2 + gamma2**2 ) ) 
#lorentz_dis_beta(x) = (b**2 + (4*b*Lambda_f*Lambda_p - 4*b*sqrt(Lambda_f*Lambda_p)*cos(phi))*gamma2**2/ ( 4*(x2-x)**2 + (gamma2)**2 ) - ( 8*b*sqrt(Lambda_f*Lambda_p)*sin(phi)*(x-x2)*gamma2 ) / ( 4*(x2-x)**2 + gamma2**2 ) ) / b**2

#lorentz_dis_beta(x) = 1 + (gamma2*sqrt(Lambda_p*Lambda_f/beta))**2 * 1/((x2-x)**2 + (gamma2/2)**2) + 2 *(gamma2*sqrt(Lambda_p*Lambda_f/beta))* 1/((x2-x)**2 + (gamma2/2)**2)* ((x2-x)*sin(phi)-gamma2/2*cos(phi)) 

gamma2 = 6.9
x2 = 46.4
coupling = 0.7
b = sqrt(coupling)
beta = b
phi = 0.01
a = 0.07
Lambda_p = 0.047
Lambda_f = 0.0354

fit lorentz(x) f1 u (($1-190000)/1000*2):3:4 via tx,x0,gamma0
fit lorentz_dis_beta(x) f1 u (($1-190000)/1000*2):3:4 via  x2, gamma2, Lambda_p, phi
#set title 'transmission extinction'
set xlabel 'probe detuning (MHz)'
set ylabel 'transmission (%)'
set key right bottom
set yrange [81:101]
set xrange [30:66]
set xtic 10
set ytics 4
unset key

plot f1 u (($1-190000)/1000*2):($3*100):($4*100) w yerr ps 1 pt 7 lw 3, lorentz_dis_beta(x)*100 lc 1 lw 3 

unset multiplot   
unset output
pause -1