Abstract
Exclusive \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) photoproduction is measured for the first time in ultraperipheral pPb collisions at \(\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 5.02\,\text {Te}\text {V} \) with the CMS detector. The cross section \(\sigma ({\upgamma _{}^{}} \mathrm{p}\rightarrow {{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\mathrm{p})\) is \(11.0 \pm 1.4\,\text {(stat)} \pm 1.0\,\text {(syst)} \) \(\mu \)b at \(\langle W_{{\upgamma _{}^{}} \mathrm{p}}\rangle = 92.6\,\text {Ge}\text {V} \) for photon–proton centre-of-mass energies \(W_{{\upgamma _{}^{}} \mathrm{p}}\) between 29 and \(213\,\text {Ge}\text {V} \). The differential cross section \(\mathrm {d}\sigma /\mathrm {d}|t |\) is measured in the interval \(0.025< |t | < 1\,\text {Ge}\text {V} ^{2}\) as a function of \(W_{{\upgamma _{}^{}} \mathrm{p}}\), where t is the squared four-momentum transfer at the proton vertex. The results are compared with previous measurements and theoretical predictions. The measured cross section \(\sigma ({\upgamma _{}^{}} \mathrm{p}\rightarrow {{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\mathrm{p})\) has a power-law dependence on the photon–proton centre-of-mass, consistent with electron–proton collision measurements performed at HERA. The \(W_{{\upgamma _{}^{}} \mathrm{p}}\) dependence of the exponential slope of the differential cross section \(\mathrm {d}\sigma /\mathrm {d}|t |\) is also measured.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
1 Introduction
Exclusive vector meson (VM) photoproduction, \({\upgamma _{}^{}} \mathrm{p}\rightarrow \mathrm {VM}\mathrm{p}\), has received renewed interest following recent studies of ultraperipheral collisions involving ions and protons at the CERN LHC [1, 2]. In such collisions, photon-induced interactions predominantly occur when the colliding hadrons are separated by a distance larger than the sum of their radii. In this case, one of the hadrons may emit a quasi-real photon that fluctuates into a quark-antiquark pair with the quantum numbers of the photon, which can then turn into a VM upon interacting with the other hadron. The interaction of the VM with the hadron proceeds via the exchange of the vacuum quantum numbers, the so-called pomeron exchange. Proton–lead (pPb) collisions are particularly interesting for studying photon–proton interactions [3, 4] because the large electric charge of the Pb nucleus strongly enhances photon emission. Also, in these events, one can determine the photon direction and hence the photon–proton centre-of-mass energy \(W_{{\upgamma _{}^{}} \mathrm{p}}\) unambiguously. This advantage is not present in symmetric colliding systems such as pp interactions. Exclusive VM photoproduction is interesting because the Fourier transform of the t distribution, with t being the squared four-momentum transfer at the proton vertex, is related to the two-dimensional spatial distribution of the struck partons in the plane transverse to the beam direction. Furthermore, some models suggest that the energy dependence of the integrated cross section and that of the t distribution may provide evidence of gluon saturation, as discussed in Refs. [5,6,7,8,9,10].
By using ultraperipheral pPb collisions at \(\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 5.02\,\text {Te}\text {V} \) at the LHC, the ALICE Collaboration has measured the exclusive photoproduction of \({{J/\uppsi _{}^{}}}{{\left( {1S}\right) }{}_{}^{}}\) mesons in the centre-of-mass energy interval \(20< W_{{\upgamma _{}^{}} \mathrm{p}} < 700\,\text {Ge}\text {V} \) [11, 12]. The LHCb Collaboration has studied exclusive \({{J/\uppsi _{}^{}}}{{\left( {1S}\right) }{}_{}^{}}\), \({{\uppsi _{}^{}}}{{\left( {2S}\right) }{}_{}^{}}\), and \(\Upsilon _{}^{}\) (nS) photoproduction in pp collisions at \(\sqrt{s} = 7\) and \(8\,\text {Te}\text {V} \) [13, 14]. Exclusive photoproduction of \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) mesons was first studied in fixed-target experiments at \(W_{{\upgamma _{}^{}} \mathrm{p}}\) values up to 20\(\,\text {Ge}\text {V}\) [15, 16]. Experiments at the HERA electron–proton collider at DESY have studied this process at \(W_{{\upgamma _{}^{}} \mathrm{p}}\) values ranging from 50 to 187\(\,\text {Ge}\text {V}\), both with quasi-real photons and for photons with larger virtualities [17, 18]. The HERA data have provided clear experimental evidence for the transition from the soft to the hard diffractive regime [19, 20]. More recently, exclusive photoproduction of \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) mesons has been studied by the STAR Collaboration in ultraperipheral AuAu collisions at the BNL RHIC collider [21,22,23], and by the ALICE Collaboration in PbPb collisions [24]. The cross sections measured by the ALICE and STAR Collaborations in photon-nucleus interactions are 40% lower than both the prediction from the Glauber approach and the corresponding measurements in photon–proton interactions [24, 25]. However, the Glauber approach reproduces the measured cross sections well at lower energies. This is an indication that nuclei do not behave as a collection of independent nucleons at high energies. In the present analysis, exclusive photoproduction of \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) mesons in the \({\uppi _{}^{}} _{}^{+}\) \({\uppi _{}^{}} _{}^{-}\) decay channel in ultraperipheral pPb collisions at \(\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 5.02\) \(\,\text {Te}\text {V}\) is measured. The cross section is measured as a function of \(W_{{\upgamma _{}^{}} \mathrm{p}}\) and t. In this paper \(|t |\) is defined as the squared transverse momentum of the \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^0\) meson, \(|t | \approx p_{\mathrm {T}} ^{2}\).
This paper is organized as follows. Section 2 describes the experimental apparatus and Sect. 3 the data and simulated Monte Carlo samples. The event selection procedure is illustrated in Sect. 4. Section 5 discusses the background contributions and Sect. 6 the strategy used to extract the signal; the systematic uncertainties are summarized in Sect. 7. The total and differential cross sections are presented in Sect. 8. The results are summarized in Sect. 9.
2 The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6\(\,\text {m}\) internal diameter, providing a magnetic field of 3.8\(\,\text {T}\). Within the solenoid volume are a silicon pixel and strip tracker, a lead tungsten crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections. The silicon tracker measures charged particles within the range \(|\eta |< 2.5\). It consists of 1440 silicon pixel and 15,148 silicon-strip detector modules and is located in the field of the superconducting solenoid. For nonisolated particles of \(1< p_{\mathrm {T}} < 10\,\text {Ge}\text {V} \) and \(|\eta | < 1.4\), the track resolutions are typically 1.5% in \(p_{\mathrm {T}}\) and 25–90 (45–150)\(\,\mu \text {m}\) in the transverse (longitudinal) direction [26].
The pseudorapidity coverage for the ECAL and HCAL detectors is \(|\eta |<3.0\). The ECAL provides coverage in the pseudorapidity range \(|\eta |<1.5\) in the barrel (EB) region and \(1.5<|\eta |<3.0\) in the two endcap (EE) regions. The HCAL provides coverage for \(|\eta |<1.3\) in the barrel (HB) region and \(1.3<|\eta |< 3.0\) in the two endcap (HE) regions. The hadron forward (HF) calorimeters (\(3.0<|\eta | < 5.2 \)) complement the coverage provided by the barrel and endcap detectors. The zero-degree calorimeters (ZDCs) are two Čerenkov calorimeters composed of alternating layers of tungsten and quartz fibers that cover the region \(|\eta |>8.3\). Both the HF and ZDC detectors are divided into two halves, one covering positive pseudorapidities, the other negative, and referred to as HF+ and ZDC+ (and HF- and ZDC-), respectively. Another calorimeter, CASTOR, also a Čerenkov sampling calorimeter, consists of quartz and tungsten plates and is located only at negative pseudorapidities with coverage of \(-6.6<\eta <-5.2\).
A more detailed description of the CMS detector, together with the definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [27].
3 Data and Monte Carlo simulation
This analysis uses data from pPb collisions at \(\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 5.02\,\text {Te}\text {V} \) collected with the CMS detector in February 2013. The beam energies are 4\(\,\text {Te}\text {V}\) for the protons and 1.58\(\,\text {Te}\text {V}\) per nucleon for the lead nuclei. The integrated luminosity is \({\mathcal {L}} = 7.4\,\mu \text {b}^{-1} \) for the pPb data set (protons circulating in the negative z direction) and \({\mathcal {L}} = 9.6\,\mu \text {b}^{-1} \) for the Pbp data set (protons circulating in the positive z direction). Since the events are asymmetric in rapidity, the pPb and Pbp samples are merged after changing the sign of the rapidity in the Pbp sample.
The starlight (version 2.2.0) Monte Carlo (MC) event generator [28] is used to simulate exclusive \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) photoproduction followed by the \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\rightarrow {{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) decay. The starlight generator models two-photon and photon-hadron interactions at ultrarelativistic energies. Two processes contribute to the exclusive \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) channel: resonant \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\rightarrow {\uppi _{}^{}} ^{+}{\uppi _{}^{}} ^{-}\) production, and nonresonant \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) production, including the interference term. Both processes are generated in order to calculate the signal acceptance and efficiency, and to extract the corrected signal yield. starlight is also used to generate exclusive \({{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}\) events. The pPb and Pbp samples are produced separately. The events are passed through a detailed Geant4 [29] simulation of the CMS detector in order to model the detector response, and are reconstructed with the same software used for the data.
4 Event selection
Table 1 presents the number of events after each selection requirement is applied. Events were selected online [30] by requiring the simultaneous presence of the two beams at the interaction point, as measured by the beam monitor timing system, in conjunction with at least one track in the pixel tracker. Offline, events are discarded if they have an energy deposit in any of the HF towers above the noise threshold of 3\(\,\text {Ge}\text {V}\). Events are also required to have exactly two tracks that pass the selection criteria defined in Ref. [31], and to be associated with a single vertex located within 15\(\,\text {cm}\) of the nominal interaction point along the beam direction. The pion mass is assigned to each track. In order to minimize the effect of the uncertainty in the low-\(p_{\mathrm {T}}\) track efficiency, one of the tracks should have a \(p_{\mathrm {T}}\) larger than 0.4\(\,\text {Ge}\text {V}\), and the other larger than 0.2\(\,\text {Ge}\text {V}\). Both tracks are selected in the interval \(|\eta | < 2.0\). The rapidity of the \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) system is required to be in the interval \(|y_{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}} |<2.0\). To reject the photoproduction of \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) mesons from \({\upgamma _{}^{}} \mathrm {Pb}\) interactions where the proton radiates a quasi-real photon, the \(p_{\mathrm {T}}\) of the \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) system is required to be larger than 0.15\(\,\text {Ge}\text {V}\) (as discussed in Sect. 5).
A sizable background contribution comes from proton dissociative events, \({\upgamma _{}^{}} \mathrm{p}\rightarrow {{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\mathrm{p}^{*}\), where \(\mathrm{p}^{*}\) indicates a low-mass hadronic state. In these events the scattered proton is excited and then dissociates. The \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) is measured in the central region, whereas the low-mass state usually escapes undetected. To suppress this contribution, events with activity above noise thresholds in the CASTOR, HE, HF, and ZDC detectors are rejected. The signal-to-noise ratio in \(\hbox {ZDC}^{+}\) is better than in \(\hbox {ZDC}^{-}\) because of differences in radiation damage to the two detectors. For this reason, the ZDC energy thresholds shown in Table 1 are asymmetric. CASTOR is used for only the pPb sample because of its location, as discussed in Sect. 2. The final selection requires the two tracks to have opposite charges. A total of 20,060 opposite-sign pair events and 1514 same-sign pair events are selected in this analysis.
5 Background
The main background sources are listed below.
-
Nonresonant \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) production. This contributes mainly through an interference term. It is included when fitting the invariant mass distribution, as discussed in Sect. 6.
-
Exclusive photoproduction of \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \) and \({{{\upphi _{}^{}}}{{\left( {1020}\right) }{}_{}^{}}} \) mesons. Contamination from the decay \({{{\upphi _{}^{}}}{{\left( {1020}\right) }{}_{}^{}}} \rightarrow \hbox {K}^{+}\hbox {K}^{-}\) is removed by assigning the kaon mass to the tracks and rejecting events with invariant mass values of the \(\hbox {K}^{+}\hbox {K}^{-}\) system larger than 1.04\(\,\text {Ge}\text {V}\). In addition, contamination is expected from the \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \rightarrow {{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} {{\uppi _{}^{}} _{}^{0}} \) and \({{{\upphi _{}^{}}}{{\left( {1020}\right) }{}_{}^{}}} \rightarrow {{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} {{\uppi _{}^{}} _{}^{0}} \) decays when the photons from the \({\uppi _{}^{}} _{}^{0}\) decay are undetected. Although the \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) invariant mass in these cases is mostly below the \({{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}\) mass, the rate of \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \) and \({{{\upphi _{}^{}}}{{\left( {1020}\right) }{}_{}^{}}} \) meson production increases with \(|t |\). As observed in this analysis and at HERA [32], undetected photons lead to an overestimate of the \(p_{\mathrm {T}}\) imbalance in the event, mimicking large \(|t |\) events. Since these processes cannot be modeled by starlight, their contribution is estimated from the fits of the unfolded invariant mass distributions described in Sect. 6. The \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \rightarrow {{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) amplitude is small, but is clearly visible through its interference with the \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\), which produces the small kink in the invariant mass spectrum near 800\(\,\text {Me}\text {V}\). This contribution is included in the invariant mass fit, as discussed in Sect. 6.
-
Exclusive photoproduction of \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) mesons.Footnote 1 The \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) decays mostly into a \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) meson and a pion pair, leading to final states with four charged pions, or with two charged pions and two neutral pions. The \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \rightarrow {{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} {{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) decay may also result in opposite-sign events when only two opposite-sign pions are detected because of the limited rapidity coverage of the detector. Such events will appear to have a \(p_{\mathrm {T}}\) imbalance, causing them to be incorrectly identified as large \(|t |~{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^0\) events, thereby resulting in a distortion of the \(|t |\) distribution. To validate the use of starlight for \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) photoproduction, exclusive \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} {{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) events are selected in the data. The data sample and the starlight simulation for \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) exclusive photoproduction are studied by applying the same selection criteria as for the \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\), except that on the number of tracks. Figure 1 shows a comparison of the \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) distributions of the reconstructed \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) mesons in the four pion event samples obtained from the data and the starlight simulation. All combinations of two oppositely charged pions are plotted in Fig. 1 if they have an invariant mass \(0.5< M_{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}} < 1.2\,\text {Ge}\text {V} \). In addition, the distribution of the same-sign events in the data is shown; they come mostly from \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) decays with two missing pions. Figure 1 shows that the data and the starlight results are in agreement, lending confidence to the performance of this generator. These distributions provide a template for the \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) distribution of the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) background used to estimate its contribution, as described in Sect. 6.
-
Proton dissociative \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) photoproduction. This contribution is suppressed by rejecting events with activity in the CASTOR, HE, HF, and ZDC detectors. In order to determine the residual contribution, a sample of dissociative events is selected by requiring activity in at least one of the forward detectors (CASTOR, HF, or ZDC). This sample provides a template for the \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) distribution of the dissociative events, under the assumption that the \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) distribution is independent of the mass of the dissociative system (the more forward the detector, the smaller the masses to which it is sensitive). Finally, this template is used to estimate the remaining dissociative background contributions, as discussed in Sect. 6.
-
Double pomeron exchange processes and photoproduction processes from \({\upgamma _{}^{}} \mathrm {Pb}\) interactions. Since the strong force has short range, only the nucleons on the surface of the nucleus may contribute to double pomeron exchange interactions; the corresponding cross section is therefore negligible [35]. For coherent processes in \({\upgamma _{}^{}} \mathrm {Pb}\) interactions, the size of the lead ion restricts the mean \(p_{\mathrm {T}}\) of the VM to be about \(60\,\text {Me}\text {V} \), corresponding to a de Broglie wavelength of the order of the nucleus size. Taking into account the detector resolution, all coherent \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) events have \(p_{\mathrm {T}}\) less than \(0.15\,\text {Ge}\text {V} \). Thus, events from \({\upgamma _{}^{}} \mathrm {Pb}\) interactions contribute to the lowest \(|t |\) region, which is not included in this analysis.
6 Signal extraction
The extraction of the signal is carried out in two steps. First, the proton dissociative and the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) contributions are estimated by performing a fit to the data as a function of \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\). This method relies on the fact that exclusive \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) events contribute mainly to the low-\(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) region (\(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}<0.7\,\text {Ge}\text {V} \)), whereas nonexclusive events dominate the high-\(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) region (\(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}>1.2\,\text {Ge}\text {V} \)), and the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) contribution is mostly at intermediate \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) values (\(0.7<p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}<1.2\,\text {Ge}\text {V} \)). This makes the identification of the proton dissociative and the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) contributions robust. Second, the yield of exclusive \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) candidates is extracted by performing a fit to the unfolded invariant mass distribution. Since the events from exclusive \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) production have a different invariant mass distribution from the signal events, they are subtracted before correcting the data for acceptance and efficiency. Conversely, the proton dissociative background has the same invariant mass and angular distributions as the signal, and its effect is corrected after unfolding by scaling the observed yields according to the fit performed in the first step.
To extract the normalizations of the proton dissociative and the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) backgrounds, an unbinned maximum likelihood fit is performed to the data as a function of \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) in the rapidity interval \(|y_{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}} | < 2\). The sum of the following distributions is fitted to the data at the reconstructed level: the signal distribution and the \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) continuum, as simulated by starlight, the distribution of the proton dissociative background, which is extracted from the data control sample, and the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) fitting template, which is simulated using starlight. The normalization of each of these components is determined from the fit. The signal \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) distribution generated by starlight is reweighted to describe the data using the theory-inspired expression \(\mathrm {e}^{- b|t |}\) [15]. The initial b value of starlight is \(12\,\text {Ge}\text {V} ^{-2}\) and the reweighted b is \(13.1^{+0.4}_{-0.3}\,\text {(stat)} \,\text {Ge}\text {V} ^{-2}\).
The result of the fit of the \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) distributions is shown in Fig. 2, including the systematic uncertainties associated with the fitting procedure that are discussed in Sect. 7. The resulting residual proton-dissociative and \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) contributions, over the whole rapidity interval, are \(18\pm 2\%\) \(\,\text {(stat)}\) and \(20\pm 2\%\,\text {(stat)} \), respectively. Similar fractions of proton dissociative and \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) contributions are obtained in the four rapidity intervals used in the differential cross section measurement as a function of rapidity. This is consistent with the small energy dependence of these processes in the energy range of this analysis. As seen in Fig. 2, the signal and both background contributions are of the same order of magnitude around \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}=1\,\text {Ge}\text {V} \), corresponding to a signal-to-background ratio of about 30%. For this reason, only the region \(|t | < 1\) \(\text {Ge}\text {V} ^{2}\) is used in this measurement.
The \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) background is subtracted in bins of invariant mass using the normalization obtained from the \(p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) fitting templates. The invariant mass distribution is then unfolded using the iterative D’Agostini method [36], which is regularized by four iterations. In particular, the Bayesian iterative unfolding technique is used, as implemented in the roounfold package [37]. This procedure leads to corrections for experimental effects including possible data migration between bins. The response matrix is obtained from starlight. The average of the combined acceptance and efficiency is 0.13 and is almost independent of both \(p_{\mathrm {T}}\) and \(\eta \), whereas it is sensitive to the invariant mass.
The invariant mass shape of the \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) in photoproduction deviates from that of a pure Breit–Wigner resonance [38]. Several parameterizations of the shape exist. One of the most often used is the Söding formula [39], where a continuum amplitude B is added to a Breit–Wigner distribution. Following the recent results by the STAR Collaboration [23] and the earlier ones by the DESY-MIT Collaboration [40], a further relativistic Breit–Wigner component is added to account for \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \) photoproduction, followed by the decay \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \rightarrow {{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \). This leads to the following fitting function:
Here A is the amplitude of the \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) Breit–Wigner function, B is the amplitude of the direct nonresonant \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) production, C is the amplitude of the \({{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}\) contribution, and the mass-dependent widths are given by
and
where \(\varGamma _{0}\) is the pole width for each meson and \(m_{{{\uppi _{}^{}} _{}^{\pm }}}\) is the charged pion mass. Since the branching fraction (\({\mathcal {B}}\)) for \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \rightarrow {{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) is small, only the first order term in the \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}}-{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) mass mixing theory is considered [40], leading to
with \({\mathcal {B}}({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \rightarrow {\uppi _{}^{}} {\uppi _{}^{}})=0.0153^{+0.0011}_{-0.0013}\) [33]. The H1 and ZEUS measurements did not include the \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}}-{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) interference component, although the ZEUS data seem to indicate its effect in the mass spectrum near 800\(\,\text {Me}\text {V}\) [17].
Figure 3 shows the fit of the unfolded distribution with the modifed Söding model. A least squares fit is performed for the interval \(0.6< M_{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}} < 1.1\,\text {Ge}\text {V} \), with the quantities \(M_{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}}\), \(M_{{{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}}}\), \(\varGamma _{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}}\), \(\varGamma _{{{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}}}\), A, B, C, and \(\phi _{{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \) treated as free parameters. This model includes the interference between resonant \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) and direct \({\uppi _{}^{}} ^{+}{\uppi _{}^{}} ^{-}\) production, as well as between \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) and \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \) production. To correct for the \({{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}\) reflection in the \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) mass spectrum, a Gaussian function peaking around 500\(\,\text {Me}\text {V}\) [18] is added as a further component of the invariant mass fit. This is only visible at high \(|t |\) values, as shown in Fig. 4. The fit yields \(M_{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}} = 773\pm 1\) \(\,\text {(stat)}\) \(\,\text {Me}\text {V}\) and \(\varGamma _{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}} = 148\pm 3\,\text {(stat)} \,\text {Me}\text {V} \), and \(M_{{{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}}} = 776\pm 2\) \(\,\text {(stat)}\) \(\,\text {Me}\text {V}\), consistent with the world average values [33]. The fitted value of the \({{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}\) width, \(\varGamma _{{{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}}} = 30\pm 5\) \(\,\text {(stat)}\) \(\,\text {Me}\text {V}\), is instead larger than the world average because of the detector resolution.
The \(|B/A |\) and C / A fractions are also determined; they measure the ratios of the nonresonant and \({{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}\) contributions to the resonant \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) production, respectively. Since the ZEUS Collaboration found that \(|B/A |\) decreases as \(|t |\) increases, the fit is repeated for \(|t |<0.5\) \(\text {Ge}\text {V} ^{2}\) resulting in \(0.50 \pm 0.06\) \(\,\text {(stat)}\) \(\text {Ge}\text {V} ^{-1/2}\). For this kinematic region H1 measured \(|B/A |=0.57 \pm 0.09\) \(\,\text {(stat)}\) \(\text {Ge}\text {V} ^{-1/2}\) and ZEUS \(|B/A |=0.70 \pm 0.04\) \(\,\text {(stat)}\) \(\text {Ge}\text {V} ^{-1/2}\). If the fit is repeated without the \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}}-{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) interference component, the result for \(|B/A |\) changes by less than its statistical uncertainty. The measured ratio of the \({{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}\) to \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) amplitudes is \(C/A=0.40\pm 0.06\) \(\,\text {(stat)}\), consistent with the prediction of starlight, \(C/A=0.32\), and the measurements of the STAR [23] and the DESY-MIT [40] experiments, which report \(C/A=0.36\pm 0.03\) \(\,\text {(stat)}\) and \(C/A=0.36\pm 0.04\) \(\,\text {(stat)}\), respectively. The present fit gives a nonzero \({{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}\) phase angle, \(\phi _{{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} =1.8\pm 0.3\) \(\,\text {(stat)}\), also in agreement with the previous measurements [23, 40].
Additionally, the fit is performed in \(|t |\) and y bins as shown in Fig. 4. To ensure fit stability, the \(M_{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}}\), \(M_{{{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}}}\), \(\varGamma _{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}}\), \(\varGamma _{{{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}}}\), \(\phi _{{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \) and \(|C/A |\) parameters are fixed to the values obtained for the full rapidity interval. The \({{{\upomega _{}^{}}}{{\left( {783}\right) }{}_{}^{}}} \rightarrow {{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} {{\uppi _{}^{}} _{}^{0}} \) contribution increases with \(|t |\), as reported by the H1 Collaboration [18] and as seen in Fig. 4. The \(|B/A |\) ratio is found to be independent of \(W_{{\upgamma _{}^{}} \mathrm{p}}\) and decreases with \(|t |\), in agreement with results reported by ZEUS [17].
7 Systematic uncertainties
The following sources of systematic uncertainty are considered.
Integrated luminosity determination: The uncertainty in the integrated luminosity is 4% for both the pPb and Pbp samples [41].
Track reconstruction: The contribution of the tracking efficiency to the systematic uncertainty is studied with the method described in Ref. [26], where the ratio of yields of neutral charm mesons decaying to two-body and four-body final states is compared with data and simulation for pion momenta above 300\(\,\text {Me}\text {V}\). The accuracy of the detector simulation to reproduce the single-pion tracking efficiency is 3.9%. For the present measurement, this yields a 7.8% uncertainty.
Unfolding: The uncertainty associated with the unfolding procedure is determined by modifying the number of iterations used for the Bayesian unfolding from the nominal value of 4 to 3 and 5. The resulting uncertainty is smaller than that found when changing the model for building the response matrix. The latter is estimated by comparing two different starlight samples: resonant \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) meson production, and combined resonant and nonresonant \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) production. The resulting effect on the integrated cross section is 3%.
Uncertainty in the photon flux: The uncertainty in the photon flux is 9% for the high-\(W_{{\upgamma _{}^{}} \mathrm{p}}\) data point and 2% at low \(W_{{\upgamma _{}^{}} \mathrm{p}}\), as discussed in Ref. [11]. The flux is computed in impact parameter space, convolved with the probability of no hadronic interactions. The radius of the lead nucleus is varied by the nuclear skin thickness (\(\pm 0.5\,\text {fm} \)). In addition, in the calculation of the photon flux, the \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) pole mass in Eq. (1) is replaced by the reconstructed \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) mass on an event-by-event basis. The effect of this variation is negligible.
Calorimeter exclusivity: The uncertainty related to the exclusivity requirements is evaluated by varying the calorimeter energy thresholds. Increasing (or decreasing) the energy scale of the HF calorimeter towers by 5% results in a 1.0% variation of the exclusive \({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}} \) yields. The CASTOR energy scale is varied by \(17\%\) [42], resulting in a difference of 1% in the extracted \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) yield. The variations of the energy thresholds for HE and ZDC within their respective energy scale uncertainties have a negligible effect.
Background estimation: The uncertainty in the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) subtraction is evaluated by varying the normalization of the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) contribution by 20% with respect to that obtained from the fit shown in Fig. 2. As mentioned in Sect. 5, the proton dissociative background template is obtained by requiring a signal in at least one of the forward detectors: HF, CASTOR, or ZDC. To calculate the systematic uncertainty related to the estimation of this background, the analysis is repeated five times and each time alternative combinations of forward detectors are used to obtain the proton dissociative template. The following variations are studied: (i) HF alone; (ii) CASTOR alone; (iii) ZDC alone; (iv) HF or CASTOR; (v) HF or ZDC. For each of these combinations the proton dissociative contributions are obtained in each \(|t |\) and rapidity bin. The maximum deviations from the nominal results are taken as conservative estimates of the systematic uncertainty. The resulting effect on the integrated exclusive \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) photoproduction cross section is smaller than 10%.
Model dependence: In order to assess the uncertainty due to the model used to fit the invariant mass distribution, the Ross–Stodolsky model [43] is used instead of the Söding model. The resulting cross section changes by up to 8%, depending on the rapidity and \(|t |\) interval studied. Another contribution to the model dependence uncertainty comes from the reweighting procedure of the starlight MC described in Sect. 6. This uncertainty is evaluated by varying the reweighting parameter b within its uncertainty; it is found to increase as a function of \(|t |\), and reaches 32% for the highest \(|t |\) bin. The second contribution turns out to be dominant for all the rapidity and \(|t |\) intervals studied. The uncertainty in the extrapolation to the region \(|t | < 0.025\,\text {Ge}\text {V} ^{2}\) is model dependent. We estimated this uncertainty by studying different fitting functions to the differential cross section measurements. In particular, we studied a dipole form [28], a pure exponential \(\mathrm {e}^{-bt}\), and a modified exponential \(\mathrm {e}^{-bt+ct^2}\). The difference between the two most extreme extrapolated values is used as an estimate of the model dependence uncertainty.
The values of the systematic uncertainties for all \(y_{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}}\) and \(|t |\) intervals are summarized in Table 2. The systematic uncertainties are added in quadrature for the integrated photoproduction cross section. For the differential cross section results, the systematic uncertainties in Table 2 are treated as correlated between bins.
8 Results
The differential cross section for exclusive photoproduction of \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) mesons is given by
where \(N^{\text {exc}}_{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}}\) is the corrected number of exclusive \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) events obtained from the fits described in Sect. 6 by integrating the resonant component in the interval \(0.28< M_{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}}<1.50\,\text {Ge}\text {V} \) (\(2M_{{{\uppi _{}^{}} _{}^{\pm }}}< M_{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}} < M_{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}} + 5\varGamma _{{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}}\)); \({\mathcal {B}}\) is the branching fraction, which equals about 0.99 for the \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\rightarrow {\uppi _{}^{}} ^{+}{\uppi _{}^{}} ^{-}\) decay [33], \(\varDelta y\) is the rapidity interval, and L is the integrated luminosity of the data sample. The cross section \(\mathrm {d}\sigma /\mathrm {d}y(\mathrm{p}\mathrm {Pb}\rightarrow \mathrm{p}\mathrm {Pb}{{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0})\) is related to the photon–proton cross section, \(\sigma ({\upgamma _{}^{}} \mathrm{p}\rightarrow {{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\mathrm{p}) \equiv \sigma (W_{{\upgamma _{}^{}} \mathrm{p}})\), through the photon flux, \(\mathrm {d}n/\mathrm {d}k\):
Here, k is the photon energy, which is determined from the \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) mass and rapidity, according to the formula
The average photon flux and the average centre-of-mass energy (\(\langle W_{{\upgamma _{}^{}} \mathrm{p}} \rangle \)) values in each rapidity interval are calculated using starlight.
The unfolded invariant mass distribution is studied in different \(|t |\) bins, and the extraction of the \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) photoproduction cross section is performed in each bin. In order to compare with the HERA results, the \(p_{\mathrm {T}}\)-related measurements are presented in terms of \(|t |\), which is approximated as \(|t | \approx (p_{\mathrm {T}} ^{{{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}})^{2}\). Figure 5 shows the differential cross sections as a function of \(|t |\), together with the unweighted starlight prediction, whose slope parameter is independent of \(W_{{\upgamma _{}^{}} \mathrm{p}}\). The starlight prediction is systematically higher than the data in the high-\(|t |\) region. This trend becomes more significant as \(W_{{\upgamma _{}^{}} \mathrm{p}}\) increases.
Figure 6 shows the differential cross section \(\mathrm {d}\sigma /\mathrm {d}|t |\) in the rapidity interval \(-1.2< y({{\uppi _{}^{}} _{}^{+}} {{\uppi _{}^{}} _{}^{-}}) < 0\) compared with the H1 and ZEUS results [17, 18] in a similar \(W_{{\upgamma _{}^{}} \mathrm{p}}\) range.
The differential cross section as a function of \(|t |\) is fitted with the form \(A\mathrm {e}^{-bt+ct{^2}}\) in the region \(0.025<|t |<0.5\) \(\text {Ge}\text {V} ^{2}\). For the integrated rapidity bin the fit gives \(b=9.2\pm 0.7\) \(\,\text {(stat)}\) \(\text {Ge}\text {V} ^{-2}\) and \(c=4.6\pm 1.6\) \(\,\text {(stat)}\) \(\text {Ge}\text {V} ^{-4}\). The resulting values of the slope b are shown in Fig. 7 as a function of \(W_{{\upgamma _{}^{}} \mathrm{p}}\), together with those measured by H1 and ZEUS [17, 18]. The values of the parameter c are found to be constant within the fit uncertainties. The Regge formula [44] \(b = b_{0} + 2\alpha '\ln (W_{{\upgamma _{}^{}} \mathrm{p}}/W_{0})^{2}\), which parametrizes the dependence of b on the collision energy, is fitted to the data using \(W_{0} = 92.6\) \(\,\text {Ge}\text {V}\), the average centre-of-mass energy of the present data. The fit to the CMS data alone gives a pomeron slope of \(\alpha ' = 0.28 \pm 0.11\,\text {(stat)} \pm 0.12\) \(\,\text {(syst)}\) \(\text {Ge}\text {V} ^{-2}\), consistent with the ZEUS [17] value and the Regge expectation of \(0.25~\text {Ge}\text {V} ^{-2}\).
The resulting photon–proton cross section, obtained for \(W_{{\upgamma _{}^{}} \mathrm{p}}\) between 29 and 213\(\,\text {Ge}\text {V}\) (\(\langle W_{{\upgamma _{}^{}} \mathrm{p}} \rangle = 92.6\) \(\,\text {Ge}\text {V}\)) is extrapolated to the range \(0<|t |<0.5~\text {Ge}\text {V} ^{2}\) using the exponential fits just discussed and the starlight predictions in order to allow direct comparison with previous experiments. The resulting value is \(\sigma = 11.0 \pm 1.4\) \(\,\text {(stat)}\) \(\pm 1.0\) \(\,\text {(syst)}\) \(\mu \)b. The photon–proton cross section values, \(\sigma ({\upgamma _{}^{}} \mathrm{p}\rightarrow {{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\mathrm{p})\), for all rapidity bins are presented in Table 3 and Fig. 8. Figure 8 also shows a compilation of fixed-target [45,46,47,48] and HERA results [17, 18]. The results of two fits are shown in Fig. 8. The dashed line indicates the result of a fit to all the plotted data with the formula \(\sigma = \alpha _{1} W^{\delta _{1}}_{{\upgamma _{}^{}} \mathrm{p}} + \alpha _{2} W^{\delta _{2}}_{{\upgamma _{}^{}} \mathrm{p}}\) (see e.g. [19, 20]). The fit describes the data well and yields the values \(\delta _{1} = -0.81 \pm 0.04\) \(\,\text {(stat)}\) \(\pm 0.09\) \(\,\text {(syst)}\), \(\delta _{2} = 0.36 \pm 0.07\) \(\,\text {(stat)}\) \(\pm 0.05\) \(\,\text {(syst)}\). The CMS and HERA data are also fitted with the function \(\sigma = \alpha W^{\delta }_{{\upgamma _{}^{}} \mathrm{p}}\) as shown in Fig. 8. The fit yields \(\delta = 0.24 \pm 0.13\) \(\,\text {(stat)}\) \(\pm 0.04\) \(\,\text {(syst)}\). Only statistical and uncorrelated systematic uncertainties are considered in these fits.
9 Summary
The CMS Collaboration has made the first measurement of exclusive \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) photoproduction off protons in ultraperipheral pPb collisions at \(\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 5.02\,\text {Te}\text {V} \). The cross section for this process is measured in the photon–proton centre-of-mass energy interval \(29< W_{{\upgamma _{}^{}} \mathrm{p}} < 213\,\text {Ge}\text {V} \). The results are consistent with those of the H1 and ZEUS Collaborations at HERA, indicating that ion–proton collisions can be used in the same way as electron–proton ones, with ions acting as a source of quasi-real photons. The combination of the present data and the earlier, lower energy results agrees with theory-inspired fits. The differential cross section \(\mathrm {d}\sigma /\mathrm {d}|t |\) for \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) photoproduction is measured as a function of \(W_{{\upgamma _{}^{}} \mathrm{p}}\). The starlight prediction is systematically higher than the data in the high-\(|t |\) region. This trend becomes more significant as \(W_{{\upgamma _{}^{}} \mathrm{p}}\) increases.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: Release and preservation of data used by the CMS Collaboration as the basis for publications is guided by the CMS policy as written in its document “CMS data preservation, re-use and open access policy” (https://cmsdocdb.cern.ch/cgi-bin/PublicDocDB/RetrieveFile?docid=6032&filename=CMSDataPolicyV1.2.pdf&version=2).]
Notes
The data on the photoproduction of excited \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) states in the four-pion decay channel are currently limited. A resonance structure with a broad invariant mass distribution around 1600\(\,\text {Me}\text {V}\) is reported in the literature. According to the Particle Data Group this resonance has two components: the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1450}\right) }{}_{}^{}}} \) and the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) [33]. The nature of these states is still under investigation. Recently, the STAR Collaboration reported a measurement of exclusive photoproduction of four charged pions [34]. Their data are consistent with the \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \) assuming that the peak is dominated by spin states with \(J^{PC} = 1^{--}\). In order to reproduce these data, starlight assumes a single resonance with a mass of 1540\(\,\text {Me}\text {V}\) and a width of 570\(\,\text {Me}\text {V}\) [28]. In the present paper, this state is referred to as \({{{\uprho _{}^{}} _{}^{}}{{\left( {1700}\right) }{}_{}^{}}} \).
References
A.J. Baltz, The physics of ultraperipheral collisions at the LHC. Phys. Rept. 458, 1 (2008). https://doi.org/10.1016/j.physrep.2007.12.001. arXiv:0706.3356
J.G. Contreras, J.D.Tapia Takaki, Ultra-peripheral heavy-ion collisions at the LHC. Int. J. Mod. Phys. A 30, 1542012 (2015). https://doi.org/10.1142/S0217751X15420129
L. Frankfurt, M. Strikman, M. Zhalov, Elastic and large t rapidity gap vector meson production in ultraperipheral proton–ion collisions. Phys. Lett. B 640, 162 (2006). https://doi.org/10.1016/j.physletb.2006.07.059. arXiv:hep-ph/0605160
V. Guzey, M. Zhalov, Rapidity and momentum transfer distributions of coherent \(\rm J/\psi \) photoproduction in ultraperipheral pPb collisions at the LHC. JHEP 02, 046 (2014). https://doi.org/10.1007/JHEP02(2014)046. arXiv:1307.6689
N. Armesto, A.H. Rezaeian, Exclusive vector meson production at high energies and gluon saturation. Phys. Rev. D 90, 054003 (2014). https://doi.org/10.1103/PhysRevD.90.054003. arXiv:1402.4831
T. Toll, T. Ullrich, Exclusive diffractive processes in electron–ion collisions. Phys. Rev. C 87, 024913 (2013). https://doi.org/10.1103/PhysRevC.87.024913. arXiv:1211.3048
S.P. Jones, A.D. Martin, M.G. Ryskin, T. Teubner, Probes of the small \(x\) gluon via exclusive \(\rm J/\psi \) and \(\varUpsilon \) production at HERA and the LHC. JHEP 11, 085 (2013). https://doi.org/10.1007/JHEP11(2013)085. arXiv:1307.7099
V.P. Goncalves, F.S. Navarra, D. Spiering, Exclusive \(\uprho (770)\) and \(J/\varPsi \) photoproduction in ultraperipheral \(pA\) collisions: predictions of the gluon saturation models for the momentum transfer distributions. Phys. Lett. B 791, 299 (2019). https://doi.org/10.1016/j.physletb.2019.03.007. arXiv:1811.09124
J. Cepila, J.G. Contreras, M. Krelina, J.D.Tapia Takaki, Mass dependence of vector meson photoproduction off protons and nuclei within the energy-dependent hot-spot model. Nucl. Phys. B 934, 330 (2018). https://doi.org/10.1016/j.nuclphysb.2018.07.010. arXiv:1804.05508
J. Cepila, J.G. Contreras, J.D.Tapia Takaki, Energy dependence of dissociative \({\rm J}/ \psi \) photoproduction as a signature of gluon saturation at the LHC. Phys. Lett. B 766, 186 (2017). https://doi.org/10.1016/j.physletb.2016.12.063. arXiv:1608.07559
ALICE Collaboration, Exclusive \({\rm J}/ \psi \) photoproduction off protons in ultraperipheral p–Pb collisions at \(\sqrt{s_{{\rm NN}}}=5.02~{\rm TeV}\). Phys. Rev. Lett. 113, 232504 (2014). https://doi.org/10.1103/PhysRevLett.113.232504. arXiv:1406.7819
ALICE Collaboration, Coherent \({\rm J}/ \psi \) photoproduction in ultra-peripheral Pb–Pb collisions at \(\sqrt{s_{NN}} = 2.76~{\rm TeV}\). Phys. Lett. B 718, 1273 (2013). https://doi.org/10.1016/j.physletb.2012.11.059. arXiv:1209.3715
LHCb Collaboration, Updated measurements of exclusive \({\rm J}/ \psi \) and \(\psi \)(2s) production cross-sections in pp collisions at \(\sqrt{s}=7~{\rm TeV} \). J. Phys. G 41, 055002 (2014). https://doi.org/10.1088/0954-3899/41/5/055002. arXiv:1401.3288
LHCb Collaboration, Measurement of the exclusive Y production cross-section in pp collisions at \(\sqrt{s}=7~{\rm TeV} ~{\rm and} ~8~{\rm TeV} \). JHEP 09, 084 (2015). https://doi.org/10.1007/JHEP09(2015)084. arXiv:1505.08139
T.H. Bauer, R.D. Spital, D.R. Yennie, F.M. Pipkin, The hadronic properties of the photon in high-energy interactions. Rev. Mod. Phys. 50, 261 (1978). https://doi.org/10.1103/RevModPhys.50.261. [Erratum: https://doi.org/10.1103/RevModPhys.51.407]
J.A. Crittenden, Exclusive production of neutral vector mesons at the electron–proton collider HERA. PhD thesis, Physikalisches Institut der Universitaet Bonn (1997). arXiv:hep-ex/9704009. DESY-97-068, BONN-IR-97-01
ZEUS Collaboration, Elastic and proton dissociative \(\rho ^0\) photoproduction at HERA. Eur. Phys. J. C 2, 247 (1998). https://doi.org/10.1007/s100529800834. arXiv:hep-ex/9712020
H1 Collaboration, Elastic photoproduction of \(\rho ^0\) mesons at HERA. Nucl. Phys. B 463, 3 (1996). https://doi.org/10.1016/0550-3213(96)00045-4. arXiv:hep-ex/9601004
P. Newman, M. Wing, The hadronic final state at HERA. Rev. Mod. Phys. 86, 1037 (2014). https://doi.org/10.1103/RevModPhys.86.1037. arXiv:1308.3368
L. Favart, M. Guidal, T. Horn, P. Kroll, Deeply virtual meson production on the nucleon. Eur. Phys. J. A 52, 158 (2016). https://doi.org/10.1140/epja/i2016-16158-2. arXiv:1511.04535
STAR Collaboration, Coherent \(\rho ^0\) production in ultraperipheral heavy ion collisions. Phys. Rev. Lett. 89, 272302 (2002). https://doi.org/10.1103/PhysRevLett.89.272302. arXiv:nucl-ex/0206004
STAR Collaboration, \(\rho ^0\) photoproduction in ultraperipheral relativistic heavy ion collisions at \(\sqrt{s_{{\rm NN}} = 200~{\rm GeV}}\). Phys. Rev. C 77, 034910 (2008). https://doi.org/10.1103/PhysRevC.77.034910. arXiv:0712.3320
STAR Collaboration, Coherent diffractive photoproduction of \(\rho ^0\) mesons on gold nuclei at \(200~\rm GeV\)/nucleon-pair at the relativistic heavy ion collider. Phys. Rev. C 96, 054904 (2017). https://doi.org/10.1103/PhysRevC.96.054904. arXiv:1702.07705
ALICE Collaboration, Coherent \(\rho ^{0}\) photoproduction in ultra-peripheral Pb–Pb collisions at \( \sqrt{s_{\rm NN}}=2.76~{\rm TeV} \). JHEP 09, 095 (2015). https://doi.org/10.1007/JHEP09(2015)095. arXiv:1503.09177
L. Frankfurt, V. Guzey, M. Strikman, M. Zhalov, Nuclear shadowing in photoproduction of \(\rho (770)^{0}\) mesons in ultraperipheral nucleus collisions at RHIC and the LHC. Phys. Lett. B 752, 51 (2016). https://doi.org/10.1016/j.physletb.2015.11.012. arXiv:1506.07150
CMS Collaboration, Description and performance of track and primary-vertex reconstruction with the CMS tracker. JINST 9, P10009 (2014). https://doi.org/10.1088/1748-0221/9/10/P10009. arXiv:1405.6569
CMS Collaboration, The CMS experiment at the CERN LHC. JINST 3, S08004 (2008). https://doi.org/10.1088/1748-0221/3/08/S08004
S.R. Klein et al., STARLIGHT: a Monte Carlo simulation program for ultra-peripheral collisions of relativistic ions. Comput. Phys. Commun. 212, 258 (2017). https://doi.org/10.1016/j.cpc.2016.10.016. arXiv:1607.03838
GEANT4 Collaboration, Geant4—a simulation toolkit. Nucl. Instrum. Methods A 506, 250 (2003). https://doi.org/10.1016/S0168-9002(03)01368-8
CMS Collaboration, The CMS trigger system. JINST 12, P01020 (2017). https://doi.org/10.1088/1748-0221/12/01/P01020. arXiv:1609.02366
CMS Collaboration, CMS tracking performance results from early LHC operation. Eur. Phys. J. C 70, 1165 (2010). https://doi.org/10.1140/epjc/s10052-010-1491-3. arXiv:1007.1988
H1 Collaboration, Diffractive electroproduction of \(\rho \) and \(\phi \) mesons at HERA. JHEP 05, 032 (2010). https://doi.org/10.1007/JHEP05(2010)032. arXiv:0910.5831
Particle Data Group, M. Tanabashi et al., Review of particle physics. Phys. Rev. D 98 030001 (2018). https://doi.org/10.1103/PhysRevD.98.030001
STAR Collaboration, Observation of \(\pi ^{+} \pi ^{-} \pi ^{+} \pi ^{-}\) photoproduction in ultraperipheral heavy-ion Collisions at STAR. Phys. Rev. C 81, 044901 (2010). https://doi.org/10.1103/PhysRevC.81.044901. arXiv:0912.0604
A.J. Schramm, D.H. Reeves, Production of \(\eta \) mesons in double Pomeron exchange. Phys. Rev. D 55, 7312 (1997). https://doi.org/10.1103/PhysRevD.55.7312. arXiv:hep-ph/9611330
G. D’Agostini, A multidimensional unfolding method based on Bayes theorem. Nucl. Instrum. Methods A 362, 487 (1995). https://doi.org/10.1016/0168-9002(95)00274-X
T. Adye, Unfolding algorithms and tests using RooUnfold. In Proceedings, PHYSTAT 2011 Workshop on Statistical Issues Related to Discovery Claims in Search Experiments and Unfolding, CERN, Geneva, Switzerland 17-20 January 2011, p. 313, CERN. CERN, Geneva, (2011). arXiv:1105.1160. https://doi.org/10.5170/CERN-2011-006.313
G. McClellan et al., Photoproduction of neutral \(\rho ^{0}\) mesons. Phys. Rev. D 4, 2683 (1971). https://doi.org/10.1103/PhysRevD.4.2683
P. Söding, On the apparent shift of the \(\rho \) meson mass in photoproduction. Phys. Lett. 19, 702 (1966). https://doi.org/10.1016/0031-9163(66)90451-3
H. Alvensleben et al., Precise determination of rho-omega interference parameters from photoproduction of vector mesons off nucleon and nuclei. Phys. Rev. Lett. 27, 888 (1971). https://doi.org/10.1103/PhysRevLett.27.888
CMS Collaboration, Luminosity calibration for the 2013 proton–lead and proton–proton data taking. CMS Physics Analysis Summary CMS-PAS-LUM-13-002 (2013)
CMS Collaboration, Measurement of the inclusive energy spectrum in the very forward direction in proton–proton collisions at \(\sqrt{s}=13~{\rm TeV}\). JHEP 08, 046 (2017). https://doi.org/10.1007/JHEP08(2017)046. arXiv:1701.08695
M.H. Ross, L. Stodolsky, Photon dissociation model for vector meson photoproduction. Phys. Rev. 149, 1172 (1966). https://doi.org/10.1103/PhysRev.149.1172
K.A. Goulianos, Diffractive interactions of hadrons at high energies. Phys. Rept. 101, 169 (1983). https://doi.org/10.1016/0370-1573(83)90010-8
E665 Collaboration, Diffractive production of \(\rho ^0(770)\) mesons in muon–proton interactions at 470-GeV. Z. Phys. C 74, 237 (1997). https://doi.org/10.1007/s002880050386
D.G. Cassel et al., Exclusive \(\rho ^0\), \(\omega \) and \(\phi \) electroproduction. Phys. Rev. D 24, 2787 (1981). https://doi.org/10.1103/PhysRevD.24.2787
CLAS Collaboration, Exclusive \(\rho ^{0}\) meson electroproduction from hydrogen at CLAS. Phys. Lett. B 605, 256 (2005). https://doi.org/10.1016/j.physletb.2004.11.019. arXiv:hep-ex/0408005
CLAS Collaboration, Exclusive \(\rho ^{0}\) electroproduction on the proton at CLAS. Eur. Phys. J. A 39, 5 (2009). https://doi.org/10.1140/epja/i2008-10683-5. arXiv:0807.3834
Acknowledgements
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centres and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMWFW and FWF (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC (China); COLCIENCIAS (Colombia); MSES and CSF (Croatia); RPF (Cyprus); SENESCYT (Ecuador); MoER, ERC IUT and ERDF (Estonia); Academy of Finland, MEC, and HIP (Finland); CEA and CNRS/IN2P3 (France); BMBF, DFG, and HGF (Germany); GSRT (Greece); OTKA and NIH (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); MSIP and NRF (Republic of Korea); LAS (Lithuania); MOE and UM (Malaysia); BUAP, CINVESTAV, CONACYT, LNS, SEP, and UASLP-FAI (Mexico); MBIE (New Zealand); PAEC (Pakistan); MSHE and NSC (Poland); FCT (Portugal); JINR (Dubna); MON, RosAtom, RAS and RFBR (Russia); MESTD (Serbia); SEIDI and CPAN (Spain); Swiss Funding Agencies (Switzerland); MST (Taipei); ThEPCenter, IPST, STAR and NSTDA (Thailand); TUBITAK and TAEK (Turkey); NASU and SFFR (Ukraine); STFC (UK); DOE and NSF (USA). Individuals have received support from the Marie-Curie programme and the European Research Council and Horizon 2020 Grant, contract No. 675440 (European Union); the Leventis Foundation; the A.P. Sloan Foundation; the Alexander von Humboldt Foundation; the Belgian Federal Science Policy Office; the Fonds pour la Formation à la Recherche dans l’Industrie et dans l’Agriculture (FRIA-Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); the F.R.S.-FNRS and FWO (Belgium) under the “Excellence of Science – EOS” – be.h project n. 30820817; the Ministry of Education, Youth and Sports (MEYS) of the Czech Republic; the Lendület (“Momentum”) Programme and the János Bolyai Research Scholarship of the Hungarian Academy of Sciences, the New National Excellence Program ÚNKP, the NKFIA research Grants 123842, 123959, 124845, 124850, and 125105 (Hungary); the Council of Science and Industrial Research, India; the HOMING PLUS programme of the Foundation for Polish Science, cofinanced from European Union, Regional Development Fund, the Mobility Plus programme of the Ministry of Science and Higher Education, the National Science Center (Poland), contracts Harmonia 2014/14/M/ST2/00428, Opus 2014/13/B/ST2/02543, 2014/15/B/ST2/03998, and 2015/19/B/ST2/02861, Sonata-bis 2012/07/E/ST2/01406; the National Priorities Research Program by Qatar National Research Fund; the Programa Estatal de Fomento de la Investigación Científica y Técnica de Excelencia María de Maeztu, Grant MDM-2015-0509 and the Programa Severo Ochoa del Principado de Asturias; the Thalis and Aristeia programmes cofinanced by EU-ESF and the Greek NSRF; the Rachadapisek Sompot Fund for Postdoctoral Fellowship, Chulalongkorn University and the Chulalongkorn Academic into Its 2nd Century Project Advancement Project (Thailand); the Welch Foundation, contract C-1845; and the Weston Havens Foundation (USA).
Author information
Authors and Affiliations
Consortia
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funded by SCOAP3
About this article
Cite this article
Sirunyan, A.M., Tumasyan, A., Adam, W. et al. Measurement of exclusive \({{{\uprho _{}^{}} _{}^{}}{{\left( {770}\right) }{}_{}^{}}} ^{0}\) photoproduction in ultraperipheral pPb collisions at \(\sqrt{\smash [b]{s_{_{\mathrm {NN}}}}} = 5.02\,\text {Te}\text {V} \). Eur. Phys. J. C 79, 702 (2019). https://doi.org/10.1140/epjc/s10052-019-7202-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-019-7202-9