double EW_ABC::sigma0_had(const bool bAlternative) const { if (!bAlternative) return ( 12.0 * M_PI / NPE.getMz() / NPE.getMz() * Gamma_l(NPE.ELECTRON) * Gamma_had() / GammaZ(false) / GammaZ(false)); else { double delta_als = (NPE.Als(NPE.getMz(), FULLNNLO) - 0.119) / M_PI; double delta_alpha = (NPE.alphaMz() - 1.0 / 128.90) / NPE.getAle(); double sigma_h0 = 41.420 * (1.0 - 0.41 * delta_als + 0.03 * delta_alpha); return ( sigma_h0 * (1.0 - 0.03 * eps1() + 0.04 * eps3() - 0.20 * epsb())); } }
double EW_CHMN::sigma0_had() const { return ( 12.0 * M_PI / SM.getMz() / SM.getMz() * GammaZ_l(StandardModel::ELECTRON) * GammaZ_had() / GammaZ() / GammaZ()); }