The Practical Pomeron for High Energy Proton Collimation

We present a model which describes proton scattering data from ISR to Tevatron energies, and which can be applied to collimation n high energy accelerators, such as the LHC and FCC. Collimators remove beam halo particles, so that they do not impinge on vulnerable regions of the machine, such as the superconducting magnets and the experimental areas. In simulating the effect of the collimator jaws it is crucial to model the scattering of protons at small momentum transfer~$t$,as these protons can subsequently survive several turns of the ring before being lost. At high energies these soft processes are well described by Pomeron exchange models. We study the behaviour of elastic and single-diffractive dissociation cross sections over a wide range of energy, and show that the model can be used as a global description of the wide variety of high energy elastic and diffractive data presently available. In particular it models low mass diffraction dissociation, where a rich resonance structure is present, and thus predicts the differential and integrated cross sections in the kinematical range appropriate to the LHC.We incorporate the physics of this model into the beam tracking code MERLIN and use it to simulate the resulting loss maps of the beam halo lost in the collimators in the LHC.


Introduction and motivation
The world's highest energy particle accelerator, the Large Hadron Collider (LHC), contains two high-energy proton beams travelling in opposite directions, guided around the accelerator ring by superconducting (SC) magnets. Its nominal stored beam energy of 360 MJ is orders of magnitude greater than previous accelerators, such as the Tevatron. This high energy stored beam passes in the machine aperture close to its magnet SC coils with a quench limit of about a e-mail: robert.appleby@manchester.ac.uk 40 mW/cm 3 at the operational current of 80 % of the magnet critical current [1]. The limit corresponds to the design limit of 13-15 mW/cm 3 including an assumed safety factor of 3 [2]. A powerful cleaning system is vital to the machine protection in order to operate below the quench limit, with a highly efficient collimation system necessary in order to remove any stray halo protons. The halo is generated by various effects [2] and it is characterised as an off-momentum halo (in which particle energies deviate from the reference) and a betatron halo (in which particles have large transverse amplitudes). Although the collimation system is adequate for the current configuration of the LHC, for the future High-Luminosity (Hi-Lumi) machine [3] upgrade the physics of the scattering of protons in the collimators must be accurately simulated, to avoid any quench of the SC magnets and to protect the vulnerable parts of the machine such as the detectors.
The tracking of protons around the ring and inside the collimator material is based on complex simulations where many different physics effects are involved. Here we focus on the scattering. Protons interact with both electrons and nuclei in the collimator material, with the former giving ionisation energy loss. The latter can be divided into elastic ( pp → pp), Single-Diffractive (SD) ( pp → pX or pp → X p), double diffractive ( pp → XY ) and inelastic scatters. Note that we ignore nuclear effects and consider a nucleus as a collection of protons and neutrons, and interactions with neutrons are treated similarly to those with protons. Experimentally pp and pn cross sections are within less than 2 % at the highestenergy pn data available, √ s = 30 GeV, and theoretically the agreement is expected to be even less than this, as the relevant processes are dominated by Pomeron and f 2 exchange. This approximation has also been used in previous studies of the LHC collimation system ( [4,5] and references therein), which agree with measured losses.
To study the beam halo we do not consider inelastic scatters, double-diffractive scatters, or SD interactions pp → X p, in which the beam proton breaks up: for such events all the energy is lost locally, within 50 m, according to studies by the CERN collimator group [6,7]. With elastic and singlediffractive scattering ( pp → pX) the emerging protons are only slightly affected and may survive several turns before being lost. The elastic scattering contributes to the betatron halo creation, and SD to the off-momentum halo.
The LHC ring is divided into 8 regions. For the nominal layout, as described in the design report [2], there are two collimation regions. In the third Interaction Region (IR3), the removal of off-momentum halo particles, known as momentum cleaning, takes place in a dispersive region. In IR7, particles with large transverse amplitude are removed; this is known as betatron cleaning. There is also an accelerating region in IR4, and a beam dump region in IR6. The remaining four regions are dedicated to the detector insertions: there are two at low β * in IR1 (ATLAS) and IR5 (CMS), and two at high β * in IR2 (ALICE) and IR8 (LHCb), where β * is the betatron function of the magnetic lattice at the interaction point. In each collimation region there is a cleaning hierarchy, and the primary collimators (TCP) in IR7 have the tightest apertures of the machine. In addition, tertiary collimators (TCT) are installed at both sides of the detector insertions to protect the final focus SC magnets and detectors.
In Table 1, the equivalent centre-of-momentum energy is given for various LHC proton energies on a 'fixed target' proton in the collimator. It varies from 29 GeV at injection (E beam = 450 GeV) to 115 GeV for the nominal beam energy, 176 GeV for the LHC energy upgrade and 306 GeV for the FCC-hh.
Experimental data for pp and pp reactions exist for many energies from different experiments and accelerators, principally the Intersecting Storage Rings (ISR) at √ s =23-63 GeV and the Tevatron at 2 TeV. There are also data from the SPPS. With plentiful data both above and below the range required, our model parameters are obtained by interpolation, rather than extrapolation.
In this paper we create a model within the Pomeron and Reggeon exchange framework of Donnachie and Landshoff [8,9]. The model is an elegant description of the strong interaction at high energies, and describes the experimental data for total, elastic and SD scattering with minimal assumptions. The fit uses a small number of parameters to describe data for 21 energies and 11 experiments, aiming to achieve the best possible fit. We use an extension of the model which we fit to most of the available elastic and SD data, in order to obtain a parametrisation which covers the required protontarget kinematical range at LHC energies. The extended model, which we simply call the DL model, is implemented into the beam tracking library MERLIN [10][11][12][13], which is then used to simulate the loss maps for the nominal LHC. We use this model to simulate the LHC loss maps, demonstrating the cleaning performance of the collimation system. This performance determines whether the accelerator can safely run at higher intensity, or whether additional shielding or collimators will be required. Realistic simulations of particle loss maps are fundamental to our ability to predict eventual quenching locations, for the nominal LHC and possible upgraded collimation systems, new materials and advanced collimation concepts such as hollow electron lenses [14] and crystal collimation [15].
The layout of this paper is as follows. In Sect. 2 we introduce the kinematics and discuss the requirements for the simulation of proton scattering within the collimator materials. In Sect. 3 we model the elastic scattering, performing a fit which achieves a good description of the available data. Then in Sect. 4 we describe the single diffractive model and obtain a fit for the double differential cross section for low and high missing mass regions, producing a good description of a wide range of data. We illustrate the fitting procedure and present the results at LHC energies and the prediction of the total SD cross section as a function of the centre-of-mass energy √ s. We show that it is possible to use the DL fit approach for elastic and SD scattering to cover the required range of kinematical variables for the LHC. In Sect. 5 we introduce the MERLIN code and the implementation of the model. The resulting loss maps for the nominal LHC at 7 TeV are presented with a detailed examination of the betatron cleaning region and the losses in the dispersion regions.
The data sources for elastic and single diffraction dissociation at different energies are reported in 1 and 1, along with references and the fit of the model to data.

