Luminosity determination using Z boson production at the CMS experiment

The measurement of Z boson production is presented as a method to determine the integrated luminosity of CMS data sets. The analysis uses proton–proton collision data, recorded by the CMS experiment at the CERN LHC in 2017 at a center-of-mass energy of 13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {Te\hspace{-.08em}V}$$\end{document}TeV. Events with Z bosons decaying into a pair of muons are selected. The total number of Z bosons produced in a fiducial volume is determined, together with the identification efficiencies and correlations from the same data set, in small intervals of 20\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {pb}^{-1}$$\end{document}pb-1 of integrated luminosity, thus facilitating the efficiency and rate measurement as a function of time and instantaneous luminosity. Using the ratio of the efficiency-corrected numbers of Z bosons, the precisely measured integrated luminosity of one data set is used to determine the luminosity of another. For the first time, a full quantitative uncertainty analysis of the use of Z bosons for the integrated luminosity measurement is performed. The uncertainty in the extrapolation between two data sets, recorded in 2017 at low and high instantaneous luminosity, is less than 0.5%. We show that the Z boson rate measurement constitutes a precise method, complementary to traditional methods, with the potential to improve the measurement of the integrated luminosity.

Events with a Z boson decaying into a pair of muons have a remarkably clean experimental signature and a large cross section that facilitates high-precision measurements.Samples of Z bosons are also used as standard tools for detector calibrations and efficiency studies.The precisely known Z boson mass and width [15] are used to calibrate energy scales and momenta and to determine the detector resolution [16,17].Efficiencies for lepton triggering, reconstruction, and identification are determined using the "tag-and-probe" method [1,[16][17][18].
The large Drell-Yan (DY) cross section for the production of Z bosons, and the possibility of simultaneously determining both the yield and the detection efficiency in situ, i.e., from the same event sample, make the process useful for precision measurements of the integrated luminosity.This was discussed before the start of the LHC [19].During LHC operation, measurements of the Z boson rate already proved to be a useful and independent method for the LHC machine operators and experiments to monitor the relative instantaneous luminosity delivered to the ATLAS and CMS experiments [20].The use of Z boson production as a measure of relative luminosities was also explored by the ATLAS experiment [21].
Both muons from the Z boson decay are detectable within the fiducial volume of the CMS detector in about one third of the Z boson events.The fiducial Z boson cross section in protonproton (pp) collisions at 13 TeV has been measured to be [3] σ Z B(Z → µ + µ − ) = 694 ± 6 (syst) ± 17 (lumi) pb. ( Theoretical predictions are available up to next-to-next-to-next-to-leading order (N 3 LO) [9] in quantum chromodynamics (QCD).Electroweak corrections, including mixed QCD-electroweak corrections, are also available [7,22,23].The current uncertainty in the prediction of the fiducial cross section is about 3%, and mainly originates from limited knowledge of proton PDFs and higher-order corrections [8].Within this uncertainty, the integrated luminosity can be directly determined from the measured number of Z bosons corrected for efficiencies.
In practice, precision luminosity calibrations at the LHC are obtained from van der Meer (vdM) scan data [21,[24][25][26][27][28][29], which are more precise than the theory predictions for the Z boson cross section.In vdM scans, which are performed at low instantaneous luminosity with zero crossing angle between the two beams, the two beams are separated in two orthogonal directions transverse to the parallel beam axes.In each scan step, for a given beam separation, the event rate measured in the luminosity detectors is recorded to determine the beam overlap area.Together with the beam currents and the measured head-on collision rate, a luminosity calibration constant, referred to as the visible cross section, is determined.A full vdM scan campaign takes about six hours per experiment and is usually performed once per year, with specifically configured beams to maximize the accuracy and precision of the measurement.A detailed description of vdM scans is reported in Ref. [29].
The most precise integrated luminosity measurement in CMS to date, achieved for the 2016 data-taking period, has a total uncertainty of 1.2% [29].Roughly half of the total uncertainty is due to the luminosity integration over the full year of data taking.This uncertainty, in turn, is composed of the uncertainty in the extrapolation of the visible cross section obtained in the vdM scan to standard data-taking conditions at high instantaneous luminosity, and the uncertainty in the integration of the instantaneous luminosity over time, obtained from comparisons between different luminosity measurements.In the 2017 data, presented in this paper, the average number of collisions per bunch crossing, usually referred to as pileup, was 32 [30].In Run 2, peak instantaneous luminosities as high as 20 nb −1 s −1 were reached, corresponding to a pileup of more than 50.For the high-luminosity LHC (HL-LHC), a pileup of up to 200 is expected [31] likely leading to an increase in uncertainty with the conventional methods due to the larger extrapolation.
In this paper, we explore an approach originally proposed in Ref. [32].The measurement of the Z boson rate is used as an alternative method for the extrapolation and integration of the luminosity calibration.The Z boson counting complements conventional luminosity measurements obtained from the CMS luminosity systems, which are taken as reference luminosity.
The fiducial Z boson production cross section is defined as σ Z fid = N Z /L, where N Z stands for the efficiency-corrected number of reconstructed Z boson events and L for the integrated luminosity.Since σ Z fid is identical for all data sets of the same center-of-mass energy, the ratio of N Z for two data sets can be used to transfer the luminosity calibration from one data set to another, without input from theoretical predictions or precise knowledge of σ Z fid .For the first time, a full quantitative uncertainty analysis of the use of Z bosons for the integrated luminosity measurement is performed.
We choose two independent data sets of Z boson events, both recorded by the CMS experiment in 2017: a data set with a bunch luminosity corresponding to about three pp collisions per bunch crossing, referred to in the following as "lowPU"; and the bulk of CMS pp collision data recorded in 2017 with a typical pileup of 32, denoted as "highPU".The luminosity of the lowPU data is used to determine that of the highPU data via the relation For both sets of data, the individual trigger and selection efficiencies are determined in situ, in intervals of 20 pb −1 , thus enhancing the sensitivity to possible variations due to changes in beam conditions or detector response as a function of time.Using the integrated luminosity for the lowPU data, which has an uncertainty of 1.7% [33], the integrated luminosity L highPU and its uncertainty are determined from Eq. ( 2), and compared with the result from the conventional integrated luminosity measurement.Due to the cleaner signature, better resolution, and smaller backgrounds, the analysis of Z boson decays into muons is more accurate than electrons.In this paper, only the decays of Z bosons into muons are used.
The paper is structured as follows.After a brief outline of the CMS detector in Section 2, the analysis of the Z boson event sample is described in Section 3. The reconstructed number of Z bosons and the trigger and selection efficiencies are extracted from fits to the data.Acceptance corrections and correlations between the efficiencies for the muon track components and among the two muon tracks are studied as a function of pileup.Subsequently, in Section 4, the luminosity information obtained from Z boson counting is compared with the results from conventional luminosity measurements.In Section 5, the benefits and advantages of Z boson counting for luminosity measurements are discussed.The paper concludes with a summary in Section 6.

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the magnet volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter, each composed of a barrel and two endcap sections.Forward calorimeters extend the pseudorapidity (η) coverage provided by the barrel and endcap detectors.The muon system consists of gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid.A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, is reported in Ref. [34].
An iterative approach is used to build tracker tracks, executing a sequence of tracking algorithms, each with slightly distinct logic [17].Muons are measured in the range |η| < 2.4, with detection planes made using three technologies: drift tubes, cathode strip chambers, and resistive plate chambers.Matching muons to tracks measured in the silicon tracker results in a relative transverse momentum (p T ) resolution of 1% in the barrel and 3% in the endcaps, for muons with p T of about 100 GeV [17].The particle-flow (PF) algorithm [37] reconstructs and identifies each individual particle in an event, combining information from the various CMS detector components.Jets are clustered using the anti-k T jet finding algorithm [38,39] with the tracks assigned to candidate vertices as inputs, and the associated missing transverse momentum p miss T , taken as the negative vector p T sum of those jets [40].The primary vertex (PV) is taken to be the vertex with the largest ∑ p 2 T of its associated tracks, as described in Section 9.4 of Ref. [41].
Events of interest are selected using a two-tiered trigger system.The first level (L1), comprised of custom hardware processors, uses information from the calorimeters and muon detectors to select events at a rate of around 100 kHz within a fixed latency of 4 µs [42].The second level, known as the high-level trigger (HLT), consists of a farm of processors running a version of the full event reconstruction software optimized for fast processing, and reduces the event rate to around 1 kHz before data storage [43].
During LHC Run 2, the main CMS luminosity subdetectors (luminometers) were the silicon pixel detector, the hadron forward calorimeter (HF), the pixel luminosity telescope (PLT) [44], and the fast beam conditions monitor (BCM1F) [45].A separate data acquisition system is used to collect and store HF, PLT, and BCM1F data, as well as LHC beam-related data.A more detailed description of the CMS luminosity system is reported in Ref. [29].For all comparisons in this paper, the reference integrated luminosity is obtained with the CMS luminometers, calibrated as described in Ref. [33] and using offline-calibrated corrections for the afterglow effects in the HF luminosity measurement.
The analysis described in this paper is largely independent of Monte Carlo (MC) simulations.However, MC simulations are used for two purposes: to determine the expected DY invariant mass distribution of the signal measured in the CMS detector; and to study possible biases in the pileup-dependent measurement of the muon track-finding efficiencies.Simulated event samples of the DY process, Z/γ * → ℓℓ, are produced at leading order using the MAD-GRAPH5 aMC@NLO (v2.6.5)[46] generator, interfaced with PYTHIA (v8.240) [47] for the parton shower simulation.The parameters describing the modeling of the parton shower and underlying event are based on the CP5 tune [48].The generated MC events are passed through a full simulation of the detector using GEANT4 [49].

The Z boson candidate selection and efficiency determination
The events were recorded using a single-muon trigger (HLT muon) that requires at least one muon candidate with p T > 24 GeV and loose isolation criteria [50].The lowPU data were recorded using different, looser trigger configurations than those used for the highPU data.To obtain identical trigger configurations for the two data sets, the trigger decision in the lowPU was recalculated from raw data using the trigger configuration of the highPU data.
Based on the offline reconstruction, selected muon candidates consist of an "outer" standalone track in the muon system, matched to an "inner" track reconstructed in the silicon tracker [35].The outer track is required to have signals in at least two muon detector planes.The inner track must have at least one valid hit in the silicon pixel detector and hits in more than five strip tracker layers.The matching is done by comparing parameters of the two tracks propagated onto a common surface.A combined Kalman filter fit [51] is performed in which the information from the inner and outer tracks is used to obtain a "global" muon track.For global muons, the inner and outer tracks are required to have p T > 20 GeV, lie within |η| < 2.4, and to be matched within ∆R = √ (∆η) 2 + (∆ϕ) 2 < 0.3.Quality criteria on the global muon track fit are imposed, and it is required that the muon candidate is also reconstructed with the PF algorithm [37].No requirements are imposed on the impact parameters of the muon track.Isolation criteria are omitted to maintain efficiency also at high pileup.For muons with p T < 200 GeV, i.e., about 99% of identified muon candidates, the track parameters are taken from the inner track.In other cases, the track parameters are determined by combining information from the inner and outer tracks.For all muon tracks, p T > 25 GeV is required to ensure that the trigger efficiency reaches a plateau.
A Z boson candidate is identified as a pair of opposite-charge muons with an invariant mass of 60 < m µµ < 120 GeV.At least one of the two muon candidates is required to be matched with an HLT muon within ∆R < 0.1.To obtain the actual number of produced Z bosons, the number of reconstructed and selected Z boson candidates, the trigger efficiency, the muonidentification efficiency, as well as the background arising from nonresonant production, are determined from dedicated fits to the data, as explained in the following.

Trigger efficiency and signal extraction
The trigger efficiency and the number of Z boson candidates are determined from fits to the invariant dimuon mass distributions of mutually exclusive sets of events with exactly one (N 1 ) or exactly two (N 2 ) selected muons matched to an HLT muon.The observables N 1 and N 2 follow the relations (3) Here, the quantity ϵ µ HLT refers to the HLT muon trigger efficiency.The correction factor C HLT accounts for the correlation between the HLT efficiencies of the two muons.A value of C HLT > 1 indicates a positive correlation between the two muons, i.e., an increased probability for the second muon to pass the HLT if the first muon passes it.The determination of C HLT is presented in Section 3.2.The terms N , respectively.Also indicated are the χ 2 values per degree of freedom (dof).The lower panels contain the pulls of the distributions, defined as the difference between the data and the fit model in each bin, divided by the statistical uncertainty estimated from the expected number of entries given by the model.template generated from simulated Z → µµ events, convolved with a Gaussian function to take into account muon momentum scale and resolution differences between data and simulation.A falling exponential function is used to describe the nonresonant background.In Fig. 1, examples of two distributions and the results of the fits are presented.The sample shown here corresponds to an integrated luminosity of 20 pb −1 , yielding about 12 000 Z boson candidates.

Muon trigger correlation
The correlation between the trigger efficiencies of the two HLT muons is described by the correction factor C HLT , as introduced in Eq. ( 3).The dependence of C HLT on the pileup is of particular interest in this analysis because it does not cancel in the ratio in Eq. ( 2), and thus constitutes an important source of systematic uncertainty.The correlation was investigated in simulation, and it is largely understood to originate from isolation requirements in the trigger selection.
We determine C HLT from an MC simulation sample of Z → µµ events.As a proxy to the amount of pileup in a given event, we use the number of reconstructed PVs, N PV , an observable that is directly accessible event-by-event in both data and simulation.At fixed pileup, the distribution of N PV approximately follows a Poisson distribution with a mean at about 80% of the true pileup, as determined from DY simulation.
In the simulation, C HLT is obtained directly, by rearranging Eq. (3), as where N   We use data to validate the result for C HLT obtained in the simulation.To this end, events are analyzed that are triggered independently of the muon trigger, namely by using the trigger condition p miss T > 120 GeV in which the contribution from muons is not included.This p miss T trigger also records Z boson candidates for which the number of HLT muons is zero, and, thus, an additional relation for the number of reconstructed Z boson candidates with no HLT muons, denoted as N 0 , is obtained, (5) Together with Eq. (3), we obtain three equations for N 0 , N 1 , and N 2 with three unknowns, ϵ µ HLT , C HLT , and ϵ Z ID N Z .The correction factor C HLT can thus be determined from the number of signal events in the three categories, each obtained from a fit.The fits are performed separately in six bins of N PV where the number of bins and their boundaries are chosen such that the number of events per bin are similar.
The result is presented in Fig. 2. The red lines indicate the expectation from the simulation in which C HLT is at the level of 0.1-0.2%above unity for N PV ∼ 30.Within the limited statistical precision of the data, good agreement of the simulation with the data is observed.We assign a systematic uncertainty of 100% of the correction, which is represented by the gray band in the figure.

Muon identification and reconstruction efficiency
The efficiency to reconstruct a Z boson, ϵ Z ID , depends on the muon identification and reconstruction efficiency ϵ µ ID for each of the two muons.In the simulation, the pileup-dependent correlation between the two identified muons is of the order of 0.01%, and thus ϵ The value for C ID ≈ 1.0001 is taken from simulation and applied as a function of N PV .The muon efficiency ϵ µ ID is defined independently of the HLT muon efficiency, such that the total number of produced Z bosons is obtained from Eq. (3).
To determine ϵ µ ID , the following factorization ansatz is used: where the efficiency ϵ µ ID|Glo is the fraction of global muons that fulfill the full set of muon identification requirements; the efficiency ϵ µ Glo|Sta is the global muon efficiency, given by the fraction of standalone muons that also qualify as global muon; and the efficiency ϵ µ Sta|Trk is the standalone muon efficiency, defined as the fraction of muons with good inner tracks that are matched within ∆R < 0.3 to outer standalone muon tracks with p T > 20 GeV and |η| < 2.4.To obtain an unbiased set of inner tracks for the measurement of the efficiency ϵ µ Sta|Trk , inner tracks that are seeded from the extrapolation of outer standalone muon tracks are excluded.The term c T&P accounts for the correlations between the efficiency terms in Eq. ( 6).The pileup dependence of the correction from c T&P between the lowPU and the highPU data sets is estimated from simulation to be about 0.01%.
The efficiencies are determined from the data using a "tag-and-probe" methodology [1].Identified muon candidates that are matched to the HLT muon are selected as "tag".For each tag, a probe muon candidate of opposite charge is selected under the condition that the muon candidate pair has an invariant mass between 60 and 120 GeV.The efficiency ϵ µ x|y is then measured as where y denotes the reference sample of muon candidates and x is the probe criterion.The numbers n p and n f correspond to the number of events that pass and fail the test criterion, respectively.
For each of the efficiencies, and in bins of 20 pb −1 , fits to the m µµ distributions of the passing and failing distributions are performed.In the fits, the same shapes as described in Section 3.1 are used to describe the signal.In the histograms with passing probes, the background contribution is low and a falling exponential is used.In the case of failing probes, the nonresonant background is much larger and a more complex analytic function, comprising an exponential at high mass above the Z boson resonance and an error function at low mass, is fit.To ensure a bias-free measurement of ϵ µ Glo|Sta , the outer standalone muon track parameters are used to determine m µµ for the passing and failing probes.Since the resolution of these tracks is much worse, the invariant mass requirement is widened to 50-130 GeV.In the case that, in a given event, the probe muon also fulfills the tag muon requirements, the tag-and-probe muons are indistinguishable and both muons are used as probes.Quantitative results for the measurement of the efficiencies are presented in Section 4.

Acceptance correction
To determine the true number of Z bosons in the visible phase space, an acceptance correction for losses, or gains, due to the finite resolution of the reconstructed muon tracks is required.The correction affects the number of reconstructed Z bosons itself.The efficiencies are also affected, primarily in the matching of inner and outer tracks, and, to a lesser extent, if muon tracks for passing and for failing probes have different resolutions.The size of the correction is determined from the simulation by comparing the efficiency-corrected number of Z bosons as obtained from the measurement with the generated number of Z bosons in the visible phase space, as defined for bare leptons after final-state radiation (FSR), but before detector simulation.
For outer muon tracks, resolution effects lead to an acceptance correction of about 1.35%, which is independent of pileup and constant over the full year of data taking.For inner tracks, the acceptance correction is 0.15% at low pileup, and it is negligibly small for the highPU data set.This pileup dependency of 0.15% is applied as an additional correction, and an uncertainty of 100% of the correction is assigned.For a direct cross section measurement, the size of the bias could be further reduced through optimized track selection criteria.However, for this analysis it suffices that the pileup-independent components of the acceptance correction cancel in the cross section ratio.

The L1 trigger corrections
The term "prefiring" describes the effect that a trigger decision is assigned to a bunch crossing preceding the one in which the collision actually took place.In the CMS experiment, the triggering and readout of events in adjacent bunch crossings is vetoed in the trigger logic.However, due to the limited time resolution of the muon system, the assignment of muon candidates to bunch crossings can be wrong, and thus lead to a loss of good events, i.e., a trigger inefficiency.Since the tag-and-probe efficiency measurement is insensitive to this effect, the inefficiency due to prefiring is measured in a dedicated analysis.During the 2017 data taking, measurable prefiring occurred at nonnegligible rates for the L1 muon and ECAL triggers [42,50].For the L1 muon trigger, a correction for trigger inefficiency of 0.6% was found, independent of pileup and time.In contrast, losses due to prefiring of the ECAL require a pileup-dependent correction of 0.05-0.2%for the pileup range 0-50.The prefiring from ECAL triggers is caused mainly by initial or final state radiation, pileup jets, or the underlying event.
The impact on the lowPU data was somewhat larger due to the lower ECAL trigger thresholds, and for the lowPU data a correction of 0.6% is applied.

Results and uncertainties
The procedures described above to measure the number of reconstructed Z bosons and their efficiencies are applied to the data in bins of 20 pb −1 .Since the amount of data at the end of a fill does not usually add up to 20 pb −1 , the last bin is included as long as it contains more than 10 pb −1 .In case the last bin contains <10 pb −1 it is merged with the second to last.Altogether, in 2017 about 2000 such bins are defined.

Normalized Z boson rate
In Fig. 3, the measured Z boson rate and efficiencies are shown for the data recorded during a typical LHC fill of the highPU data-taking period in 2017.In this fill, pp collision data were recorded continuously for about 16 hours.An integrated luminosity of about 515 pb −1 was accumulated, corresponding to 25 bins of 20 pb −1 each, and the last bin contains the remaining 15 pb −1 of data.The instantaneous luminosity decreased from initially 15 nb −1 s −1 , corresponding to a pileup of about 50, to about one third of the initial value.In Fig. 3 (left), a comparison between the conventional measurement of the recorded luminosity and the measurement using the Z boson rate is shown.The integral of the measured Z boson rate is normalized to the integral of the reference luminosity.The shapes of the two independent measurements agree very well.In Fig. 3 (right), the muon trigger and identification efficiencies, ϵ µ HLT and ϵ µ ID , separated into its different components, as applied to the respective time intervals, are presented.In particular, a significant dependence on time, and thus on pileup, is seen for the HLT muon efficiency for which a rise by about 3% is measured as the pileup decreases in the course of the fill.To compare the relative linearity between the measurement of the Z boson rate and the CMS reference luminosity, the fiducial cross section for Z boson production, normalized to the average Z boson cross section, is studied as a function of the instantaneous luminosity.The result is shown in Fig. 4, where the average instantaneous luminosity in each 20 pb −1 bin is used to assign an instantaneous luminosity bin from which the average cross section is obtained.The straight-line fit to the data yields a value of 0.2% below unity for the intercept with the y-axis at low pileup.This value gives an estimate of the agreement between Z boson counting and reference luminosity measurement in the extrapolation from the low to the high pileup data.

Measurement of the absolute luminosity
Using Eq. ( 2), the integrated luminosity of the highPU data, referred to in the following as "Z luminosity", is determined from the integrated luminosity of the lowPU data, using the ratio of the number of Z bosons recorded during the two periods, corrected for reconstruction and trigger efficiencies as determined in intervals of 20 pb −1 .In the ratio, all correlated uncertainties cancel, as detailed in the following section.
In Fig. 5, the distribution of the ratios between the Z luminosity and the reference luminosity as obtained from the CMS luminosity systems is shown.Each entry in the histogram corresponds to an interval of 20 pb −1 in the highPU data recorded in 2017.The central values of both measurements are in good agreement with a difference of 0.3%.The standard deviation of about 1.2% is predominantly of statistical nature, and close to the expectation for the pure statistical uncertainty of about 12 000 Z boson candidates reconstructed in intervals of 20 pb −1 each.The ratio of Z luminosity and reference luminosity as a function of the integrated luminosity is shown in Fig. 6.This figure shows a good stability of the Z luminosity measurement over the full year.No significant patterns in time are observed.

Statistical and systematic uncertainties, and additional cross checks
The uncertainties in the analysis were studied with the focus on the ratio r = N Z highPU /N Z lowPU of the Z boson counts between two data samples in 2017 as presented in Eq. ( 2).The full list of considered sources of uncertainty in the cross sections and their ratio is given in Table 1, and described in the following.
Statistical uncertainties are driven by the number of available Z bosons and also include the statistical uncertainty in the efficiencies.As mentioned above, in one interval of 20 pb −1 , about 12 000 Z bosons candidates with two muons in the final state are available, leading to an average statistical uncertainty of 1.17%.For all intervals combined, the statistical uncertainty for the full 2017 highPU data is negligibly small.The lowPU data set corresponds to an integrated luminosity of about 200 pb −1 , and this contributes a statistical uncertainty of about 0.35%.The luminosity as measured from Z bosons divided by the reference luminosity as a function of the integrated luminosity for the 2017 highPU data.Each green point represents the ratio from one measurement of the number of Z bosons.The blue lines show the averages of 50 consecutive measurements, corresponding to an average of 1 fb −1 of data.The gray band has a width of 1.5%, corresponding to the uncertainty in the ratio of the integrated reference luminosities from the lowPU to the one of highPU [33].As discussed in Section 3.2, the correction factor for correlations in the trigger efficiencies of the two muons C HLT is determined from data and simulation; it is about 0.1% above unity for the highPU sample, consistently for data and MC simulation.The uncertainty in C HLT is assigned to be 100% of the correction.
Possible correlations between the two identified muons and imperfect factorization of muon identification and reconstruction efficiencies were discussed in Section 3.3.The simulation shows negligible effects, and corrections at the level of 0.01% are applied.The corresponding uncertainties are estimated to be 100% of the correction.
The limited resolution of the reconstructed muon tracks leads to a bias in the measurement, as described in Section 3.4.The bias from the inner track resolution is smaller, but pileup dependent, and remains in the ratio with a magnitude of 0.15%.The outer track resolution leads to a large bias, but is mostly pileup independent and cancels in the ratio.A correction is derived from simulation and two independent sources of uncertainty, estimated to be 100% of the correction each, are assigned each for the inner and outer tracks, respectively.
Systematic uncertainties in the L1 muon prefiring corrections, described in Section 3.5, cancel completely in the ratio, whereas the uncertainties due to ECAL prefiring have a different magnitude between the two data sets and cancel only partially.The remaining uncertainty is estimated to be 20% of the nominal correction [42].
The extraction of the signal and background contributions was studied using alternative fit models.For the signal model, the Gaussian resolution function convolved with the histogram template is varied.First, the histogram template is used alone, i.e., fully relying on the simulation and leaving no further degrees of freedom to the fit.Secondly, the histogram template is convolved with a Crystal Ball function [52], which has four free parameters and gives the fit more freedom to incorporate possible differences between data and simulation.Thirdly, the histogram template is constructed from generator-level post-FSR leptons, instead of mirroring the selection at detector level.This template is then convolved with a Crystal Ball function.
The three variations lead to changes in the extracted numbers of Z bosons, and the efficiencies, in both the highPU and lowPU data sets.While the trends are correlated, the relative magnitudes are different, and this leads to a significant residual uncertainty.The envelope of the three variations is taken to quantify this uncertainty.
The two types of background models are varied independently.For the categories with major background contributions, the Das function [53], a wide Gaussian distribution with exponential tails is used as an alternative function, which has four free parameters, as opposed to three for the nominal model.In the other cases, the falling exponential is substituted by a uniform distribution.
The total systematic uncertainty is obtained by adding the systematic uncertainties listed in Table 1 in quadrature.In combination with the statistical uncertainty of 0.35%, the total uncertainty to transfer the luminosity from the lowPU data to the highPU data in 2017 is δr = +0.31%−0.28% (syst) ± 0.35% (stat) = +0.47%−0.45% , where the statistical uncertainty is due to the limited size of the lowPU data set.The systematic uncertainty is driven by the uncertainties in the signal modeling, followed by the background modeling and acceptance corrections.Overall, a total uncertainty of about 0.5% is obtained.
Multiple cross-checks were performed to test the robustness of the result.The size of the luminosity bin was varied from 20 down to 15 and up to 30 pb −1 , and negligible differences with respect to the nominal measurement were found.It was further verified that the measurement is independent of the choice of the bin width chosen for the m µµ distribution, by varying it by factors of 1/2 and 2. Independence of the result on the chosen fit interval was tested using two alternative ranges: a more narrow interval from m µµ ∈ [70, 110] GeV and a wider interval from GeV.Both variations have a strong impact on N Z since the phase space of the measurement changes, but, as expected, the effect cancels almost completely in the ratio.The results of these cross checks are summarized in Table 2.

Discussion and outlook
With an uncertainty in the transfer factor of about 0.5% for the 2017 data, this analysis shows that Z boson counting can provide an independent and competitive method to extrapolate and integrate luminosity calibrations.The results from Z boson counting are independent of the conventional luminosity measurements.They can be treated as uncorrelated in combinations, which can lead to significant improvements in the combined uncertainty.
Taking the current precision of 1.7% for the integrated luminosity in the lowPU data [33], the integrated luminosity in the highPU 2017 data could potentially be determined to a precision of better than 1.8%, in contrast to the preliminary uncertainty of the reference luminosity measurement of 2.3% [33].
A unique aspect of Z boson counting is that the relevant efficiency corrections as a function of time can be calibrated from the same event sample.This feature makes the method robust not only against small changes in detector response, but also across different detector configurations.In general, once a precision measurement of the integrated luminosity is available, such as that for the lowPU data in 2017, the integrated luminosity for all data recorded at the same center-of-mass energy can be determined using the Z boson counting.However, each transfer between data sets requires detailed studies of the correlations of the muon trigger and the reconstruction efficiencies.
In this paper, the full analysis was presented for the data from 2017, when a dedicated and sufficiently large sample of lowPU data was recorded.Under such conditions, a large fraction of the systematic uncertainties cancels in the ratio.For the most precise CMS measurement of the luminosity to date [29], published for 2016, an extrapolation and integration uncertainty of 0.7% was reported.For 2016, no lowPU data set was recorded.Further studies on the impact of different detector conditions would be required to extrapolate from the 2016 data set.If, hypothetically, an extrapolation uncertainty of 0.5% for Z boson counting were achievable also in the 2016 data, the uncertainty of 1.2% in the total integrated luminosity for 2016 could be improved to 1.1%.
The dominant contribution to the uncertainty comes from the statistical uncertainty, which is driven by the size of the lowPU data sample.The lowPU data recorded in 2017 correspond to an integrated luminosity of about 200 pb −1 .A significant increase of the sample size, e.g., by a factor 3 or 4, would make the statistical uncertainty negligible.
In the coming years, during the ongoing LHC Run 3 and beyond, additional measurements and studies on the main systematic uncertainties will be performed, and that is expected to improve the precision of the method further.Furthermore, the method is expected to contribute substantially to the combination of integrated luminosity measurements for different data sets.
In the longer term, pileup conditions of up to 200 pp collisions per bunch crossing are expected at the HL-LHC [31].In both Run 3 and at the HL-LHC, the uncertainties due to extrapolation from vdM conditions to standard data taking are expected to remain substantial.In such conditions, the method of Z boson counting has the potential to provide significant improvements.

Summary
The precision measurement of the Z boson production rate provides a complementary method to transfer integrated luminosity measurements between data sets.This study makes use of events with Z bosons decaying into a pair of muons.The data were recorded with the CMS experiment at the CERN LHC in 2017, at a proton-proton center-of-mass energy of 13 TeV.
The integrated luminosity of a larger data sample recorded in 2017 is obtained from that of a smaller data set recorded at lower pileup using the ratio of the efficiency-corrected numbers of Z bosons counted in the two data sets.The full set of efficiencies and correlation correction factors for triggering, reconstruction, and selection are determined in intervals of 20 pb −1 from the same Z boson data samples.Monte Carlo simulations are used only to describe the shape of the resonant Z boson signal and for the study of possible biases of the method.A detailed quantitative study of the systematic uncertainties and their dependencies on pileup is performed for the first time.In the integrated luminosity ratio, the systematic uncertainties cancel almost completely, with the exception of the pileup-dependent effects.The resulting uncertainty in the ratio is 0.5%.With its high precision, the Z boson counting is competitive with and independent of conventional methods for the extrapolation and integration of luminosity.

Figure 1 :
Figure 1: The upper panels show the reconstructed invariant mass distributions of Z boson candidates in a 20 pb −1 sample of data for events where one (left) or two (right) muons pass the single-muon trigger selection.The blue curve shows the fitted background contribution and the red curve illustrates the modeled signal-plus-background contribution.The error bars indicate the statistical uncertainties.The numbers of signal and background candidates are given by N sig i = N i − N bkg i and N bkg i sig 1 and N sig 2 are the number of signal events, corresponding to N 1 − N bkg 1 and N 2 − N bkg 2 in the data.

Figure 2 :
Figure 2: Correction factor C HLT for the correlation between the measured muon trigger efficiencies of the two muons as a function of the number of reconstructed primary vertices, N PV , in the simulation (lines) and the data (points).The data points are drawn at the mean value of N PV in each bin of the measurement.The horizontal error bars on the points show the bin width, and the vertical error bars show the statistical uncertainty.The gray band indicates the ±100% uncertainty in the correction factor.

Figure 3 :
Figure3: Left: the efficiency-corrected Z boson rate, compared to the reference luminosity measurement, in the LHC fill 6255, recorded on September 29, 2017[33].Each bin corresponds to about 20 pb −1 , as determined by the reference measurement.For shape comparison, the integrated Z boson rate is normalized to the reference integrated luminosity.The panel at the bottom shows the ratio of the two measurements.The vertical error bars show the statistical uncertainty in the Z boson rate.Right: the measured single-muon efficiencies as functions of time for the same LHC fill.The vertical error bars show the statistical uncertainty in the efficiency.

Figure 4 :
Figure 4: Fiducial Z boson production cross section as a function of the instantaneous recorded luminosity, normalized to the average measured cross section.In each point, multiple measurements of the delivered Z boson rates are combined, the error bars correspond to the statistical uncertainties of the Z boson rate measurement.The leftmost point, highlighted in red, corresponds to the lowPU data.The result of a fit to a linear function is shown as a red line and the statistical uncertainties are covered by the gray band.

Figure 5 :
Figure 5: Distribution of the ratio of integrated luminosities between Z boson counting and the reference luminometer.The entries, each corresponding to one interval of 20 pb −1 of highPU data, are weighted with the respective measured luminosity.

Figure 6 :
Figure6: The luminosity as measured from Z bosons divided by the reference luminosity as a function of the integrated luminosity for the 2017 highPU data.Each green point represents the ratio from one measurement of the number of Z bosons.The blue lines show the averages of 50 consecutive measurements, corresponding to an average of 1 fb −1 of data.The gray band has a width of 1.5%, corresponding to the uncertainty in the ratio of the integrated reference luminosities from the lowPU to the one of highPU[33].

Table 1 :
Summary of the uncertainties in the number of delivered Z bosons in the 2017 highPU and lowPU data, and their ratio.The symbol δ denotes the relative uncertainty, i.e., δx = ∆x/x.The systematic and statistical uncertainties are added in quadrature to obtain the total uncertainty.

Table 2 :
Summary of cross checks performed by varying the length of the luminosity interval, the bin width of the m µµ histograms, and the range of the fit.As in Table1, the resulting variations of the number of Z bosons in the 2017 highPU and lowPU data, and their ratio, are shown.The δ denotes the relative variations, i.e., δx = ∆x/x.