\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_W$$\end{document}MW helps select \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document}Z′ models for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\rightarrow s \ell \ell $$\end{document}b→sℓℓ anomalies

As shown in Allanach et al. (Global fits of third family hypercharge models to neutral current B-anomalies and electroweak precision observables. arXiv:2103.12056), the Third Family Hypercharge (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y_3$$\end{document}Y3) Model changes the Standard Model prediction for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_W$$\end{document}MW whilst simultaneously explaining anomalies in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\rightarrow s\ell \ell $$\end{document}b→sℓℓ transitions via a heavy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document}Z′ gauge boson which is spawned by a spontaneously broken gauged \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U(1)_{Y_3}$$\end{document}U(1)Y3 symmetry. The 2022 CDF II measurement of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_W$$\end{document}MW, which is far from the Standard Model prediction in the statistical sense, somewhat disfavours the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y_3$$\end{document}Y3 model. Here, we generalise the gauge charge assignments to the anomaly-free combination \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$s Y_3 + t (B_3-L_3)$$\end{document}sY3+t(B3-L3) and show that incorporating the 2022 CDF II measurement of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M_W$$\end{document}MW selects a viable domain of integers s and t. For example, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$s=1, t=-3$$\end{document}s=1,t=-3 yields a p value of .08 in a two-parameter global fit to 277 electroweak and flavour changing b data, much improving a SM p value of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\times 10^{-6}$$\end{document}1×10-6.


