rpirena/www/cgi/rpirena-uM

178 lines
5 KiB
Text
Raw Permalink Normal View History

#! /bin/bash
# -*- sh -*-
cat <<EOT
Content-type: text/html
<html>
<head>
<meta charset="UTF-8" />
<title>RPiRENA Data</title>
</head>
<body>
<h1>Live plots of µMustang on $(hostname)</h1>
<ul>
<li/><a href="rpirena?old">Show</a> the plots already processed (default).
<li/><a href="rpirena?plot">Process</a> new data and show the plots. May time out for big files.
<li/><a href="rpirena?trigger">Start processing</a> data, without waiting.
<li/><a href="rpirena-control">RPiRENA control.</a>
<li/><a href="/data/">Data Files.</a>
</ul>
EOT
exec 2>&1
cd /home/$(hostname)/stephan/rpirena || exit
lf=../../www/data/rpirena-cgi.log
lkf=../../www/data/rpirena-cgi.lock
declare -a QUERY=(${QUERY_STRING//&/ })
ACTION=${QUERY[0]}
[ -z "$ACTION" ] && TARGET="" || TARGET=${QUERY[-1]}
[ "$TARGET" = "$ACTION" ] && TARGET=""
[ -z "$ACTION" ] && ACTION=old
echo "$REMOTE_ADDR $HTTP_USER_AGENT $QUERY_STRING $ACTION $TARGET" >>$lf
[ -f "$lkf" ] && cat "$lkf" 2>/dev/null
[ -z "$TARGET" ] || cat <<EOF
<p/>Actions for the selected old data file:
<a href="rpirena?plot&amp;$TARGET">plot</a>
<a href="rpirena?old&amp;$TARGET">old</a>
<a href="rpirena?trigger&amp;$TARGET">trigger</a>
<a href="rpirena?replot&amp;$TARGET">replot</a>
<a href="rpirena?nomake&amp;$TARGET">nomake</a>
EOF
cat <<EOF
<p/>Beware of the fake hwclock of the Raspberry Pi, the times in the filenames may be fake.
<p/>You may need to reload the page in your browser to see the latest plots.
EOF
if [ "$ACTION" = "trigger" ]
then
export QUERY_STRING="replot&$TARGET"
nohup ../../www/cgi/rpirena > ../../www/data/rpirena-trigger.log 2>&1 &
ACTION="old"
fi
if [ "$ACTION" = "restart" ]
then
echo restart > auto-run
export ACTION="old"
fi
if [ "$ACTION" = "stop" ]
then
echo stop > auto-run
export ACTION="old"
fi
fn=../../www/data/$TARGET
[ -f $fn ] || fn=$(ls -1 ../../www/data/20??-??-??T??:??:??Z.EI | tail -1)
bfn=$(basename $fn .EI)
pfn=../../www/data/$bfn
sfn=${pfn}-spectra.png
hfn=${pfn}-hk.png
cfn=${pfn}-counts.png
REPLOT=$ACTION
[ -f "$fn" -a "$fn" -nt "$sfn" -a "$ACTION" = "plot" ] && REPLOT=replot
if [ -f "$fn" -a "$REPLOT" = "replot" ] && dotlockfile -r 0 $lkf
then
trap "rm $lkf" EXIT
date +"Processing $bfn since %Y-%m-%d %H:%M:%S UTC" > $lkf
date +"Plot $bfn %Y-%m-%d %H:%M:%S" >>$lf
[ "$ACTION" = "nomake" ] || make UNIT=uM ${pfn}.EmV ${pfn}.uMhist ${pfn}.HK ${pfn}.CP >>$lf 2>&1
cat <<EOT | gnuplot >>$lf 2>&1
print "plot spectra"
set title "$(hostname) RPiRENA $bfn" noenh
set xrange [-10:160]
set yrange [-10:160]
set y2range [0.9:*]
set log y2
set ytics nomirror
set y2tics
set xlab "pulse height [mV]"
set ylab "pulse height [mV]"
set y2lab "[counts/mV]"
set grid
set timestamp "%Y-%m-%d %H:%M:%S"
set term png size 1024,768
set out "${sfn}-new"
plot \
"$pfn.EmV" u "A1":"A2" tit "A1 vs A2" with d, \
"$pfn.EmV" u "B1":"B2" tit "B1 vs B2" with d, \
for [c in "A1 A2 B1 B2"] "$pfn.uMhist" u "mV":c axis x1y2 tit c w histeps
unset out
print "plot HK"
set xrange [*:*]
set xtics timedate
set format x "%d.%m.\n%Hh%M" time
set yrange [0:10]
set y2range [*:*]
unset log y2
set xlab "time [UTC]"
set ylab "[V] [10nA]"
set y2lab "[°C]"
set term png
set out "${hfn}-new"
array HN[8] = [ "V_{CORE}", "V_{CC}", "V_{IO}", "V_{PRIM}", "T_{FPGA}", "T_{CPU}", "T_{EBOX}", "I_{BIAS}" ]
plot \
for [c=8:8] "<tail -n +25 ${pfn}.HK" u 2:(column(2+c)/10) axis x1y1 tit HN[c] with l lt c, \
for [c=1:4] "<tail -n +25 ${pfn}.HK" u 2:2+c axis x1y1 tit HN[c] with l lt c lw 2, \
for [c=5:7] "<tail -n +25 ${pfn}.HK" u 2:2+c axis x1y2 tit HN[c] with l lt c lw 2
unset out
print "plot Counts"
set yrange [0.9:*]
set log y
set xlab "time [UTC]"
set ylab "[counts/min]"
unset log y2
set y2lab "p [mbar]"
set y2range [500<*:*<1100] noextend
set term png
set out "${cfn}-new"
array CN[13] = ["A_1", "A_2", "B_1", "B_2", "A", "B", "A_{12}", "B_{12}", "AB", "A_{12}B_{12}", "A_{12}B", "AB_{12}", "L_3"]
plot \
for [c=1:4] "<tail -n +3 ${pfn}.CP" u 19:(column(2+c)*60./\$17) tit "L_1[".(c-1)."]: ".CN[c+0] with l lt c+1 lw 0.5, \
for [c=5:12] "<tail -n +3 ${pfn}.CP" u 19:(column(2+c)*60./\$17) tit "L_2[".(c-5)."]: ".CN[c+0] with steps lt c-3 lw 1, \
"<tail -n +3 ${pfn}.CP" u 19:18 axis x1y2 notit "pressure" with l lt 5 lw 2, \
"<tail -n +3 ${pfn}.CP" u 19:18 axis x1y2 tit "pressure" with l lt 1 lw 0.8
unset out
EOT
mv -v $sfn-new $sfn >>$lf 2>&1
mv -v $hfn-new $hfn >>$lf 2>&1
mv -v $cfn-new $cfn >>$lf 2>&1
rm -f $lkf
trap EXIT
fi
cat <<EOT
<h3>File: $bfn</h3>
<ul>
<li><a href="/data/$bfn.CP">Count Rates</a>, <a href="/README.txt">README</a>
<li><a href="/data/$bfn.HK">Housekeeping</a>
<li><a href="/data/$bfn.EI">Raw Data</a>
</ul>
<p/><img src="/data/$bfn-spectra.png"/>
<p/><img src="/data/$bfn-hk.png"/>
<p/><img src="/data/$bfn-counts.png"/>
<h2>Older Spectra:</h2>
<ul>
EOT
cd ../../www/data
ls -r1 20??-??-??T??:??:??Z.EI | xargs -i echo '<li> <a href="rpirena?old&amp;{}">{}</a>'
cat <<EOT
</ul>
</p><img src="/data/IMG_6519-1.jpg"/>
</p><img src="/data/IMG_6518-1.jpg"/>
</body>
</html>
EOT