unset multiplot
set term pdf fontscale 0.6 size 4.95in,3.3in
outfile = 'compiled_data_collimat.pdf'
set output outfile


set timestamp offset graph 1.04,0.7 rotate font 'Verdana,6'
set label 'plotting_scripts/compiled_data_collimat.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance

set yrange [0:3]
set ylabel "captured [%]"
set xlabel 'collimator radius [mm]'
set label at graph 0.1,0.95 'Inflector tilt = -2 mrad, Injection angle = +2 mrad'

filename1 = 'compilation_collimat.dat'
set key bottom right Left font ',8' noenhance
set yrange [0:2]
set ytics 0,0.2
set grid
plot \
filename1  u 27:($11 < 0.021  ? ($18*100/$2):1/0)  i 0:0 w lp ps 0.5 pt 5  lc 1 t '204 G-circular'   ,\
''  u 27:($11 > 210.e-4  ? ($18*100/$2):1/0)  i 2:2 w lp ps 0.5 pt 6 lc 2 t '250 G-elliptical' ,\
''  u 27:($11 > 210.e-4  ? ($18*100/$2):1/0)  i 3:3 w lp ps 0.5 pt 7 lc 3 t '250 G-circular-100k'   
#''  u 27:($11 > 210.e-4  ? ($18*100/$2):1/0)  i 0:0 w lp ps 0.5 pt 8 lc 4 t '250 G'  ,\
#''  u 27:($11 > 210.e-4  ? ($18*100/$2):1/0)  i 1:1 w lp ps 0.5 pt 9 lc 5 t '250 G-elliptical' ,\




set ylabel "{/Symbol s}_p [%]"
set key Left left font ',8' noenhance
set yrange [0.0:0.12]
set ytics 0,0.02
plot \
filename1  u ($27):($11 < 210.e-4  ? ($24*100):1/0) i 0:0  w lp ps 0.5 pt 5 lc 1 t '204 G-circular'   ,\
filename1  u ($27):($11 > 210.e-4  ? ($24*100):1/0)  i 2:2 w lp ps 0.5 pt 6 lc 2 t '250 G-elliptical'  ,\
filename1  u ($27):($11 > 210.e-4  ? ($24*100):1/0)  i 3:3 w lp ps 0.5 pt 7 lc 3 t '250 G-circular-100k'   
#filename1  u ($27):($11 > 210.e-4  ? ($24*100):1/0)  i 0:0 w lp ps 0.5 pt 8 lc 4 t '250 G'   ,\
#filename1  u ($27):($11 > 210.e-4  ? ($24*100):1/0)  i 1:1 w lp ps 0.5 pt 9 lc 5 t '250 G-elliptical',\


set ylabel "<{/Symbol D}p/p> [%]"

set key right Left  font ',8' noenhance
set yrange [0:0.1]

plot \
filename1  u ($27):($11 < 210.e-4  ? ($25*100):1/0)  i 0:0 w lp ps 0.5 pt 5 lc 1 t '204 G-circular'   ,\
filename1  u ($27):($11 > 210.e-4  ? ($25*100):1/0)  i 2:2 w lp ps 0.5 pt 6 lc 2 t '250 G-elliptical' ,\
filename1  u ($27):($11 > 210.e-4  ? ($25*100):1/0)  i 3:3 w lp ps 0.5 pt 7 lc 3 t '250 G-circular-100k' 
#filename1  u ($27):($11 > 210.e-4  ? ($25*100):1/0)  i 0:0 w lp ps 0.5 pt 8 lc 4 t '250 G'   ,\
#filename1  u ($27):($11 > 210.e-4  ? ($25*100):1/0)  i 1:1 w lp ps 0.5 pt 9 lc 5 t '250 G-elliptical',\

print outfile