We report the measurement of $K^{*0}$ meson at midrapidity ($|y|<$ 1.0) in Au+Au collisions at $\sqrt{s_{\rm NN}}$~=~7.7, 11.5, 14.5, 19.6, 27 and 39 GeV collected by the STAR experiment during the RHIC beam energy scan (BES) program. The transverse momentum spectra, yield, and average transverse momentum of $K^{*0}$ are presented as functions of collision centrality and beam energy. The $K^{*0}/K$ yield ratios are presented for different collision centrality intervals and beam energies. The $K^{*0}/K$ ratio in heavy-ion collisions are observed to be smaller than that in small system collisions (e+e and p+p). The $K^{*0}/K$ ratio follows a similar centrality dependence to that observed in previous RHIC and LHC measurements. The data favor the scenario of the dominance of hadronic re-scattering over regeneration for $K^{*0}$ production in the hadronic phase of the medium.

Measurements of $ZZ$ production in the $\ell^{+}\ell^{-}\ell^{\prime +}\ell^{\prime -}$ channel in proton-proton collisions at 13 TeV center-of-mass energy at the Large Hadron Collider are presented. The data correspond to 36.1 $\mathrm{fb}^{-1}$ of collisions collected by the ATLAS experiment in 2015 and 2016. Here $\ell$ and $\ell'$ stand for electrons or muons. Integrated and differential $ZZ \to \ell^{+}\ell^{-}\ell^{\prime +}\ell^{\prime -}$ cross sections with $Z \to \ell^+\ell^-$ candidate masses in the range of 66 GeV to 116 GeV are measured in a fiducial phase space corresponding to the detector acceptance and corrected for detector effects. The differential cross sections are presented in bins of twenty observables, including several that describe the jet activity. The integrated cross section is also extrapolated to a total phase space and to all Standard-Model decays of $Z$ bosons with mass between 66 GeV and 116 GeV, resulting in a value of $17.3 \pm 0.9$ [$\pm 0.6$ (stat.) $\pm 0.5$ (syst.) $\pm 0.6$ (lumi.)] pb. The measurements are found to be in good agreement with the Standard-Model predictions. A search for neutral triple gauge couplings is performed using the transverse momentum distribution of the leading $Z$-boson candidate. No evidence for such couplings is found and exclusion limits are set on their parameters.

The production of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonances has been measured in p-Pb collisions at $\sqrt{s_{\rm NN}}$ = 8.16 TeV using the ALICE detector. Resonances are reconstructed via their hadronic decay channels in the rapidity interval $-$0.5 $<$$y$$<$ 0 and the transverse momentum spectra are measured for various multiplicity classes up to $p_{\rm T}$ = 20 GeV/$c$ for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $p_{\rm T}$ = 16 GeV/$c$ for $\mathrm{\phi(1020)}$. The $p_{\rm T}$ -integrated yields and mean transverse momenta are reported and compared with previous results in pp, p-Pb and Pb-Pb collisions. The $x_{\mathrm{T}}$ scaling for $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ resonance production is newly tested in p-Pb collisions and found to hold in the high-$p_{\rm T}$ region at LHC energies. The nuclear modification factors ($R_{\rm pPb}$) as a function of $p_{\rm T}$ for $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ at $\sqrt{s_{NN}}$ = 8.16 TeV are presented along with the new $R_{\rm pPb}$ measurements of $\mathrm{K}^{*0}$, $\mathrm{\phi}$ , $\Xi$, and $\Omega$ at $\sqrt{s_{\rm NN}}$ = 5.02 TeV. At intermediate $p_{\rm T}$ (2-8 GeV/$c$), $R_{\rm pPb}$ of $\Xi$, $\Omega$ show a Cronin-like enhancement, while $\mathrm{K}^{*0}$ and $\mathrm{\phi}$ show no or little nuclear modification. At high $p_{\rm T}$ ($>$ 8 GeV/$c$), the $R_{\rm pPb}$ values of all hadrons are consistent with unity within uncertainties. The $R_{\rm pPb}$ of $\mathrm{K}^{*}(\mathrm{892})^{0}$ and $\mathrm{\phi(1020)}$ at $\sqrt{s_{\rm NN}}$ = 8.16 and 5.02 TeV show no significant energy dependence.

