//! Six-domain empirical datasets from Al-Mayahi (2026), Tables 2-7.
//!
//! CAB-certified extraction from:
//!   HeytingLean.Bridge.AlMayahi.TauEpoch.DomainData
//!
//! All 40 data rows transcribed verbatim from the source paper.
//!
//! Provenance: Lean 4 kernel-verified -> LambdaIR -> MiniC -> C -> Rust

/// Domain I: H0 cosmology measurement (Table 2)
pub struct H0Measurement {
    pub name: &'static str,
    pub h0: f64,
    pub sigma: f64,
    pub method: &'static str,
    pub log10_tau_method: f64,
}

/// 7 H0 measurements spanning CMB to gravitational waves.
pub const H0_ENSEMBLE: &[H0Measurement] = &[
    H0Measurement { name: "Planck 2018",    h0: 67.4,  sigma: 0.5,  method: "CMB / Early Universe", log10_tau_method: 13.1 },
    H0Measurement { name: "ACT DR6",        h0: 68.3,  sigma: 1.1,  method: "CMB / Early Universe", log10_tau_method: 13.1 },
    H0Measurement { name: "Freedman CCHP",  h0: 69.96, sigma: 1.9,  method: "TRGB + SN Ia",        log10_tau_method: 0.5 },
    H0Measurement { name: "HoLiCOW",        h0: 73.3,  sigma: 1.8,  method: "Strong Lensing",       log10_tau_method: 9.5 },
    H0Measurement { name: "SH0ES",          h0: 73.04, sigma: 1.04, method: "Cepheid + SN Ia",     log10_tau_method: -2.0 },
    H0Measurement { name: "SPARC BTFR",     h0: 75.1,  sigma: 2.3,  method: "Tully-Fisher",         log10_tau_method: 8.0 },
    H0Measurement { name: "GW170817",       h0: 70.0,  sigma: 10.0, method: "Gravitational Waves",  log10_tau_method: -1.5 },
];

/// Domain II: LHC Higgs coupling offset (Table 3)
pub struct LhcKappa {
    pub channel: &'static str,
    pub kappa_atlas: f64,
    pub kappa_cms: f64,
    pub delta_kappa: f64,
}

/// 7 Higgs coupling channels; all delta_kappa > 0 (ATLAS > CMS).
pub const LHC_KAPPA_TABLE: &[LhcKappa] = &[
    LhcKappa { channel: "ggF",      kappa_atlas: 1.04, kappa_cms: 1.02, delta_kappa: 0.02 },
    LhcKappa { channel: "VBF",      kappa_atlas: 1.09, kappa_cms: 0.98, delta_kappa: 0.11 },
    LhcKappa { channel: "WH",       kappa_atlas: 1.30, kappa_cms: 1.15, delta_kappa: 0.15 },
    LhcKappa { channel: "ZH",       kappa_atlas: 1.05, kappa_cms: 1.01, delta_kappa: 0.04 },
    LhcKappa { channel: "ttH",      kappa_atlas: 1.12, kappa_cms: 1.05, delta_kappa: 0.07 },
    LhcKappa { channel: "gamma",    kappa_atlas: 1.06, kappa_cms: 1.01, delta_kappa: 0.05 },
    LhcKappa { channel: "tau-tau",  kappa_atlas: 1.02, kappa_cms: 0.99, delta_kappa: 0.03 },
];

/// Domain III: NMR chemical shift measurement (Table 4)
pub struct NmrMeasurement {
    pub field_tesla: f64,
    pub shift_ppm: f64,
    pub sigma_inter: f64,
    pub source: &'static str,
}

/// 8 NMR field strengths with 1/B0 scaling of inter-laboratory scatter.
pub const NMR_TABLE: &[NmrMeasurement] = &[
    NmrMeasurement { field_tesla: 1.5,  shift_ppm: 4.79, sigma_inter: 0.031,  source: "BMRB" },
    NmrMeasurement { field_tesla: 3.0,  shift_ppm: 4.78, sigma_inter: 0.019,  source: "BMRB" },
    NmrMeasurement { field_tesla: 7.05, shift_ppm: 4.77, sigma_inter: 0.012,  source: "BMRB/in-house" },
    NmrMeasurement { field_tesla: 9.4,  shift_ppm: 4.77, sigma_inter: 0.0092, source: "in-house" },
    NmrMeasurement { field_tesla: 11.7, shift_ppm: 4.76, sigma_inter: 0.0074, source: "SDBS" },
    NmrMeasurement { field_tesla: 14.1, shift_ppm: 4.76, sigma_inter: 0.0061, source: "SDBS" },
    NmrMeasurement { field_tesla: 18.8, shift_ppm: 4.76, sigma_inter: 0.0050, source: "in-house" },
    NmrMeasurement { field_tesla: 23.5, shift_ppm: 4.76, sigma_inter: 0.0043, source: "projection" },
];

/// Domain IV: Neutron lifetime measurement (Table 5)
pub struct NeutronMeasurement {
    pub experiment: &'static str,
    pub lifetime_sec: f64,
    pub sigma: f64,
    pub mode: &'static str,
}

