set term postscript color enhanced 'Times-Roman' 22
!rm inj_twiss.ps
set output 'inj_twiss.ps'
set xlabel 'z[m]'
set ylabel '{/Symbol b}_y [m]'
set xrange [0:5]
set yrange [0:60]
plot 'injection_line_twiss_old-inf_2m.dat' u 1:3 w l lt 1 lc 1 lw 2 t 'old inf, {/Symbol b}_y(inf)=10m',\
'injection_line_twiss_23mm_gradient015.dat' u 1:3 w l lt 1 lc 3 lw 2 t 'prop inf, {/Symbol b}_y(inf)=10m'