set term pdf fontscale 0.75 size 6in,4in
outfile = 'BetaCrossE-single.pdf'
set output outfile
set timestamp offset graph 1.04,0.7 rotate font 'Verdana,6'
set label 'plotting_scripts/BetaCrossE-single.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set xrange [0:400]
betax = 7.8
gamma = 1./betax
eta = 8.8
etap=0.
dir2 = '20180911_200122'
filename2 =  dir2.'/BetaCrossE_turn.dat
A= (0.0 - eta*0.002)**2/betax*1.e6
char_A =sprintf("%3.3f",A)
t2 = 'a_x = '.char_A.' mm-mrad'
dir1 = '20180911_201039'
filename1 =  dir1.'/BetaCrossE_turn.dat'
A= (-0.01 - eta*0.002)**2/betax*1.e6
char_A =sprintf("%3.3f",A)
t1 = 'a_x = '.char_A.' mm-mrad'
dir3='20180911_202317'
filename3 =  dir3.'/BetaCrossE_turn.dat'
A= (0.02 - eta*0.002)**2/betax *1.e6
char_A =sprintf("%3.3f",A)
t3 = 'a_x = '.char_A.' mm-mrad'


dir4='20180912_120122'
filename4 =  dir4.'/BetaCrossE_turn.dat'
A= (0.04 - eta*0.002)**2/betax*1.e6
char_A =sprintf("%3.3f",A)
t4 = 'a_x = '.char_A.' mm-mrad'

set label at graph 0.1, 0.9 '{/Symbol D}p/p= 0.2%'
set xlabel "time [us]
set ylabel "(C_e)/({/Symbol D}p/p) [ppm]"

f = 2/3.e8/1.45*1.e6
#set key noenhance
plot filename3 u ($8*1.e6):($10/$8)*f w l t t3,filename2 u ($8*1.e6):($10/$12)*f w l t t2,filename4 u ($8*1.e6):($10/$12)*f w l t t4,\
filename1 u ($8*1.e6):($10/$12)*f w l t t1
set ylabel "C_e [ppm]"

set autoscale y
plot filename3 u ($8*1.e6):($10/$8*$7)*f w l t t3,filename2 u ($8*1.e6):($10/$12*$7)*f w l t t2,filename4 u ($8*1.e6):($10/$12*$7)*f w l t t4,\
filename1 u ($8*1.e6):($10/$12*$7)*f w l t t1


print outfile