/// 9 neutron lifetime measurements (7 bottle + 2 beam).
pub const NEUTRON_TABLE: &[NeutronMeasurement] = &[
    NeutronMeasurement { experiment: "Mampe 1993",     lifetime_sec: 882.6,  sigma: 2.7,  mode: "bottle" },
    NeutronMeasurement { experiment: "Serebrov 2005",  lifetime_sec: 878.5,  sigma: 0.8,  mode: "bottle" },
    NeutronMeasurement { experiment: "Pichlmaier 2010",lifetime_sec: 880.7,  sigma: 1.8,  mode: "bottle" },
    NeutronMeasurement { experiment: "Steyerl 2012",   lifetime_sec: 882.5,  sigma: 2.1,  mode: "bottle" },
    NeutronMeasurement { experiment: "Arzumanov 2015", lifetime_sec: 880.2,  sigma: 1.2,  mode: "bottle" },
    NeutronMeasurement { experiment: "Pattie 2018",    lifetime_sec: 881.5,  sigma: 0.9,  mode: "bottle" },
    NeutronMeasurement { experiment: "Gonzalez UCNtau",lifetime_sec: 877.75, sigma: 0.34, mode: "bottle" },
    NeutronMeasurement { experiment: "Byrne 1996",     lifetime_sec: 889.2,  sigma: 4.8,  mode: "beam" },
    NeutronMeasurement { experiment: "Yue 2013",       lifetime_sec: 887.7,  sigma: 2.2,  mode: "beam" },
];

/// Domain V: Proton radius measurement (Table 6)
pub struct ProtonMeasurement {
    pub experiment: &'static str,
    pub radius_fm: f64,
    pub sigma: f64,
    pub method: &'static str,
    pub log10_tau_method: f64,
}

/// 6 proton radius measurements (2 muonic + 4 electronic).
pub const PROTON_TABLE: &[ProtonMeasurement] = &[
    ProtonMeasurement { experiment: "Pohl 2010",      radius_fm: 0.84184, sigma: 0.00067, method: "muonic",     log10_tau_method: -9.0 },
    ProtonMeasurement { experiment: "Antognini 2013", radius_fm: 0.84087, sigma: 0.00039, method: "muonic",     log10_tau_method: -9.0 },
    ProtonMeasurement { experiment: "Bernauer 2010",  radius_fm: 0.879,   sigma: 0.008,   method: "electronic", log10_tau_method: -6.0 },
    ProtonMeasurement { experiment: "Sick 2012",      radius_fm: 0.886,   sigma: 0.008,   method: "electronic", log10_tau_method: -6.0 },
    ProtonMeasurement { experiment: "Beyer 2017",     radius_fm: 0.8335,  sigma: 0.0095,  method: "electronic", log10_tau_method: -6.0 },
    ProtonMeasurement { experiment: "Xiong 2019",     radius_fm: 0.831,   sigma: 0.012,   method: "electronic", log10_tau_method: -6.0 },
];

/// Domain VI: Muon g-2 theory comparison (Table 7)
pub struct MuonMeasurement {
    pub source: &'static str,
    pub a_mu_times_1e11: f64,
    pub sigma: f64,
    pub log10_tau_method: f64,
}

/// 3 muon g-2 theory values (White Paper + 2 BMW lattice results).
pub const MUON_TABLE: &[MuonMeasurement] = &[
    MuonMeasurement { source: "White Paper 2020", a_mu_times_1e11: 11659181.0, sigma: 4.3, log10_tau_method: -23.0 },
    MuonMeasurement { source: "BMW Lattice 2020", a_mu_times_1e11: 11659207.0, sigma: 3.4, log10_tau_method: -18.0 },
    MuonMeasurement { source: "BMW Lattice 2021", a_mu_times_1e11: 11659205.0, sigma: 5.1, log10_tau_method: -17.5 },
];

#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn test_h0_count() {
        assert_eq!(H0_ENSEMBLE.len(), 7);
    }

    #[test]
    fn test_lhc_all_positive() {
        assert!(LHC_KAPPA_TABLE.iter().all(|k| k.delta_kappa > 0.0));
    }

    #[test]
    fn test_nmr_count() {
        assert_eq!(NMR_TABLE.len(), 8);
    }

    #[test]
    fn test_neutron_count() {
        assert_eq!(NEUTRON_TABLE.len(), 9);
    }

    #[test]
    fn test_proton_count() {
        assert_eq!(PROTON_TABLE.len(), 6);
    }

    #[test]
    fn test_muon_count() {
        assert_eq!(MUON_TABLE.len(), 3);
    }

    #[test]
    fn test_total_data_rows() {
        let total = H0_ENSEMBLE.len() + LHC_KAPPA_TABLE.len() + NMR_TABLE.len()
            + NEUTRON_TABLE.len() + PROTON_TABLE.len() + MUON_TABLE.len();
        assert_eq!(total, 40);
    }
}