An invariant differential cross section measurement of inclusive $\pi^{0}$ and $\eta$ meson production at mid-rapidity in pp collisions at $\sqrt{s}=8$ TeV was carried out by the ALICE experiment at the LHC. The spectra of $\pi^{0}$ and $\eta$ mesons were measured in transverse momentum ranges of $0.3<p_{\rm T}<35$ GeV/$c$ and $0.5<p_{\rm T}<35$ GeV/$c$, respectively. Next-to-leading order perturbative QCD calculations using fragmentation functions DSS14 for the $\pi^{0}$ and AESSS for the $\eta$ overestimate the cross sections of both neutral mesons, although such calculations agree with the measured $\eta/\pi^{0}$ ratio within uncertainties. The results were also compared with PYTHIA~8.2 predictions for which the Monash~2013 tune yields the best agreement with the measured neutral meson spectra. The measurements confirm a universal behavior of the $\eta/\pi^{0}$ ratio seen for NA27, PHENIX and ALICE data for pp collisions from $\sqrt{s}=27.5$ GeV to $\sqrt{s}=8$ TeV within experimental uncertainties. A relation between the $\pi^{0}$ and $\eta$ production cross sections for pp collisions at $\sqrt{s}=8$ TeV is given by $m_{\rm T}$ scaling for $p_{\rm T}>3.5$ GeV/$c$. However, a deviation from this empirical scaling rule is observed for transverse momenta below $p_{\rm T}<3.5$ GeV/$c$ in the $\eta/\pi^0$ ratio with a significance of $6.2\sigma$.

A new measurement of the proton structure function $F_2(x,Q~2)$ is reported for momentum transfers squared $Q~2$ between 1.5GeV$~2$ and 5000GeV$~2$ and for Bjorken $x$ between $3\cdot 10~{-5}$ and $0.32$ using data collected by the HERA experiment H1 in 1994. The data represent an increase in statistics by a factor of ten with respect to the analysis of the 1993 data. Substantial extension of the kinematic range towards low $Q~2$ and $x$ has been achieved using dedicated data samples and events with initial state photon radiation. The structure function is found to increase significantly with decreasing $x$, even in the lowest accessible $Q~2$ region. The data are well described by a Next to Leading Order QCD fit and the gluon density is extracted.

The production rate of charged D* mesons in jets has been measured in 1.8-TeV p¯p collisions at the Fermilab Tevatron with the Collider Detector at Fermilab. In a sample of approximately 32 300 jets with a mean transverse energy of 47 GeV obtained from an exposure of 21.1 nb−1, a signal corresponding to 25.0±7.5(stat)±2.0(syst) D*±→K∓π±π± events is seen above background. This corresponds to a ratio N(D*++D*−)/N(jet) =0.10±0.03±0.03 for D* mesons with fractional momentum z greater than 0.1.

We report our final results from the analysis of the full high statistics sample of events of the reaction ν μ + e − → μ − + ν c collected with the CHARM II detector in the CERN wide-band neutrino beam during the years 1988 to 1991. From a signal of 15758 ± 324 inverse muon decay events we derived, inthe Born approximation, a value of (16.51 ± 0.93) × 10 −42 cm 2 GeV −1 for the asymptotic cross section slope σ E ν , in goodagreement with the Standard Model prediction of 17.23 × 10 −42 cm 2 GeV −1 . The result constrains the scalar coupling of the electron and the muon to | g LL S | 2 < 0.475 at 90% CL.

We study charged particle production in proton-antiproton collisions at 300 GeV, 900 GeV, and 1.96 TeV. We use the direction of the charged particle with the largest transverse momentum in each event to define three regions of eta-phi space; toward, away, and transverse. The average number and the average scalar pT sum of charged particles in the transverse region are sensitive to the modeling of the underlying event. The transverse region is divided into a MAX and MIN transverse region, which helps separate the hard component (initial and final-state radiation) from the beam-beam remnant and multiple parton interaction components of the scattering. The center-of-mass energy dependence of the various components of the event are studied in detail. The data presented here can be used to constrain and improve QCD Monte Carlo models, resulting in more precise predictions at the LHC energies of 13 and 14 TeV.

Deep inelastic scattering (DIS) events, selected from 1993 data taken by the H1 experiment at HERA, are studied in the Breit frame of reference. The fragmentation function of the quark is compared with those of \ee data. It is shown that certain aspects of the quarks emerging from within the proton in \ep interactions are essentially the same as those of quarks pair-created from the vacuum in \ee annihilation. The measured area, peak position and widthof the fragmentation function show that the kinematic evolution variable, equivalent to the \ee squared centre of mass energy, is in the Breit frame the invariant square of the four-momentum transfer. We comment on the extent to which we have evidence for coherence effects in parton showers.

