enable redirection to specific output directory

FlexibleSUSY · Jul 14, 2016 · d53f42a · d53f42a
1 parent cb38392
commit d53f42a
Show file tree

Hide file tree

Showing 2 changed files with 92 additions and 79 deletions.
diff --git a/scan_MSSM.sh b/scan_MSSM.sh
@@ -2,11 +2,16 @@ start=91.1876
 stop=100000
 n_points=60
 
+DIR=${OUTPUT_DIR:-.}
+
 # for Xt in 0 2.44949 -2.44949
 for Xt in 0 2 -2
 do
 
-./scan.sh --parameter=MS --start=$start --stop=$stop --steps=$n_points --step-size=log --TB=5 --Xt="${Xt}" | tee scale_MSSM_TB-5_Xt-${Xt}.dat
+echo "Running with Xt = ${Xt} (saving in ${DIR})"
+
+./scan.sh --parameter=MS --start=$start --stop=$stop --steps=$n_points --step-size=log --TB=5 --Xt="${Xt}" \
+    | tee "${DIR}"/scale_MSSM_TB-5_Xt-${Xt}.dat
 
 slha_templ="
 Block MODSEL                 # Select model
@@ -65,7 +70,7 @@ echo "$slha_templ" | ./utils/scan-slha.sh \
     --scan-range=MS[]=91~100000:$n_points \
     --step-size=log \
     --output=MS[],MASS[25] \
-    > scale_MSSMMuBMu_TB-5_Xt-${Xt}.dat
+    > "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}.dat
 
 { echo "$slha_templ";
   cat <<EOF
@@ -77,7 +82,7 @@ EOF
     --scan-range=MS[]=91~100000:$n_points \
     --step-size=log \
     --output=MS[],MASS[25] \
-    > scale_MSSMMuBMu_TB-5_Xt-${Xt}_DeltaMt_high.dat
+    > "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}_DeltaMt_high.dat
 
 { echo "$slha_templ";
   cat <<EOF
@@ -89,12 +94,12 @@ EOF
     --scan-range=MS[]=91~100000:$n_points \
     --step-size=log \
     --output=MS[],MASS[25] \
-    > scale_MSSMMuBMu_TB-5_Xt-${Xt}_DeltaMt_low.dat
+    > "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}_DeltaMt_low.dat
 
-paste scale_MSSMMuBMu_TB-5_Xt-${Xt}.dat \
-      scale_MSSMMuBMu_TB-5_Xt-${Xt}_DeltaMt_low.dat \
-      scale_MSSMMuBMu_TB-5_Xt-${Xt}_DeltaMt_high.dat \
-      > scale_MSSMMuBMu_TB-5_Xt-${Xt}_DeltaMt.dat
+paste "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}.dat \
+      "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}_DeltaMt_low.dat \
+      "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}_DeltaMt_high.dat \
+      > "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}_DeltaMt.dat
 
 echo "uncertainty from Q in the tower"
 
@@ -104,7 +109,7 @@ echo "$slha_templ" | \
         --scan-range=MS[]=91~100000:$n_points \
         --step-size=log \
         --output=MS[],MASS[25] \
-        | tee scale_MSSMtower_TB-${TB}_Xt-${Xt}_scale_uncertainty.dat
+        | tee "${DIR}"/scale_MSSMtower_TB-${TB}_Xt-${Xt}_scale_uncertainty.dat
 
 echo "calculate Q uncertainty"
 
@@ -113,7 +118,7 @@ echo "$slha_templ" | ./utils/scan-slha.sh \
     --scan-range=MS[]=91~100000:$n_points \
     --step-size=log \
     --output=MS[],MASS[25] \
-    | tee scale_MSSMMuBMu_TB-5_Xt-${Xt}_scale_uncertainty.dat.$$
+    | tee "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}_scale_uncertainty.dat.$$
 
 echo "calculating uncertainty from Q_match in the tower"
 
@@ -123,11 +128,11 @@ echo "$slha_templ" | \
         --scan-range=MS[]=91~100000:$n_points \
         --step-size=log \
         --output=MS[],MASS[25] \
-    | tee scale_MSSMtower_TB-5_Xt-${Xt}_Qmatch_uncertainty.dat
+    | tee "${DIR}"/scale_MSSMtower_TB-5_Xt-${Xt}_Qmatch_uncertainty.dat
 
