Charmonium and $e^+e^-$ pair photoproduction at mid-rapidity in ultra-peripheral Pb-Pb collisions at $\sqrt{s_{\rm NN}}$=2.76 TeV

The ALICE Collaboration at the LHC has measured the J/$\psi$ and $\psi^{'}$ photoproduction at mid-rapidity in ultra-peripheral Pb-Pb collisions at $\sqrt{s_{\rm NN}}=2.76$ TeV. The charmonium is identified via its leptonic decay for events where the hadronic activity is required to be minimal. The analysis is based on an event sample corresponding to an integrated luminosity of about 23 $\mu\rm{b}^{-1}$. The cross section for coherent and incoherent J/$\psi$ production in the rapidity interval $-0.9<y<0.9$, are $\mathrm{d}\sigma_{{\rm J}/\psi}^{\rm coh}/\mathrm{d}y = 2.38^{+0.34}_{-0.24}\big(\mathrm{sta+sys}\big)$mb. and $\mathrm{d}\sigma_{{\rm J}/\psi}^{\rm inc}/\mathrm{d}y = 0.98^{+0.19}_{-0.17}\big(\mathrm{sta+sys}\big)$mb, respectively. The results are compared to theoretical models for J/$\psi$ production and the coherent cross section is found to be in good agreement with those models incorporating moderate nuclear gluon shadowing, such as EPS09 parametrization. In addition the cross section for the process $\gamma\gamma\rightarrow e^+ e^-$ has been measured and found to be in agreement with models implementing QED at leading order.


Introduction
field. The TPC acceptance covers the pseudorapidity region |η| < 0.9. Ionization measurements made along track clusters are used for particle identification [17]. Beyond the TPC, the Time-of-Flight detector (TOF) is a large cylindrical barrel of Multigap Resistive Plate Chambers (MRPCs) with about 150 000 readout channels, giving very high precision timing for tracks traversing it. Its pseudorapidity coverage matches that of the TPC. Used in combination with the tracking system, the TOF detector can be used for charged particle identification up to about 2.5 GeV/c (pions and kaons) and 4 GeV/c (protons). Still further out from the interaction region, the Electromagnetic Calorimeter (EMCAL) is a Pb-scintillator sampling calorimeter at a distance of ≈ 4.5 metres from the beam line, covering the range |η| ≤ 0.7 and a difference in azimuth ∆φ = 100 • . It has 20.1 radiation lengths and consists of 11 520 towers.
The analysis presented below also makes use of two forward detectors. The VZERO counters consist of two arrays of 32 scintillator tiles each, covering the range 2.8 < η < 5.1 (VZERO-A, on the opposite side of the muon arm) and -3.7 < η < -1.7 (VZERO-C, on the same side as the muon arm) and positioned respectively at z = 340 cm and z = -90 cm from the interaction point. The Forward Multiplicity Detector (FMD) consists of Si-strip sensors with a total of 51 240 active detection elements, arranged in five rings perpendicular to the beam direction, covering the pseudorapidity ranges −3.4 < η < −1.7 (FMD-3) and 1.7 < η < 5.1 (FMD-1 and FMD-2), a similar coverage to that of the VZERO detector. Finally, two sets of hadronic Zero Degree Calorimeters (ZDC) are located at 116 m on either side of the interaction point. The ZDCs detect neutrons emitted in the very forward region, such as neutrons produced by electromagnetic dissociation [18] (see Section 3).

Data analysis 3.1 Event selection
The present analysis is based on a sample of events collected during the 2011 Pb-Pb run, selected with a dedicated barrel ultra-peripheral collision trigger (BUPC), set up to select events containing two tracks in an otherwise empty detector. Events from two-photon production (γγ → µ + µ − , e + e − ) or from photonuclear vector meson production are selected by this trigger with the following characteristics: (i) at least two hits in the SPD detector; (ii) a number of fired pad-OR (N on ) in the TOF detector [19] in the range 2 ≤ N on ≤ 6 , with at least two of them with a difference in azimuth, ∆φ , in the range 150 • ≤ ∆φ ≤180 • ; (iii) no hits in the VZERO-A and no hits in the VZERO-C detectors. A total of about 6.5 ×10 6 events were selected by the BUPC trigger.
The integrated luminosity was measured using a trigger for the most central hadronic Pb-Pb collisions. The cross section for this process was obtained with a van der Meer scan [20], giving a cross section σ = 4.10 +0. 22 −0.13 (sys) b [21]. This gives an integrated luminosity for the BUPC trigger sample, corrected for trigger live time, of L int = 23.0 +0.7 −1.2 µb −1 . An alternative method based on using neutrons detected in the two ZDCs was also used. The ZDC trigger condition required a signal in at least one of the two calorimeters, thus selecting single electromagnetic dissociation (EMD) as well as hadronic interactions. The cross section for this trigger was also measured with a van der Meer scan, giving a cross section σ = 371.4 ± 0.6(sta) ± 24 19 (sys) b [18]. The integrated luminosity obtained for the BUPC by this method is consistent with the one quoted above within 3%.
The following selection criteria were applied in the data analysis: (i) a number of reconstructed tracks 1 ≤ N T RK ≤ 10, where a track is defined with loose criteria: more than 50% of findable clusters in the TPC fiducial volume and at least 20 TPC clusters, matching with those found in the ITS; (ii) a reconstructed primary vertex; (iii) only two good tracks passing tighter quality cuts: at least 70 TPC clusters, at least 1 SPD cluster and rejection of tracks with a kink. The tracks extrapolated to the reconstructed vertex should have a distance of closest approach (DCA) in the longitudinal beam direction DCA L ≤ 2 cm, and DCA T ≤0.0182+ 0.0350/p 1.01 T cm in the plane orthogonal to the beam direction, where p T is in (GeV/c); (iv) at least one of the two good tracks selected in (iii) with p T ≥ 1 GeV/c; this cut reduces the background while it does not affect the genuine leptons from J/ψ decay; (v) The VZERO trigger required no signal within a time window of 25 ns around the collision time in any of the scintillator tiles of both VZERO-A and VZERO-C. The size of the trigger windows are limited by the design of the VZERO front-end electronics which is operated at the frequency of the LHC clock, i.e. 40 MHz. In the offline analysis the event selection criteria consisted in an absence of a reconstructed signal in any of the VZERO scintillator tiles. The time windows in the offline analysis are enlarged to 40 ns and 60 ns around the collision time in VZERO-A and VZERO-C, respectively, and were chosen in order to maximize the vetoing efficiency; (vi) the dE/dx for the two tracks is compatible with that of electrons or muons; Figure 1 shows the TPC dE/dx of the positive lepton candidate as a function of the dE/dx of the negative lepton candidate, for J/ψ candidates in the invariant mass range 2.8 <M inv < 3.2 GeV/c 2 . It is worth noting that the TPC resolution does not allow to distinguish between muons and charged pions; (vii) the two tracks have same or opposite charges, depending on the analysis; (viii) invariant mass 2.2 < M inv < 6 GeV/c 2 .
The analysis of the γγ events is discussed in Section 5. In the remaining of this section we will focus on J/ψ analysis. The effect of the cuts on the statistics is listed in Table 1. In addition to the requirements (i) to (viii), a first sample enriched with coherent events was selected by applying a cut p T < 200 MeV/c for di-muons (p T < 300 MeV/c for di-electrons). Photoproduction of vector mesons can occur in interactions where additional photons are exchanged [22]. These additional photons can lead to break up of one or both nuclei. Since the energies of these photons are low, only a few neutrons are emitted when the nuclei break up. The exact upper limit on the number of emitted neutrons is not known, but in this analysis a cut on the neutron ZDC signal corresponding to less than 6 neutrons on each side has been applied. This cut reduces the statistics by 2.5%, which is considered as a source of systematic error +2.5% −0% . After applying all of these selections, 746 di-electron and 1301 di-muon coherent lepton-pair candidates remain. A second sample was enriched with incoherent events by applying a cut p T > 200 MeV/c on di-muons (p T > 300 MeV/c on di-electrons), giving 278 electron and 1748 muon incoherent event candidates.
As described in reference [4], during the 2011 Pb-Pb run the VZERO detector was optimized for the selection of hadronic Pb-Pb collisions, with a threshold corresponding to an energy deposit above that from a single minimum ionizing particle. Since the VZERO was used as a veto in the BUPC trigger, this setting could lead to an inefficiency in background rejection. In about 30% of the 2011 BUPC data taking sample, the FMDs were read out too. Since these detectors cover a pseudorapidity interval similar to that of the VZEROs, we have used, offline, their information to check for possible inefficiencies in the VZERO data. As expected, we found no hits in the FMD detector for the selected BUPC events, confirming that the VZERO inefficiency is very small.
A test on the electron and muon separation was applied to those tracks crossing the EMCAL. For each track we evaluated the ratio E/p between the energy released in the EMCAL and the reconstructed momentum; electrons lose their total energy in the shower generated in the EMCAL and for these a value E/p ≃ 1 is measured. Minimum ionizing particles lose only a small fraction of their energy in the EMCAL; in this case the measured E/p peaks in the region 0.1-0.2, in good agreement with the expectation. The e/µ separation was obtained by using two methods: (a) a sharp cut on Fig. 1, where all the particles beyond the dotted line are considered as electrons, and (b) using the average of the electron (muon) dE/dx and considering as electrons (muons) the particles within 3 sigma. The difference between the two methods was used as an estimate of the systematic error (see Table 2).

Acceptance and efficiency correction
The acceptance and efficiency of J/ψ reconstruction were calculated using a large sample of coherent and incoherent J/ψ events generated by STARLIGHT [23] and folded with the detector Monte Carlo simulation. STARLIGHT simulates photonuclear and two-photon interactions at hadron colliders. The simulations for exclusive vector meson production and two-photon interactions are based on the models in [9] and [24], respectively.
A separate simulation was performed for each run, in order to take into account the slight variations in run conditions during the data taking. The product of the acceptance and efficiency corrections (Acc × ε) J/ψ was calculated as the ratio of the number of the simulated events that satisfy all selections in Table 1 to the number of generated events with the J/ψ in the rapidity interval -0.9 <y< 0.9. In addition, for the coherent (incoherent) sample a reconstructed transverse momentum requirement p T < 200 MeV/c (p T > 200 MeV/c) was required for di-muons and p T < 300 MeV/c (p T > 300 MeV/c) for di-electrons. The average values for the combined acceptance and efficiency for J/ψ → e + e − (µ + µ − ) were found to be 2.71 (4.57)% and 1.8 (3.19)% for coherent and incoherent J/ψ, respectively. The STARLIGHT model predicts a dependence of the J/ψ cross section on the rapidity, giving a ≃10% variation over the rapidity range y=±0. 9. In order to evaluate the systematic error on the acceptance coming from the generator choice, we used a flat dependence of dσ J/ψ /dy in the interval -0.9 <y<0.9, as predicted by other models (see Fig. 5). The relative differences in (Acc × ε) between the methods were 2.5 (1.0)% for coherent electrons (muons), and 6.5 (3.5)% for incoherent electrons (muons), and are taken into account in the systematic error calculation. Transverse polarization is expected from helicity conservation for a quasireal photon. It is therefore assumed in these calculations that the J/ψ is transversely polarized, as found by previous experiments [25,26]. The trigger efficiency was measured relying on a data sample collected in a dedicated run triggered by the ZDCs only, as described at the beginning of this section. We analyzed events with only two tracks, showing the topology requested by the trigger selection. The resulting trigger efficiency was compared with that obtained by the Monte Carlo simulation, showing an agreement within +3.8% −9.0% . Figure 2 shows the invariant mass distribution for 2.2 < M inv < 6.0 GeV/c 2 for opposite-sign (OS) and like-sign (LS) electron and muon pairs. A J/ψ peak is clearly visible in the four spectra, on top of a continuum coming from γγ → e + e − (µ + µ − ) for the coherent enriched sample. The continuum for the incoherent enriched sample for the muon channel (bottom, left) is likely to come from misidentified π + π − pairs. To extract the J/ψ yield, the number of OS events in the interval 2.2 <M inv < 3.2 GeV/c 2 for electrons and 3.0 < M inv < 3.2 GeV/c 2 for muons were considered. In the mass intervals quoted above, 0 (3) LS electron(muon)-pairs were found for coherent enriched events, while 8 (53) LS pairs were found for incoherent enriched events. The corresponding number of OS pairs was 514 (365) for coherent enriched sample and 143 (178) for incoherent enriched events. The J/ψ yield was obtained by fitting the di-lepton invariant mass spectrum with an exponential function to describe the underlying continuum, and a Crystal Ball function [27] to extract the J/ψ signal. The Crystal Ball function tail parameters (α CB and n) were left free for the coherent enriched sample, giving a good agreement with those obtained by fitting the simulated data, and were fixed to values obtained from simulations for the incoherent enriched one. When fitting the incoherent di-muons a polynomial was also used to take into account the combinatorial background, constrained by the LS pair spectrum.

Coherent J/ψ cross section
The yield obtained for the coherent enriched sample (Fig. 2 top) is N yield = 265 ± 40(sta) ± 12(sys) for the J/ψ → e + e − channel and N yield = 291 ± 18(sta) ± 4(sys) for the J/ψ → µ + µ − channel. The systematic error on the yield is obtained by varying the bin size and by replacing the exponential with a polynomial to fit the γγ process. In addition, the Crystal Ball function parameters are obtained by fitting a simulated sample made of J/ψ and γγ event cocktail and then used to fit the coherent enriched data sample too. The difference in the yield obtained with the two Crystal Ball fit procedures is included in the systematic error. As a result we obtain a +7% −6% and +6% −5% systematic error on the signal extraction for coherent electrons and muons, respectively. For the coherent enriched sample, the central mass (standard deviation) value from the fit is 3.092 ± 0.036 GeV/c 2 (25.0±1.9 MeV/c 2 ) for the electron channel and 3.096 ± 0.002 GeV/c 2 (25±1.1 MeV/c 2 ) for the muon channel, in good agreement with the known value of the J/ψ mass and compatible with the absolute calibration accuracy of the barrel. For the incoherent sample, the central mass (standard deviation) value from the fit is 3.080 ± 0.007 GeV/c 2 (25±1.4 MeV/c 2 ) for the electron channel and 3.085 ± 0.007 GeV/c 2 (33±6 MeV/c 2 ) for the muon channel. The exponential slope parameter, λ γγ , of the continuum for the coherent enriched sample is computed at 2.2 < M inv < 2.6 GeV/c 2 (low) and 3.7 < M inv < 10 GeV/c 2 (high) for electrons with -0.9<η 1,2 <0.9, giving -1.55±0.88 GeV −1 c 2 and -0.73±0.18 GeV −1 c 2 , in good agreement with the corresponding Monte Carlo expectation, -1.07±0.16 GeV −1 c 2 and -0.81±0.01 GeV −1 c 2 , respectively. This is an additional indication that there is no important background in the invariant mass and p T range considered.
The fraction f D of coherent J/ψ mesons coming from the decay ψ ′ → J/ψ + anything, was estimated following the same prescription used in [4], with the theoretical estimates for f D ranging from 4.4% to 11.8% for electrons and 4.3% to 14.7% for muons. Alternatively, the ratio of coherent yields for ψ ′ to J/ψ can be extracted from the real data. Owing to the limited statistics, we combine the electron and muon channels to obtain N ψ ′ =17±10 and N ψ =505±48 (see Fig. 3). The fraction f D , for a given J/ψ polarization P, can be written as: The acceptance corrections are polarization dependent and give f P D ranging from 15±9% for longitudinal polarization to 11±6.5% for transverse polarization.
In what follows, we use the central value of theoretical and experimental estimates, and take the others as upper and lower limits, i.e. f D =0.10 +0.05 −0.06 . The di-electron (di-muon) p T distribution, integrated over 2.2 <M inv < 3.2 GeV/c 2 , (3.0 <M inv < 3.2 GeV/c 2 ) is shown in Fig. 4 right (left). The clear peak at low p T is mainly due to coherent interactions, while the tail extending out to 1 GeV/c comes from incoherent production. To estimate the fraction ( f I ) of incoherent over coherent events in the low p T region (p T < 300 MeV/c for di-electrons, p T < 200 MeV/c for di-muons), the ratio σ inc /σ coh , weighted by the detector acceptance and efficiency for the two processes, was calculated, giving f I =0.13 (0.06) for dielectrons (di-muons) when σ inc /σ coh was taken from STARLIGHT, and f I = 0.05 (0.03) when the model in [8] was used with leading twist contribution. For higher twist contributions the above model gives f I = 0.07 (0.03). An alternative method to extract an upper limit of f I from the data was carried out by fitting the measured p T distribution. Six different functions were used to describe the p T spectrum: (i) coherent J/ψ photoproduction; (ii) incoherent J/ψ photoproduction; (iii) J/ψ from coherent ψ ′ decay; (iv) J/ψ from incoherent ψ ′ decay; (v) two-photon production of continuum pairs; (vi) J/ψ produced in peripheral hadronic collisions.
The shapes for the first five fitting functions (Monte Carlo templates) were provided by STARLIGHT events folded with the detector simulation, while the last one is extracted from data at higher centralities [28]. The relative normalization was left free for coherent and incoherent photoproduction. The contribution from the ψ ′ was constrained from the estimate above ( f D =0.10 +0.05 −0.06 ), and the two-photon contribution was determined from the fit to the continuum in Fig. 2. The hadronic J/ψ were constrained by the fit to the region p T >1.1 GeV/c, where the ultra-peripheral J/ψ contribution is negligible. As a result of the fit, we obtain f I =0.044± 0.014 for di-muons and f I =0.15± 0.02 for di-electrons. Since these values are compatible, within the errors, with the theoretical expectations (for both models in the case of di-muons and for STARLIGHT for the di-electrons), they are used in the calculations. The fit reproduces properly the experimental data p T spectrum, clarifying the origin of the large p T J/ψs pointed out in the PHENIX paper [7]. Finally, the total number of coherent J/ψs is calculated from the yield extracted from the fit to the invariant mass distribution by resulting in N coh J/ψ (µ + µ − ) = 255 ± (sta) +14 −13 (sys) and N coh J/ψ (e + e − ) =212 ± 32(sta) +14 −13 (sys) respectively. The coherent J/ψ differential cross section is given by: where N coh J/ψ is the number of J/ψ candidates from Eq. 2 and (Acc × ε) J/ψ corresponds to the acceptance and efficiency as discussed above. BR(J/ψ → l + l − ) is the branching ratio for J/ψ decay into leptons [29], ∆y = 1.8 the rapidity interval bin size, and L int the total integrated luminosity. As a result we obtain dσ coh J/ψ /dy = 2.27 ± 0.14(sta) +0.30 −0.20 (sys) mb for the di-muon channel and dσ coh J/ψ /dy = 3.19 ± 0.50(sta) +0.45 −0.31 (sys) mb for the electron one. Since the di-electron and di-muon data are statistically separated samples, they can be combined; their weighted average gives dσ coh J/ψ /dy = 2.38 +0.34 −0.24 sta + sys mb. In addition, the fraction F n of coherent events with no neutron emission was estimated by STARLIGHT to be F n =0.68, while the model [8] predicts F n =0.76. Events with neutron emission can be efficiently tagged in ALICE by the ZDC calorimeters, taking advantage of their high efficiency (>98%). By fitting the di-electron (di-muon) invariant mass spectrum for events with and without neutron emission and with p T < 300 MeV/c (p T < 200 MeV/c), we obtain a fraction 0.70 ± 0.05(sta) in good agreement with the above estimates.

Incoherent J/ψ cross section
The incoherent cross sections are obtained in a similar way. For the incoherent enriched sample the obtained yield is (Fig.2 bottom), N yield = 61 ± 14(sta) +16 −7 (sys) for the J/ψ → e + e − channel and N yield = 91 ± 15(sta) +7 −5 (sys) for the J/ψ → µ + µ − channel. Here f D represents the fraction f D of coherent J/ψ mesons coming from the decay ψ ′ → J/ψ + anything, and was obtained only from formula (1), since the limited statistics did not allow the extraction of the ψ ′ yield from the data. The predictions for incoherent f D are calculated using both STARLIGHT and the model [8]. As a result we obtain a f D value ranging from 3.9% to 15.1% for muons and from 3.8% to 18.1% for electrons. By using the average we obtain f D = (9.5 ± 5.5)% for muons and f D = (11 ± 7)% for electrons. Using STARLIGHT, f C (the fraction of coherent J/ψ contaminating the incoherent sample), corrected by the acceptance and the efficiency, is found to be f C =0.5 for electrons and f C =0.02 for muons. By fitting the measured p T distribution (Fig. 4) we extract f C =(0.47 ± 0.09) for electrons, while f C is f C =(0.03 ± 0.03) for muons. These results are compatible with those from the models and will be used in the following. By applying the ratio 1/(1+ f D + f C ) to the N yield , the total number of incoherent muon events is N inc J/ψ (µ + µ − ) = 81 ± 13(sta) +8 −6 (sys), corresponding to dσ incoh J/ψ /dy = 1.03 ± 0.17(sta) +0.15 −0.12 (sys) mb for the di-muon channel. For electrons we obtain N inc J/ψ (e + e − ) = 39 ± 9(sta) +10 −5 (sys), corresponding to dσ incoh J/ψ /dy = 0.87 ± 0.20(sta) +0.26 −0.14 (sys) mb for the di-electron channel. Since these are statistically separate channels, their weighted average gives dσ incoh J/ψ /dy = 0.98 +0.19 −0.17 sta + sys mb.

Background and systematic error estimate
As discussed in [4], a possible loss of events might come from correlated QED pair production, i.e. interactions which produce both a J/ψ and a low mass e + e − pair (the latter process has a very large cross section), with one of the electrons hitting the VZERO detector and thus vetoing the event. This effect was studied in [4], with a control data sample where no veto at trigger level was applied. As a result, an upper limit on the inefficiency smaller than 2% was found. In the forward rapidity trigger only VZERO-A was used as a veto, and therefore we estimate, conservatively, a 4% systematic error for this study.
Another possible source of systematic error is the radiative decay J/ψ→ e + e − , neglected by the event generator used in this paper. We simulated a J/ψ → e + e − sample, where 15% of the events had a photon in the final state [30]. The Crystal Ball function fit applied to this sample provides fit parameters identical to those of the standard sample, and the (Acc × ε) is also not distinguishable from the standard value, so no correction is required in this analysis.
A possible background from hadronic J/ψ is found (Fig. 4) to be negligible for p T below around 200-300 MeV/c, and therefore it is not important for coherent production. For incoherent events this background was evaluated from the p T fit described above and gives a contribution (0.043±0.015) for di-electrons and (0.024±0.017) for di-muons. These fractions refer to events in the mass interval 2.2 <M inv < 3.2 GeV/c 2 for di-electrons and 3.0 <M inv < 3.2 GeV/c 2 for di-muons respectively, and therefore are not necessarily J/ψ only. We use these fractions as upper limits to be included in the systematic error, giving a contribution +0% −5% and +0% −3% , respectively. The hadronic combinatorial background can be estimated by LS events (see Table 2). It is negligible for coherent events and for incoherent di-electrons. For incoherent di-muons this background, possibly coming from misidentified pion pairs, was taken into account by using a polynomial function in the corresponding fit, as described at the end of Section 3.
Another source of background may come from photo-produced J/ψ by nuclei with impact parameters b<2R. According to a simulation, based on a calculation method similar to STARLIGHT, the cross section for this process (usually not included in the event generator) is 1.1 mb and 0.7 mb in the centrality bin (80-90)% and (90-100)%, respectively. The survival probability of the events in these two bins was simulated with 2.2 · 10 6 Pb-Pb minimum bias events produced by the HIJING event generator. Assuming the trigger conditions (i, ii, Section 3) and the analysis cuts (ii, iv, vi, Section 3) to be fully satisfied by di-leptons produced in UPC J/ψ decays, we find the fraction of events passing the trigger cut (iii) and the analysis cuts (iii, v) to be 0.06% and 0.3% in the two centrality bins. This process therefore gives a negligible contribution to the ultra-peripheral cross section.

Two-photon cross section
The STAR Collaboration measured the two-photon cross section with a precision of 22.5% when adding the statistical and systematic errors in quadrature [31]. This result was slightly larger than the one predicted by STARLIGHT, but within ∼ 2σ . The PHENIX Collaboration has also measured the cross section of two-photon production of di-electron pairs [7]. This measurement, which has un uncertainty of about 30%, when the statistical and the systematic errors are added in quadrature, was found in good agreement with STARLIGHT. The cross section for γγ → e + e − can be written in a similar way to Eq. 3, where N γγ was obtained by fitting the continuum in the invariant mass intervals 2.2 < M inv < 2.6 GeV/c 2 (N e + e − γγ = 186 ± 13(sta) ± 4(sys)) and 3.7 < M inv < 10 GeV/c 2 (N e + e − γγ = 93 ± 10(sta) ± 4(sys)), to avoid contamination from the J/ψ peak. In this analysis the integrated luminosity used was L int = 21.7 +0.7 −1.1 µb −1 and the cut (iv) on the track p T was removed. The cross section for the process γγ → µ + µ − was not studied due to a possible contamination (although small) from pions in the di-muons sample, suggested by the presence of LS events. The cross section for di-lepton invariant mass was computed between 2.2 < M inv < 2.6 GeV/c 2 and 3.7 < M inv < 10 GeV/c 2 , for a di-lepton rapidity in the interval -0.9 < y < 0.9, and requiring −0.9 < η 1,2 < 0.9 for each lepton. The data cuts applied to the Monte Carlo sample are the same as those applied in the analysis described above, resulting in a (Acc × ε) e + e − γγ =5.6% for 2.2 < M inv < 2.6 GeV/c 2 and (Acc× ε) e + e − γγ =4.73% for 3.7 < M inv < 10 GeV/c 2 . As a result we obtain σ e + e − γγ = 154 ± 11(sta) +16.6 −10.8 (sys) µb for the lower invariant mass interval and σ e + e − γγ = 91 ± 10(sta) +10.9 −8.0 (sys) µb for the higher invariant mass interval, to be compared with σ =128 µb and σ =77 µb given by STARLIGHT, respectively. In Fig. 6 the invariant mass distributions for 2.2 < M inv < 2.6 GeV/c 2 interval and for 3.7 < M inv < 10 GeV/c 2 are shown.

Discussion
The cross section of coherent J/ψ photoproduction is compared with calculations from six different models [8][9][10][11][12][13] in Fig. 5(a). The incoherent production cross section is compared with calculations by three different models [8,9,13]. These models calculate the photon spectrum in impact parameter space in order to exclude interactions where the nuclei interact hadronically. The differences between the models come mainly from the way the photonuclear interaction is treated. The predictions can be divided into three categories: i) those that include no nuclear effects (AB-MSTW08, see below for definition). In this approach, all nucleons contribute to the scattering, and the forward scattering differential cross section, dσ /dt at t = 0 (t is the momentum transfer from the target nucleus squared), scales with the number of nucleons squared, A 2 ; ii) models that use a Glauber approach to calculate the number of nucleons contributing to the scattering (STARLIGHT, GM, CSS and LM). The calculated cross section depends on the total J/ψ-nucleon cross section and on the nuclear geometry; iii) partonic models, where the cross section is proportional to the nuclear gluon distribution squared (AB-EPS08, AB-EPS09, AB-HKN07, and RSZ-LTA).
The rapidity region -0.9 <y< 0.9 considered here corresponds to photon-proton centre-of-mass energies, W γp , between 59 and 145 GeV. The corresponding range in Bjorken-x is between x = 5 × 10 −4 and x = 3 × 10 −3 . In this region, a rather strong shadowing is expected, and models based on perturbative QCD predict a lower value for the cross section than models using a Glauber approach to account for the nuclear effect.
The measured cross section, dσ coh J/ψ /dy = 2.38 +0.34 −0.24 sta + sys mb is in very good agreement with the calculation by Adeluyi and Bertulani using the EPS09 nuclear gluon prediction. The GM model, and the other models using a Glauber approach, predict a cross section a factor 1.5-2 larger than the data, overestimating the measured cross section by more than 3 standard deviations. So does the prediction based on the HKN07 parametrization, which includes less gluon shadowing than EPS09.
The model AB-EPS08, significantly underestimates the measured cross section by about a factor of two (about 5 standard deviations), indicating that the gluon shadowing is too strong in the EPS08 parameterization. The leading twist calculation (RSZ-LTA) is also significantly below the data, by about 2-3 sigma.
For the incoherent cross section, shown in fig. 5(b), there are three model predictions available, LM, STARLIGHT, and RSZ-LTA. The measured value deviates by about two standard deviations from the LM prediction, while STARLIGHT predicts an incoherent cross section 60% too high, and RSZ-LTA a factor 4 too low. Taking the ratio between the incoherent and coherent cross section provides further constraints on the treatment of the nuclear modifications implemented in the different models. Another advantage is that the photon spectrum is factorized out, so that the comparison directly probes the ratio of the photonuclear cross sections. The ratio obtained from data is 0.41 +0.10 −0.08 sta + sys . This can be compared with 0.21 from LM, 0.41 from STARLIGHT, and 0.17 from RSZ-LTA. Although the RSZ-LTA model is quite close for the coherent cross section at mid-rapidity, it seems to underpredict the incoherent cross section. The LM model also predicts a too low ratio. STARLIGHT, on the other hand, has about the right ratio of incoherent-to-coherent cross section, although it does not reproduce any of the cross sections individually. All three models use the Glauber model to calculate the incoherent cross section, but the implementation and the input cross section for γ + p → J/ψ + p varies. In STARLIGHT the scaling of the inelastic J/ψ + nucleus cross section, ranges from A 2/3 to A, depending on the J/ψ + nucleon cross section. In the first case, only the nucleons on the surface participate in the scattering, while in the second one all the nucleons contribute. The cross section for incoherent photoproduction is assumed in STARLIGHT to follow the same scaling, while in the other models, the reduction with respect to the A scaling is larger. The measured values for the γγ cross sections are 20% above but fully compatible within 1.0 and 1.5 sigma with the STARLIGHT prediction for the high and low invariant mass intervals, respectively, if the statistical and systematic errors are added in quadrature. This result provides important constraints on calculations that include terms of higher orders in α em . A reduction in the two-photon cross section of up to 30 % compared with leading-order calculations has been predicted [14]. The result for the two-photon cross section to di-lepton pairs, measured by ALICE with a precision of 12% and 16% for the low and high invariant mass range respectively, is thus fully consistent with STARLIGHT, and sets limits on the contribution from higher order terms [15]. This result reinforces the ALICE J/ψ photoproduction measurement in the forward rapidity region [4], where the cross section was based on σ γγ .