The Standard Model of particle physics describes the known fundamental particles and forces that make up our universe, with the exception of gravity. One of the central features of the Standard Model is a field that permeates all of space and interacts with fundamental particles. The quantum excitation of this field, known as Higgs field, manifests itself as the Higgs boson, the only fundamental particle with no spin. In 2012, a particle with properties consistent with the Higgs boson of the Standard Model was observed by the ATLAS and CMS experiments at the Large Hadron Collider at CERN. Since then, more than 30 times as many Higgs bosons have been recorded by the ATLAS experiment, allowing much more precise measurements and new tests of the theory. Here, on the basis of this larger dataset, we combine an unprecedented number of production and decay processes of the Higgs boson to scrutinize its interactions with elementary particles. Interactions with gluons, photons, and $W$ and $Z$ bosons -- the carriers of the strong, electromagnetic, and weak forces -- are studied in detail. Interactions with three third-generation matter particles (bottom ($b$) and top ($t$) quarks, and tau leptons ($\tau$)) are well measured and indications of interactions with a second-generation particle (muons, $\mu$) are emerging. These tests reveal that the Higgs boson discovered ten years ago is remarkably consistent with the predictions of the theory and provide stringent constraints on many models of new phenomena beyond the Standard Model.

Jet substructure quantities are measured using jets groomed with the soft-drop grooming procedure in dijet events from 32.9 fb$^{-1}$ of $pp$ collisions collected with the ATLAS detector at $\sqrt{s} = 13$ TeV. These observables are sensitive to a wide range of QCD phenomena. Some observables, such as the jet mass and opening angle between the two subjets which pass the soft-drop condition, can be described by a high-order (resummed) series in the strong coupling constant $\alpha_S$. Other observables, such as the momentum sharing between the two subjets, are nearly independent of $\alpha_S$. These observables can be constructed using all interacting particles or using only charged particles reconstructed in the inner tracking detectors. Track-based versions of these observables are not collinear safe, but are measured more precisely, and universal non-perturbative functions can absorb the collinear singularities. The unfolded data are directly compared with QCD calculations and hadron-level Monte Carlo simulations. The measurements are performed in different pseudorapidity regions, which are then used to extract quark and gluon jet shapes using the predicted quark and gluon fractions in each region. All of the parton shower and analytical calculations provide an excellent description of the data in most regions of phase space.

The energy spectrum and the cross section of photonuclear interactions of 180 GeV muons in iron were measured at the CERN SPS using prototype modules of the ATLAS hadron calorimeter. The differential

A search for new phenomena is performed in final states containing one or more jets and an imbalance in transverse momentum in pp collisions at a centre-of-mass energy of 13 TeV. The analysed data sample, recorded with the CMS detector at the CERN LHC, corresponds to an integrated luminosity of 2.3 inverse femtobarns. Several kinematic variables are employed to suppress the dominant background, multijet production, as well as to discriminate between other standard model and new physics processes. The search provides sensitivity to a broad range of new-physics models that yield a stable weakly interacting massive particle. The number of observed candidate events is found to agree with the expected contributions from standard model processes, and the result is interpreted in the mass parameter space of fourteen simplified supersymmetric models that assume the pair production of gluinos or squarks and a range of decay modes. For models that assume gluino pair production, masses up to 1575 and 975 GeV are excluded for gluinos and neutralinos, respectively. For models involving the pair production of top squarks and compressed mass spectra, top squark masses up to 400 GeV are excluded.

Statistical combinations of searches for charginos and neutralinos using various decay channels are performed using $139\,$fb$^{-1}$ of $pp$ collision data at $\sqrt{s}=13\,$TeV with the ATLAS detector at the Large Hadron Collider. Searches targeting pure-wino chargino pair production, pure-wino chargino-neutralino production, or higgsino production decaying via Standard Model $W$, $Z$, or $h$ bosons are combined to extend the mass reach to the produced SUSY particles by 30-100 GeV. The depth of the sensitivity of the original searches is also improved by the combinations, lowering the 95% CL cross-section upper limits by 15%-40%.

