W-Mass Constraint Method in Electroweak Physics

• W-Mass Constraint Method is a precision technique that extracts the W boson mass from constrained kinematic observables in collider experiments.
• It utilizes in situ calibration, fast parametric simulation, and maximum-likelihood fits to achieve uncertainties as low as a few MeV.
• This method underpins precision electroweak tests, constraining Standard Model parameters and probing new physics through oblique corrections and extended sectors.

The W-Mass Constraint Method denotes a family of procedures in which the mass of the charged weak boson, $M_W$, is treated as a precision-defining quantity. In hadron-collider measurements, it refers to the extraction of $M_W$ from constrained kinematic observables in $W\to \ell\nu$ events, with in situ calibration, fast parametric simulation, and correlated statistical combination. In electroweak phenomenology, the same quantity is used as an input to radiative-correction relations or to oblique-parameter fits, thereby constraining Higgs-sector parameters, extra gauge bosons, scalar multiplets, or dark-sector spectra. At lepton colliders, a closely related usage denotes kinematic reconstruction with explicit four-momentum and equal-mass constraints in $W^+W^-$ events (Sa, 2012, Azzurri, 2021).

1. Experimental definition in hadron-collider measurements

At the Tevatron, the method begins with the selection of clean leptonic $W$ samples and the construction of observables with maximal sensitivity to $M_W$. The channels used are $W\to e\nu$ and $W\to \mu\nu$. DØ, with $4.3\,\mathrm{fb}^{-1}$, uses only the central-electron channel with $|\eta|<1.05$, while CDF, with $M_W$0, uses both central electrons and muons with $M_W$1. Typical kinematic requirements include lepton $M_W$2 or $M_W$3 cuts, missing transverse energy requirements, a transverse-mass window, and recoil suppression through $M_W$4 in order to avoid high-$M_W$5 $M_W$6 events. The final sample sizes are approximately $M_W$7 for DØ $M_W$8, and for CDF approximately $M_W$9 $W\to \ell\nu$0 plus $W\to \ell\nu$1 $W\to \ell\nu$2 (Sa, 2012).

The core kinematic quantities are the neutrino transverse momentum,

$W\to \ell\nu$3

the transverse mass,

$W\to \ell\nu$4

and the single-lepton distribution $W\to \ell\nu$5. CDF also uses the $W\to \ell\nu$6 distribution, whereas DØ drops it in the final analysis. The operational logic is to isolate those one-dimensional projections whose line shape near the Jacobian region is most sensitive to changes in the trial value of $W\to \ell\nu$7 (Sa, 2012).

This formulation makes the method intrinsically differential rather than purely counting-based. The measurement is not obtained from a direct invariant-mass reconstruction, since the neutrino longitudinal momentum is not observed, but from a constrained inference in transverse kinematics.

2. Calibration, fast simulation, and likelihood extraction

A defining feature of the method is the calibration of the lepton scale in situ to $W\to \ell\nu$8. DØ uses a calorimeter-based electron scale. It fits $W\to \ell\nu$9 events simultaneously in the invariant mass $W^+W^-$0 and in the auxiliary variable

$W^+W^-$1

with $W^+W^-$2 the opening angle. A two-dimensional binned likelihood in $W^+W^-$3 yields an energy scale $W^+W^-$4 and offset $W^+W^-$5 in each luminosity bin. Materials upstream, layer-by-layer response, underlying-event dependence, and luminosity-dependent gain are tuned by shower-shape and minimum-bias overlays. The final electron-scale uncertainty is $W^+W^-$6, corresponding to approximately $W^+W^-$7 on $W^+W^-$8, and the procedure effectively measures $W^+W^-$9, cancelling many systematics (Sa, 2012).

CDF uses a tracker-based momentum scale. The COT drift chamber is aligned and calibrated with cosmic-ray muons, and weak modes are removed through $W$0 $W$1 comparisons. The momentum scale is determined from fits to the $W$2, $W$3, and $W$4 invariant-mass peaks, with corrections for $W$5-field non-uniformities and ionization energy loss. The achieved momentum-scale precision is $W$6, corresponding to approximately $W$7 on $W$8. The scale is then transferred to electrons through $W$9 fits in $M_W$0 and $M_W$1 (Sa, 2012).

