From 99a5d20d34566165f7ba874493c9a163e2e5dc73 Mon Sep 17 00:00:00 2001
From: Vesa Oikonen <vesoik@utu.fi>
Date: Tue, 27 Jun 2023 10:09:41 +0300
Subject: [PATCH] final

---
 time_frames/frames_and_delay/frames.sif |   8 +-
 time_frames/frames_and_delay/sim1.sh    | 290 +++++++-----------------
 2 files changed, 88 insertions(+), 210 deletions(-)

diff --git a/time_frames/frames_and_delay/frames.sif b/time_frames/frames_and_delay/frames.sif
index 2bc8b4f..b03d610 100644
--- a/time_frames/frames_and_delay/frames.sif
+++ b/time_frames/frames_and_delay/frames.sif
@@ -1,4 +1,4 @@
-26/5/2018 13:54:12 10 2 1
+26/5/2018 13:54:12 16 2 1
 0 5
 5 10
 10 15
@@ -9,3 +9,9 @@
 35 40
 40 50
 50 60
+60 70
+70 80
+80 90
+90 100
+100 120
+120 140
diff --git a/time_frames/frames_and_delay/sim1.sh b/time_frames/frames_and_delay/sim1.sh
index ce1a576..e2967c4 100644
--- a/time_frames/frames_and_delay/sim1.sh
+++ b/time_frames/frames_and_delay/sim1.sh
@@ -1,6 +1,6 @@
 #!/bin/bash
 #: Title      : frames_and_delay/sim1
-#: Date       : 2023-06-20
+#: Date       : 2023-06-27
 #: Author     : "Vesa Oikonen" <vesa.oikonen@utu.fi>
 #: Options    : None
 
@@ -17,7 +17,7 @@ printf "=====================================================================\n"
 printf "Making 66 s radiowater BTAC with 0.1 s sample intervals \n"
 printf "=====================================================================\n"
 
-fit2dat -c=0,66,0.1 radiowater.fit sim1/blood.tac
+fit2dat -c=0,180,0.1 radiowater.fit sim1/blood.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
 
 
@@ -75,17 +75,17 @@ if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
 tac4frpl sim1/true_idif05.tac sim1/bar_true_idif05.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
 
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif00.svg -l sim1/blood.tac sim1/bar_true_idif00.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif00.svg -l sim1/blood.tac sim1/bar_true_idif00.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif01.svg -l sim1/blood01.tac sim1/bar_true_idif01.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif01.svg -l sim1/blood01.tac sim1/bar_true_idif01.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif02.svg -l sim1/blood02.tac sim1/bar_true_idif02.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif02.svg -l sim1/blood02.tac sim1/bar_true_idif02.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif03.svg -l sim1/blood03.tac sim1/bar_true_idif03.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif03.svg -l sim1/blood03.tac sim1/bar_true_idif03.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif04.svg -l sim1/blood04.tac sim1/bar_true_idif04.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif04.svg -l sim1/blood04.tac sim1/bar_true_idif04.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif05.svg -l sim1/blood05.tac sim1/bar_true_idif05.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_true_idif05.svg -l sim1/blood05.tac sim1/bar_true_idif05.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
 
 printf "=====================================================================\n"