Infrared and collinear safe event shape distributions and their mean values are determined in e+e- collisions at centre-of-mass energies between 45 and 202 GeV. A phenomenological analysis based on power correction models including hadron mass effects for both differential distributions and mean values is presented. Using power corrections, alpha_s is extracted from the mean values and shapes. In an alternative approach, renormalisation group invariance (RGI) is used as an explicit constraint, leading to a consistent description of mean values without the need for sizeable power corrections. The QCD beta-function is precisely measured using this approach. From the DELPHI data on Thrust, including data from low energy experiments, one finds beta_0 = 7.86 +/- 0.32 for the one loop coefficient of the beta-function or, assuming QCD, n_f = 4.75 +/- 0.44 for the number of active flavours. These values agree well with the QCD expectation of beta_0=7.67 and n_f=5. A direct measurement of the full logarithmic energy slope excludes light gluinos with a mass below 5 GeV.

This paper reviews and extends searches for the direct pair production of the scalar supersymmetric partners of the top and bottom quarks in proton-proton collisions collected by the ATLAS collaboration during the LHC Run 1. Most of the analyses use 20 fb$^{-1}$ of collisions at a centre-of-mass energy of $\sqrt{s}$ = 8 TeV, although in some case an additional 4.7 fb$^{-1}$ of collision data at $\sqrt{s}$ = 7 TeV are used. New analyses are introduced to improve the sensitivity to specific regions of the model parameter space. Since no evidence of third-generation squarks is found, exclusion limits are derived by combining several analyses and are presented in both a simplified model framework, assuming simple decay chains, as well as within the context of more elaborate phenomenological supersymmetric models.

Charged-particle spectra at midrapidity are measured in Pb-Pb collisions at the centre-of-mass energy per nucleon-nucleon pair $\sqrt{s_{\rm NN}}$ = 5.02 TeV and presented in centrality classes ranging from most central (0-5%) to most peripheral (95-100%) collisions. Possible medium effects are quantified using the nuclear modification factor ($R_{\rm AA}$) by comparing the measured spectra with those from proton-proton collisions, scaled by the number of independent nucleon-nucleon collisions obtained from a Glauber model. At large transverse momenta ($8<p_{\rm T}<20$ GeV/$c$), the average $R_{\rm AA}$ is found to increase from about $0.15$ in 0-5% central to a maximum value of about $0.8$ in 75-85% peripheral collisions, beyond which it falls off strongly to below $0.2$ for the most peripheral collisions. Furthermore, $R_{\rm AA}$ initially exhibits a positive slope as a function of $p_{\rm T}$ in the $8$-$20$ GeV/$c$ interval, while for collisions beyond the 80% class the slope is negative. To reduce uncertainties related to event selection and normalization, we also provide the ratio of $R_{\rm AA}$ in adjacent centrality intervals. Our results in peripheral collisions are consistent with a PYTHIA-based model without nuclear modification, demonstrating that biases caused by the event selection and collision geometry can lead to the apparent suppression in peripheral collisions. This explains the unintuitive observation that $R_{\rm AA}$ is below unity in peripheral Pb-Pb, but equal to unity in minimum-bias p-Pb collisions despite similar charged-particle multiplicities.

We present an angular analysis of the $B^{+}\rightarrow K^{\ast+}(\rightarrow K_{S}^{0}\pi^{+})\mu^{+}\mu^{-}$ decay using 9$\,\mbox{fb}^{-1}$ of $pp$ collision data collected with the LHCb experiment. For the first time, the full set of CP-averaged angular observables is measured in intervals of the dimuon invariant mass squared. Local deviations from Standard Model predictions are observed, similar to those in previous LHCb analyses of the isospin-partner $B^{0}\rightarrow K^{\ast0}\mu^{+}\mu^{-}$ decay. The global tension is dependent on which effective couplings are considered and on the choice of theory nuisance parameters.

An angular analysis of the $B^{0}\rightarrow K^{*0}(\rightarrow K^{+}\pi^{-})\mu^{+}\mu^{-}$ decay is presented. The dataset corresponds to an integrated luminosity of $3.0\,{\mbox{fb}^{-1}}$ of $pp$ collision data collected at the LHCb experiment. The complete angular information from the decay is used to determine $C\!P$-averaged observables and $C\!P$ asymmetries, taking account of possible contamination from decays with the $K^{+}\pi^{-}$ system in an S-wave configuration. The angular observables and their correlations are reported in bins of $q^2$, the invariant mass squared of the dimuon system. The observables are determined both from an unbinned maximum likelihood fit and by using the principal moments of the angular distribution. In addition, by fitting for $q^2$-dependent decay amplitudes in the region $1.1<q^{2}<6.0\mathrm{\,Ge\kern -0.1em V}^{2}/c^{4}$, the zero-crossing points of several angular observables are computed. A global fit is performed to the complete set of $C\!P$-averaged observables obtained from the maximum likelihood fit. This fit indicates differences with predictions based on the Standard Model at the level of 3.4 standard deviations. These differences could be explained by contributions from physics beyond the Standard Model, or by an unexpectedly large hadronic effect that is not accounted for in the Standard Model predictions.

