The transverse energy cross-sectiondσ/dET has been measured in the pseudorapidity region 0.6<η<2.4 for hadron-lead collisions at 200 GeV/c incident hadron momentum. TheET distribution extends to 40 GeV, which is twice the kinematic limit forp-p collisions at the same incident beam momentum. The distribution ofET is found to shift towards low pseudorapidities with increasing total transverse energy.

This letter describes the observation of the light-by-light scattering process, $\gamma\gamma\rightarrow\gamma\gamma$, in Pb+Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. The analysis is conducted using a data sample corresponding to an integrated luminosity of 1.73 nb$^{-1}$, collected in November 2018 by the ATLAS experiment at the LHC. Light-by-light scattering candidates are selected in events with two photons produced exclusively, each with transverse energy $E_{\textrm{T}}^{\gamma} > 3$ GeV and pseudorapidity $|\eta_{\gamma}| < 2.37$, diphoton invariant mass above 6 GeV, and small diphoton transverse momentum and acoplanarity. After applying all selection criteria, 59 candidate events are observed for a background expectation of 12 $\pm$ 3 events. The observed excess of events over the expected background has a significance of 8.2 standard deviations. The measured fiducial cross section is 78 $\pm$ 13 (stat.) $\pm$ 7 (syst.) $\pm$ 3 (lumi.) nb.

A search for the production of three massive vector bosons in proton--proton collisions is performed using data at $\sqrt{s}=13\,TeV$ recorded with the ATLAS detector at the Large Hadron Collider in the years 2015--2017, corresponding to an integrated luminosity of $79.8\,\text{fb}^{-1}$. Events with two same-sign leptons $\ell$ (electrons or muons) and at least two reconstructed jets are selected to search for $WWW\to\ell\nu\ell\nu qq$. Events with three leptons without any same-flavour opposite-sign lepton pairs are used to search for $WWW\to\ell\nu\ell\nu\ell\nu$, while events with three leptons and at least one same-flavour opposite-sign lepton pair and one or more reconstructed jets are used to search for $WWZ\to\ell\nu qq \ell\ell$. Finally, events with four leptons are analysed to search for $WWZ\to\ell\nu\ell\nu\ell\ell$ and $WZZ\to qq \ell\ell\ell\ell$. Evidence for the joint production of three massive vector bosons is observed with a significance of 4.0 standard deviations, where the expectation is 3.1 standard deviations.

A search for charged leptons with large impact parameters using 139 fb$^{-1}$ of $\sqrt{s} = 13$ TeV $pp$ collision data from the ATLAS detector at the LHC is presented, addressing a long-standing gap in coverage of possible new physics signatures. Results are consistent with the background prediction. This search provides unique sensitivity to long-lived scalar supersymmetric lepton-partners (sleptons). For lifetimes of 0.1 ns, selectron, smuon and stau masses up to 720 GeV, 680 GeV, and 340 GeV are respectively excluded at 95% confidence level, drastically improving on the previous best limits from LEP.

Using 1.8 fb-1 of pp collisions at a center-of-mass energy of 7 TeV recorded by the ATLAS detector at the Large Hadron Collider, we present measurements of the production cross sections of Upsilon(1S,2S,3S) mesons. Upsilon mesons are reconstructed using the di-muon decay mode. Total production cross sections for p_T<70 GeV and in the rapidity interval |Upsilon|<2.25 are measured to be 8.01+-0.02+-0.36+-0.31 nb, 2.05+-0.01+-0.12+-0.08 nb, 0.92+-0.01+-0.07+-0.04 nb respectively, with uncertainties separated into statistical, systematic, and luminosity measurement effects. In addition, differential cross section times di-muon branching fractions for Upsilon(1S), Upsilon(2S), and Upsilon(3S) as a function of Upsilon transverse momentum p_T and rapidity are presented. These cross sections are obtained assuming unpolarized production. If the production polarization is fully transverse or longitudinal with no azimuthal dependence in the helicity frame the cross section may vary by approximately +-20%. If a non-trivial azimuthal dependence is considered, integrated cross sections may be significantly enhanced by a factor of two or more. We compare our results to several theoretical models of Upsilon meson production, finding that none provide an accurate description of our data over the full range of Upsilon transverse momenta accessible with this dataset.