Introduction
Recently, the CDF II Collaboration reported a measurement of the W boson mass M W = 80.4335 ± 0.0094 GeV [2] that disagrees by many sigma with the Standard Model (SM) prediction of the quantity. One average 1 of all current measurements yields [4] M W = 80.4133 ± 0.0080 GeV, (1.1) still in significant disagreement (with a pull of −6.4σ ) with SM predictions; we call this disagreement the M W anomaly. Note that the top quark mass m t plays an important role in the prediction of M W in the SM: throughout this paper, we impose the following constraint upon the top quark pole mass, m t = 171.79 ± 0.38 GeV, (1.2) based on a combination of the latest data [4]. We use the smelli2.3.2 defaults [5] for other electroweak constraints and parameters. There are also discrepancies between SM predictions and some measurements of flavour-changing B-meson decays, which we collectively call the b → s anomalies. For example, various lepton flavour universality (LFU) observables like the ratios of branching ratios R K ( * ) = B R(B → K ( * ) μ + μ − )/B R(B → K ( * ) e + e − ) are observed to be lower than their SM predictions in several channels and several different bins of di-lepton invariant mass squared [6][7][8]. Similar double ratio measurements in B 0 → K 0 s + − and B + → K * + + − decays [9] are consistent with a similar deficit in di-muons over di-electrons, albeit with lower statistics than in R ( * ) K . The predictions for all aforementioned double ratios have rather small theoretical uncertainties due to cancellations and the SM predictions are generically considered to be robust in the di-lepton invariant mass squared bins of interest. B R(B s → μ + μ − ) also has quite small theoretical uncertainties in its SM prediction. The combined measurements of B R(B s → μ + μ − ) are in a 2σ tension with its SM prediction [10][11][12][13][14]. Several other B-meson decay observables appear to be in tension with SM predictions even when their larger theoretical uncertainties are taken into account: for example angular distributions in B → K * μ + μ − decays [15][16][17][18][19][20], and B R(B s → φμ + μ − ) [21,22]. For these quantities though, there is room for argument about the assumed size of the theoretical uncertainties and indeed the best way of estimating the SM predictions for them.
We display the statistical tensions in the SM due to the M W and b → s anomalies in Table 1, as calculated by the computer program 2 smelli2.3.2 [5]. We see from the table that the b → s anomalies disfavour the SM: the b → s data yields a poor fit with a global p value of .0036, even when using the default smelli2.3.2 constraints upon the EWPOs (4.13), i.e. excluding the new CDF II measurement of M W (although the EWPOs have an acceptable fit in and of themselves with a p value of .17). Taking into account the CDF II measurement of M W as in (1.1) exacerbates an already poor quality of fit, lowering the global p value to 10 −6 .
Global fits find that new physics contributions to the (bγ μ P L s)(μγ μ P X μ) + H.c. vertex in the Lagrangian density can ameliorate the fit to the b → s data [23][24][25][26][27][28][29], where P L is the left-handed projection operator. P X corresponds to the helicity projection of the muon pair in the effective vertex: the fits agree that a purely right-handed projection P X = P R is disfavoured, whereas a mixture in the domain P X ≈ P L to P X ≈ P L + P R is preferred [30]. Whereas the fitting groups yield very similar results when fitting just to the 'LFU' category of observables, there are some differences observed when the 'quarks' category is included: for example, which option out of P X = P L or P X = P L + P R has a better fit. Such differences in the constraints can arise from the treatment of theoretical uncertainties in the 'quarks' category. The lesson one learns from such studies is that, from a new physics point of view in order to fit the b → s anomalies, one needs a new physics state that couples to left-handed quark fields and left-handed muon fields in a family non-universal manner (but it may or may not also couple to right-handed muon fields and/or electron fields). According to Ref. [31], there is a mild preference for a family universal coupling of new physics to leptons as well as a specific and different new physics contribution to the coupling to di-muon pairs.
In Ref. [50], the Third Family Hypercharge (Y 3 ) model was presented and shown to fit the b → s data. Here, the X charges of all fermions are zero except for the third family, which has X charge equal to its hypercharge. The model predicts Z − Z mixing because the Higgs doublet is necessarily charged under U (1) X to allow a renormalisable top Yukawa coupling, which would otherwise be forbidden by the U (1) X gauge symmetry. It was noted that this will change the SM prediction for M W [73] and in Ref. [1] it was shown that the model could simultaneously fit the b → s anomalies and EWPOs: indeed, that the SM prediction for M W was improved, since it was some 2σ too low in the SM (compared to the measurement at the time) but receiving a positive correction from the Y 3 model.
We will show below that including the new CDF II measurement of M W as in (1.1), the Y 3 model is somewhat disfavoured with a global p value of .02. Our purpose in this paper is then to build a generalisation of the Y 3 model, which is similar in construction and as simple, but which provides a simultaneously acceptable fit to both the M W anomaly and b → s anomalies. To this end, in Sect. 2, we shall propose a generalisation of the Y 3 charge assignments, which is rendered anomaly-free simply by the inclusion of right-handed neutrinos. There are two classes of charge assignment, depending upon whether one permutes the righthanded third family field's charged lepton charge with that of the second family (i.e. X e 3 ↔ X e 2 ) or not. In Sect. 3, we calculate the SM effective field theory (SMEFT) coefficients that arise from integrating out the heavy Z boson, and then matching to the SMEFT at tree level. This allows us to encode  (4.14) and the discussion of it in the surrounding text our model in a suitable form for input into smelli2.3.2 and use its calculation of EWPOs, LFU observables and the 'quarks' category of observables. We discuss the phenomenology of the models in Sect. 4, in particular the effect on M W and on b → s coefficients in the weak effective theory. The case where one permutes X e 3 ↔ X e 2 corresponds to P X = P L + P R , whereas in the case where one does not, one obtains a certain linear combination of P L and P R , depending upon s and t. In both cases, there is a purely axial coupling to electrons. We then present global fits to this latter family of model. The inclusion of the 2022 CDF II M W measurement selects a subset of models which provide acceptable fits to the collective data. We conclude in §5.