Measurements of elliptic ($v_2$) and triangular ($v_3$) flow coefficients of $\pi^{\pm}$, K$^{\pm}$, p+$\rm \overline{p}$, K$^0_{\rm S}$, and $\Lambda + \overline{\Lambda}$ obtained with the scalar product method in Xe-Xe collisions at $\sqrt{s_{\rm NN}}$ = 5.44 TeV are presented. The results are obtained in the rapidity range $\left | y \right |<0.5$ and reported as a function of transverse momentum, $p_{\rm T}$, for several collision centrality classes. The flow coefficients exhibit a particle mass dependence for $p_{\rm T}<3$ GeV/$c$, while a grouping according to particle type (i.e., meson and baryon) is found at intermediate transverse momenta (3< $p_{\rm T}$ <8 GeV/$c$). The magnitude of the baryon $v_{2}$ is larger than that of mesons up to $p_{\rm T}$ = 6 GeV/$c$. The centrality dependence of the shape evolution of the $p_{\rm T}$-differential $v_2$ is studied for the various hadron species. The $v_2$ coefficients of $\pi^{\pm}$, K$^{\pm}$, and p+$\rm \overline{p}$ are reproduced by MUSIC hydrodynamic calculations coupled to a hadronic cascade model (UrQMD) for $p_{\rm T} <1$ GeV/$c$. A comparison with $v_{\rm n}$ measurements in the corresponding centrality intervals in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV yields an enhanced $v_2$ in central collisions and diminished value in semicentral collisions.

The elliptic ($v_2$), triangular ($v_3$), and quadrangular ($v_4$) flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, ${\rm p+\overline{p}}$, ${\Lambda+\overline{\Lambda}}$, ${\rm K}^{\rm 0}_{\rm S}$, and the $\phi$-meson are measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}=5.02$ TeV. Results obtained with the scalar product method are reported for the rapidity range $\vert y \vert <$ 0.5 as a function of transverse momentum, $p_\text{T}$, at different collision centrality intervals between 0-70%, including ultra-central (0-1%) collisions for $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$. For $p_\text{T} < 3$ GeV$\kern-0.05em/\kern-0.02em c$, the flow coefficients exhibit a particle mass dependence. At intermediate transverse momenta ($3<p_\text{T}<$~8-10 GeV$\kern-0.05em/\kern-0.02em c$), particles show an approximate grouping according to their type (i.e., mesons and baryons). The $\phi$-meson $v_2$, which tests both particle mass dependence and type scaling, follows ${\rm p+\overline{p}}$ $v_2$ at low $p_\text{T}$ and $\pi^{\pm}$ $v_2$ at intermediate $p_\text{T}$. The evolution of the shape of $v_{\rm n}(p_{\mathrm{T}})$ as a function of centrality and harmonic number $n$ is studied for the various particle species. Flow coefficients of $\pi^{\pm}$, ${\rm K}^{\pm}$, and ${\rm p+\overline{p}}$ for $p_\text{T}<3$ GeV$\kern-0.05em/\kern-0.02em c$ are compared to iEBE-VISHNU and MUSIC hydrodynamical calculations coupled to a hadronic cascade model (UrQMD). The iEBE-VISHNU calculations describe the results fairly well for $p_\text{T} < 2.5$ GeV$\kern-0.05em/\kern-0.02em c$, while MUSIC calculations reproduce the measurements for $p_\text{T} < 1$ GeV$\kern-0.05em/\kern-0.02em c$. A comparison to $v_{\rm n}$ coefficients measured in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV is also provided.