Jets created in association with a photon can be used as a calibrated probe to study energy loss in the medium created in nuclear collisions. Measurements of the transverse momentum balance between isolated photons and inclusive jets are presented using integrated luminosities of 0.49 nb$^{-1}$ of Pb+Pb collision data at $\sqrt{s_\mathrm{NN}}=5.02$ TeV and 25 pb$^{-1}$ of $pp$ collision data at $\sqrt{s}=5.02$ TeV recorded with the ATLAS detector at the LHC. Photons with transverse momentum $63.1 < p_\mathrm{T}^{\gamma} < 200$ GeV and $\left|\eta^{\gamma}\right| < 2.37$ are paired inclusively with all jets in the event that have $p_\mathrm{T}^\mathrm{jet} > 31.6$ GeV and pseudorapidity $\left|\eta^\mathrm{jet}\right| < 2.8$. The transverse momentum balance given by the jet-to-photon $p_\mathrm{T}$ ratio, $x_\mathrm{J\gamma}$, is measured for pairs with azimuthal opening angle $\Delta\phi > 7\pi/8$. Distributions of the per-photon jet yield as a function of $x_\mathrm{J\gamma}$, $(1/N_\gamma)(\mathrm{d}N/\mathrm{d}x_\mathrm{J\gamma})$, are corrected for detector effects via a two-dimensional unfolding procedure and reported at the particle level. In $pp$ collisions, the distributions are well described by Monte Carlo event generators. In Pb+Pb collisions, the $x_\mathrm{J\gamma}$ distribution is modified from that observed in $pp$ collisions with increasing centrality, consistent with the picture of parton energy loss in the hot nuclear medium. The data are compared with a suite of energy-loss models and calculations.

The inclusive production rates of isolated, prompt photons in $p$+Pb collisions at $\sqrt{s_\mathrm{NN}} = 8.16$ TeV are studied with the ATLAS detector at the Large Hadron Collider using a dataset with an integrated luminosity of 165 nb$^{-1}$ recorded in 2016. The cross-section and nuclear modification factor $R_{p\mathrm{Pb}}$ are measured as a function of photon transverse energy from 20 GeV to 550 GeV and in three nucleon-nucleon centre-of-mass pseudorapidity regions, (-2.83,-2.02), (-1.84,0.91), and (1.09,1.90). The cross-section and $R_{p\mathrm{Pb}}$ values are compared with the results of a next-to-leading-order perturbative QCD calculation, with and without nuclear parton distribution function modifications, and with expectations based on a model of the energy loss of partons prior to the hard scattering. The data disfavour a large amount of energy loss and provide new constraints on the parton densities in nuclei.

A measurement of the four-lepton invariant mass spectrum is made with the ATLAS detector, using an integrated luminosity of 36.1 fb$^{-1}$ of proton-proton collisions at $\sqrt{s}$ = 13 TeV delivered by the Large Hadron Collider. The differential cross-section is measured for events containing two same-flavour opposite-sign lepton pairs. It exhibits a rich structure, with different mass regions dominated in the Standard Model by single $Z$ boson production, Higgs boson production, and $Z$ boson pair production, and non-negligible interference effects at high invariant masses. The measurement is compared with state-of-the-art Standard Model calculations, which are found to be consistent with the data. These calculations are used to interpret the data in terms of $gg\rightarrow ZZ \rightarrow 4\ell$ and $Z \rightarrow 4\ell$ subprocesses, and to place constraints on a possible contribution from physics beyond the Standard Model.

