The photon PDF in events with rapidity gaps

We consider photon-initiated events with large rapidity gaps in proton-proton collisions, where one or both protons may break up. We formulate a modified photon PDF that accounts for the specific experimental rapidity gap veto, and demonstrate how the soft survival probability for these gaps may be implemented consistently. Finally, we present some phenomenological results for the two-photon induced production of lepton and $W$ boson pairs.


Introduction
Photon-initiated processes at the LHC allow us to study γp and two-photon interactions at unprecedented collision energies, for a range of final states. In inclusive processes taking account of electroweak corrections is of increasing importance for precision phenomenology, and an essential ingredient in these is the introduction of a photon parton distribution function (PDF), where data such as the electroproduction of an isolated photon ep → eγX at HERA, and electroweak boson production at the LHC are sensitive to the size of the photon distribution (see [1,2,3] for studies by the global parton fitting groups).
In addition to the inclusive case, it also natural to consider photon-initiated exclusive and diffractive processes. The colour-singlet photon exchange can lead naturally to rapidity gaps in the final state, and in addition these modes offer some important and potentially unique advantages. For example, diffractive vector meson production provides a probe of the gluon PDF at low x and Q 2 , as well as possible gluon saturation effects, γγ → W + W − pair production provides a precise probe of potential anomalous gauge couplings [4,5,6], while the theoretically well understood case of lepton pair production, γγ → l + l − , is sensitive to the effect of soft proton interactions [7,8] as well as potentially being useful for luminosity calibration [9]. Moreover, there has recently been a renewal of interest in photon-induced processes in light of the excess of events at 750 GeV in the diphoton mass spectrum seen by ATLAS [10] and CMS [11]; if this is due to a new resonance with a sizeable coupling to photons, then the γγ-induced production mechanism may be dominant, see e.g. [12,13,14,15].
Such processes may be purely exclusive, that is with the outgoing protons remaining intact after the collision, or the interacting protons may dissociate. In the latter case, while there will be additional secondary particle production in some rapidity interval, the events can nonetheless have a diffractive topology, and thus an attractive way to select photon exchange events is to require a Large Rapidity Gap (LRG) between the centrally produced system (W + W − or l + l − pair, J/ψ or Υ, etc) and the forward outgoing secondaries. There are a range of measurement possibilities for such processes at the LHC, with a promising experimental programme underway [16]. The potential for rapidity gap vetoes to select events with a diffractive topology is in particular relevant at LHCb, for which the relatively low instantaneous luminosity and wide rapidity coverage allowed by the newly installed HER-SCHEL forward detectors [17] are highly favourable, while similar scintillation counters are also installed at ALICE [18] and CMS [17]. In addition, exclusive events may be selected at the LHC by tagging the outgoing intact protons using the approved AFP [19] and installed CT-PPS [20] forward proton spectrometers, associated with the ATLAS and CMS central detectors, respectively, see also [16,21].
In this paper we will consider γγ-induced reactions where the outgoing protons may dissociate (see [8] and references therein for discussion of the purely exclusive case), but with large rapidity gaps present between the produced object and the outgoing proton dissociation products. Provided the experimental rapidity veto region is large enough, the remaining contribution from non γγ-initiated processes (e.g. standard Drell-Yan production) will be small, and can be suppressed with further cuts and subtracted using MC simulation, see for example [7,22]. When considering these processes, there are two important effects that must be correctly accounted for. First, the secondaries produced during the DGLAP evolution of the photon PDF may populate the LRG. This means that in order to calculate the corresponding cross section we can not use the conventional inclusive PDF which describes the probability to find a photon in the proton, without any additional restrictions. Rather, we have to construct a new, modified, PDF where the evolution equation is supplemented by the condition that no s-channel partons are emitted in the LRG interval. Second, we have to include the probability that the gap will not be filled by secondaries produced by additional soft interactions of the colliding protons. This gap survival factor, S 2 , can be calculated within a given model of soft hadronic interactions, see e.g. [23,24]; although this introduces an element of model-dependence in the corresponding predictions, these can be fairly well constrained by the requirement that they give a satisfactory description of high energy proton-proton scattering data such as the differential elastic proton cross section, dσ el /dt, and the proton dissociation cross sections.
The outline of this paper is as follows. In Section 2 we describe how a rapidity gap veto may be accounted for in a modified photon PDF. In Section 3 we discuss the inclusion of the survival factor and demonstrate the effects this has on the γγ luminosity. In Section 4 we present numerical predictions for lepton and W boson pair production at √ s = 13 TeV. Finally, in Section 6 we conclude.

