The inclusive transverse momentum ($p_{\rm T}$) distributions of primary charged particles are measured in the pseudo-rapidity range $|\eta|<0.8$ as a function of event centrality in Pb-Pb collisions at $\sqrt{s_{\rm{NN}}}=2.76$ TeV with ALICE at the LHC. The data are presented in the $p_{\rm T}$ range $0.15<p_{\rm T}<50$ GeV/$c$ for nine centrality intervals from 70-80% to 0-5%. The Pb-Pb spectra are presented in terms of the nuclear modification factor $R_{\rm{AA}}$ using a pp reference spectrum measured at the same collision energy. We observe that the suppression of high-$p_{\rm T}$ particles strongly depends on event centrality. In central collisions (0-5%) the yield is most suppressed with $R_{\rm{AA}}\approx0.13$ at $p_{\rm T}=6$-7 GeV/$c$. Above $p_{\rm T}=7$ GeV/$c$, there is a significant rise in the nuclear modification factor, which reaches $R_{\rm{AA}} \approx0.4$ for $p_{\rm T}>30$ GeV/$c$. In peripheral collisions (70-80%), the suppression is weaker with $R_{\rm{AA}} \approx 0.7$ almost independently of $p_{\rm T}$. The measured nuclear modification factors are compared to other measurements and model calculations.

We present measurements of Underlying Event observables in pp collisions at $\sqrt{s}$ = 0.9 and 7 TeV. The analysis is performed as a function of the highest charged-particle transverse momentum $p_{\rm T, LT}$ in the event. Different regions are defined with respect to the azimuthal direction of the leading (highest transverse momentum) track: Toward, Transverse and Away. The Toward and Away regions collect the fragmentation products of the hardest partonic interaction. The Transverse region is expected to be most sensitive to the Underlying Event activity. The study is performed with charged particles above three different $p_{\rm T}$ thresholds: 0.15, 0.5 and 1.0 GeV/$c$. In the Transverse region we observe an increase in the multiplicity of a factor 2-3 between the lower and higher collision energies, depending on the track $p){\rm T}$ threshold considered. Data are compared to Pythia 6.4, Pythia 8.1 and Phojet. On average, all models considered underestimate the multiplicity and summed $p_{\rm T}$ in the Transverse region by about 10-30%.

We present the measurements of particle pair yields per trigger particle obtained from di-hadron azimuthal correlations in pp collisions at $\sqrt{s} = 0.9$, $2.76$, and $7$ TeV recorded with the ALICE detector. The yields are studied as a function of the charged particle multiplicity. Taken together with the single particle yields the pair yields provide information about parton fragmentation at low transverse momenta, as well as on the contribution of multiple parton interactions to particle production. Data are compared to calculations using the PYTHIA6, PYTHIA8, and PHOJET event generators.

In this paper measurements are presented of $\rm \pi$$^+$, $\rm \pi$$^-$, K$^+$, K$^-$, p and $\overline{\rm p}$ production at mid-rapidity < 0.5, in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV as a function of centrality. The measurement covers the transverse momentum ($p_{\rm T}$) range from 100, 200, 300 MeV/$c$ up to 3, 3, 4.6 GeV/$c$, for $\rm\pi$, K, and p respectively. The measured $p_{\rm T}$ distributions and yields are compared to expectations based on hydrodynamic, thermal and recombination models. The spectral shapes of central collisions show a stronger radial flow than measured at lower energies, which can be described in hydrodynamic models. In peripheral collisions, the $p_{\rm T}$ distributions are not well reproduced by hydrodynamic models. Ratios of integrated particle yields are found to be nearly independent of centrality. The yield of protons normalized to pions is a factor ~1.5 lower than the expectation from thermal models.

