#beta through injection line and once around ring
!rm twiss_mismatch.ps
set term postscript enhanced color 'Times-Roman' 22
set output 'twiss_mismatch.ps'


fx=45./7.7
fy=45./5.4
emit = 34.15 #mm-mrad
e=1
set yrange [0:60]
set xlabel 's [m]'
set ylabel '({/Symbol e b}_y)^{1/2} [mm]'
set y2label '{/Symbol (hd)} [mm]'
#set ylabel '{/Symbol s}_x [mm]'
#set key off
delta = 1.  #delta * 1000 (m to mm) (This is 0.001)

plot 'injection_line_twiss.dat' u 1:($2**.5)*fx w l lw 2 lc 1 lt 1 not,\
 'twiss_file.dat' u ($2+4.1):($4**.5)*fx w l lw 2 lc 1 lt 1 t '({/Symbol e b}_x)^{1/2}',\
 'twiss_file.dat' u ($2+4.1):($8 * 1) w l lw 2 lc 3 lt 1 t '{/Symbol (hd)}',\
 'aperture.dat' u 1:2 i 0:0 w l lt 2 lc -1 lw 2 t 'aperture'
# plot 'injection_line_twiss.dat' u 1:($2*fx**2 + ($6)**2*delta**2)**.5 w l lw 2 lc 1 lt 1 not,\
# 'twiss_file.dat' u ($2+4.1):($4*fx**2+ $8**2*delta**2)**.5 w l lw 2 lc 1 lt 1 t '({/Symbol e b}+{/Symbol (h d)^2)^{1/2}' ,\
# 'aperture.dat' u 1:2 i 0:0 w l lt 2 lc -1 lw 2 t 'aperture'

#plot 'injection_line_twiss.dat' u 1:($3**.5)*fy w l lw 2 lc 2 lt 1 not,\
# 'twiss_file.dat' u ($2+4.1):($6**.5)*fy w l lw 2 lc 2 lt 1 t '{/Symbol s}_y',\
# 'aperture.dat' u 1:2 i 1:1 w l lt 2 lc -1 lw 2 t 'aperture'





#

# 'injection_line_twiss.dat' u 1:($6)*e w l lw 2 lc 3 lt 3 not,\
# 'twiss_file.dat' u ($2+4.1):($8)*e w l lw 2 lc 3 lt 3 t '{/Symbol h}_x'