|
Go to the documentation of this file. 1Get["models/HSSUSY/HSSUSY_librarylink.m"];
2Get["model_files/HSSUSY/HSSUSY_uncertainty_estimate.m"];
5 precisionGoal -> 1.*^-5,
6 poleMassLoopOrder -> 4,
8 betaFunctionLoopOrder -> 5,
9 thresholdCorrectionsLoopOrder -> 3,
10 thresholdCorrections -> 124111421
14 alphaEmMZ -> 1/127.916, (* SMINPUTS[1] *)
15 GF -> 1.166378700*^-5, (* SMINPUTS[2] *)
16 alphaSMZ -> 0.1184, (* SMINPUTS[3] *)
17 MZ -> 91.1876, (* SMINPUTS[4] *)
18 mbmb -> 4.18, (* SMINPUTS[5] *)
19 Mt -> 173.34, (* SMINPUTS[6] *)
20 Mtau -> 1.77699, (* SMINPUTS[7] *)
21 Mv3 -> 0, (* SMINPUTS[8] *)
22 MW -> 80.385, (* SMINPUTS[9] *)
23 Me -> 0.000510998902, (* SMINPUTS[11] *)
24 Mv1 -> 0, (* SMINPUTS[12] *)
25 Mm -> 0.1056583715, (* SMINPUTS[13] *)
26 Mv2 -> 0, (* SMINPUTS[14] *)
27 md2GeV -> 0.00475, (* SMINPUTS[21] *)
28 mu2GeV -> 0.0024, (* SMINPUTS[22] *)
29 ms2GeV -> 0.104, (* SMINPUTS[23] *)
30 mcmc -> 1.27, (* SMINPUTS[24] *)
41 alphaEm0 -> 1/137.035999074,
45HSSUSYCalcMh[MS_, TB_, Xtt_] :=
47 fsSettings -> settings,
48 fsSMParameters -> smpars,
49 fsModelParameters -> {
58 AtInput -> (Xtt + 1/TB) MS,
61 msq2 -> MS^2 IdentityMatrix[3],
62 msu2 -> MS^2 IdentityMatrix[3],
63 msd2 -> MS^2 IdentityMatrix[3],
64 msl2 -> MS^2 IdentityMatrix[3],
65 mse2 -> MS^2 IdentityMatrix[3],
76LinearRange[start_, stop_, steps_] :=
77 Range[start, stop, (stop - start)/steps];
83 { N[#], Sequence @@ HSSUSYCalcMh[#, TB, Xtt] }&,
84 LinearRange[500, 10^4, 100]
87MhMin[{MS_, Mh_, DMh_}] := {MS, Mh - DMh};
88MhMax[{MS_, Mh_, DMh_}] := {MS, Mh + DMh};
89MhBest[{MS_, Mh_, DMh_}] := {MS, Mh};
96 PlotStyle -> {Red, Thick}];
98plot1 = ListLinePlot[{dataMhMax, dataMhMin},
99 PlotStyle -> LightGray,
100 Filling -> {1 -> {{2}, LightGray}},
104 BaseStyle -> {FontSize -> 16, FontFamily -> "Helvetica"},
105 PlotLabel -> Style["\*SubscriptBox[X, t] = 2.44949 \*SubscriptBox[M, S], tan\[Beta] = 20"],
106 PlotRange -> Automatic,
107 Axes -> False, Frame -> True,
108 FrameLabel -> {Style["\*SubscriptBox[M, S] / GeV"],
109 Style["\*SubscriptBox[M, h] / GeV"]}];
111Export["HSSUSY_Mh_MS.png", plot, ImageSize -> 600];
Get["models/HSSUSY/HSSUSY_librarylink.m"]
Export["HSSUSY_Mh_MS.png", plot, ImageSize -> 600]
MhBest[{MS_, Mh_, DMh_}] Mh
|