set term pdf fontscale 0.75 size 6in,4in
outfile = 'BetaDotB-single.pdf'
set output outfile
set timestamp offset graph 1.12,0.7 rotate font 'Verdana,6'
set label 'plotting_scripts/BetaDotB-single.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set xrange [0:100]
betay = 19.5
gamma = 1./betay
eta = 8.8
etap=0.
dir0='20180919_195851' #0.0
filename0 =  dir0.'/BetaCrossE_turn.dat
A= (0.0)**2/betay*1.e6
char_A =sprintf("%3.1f",A)
t0 = 'a_y = '.char_A.' mm-mrad'

dir3 = '20180919_194519' # 3cm
filename3 =  dir3.'/BetaCrossE_turn.dat
A= (0.03)**2/betay*1.e6
char_A =sprintf("%3.1f",A)
t3 = 'a_y = '.char_A.' mm-mrad'
dir4 = '20180919_194821'  #0.04
filename4 =  dir4.'/BetaCrossE_turn.dat'
A= (0.04)**2/betay*1.e6
char_A =sprintf("%3.1f",A)
t4 = 'a_y = '.char_A.' mm-mrad'
dir2='20180919_134137' #0.02
filename2 =  dir2.'/BetaCrossE_turn.dat'
A= (0.02)**2/betay *1.e6
char_A =sprintf("%3.1f",A)
t2 = 'a_y = '.char_A.' mm-mrad'


dir1='20180919_195008'  #0.01
filename1 =  dir1.'/BetaCrossE_turn.dat'
A= (0.01)**2/betay*1.e6
char_A =sprintf("%3.1f",A)
t1 = 'a_y = '.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_p(t) [ppm]"

f = 1./1.45*1.e6
#set key noenhance
plot filename0 u ($8*1.e6):($14/$12)*f w l t t0,filename1 u ($8*1.e6):($14/$12)*f w l t t1,filename2 u ($8*1.e6):($14/$12)*f w l t t2,filename3 u ($8*1.e6):($14/$12)*f w l t t3,\
filename4 u ($8*1.e6):($14/$12)*f w l t t4
set ylabel "C_p [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