-paste scale_MSSMMuBMu_TB-5_Xt-${Xt}.dat \
-      scale_MSSMMuBMu_TB-5_Xt-${Xt}_scale_uncertainty.dat.$$ \
-      > scale_MSSMMuBMu_TB-5_Xt-${Xt}_scale_uncertainty.dat
+paste "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}.dat \
+      "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}_scale_uncertainty.dat.$$ \
+      > "${DIR}"/scale_MSSMMuBMu_TB-5_Xt-${Xt}_scale_uncertainty.dat
 
 echo "calculate uncertainty in MSSMtower from varying DeltaLambda"
 
@@ -136,7 +141,7 @@ echo "$slha_templ" | ./utils/scan-slha.sh \
     --scan-range=MS[]=91~100000:$n_points \
     --step-size=log \
     --output=Ms[],MASS[25] \
-    > scale_MSSMtower_TB-5_Xt-${Xt}.dat
+    > "${DIR}"/scale_MSSMtower_TB-5_Xt-${Xt}.dat
 
 { echo "$slha_templ";
 cat <<EOF
@@ -149,7 +154,7 @@ EOF
     --scan-range=MS[]=91~100000:$n_points \
     --step-size=log \
     --output=Ms[],MASS[25] \
-    > scale_MSSMtower_TB-5_Xt-${Xt}_DeltaLambda_low.dat
+    > "${DIR}"/scale_MSSMtower_TB-5_Xt-${Xt}_DeltaLambda_low.dat
 
 { echo "$slha_templ";
 cat <<EOF
@@ -162,17 +167,17 @@ EOF
     --scan-range=MS[]=91~100000:$n_points \
     --step-size=log \
     --output=Ms[],MASS[25] \
-    > scale_MSSMtower_TB-5_Xt-${Xt}_DeltaLambda_high.dat
+    > "${DIR}"/scale_MSSMtower_TB-5_Xt-${Xt}_DeltaLambda_high.dat
 
-paste scale_MSSMtower_TB-5_Xt-${Xt}.dat \
-      scale_MSSMtower_TB-5_Xt-${Xt}_DeltaLambda_low.dat \
-      scale_MSSMtower_TB-5_Xt-${Xt}_DeltaLambda_high.dat \
-      > scale_MSSMtower_TB-5_Xt-${Xt}_DeltaLambda.dat
+paste "${DIR}"/scale_MSSMtower_TB-5_Xt-${Xt}.dat \
+      "${DIR}"/scale_MSSMtower_TB-5_Xt-${Xt}_DeltaLambda_low.dat \
+      "${DIR}"/scale_MSSMtower_TB-5_Xt-${Xt}_DeltaLambda_high.dat \
+      > "${DIR}"/scale_MSSMtower_TB-5_Xt-${Xt}_DeltaLambda.dat
 
 plot_scale="
 set terminal pdfcairo
-set output 'scale_MSSM_TB-5_Xt-${Xt}.pdf'
-set key box bottom right width -4
+set output '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.pdf'
+set key box bottom right width -2
 set logscale x
 set grid
 
@@ -195,21 +200,21 @@ set xlabel 'M_S / TeV'
 set ylabel 'M_h / GeV'
 
 plot [0.091:] [60:140] \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):2 t 'FlexibleSUSY/MSSM-tower' w lines ls 1, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):4 t 'FlexibleSUSY/MSSM' w lines ls 3, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):5 t 'FlexibleSUSY/HSSUSY' w lines ls 2, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):6 t 'SOFTSUSY 3.6.2' w lines ls 4, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):12 t 'SPheno 3.3.8' w lines ls 6, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):14 t 'FlexibleSUSY/MSSM m_t(M_S)' w points ls 12, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):15 t 'FlexibleSUSY/MSSM m_t(M_S) SPheno-like' w points ls 13, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):8 t 'FeynHiggs 2.11.3' w lines ls 8, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):(\$8-\$9):(\$8+\$9) t 'FeynHiggs uncertainty' w filledcurves ls 9 fs transparent solid 0.3, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):10 t 'SUSYHD 1.0.2' w lines ls 10, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):(\$10-\$11):(\$10+\$11) t 'SUSYHD uncertainty' w filledcurves ls 11 fs transparent solid 0.3, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):16 t 'FlexibleSUSY/MSSM-tower (0L y_t(M_S))' w points ls 14, \
