The flavio [94] and eos [95] software packages are used to fit the configuration (i) results, varying the parameter ℛ𝑒(𝒞9). For the flavio package, the default values of the SM nuisances are used, including QCD effects as in Refs. [13,27], while for the eos package, values are taken from Ref. [96]. The 3.3𝜎 discrepancy with respect to the SM value of ℛ𝑒(𝒞9) obtained with flavio in Ref. [5] becomes 3.6𝜎 ( 3.8𝜎) with flavio (eos). Including the differential branching fraction, the discrepancy increases to 4.1𝜎 ( 4.0𝜎). Neglecting the 6.0<𝑞2<8.0 GeV2/𝑐4 bin that is closest to the charmonium resonances reduces the observed significance of ℛ𝑒(𝒞9) by about 0.2𝜎. The best fit to the angular distribution is obtained with a shift in the SM value of ℛ𝑒(𝒞9) by −0.94+0.21−0.17 ( −0.92 ±0.23) with flavio (eos). When combining the information from the angular observables with the branching fraction, a shift in ℛ𝑒(𝒞9) of −0.93+0.18−0.16 ( −0.94 ±0.22) is obtained with flavio (eos). Differences between these numbers and those reported in Refs. [44,45], which fit for the Wilson coefficients directly, are due to the model inputs used, including the use of LHCb data to constrain the size of nonlocal contributions.