The models
We consider a class of Z models based on extending the SM gauge symmetry by a U (1) X factor, where  Table 2. The X charge assignment in (2.1) is the most general anomaly-free U (1) X extension 4 that couples only to a single family of SM fermions [38], including a right-handed neutrino. While we have parameterised this family of U (1) X models by two integers s and t, it is only the rational parameter t/s ∈ Q that is relevant for phenomenology. 5 This family of U (1) X extensions of the SM are known to have semi-simple gauge completions without needing any extra chiral fermions, as was shown in Ref. [74].
The Z model is designed to explain the b → s anomalies. Global fits strongly favour a lepton flavour non-universal coupling of the Z to left-handed muons, at least, as described in Sect. 1. We thus permute the non-zero left-handed lepton charge from the third to the second family, as in the original Y 3 model of [50]. Regarding right-handed leptons, acceptable fits can be obtained with or without permuting the non-zero right-handed lepton charge to the second family.
In summary, the charges of the SM fermions, together with the SM Higgs doublet H = (H + , H 0 ) T and an extra U (1) X symmetry-breaking scalar θ which is a SM singlet, are where X θ = 0, with all other U (1) X charges being zero. The original Y 3 model of Ref. [50], which is a template for the family of models we consider here, can be recovered by setting n = 3 and (s, t) = (1, 0), thus decoupling the righthanded neutrino.
Like the Y 3 model [50], the family of models we consider here allow, on the quark side, only third family Yukawa couplings to the Higgs at the renormalisable level, L ⊃ y t q 3 H c u 3 + y b q 3 Hd 3 . Thus, to zeroth order, gauging the anomaly-free chiral third-family symmetry (2.5) postdicts a heavy third family and small quark mixing angles, as observed. The light quark Yukawa couplings, responsible for the masses of the first and second generation and for the quark mixing angles, must come from higher-dimensional operators. Such operators can come from integrating out heavier fermionic representations that are vector-like under to the gauge group. These are details of the UV theory which it Table 2 Representations of fields under the SM gauge factors, which are family universal, together with the family non-universal Y 3 and (B − L) 3 symmetries on which our Z model is based.
We use the minimal integer normalisation for the charges under each U (1) factor. Note that the permutations 2 ↔ 3 (and in some cases e 2 ↔ e 3 ) are made after the assignments shown. All fields are Weyl fermions except for the complex scalar Higgs doublet H and the complex scalar flavon θ. * denotes that the charge is a non-zero number whose value does not change any of the discussion or results of this paper would be premature to specify (although see Ref. [1] for further comments and ideas).
On the other hand, the lepton sector is not so natural from the Yukawa perspective, due to the fact that we permute the non-zero X i charge into the second family. At face value, the charge assignment (2.5) with n = 3 allows only for a renormalisable off-diagonal Yukawa coupling ∼ 2 He 3 , which must be highly suppressed in order to explain the relative heaviness of the tau and the non-observation of μ − τ lepton flavour violation. This is easily explained with a little more model-building; for example, additionally gauging an anomaly-free lepton-flavoured U (1) symmetry could ban this off-diagonal Yukawa coupling. 6 In what follows we use a convention in which the covariant derivative acting on a field f is where X μ is the gauge field for U (1) X and g X is its gauge coupling. We assume that any kinetic mixing between U (1) X and U (1) Y gauge fields is negligible, leaving a study of such effects to future work. The fermion couplings to the gauge fields, in the gauge eigenbasis (indicated by 'primes'), are where the sum on ψ runs over all SM Weyl fermions and the charges X ψ are those in (2.5).

Symmetry breaking and Z − Z mixing
The U (1) X symmetry is broken predominantly by the SM singlet scalar field θ , which acquires a vacuum expectation This implies that the mass of the Z field M Z ≈ g X X θ v X . However, because the SM Higgs is also charged under U (1) X , there is mass mixing between the Z boson and the Z gauge boson. These physical gauge bosons are linear combinations of the SM combination 7 is the Weinberg angle) and the X μ field viz.
where the mixing angle α z is determined by, at tree-level, At a constant value of M Z (we shall take M Z to be an experimental input), the O(M 2 Z /M 2 Z ) correction in (2.10) translates to an upward shift in the SM prediction for M W , as we shall describe in Sect. 4.1.

Fermion mixing matrices
Following Refs. [1,50], we assume a simple ansatz for the 3by-3 unitary fermion mixing matrices describing the change from the gauge eigenbasis to the (unprimed) mass eigenbasis of the fermionic fields. The purpose of this ansatz is to characterise the main physical flavour characteristics of the model without introducing large flavour-changing neutral currents that would be subject to strong experimental constraints. We may think of the ansatz as a limit to expand around: the fairly strong assumptions might be motivated by further model-building involving additional symmetries or dynamics, but we leave such considerations to one side, for now.
For left-handed down-type quarks, we parameterise the mixing matrix as (2.11) For simplicity, and to avoid unwanted large contributions to F = 2 processes and charged lepton flavour violating processes, we choose Using these rotation matrices, the couplings of the Z boson to the physical fermion states can be obtained from (2.7). For the left-handed quarks, we work in a 'downaligned' basis in which the (unprimed) quark doublets are For the left-handed leptons, we work in a basis in which the left-handed charged leptons e L i align with the mass eigenstates, and thus the (unprimed) lepton doublets are i = (U † i j ν L , j , e L ,i ). Since V e L = V e R = 1 the charged lepton couplings remain diagonal. The down quark couplings are mixed away from the diagonal, however, as is indeed necessary to obtain a quark flavour-violating coupling of the Z to bs andsb. We have: 2 + X e n e n γ μ e n + X ν n ν n γ μ ν n ,