@@ -121,15 +121,15 @@ if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
 tac4frpl sim1/tacmove_idif05.tac sim1/bar_tacmove_idif05.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
 
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif01.svg -l sim1/blood01.tac sim1/bar_tacmove_idif01.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif01.svg -l sim1/blood01.tac sim1/bar_tacmove_idif01.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif02.svg -l sim1/blood02.tac sim1/bar_tacmove_idif02.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif02.svg -l sim1/blood02.tac sim1/bar_tacmove_idif02.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif03.svg -l sim1/blood03.tac sim1/bar_tacmove_idif03.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif03.svg -l sim1/blood03.tac sim1/bar_tacmove_idif03.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif04.svg -l sim1/blood04.tac sim1/bar_tacmove_idif04.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif04.svg -l sim1/blood04.tac sim1/bar_tacmove_idif04.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif05.svg -l sim1/blood05.tac sim1/bar_tacmove_idif05.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_tacmove_idif05.svg -l sim1/blood05.tac sim1/bar_tacmove_idif05.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
 
 
@@ -176,15 +176,15 @@ if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
 
 printf "\nPlot these with AIF BTACs\n"
 
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif01.svg -l sim1/blood01.tac sim1/trad_mid_idif01.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif01.svg -l sim1/blood01.tac sim1/trad_mid_idif01.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif02.svg -l sim1/blood02.tac sim1/trad_mid_idif02.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif02.svg -l sim1/blood02.tac sim1/trad_mid_idif02.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif03.svg -l sim1/blood03.tac sim1/trad_mid_idif03.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif03.svg -l sim1/blood03.tac sim1/trad_mid_idif03.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif04.svg -l sim1/blood04.tac sim1/trad_mid_idif04.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif04.svg -l sim1/blood04.tac sim1/trad_mid_idif04.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -x2=65 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif05.svg -l sim1/blood05.tac sim1/trad_mid_idif05.tac
+tac2svg -x2=125 -legend=n -mt="" sim1/sim_frames_and_peak_trad_mid_idif05.svg -l sim1/blood05.tac sim1/trad_mid_idif05.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
 
 printf "\nInterpolate into original frame durations (same result as above therefore not plotted) \n"
@@ -207,200 +207,72 @@ printf "Calculate AUC from zero to frame end times. \n"
 printf "Plot its error at each frame end time.\n"
 printf "=====================================================================\n"
 
