adding scan script for SUSYHD

FlexibleSUSY · Apr 28, 2016 · 6f91c34 · 6f91c34
1 parent bbdacec
commit 6f91c34
Show file tree

Hide file tree

Showing 2 changed files with 37 additions and 2 deletions.
diff --git a/run_SUSYHD.m b/run_SUSYHD.m
@@ -0,0 +1,31 @@
+Needs["SUSYHD`"];
+scale=Q;
+tanBeta=TB;
+Xtt=Xt;
+gauginoM1=scale;
+gauginoM2=scale;
+gauginoM3=scale;
+higgsMu=scale;
+higgsAt=(higgsMu/tanBeta + Xtt scale);
+mq3=scale;
+mu3=scale;
+md3=scale;
+mq2=scale;
+mu2=scale;
+md2=scale;
+mq1=scale;
+mu1=scale;
+md1=scale;
+ml3=scale;
+me3=scale;
+ml2=scale;
+me2=scale;
+ml1=scale;
+me1=scale;
+mA=scale;
+mass = MHiggs[{tanBeta, gauginoM1, gauginoM2, gauginoM3, higgsMu, higgsAt, mq3, mu3, md3, mq2, mu2, md2, mq1, mu1, md1, ml3, me3, ml2, me2, ml1, me1, mA}, Rscale->scale, scheme->"DRbar", hiOrd->{1,1,1,1}, numerical->True, split->False];
+dmass= \[CapitalDelta]MHiggs[{tanBeta, gauginoM1, gauginoM2, gauginoM3, higgsMu, higgsAt, mq3, mu3, md3, mq2, mu2, md2, mq1, mu1, md1, ml3, me3, ml2, me2, ml1, me1, mA}, Rscale->scale, scheme->"DRbar", sources->{1,1,1}, numerical->False];
+
+Write[Streams["stderr"], ToString[mass] <> " " <> ToString[dmass]];
+
+Quit[0]
diff --git a/scan.sh b/scan.sh
@@ -359,7 +359,7 @@ if test $# -gt 0 ; then
 fi
 
 printf "# MS = ${MS}, TanBeta = ${TB}, Xt = ${Xt}\n"
-printf "# %14s %16s %16s %16s %16s %16s %16s %16s %16s\n" "$parameter" "MSSMtower" "EFTtower" "MSSMMuBMu" "HSSUSY" "Softsusy" "MSSMMuBMuSPheno" "FeynHiggs" "DeltaFeynHiggs"
+printf "# %14s %16s %16s %16s %16s %16s %16s %16s %16s %16s %16s\n" "$parameter" "MSSMtower" "EFTtower" "MSSMMuBMu" "HSSUSY" "Softsusy" "MSSMMuBMuSPheno" "FeynHiggs" "DeltaFeynHiggs" "SUSYHD" "DeltaSUSYHD"
 
 for i in `seq 0 $steps`; do
     # calculate current value for the scanned variable
@@ -415,6 +415,10 @@ EOF
     MhFH=$(echo "$FHout" | awk '{ print $1 }')
     DeltaMhFH=$(echo "$FHout" | awk '{ print $2 }')
 
-    printf "%16s %16s %16s %16s %16s %16s %16s %16s %16s\n" "$value" "$MhMSSMtower" "$MhEFTtower" "$MhMSSMMuBMu" "$MhHSSUSY" "$MhSoftsusy" "$MhMSSMMuBMuSPheno" "$MhFH" "$DeltaMhFH"
+    SUSYHDout=$(math -run "Q=${MS}; TB=${TB}; Xt=${Xt}; Get[\"run_SUSYHD.m\"]" 2>&1 >/dev/null)
+    MhSUSYHD=$(echo "$SUSYHDout" | awk '{ print $1 }')
+    DeltaMhSUSYHD=$(echo "$SUSYHDout" | awk '{ print $2 }')
+
+    printf "%16s %16s %16s %16s %16s %16s %16s %16s %16s %16s %16s\n" "$value" "$MhMSSMtower" "$MhEFTtower" "$MhMSSMMuBMu" "$MhHSSUSY" "$MhSoftsusy" "$MhMSSMMuBMuSPheno" "$MhFH" "$DeltaMhFH" "$MhSUSYHD" "$DeltaMhSUSYHD"
 
 done