set term pdf enhanced fontscale 0.75 size 6.in,4.in
output_name = 'CaptureVsEmit.pdf'
set output output_name
data_file = '../scan_emittance/scan_emit_data.dat'
reset
f(x) = a*x+b
a=1.
b=1.
set title '1/(Number captured) vs Emittance'
set yrange [0:]
set xlabel 'emittance [mm-mrad]'
set ylabel '(number captured)^{-1}'
set label at graph 0.1,0.95 'Kicker at 78%'
set label at graph 0.1,0.89 'Quad at 18kV'
set label at graph 0.1,0.83 '{/Symbol e}_x =  {/Symbol e}_y'
set label '/plotting_scripts/CaptureVsEmit.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set label 'mytest/scan_emittance/' at graph 0.8,1.02 font 'Verdana,6' noenhance
set timestamp offset graph 1.04,0.7 rotate font 'Verdana,6'
fit  f(x)  data_file u ($1*2.5/1.e-6):(1/$3) via a,b
plot data_file u ($1*2.5/1.e-6):(1/$3) pt 7 ps 1.5 not,f(x) t 'constant/{/Symbol e}'
print 'write '.output_name
#clear
#unset term
#set term pdf enhanced fontscale 0.75 size 6.in,4.in
#output_name = 'CaptureVsEmit.pdf'
#set output output_name
reset
set title 'Number captured vs Emittance'
set xlabel 'emittance [mm-mrad]'
set ylabel 'number captured'
set label at graph 0.1,0.95 'Kicker at 78%'
set label at graph 0.1,0.89 'Quad at 18kV'
set label at graph 0.1,0.83 '{/Symbol e}_x =  {/Symbol e}_y'
a=1
#f(x) = a + b/x + c/x**2
#fit  f(x)  'scan_emit_data.dat' u ($1*2.5/1.e-6):($3) via a,b,c
plot data_file u ($1*2.5/1.e-6):($3) pt 7 ps 1.5 not
print 'write '.output_name
unset term