-REM interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/blood.tac sim1/auc.tac
-
-interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/blood01.tac sim1/auc01.tac
-
-tacfr2x sim1/true_idif01.tac sim1/steps.tac
-interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/steps.tac sim1/temp.tac
-taccalc sim1/temp.tac - sim1/auc01.tac sim1/aucerror01.tac
-
-tacfr2x sim1/tacmove_idif01.tac sim1/steps.tac
-interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/steps.tac sim1/temp.tac
-taccalc sim1/temp.tac - sim1/auc01.tac sim1/temp2.tac
-tacadd sim1/aucerror01.tac sim1/temp2.tac
-
-interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/trad_mid_idif01.tac sim1/temp.tac
-taccalc sim1/temp.tac - sim1/auc01.tac sim1/temp2.tac
-tacadd sim1/aucerror01.tac sim1/temp2.tac
-
-REM tacfr2x sim1/trad_frame_idif01.tac sim1/steps.tac
-REM interpol -i -x=5,10,15,20,25,30,35,40,60 sim1/steps.tac sim1/temp.tac
-REM taccalc sim1/temp.tac - sim1/auc01.tac sim1/temp2.tac
-REM tacadd sim1/aucerror01.tac sim1/temp2.tac
-
-tacren sim1/aucerror01.tac 1 True IDIF
-tacren sim1/aucerror01.tac 2 tacmove
-tacren sim1/aucerror01.tac 3 Trad
-REM tacren sim1/aucerror01.tac 4 Trad frame
-
-tac2svg -x2=65 -legend=y -mt="" -yt="s*kBq/mL" sim1/sim_frames_and_peak_aucerror01.svg -l sim1/aucerror01.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-
-
-
-
-exit 0
-
-printf "=====================================================================\n"
-printf "Assuming that BTAC is derived from PET image with frame times given in SIF \n"
-printf "=====================================================================\n"
-
-simframe blood.tac frames.sif idif00.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-simframe blood01.tac frames.sif idif01.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-simframe blood02.tac frames.sif idif02.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-simframe blood03.tac frames.sif idif03.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-simframe blood04.tac frames.sif idif04.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-simframe blood05.tac frames.sif idif05.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-
-printf "\n Plot original BTAC and its framed version with bars \n"
-tac4frpl idif00.tac blood2frpl.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -y2=460 -legend=n -mt="" sim_frames_and_delay_btac.svg -l blood2frpl.tac blood.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-
-printf "=====================================================================\n"
-printf "Move framed IDIF back in time \n"
-printf "=====================================================================\n"
-
-tactime idif01.tac -1 idif01dc.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tactime idif02.tac -2 idif02dc.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tactime idif03.tac -3 idif03dc.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tactime idif04.tac -4 idif04dc.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tactime idif05.tac -5 idif05dc.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-
-printf "=====================================================================\n"
-printf "And interpolate to original framing \n"
-printf "=====================================================================\n"
-
-tacmove 
-
-
-
-exit 0
-
-
-printf "\nPlot delay time matched original BTAC and simulated IDIF together \n"
-tac2svg -legend=n -mt="" idif00.svg -l blood.tac -b idif00.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -legend=n -mt="" idif01.svg -l blood01.tac -b idif01.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -legend=n -mt="" idif02.svg -l blood02.tac -b idif02.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -legend=n -mt="" idif03.svg -l blood03.tac -b idif03.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -legend=n -mt="" idif04.svg -l blood04.tac -b idif04.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -legend=n -mt="" idif05.svg -l blood05.tac -b idif05.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -legend=n -mt="" idif06.svg -l blood06.tac -b idif06.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -legend=n -mt="" idif07.svg -l blood07.tac -b idif07.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -legend=n -mt="" idif08.svg -l blood08.tac -b idif08.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -legend=n -mt="" idif09.svg -l blood09.tac -b idif09.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-tac2svg -legend=n -mt="" idif10.svg -l blood10.tac -b idif10.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-
-printf "\nPlot all simulated IDIFs together \n"
-tac2svg -y2=460 -legend=n -mt="" sim_frames_and_peak_idif.svg idif??.tac
+DN=01
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/blood$DN.tac sim1/auc$DN.tac
+tacfr2x sim1/tacmove_idif$DN.tac sim1/steps.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/steps.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc$DN.tac sim1/aucerror$DN.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/trad_mid_idif$DN.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc$DN.tac sim1/temp2.tac
+tacadd sim1/aucerror$DN.tac sim1/temp2.tac
+tacren sim1/aucerror$DN.tac 1 Bar
+tacren sim1/aucerror$DN.tac 2 Mid_point
+tac2svg -x2=125 -legend=y -mt="" -yt="Error in AUC (s*kBq/mL)" sim1/sim_frames_and_peak_aucerror$DN.svg -l sim1/aucerror$DN.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+DN=02
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/blood$DN.tac sim1/auc$DN.tac
+tacfr2x sim1/tacmove_idif$DN.tac sim1/steps.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/steps.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc$DN.tac sim1/aucerror$DN.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/trad_mid_idif$DN.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc$DN.tac sim1/temp2.tac
+tacadd sim1/aucerror$DN.tac sim1/temp2.tac
+tacren sim1/aucerror$DN.tac 1 Bar
+tacren sim1/aucerror$DN.tac 2 Mid_point
+tac2svg -x2=125 -legend=y -mt="" -yt="Error in AUC (s*kBq/mL)" sim1/sim_frames_and_peak_aucerror$DN.svg -l sim1/aucerror$DN.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+DN=03
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/blood$DN.tac sim1/auc$DN.tac
+tacfr2x sim1/tacmove_idif$DN.tac sim1/steps.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/steps.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc$DN.tac sim1/aucerror$DN.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/trad_mid_idif$DN.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc$DN.tac sim1/temp2.tac
+tacadd sim1/aucerror$DN.tac sim1/temp2.tac
+tacren sim1/aucerror$DN.tac 1 Bar
+tacren sim1/aucerror$DN.tac 2 Mid_point
+tac2svg -x2=125 -legend=y -mt="" -yt="Error in AUC (s*kBq/mL)" sim1/sim_frames_and_peak_aucerror$DN.svg -l sim1/aucerror$DN.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+DN=04
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/blood$DN.tac sim1/auc$DN.tac
+tacfr2x sim1/tacmove_idif$DN.tac sim1/steps.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/steps.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc$DN.tac sim1/aucerror$DN.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/trad_mid_idif$DN.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc$DN.tac sim1/temp2.tac
+tacadd sim1/aucerror$DN.tac sim1/temp2.tac
+tacren sim1/aucerror$DN.tac 1 Bar
+tacren sim1/aucerror$DN.tac 2 Mid_point
+tac2svg -x2=125 -legend=y -mt="" -yt="Error in AUC (s*kBq/mL)" sim1/sim_frames_and_peak_aucerror$DN.svg -l sim1/aucerror$DN.tac
+if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
+
+DN=05
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/blood$DN.tac sim1/auc$DN.tac
+tacfr2x sim1/tacmove_idif$DN.tac sim1/steps.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/steps.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc$DN.tac sim1/aucerror$DN.tac
+interpol -i -x=5,10,15,20,25,30,35,40,60,70,80,90,100,120 sim1/trad_mid_idif$DN.tac sim1/temp.tac
+taccalc sim1/temp.tac - sim1/auc$DN.tac sim1/temp2.tac
+tacadd sim1/aucerror$DN.tac sim1/temp2.tac
+tacren sim1/aucerror$DN.tac 1 Bar
+tacren sim1/aucerror$DN.tac 2 Mid_point
+tac2svg -x2=125 -legend=y -mt="" -yt="Error in AUC (s*kBq/mL)" sim1/sim_frames_and_peak_aucerror$DN.svg -l sim1/aucerror$DN.tac
 if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
 
 