The yields of the K*(892)$^{0}$ and $\Phi$(1020) resonances are measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV through their hadronic decays using the ALICE detector. The measurements are performed in multiple centrality intervals at mid-rapidity (|$y$|<0.5) in the transverse-momentum ranges 0.3 < $p_{\rm T}$ < 5 GeV/$c$ for the K*(892)$^{0}$ and 0.5 < $p_{\rm T}$ < 5 GeV/$c$ for the $\Phi$(1020). The yields of K*(892)$^{0}$ are suppressed in central Pb-Pb collisions with respect to pp and peripheral Pb-Pb collisions (perhaps due to rescattering of its decay products in the hadronic medium), while the longer lived $\Phi$(1020) meson is not suppressed. These particles are also used as probes to study the mechanisms of particle production. The shape of the $p_{\rm T}$ distribution of the $\Phi$(1020) meson, but not its yield, is reproduced fairly well by hydrodynamic models for central Pb-Pb collisions. In central Pb-Pb collisions at low and intermediate $p_{\rm T}$, the p/$\Phi$(1020) ratio is flat in $p_{\rm T}$, while the p/$\pi$ and $\Phi$(1020)/$\pi$ ratios show a pronounced increase and have similar shapes to each other. These results indicate that the shapes of the $p_{\rm T}$ distributions of these particles in central Pb-Pb collisions are determined predominantly by the particle masses and radial flow. Finally, $\Phi$(1020) production in Pb-Pb collisions is enhanced, with respect to the yield in pp collisions and the yield of charged pions, by an amount similar to the $\Lambda$ and $\Xi$.

The ALICE measurement of K$^0_{\rm S}$ and $\rm\Lambda$ production at mid-rapidity in Pb-Pb collisions at $\sqrt{s_{\rm NN}} = 2.76$ TeV is presented. The transverse momentum ($p_{\rm T}$) spectra are shown for several collision centrality intervals and in the $p_{\rm T}$ range from 0.4 GeV/$c$ (0.6 GeV/$c$ for $\rm\Lambda$) to 12 GeV/$c$. The $p_{\rm T}$ dependence of the $\rm \Lambda$/K$^0_{\rm S}$ ratios exhibits maxima in the vicinity of 3 GeV/$c$, and the positions of the maxima shift towards higher $p_{\rm T}$ with increasing collision centrality. The magnitude of these maxima increases by almost a factor of three between most peripheral and most central Pb-Pb collisions. This baryon excess at intermediate $p_{\rm T}$ is not observed in pp interactions at sqrt(s) = 0.9 TeV and at sqrt(s) = 7 TeV. Qualitatively, the baryon enhancement in heavy-ion collisions is expected from radial flow. However, the measured $p_{\rm T}$ spectra above 2 GeV/$c$ progressively decouple from hydrodynamical-model calculations. For higher values of $p_{\rm T}$, models that incorporate the influence of the medium on the fragmentation and hadronization processes describe qualitatively the $p_{\rm T}$ dependence of the $\rm\Lambda$/K$^0_{\rm S}$ ratio.

The average transverse momentum $\langle p_{\rm T}\rangle$ versus the charged-particle multiplicity $N_{\rm ch}$ was measured in p-Pb collisions at a collision energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}=5.02$ TeV and in pp collisions at collision energies of $\sqrt{s}=0.9$, 2.76, and 7 TeV in the kinematic range $0.15<p_{\rm T}<10.0$ GeV/$c$ and $|\eta|<0.3$ with the ALICE apparatus at the LHC. These data are compared to results in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV at similar charged-particle multiplicities. In pp and p-Pb collisions, a strong increase of $\langle p_{\rm T}\rangle$ with $N_{\rm ch}$ is observed, which is much stronger than that measured in Pb-Pb collisions. For pp collisions, this could be attributed, within a model of hadronizing strings, to multiple-parton interactions and to a final-state color reconnection mechanism. The data in p-Pb and Pb-Pb collisions cannot be described by an incoherent superposition of nucleon-nucleon collisions and pose a challenge to most of the event generators.

In high-energy collisions, the spatio-temporal size of the particle production region can be measured using the Bose-Einstein correlations of identical bosons at low relative momentum. The source radii are typically extracted using two-pion correlations, and characterize the system at the last stage of interaction, called kinetic freeze-out. In low-multiplicity collisions, unlike in high-multiplicity collisions, two-pion correlations are substantially altered by background correlations, e.g. mini-jets. Such correlations can be suppressed using three-pion cumulant correlations. We present the first measurements of the size of the system at freeze-out extracted from three-pion cumulant correlations in pp, p-Pb and Pb-Pb collisions at the LHC with ALICE. At similar multiplicity, the invariant radii extracted in p-Pb collisions are found to be 5-15% larger than those in pp, while those in Pb-Pb are 35-55% larger than those in p-Pb. Our measurements disfavor models which incorporate substantially stronger collective expansion in p-Pb as compared to pp collisions at similar multiplicity.

