set multiplot
set xlabel font 'Verdana,20'
set ylabel font 'Verdana,20'
set xtics font 'Verdana,20'
set ytics font 'Verdana,20'
set origin 0.,0.
set size 1.,0.5
set xlabel "x [mm]"
set ylabel "px [mrad]"
set xrange [-5:5]
set yrange [-0.3:0.3]
plot "outgoing.dat" using ($3*1000):($2>2600 && $2<3900 ? $4*1000 : 1/0) title "horizontal phase space"
set xlabel "turn number"
set ylabel "particles per turn"
set xrange [2600:3900]
set yrange [0:25]
set origin 0.,0.5
plot "runstats.dat" u 2 t "particles extracted per turn" with lines lt rgb "blue", "" u ($3 +$4 +$5 +$6) t "particles lost per turn" with lines lt rgb "red"