where \(\mathcal {O}_{\rm th}(C_i,C_j)\)  are the theoretical predictions of N = 168 observables used in the \(\chi ^2\)  fit, \(\mathcal {O}_{\rm exp}\)  are the corresponding central values of the experimental measurements and \(\mathcal {C}\)  is the total covariance matrix. This \(N \times N\)  matrix is constructed by adding the individual theoretical and experimental covariance matrices.
The experimental correlations are included for the angular observables in \(B^0 \rightarrow K^{*0} \mu ^+ \mu ^-\)  [1]}, \(B^+ \rightarrow K^{*+} \mu ^+ \mu ^-\)  [2]}
and \(B_s \rightarrow \phi \mu ^+ \mu ^-\)  [3]} whereas the theoretical covariance matrix includes form factors and power corrections uncertainties. These are computed using flavio [4]} where the observables are pre-implemented based on refs. [5]}, [6]}.