Template generation is performed with a fast Parametrized Monte Carlo Simulation. It generates $M_W$2, $M_W$3, and $M_W$4 distributions as functions of trial $M_W$5, including lepton resolution, radiative tails from QED radiation, a hadronic recoil model tuned on $M_W$6 events, and a parton-level $M_W$7 spectrum from resummed NLO QCD. Backgrounds include QCD multijet events faking leptons, $M_W$8, $M_W$9, and cosmic muons. The extraction is based on binned maximum-likelihood fits of the form

$W\to e\nu$0

DØ uses two observables, $W\to e\nu$1 and $W\to e\nu$2, while CDF uses six fits: $W\to e\nu$3, $W\to e\nu$4, and $W\to e\nu$5 in both the electron and muon channels. The summary of the method characterizes these as blinded maximum-likelihood fits built on fast parametric Monte Carlo templates (Sa, 2012).

Systematic uncertainties are evaluated by varying detector and production-model inputs in the fast simulation and refitting $W\to e\nu$6. Experimental systematics include lepton energy or momentum scale and resolution, hadronic recoil scale and resolution, lepton identification and trigger efficiencies versus $W\to e\nu$7, and background normalizations and shapes. Production-model systematics include PDFs, the $W\to e\nu$8 model, and higher-order QED radiation. The total systematic error is approximately $W\to e\nu$9 for DØ and $W\to \mu\nu$0 for CDF, with CDF also quoting a statistical uncertainty of approximately $W\to \mu\nu$1 and a total uncertainty of approximately $W\to \mu\nu$2 (Sa, 2012).

3. BLUE combination and indirect Higgs-mass inference

Once individual measurements are obtained, the method proceeds to a global combination with the BLUE estimator. If $W\to \mu\nu$3 are the individual measurements and $W\to \mu\nu$4 is the full covariance matrix, the combined value is

$W\to \mu\nu$5

with variance

$W\to \mu\nu$6

Experimental systematics such as detector effects and backgrounds are taken as uncorrelated between CDF and DØ, whereas production-model systematics such as PDFs, QED effects, and the $W\to \mu\nu$7 model are treated as partially correlated, with a common theory component taken as $W\to \mu\nu$8 correlated and the experiment-specific remainder uncorrelated. All measurements are corrected to a common $W\to \mu\nu$9 width $4.3\,\mathrm{fb}^{-1}$0 (Sa, 2012).

Using this procedure, the DØ combined value is reported as

$4.3\,\mathrm{fb}^{-1}$1

and the CDF value as

$4.3\,\mathrm{fb}^{-1}$2

The Tevatron combination over Run 0, I, and II yields

$4.3\,\mathrm{fb}^{-1}$3

with $4.3\,\mathrm{fb}^{-1}$4 and $4.3\,\mathrm{fb}^{-1}$5, and the world average obtained from Tevatron and LEP is

$4.3\,\mathrm{fb}^{-1}$6

These values define the precision benchmark against which later phenomenological studies formulate their constraints (Sa, 2012).

The same measured quantity enters the electroweak radiative-correction relation

$4.3\,\mathrm{fb}^{-1}$7

Here $4.3\,\mathrm{fb}^{-1}$8 encapsulates loop contributions and has logarithmic sensitivity to the Higgs mass $4.3\,\mathrm{fb}^{-1}$9. The global electroweak fit quoted in the Tevatron study gives

$|\eta|<1.05$0

before the new Tevatron $|\eta|<1.05$1 results and

$|\eta|<1.05$2

after adding the Tevatron Run II measurements. The updated fit strongly favors a light Higgs in the $|\eta|<1.05$3–$|\eta|<1.05$4 window and excludes high values above approximately $|\eta|<1.05$5 (Sa, 2012).

In this sense, the method is not limited to measuring $|\eta|<1.05$6; it uses $|\eta|<1.05$7 as a high-leverage electroweak input.

4. Kinematic-fit realization at lepton colliders

At FCC-ee, the phrase denotes a more explicit constraint fit in reconstructed $|\eta|<1.05$8 events. Measured objects are represented by four-vectors