Results from a search for supersymmetry in events with four or more charged leptons (electrons, muons and taus) are presented. The analysis uses a data sample corresponding to 36.1 fb$^{-1}$ of proton-proton collisions delivered by the Large Hadron Collider at $\sqrt{s}=13$ TeV and recorded by the ATLAS detector. Four-lepton signal regions with up to two hadronically decaying taus are designed to target a range of supersymmetric scenarios that can be either enriched in or depleted of events involving the production and decay of a $Z$ boson. Data yields are consistent with Standard Model expectations and results are used to set upper limits on the event yields from processes beyond the Standard Model. Exclusion limits are set at the 95% confidence level in simplified models of General Gauge Mediated supersymmetry, where higgsino masses are excluded up to 295 GeV. In $R$-parity-violating simplified models with decays of the lightest supersymmetric particle to charged leptons, lower limits of 1.46 TeV, 1.06 TeV, and 2.25 TeV are placed on wino, slepton and gluino masses, respectively.

The prevalence of hadronic jets at the LHC requires that a deep understanding of jet formation and structure is achieved in order to reach the highest levels of experimental and theoretical precision. There have been many measurements of jet substructure at the LHC and previous colliders, but the targeted observables mix physical effects from various origins. Based on a recent proposal to factorize physical effects, this Letter presents a double-differential cross-section measurement of the Lund jet plane using 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV proton-proton collision data collected with the ATLAS detector using jets with transverse momentum above 675 GeV. The measurement uses charged particles to achieve a fine angular resolution and is corrected for acceptance and detector effects. Several parton shower Monte Carlo models are compared with the data. No single model is found to be in agreement with the measured data across the entire plane.

This letter presents a combination of searches for Higgs boson pair production using up to 36.1 fb$^{-1}$ of proton-proton collision data at a centre-of-mass energy $\sqrt{s} = 13$ TeV recorded with the ATLAS detector at the LHC. The combination is performed using six analyses searching for Higgs boson pairs decaying into the bbbb, bbWW, bb$\tau\tau$, WWWW, bb$\gamma \gamma$ and WW$\gamma\gamma$ final states. Results are presented for non-resonant and resonant Higgs boson pair production modes. No statistically significant excess in data above the Standard Model predictions is found. The combined observed (expected) limit at 95% confidence level on the non-resonant Higgs boson pair production cross-section is 6.9 (10) times the predicted Standard Model cross-section. Limits are also set on the ratio ($ \kappa_{\lambda} $) of the Higgs boson self-coupling to its Standard Model value. This ratio is constrained at 95% confidence level in observation (expectation) to $ -5.0 < \kappa_{\lambda} <12.0 $ ($ -5.8 < \kappa_{\lambda} <12.0 $). In addition, limits are set on the production of narrow scalar resonances and spin-2 Kaluza-Klein Randall-Sundrum gravitons. Exclusion regions are also provided in the parameter space of the habemus Minimal Supersymmetric Standard Model and the Electroweak Singlet 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.

A search for direct pair production of top squarks in final states with two tau leptons, $b$-jets, and missing transverse momentum is presented. The analysis is based on proton-proton collision data at $\sqrt{s} = 13$ TeV corresponding to an integrated luminosity of 36.1 fb$^{-1}$ recorded with the ATLAS detector at the Large Hadron Collider in 2015 and 2016. Two exclusive channels with either two hadronically decaying tau leptons or one hadronically and one leptonically decaying tau lepton are considered. No significant deviation from the Standard Model predictions is observed in the data. The analysis results are interpreted in terms of model-independent limits and used to derive exclusion limits on the masses of the top squark $\tilde t_1$ and the tau slepton $\tilde \tau_1$ in a simplified model of supersymmetry with a nearly massless gravitino. In this model, masses up to $m(\tilde t_1) = 1.16$ TeV and $m(\tilde \tau_1) = 1.00$ TeV are excluded at 95% confidence level.

A search for magnetic monopoles and high-electric-charge objects is presented using 34.4 fb$^{-1}$ of 13 TeV $pp$ collision data collected by the ATLAS detector at the LHC during 2015 and 2016. The considered signature is based upon high ionization in the transition radiation tracker of the inner detector associated with a pencil-shape energy deposit in the electromagnetic calorimeter. The data were collected by a dedicated trigger based on the tracker high-threshold hit capability. The results are interpreted in models of Drell-Yan pair production of stable particles with two spin hypotheses (0 and 1/2) and masses ranging from 200 GeV to 4000 GeV. The search improves by approximately a factor of five the constraints on the direct production of magnetic monopoles carrying one or two Dirac magnetic charges and stable objects with electric charge in the range $20\le|z|\le60$ and extends the charge range to $60<|z|\le100$.