Particle beam dynamics and dispersion
The horizontal transverse motion of a particle in an accelerator is given by the Courant-Snyder parameterization of the solution to Hill's equation [16] Fig. 1 Horizontal dispersion (green line) and horizontal and vertical β functions (red and blue lines) for beam 1 around the IR7 region in the LHC as generated by MERLIN using the V6.503 optics layout. The horizontal axis represents the distance from ATLAS interaction point Here s is the longitudinal position along the accelerator lattice, J is the particle action, β x the betatron function of the accelerator magnetic lattice, and μ x the betatron phase. p is the reference momentum for the lattice and p is the deviation of the particle from this reference momentum. D x is the dispersion function, describing the motion of particles with such a deviation. Protons that have lost momentum in diffractive interactions may have large transverse displacements from the reference orbit in regions where |D x | is large. Figure 1 shows the β-functions and the horizontal dispersion in the betatron collimation region IR7. There, the dispersion is small but the β-functions are large, so the collimators placed here remove protons in the betatron halo but not the energy halo. Figure 2 shows the β-functions and horizontal dispersion in the IR5 region, where the CMS detector is located. The magnetic elements are shown above the plot, including the quadrupole triplets on both sides of the detector that squeeze the beam at the interaction point. For the nominal LHC, the value of β at the interaction point is 55 cm in IR1 (ATLAS) and IR5 (CMS) and 10 m in IR2 (ALICE) and IR8 (LHCb).

Kinematics and the relevant range of t and M X
If a proton with mass M p and 4-momentum p μ i = (E i , p i ) in the lab frame interacts with a stationary proton P μ = (M p , 0), and scatters to p μ f = (E f , p f ), the invariants s and t are defined as s is the energy in the centre of momentum system. The expression for t can usefully be rewritten in terms of the proton scattering angle θ , The invariant mass of the diffracted proton-target, M X also called the missing mass, is given by energy-momentum conservation It is convenient to define the dimensionless variable The energy loss by the scattered proton is given by combining Eqs. 4 and 5 (For elastic scattering M X is equal to M p and this simplifies further.) Thus the quantities relevant for the simulation: E the energy change, and θ , the direction of the outgoing particle (apart from a random azimuthal angle) are determined by the quantities t and M X or ξ , and it is the distributions for these two quantities that are predicted by the model. If the scattering angle is significantly larger than the beam divergence, the scattered proton will be lost immediately, or in the nearby downstream region of the machine. Thus our model of elastic scattering must be accurate at small |t|/small θ but need not model large |t|/large θ , where 'large' and 'small' refer to comparison of the scattering angle with the angular beam divergence at that location. Table 2 shows the LHC V6.503 machine optics, characterised by the Twiss parameters α and β, at three typical collimators in IR7 (defined earlier) and, assuming a normalised nominal beam emittance of 3.75 mm.mrad, shows the angular The table shows that, for the collimation optics in IR7, modelling elastic scattering events up to |t| = 0.34 GeV 2 is sufficient to correctly model scattering events which could change the collimator-induced loss map beyond the immediately vicinity of the primary collimator. The table also includes an extreme case of collimation optics (arbitrarily chosen) which shows that for larger values of α at collimator locations we need to model an approximately double |t| range of elastic events. The elastic fits in this paper are valid over the range of available data, and extend to t = − 14.2 GeV 2 , which is more than ample for our simulations.
Detailed modelling of scattering at very small |t| is also unnecessary as very small angle scatters do not lead to beam loss. To investigate this we have used MERLIN to perform a full phase space aperture scan in both planes, injecting a beam filling one plane of phase space, i.e. a grid in x and x , with the other coordinates matched to the optical lattice at each collimator. These particles were then tracked for 100 turns, with particles removed if they touch any aperture restrictions, and the surviving particles' initial angles at the collimator jaws recorded. The smallest possible angle for which a particle is lost gives the minimum t value required. For the collimator jaw around the experimental regions, which have the minimum aperture available for scattering, an appropriate minimum value of |t| is 0.0001 GeV 2 .
The ranges of M X and |t| required for modelling diffractive scattering depend on the beam's angular divergence and its intrinsic energy spread. For the former, referring to Table 2, we take a conservative value of σ = 15 μrad to cover all possible current and future cases including deviation from the specified normalised emittance value. For the latter, the LHC beam energy spread σ e is the nominal 1.1 · 10 −4 at 7 TeV, and has been measured to be 1.36 ± 0.04 · 10 −4 at 3.5 TeV. The dependence of the scattering angle on ξ is weak, and a |t| limit of 4 GeV 2 corresponds to 20σ over all relevant ξ . For a 3.5 TeV beam energy, ξ = 0.12 corresponds to M X = 28 GeV and, even at our maximum |t| of 4 GeV 2 , this gives an energy deviation 420 GeV, which is 856σ e , and also 41σ . At 7 TeV it corresponds to 1109σ e , The conclusion for all energies is that a kinematical range of ξ up to 0.12 and |t| up to 4 GeV 2 is sufficient and conservative for the single diffractive fit. The fits are not sensitive to the minimum values, and we take the fits down to threshold for ξ and down to t = − 0.0001 GeV 2 .