SMEFT matching
Integrating out the heavy X boson and matching onto the SM Effective Field Theory (SMEFT) at a scale M X := g X X θ v X , we obtain the Wilson coefficients (WCs) for dimension-6 SMEFT operators written in We write the coefficients as functions of the charges X f , which are explicitly parameterised in (2.5). The integer n = 2 or n = 3 corresponds to two variations of the model, as explained in Sect. 2. All Wilson coefficients are in units of g 2 X /M 2

X 4 Phenomenology
Starting from the SMEFT matching of Sect. 3, we will use the smelli2.3.2 program to evaluate the likelihood of the model given hundreds of observables in the electroweak and flavour sectors. Before we do so, however, we think it important to highlight the most important observables that are sensitive to our model, and how these depend upon the SMEFT coefficients and hence upon the X f charges.

Electroweak
In light of the recent CDF II measurement [2], the M W prediction is deserving of special attention. We can parameterise its deviation from the SM prediction in our Z model via the parameter where the parameterρ = 1.01019 ± 0.00009 includes the custodial-violating top contributions to the gauge boson masses (see Ref. [75]). In (4.1), we use the conventional notationĉ 2 to denote the cosine squared of the renormalised Weinberg angle in the MS scheme. The ρ 0 parameter is defined so as to equal precisely 1 in the SM using the MS scheme. Its deviation from unity in our heavy Z model, which recall is due to Z − Z mixing, is Importantly, ρ 0 unavoidably shifts upwards [1,77], easing the tension due to the CDF II M W measurement irrespective of the sign of the Higgs charge X H . At the level of the SMEFT, this shift is captured by the Wilson coefficients C φ D and C φ , which are (as in Table 3) Global fits of the SMEFT to electroweak data that turn on only these two operators have been found to give a good fit to electroweak data in light of the CDF measurement (e.g. [78]). The status of the electroweak fit is of course much more delicate in our model for the b → s anomalies, because we have a plethora of new physics effects in Z pole observables due to modified Z couplings to fermions. Important among these are modified Z couplings to leptons, especially muons, and the forward-backward asymmetry variable A F B b for Z decays to bb, because these observables already exhibit small tensions with the SM and because they receive large-ish corrections in our model. These effects are strictly correlated to the shifts in C φ D and C φ , since they arise from the Z − Z mixing.

b → s
Regarding the b → s anomalies and related observables, the new physics effects can be parameterised by contributions to Wilson coefficients in the weak effective theory (WET). The WET Hamiltonian is, using a conventional normalisation, where we emphasise that in the present paper, C i denotes the beyond the SM (BSM) contribution to the Wilson coefficient. When coupling tobs/sb currents, the Z only couples to the left-handed component and so there are contributions to the dimension-6 semi-leptonic operators O 10 = e 2 16π 2 (sγ μ P L b)(lγ μ γ 5 ). (4.5) Because of the Z − Z mixing, the Z boson picks up a small quark flavour-violating coupling tobs, viz. L ⊃ g sb Z Z μbL γ μ s L + H.c., where the coupling g sb Z is proportional to the Z − Z mixing angle (2.9). Specifically, There are thus BSM contributions to the 4-fermion operators O 9,10 from both Z exchange and from Z exchange, and these are of the same order. The Z contributions are lepton flavour universal (LFU), whereas the Z contributions introduce lepton flavour universality violation (LFUV). Accounting for both contributions 8 we have where recall that k = 1/(1 − 4 sin 2 θ w ) ≈ 9.23. We have defined the common pre-factor where e is the electromagnetic gauge coupling. The Zinduced LFU pieces in (4.8) clearly vanish in the limit X H → 0, in which there is no Z − Z mixing. The LFU contributions to the b → s observables are therefore correlated to their effects on EWPOs, including the shift in M W . This correlation was explored, prior to the updated CDF II measurement, in Ref. [77]. Substituting in the parametrisation (2.2-2.5) of charges in our models, and dropping the 1/k suppressed C U 9 contributions, we have This is equivalent to writing the chirality of the coupling to muons in terms of left and right projection operators as in 8 Of course these contributions to the WET coefficients can also be derived at the level of the SMEFT matching. From this perspective, the SMEFT-WET matching formulae are [79,80]: (4.7) Substituting in the SMEFT Wilson coefficients from Table 3 reproduces (4.8).
Sect. 1, i.e. via an effective operator ∝ (sγ μ P L b)(μγ μ P X μ). For a general s and t, this projection operator is For the variation of the model with n = 2, we have (4.12) thus P X ≈ P L + P R , a vectorial coupling of the Z to μ + μ − .