The modified photon distribution
The photon PDF is given in terms of an input term γ(x, Q 2 0 ) at the starting scale Q 0 , and a term due to photon emission from quarks during the DGLAP evolution up to the hard scale µ. Since the QED coupling α is very small it is sufficient to consider just the leading O(α) contribution to the evolution, while the appropriate splitting functions which allow the evolution to be evaluated at NLO in the strong coupling α S have recently been calculated in [25], and are included here 1 . The photon PDF is thus given by where γ(x, Q 2 0 ) is the input photon distribution at the scale Q 0 . This may be written in terms of a coherent component due to the elastic process, p → p + γ, and p → N * + γ excitation, see [26], as well as a component due to emission from the individual quarks within the proton (i.e. the direct analogue of perturbative emission in the QCD case). The P γq (z) and P γg (z) are the NLO (in α S ) splitting functions. At LO we have where the indices q and l denote the light quark and the lepton flavours respectively 2 , see [25] for the full NLO results. In fact, if we ignore the small corrections that the photon PDF will give to the evolution of the quark and gluons then the equation (1) for the DGLAP evolution 1 In general, to be consistent the NLO correction to the γγ → X matrix element should also be included. However for the production of colourless particles that we will consider in this paper, these are zero.
2 In [26] the lepton contribution to P γγ was mistakenly omitted.
of the photon PDF can be solved exactly, giving where the photon Sudakov factor corresponds to the probability for the photon PDF to evolve from scales Q 0 to µ without further branching; here P q(l)γ (z) is the γ to quark (lepton) splitting function at NLO in α s . At LO it is given by where N a = N c e 2 q for quarks and N a = e 2 l for leptons, while the factor of 1/2 in (7) is present to avoid double counting over the quark/anti-quarks (lepton/anti-leptons). The Sudakov factor is generated by resumming the term proportional to P γγ , due to virtual corrections to the photon propagator, which is a relatively small correction to the photon evolution. However this correction is not negligible, in particular for larger masses; we have S γ ∼ 0.97 − 0.93 for M X = 20 − 500 GeV.
As described above, the solution (5) is only exact if we neglect the dependence of the quark and gluon PDFs on the photon PDF, through the P qγ and P gγ terms in their evolution, respectively. These correspond to O(α 2 ) corrections to the photon evolution, and are therefore formally higher-order in α, so that they can be safely neglected. To confirm this expectation, we have compared (5) with the result of solving (1) numerically with the P qγ term included in the quark evolution, at LO in α S and only considering QED evolution (i.e. using the QECDS scheme [27] described below) for concreteness; the contribution from P gγ only enters at NLO in α S and so will be further suppressed. As expected, the difference is very small, and the results are found to coincide to within less than 0.1%. We have also confirmed this by using the APFEL evolution code [27], with the results with and without the P qγ term in the quark evolution coinciding to a very similar level, irrespective of the evolution scheme used.
The above equations correspond to the fully inclusive distribution, that is without any gap survival conditions. To include these, we note that as shown in (6) the photon PDF at a scale µ may be expressed as a sum of a term, γ in (x, µ 2 ), due to the input PDF, i.e. generated by coherent and incoherent photon emission up to the scale Q 0 , multiplied by the probability of no further emission up to the hard scale µ, and a second term, γ evol (x, µ 2 ), due purely to DGLAP emission from the quark/gluons, which is independent of the input photon PDF. For the coherent input component, there is naturally a large rapidity gap between y LRG y X y q y p emission | {z } Figure 1: Schematic diagram corresponding diffractive topology described in text, corresponding to case that quark of rapidity y q is emitted beyond the edge of a LRG region.
Due to strong q t ordering, the transverse momentum of the recoiled quark is given by −q t , that is equal and opposite to that of the final-state photon. The rapidity of this quark is We require that the quark be produced with rapidity greater than some y LRG , corresponding to the end of the experimentally defined gap 2 : in this case, it is convenient to work in terms of the rapidity interval, δ = y p − y LRG between the edge of the gap and outgoing proton in which the quark may be emitted, see Fig. 1. The condition y q > y LRG in this notations takes the form and thus to obtain the modified photon PDF, corresponding to the kinematics with a LRG present, we have to supplement the integrand in (1) by a Θ function which ensure that the condition (7) is satisfied. This gives 3 where q t = (1 − z)Q ′ 2 . Due to strong q t ordering all the previous partons emitted during the evolution will have larger rapidities, y > y q , and therefore evidently do not spoil the rapidity gap; this condition is therefore sufficient for a LRG to be present. the intact proton or the N * excitation and the central system. This is also the case for the incoherent input term: since the input value of Q 2 0 ∼ 1 GeV 2 is small and we have the kinematic requirement on the quark q t < Q 0 , this implies that the transverse momentum of the final-state quark produced in the incoherent emission is small, and the rapidity of the produced secondaries is large, that is, similar to the outgoing proton/N * rapidity in the coherent case 3 .
Next, we consider the second term in (1), that is due to the DGLAP evolution. At LO, this corresponds to the splitting of a quark with a fraction x/z of the proton momentum p to a photon with longitudinal momentum xp (where p = | p|) and squared transverse momentum and to a s-channel quark with longitudinal momentum Due to strong q t ordering, the transverse momentum of the recoiled quark is given by −q t , that is equal and opposite to that of the final-state photon. The rapidity of this quark is We require that the quark be produced with rapidity greater than some y LRG , corresponding to the end of the experimentally defined gap 4 : in this case, it is convenient to work in terms of the rapidity interval, δ = y p − y LRG between the edge of the gap and outgoing proton in which the quark may be emitted, see Fig. 1. The condition y q > y LRG in this notation takes the form and thus to obtain the modified photon PDF, corresponding to the kinematics with a LRG present, we have to simply supplement the integrand in (5) by a Θ function which ensures that the condition (12) is satisfied. This gives 5 where q t = (1 − z)Q 2 and in the final expression serves to define the δ-dependent evolution component γ evol (x, µ 2 ; δ). Due to strong q t ordering all the previous partons emitted during the evolution will have larger rapidities, y > y q , and therefore evidently do not spoil the rapidity gap; this condition is therefore sufficient for a LRG to be present. We now consider some numerical results. As described above, for the input photon PDF, following [26] we include a coherent component due to purely elastic photon emission and an incoherent component due to emission from the individual quark lines, such that with where q t is the transverse momentum of the emitted photon, and Q 2 is the modulus of the photon virtuality, given by The functions F E and F M are the usual proton electric and magnetic form factors with 5 If we consider the evolution equation (5) in terms of the scale q t then the limit (12) corresponds to a straightforward upper limit on the momentum fraction, z. in the dipole approximation, where G E and G M are the 'Sachs' form factors. The incoherent input term is given by 6 accounts for the fact that the lowest proton excitation is the ∆-isobar, and the final factor (1 − G 2 E (Q 2 )) corresponds to the probability to have no intact proton in the final state (which is already included in the coherent component). Here m q = m d (m u ) when convoluted with d 0 (u 0 ), and the current quark masses are taken. As the quark distributions are frozen for Q < Q 0 , this represents an upper bound on the incoherent contribution. Although other models for this incoherent component may also be taken, the conclusions which follow are relatively insensitive to the specific choice, and so for simplicity we will not consider them here. We also note that it is possible to account explicitly for the first ∆-isobar excitation in the coherent component, see [26], however this does not have a noticeable effect on the results which follow, and is omitted here.
In Fig. 2 we show the effect of including the rapidity gap constraint (12) on the photon PDF, for two choice of scale and for different values of δ. Here, and in all numerical results which follow, we for concreteness use MMHT2014 NLO PDFs [28] for the quark term in (13). The suppression in the PDFs relative to the inclusive case, which becomes stronger as δ decreases, is clear. In addition, we can see that the suppression is stronger at lower x and higher µ 2 , as expected from (12): in the former case, the outgoing quark in the xγ(x, µ 2 = 10 4 GeV 2 ) , NNPDF2.3QED x Inclusive δ = 2 δ = 5 δ = 7 Figure 3: As in Fig. 2, but with the NNPDF2.3QED [2] set taken for the input PDF at Q 2 0 = 2 GeV 2 . The 68% confidence error bands are shown in the inclusive case.
q → qγ splitting has on average lower longitudinal momentum, while in the latter the quark transverse momentum is higher, such that in both cases the quark tends to be produced more centrally. These effects are not limited to the particular approach to modelling the photon PDF described above: in Fig. 3 we show the same plots as before, but using the NNPDF2.3QED [2] photon PDF as input 7 . The increased suppression with decreasing x and increasing Q 2 is again clear, with the resulting PDFs generally lying outside the uncertainty band in the inclusive PDF. We end this section with some comments. First, we note that qualitatively speaking the inclusion of the Θ function in the integral (13) plays the role of the Sudakov factor in gluonmediated central exclusive production (CEP) processes, see e.g. [29], that is, it accounts for the probability for no secondary partons emission. In the case of pure CEP processes, such emission is entirely forbidden, whereas here we only require that no secondaries are emitted into the veto region. Second, in accounting for the veto condition (12) in the case of the NLO splitting functions we should consider vetoes on the two emitted partons individually, i.e. qg(qq) for P γq (P γg ). However since the effect of the NLO correction is rather small (∼ 5% ) here we for simplicity use the same veto as in the LO case. This corresponds to a veto on the kinematics of the parton pair and so only gives an approximate indication of the effect to the NLO contribution. In addition, we emphasise that (13) corresponds to the survival of the LRG in terms of the secondary partons only. A complete evaluation, in which the probability that no secondary hadrons spoil the gap would require a Monte Carlo simulation which accounts for the fluctuations during the fragmentation and hadronization 7 We note that the PDF evolution for the NNPDF set is performed in the so-called QECDS scheme [27], where the QED and QCD factorization scales are treated seperately, with the QED evolution performed first. In the context of our approach, this corresponds to evaluating the quark PDFs in (1) at fixed scale Q 2 0 (we treat the QED evolution here at LO in α s , consistently with NNPDF, and hence no gluon term is present). However, as discussed in [27], this QECDS scheme leads to potentially large unresummed logarithms at higher scales, and we use it here and in all NNPDF results which follow only for the sake of comparison.
process. However, the results presented above should not be too sensitive to these effects, in particular if δ is large enough. Finally, we note that the photon PDF is formally defined for fully inclusive processes, i.e. where a complete sum over all final states is performed. For the distribution in the semi-exclusive case, only a subset of final states is summed over, and so the PDF no longer has this formal definition. Moreover, as discussed in the following section, factorization is explicitly violated by soft survival effects. Nonetheless, we may consider these PDFs as phenomenological objects which capture the important physics of the semi-exclusive process, i.e. the constraint (12), for which the Q 2 -evolution can be described within the same leading logarithmic approximation as in the standard DGLAP approach.

