#!/usr/bin/env python3
"""
verify_P118.py -- Addendum 118: rho-omega isospin mixing.

This verifier checks the omega VMD correction in
118_Addendum_RhoOmegaMixing.tex: SU(3) baseline 72*pi, charge-factor ratio
C_rho/C_omega=3, epsilon=-alpha_s/3, f_omega^2=72*pi/(1-alpha_s)^2,
residuals, and the comparison to the neutral-vector table.

The omega arithmetic reproduces. Flagged issues are proof/status and sign
convention problems: the rho row's -2.3% sign is opposite the TOE-minus-PDG
convention, the "complete f_V^2 sector" wording conflicts with the paper's
own rho/open-items section, and the epsilon=-alpha_s/3 and MZ-scale choices
are asserted/imported rather than derived here.
"""

from __future__ import annotations

import math
import sys
from pathlib import Path

sys.path.insert(0, str(Path(__file__).resolve().parent))
from verify_common import CheckResult, Verifier


class ModernVerifier(Verifier):
    """Local adapter: tolerance logic inherited byte-identical from
    verify_common.Verifier; only the output layer is modernised."""

    def __init__(self, name: str) -> None:
        self.name = name
        self.results = []
        print(name)

    def record(self, label, ok, computed="", claimed="", detail=""):
        self.results.append(CheckResult(label, ok, computed, claimed, detail))
        n = len(self.results)
        desc = f"{label} -- {detail}" if detail else label
        print(f"  [{'PASS' if ok else 'FAIL'}] {n:>2}. {desc}")
        if computed != "" or claimed != "":
            print(f"        computed: {computed}")
            print(f"        claimed : {claimed}")
        return ok

    def summary(self):
        passed = sum(r.ok for r in self.results)
        failed = len(self.results) - passed
        print(f"\n{'='*60}\nRESULT: {passed} PASS / {failed} FAIL")
        return 1 if failed else 0


v = ModernVerifier("P118 -- rho-omega Mixing")

ROOT = Path(__file__).resolve().parents[1]
TEX = (ROOT / "118_Addendum_RhoOmegaMixing.tex").read_text()

PI = math.pi
MU0 = 4.0 * PI**3 + PI**2 + PI
MU1 = 16.0 * PI**3 / 5.0 + 3.0 * PI**2 / 4.0 + 2.0 * PI / 3.0
MU = MU1 / MU0
ALPHA_S = 0.1186

PDG_RHO = 24.56
PDG_OMEGA = 290.97
PDG_PHI = 179.1

C_RHO2 = 1.0 / 2.0
C_OMEGA2 = 1.0 / 18.0
C_PHI2 = 1.0 / 9.0
C_RATIO = math.sqrt(C_RHO2 / C_OMEGA2)

F_RHO = 8.0 * PI
F_OMEGA_SU3 = 72.0 * PI
EPSILON = -ALPHA_S / C_RATIO
F_OMEGA = F_OMEGA_SU3 / (1.0 + C_RATIO * EPSILON) ** 2
F_PHI = 36.0 * PI / MU**2
OMEGA_GAP_REL_SU3 = 100.0 * (PDG_OMEGA - F_OMEGA_SU3) / F_OMEGA_SU3
OMEGA_GAP_REL_EXP = 100.0 * (F_OMEGA_SU3 - PDG_OMEGA) / PDG_OMEGA
OMEGA_RESIDUAL = 100.0 * (F_OMEGA - PDG_OMEGA) / PDG_OMEGA
OMEGA_GAP_CLOSED_REDUCTION = 100.0 - 100.0 * abs(F_OMEGA - PDG_OMEGA) / abs(PDG_OMEGA - F_OMEGA_SU3)
RHO_RESIDUAL_TOE_MINUS = 100.0 * (F_RHO - PDG_RHO) / PDG_RHO
PHI_RESIDUAL = 100.0 * (F_PHI - PDG_PHI) / PDG_PHI
P114_DELTA = (8.0 / 9.0) * ALPHA_S / PI


