forked from Expander/FlexibleSUSY
/
scan_MRSSM.sh
287 lines (256 loc) · 12.4 KB
/
scan_MRSSM.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
n_points=60
parameter_point="
Block MINPAR
3 5 # TanBeta
Block HMIXIN
301 -0.01 # LSD
302 -0.01 # LSU
303 -0.01 # LTD
304 -0.01 # LTU
201 1e3 # MuD
202 1e3 # MuU
"
slha_templ="
Block MODSEL
1 1 # 1/0: High/low scale input
2 1 # Boundary Condition
6 1 # Generation Mixing
Block FlexibleSUSY
0 1.000000000e-04 # precision goal
1 0 # max. iterations (0 = automatic)
2 0 # algorithm (0 = two_scale, 1 = lattice)
3 0 # calculate SM pole masses
4 2 # pole mass loop order
5 2 # EWSB loop order
6 2 # beta-functions loop order
7 2 # 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 # Higgs logarithmic resummation
15 1.000000000e-11 # beta-function zero threshold
16 0 # calculate observables (a_muon, ...)
17 0 # Mt method (0 = FS, 1 = SPheno)
18 0 # write EFT output
Block SMINPUTS # Standard Model inputs
1 1.279440000e+02 # alpha^(-1) SM MSbar(MZ)
2 1.166380000e-05 # G_Fermi
3 1.184000000e-01 # alpha_s(MZ) SM MSbar
4 9.118760000e+01 # MZ(pole)
5 4.180000000e+00 # mb(mb) SM MSbar
6 1.733400000e+02 # mtop(pole)
7 1.777000000e+00 # mtau(pole)
8 0.000000000e+00 # mnu3(pole)
9 80.384 # 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
${parameter_point}
"
slha_templ_fs="
${slha_templ}
Block EXTPAR
0 1e4 # Ms
"
slha_templ_spheno_like="
${slha_templ}
Block FlexibleSUSY
17 1 # Mt method (0 = FS, 1 = SPheno)
"
slha_templ_spheno_1L="
${slha_templ}
Block SPhenoInput # SPheno specific input
1 -1 # error level
2 0 # SPA conventions
7 1 # Skip 2-loop Higgs corrections
8 3 # Method used for two-loop calculation
9 1 # Gaugeless limit used at two-loop
10 1 # safe-mode used at two-loop
11 0 # calculate branching ratios
13 0 # 3-Body decays: none (0), fermion (1), scalar (2), both (3)
14 0 # Run couplings to scale of decaying particle
12 1.000E-04 # write only branching ratios larger than this value
15 1.000E-30 # write only decay if width larger than this value
31 -1 # fixed GUT scale (-1: dynamical GUT scale)
32 0 # Strict unification
34 1.000E-04 # Precision of mass calculation
35 40 # Maximal number of iterations
36 5 # Minimal number of iterations before discarding points
37 1 # Set Yukawa scheme
38 2 # 1- or 2-Loop RGEs
50 1 # Majorana phases: use only positive masses (put 0 to use file with CalcHep/Micromegas!)
51 0 # Write Output in CKM basis
52 0 # Write spectrum in case of tachyonic states
55 1 # Calculate loop corrected masses
57 0 # Calculate low energy constraints
65 1 # Solution tadpole equation
75 1 # Write WHIZARD files
76 1 # Write HiggsBounds file
86 0 # Maximal width to be counted as invisible in Higgs decays; -1: only LSP
510 0 # Write tree level values for tadpole solutions
515 0 # Write parameter values at GUT scale
520 0 # Write effective Higgs couplings (HiggsBounds blocks): put 0 to use file with MadGraph!
521 0 # Diphoton/Digluon widths including higher order
525 0 # Write loop contributions to diphoton decay of Higgs
530 0 # Write Blocks for Vevacious
"
slha_templ_spheno_2L="
${slha_templ_spheno_1L}
Block SPhenoInput # SPheno specific input
7 0 # Skip 2-loop Higgs corrections
"
echo "$slha_templ" | ./utils/scan-slha.sh \
--spectrum-generator=models/MRSSMtower/run_MRSSMtower.x \
--scan-range=EXTPAR[0]=91~100000:$n_points \
--step-size=log \
--output=EXTPAR[0],MASS[25] \
> scale_MRSSMtower_TB-5.dat
echo "$slha_templ" | ./utils/scan-slha.sh \
--spectrum-generator=models/MRSSMMSUSY/run_MRSSMMSUSY.x \
--scan-range=EXTPAR[0]=91~100000:$n_points \
--step-size=log \
--output=EXTPAR[0],MASS[25] \
> scale_MRSSMMSUSY_TB-5.dat
echo "$slha_templ_spheno_like" | ./utils/scan-slha.sh \
--spectrum-generator=models/MRSSMMSUSY/run_MRSSMMSUSY.x \
--scan-range=EXTPAR[0]=91~100000:$n_points \
--step-size=log \
--output=EXTPAR[0],MASS[25] \
> scale_MRSSMMSUSY_TB-5_SPheno-like.dat
echo "$slha_templ" | ./utils/scan-slha.sh \
--spectrum-generator=models/MRSSMMSUSYYuatMS/run_MRSSMMSUSYYuatMS.x \
--scan-range=EXTPAR[0]=91~100000:$n_points \
--step-size=log \
--output=EXTPAR[0],MASS[25] \
> scale_MRSSMMSUSYYuatMS_TB-5.dat
echo "$slha_templ_spheno_like" | ./utils/scan-slha.sh \
--spectrum-generator=models/MRSSMMSUSYYuatMS/run_MRSSMMSUSYYuatMS.x \
--scan-range=EXTPAR[0]=91~100000:$n_points \
--step-size=log \
--output=EXTPAR[0],MASS[25] \
> scale_MRSSMMSUSYYuatMS_TB-5_SPheno-like.dat
echo "$slha_templ_spheno_1L" | ./utils/scan-slha.sh \
--spectrum-generator=./SPhenoMRSSM2 \
--scan-range=MINPAR[1]=91~100000:$n_points \
--step-size=log \
--output=MINPAR[1],MASS[25] \
--type=SPheno \
> scale_SPhenoMRSSM_TB-5_1L.dat
echo "$slha_templ_spheno_2L" | ./utils/scan-slha.sh \
--spectrum-generator=./SPhenoMRSSM2 \
--scan-range=MINPAR[1]=91~100000:$n_points \
--step-size=log \
--output=MINPAR[1],MASS[25] \
--type=SPheno \
> scale_SPhenoMRSSM_TB-5_2L.dat
echo "$slha_templ_spheno_1L" | ./utils/scan-slha.sh \
--spectrum-generator=./SPhenoMRSSM2_FlexibleSUSY_like \
--scan-range=MINPAR[1]=91~100000:$n_points \
--step-size=log \
--output=MINPAR[1],MASS[25] \
--type=SPheno \
> scale_SPhenoMRSSM_TB-5_1L_FSlike.dat
echo "$slha_templ_spheno_2L" | ./utils/scan-slha.sh \
--spectrum-generator=./SPhenoMRSSM2_FlexibleSUSY_like \
--scan-range=MINPAR[1]=91~100000:$n_points \
--step-size=log \
--output=MINPAR[1],MASS[25] \
--type=SPheno \
> scale_SPhenoMRSSM_TB-5_2L_FSlike.dat
echo "calculating parametric uncertainty from Q in FlexibleSUSY/MRSSM"
echo "$slha_templ" | ./utils/scan-slha.sh \
--spectrum-generator=./MRSSMMSUSY_uncertainty.sh \
--scan-range=EXTPAR[0]=91~100000:$n_points \
--step-size=log \
--output=EXTPAR[0],MASS[25] \
> scale_MRSSMMSUSY_TB-5_scale_uncertainty.dat
# calculate parametric uncertainty from Q
echo "calculating parametric uncertainty from Q in SPheno"
echo "$slha_templ_spheno_2L" | ./utils/scan-slha.sh \
--spectrum-generator=./SPhenoMRSSM2_uncertainty.sh \
--scan-range=MINPAR[1]=91~100000:$n_points \
--step-size=log \
--output=MINPAR[1],MASS[25] \
--type=SPheno \
> scale_SPhenoMRSSM_TB-5_2L_scale_uncertainty.dat
echo "calculating parametric uncertainty from Q in the tower"
echo "$slha_templ" | ./utils/scan-slha.sh \
--spectrum-generator=./MRSSMtower_uncertainty.sh \
--scan-range=EXTPAR[0]=91~100000:$n_points \
--step-size=log \
--output=EXTPAR[0],MASS[25] \
> scale_MRSSMtower_TB-5_scale_uncertainty.dat
echo "calculating parametric uncertainty from Q_match in the tower"
echo "$slha_templ" | ./utils/scan-slha.sh \
--spectrum-generator=./MRSSMtower_Qmatch_uncertainty.sh \
--scan-range=EXTPAR[0]=91~100000:$n_points \
--step-size=log \
--output=EXTPAR[0],MASS[25] \
> scale_MRSSMtower_TB-5_Qmatch_uncertainty.dat
echo "run MRSSM-tower with yt(0L)"
{ echo "$slha_templ";
cat <<EOF
Block FlexibleSUSY
19 1 # mf tree-level matching
EOF
} | ./utils/scan-slha.sh \
--spectrum-generator=models/MRSSMtower/run_MRSSMtower.x \
--scan-range=EXTPAR[0]=91~100000:$n_points \
--step-size=log \
--output=EXTPAR[0],MASS[25] \
> scale_MRSSMtower_TB-5_yt-0L.dat
plot_scale="
set terminal pdfcairo size 5in,4in
set tmargin 7
set border back
set output 'scale_MRSSM_uncertainty.pdf'
#set key box top left width -2
set key box top left width -2 at graph 0.01, graph 1.35 opaque
set logscale x
set grid
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'
set style line 3 lt 1 dt 4 lw 2 lc rgb '#45AD53'
set style line 4 lt 1 dt 3 lw 2 lc rgb '#FFBF00'
set style line 5 lt 1 dt 5 lw 2 lc rgb '#FF00FF'
set style line 6 lt 1 dt 6 lw 2 lc rgb '#00FFFF'
set style line 7 lt 1 dt 7 lw 2 lc rgb '#000000'
set style line 8 lt 1 dt 4 lw 2 lc rgb '#00FF00'
set style line 9 lt 1 dt 1 lw 0 lc rgb '#00FF00'
set style line 10 lt 1 dt 2 lw 2 lc rgb '#9C4C17'
set style line 12 lt 1 dt 8 lw 1 lc rgb '#000000' pt 1
set style line 13 lt 1 dt 9 lw 1 lc rgb '#000000' pt 2
set xlabel 'M_S / TeV'
set ylabel 'M_h / GeV'
min(x,y) = x < y ? x : y
max(x,y) = x < y ? y : x
plot [0.1:] [:] \
'scale_MRSSMtower_TB-5.dat' u (\$1/1000):2 t 'FS/MRSSM-tower' w lines ls 1, \
'scale_MRSSMMSUSY_TB-5.dat' u (\$1/1000):2 t 'FS/MRSSM 1L' w lines ls 3, \
'scale_MRSSMMSUSY_TB-5_SPheno-like.dat' u (\$1/1000):2 t 'FS/MRSSM 1L SPheno-like' w lines ls 5, \
'scale_MRSSMMSUSYYuatMS_TB-5.dat' u (\$1/1000):2 t 'FS/MRSSM m_t(M_S)' w points ls 12, \
'scale_MRSSMMSUSYYuatMS_TB-5_SPheno-like.dat' u (\$1/1000):2 t 'FS/MRSSM m_t(M_S) SPheno-like' w points ls 13, \
'scale_SPhenoMRSSM_TB-5_1L.dat' u (\$1/1000):2 t 'SPheno/MRSSM 1L' w lines ls 2, \
'scale_SPhenoMRSSM_TB-5_2L.dat' u (\$1/1000):2 t 'SPheno/MRSSM 2L' w lines ls 4, \
'scale_SPhenoMRSSM_TB-5_1L_FSlike.dat' u (\$1/1000):2 t 'SPheno/MRSSM 1L FS-like' w lines ls 6, \
'scale_SPhenoMRSSM_TB-5_2L_FSlike.dat' u (\$1/1000):2 t 'SPheno/MRSSM 2L FS-like' w lines ls 7, \
'scale_SPhenoMRSSM_TB-5_2L.dat' u (\$1/1000):(min(\$4,\$6)):(max(\$4,\$6)) t 'alpha_s uncertainty' w filledcurves ls 4 dt 1 lw 0 fs transparent solid 0.3, \
'scale_SPhenoMRSSM_TB-5_2L.dat' u (\$1/1000):(min(\$8,\$10)):(max(\$8,\$10)) t 'M_t uncertainty' w filledcurves ls 5 dt 1 lw 0 fs transparent solid 0.3, \
'scale_SPhenoMRSSM_TB-5_2L.dat' u (\$1/1000):(\$2-\$12/2):(\$2+\$12/2) t 'scale uncertainty' w filledcurves ls 6 dt 1 lw 0 fs transparent solid 0.3, \
'scale_MRSSMtower_TB-5.dat' u (\$1/1000):(min(\$4,\$6)):(max(\$4,\$6)) t '' w filledcurves ls 4 dt 1 lw 0 fs transparent solid 0.3, \
'scale_MRSSMtower_TB-5.dat' u (\$1/1000):(min(\$8,\$10)):(max(\$8,\$10)) t '' w filledcurves ls 5 dt 1 lw 0 fs transparent solid 0.3, \
'scale_MRSSMtower_TB-5.dat' u (\$1/1000):(min(\$12,\$14)):(max(\$12,\$14)) t '{/Symbol D}{/Symbol l}^{(2)} uncertainty' w filledcurves ls 1 dt 1 lw 0 fs transparent solid 0.3, \
'scale_MRSSMtower_TB-5.dat' u (\$1/1000):(\$2-\$16/2):(\$2+\$16/2) t '' w filledcurves ls 6 dt 1 lw 0 fs transparent solid 0.3, \
'scale_MRSSMtower_TB-5.dat' u (\$1/1000):(\$2-\$18/2):(\$2+\$18/2) t 'Q_{match} uncertainty' w filledcurves ls 2 dt 1 lw 0 fs transparent solid 0.3, \
'scale_MRSSMMSUSY_TB-5.dat' u (\$1/1000):(min(\$4,\$6)):(max(\$4,\$6)) t '{/Symbol D}M_t uncertainty' w filledcurves ls 3 dt 1 lw 0 fs transparent solid 0.3, \
'scale_MRSSMMSUSY_TB-5_SPheno-like.dat' u (\$1/1000):(min(\$4,\$6)):(max(\$4,\$6)) t '' w filledcurves ls 3 dt 1 lw 0 fs transparent solid 0.3
"
echo "$plot_scale" | gnuplot