Azimuthal anisotropy of produced particles is one of the most important observables used to access the collective properties of the expanding medium created in relativistic heavy-ion collisions. In this paper, we present second ($v_{2}$) and third ($v_{3}$) order azimuthal anisotropies of $K_{S}^{0}$, $\phi$, $\Lambda$, $\Xi$ and $\Omega$ at mid-rapidity ($|y|<$1) in Au+Au collisions at $\sqrt{s_{\text{NN}}}$ = 54.4 GeV measured by the STAR detector. The $v_{2}$ and $v_{3}$ are measured as a function of transverse momentum and centrality. Their energy dependence is also studied. $v_{3}$ is found to be more sensitive to the change in the center-of-mass energy than $v_{2}$. Scaling by constituent quark number is found to hold for $v_{2}$ within 10%. This observation could be evidence for the development of partonic collectivity in 54.4 GeV Au+Au collisions. Differences in $v_{2}$ and $v_{3}$ between baryons and anti-baryons are presented, and ratios of $v_{3}$/$v_{2}^{3/2}$ are studied and motivated by hydrodynamical calculations. The ratio of $v_{2}$ of $\phi$ mesons to that of anti-protons ($v_{2}(\phi)/v_{2}(\bar{p})$) shows centrality dependence at low transverse momentum, presumably resulting from the larger effects from hadronic interactions on anti-proton $v_{2}$.

We present systematic measurements of azimuthal anisotropy for strange and multistrange hadrons ($K^{0}_{s}$, $\Lambda$, $\Xi$, and $\Omega$) and $\phi$ mesons at midrapidity ($|y| <$ 1.0) in collisions of U + U nuclei at $\sqrt{s_{NN}} = 193$ GeV, recorded by the STAR detector at the Relativistic Heavy Ion Collider. Transverse momentum ($p_{\text{T}}$) dependence of flow coefficients ($v_{2}$, $v_{3}$, and $v_{4}$) is presented for minimum bias collisions and three different centrality intervals. Number of constituent quark scaling of the measured flow coefficients in U + U collisions is discussed. We also present the ratio of $v_{n}$ scaled by the participant eccentricity ($\varepsilon_{n}\left\lbrace 2 \right\rbrace$) to explore system size dependence and collectivity in U + U collisions. The magnitude of $v_{2}/\varepsilon_{2}$ is found to be smaller in U + U collisions than that in central Au + Au collisions contradicting naive eccentricity scaling. Furthermore, the ratios between various flow harmonics ($v_{3}/v_{2}^{3/2}$, $v_{4}/v_{2}^{4/2}$) are studied and compared with hydrodynamic and transport model calculations.

The Fourier coefficients v[2] and v[3] characterizing the anisotropy of the azimuthal distribution of charged particles produced in PbPb collisions at sqrt(s[NN]) = 5.02 TeV are measured with data collected by the CMS experiment. The measurements cover a broad transverse momentum range, 1 < pT < 100 GeV. The analysis focuses on pT > 10 GeV range, where anisotropic azimuthal distributions should reflect the path-length dependence of parton energy loss in the created medium. Results are presented in several bins of PbPb collision centrality, spanning the 60% most central events. The v[2] coefficient is measured with the scalar product and the multiparticle cumulant methods, which have different sensitivities to the initial-state fluctuations. The values of both methods remain positive up to pT of about 60-80 GeV, in all examined centrality classes. The v[3] coefficient, only measured with the scalar product method, tends to zero for pT greater than or equal to 20 GeV. Comparisons between theoretical calculations and data provide new constraints on the path-length dependence of parton energy loss in heavy ion collisions and highlight the importance of the initial-state fluctuations.

Azimuthal correlations between the two jets with the largest transverse momenta $ {p_{\mathrm{T}}} $ in inclusive 2-, 3-, and 4-jet events are presented for several regions of the leading jet $ {p_{\mathrm{T}}} $ up to 4 TeV. For 3- and 4-jet scenarios, measurements of the minimum azimuthal angles between any two of the three or four leading $ {p_{\mathrm{T}}} $ jets are also presented. The analysis is based on data from proton-proton collisions collected by the CMS Collaboration at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. Calculations based on leading-order matrix elements supplemented with parton showering and hadronization do not fully describe the data, so next-to-leading-order calculations matched with parton shower and hadronization models are needed to better describe the measured distributions. Furthermore, we show that azimuthal jet correlations are sensitive to details of the parton showering, hadronization, and multiparton interactions. A next-to-leading-order calculation matched with parton showers in the MC@NLO method, as implemented in HERWIG 7, gives a better overall description of the measurements than the POWHEG method.