#!/usr/bin/gnuplot
#set terminal qt size 550,362 enhanced font 'Arial,18' persist
# Line width of the axes
set border linewidth 1.5

# solution for epsilon:
lambda = 0.03
tx_fit(x) = lambda*(1-lambda) * 4 / (1+x)
lambda_fix = 0.033
tx(x) = lambda_fix*(1-lambda_fix) * 4 / (1+x)
## input file ##

infile_extinction_data_decaying='tx_safwm_full_analysis/tx_safwm_summary.dat'

cw_extinction='cw_extinction.dat'


#### Fit lambda ####

#fit tx_fit(x) infile_extinction_data_decaying u (26.23/$1):($3) via lambda
#################
unset log                              # remove any log-scaling
unset label                            # remove any previous labels

set terminal postscript eps color enhanced solid font "Helvetica,22" size 12.6cm,8cm
set output "fig_tx.eps"

set ytics 3 nomirror format "%g" 
set yr [0.0:13]
set xtics autofreq nomirror
set xr [-0.2:6.5]
#set xr [4:14]
set xlabel "Relative bandwidth {/Symbol G}_p / {/Symbol G}_0"
set ylabel "Extinction {/Symbol e} (%)"
set border 3
set key right top
unset grid
unset bars
#unset key
set key invert
unset key
 
plot tx(x)*100 lw 3 lc rgb "blue" t "Eq.7 with {/Symbol L}=3.3%", \
    infile_extinction_data_decaying u (26.23/$1):(100*$3):(($2*26.23/($1*$1))):(100*$4)  w xyerrorbars pt 7 ps 2 lc rgb "red" lw 3  t "experiment",\
    cw_extinction u 1:(100*$2):(100*$3) w yerrorbars pt 13 ps 2.5 lc rgb "black" lw 3

unset output

#pause -1