Code
- SMPoleMatching.m
- Mathematica code for the matching of the pole masses to the MS parameters in the Standard Model
- SMrunning.m
- Mathematica code for evolving the MS coupling
constants above the top quark mass using teh Standrd Model beta functions
up to 3 loops.
- Example.nb
- Examples of usage.
- results-MAPLE.mpl
- Maple code for the pole mass matching from M. Kalmykov
Examples
Mathematica code
In[1]:= <<"SMPoleMatching.m"
In[2]:= mh = 130;
In[3]:= mt = 172.9;
In[4]:= aSMZ = 0.1184;
In[5]:= crule = constantsQCD3top2higgs2[mh, mt, aSMZ, mt]
Out[6]= {g10 -> 0.358729, g20 -> 0.648382, g30 -> 1.1646,
yt0 -> 0.934456, lambda0 -> 0.136494, mu0 -> 172.9}
In[7]:= lambda0/.crule
Out[7]= 0.136494
In[8]:= yt0/.crule
Out[8]= 0.934456
In[9]:= sol1 = solveSMRG[beta3loop, constantsQCD3top2higgs2[mh0, mt0, aS0, mt0]]
Out[9]= {g1 -> InterpolatingFunction[{{0., 50.}], <>],
g2 -> InterpolatingFunction[{{0., 50.}], <>],
g3 -> InterpolatingFunction[{{0., 50.}], <>],
yt -> InterpolatingFunction[{{0., 50.}], <>],
λ -> InterpolatingFunction[{{0., 50.}], <>],
mu0 ->172.9}
In[10]:= yt[Log[500/mu0]] /. sol1
Out[10]= 0.882464
In[11]:= LogLinearPlot[λ[Log[mu/mu0]] /. sol1, {mu, 200, 1000 mP}]
Out[11]=