-printf "=====================================================================\n"
-printf "Compute bias in IDIFs caused by framing but time delay is accounted for \n"
-printf "=====================================================================\n"
-
-tacformat -mid -f=PMOD idif00.tac midtime.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-
-interpol -f=midtime.tac blood.tac bl00.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-interpol -f=midtime.tac blood01.tac bl01.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-interpol -f=midtime.tac blood02.tac bl02.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-interpol -f=midtime.tac blood03.tac bl03.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-interpol -f=midtime.tac blood04.tac bl04.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-interpol -f=midtime.tac blood05.tac bl05.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-interpol -f=midtime.tac blood06.tac bl06.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-interpol -f=midtime.tac blood07.tac bl07.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-interpol -f=midtime.tac blood08.tac bl08.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-interpol -f=midtime.tac blood09.tac bl09.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-interpol -f=midtime.tac blood10.tac bl10.tac 
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-
-taccalc idif00.tac - bl00.tac dif00.tac
-taccalc idif01.tac - bl01.tac dif01.tac
-taccalc idif02.tac - bl02.tac dif02.tac
-taccalc idif03.tac - bl03.tac dif03.tac
-taccalc idif04.tac - bl04.tac dif04.tac
-taccalc idif05.tac - bl05.tac dif05.tac
-taccalc idif06.tac - bl06.tac dif06.tac
-taccalc idif07.tac - bl07.tac dif07.tac
-taccalc idif08.tac - bl08.tac dif08.tac
-taccalc idif09.tac - bl09.tac dif09.tac
-taccalc idif10.tac - bl10.tac dif10.tac
-
-taccalc dif00.tac div bl00.tac rdif00.tac
-taccalc dif01.tac div bl01.tac rdif01.tac
-taccalc dif02.tac div bl02.tac rdif02.tac
-taccalc dif03.tac div bl03.tac rdif03.tac
-taccalc dif04.tac div bl04.tac rdif04.tac
-taccalc dif05.tac div bl05.tac rdif05.tac
-taccalc dif06.tac div bl06.tac rdif06.tac
-taccalc dif07.tac div bl07.tac rdif07.tac
-taccalc dif08.tac div bl08.tac rdif08.tac
-taccalc dif09.tac div bl09.tac rdif09.tac
-taccalc dif10.tac div bl10.tac rdif10.tac
-
-taccalc rdif00.tac x 100 pdif00.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-taccalc rdif01.tac x 100 pdif01.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-taccalc rdif02.tac x 100 pdif02.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-taccalc rdif03.tac x 100 pdif03.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-taccalc rdif04.tac x 100 pdif04.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-taccalc rdif05.tac x 100 pdif05.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-taccalc rdif06.tac x 100 pdif06.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-taccalc rdif07.tac x 100 pdif07.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-taccalc rdif08.tac x 100 pdif08.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-taccalc rdif09.tac x 100 pdif09.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-taccalc rdif10.tac x 100 pdif10.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-
-printf "\nPlot percent bias curves \n"
-tac2svg -legend=n -yt="% bias" -mt="" sim_frames_and_peak_idif_bias.svg -l pdif??.tac
-if [ $? -ne 0 ] ; then printf "Failed!\n" ; exit 1 ; fi
-
 
 
 printf "\n===================================================================\n"
-- 
GitLab