Soft gap survival factor
In addition to accounting for the rapidity gap veto in the q → qγ splitting associated with the DGLAP evolution of the photon PDF, we must also consider the possibility of additional soft proton-proton rescattering, that is underlying event activity, which can fill the rapidity gaps with secondary particles. The probability that this does not occur is encoded in the eikonal survival factor, S 2 elk , see [23,24] for some more recent theoretical work, and [8,29] for further discussion and references 8 .
Here we follow the approach described in detail in [8], where it is emphasised that the impact of survival effects depends sensitively on the subprocess, through the specific proton impact parameter dependence. For example, for exclusive photon exchange processes, the low virtuality (and hence transverse momentum) of the quasi-real photon exchange, which corresponds to relatively large impact parameters between the colliding protons, leads to an average survival factor that is close to unity, while for the less peripheral QCD-induced exclusive processes the suppression is much larger. This fact has important consequences in the current situation. In particular, the expression (13) for the photon PDF evaluated at a scale µ 2 is given in terms of an input distribution corresponding to low photon q 2 ⊥ < Q 2 0 , and a term generated by DGLAP evolution, for which we have q 2 ⊥ Q 2 0 . This larger scale (and hence smaller average impact parameter) suggests that the survival factor in the latter case will be much smaller.
To demonstrate this, we must determine precisely how survival effects are to be included in this semi-exclusive case. First, we recall how such corrections are included in the exclusive case, relevant for the input component of the photon PDF. Here, we can work at the amplitude level, as described in detail in [8], where purely exclusive two-photon initiated processes are considered. Survival effects are generated by the so-called 'screened' amplitude shown in Fig. 4, for the representative case of lepton pair production, where the grey oval represents the exchange of a pomeron with transverse momentum k t transferred through the loop. The where q 1t = q 1t +k t and q 2t = q 1t −k t are the incoming photon transverse momenta, as show in Fig. 4 (b), and other kinematic arguments of the amplitudes are omitted for simplicity. Here, T (q 1t , q 2t ) is the production amplitude; for q it = q it this corresponds to Fig. 4 (a). While the equivalent photon approximation and the corresponding expression (16) for the coherent component of the photon PDF are defined at the cross section level, as discussed in [8] the coherent amplitude may be unambiguously defined for the term proportional to the proton electromagnetic form factor, and thus included inside the k t integral (22). After adding to the 'bare' amplitude (i.e. without survival effects) and squaring, the average survival factor may be evaluated using where for illustration we have considered only the simplest, so-called 'one-channel' approach, which ignores any internal structure of the proton: see [31,32] for discussion of how this can be generalized to the more realistic 'multi-channel' case. As discussed in [8], the inclusion of survival effects requires a non-trivial vector combination of the incoming photon transverse momenta q it , which only after squaring and angular averaging allows the q it dependence to be factorized as in (16). For example, for the production of a spin-0 object the production amplitude in (22) should be decomposed as where the T λ 1 λ 2 are now the γ(λ 1 )γ(λ 2 ) → X helicity amplitudes and n 0 is a unit vector along the (beam) z-axis; see [8] for the full case, including the |J z | = 2 amplitudes. Using this, it can readily be shown that the bare amplitude squared reduces to the correct cross section level expressions corresponding to (16) or (20), while in the screened case it is again crucial to include this correct vector form of the amplitude, i.e. with the q it included inside the integral (22). We therefore do this here, following the approach of [8], for both the coherent and incoherent contributions. As the relative contributions of the amplitudes T λ 1 λ 2 affect the q ⊥ dependence in (24), the survival factor will in general depend on the helicity structure of the γγ → X process. In all numerical calculations that follow we consider for simplicity the case of a scalar object (T ++ = T −− , T ±∓ = 0). For other quantum numbers, the survival factor can vary by up to ∼ 10%. Next, we must consider the component of the photon PDF due to the DGLAP evolution term in (13), generated by the q → qγ splitting. We will first consider the case that the photon PDFs from both protons correspond to this evolution component, before discussing the mixed case. As mentioned above, in contrast to the input component, for which q 2 t < Q 2 0 , we now have q 2 t Q 2 0 , due to strong DGLAP q t ordering. The transverse momentum transferred through the pomeron loop corresponds to a soft physics scale, and is generally very low, being set by k 2 t ∼ 2/B el , where B el is the t-slope for elastic pp scattering, with at the LHC B el ∼ 20 GeV −2 [33] and so k 2 t ∼ 0.1 GeV 2 . Thus, when we consider a screened diagram of the type shown in Fig. 5 (top right), we have q it k t , and so when considering the screened amplitude (22), the incoming photon transverse momenta are simply given by q it ∼ q it , and the production amplitude T (q 1t , q 2t ) factorizes from the k t integral.
More precisely, we note that in fact the evolution component can no longer be defined at the amplitude level, as it is the cross section which is written in terms of the photon PDF. This is shown schematically in Fig. 5 (bottom left) for the unscreened case. We must therefore define screening corrections at the cross section level; this can be achieved by observing that for the screened contribution to the cross section of the type shown in Fig. 5 (bottom left), instead of dealing with the usual 'diagonal' photon PDF, we should instead consider the generalised PDF (GPDF), see [34,35] for reviews and references, that is where the skewedness is provided by the non-zero squared transverse momentum t ≈ −k t transferred through the t-channel exchange (i.e. the hard process and the DGLAP ladder generating the evolution of the GPDF) by the screening pomeron 9 . We can then to good approximation neglect this k t dependence in the evolution of the GPDF: this is well justified at LO as due to the strong q t ordering of the DGLAP evolution, the k t Q 0 transferred through the pomeron loop can be neglected in every rung of the evolution ladder for H γ except that nearest the proton. This allows us to write the above expression in the factorized form where γ evol is defined in (14) and F 1 is the proton Dirac form factor; this choice is motivated by the sum rule derived in [36] for the quark GPDF, which is appropriate here as it is the quark PDF which is driving the photon evolution. This allows the k t dependence induced by the screening corrections to be properly accounted for at the cross section level, namely by expanding the squared amplitude in the numerator of (23), giving where overall factors have been omitted for simplicity. This expansion is shown schematically in Fig. 5 (bottom). In fact, the calculation may be more easily performed in impact parameter , for different systems of mass M X produced with Y X = 0. The coherent, incoherent and evolution components are shown for (proton 1, proton 2).
space, in which case we simply have where b t is the impact parameter separation of the two protons, F 1 (b t ) is the Fourier transform of the Dirac form factor, and Ω(b t ) is the so-called proton opacity; physically, exp (−Ω(b t )) represents the probability of no inelastic scattering at impact parameter b t . Here, as above, we work in the single-channel approximation for the sake of clarity, but in actual calculations we use the two-channel approach described in [37]. The above discussion only applies for the case that the evolution components are probed from both protons, however we must also consider the mixed case, where we probe the evolution component from one proton and the input from the other. Such a process will induce a distinct proton impact parameter distribution from the cases above, leading to a different survival factor. The precise size of this can be readily calculated by adjusting (27) to include an integral over the photon q t emitted from the input side. However here care must again be taken as in (24) to include the correct vector dependence on the incoming photon q t on the input side. As the transverse momentum q 2t (say) on the evolution side is much larger than the k t transferred through the exchanged pomeron, the q 2t dependence factorizes from the k t integral in (22), and q 2t can be averaged over at the cross section level, leaving only a dependence on q 1t in the screened amplitude. We find that this averaging washes out the dependence (24) seen in the purely coherent/incoherent contributions, such that the survival factor here does not in fact depend on the quantum numbers of the produced object.
We now consider some numerical results for the average survival factors. To calculate these, we use model 4 of [37], which applies a two-channel eikonal approach, in which the incoming proton is considered to be a coherent superposition of two diffractive 'Good-Walker' (GW) eigenstates [38], each of which may scatter elastically. In this situation, we have some freedom as to how the expression (26) may be suitably generalised, that is how the photon GPDF couples to the individual GW eigenstates. By default, we assume that these are the same as the coupling to the pomeron taken in [37], that is H γ i ∼ γ i · F 1 (t), for eigenstate i = 1, 2, where γ i is coupling to the pomeron. However, this is not the only possibility: for larger x where the quark contribution to H γ is more important, it may be more sensible to instead assume that this coupling is universal, i.e. simply H γ i ∼ F 1 (t). A further question is whether the proton form factor F 1 is the appropriate choice: it may be be more suitable, in particular at low x, to take the same form factors as in [37] for the coupling of the pomeron to the GW eigenstates. In fact, it turns out that these different choices generally have a small effect on the observable predictions; we will comment on this further below.
The corresponding average survival factors for all combinations of photon PDF components from each proton are given in Table 1. A large range of expected suppression factors is evident, with as anticipated S 2 for the lower scale (and hence more peripheral) coherent production process being higher than for the higher scale evolution component. The survival factor for the incoherent component of the input PDF is seen to be particularly small: this is due to the (1 − G 2 E (t)) factor in (20), which accounts for probability to have no intact proton in the final state, and is therefore peaked towards larger t, i.e. less peripheral interactions, where it is less likely to produce an intact proton.
These results have important implications for the standard factorisation formula as, using (14) and (15), we have Crucially, each of these six independent contributions now has a distinct (and in principle, x and/or µ dependent) survival factor associated with it, and therefore the simple factorisation implied by (29), where the photon flux associated with each proton i can be factorised in terms on an independent PDF γ(x i , µ 2 ), no longer holds; instead, this now depends on the state of the other interacting proton, through the influence this has on the survival factor. Physically, this is to be expected, as the survival factor is generated by additional soft proton-proton interactions, which then prevent all of the physics associated with the initial-state photon produced by a given proton being considered independently from the other proton. Such factorisation breaking effects have already been seen in, for example, diffractive production at HERA, where the predictions using the so-called diffractive PDFs are known to dramatically overshoot the data when naively applied in hadron-hadron collisions [39]. It is therefore not possible to show equivalent plots to Fig. 2 demonstrating the impact of survival effects on the individual photon PDF. Instead, we can consider the γγ luminosity, given by dL dM 2 where τ = M 2 X /s and we take µ 2 = M 2 X as the scale of the PDFs. Such a variable also gives a clearer picture of the suppression we can expect in physical cross sections, as comparing with (29) we have The photon-photon luminosity at √ s = 13 TeV is shown in Fig. 6 (left), for inclusive and semi-exclusive (δ = 5 for both protons) production. For demonstration we show the latter case both with and without survival effects included. We can see that the inclusion of the condition (12) leads to a factor of ∼ 2 reduction in the luminosity, roughly consistently with Fig. 2, while the inclusion of survival effects leads to a further suppression of a similar size. That the suppression due to both effects is similar in size is not necessarily to be expected, and indeed for different choices of δ and/or √ s, the relative contribution of these effects will differ. It is also interesting to consider how the suppression varies with the central system mass, M X . This is shown in Fig. 6 (right), and in both cases the dependence is seen to relatively mild. The suppression due to introducing the δ cut decreases at both low and high M X , due to the counteracting effects seen in Fig. 2: while increasing M X leads to a generally larger suppression due to the higher scale at which the PDF is evaluated, this also leads to a larger average x value probed, for which the suppression is less, with similar, but opposite,  effects for decreasing M X . Once soft survival effects are included, however, the overall trend is simply decreasing with M X .
A further observable of interest is the relative contribution to the photon-photon luminosity from the different components of the input photon PDFs, i.e. the coherent, incoherent and evolution defined above, and which roughly speaking correspond to exclusive, low and high-mass diffractive production at the corresponding proton vertex, respectively. This is shown in Fig. 7, again for the inclusive and semi-exclusive (with/without survival effects) cases. Considering the inclusive case in Fig. 7 (top left), we can see that the incoherent input component is found to be very small, consistently with the results of [26], while as M X increases, the relative contribution from the PDF evolution increases, due to the larger evolution length, although this effect is somewhat softened due to the fact that at higher M X the average x probed increases, where the impact of PDF evolution is less significant. Next, we consider the semi-exclusive case in Fig. 7 (top right). The most dramatic effect is 12.2 2.4 0.015 σ δ=3 /σ inc. 0.18 0.16 0.14 σ δ=7 /σ inc. 0.39 0.36 0.31 Table 2: Cross section predictions for photon-induced muon and W boson pair production at the √ s = 13 TeV LHC. The inclusive and ratio of semi-exclusive to inclusive cross sections are shown, with two choices of emission region δ. In the muon pair case, results for two values of cut on the transverse momenta are shown. For the W boson case, the leptonic decay (l = e, µ) is considered, with lepton transverse momenta p l ⊥ > 25 GeV, and missing transverse energy E miss ⊥ > 20 GeV. In all cases the lepton pseudorapidity is required to satisfy |η| µ < 2.5.
that the relative contribution from the evolution components is significantly reduced: this is a simple consequence of the fact that the constraint (12) only applies to the PDF evolution. Finally, considering the case with survival effects in Fig. 7 (bottom), we can see that the relative contributions change consistently with the results of Table 1.
The most important consequence of the above results is the suppression seen in the contribution from the evolution components, due both to the impact of (12) on the evolution and the smaller survival factors for this less peripheral interaction. While in the inclusive case the evolution component is larger than the purely coherent one, for the semi-exclusive the result is dramatically different, with the contribution due to the evolution components of both protons being completely negligible. Apart from at the highest M X , it is only the 'coherent, evolution' contribution, due to the larger survival factor and smaller impact of (12), which remains significant. Interpreting this result physically, for semi-exclusive production we expect some M X dependent contribution from exclusive and single proton dissociation, while the double proton dissociation contribution is only noticeable at very high M X .

