research30 min studyLesson 620 in Physics

Electroweak fit inputs and correlated tests

Construct a reproducible precision comparison with declared inputs, correlated residuals, nuisance parameters, and an honest separation of data from predictions.

On this page

What you’ll understand

  • Separate independent input observables from genuine predicted tests.
  • Calculate a correlated Gaussian chi-squared and interpret its eigen-directions.
  • Profile a simple nuisance parameter and state the statistical assumptions of the comparison.

A successful electroweak fit compares several measurements with predictions from a shared parameter set. Its evidential strength comes from redundant constraints, not from treating the same input several times as independent data. A reported tension also depends on experimental correlations, theoretical uncertainties, mass and coupling definitions, and which parameters were fitted. We build a transparent Gaussian classroom model to expose these issues. Every measurement, error and correlation in the numerical exercises is invented; none represents a current global electroweak fit.

Choose inputs before declaring a prediction

At leading order a conventional electroweak input set contains the electromagnetic coupling, the Fermi constant and the Z mass. These determine combinations of the gauge couplings and vacuum scale, allowing a W-mass prediction. At higher order additional parameters such as top and Higgs masses, strong coupling and hadronic vacuum polarization enter. The chosen renormalization and unstable-particle mass conventions determine how those measured inputs map to Lagrangian parameters.

If the measured W mass is instead included as an input, reproducing that same number cannot be counted as an independent successful W-mass test. One can still test other observables or examine an indirect prediction obtained by excluding that datum. A leave-one-out result and a full-fit residual answer related but different questions. Experimental combinations also need their shared calibration and modeling systematics retained, rather than replacing several correlated analyses with an artificially small independent error.

ri=Oiobs−Oith(ϑ),χ2(ϑ)=rTC−1r,Cij=Cov⁡(Oi,Oj)r_i=O_i^{\rm obs}-O_i^{\rm th}(\vartheta),\qquad\chi^2(\vartheta)=r^TC^{-1}r,\qquad C_{ij}=\operatorname{Cov}(O_i,O_j)
The vector ϑ contains fitted theory parameters. Covariance entries carry the product of their observable units, making χ² dimensionless. The simple Gaussian form assumes a suitable likelihood approximation and a consistent positive-definite covariance.↔ Scroll wide equations horizontally.

Correlations change the meaningful residual directions

For two observables, divide each residual by its own standard deviation and denote their correlation by c. Inverting the two-by-two correlation matrix yields a quadratic form with a cross term. Equal-sign standardized residuals lie along a common direction whose variance grows for positive correlation. Opposite-sign residuals lie along a difference direction whose variance shrinks. Correlation can therefore reduce one discrepancy and strengthen another; ignoring it has no uniformly conservative effect.

A covariance built from experimental errors is not automatically a complete theory-error model. Common missing corrections may be correlated across observables, and some theory uncertainties are better represented by explicit nuisance parameters or non-Gaussian likelihoods. If input uncertainties have already been propagated into a covariance, also fitting the same fluctuations independently can double count them. The fit specification should make every uncertainty source and dependence reproducible rather than hiding it inside a single quoted total error.

χ2=z12+z22−2cz1z21−c2,zi=riσi,−1<c<1\chi^2=\frac{z_1^2+z_2^2-2cz_1z_2}{1-c^2},\qquad z_i=\frac{r_i}{\sigma_i},\qquad -1<c<1
The correlation c is distinct from the weak-angle cosine cW. This expression assumes known Gaussian covariance and two observables with positive standard deviations; the singular limits require a reduced description.↔ Scroll wide equations horizontally.

Fit parameters and nuisance assumptions affect the inference

A nuisance parameter can represent a shared calibration shift, a hadronic input, or a theoretical normalization. Profiling minimizes the objective over its allowed values, including any external constraint term; marginalization instead integrates over it with specified probability distributions. These procedures need not agree beyond simple Gaussian settings. A change in chi-squared should identify which parameters were reoptimized and which external constraints were retained.

For a regular Gaussian model with known covariance and a correctly specified linear mean, a minimized chi-squared may be compared with a distribution whose degrees of freedom subtract the independently fitted parameters. Boundaries, nonlinearities, estimated covariance and model selection can change that calibration. A p-value measures compatibility with a specified test construction, not the probability that the Standard Model is true. Precision information is most useful when a published input list and likelihood allow other researchers to reproduce the claimed constraints and test alternative explanations.

