corrections to m_z and m_w about theoretical-physics HOT 1 OPEN

A better script is:

from math import pi, log, sqrt, sin, cos, asin

eV = 1.
KeV = 1e3 * eV
MeV = 1e6 * eV
GeV = 1e9 * eV

M = 76500 * MeV
M0 = 87640 * MeV
m_u = 250 * MeV
m_d = 250 * MeV
m_e = 0.511 * MeV
m_mu = 105.6 * MeV
m_tau = 1777 * MeV
g = sqrt(0.039 * pi**2)
theta = asin(sqrt(0.238))
s_t = sin(theta)
c_t = cos(theta)
#m_H = 3 * GeV
#m_H = 200 * GeV
m_H = 1000 * GeV

#print g**2/(384*pi**2) * M * log(m_H**2) / MeV
#print g**2/(384*pi**2) * M0 * (1 + 10*s_t**2/c_t**2)* log(m_H**2) / MeV

delta_M_e = - g**2 * M**2 / pi**2 * (5./36 + 1./48 * \
                (2*log(m_u**2/M**2) + log(m_d**2/M**2) + log(m_e**2/M**2)))
#delta_M_mu = - g**2 * M**2 / pi**2 * (5./36 + 1./48 * \
#                (2*log(m_u**2/M**2) + log(m_d**2/M**2) + log(m_mu**2/M**2)))
#delta_M_tau = - g**2 * M**2 / pi**2 * (5./36 + 1./48 * \
#                (2*log(m_u**2/M**2) + log(m_d**2/M**2) + log(m_tau**2/M**2)))
delta_M0 = - g**2 * M0**2 / (pi**2 * c_t**2) *  ( \
        1./108 * (40*c_t**4 - 50*c_t**2 + 25  ) + \
        1./72  * (8 *c_t**4 - 10*c_t**2 + 5   ) * log(m_u**2/M0**2) + \
        1./72  * (2 *c_t**4 -    c_t**2 + 1./2) * log(m_d**2/M0**2) + \
        1./48  * (4 *c_t**4 - 6 *c_t**2 + 3   ) * log(m_e**2/M0**2))

M_higgs = 120 * MeV
M0_higgs = 140 * MeV

print 0.5 * delta_M_e / M / MeV
print 0.5 * delta_M0 / M0 / MeV

print "bare masses:"
print M / GeV
print M0 / GeV
print "shifted masses:"
print (M + 0.5 * delta_M_e / M + M_higgs) / GeV
print (M0 + 0.5 * delta_M0 / M0 + M0_higgs) / GeV

print "shifted masses (correct):"
print (M + 3080 * MeV) / GeV
print (M0 + 3310 * MeV) / GeV

print "experimental values:"
print 80.399
print 91.1876

Which produces:

1600.76296787
1702.62148258
bare masses:
76.5
87.64
shifted masses:
78.2207629679
89.4826214826
shifted masses (correct):
79.58
90.95
experimental values:
80.399
91.1876

from theoretical-physics.