Cross section predictions at the LHC
In this section we present a short selection of cross section predictions at the LHC. We will consider for concreteness the cases of muon and W boson pair production at √ s = 13 TeV, although the approach described above may be readily applied to other two-photon induced processes. In principle it would also be possible to apply this approach to the semi-exclusive photoproduction of e.g. vector meson (J/ψ, Υ...) states, however the situation is then greatly complicated by the presence of the QCD-induced vertex, due in the language of Regge theory to pomeron, rather than photon, exchange; for this reason we do not consider it here.
In Table 2 we show cross section predictions for photon-induced muon and W boson pair production at the √ s = 13 TeV LHC. We show the inclusive cross section and ratio of semiexclusive to inclusive cases, with survival effects included and for two values of δ (applied to  Table 3: Predicted relative contribution of exclusive, single (SD) and double (DD) dissociative components for the same processes as in Table 2. Results are rounded to two significant figures, so may not sum to unity precisely.
both proton sides). These choices are motivated by the experimental situation at the LHC. Considering just the central tracking detector at ATLAS/CMS, we have an uninstrumented region beyond |η| ∼ 2.5, which for √ s = 13 TeV corresponds to roughly δ ≈ 7. On the other hand, a much larger rapidity gap may be vetoed on using forward shower counters, which are currently installed at LHCb (the HERSCHEL forward detectors [17]), CMS [17] and ALICE [18]. Roughly speaking, these extend rapidity coverage out to |η| 8, corresponding to δ ≈ 2 is representative in this case. These can be considered as lower and upper bounds on experimentally realistic values of δ: for other specific experimental configurations, the appropriate value may lie somewhere in between. However, δ = 2 represents a relatively small region for non-vetoed emission, which may be sensitive in particular to fluctuations due to fragmentation and hadronisation, see e.g. [40], and for which the transverse momentum of the outgoing quark due to DGLAP evolution satisfying (12) may not be sufficiently high that the simple factorisation (26) holds; for example, if we take a characteristic z ∼ 0.2 x in (12) we have q t 1.5 GeV. As well as potentially spoiling this factorisation, this relatively low scale indicates that such a veto may be sensitive to the incoherent input component of the photon PDF, for which the recoil quarks may be produced with sufficiently large q t to fill the rapidity gap. We thus present results for a somewhat higher value of δ = 3, in order to avoid too great a sensitivity to these effects.
We can see from Table 2 that even for the relatively large δ = 7 there is expected to be a significant reduction, by a factor of ∼ 3, in the predicted cross section relative to the inclusive case; for δ = 3 this reduces by a further factor of ∼ 2, depending on the process. The greater suppression with increasing M X , already seen in Fig. 6, is also evident. It would be interesting to test these predicted trends, which we recall is a non-trivial result of the formalism described above, in particular that used to model survival effects.
In addition, we may investigate the relative contributions of the pure exclusive, single dissociative, and double dissociative components of the photon PDFs to these processes, described in the preceding section. This is shown in Table 3: we can see the clear relative enhancement of the pure exclusive component which becomes more pronounced as δ is decreased. While such a qualitative enhancement in the exclusive component is to be trivially expected when such a rapidity gap veto is applied, the quantitative predictions depend on this particular approach. It is also observed that an increase in the double dissociative component is predicted as δ increases. We end the section with some comments about the theoretical uncertainty in the above predictions. First, it should be emphasised that the coherent component of the input photon PDF, due to purely QED emission from the proton, is theoretically very well understood. Moreover, for purely exclusive processes, i.e. due to coherent emission from both protons, it was found in [8] that the survival factor is almost independent of the specific model taken, as the main model dependence in the evaluation of the soft survival factor lies in the region of small impact parameter b t R p (where R p is the proton radius), whereas for coherent photon emission the average impact parameter is much larger, and so it is relatively insensitive to this region.
For the incoherent input component, which depends on such issues as the treatment of the quark masses, the uncertainty is larger. We recall that freezing the input u and d quark distributions in (20) at the starting scale, as is done by default, represents an upper limit on the incoherent component of the input distribution. We can instead take the non-relativistic quark model expectations, see [26], which represent a lower limit on the incoherent component. In this case, for example, the incoherent-coherent component of the γγ luminosity in Fig. 7 (bottom) is a factor of ∼ 2 smaller at lower M X . However, this effect is typically washed out in observable predictions: the inclusive cross sections in Table 2 are a factor of ∼ 10% smaller, with the ratio of semi-exclusive to inclusive almost unchanged, while the inclusive and δ = 7 predictions in Table 3 are similarly effected at the 10% level in the case of muon pair production, with the exclusive component being somewhat larger and the single and double dissociation somewhat smaller, while for the higher mass W + W − cross section, the effect is even smaller. The most sensitive case are the δ = 3 predictions for the lower M X muon pair case, for which the exclusive component is ∼ 20% bigger, the single dissociation ∼ 10% bigger, and the double dissociation is a factor of ∼ 2 smaller; the effect is similar, but less significant for the W + W − cross section. Thus, generally the uncertainty on the incoherent contribution to the input photon PDF affects the cross section prediction at the level of 20% or lower, with the exception of the double dissociation cross section for lower δ, which displays a larger sensitivity. However, this uncertainty may be further reduced by refining the treatment of the input distribution to include higher mass (∆...) resonances, as well as potentially by constraining it with inclusive measurements that are sensitive to the photon PDF.
Finally, we must consider the evolution component, where as discussed in the previous section there is some uncertainty in the choice of form factor for the coupling of the photon GPDF to the proton, as in (26). If instead of assuming the same coupling of the GPDF to the GW eigenstates as to the pomeron, we take a universal coupling, this leads to a survival factor of 0.15 for the evolution-evolution component, i.e. roughly 50% larger than the default value in Table 1 10 . However, if we instead consider the coherent-evolution component, then the survival factor is insensitive, at the level of a few percent, to this choice: the interaction is automatically highly peripheral, due to the coherent photon-proton vertex, so that any contribution from the vertex for the evolution component is largely washed out. Thus, while it might appear at first sight that the uncertainty associated with the choice of coupling may be quite large, this is in fact not the case. From Fig. 7 (top right) we can see that the contribution from the evolution-evolution component is small, ∼ 10%, even prior to including survival effects, for reasonable values of M X . After including the survival factor, this contribution becomes much smaller still, so that the final cross section prediction will be relatively insensitive to its precise value. In fact, the evolution component enters dominantly through the mixed evolution-coherent term in the γγ luminosity, which as discussed above carries little uncertainty due to this choice.