Global fits
We wrote a computer program (included in the ancillary information attached to the arXiv preprint of this paper) that uses smelli2.3.2 to calculate the likelihoods of 277 observables for our model by first calculating the WCs in Table 3  For a given s and t the likelihoods are a function of two effective parameters: α := g X × 3 TeV/M X and the quark mixing angle θ sb . Throughout, we shall illustrate with M X = 3 TeV as an example (similar third-family type models with M X ≥ 1 TeV were not ruled out by search constraints in the parameter region where they fit the b → s anomalies [81,82] and so we expect M X = 3 TeV to be allowed by direct Z search constraints). We do not expect a significant change if we were to change M X to a different value M X as long as one scales g X by the same factor, since the boundary conditions upon the Wilson coefficients (shown in Table 3) depend only upon the ratio g X /M X . This approximation is good up to small loop suppressed corrections from the renormalisation group running between M X and M X , which induce relative multiplicative changes of O[log(M X /M X )/(16π 2 )] in corrections to predictions of observables 9 .
We consider each different value of the ratio t/s to constitute a different model. smelli2.3.2 then calculates χ 2 via the predicted values of observables. We minimise χ 2 by varying α and θ sb using the Nelder-Mead algorithm, given a guess at a starting point. This was obtained by doing a rough initial scan for one value of t = 0 and s = 0 and roughly estimating how electroweak and b-observable χ 2 values are expected to scale with s and t. The result, for a given value of t/s, is a best-fit point, where the fit in the EWPOs is balanced against those of the b data. This balance crucially depends upon the experimental constraint that is taken on M W , as we illustrate in the top panel of Fig. 1, where we display p values using (1.1), which includes the CDF II measurement. Here, we can see that, for example, the Third Family Hypercharge Model model, denoted Y 3 , is somewhat disfavoured: its global p value is around .02. Models at larger values of |t/s| approximate the B 3 − L 3 model and so the p values all asymptote towards the left-hand side and the right-hand side of the plot. The B 3 − L 3 model has no Z −Z mixing because the Higgs doublet field is not charged under B 3 − L 3 , meaning that in the limit of large |t/s|, the CDF II M W measurement strongly disfavours each model. We see that models with −5 < t/s < −2 all fare well with global p values above the canonical .05 bound shown by the red line. We notice by examining the EWPO p values that the M W constraint prefers a localised value of t/s. The other effects that are relevant are in the LFU and 'quarks' categories: both have a valley around t/s = −1. Here, at s = −t, one can see from (2.2) that X q 3 = X 2 = 0, meaning that there is no coupling at tree-level between the Z and lefthanded quarks or left-handed muons; the Z does therefore not help with the b → s anomalies and we revert to the poor fit of the SM both for the 'LFU' and for the 'quarks' category of observable. It is of interest that the p value is suppressed somewhat for all t/s by the 'quarks' category of observable, in which there is room for disagreement with the theoretical uncertainty budget of the prediction.
We see a rather weak M W selection effect in the bottom panel of Fig. 1  constraint leads to weaker selection than the one including the CDF II M W measurement because it has much larger uncertainties and less need for a large contribution from the Z . We see here that only the region −2 < t/s < 0 has a global p value of less than .05, and this is clearly driven by the 'quarks' and 'LFU' categories, not by EWPOs. We summarise the χ 2 and p values for both options of experimental M W constraint in Table 4.
We now go on to examine the pulls at each best-fit point for t/s = −3. We define the pull for an observable with theoretical prediction P, measured central value M and uncertainty C to be pull = P − M C , (4.14) where C does not include correlations with other observables, but may include theoretical uncertainties added in quadrature in some cases. In particular, for M W , we have added an estimated uncertainty in the prediction of 5.6 MeV [4] in quadrature. We display the pulls for the combination of M W measurements that include the recent CDF II determination in Fig. 2. Several notable effects are evident: for example, unsurprisingly M W itself is better fit, with a pull of −1. The observable R μ , the branching ratio of the Z 0 boson to μ + μ − , has been increased by a fair amount; in the SM the prediction was high by over 1σ , but in the Z model with t/s = −3 it is less than one sigma too low 10 . We see that the forward-backward asymmetry measured in e + e − → bb, A b F B , has a worse fit in the Z model. A e , the left-right asymmetry in e + e − → e + e − , has a smaller pull in the Z model. Overall, the quality of fit to 10 Observables such as R μ receive contributions from (C (1) φl ) 22 , C 22 φe . As a reference to Table 3 and (2.2-2.5) shows, the sign of these contributions depends upon the sign of t/s. We have checked that the predictions for such observables do indeed go in the 'wrong' direction for positive t/s, resulting in a less favourable fit. the EWPOs is fine: Table 4 shows that the p-value is .29. The right-hand panel of Fig. 2 shows that many of the selected observables of interest to b-meson decays are fit better than in the SM save for the B s − B s mixing observable m s , which receives (fairly mild) positive corrections from the Z contribution.
We compare and contrast the fits including the CDF II M W measurement from Fig. 2 with fits excluding it in Fig. 3. The most obvious effect of excluding the CDF II M W measurement is that the SM pull of M W of is only 2σ when the CDF II M W measurement is excluded. Because the required shift in M W is smaller, the relative effect upon the other EWPOs is smaller and the result is a good fit to EWPOs for t/s = −3: Table 4 reveals the p value in the fit to be .49. The fits in the b-observables on the right-hand panel show a very similar pattern between Figs. 2 and 3. Essentially, g X is being fixed by the EWPOs (and is being pulled by M W in particular), and then θ sb is fit to a value which fits the b → s anomalies. Here, many of the SM-discrepant observables relevant to the b-anomalies receive a relative contribution from the Z [50] ∝ g 2 X sin 2θ sb /M 2 X . We see that the pull of m s , which is a measure of the B s − B s mixing, decreases when one includes the CDF II M W measurement. This is because the fit is pushed to a larger value of g X in order to fit the larger needed new physics contribution in M W . In turn, in order to fit the b-anomalies, one requires a smaller value of θ sb in order to keep the Z contributions to them (4.9) constant, in turn reducing the Z contribution to B s − B s mixing.
The p values and pulls for the n = 2 fits are qualitatively similar to those for n = 3 and we neglect to present them here, noting that they are presented in the ancillary information attached to the arXiv preprint version of this paper. Excluding all of the M W measurements except for the CDF II one, the fits do not differ in significant details (the global p values differ by less than .02, for example) and we also relegate plots for them to the ancillary information.