This paper presents a measurement of jet fragmentation functions in 0.49 nb$^{-1}$ of Pb+Pb collisions and 25 pb$^{-1}$ of $pp$ collisions at $\sqrt{s_{NN}} = 5.02$ TeV collected in 2015 with the ATLAS detector at the LHC. These measurements provide insight into the jet quenching process in the quark-gluon plasma created in the aftermath of ultra-relativistic collisions between two nuclei. The modifications to the jet fragmentation functions are quantified by dividing the measurements in Pb+Pb collisions by baseline measurements in $pp$ collisions. This ratio is studied as a function of the transverse momentum of the jet, the jet rapidity, and the centrality of the collision. In both collision systems, the jet fragmentation functions are measured for jets with transverse momentum between 126 GeV and 398 GeV and with an absolute value of jet rapidity less than 2.1. An enhancement of particles carrying a small fraction of the jet momentum is observed, which increases with centrality and with increasing jet transverse momentum. Yields of particles carrying a very large fraction of the jet momentum are also observed to be enhanced. Between these two enhancements of the fragmentation functions a suppression of particles carrying an intermediate fraction of the jet momentum is observed in Pb+Pb collisions. A small dependence of the modifications on jet rapidity is observed.

Multi-particle azimuthal cumulants are measured as a function of centrality and transverse momentum using 470 $\mu$b$^{-1}$ of Pb+Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV with the ATLAS detector at the LHC. These cumulants provide information on the event-by-event fluctuations of harmonic flow coefficients $v_n$ and correlated fluctuations between two harmonics $v_n$ and $v_m$. For the first time, a non-zero four-particle cumulant is observed for dipolar flow, $v_1$. The four-particle cumulants for elliptic flow, $v_2$, and triangular flow, $v_3$, exhibit a strong centrality dependence and change sign in ultra-central collisions. This sign change is consistent with significant non-Gaussian fluctuations in $v_2$ and $v_3$. The four-particle cumulant for quadrangular flow, $v_4$, is found to change sign in mid-central collisions. Correlations between two harmonics are studied with three- and four-particle mixed-harmonic cumulants, which indicate an anti-correlation between $v_2$ and $v_3$, and a positive correlation between $v_2$ and $v_4$. These correlations decrease in strength towards central collisions and either approach zero or change sign in ultra-central collisions. To investigate the possible flow fluctuations arising from intrinsic centrality or volume fluctuations, the results are compared between two different event classes used for centrality definitions. In peripheral and mid-central collisions where the cumulant signals are large, only small differences are observed. In ultra-central collisions, the differences are much larger and transverse momentum dependent. These results provide new information to disentangle flow fluctuations from the initial and final states, as well as new insights on the influence of centrality fluctuations.

The dynamics of isolated-photon plus two-jet production in $pp$ collisions at a centre-of-mass energy of 13 TeV are studied with the ATLAS detector at the LHC using a dataset corresponding to an integrated luminosity of 36.1 fb$^{-1}$. Cross sections are measured as functions of a variety of observables, including angular correlations and invariant masses of the objects in the final state, $\gamma+jet+jet$. Measurements are also performed in phase-space regions enriched in each of the two underlying physical mechanisms, namely direct and fragmentation processes. The measurements cover the range of photon (jet) transverse momenta from 150 GeV (100 GeV) to 2 TeV. The tree-level plus parton-shower predictions from SHERPA and PYTHIA as well as the next-to-leading-order QCD predictions from SHERPA are compared with the measurements. The next-to-leading-order QCD predictions describe the data adequately in shape and normalisation except for regions of phase space such as those with high values of the invariant mass or rapidity separation of the two jets, where the predictions overestimate the data.

