Skip to content
Snippets Groups Projects
Commit 99a5d20d authored by Vesa Oikonen's avatar Vesa Oikonen
Browse files

final

parent 1f535e01
Branches
No related tags found
No related merge requests found
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
#!/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"
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment