measurements are made using the combination of events
from the three channels, where the fraction of signal events
in the data sample, estimated from simulation, is 79%.
In Figs. 2 and 3, distributions of all the reconstructed
angular observables (defined in Sec. II) are shown. There is
reasonable agreement between the data and the sum of the
expected signal and background contributions given the
systematic uncertainties. In addition to the systematic
uncertainties discussed in Sec. VII, two uncertainties that
affect only the normalization of the measured differential
cross section are considered: the 2.5% uncertainty in the
integrated luminosity of the data sample [5] is applied to the
normalization of all simulated predictions, and a 1.5%
normalization uncertainty is applied to the t
¯
t prediction to
account for the uncertainty in the dileptonic branching
fraction (BF) [1]. The shapes of the reconstructed distri-
butions differ substantially from the expected parton-level
functional forms of Eqs. (8)–(10) owing to the effects of bin
migration, detector acceptance and efficiency, and back-
ground events, which are described in Sec. VI.
VI. UNFOLDING THE DIFFERENTIAL
CROSS SECTIONS
The effects of detector acceptance and efficiency sculpt
the reconstructed distributions, and the smearing intro-
duced by the detector response, kinematic reconstruction
algorithm, PS, and hadronization leads to the migration of
events across bins. In order to measure the differential cross
sections at the parton level in the full phase space, these
effects are accounted for by using the TUnfold regularized
unfolding method [51]. The response matrix used in the
unfolding is calculated for each measured distribution using
the default t
¯
t simulation, where the momenta of the parton-
level top quarks are defined after QCD radiation has been
simulated but before the top quark decays.
To keep the bin-to-bin migrations small (to avoid strong
bin-to-bin correlations in the unfolded distributions), the
widths of the measurement bins are chosen according to the
reconstruction resolution of the observable. This is quanti-
fied both directly by comparing the generator level and
detector level observables in simulation, and by measuring
the purity and stability. Purity is defined as the fraction of
events in a given bin at the detector level that originate from
the same bin at the generator level, and stability is defined
as the fraction of events in a given bin at the generator level
that are reconstructed in the same bin at the detector level.
For all observables measured in the top quark rest frame,
the use of six bins of uniform width is found to be well
matched to the reconstruction resolution. The purities and
stabilities are typically 40%. For the observables measured
in the laboratory frame (cos φ
lab
and jΔϕ
ll
j), six uniform-
width bins are also used. These observables have excellent
experimental resolution, and the purities and stabilities
are >99%.
The presence of background events is accounted for prior
to performing the unfolding. After subtracting all other
background components, the background from dileptonic t
¯
t
events with leptonically decaying τ leptons is subtracted as
a fraction of the total remaining events. The fraction is
evaluated per bin as the ratio of the background to the total
dileptonic t
¯
t events in simulation. Thus, the shapes of the
distributions for dileptonic t
¯
t events are taken from data,
and any dependence on the total cross section used in the
normalization of the simulated t
¯
t sample is avoided.
In TUnfold, a procedure based on matrix inversion is
used to obtain an unfolded distribution from the measured
distribution by applying a χ
2
minimization technique. The
potential large statistical fluctuations and strong anticorre-
lations between adjacent bins arising from the matrix
inversion are suppressed by introducing a term in the χ
2
expression that smooths (regularizes) the shape of the
unfolded distribution [51]. The regularization term penal-
izes the curvature of a vector constructed from the product
of the difference between the unfolded and simulated bin
values and a factor calculated using the expected functional
form [Eqs. (8)–(10)] such that a deviation in the coefficient
corresponds to a linear change in the vector. Since linear
changes are unconstrained by regularization of the curva-
ture, and the functional forms at the parton level (which are
unaffected by BSM phenomena in t
¯
t production) depend
only on the coefficient, this ensures that the regularization
cannot introduce a bias in the unfolded distribution. For the
laboratory-frame distributions there are no such simple
functional forms, and no factor is applied to the difference
vector. However, this choice is of little consequence
because the regularization is very weak owing to the low
level of bin migration.
The use of wide bins for the response matrix loses
information about its dependence inside each bin, meaning
the unfolding can be biased if the physical process density
differs from the simulation. Since the curvature regulari-
zation is unbiased, we make use of narrower bins in the
TUnfold χ
2
minimization; a factor 4 narrower is found to be
sufficient to reduce the bias from binning to a negligible
level. We have thus replaced the biased implicit regulari-
zation from binning with an unbiased regularization of the
curvature within each of the original bins.
The regularization level is determined for each distribu-
tion by minimizing the average global correlation coef-
ficient (ρ
avg
) [51], where ρ
avg
is determined after rebinning
to the original six bins.
For each measured bin, we perform tests using pseudo-
data to confirm a linear response of the method to variations
in the coefficient, and confirm that the distribution of the
difference between the nominal bin value and that mea-
sured in pseudodata, normalized to the measured uncer-
tainty, is consistent with having zero mean and unit width.
The data in the three channels are combined before
unfolding in order to model correlations between channels
MEASUREMENT OF THE TOP QUARK POLARIZATION AND t
¯
t … PHYS. REV. D 100, 072002 (2019)
072002-7