Results of a search for the pair production of photon-jets$-$collimated groupings of photons$-$in the ATLAS detector at the Large Hadron Collider are reported. Highly collimated photon-jets can arise from the decay of new, highly boosted particles that can decay to multiple photons collimated enough to be identified in the electromagnetic calorimeter as a single, photonlike energy cluster. Data from proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 36.7 fb$^{-1}$, were collected in 2015 and 2016. Candidate photon-jet pair production events are selected from those containing two reconstructed photons using a set of identification criteria much less stringent than that typically used for the selection of photons, with additional criteria applied to provide improved sensitivity to photon-jets. Narrow excesses in the reconstructed diphoton mass spectra are searched for. The observed mass spectra are consistent with the Standard Model background expectation. The results are interpreted in the context of a model containing a new, high-mass scalar particle with narrow width, $X$, that decays into pairs of photon-jets via new, light particles, $a$. Upper limits are placed on the cross section times the product of branching ratios $\sigma \times \mathcal{B}(X \rightarrow aa) \times \mathcal {B}(a \rightarrow \gamma \gamma)^{2}$ for 200 GeV $< m_{X} <$ 2 TeV and for ranges of $ m_a $ from a lower mass of 100 MeV up to between 2 and 10 GeV, depending upon $ m_X $. Upper limits are also placed on $\sigma \times \mathcal{B}(X \rightarrow aa) \times \mathcal {B}(a \rightarrow 3\pi^{0})^{2}$ for the same range of $ m_X $ and for ranges of $ m_a $ from a lower mass of 500 MeV up to between 2 and 10 GeV.

This paper presents a measurement of quantities related to the formation of jets from high-energy quarks and gluons (fragmentation). Jets with transverse momentum 100 GeV $<p_T<$ 2.5 TeV and pseudorapidity $|\eta| < 2.1$ from an integrated luminosity of 33 fb$^{-1}$ of $\sqrt{s}=13$ TeV proton-proton collisions are reconstructed with the ATLAS detector at the Large Hadron Collider. Charged-particle tracks with $p_T > 500$ MeV and $|\eta| < 2.5$ are used to probe the detailed structure of the jet. The fragmentation properties of the more forward and the more central of the two leading jets from each event are studied. The data are unfolded to correct for detector resolution and acceptance effects. Comparisons with parton shower Monte Carlo generators indicate that existing models provide a reasonable description of the data across a wide range of phase space, but there are also significant differences. Furthermore, the data are interpreted in the context of quark- and gluon-initiated jets by exploiting the rapidity dependence of the jet flavor fraction. A first measurement of the charged-particle multiplicity using model-independent jet labels (topic modeling) provides a promising alternative to traditional quark and gluon extractions using input from simulation. The simulations provide a reasonable description of the quark-like data across the jet $p_T$ range presented in this measurement, but the gluon-like data have systematically fewer charged particles than the simulations.

Azimuthal decorrelations between the two central jets with the largest transverse momenta are sensitive to the dynamics of events with multiple jets. We present a measurement of the normalized differential cross section based on the full dataset (L=36/pb) acquired by the ATLAS detector during the 2010 sqrt(s)=7 TeV proton-proton run of the LHC. The measured distributions include jets with transverse momenta up to 1.3 TeV, probing perturbative QCD in a high energy regime.

Correlations of two flow harmonics $v_n$ and $v_m$ via three- and four-particle cumulants are measured in 13 TeV $pp$, 5.02 TeV $p$+Pb, and 2.76 TeV peripheral Pb+Pb collisions with the ATLAS detector at the LHC. The goal is to understand the multi-particle nature of the long-range collective phenomenon in these collision systems. The large non-flow background from dijet production present in the standard cumulant method is suppressed using a method of subevent cumulants involving two, three and four subevents separated in pseudorapidity. The results show a negative correlation between $v_2$ and $v_3$ and a positive correlation between $v_2$ and $v_4$ for all collision systems and over the full multiplicity range. However, the magnitudes of the correlations are found to depend strongly on the event multiplicity, the choice of transverse momentum range and collision system. The relative correlation strength, obtained by normalisation of the cumulants with the $\langle v_n^2\rangle$ from a two-particle correlation analysis, is similar in the three collision systems and depends weakly on the event multiplicity and transverse momentum. These results based on the subevent methods provide strong evidence of a similar long-range multi-particle collectivity in $pp$, $p$+Pb and peripheral Pb+Pb collisions.

