#!/usr/bin/gnuplot

set terminal postscript eps color enhanced solid font "Helvetica,18" size 8.6cm, 6cm
#set terminal epslatex color
set output "fig3_label.eps"

#################################
#######
####GENERAL SETTINGS####
set border linewidth 1.5
set samples 5000
set multiplot 
#set bar 2
unset bar
unset key
#set style fill solid 0.3 noborder

###########################
#####            #####

unset key

#unset xtics
#unset xlabel
#unset ytics
#unset ylabel
unset key

set border 15
#set tmargin at screen 0.53
#set bmargin at screen 0.18
set lmargin at screen 0.20
#set rmargin at screen 0.95

f1 = 'tms.dat'


set ytics 40 nomirror format "%g" 
set xtics 0.5 nomirror

set yrange [0:160]
set xlabel 'critical distance d ({/Symbol m}m)'
set ylabel 'spacing {/Symbol D}{/Symbol n}_{tr}  (MHz)'
unset label

unset title
unset grid 
unset key



#set label 1 "b" at graph 0.08,0.23
#set label 2 "a" at graph 0.18,0.35
#lorentz_T(x) = T0 / ( 4*( (x-x0) / (gamma0) )**2 + 1 )
#fit [-200:200] lorentz_T(x) f3 u 1:($2/max_T_cts*max_T):($3/max_T_cts*max_T) via x0, gamma0, T0

n0=28205.0
R=5.5e-3
lambda=780*10**-9/1.00029
c=299792458
pi=3.141592653589793238
tms(x)=-c/(2*(n0-x)*lambda/2)*(1/pi*acos(1-(n0-x)*lambda/2/R)-1)
fit tms(x) f1 using 1:($2*10**6):($3*10**6) via n0,R
tms_l(x)=-c/(2*(-x+2*R))*(1/pi*acos(1-(-x+2*R)/R)-1)/10**6
set xrange [0:3.5]

plot f1 using ((2*R-(n0-$1)*lambda/2)*10**6):2:3 w yerr  pt 7 lc rgb "red", tms_l(x/10**6) lc rgb "red" 

#################################