rpirena/hour_rate.gpt

# -*- gnuplot -*-

set ylabel "[counts/hour]" 
set yrange [ 0:100 ]
set locale "de_DE.UTF-8"

t0 = `date +%s -d 0`
t(x) = (x-t0)/3600.0

U=""
if (ARG1 eq "MDDM") { U="-v U=MDDM"; ARG1=ARG2; ARG2=ARG3; ARG3=ARG4 }

if (ARG1 ne "") { fn=ARG1 }
T=3600
if (ARG2 ne "") { T=ARG2 }
S = 3600
if (ARG3 ne "") { S=ARG3 }

plot "<./hour_rate.awk ".U." T=".T." ".fn \
     u (t($1)):($2*S):($3*S) tit "Center cm^2" w errorl pt 7, \
  "" u (t($1)):($4*S):($5*S) tit "Vertical" w errorl pt 7, \
  "" u (t($1)):($6*S):($7*S) tit "Horizontal"w errorl pt 7