Comparison to CMS semi-exclusive µ + µ − data
In [22] the CMS collaboration have presented a measurement of semi-exclusive W + W − production via leptonic decays at √ s = 7 TeV. That is, events are selected by demanding a dilepton vertex with no additional associated charged tracks within the tracker acceptance (η < 2.4). In this way the contribution from exclusive γγ → W + W − is enhanced, and the corresponding cross section measurement can be used to set limits on anomalous quartic gauge couplings. To compare the measured semi-exclusive data with the purely exclusive prediction, a correction factor derived from a larger data sample of µ + µ − events selected with the same track veto and in the same invariant mass region as the W + W − events was applied. Namely, the ratio F of the total number of observed µ + µ − events, after correcting for the remaining Drell-Yan background, to the prediction found with the LPAIR MC [41,42], corresponding to the equivalent photon approximation without any survival factor applied, is taken. Assuming the same factor applies in the W + W − case this allows the measured semi-exclusive cross section to be converted to an exclusive one and for corresponding limits on anomalous coupling to be set.
In this paper we are not interested in such limits, but rather on this ratio F : assuming that the CMS track requirement corresponds to good approximation to a veto on all additional particle in the |η| < 2.4 region, this is precisely the experimental situation considered in this paper. In particular, for √ s = 7 TeV this corresponds to a rapidity veto with δ ≈ 6.5, as defined in Section 2, for µ + µ − production. In Table 4 we show our prediction for this scenario 11 and with the same cuts applied on the final-state muons as in [22]. We can see that F Inclusive 10.9 δ = 6.5 3.6 δ = 6.5, γ incoh = 0 3.0 CMS [22] 3.23 ± 0.53 Table 4: Predictions for ratio F of semi-exclusive to exclusive µ + µ − cross sections, with M µµ > 160 GeV, p µ ⊥ > 20 GeV, |η µ | < 2.4 at √ s = 7 TeV, for a rapidity gap veto with δ = 6.5, compared to the measurement of [22] which closely corresponds to this scenario. Consistently with [22], no survival factor is applied for the purely exclusive cross section. the prediction is in excellent agreement with the CMS measurement of F = 3.23 ± 0.53. Such an encouraging level of agreement, given the various ingredients that enter in the calculation, is not trivial; we recall in addition that the uncertainty on such a result is expected to be at the ∼ 10% level, see the discussion in the preceding section, and hence is of the same order as the experimental uncertainty. We also show for comparison the prediction with the incoherent contribution (20) turned off; this result, which is about ∼ 20% lower, is also completely consistent with the data within current uncertainties. With further data and a reduction in the (dominantly statistical) uncertainty on the measured F it may become possible for the incoherent contribution in our approach to be pinned down, although given that the theoretical uncertainties are of the same order as this 20% difference, this may be challenging.
It is also interesting to consider the predicted ratios F using other available photon PDFs, namely the NNPDF2.3QED [2], CT14QED [3] and the older MRST2004QED [1] sets. This is achieved by applying the same procedure as above, but instead using these sets for the input PDF in (13) at the corresponding staring scale Q 0 of the set. Unfortunately, for the NNPDF set one added complication is that as no separation is made between the coherent and incoherent inputs, there is a certain amount of freedom in how to treat the survival factor in these cases, which from Table 1 we can see are generally quite different. As the coherent contribution is expected on general grounds to be dominant, we assume the inputs to be purely coherent when calculating the survival factor, however a more complete treatment would give a somewhat smaller value for F ; making the unphysical assumption that the input is completely incoherent, the predicted F is ∼ 30 − 40% smaller. More realistically, the correct predictions for F could be ∼ 10% smaller than the quoted value. For simplicity, we make the same assumption for the MRST2004QED set, with the results of a more precise treatment being well within the PDF uncertainty range. The CT14QED set is constrained using ZEUS data on isolated photon production in DIS [43], and is interpreted in [3] as being due to the purely inelastic production 12 . In the context of our approach, this corresponds CMS [22] 3.23 ± 0.53 Table 5: Predictions for ratio F as in Table 4, but with the input photon PDF calculated using the NNPDF2.3QED [2] (NLO in QCD, α s (M Z ) = 0.118), CT14QED [3] and the older MRST2004QED [1] sets. The MRST2004 range corresponds to the constituent and current quark mass results, the CT14 range to the results with the photon momentum fraction p γ 0 between 0 -0.14%, and the NNPDF2.3 uncertainties correspond to a 68% confidence envelope.
to a purely incoherent input component at Q 0 : we therefore add this to the coherent input (16), with the corresponding survival factors included as described in the previous sections.
In Table 5 we show predictions for these three PDF sets, with uncertainties calculated as described in the Table caption. In all cases these are completely consistent with the CMS data, within the very large PDF uncertainties. The NNPDF and MRST predictions, which include no explicit coherent input, span a large range above and below the measured value of F . On the other hand the CT prediction, for which we have included the coherent contribution explicitly, predicts F 3, consistently with our prediction with γ incoh = 0 in Table 4 (the values are not precisely the same as the evolution is now performed at LO in α s , consistently with [3]), but extends to significantly higher values of F than in our approach, well beyond the measurement. Given the size of the PDF uncertainties, it is clear that including constraints of this type, even allowing for some conservative theoretical uncertainty, could have a dramatic effect in constraining the photon PDF within these approaches. It is also worth emphasising that the data are in excellent agreement with a dominantly coherent contribution, as calculated in this paper.

