filename = 'FiberHarp.dat'

#set terminal pdf enhanced fontscale 0.5 size 6.0in, 4.0in
#harp = 1
harp_num=sprintf("%i",harp)
#set output 'FiberHarp'.harp_num.'-Position.pdf'
print 'Harp '.harp_num
set label 'Harp '.harp_num at graph 0.1, 1.04
#set multiplot #title filename noenhance
set label 'plotting_scripts/harp_centroid.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
#set size 0.5,0.97
#set origin 0.,0.
set xrange [0:50]
set xlabel 'Time [us]'
set ylabel 'Horizontal displacement [mm]'
plot filename u ($2/100):($1 == harp ? ($3*1000):1/0) w l t '<x>'
#set size 0.5,0.5
unset label
unset ylabel
set xlabel 'frequency [MHz]'
#set origin 0.5,0.
set xrange [0.1:10]
plot filename u 6:($1 == harp ? ($9*2+$10**2)**.5:1/0) w l t '<x> fft'

#plot filename u 6:($1 == harp ? ($11**2+$12**2)**.5:1/0) w l t '<x-width> fft'