The differential charged jet cross sections, jet fragmentation distributions, and jet shapes are measured in minimum bias proton-proton collisions at centre-of-mass energy $\sqrt{s}=7$ TeV using the ALICE detector at the LHC. Jets are reconstructed from charged particle momenta in the mid-rapidity region using the sequential recombination $k_{\rm T}$ and anti-$k_{\rm T}$ as well as the SISCone jet finding algorithms with several resolution parameters in the range $R=0.2$ to $0.6$. Differential jet production cross sections measured with the three jet finders are in agreement in the transverse momentum ($p_{\rm T}$) interval $20<p_{\rm T}^{\rm jet,ch}<100$ GeV/$c$. They are also consistent with prior measurements carried out at the LHC by the ATLAS collaboration. The jet charged particle multiplicity rises monotonically with increasing jet $p_{\rm T}$, in qualitative agreement with prior observations at lower energies. The transverse profiles of leading jets are investigated using radial momentum density distributions as well as distributions of the average radius containing 80% ($\langle R_{\rm 80} \rangle$) of the reconstructed jet $p_{\rm T}$. The fragmentation of leading jets with $R=0.4$ using scaled $p_{\rm T}$ spectra of the jet constituents is studied. The measurements are compared to model calculations from event generators (PYTHIA, PHOJET, HERWIG). The measured radial density distributions and $\langle R_{\rm 80} \rangle$ distributions are well described by the PYTHIA model (tune Perugia-2011). The fragmentation distributions are better described by HERWIG.

In this Letter, comprehensive results on ${\rm\pi}^\pm$, K$^\pm$, K$^0_S$, p, $\rm\bar{p}$, $\rm \Lambda$ and $\rm \bar{\Lambda}$ production at mid-rapidity ($0 < y_{\rm cms} < 0.5$) in p-Pb collisions at $\sqrt{s_{\rm NN}} = 5.02$ TeV, measured by the ALICE detector at the LHC, are reported. The transverse momentum distributions exhibit a hardening as a function of event multiplicity, which is stronger for heavier particles. This behavior is similar to what has been observed in pp and Pb-Pb collisions at the LHC. The measured $p_{\rm T}$ distributions are compared to results at lower energy and with predictions based on QCD-inspired and hydrodynamic models.

A measurement of the transverse momentum spectra of jets in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV is reported. Jets are reconstructed from charged particles using the anti-$k_{\rm T}$ jet algorithm with jet resolution parameters $R$ of $0.2$ and $0.3$ in pseudo-rapidity $|\eta|<0.5$. The transverse momentum $p_{\rm T}$ of charged particles is measured down to $0.15$ GeV/$c$ which gives access to the low $p_{\rm T}$ fragments of the jet. Jets found in heavy-ion collisions are corrected event-by-event for average background density and on an inclusive basis (via unfolding) for residual background fluctuations and detector effects. A strong suppression of jet production in central events with respect to peripheral events is observed. The suppression is found to be similar to the suppression of charged hadrons, which suggests that substantial energy is radiated at angles larger than the jet resolution parameter $R=0.3$ considered in the analysis. The fragmentation bias introduced by selecting jets with a high $p_{\rm T}$ leading particle, which rejects jets with a soft fragmentation pattern, has a similar effect on the jet yield for central and peripheral events. The ratio of jet spectra with $R=0.2$ and $R=0.3$ is found to be similar in Pb-Pb and simulated PYTHIA pp events, indicating no strong broadening of the radial jet structure in the reconstructed jets with $R<0.3$.

We report the measurement of charged $D^*$ mesons in inclusive jets produced in proton-proton collisions at a center of mass energy $\sqrt{s}$ = 200 GeV with the STAR experiment at RHIC. For $D^{*}$ mesons with fractional momenta $0.2 < z < 0.5$ in inclusive jets with 11.5 GeV mean transverse energy, the production rate is found to be $N(D^{*+}+D^{*-})/N(\mathrm{jet}) = 0.015 \pm 0.008 (\mathrm{stat}) \pm 0.007 (\mathrm{sys})$. This rate is consistent with perturbative QCD evaluation of gluon splitting into a pair of charm quarks and subsequent hadronization.

The first measurements of the invariant differential cross sections of inclusive $\pi^0$ and $\eta$ meson production at mid-rapidity in proton-proton collisions at $\sqrt{s}=0.9$ TeV and $\sqrt{s}=7$ TeV are reported. The $\pi^0$ measurement covers the ranges $0.4<p_T<7$ GeV/$c$ and $0.3<p_T<25$ GeV/$c$ for these two energies, respectively. The production of $\eta$ mesons was measured at $\sqrt{s}=7$ TeV in the range $0.4<p_T<15$ GeV/$c$. Next-to-Leading Order perturbative QCD calculations, which are consistent with the $\pi^0$ spectrum at $\sqrt{s}=0.9$ TeV, overestimate those of $\pi^0$ and $\eta$ mesons at $\sqrt{s}=7$ TeV, but agree with the measured $\eta/\pi^0$ ratio at $\sqrt{s}=7$ TeV.

Invariant mass spectra for jets reconstructed using the anti-kt and Cambridge-Aachen algorithms are studied for different jet "grooming" techniques in data corresponding to an integrated luminosity of 5 inverse femtobarns, recorded with the CMS detector in proton-proton collisions at the LHC at a center-of-mass energy of 7 TeV. Leading-order QCD predictions for inclusive dijet and W/Z+jet production combined with parton-shower Monte Carlo models are found to agree overall with the data, and the agreement improves with the implementation of jet grooming methods used to distinguish merged jets of large transverse momentum from softer QCD gluon radiation.

Measurements of the differential cross sections for the production of exactly four jets in proton-proton collisions are presented as a function of the transverse momentum pt and pseudorapidity eta, together with the correlations in azimuthal angle and the pt balance among the jets. The data sample was collected in 2010 at a center-of-mass energy of 7 TeV with the CMS detector at the LHC, with an integrated luminosity of 36 inverse picobarns. The cross section for a final state with a pair of hard jets with pt > 50 GeV and another pair with pt > 20 GeV within abs(eta) < 4.7 is measured to be sigma = 330 +- 5 (stat.) +- 45 (syst.) nb. It is found that fixed-order matrix element calculations including parton showers describe the measured differential cross sections in some regions of phase space only, and that adding contributions from double parton scattering brings the Monte Carlo predictions closer to the data.

The underlying event activity in proton-proton collisions at forward pseudorapidity (-6.6 < eta < -5.2) is studied with the CMS detector at the LHC, using a novel observable: the ratio of the forward energy density, dE/d(eta), for events with a charged-particle jet produced at central pseudorapidity (abs(eta[jet]) < 2) to the forward energy density for inclusive events. This forward energy density ratio is measured as a function of the central jet transverse momentum, pt, at three different pp centre-of-mass energies (sqrt(s) = 0.9, 2.76, and 7 TeV). In addition, the sqrt(s) evolution of the forward energy density is studied in inclusive events and in events with a central jet. The results are compared to those of Monte Carlo event generators for pp collisions and are discussed in terms of the underlying event. Whereas the dependence of the forward energy density ratio on jet pt at each sqrt(s) separately can be well reproduced by some models, all models fail to simultaneously describe the increase of the forward energy density with sqrt(s) in both inclusive events and in events with a central jet.

A search for physics beyond the standard model in events with at least three leptons is presented. The data sample, corresponding to an integrated luminosity of 19.5 inverse femtobarns of proton-proton collisions with center-of-mass energy sqrt(s) = 8 TeV, was collected by the CMS experiment at the LHC during 2012. The data are divided into exclusive categories based on the number of leptons and their flavor, the presence or absence of an opposite-sign, same-flavor lepton pair (OSSF), the invariant mass of the OSSF pair, the presence or absence of a tagged bottom-quark jet, the number of identified hadronically decaying tau leptons, and the magnitude of the missing transverse energy and of the scalar sum of jet transverse momenta. The numbers of observed events are found to be consistent with the expected numbers from standard model processes, and limits are placed on new-physics scenarios that yield multilepton final states. In particular, scenarios that predict Higgs boson production in the context of supersymmetric decay chains are examined. We also place a 95% confidence level upper limit of 1.3% on the branching fraction for the decay of a top quark to a charm quark and a Higgs boson (t to c H), which translates to a bound on the left- and right-handed top-charm flavor-violating Higgs Yukawa couplings, lambda[H, tc] and lambda[H, ct], respectively, of sqrt(abs(lambda[H, tc])^2 + abs(lambda[H, ct])^2) < 0.21.

The polarizations of prompt J/psi and psi(2S) mesons are measured in proton-proton collisions at sqrt(s) = 7 TeV, using a dimuon data sample collected by the CMS experiment at the LHC, corresponding to an integrated luminosity of 4.9 inverse femtobarns. The prompt J/psi and psi(2S) polarization parameters lambda[theta], lambda[phi], and lambda[theta,phi], as well as the frame-invariant quantity lambda(tilde), are measured from the dimuon decay angular distributions in three different polarization frames. The J/psi results are obtained in the transverse momentum range 14 < pt < 70 GeV, in the rapidity intervals abs(y) < 0.6 and 0.6 < abs(y) < 1.2. The corresponding psi(2S) results cover 14 < pt < 50 GeV and include a third rapidity bin, 1.2 < abs(y) < 1.5. No evidence of large transverse or longitudinal polarizations is seen in these kinematic regions, which extend much beyond those previously explored.

A measurement is presented of the ratio of the inclusive 3-jet cross section to the inclusive 2-jet cross section as a function of the average transverse momentum, <pT[1,2]>, of the two leading jets in the event. The data sample was collected during 2011 at a proton-proton centre-of-mass energy of 7 TeV with the CMS detector at the LHC, corresponding to an integrated luminosity of 5.0 inverse femtobarns. The strong coupling constant at the scale of the Z boson mass is determined to be alphaS[MZ] = 0.1148 +/- 0.0014 (exp.) +/- 0.0018 (PDF) +/- 0.0050 (theory), by comparing the ratio in the range 0.42 < <pT[1,2]> < 1.39 TeV to the predictions of perturbative QCD at next-to-leading order. This is the first determination of alphaS[MZ] from measurements at momentum scales beyond 0.6 TeV. The predicted ratio depends only indirectly on the evolution of the parton distribution functions of the proton such that this measurement also serves as a test of the evolution of the strong coupling constant. No deviation from the expected behaviour is observed.

The polarizations of the Y(1S), Y(2S), and Y(3S) mesons are measured in proton-proton collisions at sqrt(s) = 7 TeV, using a data sample of Y(nS) to oppositely charged muon pair decays collected by the CMS experiment, corresponding to an integrated luminosity of 4.9 inverse femtobarns. The dimuon decay angular distributions are analyzed in three different polarization frames. The polarization parameters lambda[theta], lambda[phi], and lambda[theta,phi], as well as the frame-invariant quantity lambda-tilde, are presented as a function of the Y(nS) transverse momentum between 10 and 50 GeV, in the rapidity ranges abs(y) < 0.6 and 0.6 < abs(y) < 1.2. No evidence of large transverse or longitudinal polarizations has been seen in the explored kinematic region.

The $\Upsilon$(1S), $\Upsilon$(2S), and $\Upsilon$(3S) production cross sections are measured using a data sample corresponding to an integrated luminosity of 35.8 $\pm$ 1.4 inverse picobarns of proton-proton collisions at $\sqrt{s}$ = 7 TeV, collected with the CMS detector at the LHC. The Upsilon resonances are identified through their decays to dimuons. Integrated over the $\Upsilon$ transverse momentum range $p_{t}^{\Upsilon} \lt$ 50GeV and rapidity range |$y^\Upsilon$| $\lt$ 2.4, and assuming unpolarized Upsilon production, the products of the Upsilon production cross sections and dimuon branching fractions are \begin{equation*}\sigma(pp \to \Upsilon(1S) X) . B(\Upsilon(1S) \to \mu^+ \mu^-) = (8.55 \pm 0.05^{+0.56}_{-0.50} \pm 0.34) nb,\end{equation*} \begin{equation*}\sigma(pp \to \Upsilon(2S) X) . B(\Upsilon(2S) \to \mu^+ \mu^-) = (2.21 \pm 0.03^{+0.16}_{-0.14} \pm 0.09) nb,\end{equation*} \begin{equation*}\sigma(pp \to \Upsilon(3S) X) . B(\Upsilon(3S) \to \mu^+ \mu^-) = (1.11 \pm 0.02^{+0.10}_{-0.08} \pm 0.04) nb, \end{equation*} where the first uncertainty is statistical, the second is systematic, and the third is from the uncertainty in the integrated luminosity. The differential cross sections in bins of transverse momentum and rapidity, and the cross section ratios are presented. Cross section measurements performed within a restricted muon kinematic range and not corrected for acceptance are also provided. These latter measurements are independent of Upsilon polarization assumptions. The results are compared to theoretical predictions and previous measurements.