Summary
In summary, the first measurement of coherent and incoherent J/ψ photoproduction and two-photon production of di-lepton pairs in Pb-Pb collisions at mid-rapidity at the LHC has been presented and compared with model calculations. The J/ψ photoproduction cross sections provide a powerful tool to constrain the nuclear gluon shadowing in the region x ≈ 10 −3 . The coherent J/ψ cross section is found to be in good agreement with the model which incorporates the nuclear gluon shadowing according to the EPS09 parameterization (AB-EPS09).
Models which include no nuclear gluon shadowing are inconsistent with the measured results, as those which use the Glauber model to incorporate nuclear effects. The AB-HKN07 and AB-EPS08 models contain too little or too much shadowing, respectively, to match the data. Our results are about 3 sigma higher than the RSZ-LTA model prediction, although a deviation of just 1.5 sigma is found from the model upper limit. Nevertheless the above predictions may have large uncertainties coming not only from the parametrization of the nuclear gluon distribution but also from the selection of the hard scale, the contributions from the higher order terms and the treatment of the photon fluctuation to a quarkantiquark pair. The current measurement will contribute to resolve these uncertainties.
None of the three existing models predicts the incoherent photoproduction cross section correctly, but STARLIGHT predicts a correct incoherent-to-coherent ratio.
Finally, the measured two-photon cross section for di-electron production is consistent with the STARLIGHT model. This implies the models predicting a strong contribution of higher-order terms (not included in STARLIGHT) to the cross section are not favored.    Table 1, in the invariant mass range 2.2 < M inv < 6 GeV/c 2 . Coherent enriched sample (top) and incoherent enriched sample (bottom) for muons (left) and electrons (right). Blue (red) circles (triangles) are opposite-sign (like-sign) pairs. For like-sign pair the penultimate cut in Table 1 is replaced by the request of a same-sign pair. No LS events were found for coherent di-electron events.  Table 1, in the invariant mass interval 2.2 < M inv < 6 GeV/c 2 . Coherent dielectron and di-muon candidates are summed together. For like-sign pair the penultimate cut in Table 1 is replaced by the request of a same-sign pair data sum γ γ ψ coherent J/ ψ incoherent J/ ' decay ψ from ψ coherent J/ ' decay ψ from ψ incoherent J/ hadronic Fig. 4: Di-muon (left) and di-electron (right) p T distribution for ultra-peripheral Pb-Pb collisions at √ s NN = 2.76 TeV and -0.9 < y < 0.9 for events satisfying the event selection in the invariant mass interval 3.0 < M inv < 3.2 GeV/c 2 and 2.2 < M inv < 3.2 GeV/c 2 respectively, with the p T -range extended to p T < 1 GeV/c (top) and to p T < 5 GeV/c (bottom). The data points are fitted summing six different Monte Carlo templates: coherent J/ψ production (black), incoherent J/ψ production (red), J/ψs from coherent ψ ′ decay (light blue), J/ψs from incoherent ψ ′ decay (violet), γγ (green), and J/ψ produced in peripheral hadronic collisions (grey). The solid histogram (blue) is the sum.