$|\eta|<1.05$9

with covariance matrices $M_W$00, and one introduces fitted four-vectors

$M_W$01

that are adjusted within their uncertainties to satisfy physics constraints. The standard constraints are the 4C conditions of total energy equal to $M_W$02 and total three-momentum equal to zero, plus a fifth constraint in the 5C fit requiring the two reconstructed $M_W$03 bosons to have equal invariant mass (Azzurri, 2021).

A Lagrange-multiplier formulation of the fitted objective is

$M_W$04

For the fully hadronic channel, the equal-mass constraint is

$M_W$05

with $M_W$06 and $M_W$07. After the fit, the reconstructed mass is taken as

$M_W$08

In the $M_W$09 channel all three jet pairings are tried and the one with the smallest $M_W$10 is retained; in the $M_W$11 channel the neutrino is treated as an invisible object with measured $M_W$12 and large covariance (Azzurri, 2021).

The projected statistical precision of this kinematic-reconstruction method is similar to the threshold-scan method: approximately $M_W$13 for the $M_W$14 mass and approximately $M_W$15 for the width, using $M_W$16-pair data collected at threshold and at $M_W$17–$M_W$18. The threshold scan itself, with $M_W$19 shared on energy points between $M_W$20 and $M_W$21, is projected to yield a statistical uncertainty of $M_W$22 on the mass and $M_W$23 on the width. Uncertainty propagation is summarized by relations such as $M_W$24 and $M_W$25, while the jet-energy-scale response is approximated by $M_W$26 (Azzurri, 2021).

The dominant systematic issues differ from the Tevatron template method. Beam-energy calibration enters directly through the energy constraint, and hadronization and fragmentation modeling affect jet response and the fitted $M_W$27 factors. The control-sample strategy relies on $M_W$28 and $M_W$29 events reconstructed and fitted with the same techniques as the $M_W$30 events, together with $M_W$31–$M_W$32 $M_W$33 hadrons events at the $M_W$34 pole (Azzurri, 2021).

5. Precision-constraint use in extended electroweak sectors

In beyond-the-Standard-Model studies, the method shifts from direct extraction to indirect exclusion or accommodation. One class of realizations introduces tree-level $M_W$35–$M_W$36 mixing. In the $M_W$37 basis, the neutral-boson mass matrix is

$M_W$38

with eigenvalues

$M_W$39

and mixing angle

$M_W$40

The lighter eigenvalue is identified with the physical $M_W$41 mass, and solving for $M_W$42 in terms of $M_W$43 and $M_W$44 fixes $M_W$45 and hence $M_W$46. In the Derivative Portal Dark Matter model, fitting the shifted $M_W$47 mass together with $M_W$48, $M_W$49, and $M_W$50 yields a best fit at $M_W$51, $M_W$52, and $M_W$53, with all constraints overlapping in a broad region $M_W$54 and $M_W$55–$M_W$56. In the simple $M_W$57 extension, the best compromise gives only $M_W$58, with $M_W$59, $M_W$60, and $M_W$61 (Zeng et al., 2022).

A second class uses oblique corrections from extended scalar sectors. In the 2HDM+$M_W$62 model, the leading precision effect is encoded in the Peskin–Takeuchi parameters $M_W$63, $M_W$64, and $M_W$65, with

$M_W$66

and

$M_W$67

ScannerS is used to compute the one-loop $M_W$68, $M_W$69, and $M_W$70 numerically from scalar loops. The fit is performed through

$M_W$71

with benchmarks taken from CDF II, ATLAS, and the world average. The study requires $M_W$72 relative to the world average and separately relative to the ATLAS result, and finds that no scalar-only point reaches the high CDF value within $M_W$73, whereas many lie in the ATLAS and world-average bands. The surviving points populate a region with $M_W$74 and $M_W$75 (Mulaudzi et al., 2023).

A third class introduces high-dimensional scalar multiplets. For real $M_W$76 multiplets of odd dimension $M_W$77 and hypercharge $M_W$78, one has

$M_W$79