A search for heavy long-lived multi-charged particles is performed using the ATLAS detector at the LHC. Data with an integrated luminosity of 36.1 fb$^{-1}$ collected in 2015 and 2016 from proton-proton collisions at $\sqrt{s}$ = 13 TeV are examined. Particles producing anomalously high ionization, consistent with long-lived massive particles with electric charges from |q|=2e to |q|=7e, are searched for. No events are observed, and 95% confidence level cross-section upper limits are interpreted as lower mass limits for a Drell-Yan production model. Multi-charged particles with masses between 50 GeV and 980-1220 GeV (depending on their electric charge) are excluded.

This paper reports on a search for heavy resonances decaying into $WW$, $ZZ$ or $WZ$ using proton-proton collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV. The data, corresponding to an integrated luminosity of 139 $\mathrm{fb^{-1}}$, were recorded with the ATLAS detector from 2015 to 2018 at the Large Hadron Collider. The search is performed for final states in which one $W$ or $Z$ boson decays leptonically, and the other $W$ boson or $Z$ boson decays hadronically. The data are found to be described well by expected backgrounds. Upper bounds on the production cross sections of heavy scalar, vector or tensor resonances are derived in the mass range 300-5000 GeV within the context of Standard Model extensions with warped extra dimensions or including a heavy vector triplet. Production through gluon-gluon fusion, Drell-Yan or vector-boson fusion are considered, depending on the assumed model.

A dedicated sample of Large Hadron Collider proton-proton collision data at centre-of-mass energy $\sqrt{s}=8$ TeV is used to study inclusive single diffractive dissociation, $pp \rightarrow Xp$. The intact final-state proton is reconstructed in the ATLAS ALFA forward spectrometer, while charged particles from the dissociated system $X$ are measured in the central detector components. The fiducial range of the measurement is $-4.0 < \log_{10} \xi < -1.6$ and $0.016 < |t| < 0.43 \ {\rm GeV^2}$, where $\xi$ is the proton fractional energy loss and $t$ is the squared four-momentum transfer. The total cross section integrated across the fiducial range is $1.59 \pm 0.13 \ {\rm mb}$. Cross sections are also measured differentially as functions of $\xi$, $t$, and $\Delta \eta$, a variable that characterises the rapidity gap separating the proton and the system $X$. The data are consistent with an exponential $t$ dependence, ${\rm d} \sigma / {\rm d} t \propto \text{e}^{Bt}$ with slope parameter $B = 7.65 \pm 0.34 \ {\rm GeV^{-2}}$. Interpreted in the framework of triple Regge phenomenology, the $\xi$ dependence leads to a pomeron intercept of $\alpha(0) = 1.07 \pm 0.09$.

A measurement of $W^\pm$ boson production in lead-lead collisions at $\sqrt{s_\mathrm{NN}} = 5.02$ TeV is reported using data recorded by the ATLAS experiment at the LHC in 2015, corresponding to a total integrated luminosity of $0.49\;\mathrm{nb^{-1}}$. The $W^\pm$ bosons are reconstructed in the electron or muon leptonic decay channels. Production yields of leptonically decaying $W^\pm$ bosons, normalised by the total number of minimum-bias events and the nuclear thickness function, are measured within a fiducial region defined by the detector acceptance and the main kinematic requirements. These normalised yields are measured separately for $W^+$ and $W^-$ bosons, and are presented as a function of the absolute value of pseudorapidity of the charged lepton and of the collision centrality. The lepton charge asymmetry is also measured as a function of the absolute value of lepton pseudorapidity. In addition, nuclear modification factors are calculated using the $W^\pm$ boson production cross-sections measured in $pp$ collisions. The results are compared with predictions based on next-to-leading-order calculations with CT14 parton distribution functions as well as with predictions obtained with the EPPS16 and nCTEQ15 nuclear parton distribution functions. No dependence of normalised production yields on centrality and a good agreement with predictions are observed for mid-central and central collisions. For peripheral collisions, the data agree with predictions within 1.7 (0.9) standard deviations for $W^-$ ($W^+$) bosons.