|
Go to the documentation of this file. 1Get["models/MSSMEFTHiggs/MSSMEFTHiggs_librarylink.m"];
2Get["model_files/MSSMEFTHiggs/MSSMEFTHiggs_uncertainty_estimate.m"];
5 precisionGoal -> 1.*^-5,
6 poleMassLoopOrder -> 3,
8 betaFunctionLoopOrder -> 3,
9 thresholdCorrectionsLoopOrder -> 3,
10 thresholdCorrections -> 123111321
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,
45MSSMEFTHiggsCalcMh[MS_, TB_, Xtt_] :=
47 fsSettings -> settings,
48 fsSMParameters -> smpars,
49 fsModelParameters -> {
57 mq2Input -> MS^2 IdentityMatrix[3],
58 mu2Input -> MS^2 IdentityMatrix[3],
59 md2Input -> MS^2 IdentityMatrix[3],
60 ml2Input -> MS^2 IdentityMatrix[3],
61 me2Input -> MS^2 IdentityMatrix[3],
62 AuInput -> {{MS/TB, 0 , 0},
64 {0 , 0 , MS/TB + Xtt MS}},
65 AdInput -> MS TB IdentityMatrix[3],
66 AeInput -> MS TB IdentityMatrix[3]
70LinearRange[start_, stop_, steps_] :=
71 Range[start, stop, (stop - start)/steps];
77 { N[#], Sequence @@ MSSMEFTHiggsCalcMh[#, TB, Xtt] }&,
78 LinearRange[500, 10^4, 100]
81MhMin[{MS_, Mh_, DMh_}] := {MS, Mh - DMh};
82MhMax[{MS_, Mh_, DMh_}] := {MS, Mh + DMh};
83MhBest[{MS_, Mh_, DMh_}] := {MS, Mh};
90 PlotStyle -> {Red, Thick}];
92plot1 = ListLinePlot[{dataMhMax, dataMhMin},
93 PlotStyle -> LightGray,
94 Filling -> {1 -> {{2}, LightGray}},
98 BaseStyle -> {FontSize -> 16, FontFamily -> "Helvetica"},
99 PlotLabel -> Style["\*SubscriptBox[X, t] = 2.44949 \*SubscriptBox[M, S], tan\[Beta] = 5"],
100 PlotRange -> Automatic,
101 Axes -> False, Frame -> True,
102 FrameLabel -> {Style["\*SubscriptBox[M, S] / GeV"],
103 Style["\*SubscriptBox[M, h] / GeV"]}];
105Export["MSSMEFTHiggs_Mh_MS.png", plot, ImageSize -> 600];
Get["models/MSSMEFTHiggs/MSSMEFTHiggs_librarylink.m"]
MhBest[{MS_, Mh_, DMh_}] Mh
Export["MSSMEFTHiggs_Mh_MS.png", plot, ImageSize -> 600]
|