-     'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):17 t 'SuSpect 2.43' w lines ls 15
-#    'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):7 t 'FlexibleSUSY/MSSM SPheno-like' w lines ls 5, \
-#    'scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):13 t 'SPheno/MSSM FS-like' w lines ls 7, \
+     '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):2 t 'FlexibleSUSY/MSSM-tower' w lines ls 1, \
+     '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):4 t 'FlexibleSUSY/MSSM' w lines ls 3, \
+     '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):5 t 'FlexibleSUSY/HSSUSY' w lines ls 2, \
+     '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):6 t 'SOFTSUSY 3.6.2' w lines ls 4, \
+     '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):12 t 'SPheno 3.3.8' w lines ls 6, \
+     '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):8 t 'FeynHiggs 2.11.3' w lines ls 8, \
+     '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):(\$8-\$9):(\$8+\$9) t '' w filledcurves ls 9 fs transparent solid 0.3, \
+     '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):10 t 'SUSYHD 1.0.2' w lines ls 10, \
+     '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):(\$10-\$11):(\$10+\$11) t '' w filledcurves ls 11 fs transparent solid 0.3, \
+     '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):17 t 'SuSpect 2.43' w lines ls 5, \
+#    '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):14 t 'FlexibleSUSY/MSSM m_t(M_S)' w points ls 12, \
+#    '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):15 t 'FlexibleSUSY/MSSM m_t(M_S) SPheno-like' w points ls 13, \
+#    '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):16 t 'FlexibleSUSY/MSSM-tower (0L y_t(M_S))' w points ls 14, \
+#    '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):7 t 'FlexibleSUSY/MSSM SPheno-like' w lines ls 5, \
+#    '${DIR}/scale_MSSM_TB-5_Xt-${Xt}.dat' u (\$1/1000):13 t 'SPheno/MSSM FS-like' w lines ls 7, \
 "
 
 echo "$plot_scale" | gnuplot

diff --git a/scan_MSSM_xt.sh b/scan_MSSM_xt.sh
@@ -3,8 +3,13 @@ stop=3.5
 n_points=60
 TB=5
 
-for MS in 1000 2000 10000
+DIR=${OUTPUT_DIR:-.}
+
+for MS in 1000 2000 10000 30000
 do
+
+echo "Running with MS = ${MS} (saving in ${DIR})"
+
 slha_templ="
 Block FlexibleSUSY
     0   1.000000000e-05      # precision goal
@@ -80,7 +85,7 @@ echo "calculating Mh(Xt)"
           --step-size=linear \
           --TB=${TB} \
           --MS=${MS} \
-    | tee xt_TB-${TB}_MS-${MS}.dat
+    | tee "${DIR}"/xt_TB-${TB}_MS-${MS}.dat
 
 
 echo "calculate DeltaMt uncertainty in the MSSM"
@@ -94,7 +99,7 @@ EOF
         --spectrum-generator=models/MSSMMuBMu/run_MSSMMuBMu.x \
         --scan-range=Xtt[]=${start}:${stop}:${n_points} \
         --output=Xtt[],MASS[25] \