so that the $M_W$80 mass is unchanged by the new vacuum expectation values. For a single real septuplet, the required value is $M_W$81–$M_W$82. The same paper studies the one-loop mechanism from a complex scalar octuplet with $M_W$83, expressing the shift through $M_W$84 and $M_W$85; the viable Type A region is approximately $M_W$86 with $M_W$87–$M_W$88 (Wu et al., 2023).

Across these examples, the method has a common structure: a model modifies either the tree-level relation or the self-energies entering $M_W$89, $M_W$90, $M_W$91, and the measured or benchmarked $M_W$92 is then converted into a sharply delimited region of parameter space.

6. Dark-sector and neutrino-mass implementations

The same logic is applied in models where the $M_W$93-mass shift is tied to dark matter, neutrino mass, or both. In the singlet-triplet scotogenic model, a hyperchargeless real $M_W$94 triplet scalar $M_W$95 with vacuum expectation value $M_W$96 modifies the gauge-boson masses as

$M_W$97

For small shifts,

$M_W$98

Matching the difference between the CDF-II central value and the Standard Model prediction gives the “naïve” tree-level bound

$M_W$99

The model also includes one-loop contributions to $W\to \ell\nu$00, $W\to \ell\nu$01, and $W\to \ell\nu$02 from the $W\to \ell\nu$03-odd scalar doublet $W\to \ell\nu$04, and a covariance-matrix fit using

$W\to \ell\nu$05

selects the viable parameter region. In the loop-dominated regime $W\to \ell\nu$06, the allowed region in the $W\to \ell\nu$07 plane is a narrow diagonal band, and the maximal splitting allowed at $W\to \ell\nu$08 is approximately $W\to \ell\nu$09 in the $W\to \ell\nu$10 fit (Batra et al., 2022).

In the singlet-doublet Majorana-fermion model, the relevant quantity is the new-physics contribution to $W\to \ell\nu$11,

$W\to \ell\nu$12

which is translated into the $W\to \ell\nu$13-mass shift through the on-shell relation, with leading-order approximation

$W\to \ell\nu$14

To reproduce the CDF central value, the required oblique correction is

$W\to \ell\nu$15

and the $W\to \ell\nu$16 band is

$W\to \ell\nu$17

A single generation of singlet-doublet fermions does not leave overlap between the $W\to \ell\nu$18-preferred region and the dark-matter-allowed region, whereas two or three generations can do so, with the heavier generation or generations driving the $W\to \ell\nu$19-mass correction and the lighter generation accounting for the dark-matter phenomenology (Borah et al., 2022).

These implementations illustrate that the method can be either tree-level or loop-level, but its statistical endpoint is the same: a narrow numerical band in a space of masses, splittings, mixing angles, or vacuum expectation values.

7. Terminological extensions outside electroweak phenomenology

The expression also appears in mathematical contexts unrelated to the $W\to \ell\nu$20 boson. In constrained optimal transport, the “W–Mass Constraint Method” refers to a constrained Wasserstein problem built on the Benamou–Brenier dynamic formulation,

$W\to \ell\nu$21

supplemented by a hard or soft mass-control constraint

$W\to \ell\nu$22

The resulting saddle-point problem is solved by a primal-dual proximal splitting scheme, and a convergence theorem is proved under the step-size condition

$W\to \ell\nu$23

The paper explicitly states that the name reflects the computation of a constrained Wasserstein geodesic under an additional mass or flux constraint (Kerrache et al., 2022).

In the Cahn–Hilliard equation with dynamic boundary conditions, the same wording designates a prescribed weighted boundary-mass constraint rather than a gauge-boson observable. The system imposes

$W\to \ell\nu$24

and introduces two Lagrange multipliers, $W\to \ell\nu$25 and $W\to \ell\nu$26, yielding

$W\to \ell\nu$27

in the bulk and

$W\to \ell\nu$28

on the boundary. Well-posedness is established in a subdifferential-evolution framework, with existence, uniqueness, and characterization of the two multipliers (Colli et al., 2014).

These usages are mathematically unrelated to electroweak precision physics, but they show that the phrase “W-Mass Constraint Method” is not unique to collider phenomenology. In current high-energy usage, however, the dominant meaning is the constrained extraction or precision use of the $W\to \ell\nu$29-boson mass as a discriminator of Standard Model consistency and of new-physics parameter space.