v.check("omega SU3 baseline 72*pi", F_OMEGA_SU3, 226.2, rel=3e-5)
v.check("omega exp/SU3 gap ratio", PDG_OMEGA / F_OMEGA_SU3, 1.286, rel=3e-4)
v.check("omega original gap percent relative SU3", OMEGA_GAP_REL_SU3, 28.6, rel=2e-3)
v.check("omega original gap percent relative exp", OMEGA_GAP_REL_EXP, -22.0, rel=2e-2)
v.check("rho charge factor squared", C_RHO2, 1.0 / 2.0, rel=1e-12)
v.check("omega charge factor squared", C_OMEGA2, 1.0 / 18.0, rel=1e-12)
v.check("phi charge factor squared", C_PHI2, 1.0 / 9.0, rel=1e-12)
v.check("C_rho/C_omega", C_RATIO, 3.0, rel=1e-12)
v.check("epsilon", EPSILON, -0.03953, rel=1e-4)
v.check("3 epsilon equals -alpha_s", 3.0 * EPSILON, -ALPHA_S, rel=1e-12)
v.check("suppression factor 1-alpha_s", 1.0 + C_RATIO * EPSILON, 0.8814, rel=1e-12)
v.check("denominator (1-alpha_s)^2", (1.0 - ALPHA_S) ** 2, 0.7769, rel=5e-5)
v.check("omega f_V^2 corrected", F_OMEGA, 291.2, rel=2e-4)
v.check("omega residual percent", OMEGA_RESIDUAL, 0.08, abs_tol=0.02)
v.check("omega gap closed by residual reduction", OMEGA_GAP_CLOSED_REDUCTION, 99.7, rel=2e-4)
v.record(
    "mixing magnitude is in the quoted experimental range",
    0.04 <= abs(EPSILON) <= 0.05 or abs(abs(EPSILON) - 0.04) < 0.001,
    computed=abs(EPSILON),
    claimed="|epsilon| approx 0.04--0.05",
)
v.check("P114 Peirce QCD correction percent", 100.0 * P114_DELTA, 3.356, rel=3e-4)
v.check("alpha_s/3 percent", 100.0 * ALPHA_S / 3.0, 3.953, rel=9e-5)
v.check("phi f_V^2 table value", F_PHI, 179.7, rel=7e-4)
v.check("phi residual percent", PHI_RESIDUAL, 0.36, rel=1e-1)

v.check(
    "rho table residual sign under TOE-minus-PDG convention",
    RHO_RESIDUAL_TOE_MINUS,
    -2.3,
    rel=5e-2,
    detail="Expected fail: 8*pi is above 24.56, so TOE-minus-PDG is about +2.33%, not -2.3%.",
)
v.record(
    "complete f_V^2 sector is internally closed in P118",
    False,
    computed="the paper's open-items section says the rho residual has not yet been closed and OP1-NLO remains open",
    claimed="complete neutral vector meson f_V^2 sector fully derived",
    detail="Expected status fail.",
)
v.record(
    "epsilon = -alpha_s/3 is derived from a mixing Hamiltonian",
    False,
    computed="P118 gives a proportionality/sign argument but no explicit rho-omega mass/mixing matrix or derivation of the alpha_s/3 normalization",
    claimed="rho-omega mixing from J1 QCD coupling",
    detail="Expected proof-status fail.",
)
v.record(
    "MZ scale choice is derived in this addendum",
    False,
    computed="MZ is imported as the P106/P114 canonical anchor; no local derivation selects MZ over the hadronic scale",
    claimed="use alpha_s(MZ), not alpha_s(m_rho)",
    detail="Expected scale-selection fail.",
)
v.record(
    "rho residual open item is present",
    "has not yet been closed" in TEX,
    computed="open rho residual caveat found",
    claimed="rho correction still open at P118 stage",
)

sys.exit(v.summary())