From ae6a36eb8b99a62fa2ad9102f968e39af31685d4 Mon Sep 17 00:00:00 2001 From: Peter Stangl Date: Sat, 16 Dec 2023 15:19:50 +0100 Subject: [PATCH 1/3] add beta parameter --- flavio/data/parameters_metadata.yml | 3 +++ flavio/data/parameters_uncorrelated.yml | 2 ++ 2 files changed, 5 insertions(+) diff --git a/flavio/data/parameters_metadata.yml b/flavio/data/parameters_metadata.yml index a12bb464..e767fef0 100644 --- a/flavio/data/parameters_metadata.yml +++ b/flavio/data/parameters_metadata.yml @@ -42,6 +42,9 @@ Vcb: gamma: tex: $\gamma$ description: CKM Unitarity Triangle angle $\gamma$ in radians +beta: + tex: $\beta$ + description: CKM Unitarity Triangle angle $\beta$ in radians # Quark masses diff --git a/flavio/data/parameters_uncorrelated.yml b/flavio/data/parameters_uncorrelated.yml index e5720f5b..8bdf1c41 100644 --- a/flavio/data/parameters_uncorrelated.yml +++ b/flavio/data/parameters_uncorrelated.yml @@ -42,6 +42,8 @@ Vus: 0.2248(6) # FLAG 2023 Nf=2+1+1 Vub: 3.64(16)e-3 # FLAG 2023 Nf=2+1 B->pilnu Vcb: 41.97(48) 1e-3 # 2310.20324 gamma: 1.15(6) # HFLAV 2022 B -> DK +#beta: 0.387(12) # HFLAV 2022 B -> J/psi K +beta: 0.405(7) # LHCb 2309.09728 (Run 1 & 2)) # Meson decay constants From 5ab22f6fbd9510b7e87fcd7806d231caf67f136b Mon Sep 17 00:00:00 2001 From: Peter Stangl Date: Mon, 15 Jan 2024 19:42:04 +0100 Subject: [PATCH 2/3] update value of beta and use beta gamma CKM scheme as default --- flavio/data/config.yml | 2 +- flavio/data/parameters_uncorrelated.yml | 3 +-- 2 files changed, 2 insertions(+), 3 deletions(-) diff --git a/flavio/data/config.yml b/flavio/data/config.yml index 4aa6efaa..4a544a7c 100644 --- a/flavio/data/config.yml +++ b/flavio/data/config.yml @@ -93,7 +93,7 @@ RGE thresholds: implementation: # CKM matrix parametrization - CKM matrix: Tree + CKM matrix: beta gamma # Unitarity relation in Kaon mixing (either 'ut' or 'ct') K mixing unitarity: ut diff --git a/flavio/data/parameters_uncorrelated.yml b/flavio/data/parameters_uncorrelated.yml index 8bdf1c41..2fad37ce 100644 --- a/flavio/data/parameters_uncorrelated.yml +++ b/flavio/data/parameters_uncorrelated.yml @@ -42,8 +42,7 @@ Vus: 0.2248(6) # FLAG 2023 Nf=2+1+1 Vub: 3.64(16)e-3 # FLAG 2023 Nf=2+1 B->pilnu Vcb: 41.97(48) 1e-3 # 2310.20324 gamma: 1.15(6) # HFLAV 2022 B -> DK -#beta: 0.387(12) # HFLAV 2022 B -> J/psi K -beta: 0.405(7) # LHCb 2309.09728 (Run 1 & 2)) +beta: 0.3937(81) # HFLAV 2022 all charmonium w/o LHCb + LHCb 2309.09728 psiKS Run 1&2 # Meson decay constants From ec9312b6cdecb73def7b62cf29885bbd0b2b1283 Mon Sep 17 00:00:00 2001 From: Peter Stangl Date: Thu, 15 Feb 2024 22:16:46 +0100 Subject: [PATCH 3/3] update tests --- flavio/physics/test_ckm.py | 7 ++++--- flavio/test_functions.py | 2 +- 2 files changed, 5 insertions(+), 4 deletions(-) diff --git a/flavio/physics/test_ckm.py b/flavio/physics/test_ckm.py index 88e6bf55..06931689 100644 --- a/flavio/physics/test_ckm.py +++ b/flavio/physics/test_ckm.py @@ -45,6 +45,7 @@ class TestCKM(unittest.TestCase): par_s = dict(t12=t12,t13=t13,t23=t23,delta=delta) par_w = dict(laC=laC,A=A,rhobar=rhobar,etabar=etabar) par_t = dict(Vus=Vus,Vub=Vub,Vcb=Vcb,gamma=gamma) + par_b = dict(Vus=Vus,Vcb=Vcb,beta=beta,gamma=gamma) def test_ckm_parametrizations(self): np.testing.assert_almost_equal(self.v_t/self.v_s, np.ones((3,3)), decimal=5) @@ -79,11 +80,11 @@ def test_ckm_xi(self): xi('x','bs') def test_ckm_angles(self): - c_gamma = get_ckmangle_gamma(self.par_t) + c_gamma = get_ckmangle_gamma(self.par_b) # angle gamma should be equal to input self.assertAlmostEqual(c_gamma/gamma, 1., places=3) - c_beta = get_ckmangle_beta(self.par_t) - c_alpha = get_ckmangle_alpha(self.par_t) + c_beta = get_ckmangle_beta(self.par_b) + c_alpha = get_ckmangle_alpha(self.par_b) # some of angles should be 180° self.assertEqual( degrees(c_alpha) + degrees(c_beta) + degrees(c_gamma), 180.) diff --git a/flavio/test_functions.py b/flavio/test_functions.py index 7aea30c3..11797ce1 100644 --- a/flavio/test_functions.py +++ b/flavio/test_functions.py @@ -48,7 +48,7 @@ def test_exp_combo(self): def test_get_dep_par(self): self.assertEqual( get_dependent_parameters_sm('BR(Bs->mumu)'), - {'DeltaGamma/Gamma_Bs', 'GF', 'Vcb', 'Vub', 'Vus', 'alpha_e', 'alpha_s', 'f_Bs', 'gamma', 'm_Bs', 'm_b', 'm_mu', 'm_s', 'tau_Bs', 'm_t'} + {'DeltaGamma/Gamma_Bs', 'GF', 'Vcb', 'beta', 'Vus', 'alpha_e', 'alpha_s', 'f_Bs', 'gamma', 'm_Bs', 'm_b', 'm_mu', 'm_s', 'tau_Bs', 'm_t'} ) self.assertEqual( get_dependent_parameters_sm('BR(B0->ee)'),