#gnuplot -e 'files = "list.dat"; filename = "scraping" '
#'filenames = "20160922_185803 20160922_212434 ..."
set term x11
print 'files = ', files
load files


columns=17
print filenames
nwords = words(filenames)
print 'nwords =', nwords
if(nwords > 13)print 'More than 11 files'
#do for [n=1:nwords] {item(n) = word(filenames,n).'/Beam_moments_tbt.dat'}
#do for [n=1:nwords] {print item(n)}
set key noenhance
#plot item(1) u ($5/1.e-6):16 w l t item(1)
#replot item(2) u ($5/1.e-6):16 w l t item(2)
print filenames
c1=5
c2=16
set ylabel 'Muons remaining'
set xlabel 'Time after injection {/Symbol m}s' 
if(nwords-2 == 11) {plot filenames u (column(c1)/1.e-6):(column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118) +column(135)+column(152)+column(16+8*17)+column(169)) w l not}
if(nwords-2 == 10) {plot filenames u (column(c1)/1.e-6):(column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118) +column(135)+column(152)+column(16+8*17)) w l not}
if(nwords-2 == 9) {plot filenames u (column(c1)/1.e-6):(column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118) +column(135)+column(152)) w l not}
if(nwords-2 == 8) {plot filenames u (column(c1)/1.e-6):(column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118)+column(135)) w l not }
if(nwords-2 == 7) {plot filenames u (column(c1)/1.e-6):(column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118)) w l not} 
if(nwords-2 == 6) {plot filenames u (column(c1)/1.e-6):(column(16) +column(33)+column(50)+column(67)+column(84)+column(101)) w l not} 
if(nwords-2 == 5) {plot filenames u (column(c1)/1.e-6):(column(16) +column(33)+column(50)+column(67)+column(84)) w l not}
if(nwords-2 == 4) {plot filenames u (column(c1)/1.e-6):(column(16) +column(33)+column(50)+column(67)) w l not}
if(nwords-2 == 3) {plot filenames u (column(c1)/1.e-6):(column(16) +column(33)+column(50)) w l not}
if(nwords-2 == 2) {plot filenames u (column(c1)/1.e-6):(column(16) +column(33)) w l not}
if(nwords-2 == 1) {plot filenames u (column(c1)/1.e-6):(column(16)) w l not}
pause -1
clear
reset
if(nwords-2 == 11) {stats filenames u (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118) +column(135)+column(152)+column(16+8*17)+column(16+9*17))}
if(nwords-2 == 10) {stats filenames u (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118) +column(135)+column(152)+column(16+8*17))}
if(nwords-2 == 9) {stats filenames u (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118) +column(135)+column(152))}
if(nwords-2 == 8) {stats filenames u (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118)+column(135))}
if(nwords-2 == 7) {stats filenames u (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118))} 
if(nwords-2 == 6) {stats filenames u (column(16) +column(33)+column(50)+column(67)+column(84)+column(101))} 
if(nwords-2 == 5) {stats filenames u (column(16) +column(33)+column(50)+column(67)+column(84))}
if(nwords-2 == 4) {stats filenames u (column(16) +column(33)+column(50)+column(67))}
if(nwords-2 == 3) {stats filenames u (column(16) +column(33)+column(50))}
if(nwords-2 == 2) {stats filenames u (column(16) +column(33))}
if(nwords-2 == 1) {stats filenames u (column(16))} 


print STATS_min


YMIN=STATS_min
set yrange [YMIN-1:YMIN+100]
set ylabel 'Muons remaining'
set xlabel 'Time after injection [{/Symbol m}s]'
effic = STATS_min/STATS_max
ratio = sprintf("%1.5f",effic)
filename = 'lost muons noscrape norf'
set label filename at graph 0.4,0.9 noenhance
set label 'Fraction stored = '.ratio at graph 0.4,0.85
set label 'plotting_scripts/Stragglers.gnu' at graph 1.02,0.02 rotate left font 'Verdana,8' noenhance
set terminal pdf enhanced fontscale 0.75 size 6.0in, 6.0in
outputfile = filename.'.pdf' 
print outputfile
set output outputfile

if(nwords-2 == 11) {stats filenames u ($5 > 30.e-6 ? (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118) +column(135)+column(152)+column(16+8*17)+column(16+9*17)):1/0)}
if(nwords-2 == 10) {stats filenames u ($5 > 30.e-6 ? (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118) +column(135)+column(152)+column(16+8*17)):1/0)}
if(nwords-2 == 9) {stats filenames u ($5 > 30.e-6 ? (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118) +column(135)+column(152)):1/0)}
if(nwords-2 == 8) {stats filenames u ($5 > 30.e-6 ? (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118)+column(135)):1/0)}
if(nwords-2 == 7) {stats filenames u ($5 > 30.e-6 ? (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)+column(118)):1/0)} 
if(nwords-2 == 6) {stats filenames u ($5 > 30.e-6 ? (column(16) +column(33)+column(50)+column(67)+column(84)+column(101)):1/0)} 
if(nwords-2 == 5) {stats filenames u ($5 > 30.e-6 ? (column(16) +column(33)+column(50)+column(67)+column(84)):1/0)}
if(nwords-2 == 4) {stats filenames u ($5 > 30.e-6 ? (column(16) +column(33)+column(50)+column(67)):1/0)}
if(nwords-2 == 3) {stats filenames u ($5 > 30.e-6 ? (column(16) +column(33)+column(50)):1/0)}
if(nwords-2 == 2) {stats filenames u ($5 > 30.e-6 ? (column(16) +column(33)):1/0)}
if(nwords-2 == 1) {stats filenames u ($5 > 30.e-6 ? $16 : 1/0) } 

diff = STATS_max-STATS_min
if(diff<1) diff=1
print 'STATS_max = ',STATS_max,' STATS_min = ', STATS_min
print 'fraction lost = ', (diff)/STATS_min
set label 'Num Lost after 30{/Symbol m}s < '.sprintf("%1.6f",(diff)*100./STATS_min).'%' at graph 0.4,0.8
replot

#plot for [file in filenames] file.'/Beam_moments_tbt.dat' u ($5/1.e-6):($16) w l t file

#plot for [n=1:2] item(n).'/Beam_moments_tbt.dat' u ($5/1.e-6):16 w l t item(n)