#!/bin/bash
BASE="$1"
BASE=${BASE%.img}
BASE=${BASE%MTT}
BASE=${BASE%Tmax}
BASE=${BASE%rCBF}
BASE=${BASE%rCBV}

for a in MTT Tmax rCBF rCBV; do
	NEWFILE="${BASE}${a}.img_angle.csv"
	[ -e "${NEWFILE}" ] || {
		xmessage "Fehler - Datei ${NEWFILE} nicht gefunden."
		exit 1
	}
done

paste "${BASE}MTT.img_angle.csv" \
      "${BASE}Tmax.img_angle.csv" \
      "${BASE}rCBF.img_angle.csv" \
      "${BASE}rCBV.img_angle.csv" | \
      tr ";" "	" | \
      cut -d"	" -f1,2,4,6,8 | \
      grep -v "^BG" >"${BASE}all.tsv"
{
	echo "set terminal png size 1280,1024 giant x888888 x000000 x222222 xaa0000 x000000 xffff00 xffffff"
	echo "set output '${BASE}all.png'"
	echo "set title \"${BASE%-}\""	
	echo "set ylabel \"MTT/Tmax\""
	echo "set yrange [0:50]"
	echo "set ytics 10"
	echo "set y2label \"rCBV/F\""
	echo "set y2range [0:350]"
	echo "set y2tics 50"
	echo "set bmargin 4"
	echo "set grid ytics"
	echo "set samples 720"
	echo "set xtics rotate 10"
	b=5
	for a in ACPR ACPR ACPR HGZR HGZR ACMR ACMR ACMR ACMR ACMR ACMR ACMR ACMR ACMR VGZR VGZR ACAR ACAR \
	         ACAL ACAL VGZL VGZL ACML ACML ACML ACML ACML ACML ACML ACML ACML HGZL HGZL ACPL ACPL ACPL; do
		echo "set label \"$a\" at $b,1 rotate"
		b=$((b+10))
	done
	echo -n "plot "
	echo -n      "'${BASE}all.tsv' using 1:2 w l lw 5 title \"MTT\" smooth csplines,"
	echo -n      "'${BASE}all.tsv' using 1:3 w l lw 5 title \"Tmax\" smooth csplines,"
	echo -n      "'${BASE}all.tsv' using 1:4 w l lw 5 title \"rCBF\" smooth csplines axes x1y2,"
	echo -n      "'${BASE}all.tsv' using 1:5 w l lw 5 title \"rCBV\" smooth csplines axes x1y2"
	echo
} | gnuplot
