set term pdf fontscale 0.75 size 6in,4in
outfile = 'EfieldCorrectionSummary.pdf'
set output outfile
set timestamp offset graph 1.04,0.7 rotate font 'Verdana,6'
set label  'plotting_scripts/EfieldCorrectionSummary.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
filename='summary.txt'
gamma=29.3
r0=7112.
b = (1.-1./gamma**2)**.5
f=2*b/r0**2*1.e6
set key bottom right Left font ',8' box
set xlabel '2n(1-n){/Symbol b}^2<x_e^2>/R_0^2 [X 10^6]
set ylabel '<C_e> [ppm]'
set xrange [0.26:0.5]
set yrange [0.26:0.5]
set xtics 0.26,0.04
set ytics 0.26,0.04
plot filename u ($1*($10*(1-$10))*f):3 i 1:1 w p ps 1 pt 5 t 'inj n=0.1058', \
'' u ($1*($10*(1-$10))*f):3 i 2:2 w p ps 1 pt 6 t 'inj n=0.1445',\
'' u ($1*($10*(1-$10))*f):3 i 3:3 w p ps 1 pt 7 t 'inj n=0.185',\
'' u ($1*($10*(1-$10))*f):3 i 0:0 w p ps 1 pt 8 t 'dir n=0.1158',\
'' u ($1*($10*(1-$10))*f):3 i 5:5 w p ps 1 pt 11 t 'dir n=0.1058 - scraping',\
'' u ($1*($10*(1-$10))*f):3 i 4:4 w p ps 1 pt 9 t 'dir n=0.1158 b1'

#pause -1

reset
g=1./2./r0**2 * 1.e6
set key bottom right Left font ',8' box
set timestamp offset graph 1.04,0.7 rotate font 'Verdana,6'
set label  'plotting_scripts/EfieldCorrectionSummary.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set xlabel 'n<y^2>/(2 R_0^2) [X 10^6]'
set ylabel '<C_p> [ppm]'
set xtics 0.14,0.04
set ytics 0.14,0.04
plot filename u ($8*g*($10)):6 i 1:1 w p ps 1 pt 5 t 'inj - 0.1058', \
'' u ($8*g*($10)):6 i 2:2 w p ps 1 pt 6 t 'inj-0.1445' ,\
'' u ($8*g*($10)):6 i 3:3 w p ps 1 pt 7 t 'inj-0.185' ,\
'' u ($8*g*($10)):6 i 0:0 w p ps 1 pt 8 t 'dir-0.1445',\
'' u ($8*g*($10)):6  i 5:5 w p ps 1 pt 11 t 'dir n=0.1058 - scraping',\
'' u ($8*g*($10)):6 i 4:4 w p ps 1 pt 9 t 'dir -b1'

print outfile