Two correlated toy residuals

A toy theory predicts mW=80.0 GeV and an effective angle-squared observable 0.230. Invented observations are 80.1±0.1 GeV and 0.231±0.001, with correlation c=0.8. Treat the theory inputs as fixed. Compute χ² with and without the correlation.

  1. The standardized residuals are z₁=(80.1−80.0)/0.1=1 and z₂=(0.231−0.230)/0.001=1. The original observables have different units, which disappear after standardization.

  2. With c=0.8, the numerator is 1+1−2×0.8=0.4 and denominator 1−0.8²=0.36. Thus χ²=0.4/0.36=10/9≈1.111111.

  3. Ignoring correlation gives χ²=1²+1²=2. The positive correlation makes a shared same-sign displacement less surprising; it does not erase either individual standardized residual.

The result

The fixed-input toy comparison gives correlated χ²=10/9, compared with 2 under an independent-error assumption. Neither calculation is an actual electroweak fit.

Explore the model

Interactive lab

Interactive lab · sampling

Sampling spread depends on sample size

For each fresh comparison, Reset and set both controls. Reset and set success probability p=0.5 and independent sample count n=25. Predict sample-mean variance 0.01. Compare a fresh run with both controls p=0.5 and n=100, giving variance 0.0025 and half the standard deviation. These are exact independent Bernoulli sampling distributions, not resampled measurements. The lab has no correlated electroweak likelihood, nuisance parameters, global fit, or experimental precision prediction.

Sampling spread depends on sample sizeExact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Current readings: Mean sample proportion: 0.5; Sample-mean variance: 0.0125; Uniform CDF error bound: ≤ 0.1062. All parameters are adjustable with the labeled controls below.Exact distribution of the mean of 20 independent Bernoulli variables with success probability 0.5; playback highlights outcome 0.00.250.50.75100.0580.1150.1730.23Sample mean · k / nProbability massNormal bin masses are an approximationk = 0 · exact probability = 9.54e-7
Exact binomial massNormal interval massHighlighted outcome / mean
Mean sample proportion
0.5
Sample-mean variance
0.0125
Uniform CDF error bound
≤ 0.1062
t = 0 s

Exact binomial masses for the mean of n independent Bernoulli trials, compared with normal mass between adjacent half-count boundaries; normal tails outside the finite support remain outside the bars. Playback highlights outcomes, not random samples. The Berry–Esseen bound applies to standardized CDFs, not individual bars or continuity-corrected errors. At p = 0 or 1, normal comparison and that bound are disabled. Normal-CDF evaluation uses a numerical approximation with absolute error below 7.5×10⁻⁸ (at most 1.5×10⁻⁷ per interval, before rounding). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forIndependent statistical error can shrink with sample size; shared systematic or theory correlations require a separate model.

Check your understanding

A measured observable was used to fix an input parameter. What is required before claiming an independent prediction for that same observable?

Try it yourself

Two invented dimensionless measurements are 2 and 4, each with independent Gaussian standard deviation 1. Both measure the same unknown parameter θ plus a shared nuisance shift b constrained by b=0±2. Profile the objective (2−θ−b)²+(4−θ−b)²+b²/4.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

First minimize over the common predicted value θ+b, then use the external nuisance constraint.

Compare with the solution

For a fixed b, the best common prediction is θ+b=3, leaving residuals −1 and +1 and a data contribution 2. The remaining objective is 2+b²/4, minimized at b=0 and θ=3, with χ²min=2. The free θ absorbs a common shift, so the nuisance prior does not improve the mismatch between the two measurements. In this linear Gaussian model the uncertainty of θ is √(1/2+4)=√4.5≈2.12132; shared uncertainty affects the inferred parameter even though the residual difference remains fixed.

Take it with you

  • A genuine precision test distinguishes externally supplied inputs from independently predicted observables.
  • Correlations define common and difference residual directions with different uncertainties.
  • Nuisance handling and likelihood assumptions belong to the result, alongside the central values and fitted parameters.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.