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_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..2fad37ce 100644
--- a/flavio/data/parameters_uncorrelated.yml
+++ b/flavio/data/parameters_uncorrelated.yml
@@ -42,6 +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.3937(81) # HFLAV 2022 all charmonium w/o LHCb + LHCb 2309.09728 psiKS Run 1&2
 
 
 # Meson decay constants
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)'),