Replication Challenge - ICHEP 2026 MH Validation

An independent verifier can reproduce the core numerical results of this validation without access to any software beyond Python 3.7+ standard library.

Step 1: Run the scorer

cd 04_EXECUTABLE_SCORER/
python residual_calc.py

Compare all values to expected-output.txt in this directory.

Step 2: Verify four key numerical claims independently

Claim 1: kappa_lambda,MH ~ 0.993

phi = (1+5**0.5)/2
lambda_MH = phi**(-2*(2**0.5)) / 2           # 0.12819...
lambda_SM = 125.08**2 / (2 * 246.22**2)      # 0.12903...
kappa = lambda_MH / lambda_SM                 # 0.9932...

Claim 2: 172.2 GeV inside [172.19, 172.85] GeV

m_t = 172.52; unc = 0.33
lo = m_t - unc; hi = m_t + unc
print(lo, "<", 172.2, "<", hi, ":", lo <= 172.2 <= hi)

Claim 3: MH lightest mass inside JUNO substrate-consistent bound

import math
phi = (1+5**0.5)/2; alpha = 1/137.035989
me_meV = 0.5109989e9
mv3 = me_meV * alpha**3 * phi**(-9)           # 2.612 meV
mv2 = me_meV * alpha**3 * phi**(-6)           # 11.07 meV
dm2_21 = (mv2/1000)**2 - (mv3/1000)**2        # eV^2
bound = math.sqrt(61 * dm2_21) / 30 * 1000   # meV
print("mv3:", round(mv3,3), "< bound:", round(bound,4))

Claim 4: Smv = 60.47 meV inside [50, 140] meV 1-sigma detection band

mv1 = me_meV * alpha**3 * phi**(-3)
mv2 = mv1 * phi**(-3)
mv3 = mv2 * phi**(-3)
Smv = mv1 + mv2 + mv3
print("Smv:", round(Smv,2), "in [50, 140]:", 50 <= Smv <= 140)

Step 3: Verify ATLAS di-Higgs interval claim

From arXiv:2607.26879 abstract:

kappa_lambda,MH = 0.993
ATLAS 95% CL observed: [-3.4, 1.6]
print(-3.4 <= 0.993 <= 1.6)   # True