Conclusions
There are two reasons to include the photon as a parton in the proton. First, in inclusive production the influence of electroweak corrections is increasingly relevant as we enter the era of precision LHC phenomenology, where NNLO QCD calculations are becoming the standard reconstructed. Such a requirement will remove the contribution from the coherent component, and from MC is found to correspond to a constraint of W X > 5 GeV on the mass of the produced hadronic system [44]. As at least part of the input incoherent component is expected to fail this extra track requirement, some care may be needed in using these data to constrain the photon PDF, although a study of this is beyond the scope of the current work. for many processes. In this case it is necessary to introduce the photon PDF in exact analogy to the conventional PDFs of the quarks and gluons, in order to include such corrections consistently. However, the importance of the photon PDF is also seen in more specific photon-induced processes, such as vector meson photoproduction and γγ-initiated lepton or W boson pair production. As a rule to select these relatively rare events experimentally some additional cuts, including generally a requirement of large rapidity gaps on either side of the produced object, must be imposed. These cuts affect the incoming photon luminosity and require us to modify the corresponding photon distribution, which no longer corresponds to the usual inclusive one.
In this paper we have considered photon-initiated processes in proton-proton collisions, where one or both protons break up following the interaction. Due to the colour-singlet nature of the photon exchange, these can lead naturally to rapidity gaps in the final-state between the centrally produced system and the proton decay products. The question, which we have attempted to address, is then what the cross section is for a two-photon initiated process to pass a specific experimental rapidity gap veto. We have shown here how such a veto may be accounted for at leading order by a relatively straightforward adjustment to the usual DGLAP evolution of the photon PDF to account for the kinematic condition implied by the rapidity gap requirement. In addition to this, it is necessary to include the probability for no extra particle production in the gaps due to underlying event activity: the so-called 'survival factor'. This object is dependent on soft physics and so introduces an additional model dependence to the calculation. However, we have shown explicitly that provided the rapidity region for allowed emission is not too restrictive, then the uncertainty due to this is quite limited. This is in fact not surprising, and is a result of a more general point: it is the relatively simple nature of the photon emission (in particular, the lack of colour exchange induced) which allows for the rapidity gap veto to be included in the DGLAP evolution in a straightforward way, while the peripheral nature of the photon exchange leads to the uncertainty due to the survival factor in the observable cross section generally being small.
As a result of the effects described above we expect the 'effective' γγ luminosity and hence predicted cross section to be reduced when an experimentally realistic rapidity gap veto is applied to both sides of the produced state, with the impact of vetoing on quark emission in the PDF DGLAP evolution and the survival factor being comparable in size. The level of suppression depends on the exact gap size, as well as the hard process and event selection criteria; some representative predictions for realistic experimental scenarios have been presented in Table 2, where the overall suppression is seen to lie in the range ∼ 2 − 7. In Table 3 the breakdown between the fraction of purely exclusive events and those with one or both protons dissociating for a selection of scenarios has been given; while demanding rapidity gaps will naturally lead to an enhancement of exclusive events, the precise predictions for these fractions and their dependence on the gap size is a non-trivial result of the theoretical approach, in particular soft survival effects, which do not enter uniformly and are expected to suppress the dissociative contributions beyond naive estimates. Depending on the gap size, we expect the fraction of exclusive events to be ∼ 20 − 50 %, with the remaining events being principally due to single proton dissociation, while the contribution from double dissociation is below ∼ 15%, and may be significantly lower. In all cases, we have demonstrated that it is crucial to fully account for the experimental cuts, and in particular the rapidity gap veto, when comparing data to theoretical predictions for photon-initiated production; it is not sufficient to use the conventional inclusive photon PDF.
Such rapidity gap events are of much phenomenological interest at the LHC, with a range of measurement possibilities, using just the central tracking detector at ATLAS/CMS, or across the wider rapidity region allowed by the forward shower counters currently installed at LHCb (the HERSCHEL forward detectors [17]), CMS [17] and ALICE [18]. We have therefore presented a selection of cross section predictions for semi-exclusive photon-initiated lepton and W boson pair production at the √ s = 13 TeV LHC, with rapidity gap veto regions relevant to both these scenarios. The cross sections are quite large, and the relatively simple formalism outlined in this paper leads to some clear predictions for such observables as the ratios of the single and double proton dissociative contributions to the total cross section, and their dependence on the central system mass M X and rapidity gap veto size. By measuring these at the LHC, a better understanding of the photon PDF and models of soft physics may be possible.