From f3c4ff0d1a6392b7ef82b673901d50076716faf9 Mon Sep 17 00:00:00 2001
From: Martin Craig <martin.craig@eng.ox.ac.uk>
Date: Tue, 9 Oct 2018 13:05:03 +0100
Subject: [PATCH] Added some useful output

---
 asl_calib.in | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/asl_calib.in b/asl_calib.in
index 180a330..b0b6949 100755
--- a/asl_calib.in
+++ b/asl_calib.in
@@ -587,6 +587,9 @@ if [ $mode = longtr ]; then
 
     # calculate M0_ref value
     Moval=`fslstats $temp_calib/calib -M` # this is Mz of CSF
+    echo "Signal in reference tissue: $Moval" >> $log
+    corr=`echo "(1 - e(- ( $tr - $taq ) / $T1r) )" | bc -l`
+    echo "Correction factor: $corr (tr=$tr, tsq=$taq, t1=$T1r)" >> $log
     Moval=`echo "$Moval / (1 - e(- ( $tr - $taq ) / $T1r) )" | bc -l` #this is now M0 of the reference
     echo "Mz of reference tissue: $Moval" >> $log