#!/usr/bin/gnuplot

infile1 = 'tmp'
#set terminal postscript eps color enhanced solid font "Helvetica,26" size 16cm,10cm
#set term pdf
#set term pngcairo
#set output "tmp.pdf"

tx=0.15
gamma0 = 7
x0 = 40

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

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 ) ) 
 

gamma2 = 6.9
x2 = 40
phi = 0.01
a = 0.07
Lambda_p = 0.047


fit lorentz(x) infile1 u (($1-190000)/1000*2):3:4 via tx,x0,gamma0
#fit lorentz_dis_beta(x) infile1 u (($1-190000)/1000*2):3:4 via  x2, gamma2, Lambda_p, phi


set title "transmission of a faint coherent field "
set xlabel "detuning from bare atomic resonance (MHz)"
set ylabel 'transmission (%)'
set key right bottom
set pointsize 1.0
set xrange [-4:80]
set yrange [75:101]
#set label 1 sprintf("Lorentz+Dispersion fit:") at graph 0.02,0.45
set label 2 sprintf("extinction = %.3g\n", tx*100) at graph 0.02,0.4
#set label 2 sprintf("Lambda = %.3g\n", Lambda_p*100) at graph 0.02,0.4
#set label 3 sprintf("phi = %.3g\n", phi) at graph 0.02,0.35
set label 4 sprintf("gamma = %.3g MHz\n", gamma0) at graph 0.02,0.3
set label 5 sprintf("res = %.3g MHz\n", x0) at graph 0.02,0.25


plot infile1 u (($1-190000)/1000*2):($3*100):($4*100) w yerr ps 0.8 pt 7 lw 2 not ,  lorentz(x)*100 lc 3 lw 2 t 'lorentz' 

unset output
pause -1