-        | tee xt_MSSMMuBMu_TB-${TB}_MS-${MS}_DeltaMt_high.dat
+        | tee "${DIR}"/xt_MSSMMuBMu_TB-${TB}_MS-${MS}_DeltaMt_high.dat
 
 { echo "$slha_templ";
   cat <<EOF
@@ -105,11 +110,11 @@ EOF
         --spectrum-generator=models/MSSMMuBMu/run_MSSMMuBMu.x \
         --scan-range=Xtt[]=${start}:${stop}:${n_points} \
         --output=Xtt[],MASS[25] \
-        | tee xt_MSSMMuBMu_TB-${TB}_MS-${MS}_DeltaMt_low.dat
+        | tee "${DIR}"/xt_MSSMMuBMu_TB-${TB}_MS-${MS}_DeltaMt_low.dat
 
-paste xt_MSSMMuBMu_TB-${TB}_MS-${MS}_DeltaMt_low.dat \
-      xt_MSSMMuBMu_TB-${TB}_MS-${MS}_DeltaMt_high.dat \
-      > xt_MSSMMuBMu_TB-${TB}_MS-${MS}_DeltaMt.dat
+paste "${DIR}"/xt_MSSMMuBMu_TB-${TB}_MS-${MS}_DeltaMt_low.dat \
+      "${DIR}"/xt_MSSMMuBMu_TB-${TB}_MS-${MS}_DeltaMt_high.dat \
+      > "${DIR}"/xt_MSSMMuBMu_TB-${TB}_MS-${MS}_DeltaMt.dat
 
 echo "calculate Q uncertainty in the MSSM"
 
@@ -118,7 +123,7 @@ echo "$slha_templ" | \
         --spectrum-generator=./MSSMMuBMu_uncertainty.sh \
         --scan-range=Xtt[]=${start}:${stop}:${n_points} \
         --output=Xtt[],MASS[25] \
-        | tee xt_MSSMMuBMu_TB-${TB}_MS-${MS}_scale_uncertainty.dat
+        | tee "${DIR}"/xt_MSSMMuBMu_TB-${TB}_MS-${MS}_scale_uncertainty.dat
 
 echo "calculating parametric uncertainty from Q in the tower"
 
@@ -127,7 +132,7 @@ echo "$slha_templ" | \
         --spectrum-generator=./MSSMtower_uncertainty.sh \
         --scan-range=Xtt[]=${start}:${stop}:${n_points} \
         --output=Xtt[],MASS[25] \
-        | tee xt_MSSMtower_TB-${TB}_MS-${MS}_scale_uncertainty.dat
+        | tee "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_scale_uncertainty.dat
 
 echo "calculating parametric uncertainty from Q_match in the tower"
 
@@ -136,7 +141,7 @@ echo "$slha_templ" | \
         --spectrum-generator=./MSSMtower_Qmatch_uncertainty.sh \
         --scan-range=Xtt[]=${start}:${stop}:${n_points} \
         --output=Xtt[],MASS[25] \
-        | tee xt_MSSMtower_TB-${TB}_MS-${MS}_Qmatch_uncertainty.dat
+        | tee "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_Qmatch_uncertainty.dat
 
 echo "calculating parametric uncertainty from delta in the tower"
 
@@ -145,46 +150,47 @@ echo "$slha_templ_delta_low" | \
         --spectrum-generator=models/MSSMtower/run_MSSMtower.x \
         --scan-range=Xtt[]=${start}:${stop}:${n_points} \
         --output=Xtt[],MASS[25] \
-        | tee xt_MSSMtower_TB-${TB}_MS-${MS}_delta_low.dat
+        | tee "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_delta_low.dat
 
 echo "$slha_templ_delta_high" | \
     ./utils/scan-slha.sh \
         --spectrum-generator=models/MSSMtower/run_MSSMtower.x \
         --scan-range=Xtt[]=${start}:${stop}:${n_points} \
         --output=Xtt[],MASS[25] \
-        | tee xt_MSSMtower_TB-${TB}_MS-${MS}_delta_high.dat
+        | tee "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_delta_high.dat
 
 printf "# %14s %16s\n" \
        "MS" "Mh" \
-       > xt_MSSMtower_TB-${TB}_MS-${MS}_Mh.dat
+       > "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_Mh.dat
 
-awk '{ if ($1 !~ "#") { printf "%16s %16s\n", $1, $2 } }' xt_TB-${TB}_MS-${MS}.dat \
-    >> xt_MSSMtower_TB-${TB}_MS-${MS}_Mh.dat
+awk '{ if ($1 !~ "#") { printf "%16s %16s\n", $1, $2 } }' \
+    "${DIR}"/xt_TB-${TB}_MS-${MS}.dat \
+    >> "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_Mh.dat
 
-paste xt_MSSMtower_TB-${TB}_MS-${MS}_Mh.dat \
-      xt_MSSMtower_TB-${TB}_MS-${MS}_scale_uncertainty.dat \
-      xt_MSSMtower_TB-${TB}_MS-${MS}_Qmatch_uncertainty.dat \
-      xt_MSSMtower_TB-${TB}_MS-${MS}_delta_low.dat \
-      xt_MSSMtower_TB-${TB}_MS-${MS}_delta_high.dat \
-      > xt_MSSMtower_TB-${TB}_MS-${MS}.dat.$$
+paste "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_Mh.dat \
+      "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_scale_uncertainty.dat \
+      "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_Qmatch_uncertainty.dat \
+      "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_delta_low.dat \
+      "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_delta_high.dat \
+      > "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}.dat.$$
 
 printf "# %14s %16s %16s %16s %16s %16s %16s %16s %16s %16s\n" \
        "MS" "Mh" \
        "MS" "Q uncert." \
        "MS" "Qmatch uncert." \
        "MS" "lambda(2L) low" \
        "MS" "lambda(2L) high" \
-       > xt_MSSMtower_TB-${TB}_MS-${MS}_uncertainties.dat
+       > "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_uncertainties.dat
 
-cat xt_MSSMtower_TB-${TB}_MS-${MS}.dat.$$ >> xt_MSSMtower_TB-${TB}_MS-${MS}_uncertainties.dat
+cat "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}.dat.$$ \
+    >> "${DIR}"/xt_MSSMtower_TB-${TB}_MS-${MS}_uncertainties.dat
 
     plot_scale="