Conclusions
A Z model where the SM is augmented by an additional, spontaneously broken U (1) gauge group can simultaneously fit both the CDF M W anomaly and the b → s anomalies. The model retains the other desirable properties of the Y 3 model on which it is based [50]: namely that it qualitatively explains a hierarchically heavy third generation of quarks and small CKM angles. As Fig. 1 demonstrates, the log likelihood contribution from CDF's M W measurement is instrumental in picking out favoured values for the U (1) X quantum numbers of the fields. A particularly simple anomaly-free combination Y 3 − 3(B 3 − L 3 ) has a high quality of fit (where the 3 field is subsequently identified with the left-handed muon lepton  In the left-hand panel, we display the EWPOs, whereas in the right-hand panel, we display selected observables of interest to b → sl + l − anomalies doublet 11 ). We note that for the combination Y 3 −3(B 3 − L 3 ), the 'LFU' and 'EWPO' classes of observable (both of which have small theoretical uncertainties) each separately have a better quality of fit than the 'quarks' class, where there is more room for argument about the prediction and the size of the theoretical uncertainty assigned.  .13). In the left-hand panel, we display the EWPOs, whereas in the right-hand panel, we display selected observables of interest to b → sl + l − anomalies tations of the gauge group, for example. We argue, following Ref. [86], that it is premature to set the details of the model in the ultra-violet more than we have, preferring instead to allow the data (such as the new measurement of M W ) to inform the model building in a more fundamental and vital way. Indeed, we have used it in the present paper to select U (1) X charges under the symmetry group.
In the coming years, the LHC experiments will provide valuable further empirical measurements of M W and observables pertinent to the b → s anomalies and Belle II data will also weigh in. In the meantime, an obvious avenue of interest is from direct searches for the Z . The most promising channels [81] are likely to be pp → Z → μ + μ − with or without additional b-jets. LHC and HL-LHC sensitivity estimates are an obvious target for our models for future study. 13 Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theoretical study and no experimental data has been listed.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.