The charge asymmetry in t t-bar events is measured using dilepton final states produced in pp collisions at the LHC at sqrt(s) = 8 TeV. The data sample, collected with the CMS detector, corresponds to an integrated luminosity of 19.5 inverse femtobarns. The measurements are performed using events with two oppositely charged leptons (electrons or muons) and two or more jets, where at least one of the jets is identified as originating from a bottom quark. The charge asymmetry is measured from differences in kinematic distributions, unfolded to the parton level, of positively and negatively charged top quarks and leptons. The t t-bar and leptonic charge asymmetries are found to be 0.011 +/- 0.011 (stat) +/- 0.007 (syst) and 0.003 +/- 0.006 (stat) +/- 0.003 (syst), respectively. These results, as well as charge asymmetry measurements made as a function of the invariant mass, rapidity, and transverse momentum of the t t-bar system, are in agreement with predictions of the standard model.

Measurements of the top quark-antiquark (t t-bar) spin correlations and the top quark polarization are presented for t t-bar pairs produced in pp collisions at sqrt(s) = 8 TeV. The data correspond to an integrated luminosity of 19.5 inverse femtobarns collected with the CMS detector at the LHC. The measurements are performed using events with two oppositely charged leptons (electrons or muons) and two or more jets, where at least one of the jets is identified as originating from a bottom quark. The spin correlations and polarization are measured from the angular distributions of the two selected leptons, both inclusively and differentially, with respect to the invariant mass, rapidity, and transverse momentum of the t t-bar system. The measurements are unfolded to the parton level and found to be in agreement with predictions of the standard model. A search for new physics in the form of anomalous top quark chromo moments is performed. No evidence of new physics is observed, and exclusion limits on the real part of the chromo-magnetic dipole moment and the imaginary part of the chromo-electric dipole moment are evaluated.

Jet multiplicity distributions in top quark pair (t t-bar) events are measured in pp collisions at a centre-of-mass energy of 8 TeV with the CMS detector at the LHC using a data set corresponding to an integrated luminosity of 19.7 inverse femtobarns. The measurement is performed in the dilepton decay channels (e+ e-, mu+ mu-, and e+/- mu-/+). The absolute and normalized differential cross sections for t t-bar production are measured as a function of the jet multiplicity in the event for different jet transverse momentum thresholds and the kinematic properties of the leading additional jets. The differential t t-bar b and t t-bar b b-bar cross sections are presented for the first time as a function of the kinematic properties of the leading additional b jets. Furthermore, the fraction of events without additional jets above a threshold is measured as a function of the transverse momenta of the leading additional jets and the scalar sum of the transverse momenta of all additional jets. The data are compared and found to be consistent with predictions from several perturbative quantum chromodynamics event generators and a next-to-leading order calculation.

Measurements are presented of the cross section for the production of at least four jets, of which at least two originate from b quarks, in proton-proton collisions. Data collected with the CMS detector at the LHC at a center-of-mass energy of 7 TeV are used, corresponding to an integrated luminosity of 3 inverse picobarns. The cross section is measured as a function of the jet transverse momentum for pt > 20 GeV, and of the jet pseudorapidity for abs(eta) < 2.4 (b jets), 4.7 (untagged jets). The correlations in azimuthal angle and pt between the jets are also studied. The inclusive cross section is measured to be sigma(pp to 2 b + 2 j + X) = 69 +/- 3 (stat) +/- 24 (syst) nb. The eta and pt distributions of the four jets and the correlations between them are well reproduced by event generators that combine perturbative QCD calculations at next-to-leading-order accuracy with contributions from parton showers and multiparton interactions.