-set terminal pdfcairo enhanced size 5in,4in
-set output 'xt_TB-${TB}_MS-${MS}.pdf'
-set key box top center width -2 at graph 0.5, graph 1.4 opaque
+set terminal pdfcairo enhanced size 5in,3.4in
+set output '${DIR}/xt_TB-${TB}_MS-${MS}.pdf'
+set key box bottom right width -2
 set border back
 set grid
-set tmargin 8
 
 set style line 1 lt 1 dt 1 lw 2 lc rgb '#FF0000'
 set style line 2 lt 1 dt 2 lw 2 lc rgb '#0000FF'
@@ -207,25 +213,27 @@ set ylabel 'M_h / GeV'
 min(x,y) = x < y ? x : y
 max(x,y) = x < y ? y : x
 
-filename = 'xt_TB-${TB}_MS-${MS}.dat'
-filename2 = 'xt_MSSMtower_TB-${TB}_MS-${MS}_uncertainties.dat'
+filename = '${DIR}/xt_TB-${TB}_MS-${MS}.dat'
+filename2 = '${DIR}/xt_MSSMtower_TB-${TB}_MS-${MS}_uncertainties.dat'
 
-plot [:] [:] \
-     filename u 1:2 t 'FS/MSSM-tower' w lines ls 1, \
-     filename u 1:4 t 'FS/MSSM' w lines ls 3, \
-     filename u 1:5 t 'FS/HSSUSY' w lines ls 2, \
+plot [-3.5:3.5] [:] \
+     filename u 1:2 t 'FlexibleSUSY/MSSM-tower' w lines ls 1, \
+     filename u 1:4 t 'FlexibleSUSY/MSSM' w lines ls 3, \
+     filename u 1:5 t 'FlexibleSUSY/HSSUSY' w lines ls 2, \
      filename u 1:6 t 'SOFTSUSY 3.6.2' w lines ls 4, \
-     filename u 1:7 t 'FS/MSSM SPheno-like' w lines ls 5, \
      filename u 1:8 t 'FeynHiggs 2.11.3' w lines ls 8, \
+     filename u 1:(\$8-\$9):(\$8+\$9) t '' w filledcurves ls 9 fs transparent solid 0.3, \
      filename u 1:10 t 'SUSYHD 1.0.2' w lines ls 10, \
+     filename u 1:(\$10-\$11):(\$10+\$11) t '' w filledcurves ls 11 fs transparent solid 0.3, \
      filename u 1:12 t 'SPheno 3.3.8' w lines ls 6, \
-     filename u 1:14 t 'FS/MSSM m_t(M_S)' w points ls 12, \
-     filename u 1:15 t 'FS/MSSM m_t(M_S) SPheno-like' w points ls 13, \
-     filename u 1:16 t 'FS/MSSM-tower (0L y_t(M_S))' w lines ls 14, \
-     filename2 u 1:(\$2-\$4/2):(\$2+\$4/2) t '{/Symbol D} Q' w filledcurves ls 7 fs transparent solid 0.3, \
-     filename2 u 1:(\$2-\$6/2):(\$2+\$6/2) t '{/Symbol D} Q_{match}' w filledcurves ls 5 fs transparent solid 0.3, \
-     filename2 u 1:(min(\$8,\$10)):(max(\$8,\$10)) t '{/Symbol D} {/Symbol l}^{(2)}' w filledcurves ls 10 fs transparent solid 0.3
-
+     filename u 1:17 t 'SuSpect 2.43' w lines ls 5, \
+#    filename u 1:7 t 'FS/MSSM SPheno-like' w lines ls 5, \
+#    filename u 1:14 t 'FS/MSSM m_t(M_S)' w points ls 12, \
+#    filename u 1:15 t 'FS/MSSM m_t(M_S) SPheno-like' w points ls 13, \
+#    filename u 1:16 t 'FS/MSSM-tower (0L y_t(M_S))' w lines ls 14, \
+#    filename2 u 1:(\$2-\$4/2):(\$2+\$4/2) t '{/Symbol D} Q' w filledcurves ls 7 fs transparent solid 0.3, \
+#    filename2 u 1:(\$2-\$6/2):(\$2+\$6/2) t '{/Symbol D} Q_{match}' w filledcurves ls 5 fs transparent solid 0.3, \
+#    filename2 u 1:(min(\$8,\$10)):(max(\$8,\$10)) t '{/Symbol D} {/Symbol l}^{(2)}' w filledcurves ls 10 fs transparent solid 0.3
 "
 
 echo "$plot_scale" | gnuplot