The Relativistic Heavy Ion Colider (RHIC) showed that strongly interacting matter produced in collision of two heavy nuclei takes the form of quark gluon plasma (QGP) . In such medium produced color particles are expected to lose energy which leads to a phenomenon known as jet quenching. Suppression of charge hadron yields in heavy-ion collisions was already reported by experiments at RHIC and the LHC [2,3,4].
A direct way to study the initial stage of proton-proton (pp), proton-nucleus and nucleus-nucleus collisions is to measure production of colorless electroweak (EW) bosons, namely high-energy photons, W and Z. These bosons are created at the very early stage of the collision in the hard parton-parton interaction. Leptonic decays of W and Z are expected to not interact substantially with the QGP, making them useful probes of the QGP. This unique feature makes them sensitive to the initial geometry of the nuclei and nuclear modifications to parton distribution functions (nPDF) can be studied. A measurement of direct photon production  by the PHENIX experiment at RHIC showed that their rates scale with the nuclear thickness function. Very similar conclusions came from the first measurements of electroweak vector bosons in lead-lead (Pb+Pb) collisions performed by the ATLAS and CMS experiments [6-10]. In principle their production rates were found to be independent of the presence of QGP. Furthermore, the weak boson production in proton-lead (p+Pb) and Pb+Pb collisions may differ from the pp system. There are two main sources of these differences. The first one is related to the presence of neutrons in the nucleus which have different quark composition compared to protons. Secondly, EW boson productions can be affected by effects arising from partons being bound in the nucleus where the PDF of the free nucleon might be modified leading to parton depletion or enhancement .
This report covers the most recent measurements on W and Z boson production in the pp, p+Pb and Pb+Pb systems at TeV by the ATLAS Collaboration. They are based on data collected in 2013 and 2015. A detailed description of the ATLAS detector is provided in Ref. .
2. Z boson production
The Z boson production has been measured in pp, p+Pb and Pb+Pb collisions at TeV using data corresponding to an integrated luminosity of 29 nb and 0.49 nb [13,14], respectively. Figure 1 presents Z boson yields scaled by the average nuclear thickness function as a function for three selected centrality intervals compared to the measurement in the pp system. The Z bosons are reconstructed via the di-muon decay channel. Good agreement with the model predictions and measured pp data points is observed. Only slight excess is observed in the most peripheral class of events. The nuclear modification factor is also shown in the lower panel. Good agreement with unity indicates small impact of nuclear effects.
Figure 2 presents the differential cross section as a function of Z boson rapidity measured in p+Pb at TeV. Model estimations based on the CT10 PDF  set with and without EPS09  nuclear corrections and MSTW2008 are also shown. The strong asymmetry around comes from the asymmetry of the collision. A test of compatibility between data and the model shapes (irrespective of normalization) shows that the CT10+EPS09 shape of the Z boson rapidity disitrbution gives a highest p-value of 0.79 .
3. W boson production
W boson production in the muon decay mode has been measured in Pb+Pb collisions at TeV using data corresponding to integrated luminosity of 0.49 nb . The production yields are calculated within the fiducial region defined by muon pseudorapidity, , muon transverse momentum, GeV, missing transverse energy, GeV and transverse mass of the muon and missing transverse momentum, GeV. The total numbers of background subtracted and efficiency corrected events for W production are 34500 290 (syst.) and 31300 270 (syst.) for W and W respectively. Figure 3 presents the difference in production yields between positive and negative W bosons called as charge asymmetry and defined by:
where is a corrected number of W (W ) bosons. Charge asymmetry is extracted from the 0-80% centrality interval as a function of the muon absolute pseudorapidity. This observable is sensitive to PDFs and may also be sensitive to the nuclear modifications. The data are compared to Powheg scaled by using CT10 free nucleon PDF and MCFM using the most recent nuclear modifications to PDF, namely EPPS16  and nCTEQ15 . Sensitivity to nuclear modifications is not observed as all three predictions agree well. The muon charge asymmetry is roughly constant for and the data are described well by the predictions. For higher data points deviate from the expectations by no more than 3 .
Measurements of Z and W boson production have been reported based on pp, p+Pb and Pb+Pb data collected at TeV by the ATLAS experiment at the LHC. Predictions based on perturbative QCD calculations describe the data well. Differences between the data and predictions are observed for the Pb-going side, in case of p+Pb collisions, where there appears to be an excess above the model. The muon charge asymmetry for Pb+Pb collisions is in good agreement with all three predictions, except for the forward direction where there are discrepancies between the data and predictions at the level up to 3 .