filename='../data/exp_pulse/combined data'
set terminal postscript eps enhanced color font "Helvetica, 18"
set output '../pmax_vs_bandwidth_exp.eps'

timebin=1.0
lifetime=26.3
effiency=0.3
overlap=0.03
rho(x)=x/(overlap*effiency*(timebin/lifetime))

set key at 1.5,1
set title 'Exponential pulse ({/Symbol t} = 1/{/Symbol W})'
set xlabel '{/Symbol W}/{/Symbol G}'
set ylabel 'P_{e,max}'

plot filename i 0 u (1/(column(17)/lifetime)):(rho(column(15))) w lp title 'n=\~5',\
     filename i 2 u (1/(column(17)/lifetime)):(rho(column(15))) w lp title 'n=\~25',\
     filename i 1 u (1/(column(17)/lifetime)):(rho(column(15))) w lp title 'n=\~50',\
     filename i 3 u (1/(column(17)/lifetime)):(rho(column(15))) w lp title 'n=\~100',\
     filename i 4 u (1/(column(17)/lifetime)):(rho(column(15))) w lp title 'n=\~250',\
     filename i 5 u (1/(column(17)/lifetime)):(rho(column(15))) w lp lc 7 title 'n=\~500'

filename2='../data/sq_pulse/combined data'
set terminal postscript eps enhanced color font "Helvetica, 18"
set output '../pmax_vs_bandwidth_sq.eps'

set key at 10,0.5
set title 'Square pulse ({/Symbol t} = 2/{/Symbol W})'
set xlabel '{/Symbol W}/{/Symbol G}'
set ylabel 'P_{e,max}'

plot filename2 i 0 u (1/(column(17)/2/lifetime)):(rho(column(15))) w lp title 'n=\~5',\
     filename2 i 1 u (1/(column(17)/2/lifetime)):(rho(column(15))) w lp title 'n=\~31',\
     filename2 i 2 u (1/(column(17)/2/lifetime)):(rho(column(15))) w lp title 'n=\~45',\
     filename2 i 3 u (1/(column(17)/2/lifetime)):(rho(column(15))) w lp title 'n=\~250',\
     filename2 i 4 u (1/(column(17)/2/lifetime)):(rho(column(15))) w lp lc 7 title 'n=\~500'