Elastic proton scattering and the Pomeron
The differential cross-section dσ/dt of elastic pp and pp scattering is described by Coulomb scattering at very small |t| and nuclear scattering for larger |t|. Early measurements at the ISR [17,18] with energies between 23 GeV and 63 GeV revealed at low |t| an approximately exponential behaviour, e −B|t| , where B is known as the slope parameter. This is followed by a diffractive minimum at around |t| 1.4 GeV 2 , and subsequently a broad peak. The energy dependence of dσ/dt shows a shrinkage of the elastic peak, i.e. an increase in B, with increasing √ s [19]. The DL model includes Regge (ρ, ω and a 2 , f 2 trajectories) and Pomeron exchange [20], including multiple Regge and Pomeron exchanges [8,21]. At large t triple gluon exchange is also present [22,23]. Recently, in the light of the LHC data from the TOTEM experiment, a hard Pomeron term has also been added [24].
We extend the DL nuclear model to take into account the low t Coulomb peak in order to simulate elastic scattering in the energy ranges given in Table 1. The DL model has been fitted to all elastic data to obtain the fit parameters of the model. The fitting procedure is different from the one originally used in [24] where the normalisations were kept constant. The details of our approach are given here; a full account can be found in [25].
In this section we describe the elastic model general formulation and the fitting procedures. We then present the fitted differential cross section and the total elastic cross section.

The general formulation
The method used to calculate the elastic differential cross section is well established [26]. It is given by where f c and f n are the Coulomb and nuclear amplitudes, and the formalism includes the interference term. In general, there is a phase difference between the Coulomb and nuclear amplitudes, e iαφ(t) , such that To find the Coulomb phase φ, we use a fit to the cross section slope at t = 0, and the calculation by Cahn [27], The upper sign refers to pp scattering, the lower to pp, γ is Euler's constant and B is given by

Photon exchange
The Coulomb amplitude f c is given by [28] where F (t) is the proton electromagnetic form factor, given by equation 3.17 in [9], However it is well approximated by the simpler formula [9],

Hadronic exchange
In the DL model, 4 Regge trajectories are used, the hard Pomeron, the soft Pomeron, the f 2 and a 2 trajectory, and the ω and ρ trajectory. In the following these are labeled by 0 to 3 respectively. This gives, for the nuclear amplitude, where we have included a triple-gluon exchange amplitude. Its form varies like 1/t 4 at large t, is exponential at small t and the expressions are forced to match at an intermediate The amplitude is purely real and energy independent. t 0 is used as a free parameter in the fitting procedure. The exchange amplitudes A i are where s and u are the Mandelstam variables and F(t) is the form factor of the proton. The X i are real and positive; the factor i multiplying X 3 ensures the correct signature factor [8] for negative C-parity exchange. The amplitude for pp scattering is the same except that Y 3 has the opposite sign. The form of each Regge trajectory is where α i is the slope, and 1 + i is the intercept at t = 0. We extend this single scattering model to include double Pomeron exchange [24,25], which is necessary to account for the observed dip in the elastic differential cross section. The appropriate term 1 is included in our computed amplitude; it involves no new parameters apart from parameterising higher order scattering terms not included through scaling the double scattering amplitude by a factor λ.

The full model
Combining the contributions to Eq. (9) gives the differential cross section for elastic scattering as where the optical theorem has been used for the cross term and we define the ratio of the real and imaginary components of the nuclear term as ρ. The real and imaginary terms are the sum of all the corresponding terms from Regge exchange, and both single and double Pomeron exchange. We have assumed that αφ is small, so that the exponential term could be expanded. For pp scattering, the sign of the terms must be inverted for the C = −1 ω and ρ trajectory, and the triple gluon term.

