#!/bin/sh output="MASS-25" # specify output observable parameter=MS # specify scan parameter start=100 # range stop=100000 steps=30 step_size=log # log or linear Xt=0 TB=20 MS=10000 At=$(echo "scale=10; (1./${TB} + ${Xt}) * ${MS}" | bc) M3factor=1 # Mg = M3factor * MS M3=1 AS="1.18000000e-01" AI="1.28944742e+02" MT="1.733400000e+02" # Mt_pole MB="4.20000000e-00" # mb(mb) SM MSbar MTau="1.777030000e-00" # Mtau_pole MTmethod=0 GF=0.0000116639 MZ=91.1876 loop=0 # controls SE tadpole looporder yuk=0 # controls the matching(BSM-SM) loop order of Yukawas tresh=0 lam=0 # controls the matching(BSM-SM) loop order of lambda mtpol=0 # controls yukawa corrections order at Q=MZ BSM=0 BL=3 QS=0 #dump_fs_slha_input_file= #dump_fs_slha_output_file= #dump_ss_slha_input_file= #dump_ss_slha_output_file= sminputs_tmpl="\ Block SMINPUTS # Standard Model inputs 1 128.944742 #1.279440000e+02 # alpha^(-1) SM MSbar(MZ) 2 1.166390000e-05 # G_Fermi 3 1.18000000e-01 # alpha_s(MZ) SM MSbar 4 9.118760000e+01 # MZ(pole) 5 4.20000000e+00 # mb(mb) SM MSbar 6 1.733400000e+02 # mtop(pole) 7 1.777030000e+00 # mtau(pole) 8 0.000000000e+00 # mnu3(pole) 9 80.385 # MW pole 11 5.109989020e-04 # melectron(pole) 12 0.000000000e+00 # mnu1(pole) 13 1.056583570e-01 # mmuon(pole) 14 0.000000000e+00 # mnu2(pole) 21 4.750000000e-03 # md(2 GeV) MS-bar 22 2.400000000e-03 # mu(2 GeV) MS-bar 23 1.040000000e-01 # ms(2 GeV) MS-bar 24 1.270000000e+00 # mc(mc) MS-bar " slha_tmpl="\ Block MODSEL # Select model # 6 0 # flavour violation Block FlexibleSUSY 0 1.000000000e-05 # precision goal 1 0 # max. iterations (0 = automatic) 2 0 # algorithm (0 = two_scale, 1 = lattice) 3 0 # calculate SM pole masses 4 ${loop} # pole mass loop order 5 ${loop} # EWSB loop order 6 ${BL} # beta-functions loop order 7 ${tresh} # threshold corrections loop order 8 1 # Higgs 2-loop corrections O(alpha_t alpha_s) 9 1 # Higgs 2-loop corrections O(alpha_b alpha_s) 10 1 # Higgs 2-loop corrections O((alpha_t + alpha_b)^2) 11 1 # Higgs 2-loop corrections O(alpha_tau^2) 12 0 # force output 13 1 # Top quark 2-loop corrections QCD 14 1.000000000e-11 # beta-function zero threshold 15 0 # calculate observables (a_muon, ...) 16 0 # force positive majorana masses 17 0 # pole mass renormalization scale (0 = SUSY scale) 18 0 # pole mass renormalization scale in the EFT (0 = min(SUSY scale, Mt)) 19 0 # EFT matching scale (0 = SUSY scale) 20 ${yuk} # EFT loop order for upwards matching 21 ${lam} # EFT loop order for downwards matching 22 0 # EFT index of SM-like Higgs in the BSM model 23 ${BSM} # calculate BSM pole masses 24 122111121 # individual threshold correction loop orders 25 0 # ren. scheme for Higgs 3L corrections (0 = DR, 1 = MDR) 26 0 # Higgs 3-loop corrections O(alpha_t alpha_s^2) 27 0 # Higgs 3-loop corrections O(alpha_b alpha_s^2) 28 0 # Higgs 3-loop corrections O(alpha_t^2 alpha_s) 29 0 # Higgs 3-loop corrections O(alpha_t^3) 30 0 # Higgs 4-loop corrections O(alpha_t alpha_s^3) ${sminputs_tmpl} Block MINPAR # Input parameters 4 1 # SignMu #Block EXTPAR # Input parameters # 100 2 # LambdaLoopOrder (HSSUSY) Block Ms ${MS} # SUSY scale Block TanBeta ${TB} # tan(Beta) at the SUSY scale Block Xtt ${Xt} # Xt / Ms " # prints SLHA block print_slha_block_awk=' BEGIN { is_block = 0; if (block == "") { print "Error: block name not defined"; print " Please define the block name with -v block="; exit 1 } } { pattern = "^block[[:blank:]]*" tolower(block) "([^[:graph:]].*)?$"; not_pattern = "^block[[:blank:]]*.*$"; if (tolower($0) ~ pattern) { is_block = 1 } else if (tolower($0) ~ not_pattern) { is_block = 0 }; if (is_block) print $0 } ' # prints block entry # expects block entry keys in the form x or x:y or x:y:z etc. print_block_entry_awk=' { len = split(keys,k,":"); matches = 1; for (i in k) { if ($(i) != k[i]) matches = 0 } if (matches == 1) print $(len + 1) } ' run_HSSUSY() { local SG="$1" local MS2=$(echo "scale=5; ${MS}^2" | bc) local At=$(echo "scale=10; (1./${TB} + ${Xt}) * ${MS}" | bc) local slha_output= local block= local value= local slha_input= local output_block=$(echo "${output}" | cut -d'-' -f1) local output_entry=$(echo "${output}" | cut -d'-' -f2) slha_input=$( { echo "$slha_tmpl" ; \ cat </dev/null) block=$(echo "$slha_output" | awk -v block="$output_block" "$print_slha_block_awk") value=$(echo "$block" | awk -v keys="$output_entry" "$print_block_entry_awk") [ "x$value" = "x" ] && value="-" [ "x$dump_fs_slha_input_file" != "x" ] && \ echo "$slha_input" > "$dump_fs_slha_input_file" [ "x$dump_fs_slha_output_file" != "x" ] && \ echo "$slha_output" > "$dump_fs_slha_output_file" echo $value } run_sg() { local SG="$1" local MS2=$(echo "scale=5; ${MS}^2" | bc) local At=$(echo "scale=10; (1./${TB} + ${Xt}) * ${MS}" | bc) local slha_output= local block= local value= local slha_input= local output_block=$(echo "${output}" | cut -d'-' -f1) local output_entry=$(echo "${output}" | cut -d'-' -f2) slha_input=$( { echo "$slha_tmpl" ; \ cat </dev/null) block=$(echo "$slha_output" | awk -v block="$output_block" "$print_slha_block_awk") value=$(echo "$block" | awk -v keys="$output_entry" "$print_block_entry_awk") [ "x$value" = "x" ] && value="-" [ "x$dump_fs_slha_input_file" != "x" ] && \ echo "$slha_input" > "$dump_fs_slha_input_file" [ "x$dump_fs_slha_output_file" != "x" ] && \ echo "$slha_output" > "$dump_fs_slha_output_file" echo $value } help() { cat < double conversion run with bash instead sh MStwo=$(echo "scale=10; (${MS} * 2.)" | bc ) QLhalf=$(echo "scale=10; ( 173.34 * 0.5)" | bc ) # -> double conversion run with bash instead sh QLtwo=$(echo "scale=10; ( 173.34* 2.)" | bc ) QS=0 #${MStwo} QL=0 #${QLtwo} #specification of threshold corrections thresh=3 mtpol=3 yuk=2 BSM=0 BL=4 #run the spectrum generators lam=1 loop=1 feft=$(run_HSSUSY "models/MSSMEFTHiggs/run_MSSMEFTHiggs.x") loop=3 yuk=2 lam=2 lam=3 eft=$(run_HSSUSY "models/HSSUSY/run_HSSUSY.x") mtpol=2 yuk=2 loop=3 BL=3 fo=$(run_sg "models/NUHMSSMNoFVHimalaya/run_NUHMSSMNoFVHimalaya.x") printf "%16s %16s %16s %16s\n" "$value" "$feft" "$eft" "$fo" done