The elastic model fit
Using this model, we fit all suitable available elastic data. Since the electromagnetic cross section diverges as t −→ 0, a minimum t value must be defined otherwise the integrated cross section will be infinite.
Using MINUIT within ROOT [29], a global fit is performed over the data shown in Table 11 of 1, where √ s > 23 GeV. At and below this value, the fit quality of the model starts to degrade.
In the fitting, both the Regge trajectories in the model are "effective" trajectories and are initialized with values taken from a Chew-Frautschi plot. The Pomeron trajectories, X i factors, and both λ and t 0 are taken as free parameters, making a total of 14 parameters. However the Regge trajectory slopes are fixed to control the stability of the fit. Full systematic uncertainties are taken into account and correlations between experimental data sets are included. The fit is performed over all data, and over the full t range available, yielding a χ 2 /N DF = 4.00. This overall figure covers a considerable variation: for many datasets the fit quality is acceptable (χ 2 /N DF ∼ 1) but there are some features of some datasets where the model and the data systematically disagree in terms of the statistical errors, which are, particularly at low t in the peak, sometimes very small.
The resulting fit parameters and uncertainties given by MINUIT are given in Table 3. We note the two Pomeron intercepts are both soft ( 0 and 1 ) -in essence we started off with a hard and a soft Pomeron and the fit to the data pushed them together, leaving a single soft Pomeron [30,31] with a complicated t dependence.

Total and differential elastic cross section
The total elastic cross section is usually quoted as the contribution from the nuclear term only, Our values, calculated by integration over the differential cross section according to this convention, for the total elastic pp cross sections at LHC energies are listed in Table 4. In Fig. 3 we show a subset of results for the total elastic cross section fit against experimental pp and pp data. The black dashed vertical lines at 450 and 7000 GeV show the energy range of interest for the LHC collimation system. Figure 4 shows the fit for pp elastic scattering at √ s = 7000 GeV using TOTEM data from the LHC, showing the fit performance at very high energy. Figure 5 shows the fit The remaining plots for pp and pp elastic scattering are presented in 1.
All figures show the combined systematic and statistical errors. Each plot shows the differential dσ/dt distribution, with each experimental data set in a different colour. The black line is the fitted function with the parameters given in Table 3. The normalisations given in 1 have been applied to the data.

Single diffraction dissociation
Single diffraction dissociation in pp interactions is the process pp → pX, (22) in which one proton breaks up into a system X while the other scatters elastically. Diffractive kinematics are described by s, t and M X . In fixed-target SD events at LHC energies M X can vary from M p + m π to more than 50 GeV. The simplest description of high energy process is given in the diagram of Fig. 7a in which a Reggeon or a Pomeron is exchanged between the elastically-scattered proton and the system X . In the limit s M 2 X |t| and M 2 X not too small the process may be described by the triple-Regge model [8,9,[32][33][34] as illustrated in Fig. 7c and discussed in Sect. 4.1. For small values of the missing mass M X , around a few GeV, the system X is dominated by baryon resonances and requires a different treatment. A simple model [35], based on duality arguments, allows us to extend the fit to low mass where existing data are scarce. This is discussed in Sect. 4.2.
The advent of the LHC has renewed interest in diffraction dissociation [36][37][38][39]. The associated models go beyond the simple triple-Regge model, principally by the inclusion of absorptive corrections, and they are successful in describing the total single diffraction cross section and, in some cases, the double differential cross section d 2 σ/dtdξ at small t. In one sense our approach is less ambitious in that we use the standard triple-Regge model without modification. In another sense, however, it is much more ambitious as we attempt, successfully, to describe all existing single diffractive dissociation data in pp interactions.

High mass: triple-Regge formalism
The triple-Regge description, shown in Fig. 7, describes the pp SD cross section in the region of high ξ as the sum of Fig. 7 The triple-Regge description of high-mass diffractive dissociation. a The squared amplitude summed over all possible system X in the large triple-Regge limit (s M 2 X |t|). b The discontinuity across M 2 X of the scattering amplitude. c The total cross section d 2 σ/dtdξ as the sum of triple-Reggeon contributions contributions from triple-Regge exchanges [8] (and applied to the LHC in [34]). In Fig. 7c each of the upper exchanges carry momentum transfer t while the lower one carries zero momentum transfer; the f i (t), i = 1, 2, 3 are the couplings of the exchanges to the relevant hadrons and G 12 3 (t) is the triple-Reggeon vertex. In addition to the Pomeron, Reggeised f 2 , a 2 , ω, ρ exchanges are allowed so in principle we require a whole series of terms, given by PP PP RR RR RP PR RP PR P R P R P P R R .
Here P refers to the Pomeron and R to any of f 2 , a 2 , ω and ρ trajectories. A term ( 1 2 3 ) contributes to d 2 σ/dtdξ as Here the α i (t), i = 1, 2, 3 are the relevant Pomeron and Reggeon trajectories, s 0 is a scale factor that we take to be 1 GeV 2 and φ(α(t)) is the Reggeon phase. In practice only the first four terms of the series are required [33], those being ( P P P ), ( P P R ), ( R R P ) and ( R R R ). Further it is sufficient to consider the f 2 , a 2 , ω and ρ trajectories to be degenerate and use a single generic Reggeon exchange. The differential cross section may be written as with where i and k denote P or R as appropriate. The Pomeron trajectory [33] is given by and is an "effective" reggeon trajectory, which is a reasonable average of the f 2 , a 2 , ω, ρ trajectories. For the fitting procedure various parametrisations have been tested for the function g iik (t), The most effective of these parametrisations, and the one we use here, is given by Eq. (32), where A i , B i and C k are free parameters given by the fit to the data. The chosen parametrisation works perfectly at low energy but gives too high a cross section at √ s = 546/640 GeV. For this reason the triple Pomeron coupling term g PPP (t) is parametrized in a different way based on the possibility that it vanishes at where t 0 = −0.05 GeV 2 is the optimum value. At small t, where the triple-Pomeron coupling term dominates, the vanishing term in the parametrisation improves the double differential cross section fit but it reduces the agreement of the differential cross section dσ/dt with high-energy data. At the same time, at high energy and high t, the predicted differential cross section is lower than the data and this simple parametrisation fails. To improve the fit for both differential and double differential cross section the triple-Pomeron coupling parametrisation is divided into three different regions of t. For −0.25 ≤ t < −0.0001 we use for −1.15 ≤ t < −0.25 we use and for −4.00 ≤ t < −1.155 we use We use a linear parameterisation at low |t| to avoid unphysical behaviour and a modified form at high |t| to increase the integrated cross section.
One additional term is required: a "Reggeized" pion exchange term which is important at low t [33], this term (39) is kept fixed during the fitting procedure. In Regge theory the pion exchange term is given by where g 2 ππ p /4π = 14.4 [41] is the on-mass-shell coupling, m π is the pion mass, α π (t) = 0.93(t − m 2 π ) is the pion trajectory and F 2 (t) is the proton form factor. In Eq. (39), σ π 0 p (sξ) [mb], denotes the pion-proton cross section, modeled by [41] σ π 0 p (sξ) = 13.63(ξ s) 0.0808 + 31.79(ξ s) −0.4525 , so the overall SD double differential cross for high missing mass can be written as: The fitting procedure described in Sect. 4.3.

Low mass: background and resonances
The single-diffractive dissociation at low mass is a delicate issue in diffractive dissociation studies. A lot of pp → pX data in the resonance region is available at very low energy, much of which is not relevant, but some is used as a guide. Useful information to model the resonance region comes from data at s = 565 GeV 2 for t = −0.05 GeV 2 [42], where the d 2 σ/dξ dt are averaged over 1.5 ≤ M 2 X ≤ 2.5GeV 2 . Both Pomeron and Reggeon exchange conserve helicity, so resonance excitation is primarily through incremental angular momentum with no change in the quark spin. On the basis of the Gribov-Morrison rule [32,43] we expect the resonance to have spin-parity (1/2) + , (3/2) − , (5/2) + , (7/2) − etc. Also the dominant exchanges (Pomeron, f 2 , ω) are isoscalar, so we expect that the leading resonances produced are P 11 (1440), D 13 (1520), F 15 (1680) and G 17 (2190). The background to these leading resonances comes from the low-mass continuation of the high-mass model of Pomeron and Reggeon exchange.
Hadron-hadron scattering at low energies can be described by the sum of a few amplitudes for direct s-channel production; as the energy increases these resonances increasingly overlap and more amplitudes need to be considered. At high energies it can be described by the sum of a few simple Reggeon exchange amplitudes in the t-channel; as the energy falls more of these are required. The principle of duality asserts that these are two descriptions of the same physics, valid at lower and higher energies.
This can be extended to the principle of two component duality [41] in which the s-channel amplitudes comprise a smooth background which is dual to Pomeron exchange, and a set of resonances which is dual to Reggeon (non-Pomeron) exchanges.
We determine the background term first. We assume it is quadratic, unlike a previous analysis [42] which used a general polynomial, as this can match the triple-Regge form. The contribution vanishes at threshold, ξ th = (M p +m π ) 2 /s, and can therefore be written a and b are then determined by requiring that the background matches smoothly onto the high mass region at some chosen value ξ c which represents the division between 'low' and 'high' mass.
Writing the triple-Regge function ∂ 2 σ H M ∂t∂ξ (ξ, t, s) as A(ξ, t, s) the boundary conditions can be written as and the resulting background coefficients are given by To complete the low mass model we now add the resonances contribution. Each baryon resonance is parametrised by a Breit-Wigner function, with a mass m l and a width γ l . So the total contributions from resonances to the SD cross section at low-mass is given by with l = γ l q q l where q and q l are respectively the 3-momenta at M X and m l in the resonance rest frame, assuming that π p is the dominant final state. They are given by The data from [42] for t = −0.05 GeV 2 at √ s = 23.7 GeV are fitted as a sum of the background and these four leading resonances.
The t dependence in the resonance region comes from Schamberger [42]. For 1.5 ≤ M 2 X ≤ 2.5 GeV 2 , dσ/dt ≈ exp((13.2 ± 0.3)t). As there is some slight t-dependence in the background the double-differential cross section is obtained by multiplying Eq. (46) by exp(13.5(t + 0.05)) In terms of the variable ξ , this is If the resonance contribution is simply added to the background there will be a small step between high and low mass regions at ξ = ξ c . This is remedied by subtracting a small matching term R m linear in ξ , which is zero at threshold and equal to the magnitude of the resonance term at the matching point, The total resonance contribution can then be written as and the complete single-diffractive double differential cross section at low mass is given by

Fitting procedures
The large amount of data on soft diffraction dissociation is given in 1. It covers the ranges 17.2 < √ s < 546 GeV and 0.015 < |t| < 4.15 GeV 2 , thus spanning the energies and the range of momentum transfer required. However there are clear inconsistencies of normalisation between different data sets and considerable variation in quality. In some data, for example from the ISR, the experimental resolution is insufficient to delineate clearly the resonance from the triple-Regge region, so fits to these data were restricted to ξ > 0.01. The twelve parameters of the parametrisation of Eqs. (31)- (34), g iik , are obtained from a global fit over all the available data using MINUIT within ROOT [29].
Full systematic errors and correlations between experimental data sets are taken into account, and we consider two ways of doing this. The data are quoted with statistical and systematic errors. The former are due to Poisson statistics on the number of particles counted. The latter are dominated by uncertainties in the acceptance, and are common to all measurements made by a given experiment. These errors are strictly multiplicative but they are small enough in practice to be taken as additive, greatly simplifying the analysis.
If the systematic errors are not considered, the quantity to be minimised is where i is the number of the experiment, j the measurement within that experiment's dataset, and f i j , d i j and σ i j are respectively the fitted function, the measured cross section, and the quoted statistical error. If the experiment also quotes a systematic uncertainty on the acceptance of α i , so that each measurement has an error S i j = α i d i j , it is not possible to merely replace σ 2 i j in the denominator by σ 2 i j + S 2 i j , because although this correctly expresses the variance of an individual point, it does not take into account the correlation between points. One way to allow for the correlations is to include them in the χ 2 using the formula where V is the covariance matrix for the measurements. This is included in the fitting procedure and gives reasonable results, although it is not clear what the fit is doing to the normalisations of the individual experiments.
An alternative approach is therefore taken, in which each experiment's results are adjusted by a factor F i . These factors are included as parameters in the fit, with values initially set to 1, and a contribution of (F i −1) 2 α 2 i is added to the χ 2 for each such term. The results of this procedure are very similar (and, for test cases, identical) and the fitted value of F yielded useful information as to what the fit was doing to the normalisation of the individual experiments. F i values initially set to 1 have been found to vary between 0.9 and 1.15 which is within the 20-25 % of quoted systematic errors between different experiments.
The minimisation is performed over 5562 data points yielding a χ 2 /N DF = 8.61. The fit parameters are listed in Table 5 and the coefficients of the resonance contribution c l in the low mass region arising from the fit (we do not float   Table 6). The normalisation for each experimental data set is presented in Table 7.
Typical fits are shown in Figs. 8 and 9 for the double differential cross section; other results at different momentum transfer t and energies are reported in 1.
In Fig. 10 we show the fit of d 2 σ L M /d M 2 X dt to the data [42] at low mass. The background is the green line, the resonance structure is blue and the total in red. The resonance structure of the low mass region is reflected in the data with a strong peak at P 11 (1440) followed by a decreasing contribution from the remaining resonances. A comparison for the full range of the missing mass M X is presented in Fig. 11; the data at high mass are from [44].

Single differential SD cross section
Data for the single differential SD cross section, dσ/dt after integration over ξ , are available for a large range of energies and are used to compare the fit results and as a guide in the fitting procedure of the double differential cross section data. Fig. 12 The single diffraction differential cross section at √ s = 30.5 GeV. In blue the resonance, in green the background and in black the high mass contribution to the differential cross section. The total DL model is shown in red. The black points represent data from Albrow  [45]. We show in blue the resonance contribution, in green the background and in black the high mass contribution to the differential cross section. The total DL model is shown in red.
The DL model at low energies does not match the experimental data perfectly, due to an underestimated contribution from the resonances. At low t there is a strong contribution to the differential cross section from resonances and background, for higher t the high mass term dominates. The effect of the correction at low t from the triple-Pomeron term of the high mass contribution is visible, for example, in Fig. 14 where the DL model matches the data very well. The modification at high t is a reasonable compromise at low and high energies. Fig. 14 The single diffraction differential cross section at √ s = 546 GeV. In blue the resonance, in green the background and in black the high mass contribution to the differential cross section. The total DL model is shown in red. The black points represent data from UA4 collaboration

Integrated SD cross section s dependence
The energy dependence of the total SD cross section is a controversial topic. For energy ( √ s) below 25 GeV the standard Regge theory reproduces the SD cross section well, however it rises faster than the experimental observations at higher energy. This behaviour was already expected theoretically due to problems related with the violation of the unitarity at high energy, i.e. σ S D > σ tot and the Froissart bound [46]. Some different theoretical approaches have been attempted to overcome this problem, including the renormalisation of standard Pomeron flux to agree with the data [47] or decoupling of the triple Pomeron vertex [48].
Single-diffraction data at √ s = 7 TeV have been obtained by ALICE [49], CMS [50] and TOTEM [51]. The ALICE data are consistent with the integrated single-diffraction cross section increasing with energy. In contrast the CMS and TOTEM integrated cross sections, for 8 < M X < 350 GeV and 6.5 < M X < 1100 GeV respectively, appear to show a decrease with increasing energy when compared to extrapolations of conventional models. For M X < 3.4 GeV TOTEM [52] give an integrated cross section of (2.2 ± 2.17) mb and an upper limit of 6.31 mb at 95 % confidence level which is not inconsistent with the UA4 value of 3.0 ± 0.8 mb at √ s = 546 GeV and still allows, in principle, some increase with energy.
A full discussion of the problem and possible modifications to the models is given in [53]. Including sophisticated modifications in our fitting of multiple data sets is impractical, so we adopted the simple alternative of modifying the chosen matching point M c = √ sξ c . We do not include data above √ s = 546/630 GeV and ensure agreement with the UA4 integrated cross sections [45] and dσ/dt. This is achieved with a parametrisation given by where α = 0.6 GeV and s 0 = 4000 GeV 2 .
A comparison between our model and the UA4 results is presented in Table 8, with the experimentally determined SD cross section σ expt S D = 2σ S D , where σ S D is the integrated cross section, to take into account both arms of the SD diagram as usually quoted by experimentalists. The agreement between the data and the integrated model is very good. The total integrated single diffractive cross section, integrated up to ξ < 0.05, is shown in Fig. 15. The red line is the DL model and the experimental points are indicated with their normalisations. The vertical dashed blue lines represents the energy range of the LHC.
The model works well in this region, and some distance beyond. At very low energy the Regge model does not accurately describe the data, which is dominated by s-channel resonances, but this is well below our region of applicability.

Application of the model at the LHC energies
For the LHC collimation studies the range of centre-of-mass energy √ s is between 81-115 GeV for collisions. We use our model to predict the total and double differential SD cross section at LHC energies. The single diffractive cross sections are reported in Table 9 for ξ < 0.05 and ξ < 0.12, the first is given as reference because it is the upper limit used in many publications on the subject and the second to present the values that we are using in our model to predict the LHC loss maps. Figure 16 shows the contribution of the background and resonances to the double differential cross section at s = 114 2 GeV 2 for different values of the transfer momentum t (from left to right t = −0.01, t = −0.4 and t = −2 GeV 2 ). The blue line represents the resonance contribution, the red line the background contribution and the green line the total low mass fit. The black line represents the triple-Regge fit at high mass. The contributions of resonances and background are stronger at low t and decrease at medium t; for high t the contribution from resonances disappears and the cross section is dominated by the background.
For the high mass region the contribution of the leading trajectories and the pion-exchange term are shown in Fig. 17 for different value of the momentum transfer t. The plot on the left shows the contributions for low t, the triple Pomeron term, PPP in red, and the PPR term in blue, dominate at low ξ , with some contribution from RRP in green and RRR in cyan. At higher ξ the main contribution to the sum in green the total SD double differential cross section at low mass. The black line represents the triple-Regge fit at high mass, note that it goes up to ξ = 0.12

Fig. 17
Individual contributions of triple-Reggeon exchanges and pion exchange for the SD double differential cross section for the LHC at top energy, s = 114.6 2 GeV 2 , from left to right t = −0.01, t = −0.4 and t = −2 GeV 2 of all terms, in brown, is given by the RRP and the pionexchange term in black. At medium and high t the SD double differential cross section is dominated by the triple-Pomeron PPP term with some slight contribution from the remaining terms.

Simulation of LHC loss maps using MERLIN
We have incorporated the DL model, with the fitted parameters described in the previous sections, into the simulation code MERLIN, a C++ accelerator physics library [13] which has been extended to include proton collimation [10][11][12] with the aim of providing an accurate simulation of the Large Hadron Collider (LHC) collimation system. Up to now, the FORTRAN program SixTrack has been the main tool used to calculate dynamical aperture in the LHC , and, incorporating the scattering model known as K2 [54], has been used for many studies of the LHC collimation system and resulting the loss maps [4][5][6][7]55]. However a second program was felt to be desirable, as an independent check and to provide a more flexible and future-proof design to which new features, such as this model, can readily be added. The modular nature of MERLIN allows one to easily switch between the K2 scattering model and the DL model. As well as these nuclear processes, MERLIN also includes Coulomb scattering and other EM effects of protons in matter. The relative importance of these is discussed in [25] and will be presented in a future paper on the MERLIN code.
MERLIN was used to simulate the nominal optics at 7 TeV of the LHC in order to generate loss maps, using a thick lens tracking model. These can be generated for different optics configurations, e.g. the β-function at the interaction points (β*) and beam crossing angles. As a first step, we benchmark its prediction of the optical functions against the MAD-X program [56], with the excellent agreement shown in Fig. 18.
We then, using the same scattering physics models, benchmarked it against the standard SixTrack+K2 predictions for the loss map calculation [5,12,57]. These studies have demonstrated the quantitative agreement between the established tools for loss map analysis and MERLIN.
The loss maps are characterised by the local inefficiency which is plotted as a function of s, the distance around the ring. Here z is the longitudinal resolution (10 cm), N AB S is the number of particles absorbed in z and N tot coll is the total loss in the collimators along the whole machine. For the collimator z is set to the collimator length and N AB S are the total losses in the collimator.
In these simulations, if a proton undergoes an inelastic interaction inside the collimators, or if it touches the beam Fig. 18 The difference of the optical β-functions for the nominal LHC between MADX [56] and MERLIN pipe, it is considered lost: the particle is removed from the bunch and the location at which this takes place is recorded, using a default bin size of 10 cm. This is necessary as these conventions are those used by existing studies [4][5][6][7]55]; against which we benchmark the new code. They mark a sensible line between the study of long-distance (multi-turn) effects, which are the subject of simulations such as these, and short range effects which are the subject of programs like FLUKA [58] and GEANT4 [59]. When short-range effects are also needed, for example in the simulation of beam loss monitor responses, the lost particles are saved to a file and the shower evolution studied by one of these more detailed codes.
The main parameters for the loss map calculation are summarised in Table 10 for a squeezed beam with an IP beam separation and crossing angle applied. The beam impinges on the horizontal primary collimators (TCP.C6L7) in IR7, and is tracked for 200 turns. The transverse offset between the jaw surface and the impact point, called the impact parameter, is set to 1 µm.
The loss map calculated for this machine configuration and the DL scattering models is shown in Fig. 19. The plot is colour coded: black spikes represent losses in the collimator jaws, red spikes losses in warm elements of the accelerator, and most importantly blue spikes which indicate losses in the superconducting magnets. Using 64M simulated protons, MERLIN calculates a total loss inefficiency of 77.65 %, with 0.010 % lost in cold regions and 0.011 % in warm elements. The remaining protons are lost in the collimators.
The highest black peak in the map corresponds to the horizontal primary collimator in IR7, which has an aperture of 6 σ . This is the tightest aperture in the machine. The loss map has been compared to those computed with Sixtrack+K2 [3], and has been found to be consistent [11,57] with most of the losses in the dedicated collimation regions IR7 and IR3, and black peaks in the tertiary collimators used to protect the high luminosity insertion in IR1 (ATLAS ) and IR5 (CMS). The peaks downstream of IR7 and in the arcs IR7-IR1 and IR1-IR2 are particularly important because they show where halo protons touch the SC dipoles. These predictions allow us to understand where possible quenching events may occur and how to optimise the collimation system. The losses in the dispersion suppressor region of IR7 are shown in Fig. 20 along with the horizontal dispersion. The highest beam losses per unit length are in the primary collimators, followed by lower losses in the secondary and tertiary collimators. Between 19700 m and 20200 m there are not only losses on the collimators, appearing as sharp spikes with an efficiency 10 −4 m −1 or above, but also regions with efficiencies around 10 −5 to 10 −6 m −1 in the warm magnets between the collimators; these are mainly single diffracted Fig. 19 Loss map result for beam 1 with DL scattering physics. The beam is impinged on the primary collimator TCP.C6L7 in IR7 region, corresponding to the highest peak in the plot. The black peaks correspond to the losses in the collimator elements, the red ones to losses in the warm elements and the blue ones are the losses in the cold magnetic elements of the machine Fig. 20 Loss map in the betatron collimation region IR7. The local maxima of the dispersion (green line) are linked to the losses between 20300 m and 29400 m, which is in a cold region. The black peaks correspond to the losses in the collimator elements, the red ones in the warm elements and the blue one are the losses in the cold magnetic element of the machine protons with high momentum losses and scattered at high angles. MERLIN predicts more losses in this region than the K2 scattering routine [12] which has a different description of diffractive scattering. Losses downstream of IR7 in the dispersion suppressor, which are particularly sensitive areas, are shown in blue. Protons which experience single diffractive scattering in the bulk material of the collimator emerge with a transverse kick and a lower energy, so protons entering the dispersion suppressors, where the dispersion rises rapidly, can be lost in these cold areas. Most of the peaks in the cold part of the arc between IR7/IR8 are located at the local maxima of the dispersion, as shown for the first peak downstream of the dispersion suppressor in Fig. 20.
In Fig. 21 we present the distribution of energy lost |δp/ p| by protons impacting on the warm elements among the collimators in IR7. The plot shows a peak at low lost energy followed by a plateau till 10 % and a rapid drop off. Figure 22 shows the distribution of energy lost in the dispersion suppressors (DS1, DS2) and particles lost downstream of IR7. The range of energy lost is between 2-10 % for the first dispersion suppressor and 1-2 % for the second.

Conclusion
We have presented a development of the model of Donnachie and Landshoff for elastic and single-diffractive proton scattering for use in simulating collimation systems in high energy proton accelerators. The model includes a description of elastic scattering combining Coulomb with Regge exchange amplitudes, and a description of diffractive scattering that combines s-channel resonance formation with tchannel Regge exchange. It is valid over a wide range in the centre of mass energy √ s, the invariant 4-momentum We have taken elastic and diffractive scattering data from a large number of previous experiments, with different systematic errors, and fitted them with a small number of model parameters. The results have been incorporated into the MERLIN tracking code and this has been used to predict loss maps for 7 TeV running at the LHC, as a first contribution to a future study of the collimation system in the LHC and proposed High-Luminosity upgrades, where it can contribute significantly to the design of the necessary improvements to the collimation system. The model includes important physics in an effective way, which has been included in a practical simulation program with the potential of understanding and improving the performance of collimation systems at present and future high energy proton accelerators. and model fit The single diffraction dissociation data sources with their energies and references are reported in Table 12. The complete set of model fit comprise more than 300 plots. In Figs. 25 and 26 we present a selection of these plots to cover the full range of energies ( √ s) and ξ and for low and high momentum transfer t.
Elastic and inelastic data has recently been published at LHC energies [88][89][90][91] . This would provide an interesting further study, but is not relevant to the relevant kinematic ranges of t and ξ , as considered in this paper.  23 The pp elastic scattering model fit (fitted over all data) shown for a range of √ s, over the full t range, including the coulomb peak and down to the lower cut on t Fig. 24 The pp elastic scattering model fit (fitted over all data) shown for a range of √ s, over the full t range, including the coulomb peak and down to the lower cut on t