Measurement of Feynman-x Spectra of Photons and Neutrons in the

DESY 14-035 March 2014

ISSN 0418-9833

Measurement of Feynman-x Spectra of Photons and Neutrons in the Very Forward Direction in Deep-Inelastic Scattering at HERA

H1 Collaboration

Abstract Measurements of normalised cross sections for the production of photons and neutrons at very small angles with respect to the proton beam direction in deep-inelastic ep scattering at HERA are presented as a function of the Feynman variable xF and of the centre-of-mass energy of the virtual photon-proton system W . The data are taken with the H1 detector in the years 2006 and 2007 and correspond to an integrated luminosity of 131 pb−1 . The measurement is restricted to photons and neutrons in the pseudorapidity range η > 7.9 and covers the range of negative four momentum transfer squared at the positron vertex 6 < Q2 < 100 GeV2 , of inelasticity 0.05 < y < 0.6 and of 70 < W < 245 GeV. To test the Feynman scaling hypothesis the W dependence of the xF dependent cross sections is investigated. Predictions of deep-inelastic scattering models and of models for hadronic interactions of high energy cosmic rays are compared to the measured cross sections.

submitted to Eur.Phys.J.C

V. Andreev22 , A. Baghdasaryan34 , K. Begzsuren31 , A. Belousov22 , P. Belov10 , V. Boudry25 , G. Brandt46 , M. Brinkmann10 , V. Brisson24 , D. Britzger10 , A. Buniatyan13 , A. Bylinkin21,43 , L. Bystritskaya21 , A.J. Campbell10 , K.B. Cantun Avila20 , F. Ceccopieri3 , K. Cerny28 , V. Chekelian23 , J.G. Contreras20 , J.B. Dainton17 , K. Daum33,38 , C. Diaconu19 , M. Dobre4 , V. Dodonov10 , A. Dossanov11,23 , G. Eckerlin10 , S. Egli32 , E. Elsen10 , L. Favart3 , A. Fedotov21 , J. Feltesse9 , J. Ferencei15 , M. Fleischer10 , A. Fomenko22 , E. Gabathuler17 , J. Gayler10 , S. Ghazaryan10 , A. Glazov10 , L. Goerlich6 , N. Gogitidze22 , M. Gouzevitch10,39 , C. Grab36 , A. Grebenyuk3 , T. Greenshaw17 , G. Grindhammer23 , D. Haidt10 , R.C.W. Henderson16 , M. Herbst14 , J. Hladk`y27 , D. Hoffmann19 , R. Horisberger32 , T. Hreus3 , F. Huber13 , M. Jacquet24 , X. Janssen3 , H. Jung10,3 , M. Kapichine8 , C. Kiesling23 , M. Klein17 , C. Kleinwort10 , R. Kogler11 , P. Kostka35 , J. Kretzschmar17 , K. Krüger10 , M.P.J. Landon18 , W. Lange35 , P. Laycock17 , A. Lebedev22 , S. Levonian10 , K. Lipka10,42 , B. List10 , J. List10 , B. Lobodzinski10 , L. Lytkin8 , E. Malinovski22 , H.-U. Martyn1 , S.J. Maxfield17 , A. Mehta17 , A.B. Meyer10 , H. Meyer33 , J. Meyer10 , S. Mikocki6 , A. Morozov8 , K. Müller37 , Th. Naumann35 , P.R. Newman2 , C. Niebuhr10 , G. Nowak6 , K. Nowak11 , J.E. Olsson10 , D. Ozerov10 , P. Pahl10 , C. Pascaud24 , G.D. Patel17 , E. Perez9,40 , A. Petrukhin10 , I. Picuric26 , H. Pirumov10 , D. Pitzl10 , R. Plaˇcakyt˙e10,42 , B. Pokorny28 , R. Polifka28,44 , B. Povh12 , V. Radescu10,42 , N. Raicevic26 , T. Ravdandorj31 , P. Reimer27 , E. Rizvi18 , P. Robmann37 , ˇ alek28 , D.P.C. Sankey5 , R. Roosen3 , A. Rostovtsev21 , M. Rotaru4 , S. Rusakov22 , D. S´ M. Sauter13 , E. Sauvan19,45 , S. Schmitt10 , L. Schoeffel9 , A. Schöning13 , H.-C. Schultz-Coulon14 , F. Sefkow10 , S. Shushkevich10 , Y. Soloviev10,22 , P. Sopicki6 , D. South10 , V. Spaskov8 , A. Specka25 , M. Steder10 , B. Stella29 , U. Straumann37 , T. Sykora3,28 , P.D. Thompson2 , D. Traynor18 , P. Truöl37 , I. Tsakov30 , B. Tseepeldorj31,41 , J. Turnau6 , A. Valkárová28 , C. Vallée19 , P. Van Mechelen3 , Y. Vazdik22 , D. Wegener7 , E. Wünsch10 , ˇ acˇ ek28 , Z. Zhang24 , R. Zlebˇ ˇ c´ık28 , H. Zohrabyan34 , and F. Zomer24 J. Z´ 1

I. Physikalisches Institut der RWTH, Aachen, Germany School of Physics and Astronomy, University of Birmingham, Birmingham, UKb 3 Inter-University Institute for High Energies ULB-VUB, Brussels and Universiteit Antwerpen, Antwerpen, Belgiumc 4 National Institute for Physics and Nuclear Engineering (NIPNE) , Bucharest, Romaniaj 5 STFC, Rutherford Appleton Laboratory, Didcot, Oxfordshire, UKb 6 Institute for Nuclear Physics, Cracow, Polandd 7 Institut für Physik, TU Dortmund, Dortmund, Germanya 8 Joint Institute for Nuclear Research, Dubna, Russia 9 CEA, DSM/Irfu, CE-Saclay, Gif-sur-Yvette, France 10 DESY, Hamburg, Germany 11 Institut für Experimentalphysik, Universität Hamburg, Hamburg, Germanya 12 Max-Planck-Institut für Kernphysik, Heidelberg, Germany 13 Physikalisches Institut, Universität Heidelberg, Heidelberg, Germanya 14 Kirchhoff-Institut für Physik, Universität Heidelberg, Heidelberg, Germanya 15 Institute of Experimental Physics, Slovak Academy of Sciences, Koˇsice, Slovak Republice 16 Department of Physics, University of Lancaster, Lancaster, UKb 17 Department of Physics, University of Liverpool, Liverpool, UKb 2

1

School of Physics and Astronomy, Queen Mary, University of London, London, UKb 19 CPPM, Aix-Marseille Univ, CNRS/IN2P3, 13288 Marseille, France 20 Departamento de Fisica Aplicada, CINVESTAV, Mérida, Yucatán, Méxicoh 21 Institute for Theoretical and Experimental Physics, Moscow, Russiai 22 Lebedev Physical Institute, Moscow, Russia 23 Max-Planck-Institut für Physik, München, Germany 24 LAL, Université Paris-Sud, CNRS/IN2P3, Orsay, France 25 LLR, Ecole Polytechnique, CNRS/IN2P3, Palaiseau, France 26 Faculty of Science, University of Montenegro, Podgorica, Montenegrok 27 Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech Republicf 28 Faculty of Mathematics and Physics, Charles University, Praha, Czech Republicf 29 Dipartimento di Fisica Università di Roma Tre and INFN Roma 3, Roma, Italy 30 Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria 31 Institute of Physics and Technology of the Mongolian Academy of Sciences, Ulaanbaatar, Mongolia 32 Paul Scherrer Institut, Villigen, Switzerland 33 Fachbereich C, Universität Wuppertal, Wuppertal, Germany 34 Yerevan Physics Institute, Yerevan, Armenia 35 DESY, Zeuthen, Germany 36 Institut für Teilchenphysik, ETH, Zürich, Switzerlandg 37 Physik-Institut der Universität Zürich, Zürich, Switzerlandg

18

38

Also at Rechenzentrum, Universität Wuppertal, Wuppertal, Germany Also at IPNL, Université Claude Bernard Lyon 1, CNRS/IN2P3, Villeurbanne, France 40 Also at CERN, Geneva, Switzerland 41 Also at Ulaanbaatar University, Ulaanbaatar, Mongolia 42 Supported by the Initiative and Networking Fund of the Helmholtz Association (HGF) under the contract VH-NG-401 and S0-072 43 Also at Moscow Institute of Physics and Technology, Moscow, Russia 44 Also at Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 45 Also at LAPP, Université de Savoie, CNRS/IN2P3, Annecy-le-Vieux, France 46 Department of Physics, Oxford University, Oxford, UKb 39

a

Supported by the Bundesministerium für Bildung und Forschung, FRG, under contract numbers 05H09GUF, 05H09VHC, 05H09VHF, 05H16PEA b Supported by the UK Science and Technology Facilities Council, and formerly by the UK Particle Physics and Astronomy Research Council c Supported by FNRS-FWO-Vlaanderen, IISN-IIKW and IWT and by Interuniversity Attraction Poles Programme, Belgian Science Policy d Partially Supported by Polish Ministry of Science and Higher Education, grant DPN/N168/DESY/2009 e Supported by VEGA SR grant no. 2/7062/ 27 f Supported by the Ministry of Education of the Czech Republic under the projects LC527, INGO-LA09042 and MSM0021620859 g Supported by the Swiss National Science Foundation 2

h

Supported by CONACYT, México, grant 48778-F Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2 and Rosatom j Supported by the Romanian National Authority for Scientific Research under the contract PN 09370101 k Partially Supported by Ministry of Science of Montenegro, no. 05-1/3-3352 i

3

1 Introduction Measurements of particle production at very small polar angles with respect to the proton beam direction (forward direction) in positron-proton collisions are important inputs for the theoretical understanding of proton fragmentation mechanisms. Forward particle measurements are also valuable for high energy cosmic ray experiments, as they provide important new constraints for high energy air shower models [1, 2]. The H1 and ZEUS experiments at the ep collider HERA have studied the production of forward baryons (protons and neutrons) and photons, which carry a large fraction of the longitudinal momentum of the incoming proton [3–8]. These analyses have demonstrated that models of deep-inelastic scattering (DIS) are able to reproduce the forward baryon measurements if contributions from different production mechanisms are considered, such as string fragmentation, pion exchange, diffractive dissociation and elastic scattering of the proton [6, 7]. The forward photon production rate, however, is overestimated by the models by 50 to 70% [8]. The measurements also confirm the hypothesis of limiting fragmentation [9, 10], according to which, in the high-energy limit, the cross section for the inclusive production of particles in the target fragmentation region is independent of the incident projectile energy. Measurements in the DIS regime provide a possibility to investigate the process at different centre-of-mass (CM) energies of the virtual photon-proton system, W , within the same experiment. The studies of the energy dependence of particle production allow a test of the Feynman scaling [11] hypothesis, according to which particle production is expected to show a scaling behaviour, i.e. independence of the CM energy in terms of the Feynman-x variable, xF = 2p∗|| /W . Here p∗|| is the longitudinal momentum of the particle in the virtual photon-proton CM frame with respect to the direction of the beam proton. In several previous measurements Feynman scaling was found to be violated in the fragmentation process in the central rapidity region [12–21]. On the contrary, no sizable violation of Feynman scaling has been observed in the target fragmentation region in pp and p¯ p collisions by comparing the π 0 production cross sections at the SPS collider [22] with π ± measurements at the ISR [23–26]. However, these conclusions are debated [27] and the scarcity of other experimental forward particle production data motivates further studies of forward particle production. In this paper the production of forward neutrons and photons in DIS is studied as a function of xF 1 and W . This is the first direct experimental test of Feynman scaling for photons and neutrons produced in the very forward direction. Predictions from different DIS and different cosmic ray (CR) hadronic interaction Monte Carlo (MC) models are compared to the results. The simultaneous measurement of forward neutrons and photons provides a useful input for MC model development also because of the respective different production mechanisms: forward photons almost exclusively originate from decays of neutral mesons produced in the fragmentation of the proton remnant (Figure 1a), while forward neutrons are produced also via a colour singlet exchange process (Figure 1b). 1

In the kinematic range of this measurement the variable xF is numerically almost equal to the longitudinal momentum fraction xL used in the previous publications [3–8]. There, xL was defined as xL = En,γ /Ep , where Ep , En and Eγ are the energies of the proton beam, the forward neutron and the forward photon in the laboratory frame, respectively.

4

e’

(a)

e’

(b)

e (k)

e

γ* γ * (q) W

X

π+

X

p (p)

p

t

n, γ

n

Figure 1: (a) Generic diagram for forward photon or neutron production ep → e′ γX, ep → e′ nX in deep-inelastic scattering. (b) Diagram of forward neutron production via pion exchange.

The neutrons and photons studied here are produced at polar angles below 0.75 mrad and are measured in the Forward Neutron Calorimeter (FNC) of the H1 detector. The data used in this analysis were collected with the H1 detector at HERA in the years 2006 and 2007 and correspond to an integrated luminosity of 131 pb−1 . During this period HERA collided positrons and protons with energies of E √e = 27.6 GeV and Ep = 920 GeV, respectively, corresponding to a centre-of-mass energy of s = 319 GeV.

2 Experimental Procedure and Data Analysis 2.1 H1 main detector A detailed description of the H1 detector can be found elsewhere [28–33]. Only the detector components relevant to this analysis are briefly described here. The origin of the right-handed H1 coordinate system is the nominal e+ p interaction point. The direction of the proton beam defines the positive z–axis; the polar angle θ is measured with respect to this axis. Transverse momenta are measured in the x–y plane. The pseudorapidity is defined by η = − ln [tan(θ/2)] and is measured in the laboratory frame. The interaction region is surrounded by a two-layer silicon strip detector and large concentric drift chambers, operated inside a 1.16 T solenoidal magnetic field. Charged particle momenta are measured in the angular range 15◦ < θ < 165◦ . The forward tracking detector is used to supplement track reconstruction in the region 7◦ < θ < 30◦ and improves the hadronic final state reconstruction of forward going low momentum particles. The tracking system is surrounded by a finely segmented liquid argon (LAr) calorimeter, which covers the polar angle 5

range 4◦ < θ < 154◦ with full azimuthal acceptance. The LAr calorimeter consists of an electromagnetic section with lead absorber and a hadronic section with steel absorber. The total depth of the LAr calorimeter ranges from 4.5 to 8 hadronic interaction lengths λ. The absolute electromagnetic energy scale is known with a precision of 1%, while the absolute hadronic energy scale is known for the present data with a precision of 4%. The backward region (153◦ < θ < 177.8◦ ) is covered by a lead/scintillating-fibre calorimeter called the SpaCal; its main purpose is the detection of scattered positrons. The polar angle of the positron is measured with a precision of 1 mrad. The energy resolution for positrons is p σ(E)/E ≈ 7.1%/ E[GeV] ⊕ 1% [34] and the energy scale uncertainty is less than 1%. The hadronic energy scale in the SpaCal is known with a precision of 7%. The luminosity is determined from the rate of the elastic QED Compton process with the electron and the photon detected in the SpaCal calorimeter, and the rate of DIS events measured in the SpaCal calorimeter [35].

2.2 Forward detector for neutral particles Neutral particles produced at very small polar angles with respect to the proton beam direction can be detected in the FNC, which is situated at a polar angle of 0◦ at 106 m from the interaction point. A detailed description of the detector is given in [7, 8]. The FNC is a lead–scintillator sandwich calorimeter. It consists of two longitudinal sections: the Preshower Calorimeter with a length corresponding to about 60 radiation lengths or 1.6λ and the Main Calorimeter with a total length of 8.9λ. The acceptance of the FNC is defined by the aperture of the HERA beam-line magnets and is limited to scattering angles of θ < ∼ 0.8 mrad with approximately 30% azimuthal coverage. The longitudinal segmentation of the FNC allows an efficient discrimination of photons from neutrons. Photons are absorbed completely in the Preshower Calorimeter, while neutrons have a significant fraction of their energy deposited in the Main Calorimeter. Therefore, energy deposits in the FNC, which are contained in the Preshower Calorimeter with no energy deposits above the noise level in the Main Calorimeter, are classified as electromagnetic clusters. According to the Monte Carlo simulation about 98% of all reconstructed photon and neutron candidates originate from generated photons and neutrons, respectively. Due to the relatively large size of the FNC readout modules in combination with the small geometrical acceptance window, two or more particles entering the FNC are reconstructed as a single cluster. In the MC simulation about 1% of all hadronic clusters in the FNC associated with neutrons are overlapped with a photon, which was scattered within the FNC acceptance together with the neutron. At lower energies the electromagnetic clusters reconstructed in the FNC mainly originate from single photons. At higher measured energies there is a significant contribution from two photons, with the fraction of two-photon events increasing from 0.5% at 100 GeV to 10% at about 450 GeV and to 80% at 900 GeV. The two photons typically originate from the decay of the same meson. The absolute electromagnetic and hadronic energy scales of the FNC are known to 5% [8] and 2% [7] precision, respectively. Thepenergy resolution of the FNC calorimeter for electromagnetic showers is σ(E)/E ≈ 20%/ E [GeV] ⊕ 2% and for hadronic showers σ(E)/E ≈ 6

p 63%/ E [GeV] ⊕p 3%, as determined in test beam measurements [36]. The spatial resolution is σ(x, y) ≈ 10cm/ E [GeV] ⊕ 0.6 cm for hadronic showers starting in the Main Calorimeter. A better spatial resolution of about 2 mm is achieved for the electromagnetic showers and for those hadronic showers which start in the Preshower Calorimeter.

2.3 Cross Section Definition The kinematics of semi-inclusive forward photon and neutron production are shown in Figure 1a, where the four-vectors of the incoming and outgoing particles and of the exchanged virtual photon γ ∗ are indicated. This measurement is restricted to the DIS kinematic range, determined by the photon virtuality 6 < Q2 < 100 GeV2 and inelasticity 0.05 < y < 0.6. They are defined as p·q Q2 = −q 2 , y= , (1) p·k

where p, k and q are the four-momenta of the incident proton, the incident positron and the virtual photon, respectively. The CM energy of the virtual photon-proton system, W , is related p to Q2 and y as W ≈ ys − Q2 , where s is the squared total CM energy of the positron-proton system. The present analysis is restricted to the range 70 < W < 245 GeV.

The analysis is performed in the pseudorapidity range η > 7.9 for forward neutrons and photons. The pseudorapidity range η > 7.9 corresponds to polar angles θ < 0.75 mrad. In the virtual photon-proton CM frame the neutron transverse momentum p∗T and the neutron xF are restricted to the range 0 < p∗T < 0.6 GeV and 0.1 < xF < 0.94, respectively. For the forward photons measurement p∗T and xF are defined for the most energetic photon in the pseudorapidity range η > 7.9 and are restricted to the range 0 < p∗T < 0.4 GeV and 0.1 < xF < 0.7. The requirement that xF is below 0.7 for photons ensures that the electromagnetic clusters reconstructed in the FNC mainly originate from single photons, according to MC predictions. The kinematic phase space of the measurements is summarised in Table 1. Cross sections of neutrons and photons produced in the forward direction, normalised to the inclusive DIS cross section, 1/σDIS dσ/dxF , are determined differentially in xF in three ranges of W . In addition, γ,n the cross section ratios integrated over xF , σDIS /σDIS , are measured as a function of W .

2.4 Event selection The data selection and analysis procedures are similar to those described in previous publications using the FNC [7, 8]. The data sample of this analysis was collected using triggers which require the scattered positron to be measured in the SpaCal. The trigger efficiency is about 96% for the analysis phase space as determined from data using an independently triggered data sample. The selection of DIS events is based on the identification of the scattered positron as the most energetic, isolated compact calorimetric deposit in the SpaCal with an energy Ee′ > 11 GeV and a polar angle 156◦ < θe′ < 175◦ . The z-coordinate of the primary event vertex is required to be within ±35 cm of the nominal interaction point. The hadronic final state 7

NC DIS Selection 6 < Q2 < 100 GeV2 0.05 < y < 0.6 70 < W < 245 GeV Forward photons

Forward neutrons

η > 7.9

η > 7.9

0.1 < xF < 0.7

0.1 < xF < 0.94

0 < p∗T < 0.4 GeV

0 < p∗T < 0.6 GeV 1

dσ

W ranges for cross sections σ DIS dxF 70 < W < 130 GeV 130 < W < 190 GeV 190 < W < 245 GeV Table 1: Definition of the kinematic phase space of the measurements. is reconstructed using an energy flow algorithm which combines charged particles measured in the trackers with information from the SpaCal and LAr calorimeters [37,38]. To suppress events with hardP initial state QED radiation, as well as events originating from non-ep interactions, the quantity E − pz , summed over all reconstructed particles including the positron, is required to lie between 35 GeV and 70 GeV. This cut also efficiently removes events from photoproduction processes, where the positron is scattered into the backward beam-pipe and a particle from the hadronic final state fakes the positron signature in the SpaCal. The kinematic variables Q2 and y are reconstructed using a technique which optimises the resolution throughout the measured y range, exploiting information from both the scattered positron and the hadronic final state [39]. Events are restricted to the range 6 < Q2 < 100 GeV2 and 0.05 < y < 0.6. The DIS data sample contains about 9.3 million events. A subsample of events containing forward photons or neutrons is selected by requiring either an electromagnetic or a hadronic cluster in the FNC with a pseudorapidity above 7.9 and an energy above 92 GeV. The data sample, called ‘FNC sample’ in the following, contains about 83, 000 events with photons and 230, 000 events with neutrons.

2.5 Monte Carlo simulations and corrections to the data Monte Carlo simulations are used to correct the data for the effects of detector acceptance, inefficiencies, QED radiation from the positron and migrations between measurement bins due to the finite detector resolution. All generated events are passed through a GEANT3 [40] based 8

simulation of the H1 apparatus and are subject to the same reconstruction and analysis chain as the data. The DJANGOH [41] program is used to generate inclusive DIS events. It is based on leading order electroweak cross sections and takes into account QCD effects up to order αs . Higher order QCD effects are simulated using leading log parton showers as implemented in LEPTO [42], or using the Colour Dipole Model (CDM) as implemented in ARIADNE [43]. Subsequent hadronisation effects are modelled using the Lund string fragmentation model as implemented in JETSET [44, 45]. Higher order electroweak processes are simulated using an interface to HERACLES [46]. The LEPTO program optionally includes the simulation of soft colour interactions (SCI) [47], in which the production of diffraction-like configurations is enhanced via non-perturbative colour rearrangements between the outgoing partons. The SCI option in LEPTO is used for the simulation of forward photons. For the DJANGOH MC simulations the H1PDF 2009 parameterisation [48] of the parton distributions in the proton is used. In the following, the DJANGOH predictions based on LEPTO and ARIADNE are denoted as LEPTO and CDM, respectively. In all DJANGOH simulations forward particles originate exclusively from the hadronisation of the proton remnant and forward photons are therefore mainly produced from the decay of π 0 mesons. RAPGAP [49] is a general purpose event generator for inclusive and diffractive ep interactions. Higher order QCD effects are simulated using parton showers and the final state hadrons are obtained via the Lund string model. As in DJANGOH higher order electroweak processes are simulated using an interface to HERACLES [46]. In the version denoted below as RAPGAP-π, the program simulates exclusively the scattering of virtual or real photons off an exchanged pion (Figure 1b). In this model the cross section for ep scattering to the final state nX takes the form dσ(ep → e′ nX) = fπ+ /p (xL , t) · dσ(eπ + → e′ X).

(2)

Here xL is the longitudinal momentum fraction and t is the squared four-momentum transfer between the incident proton and the final state neutron; fπ+ /p (xL , t) represents the pion flux associated with the splitting of a proton into a π + n system and dσ(eπ + → e′ X) is the cross section of the positron-pion interaction. There are several parameterisations of the pion flux [50–54] and the one used here is taken from [51]. The details of the pion flux parameterisation are described in [7]. Using other parameterisations of the pion flux affects mainly the absolute normalisation by up to 30 %. As was shown in [7], the best description of the forward neutron data is achieved by a combination of events with neutrons originating from pion exchange, as simulated by RAPGAPπ, and events with neutrons from proton remnant fragmentation, simulated by DJANGOH. RAPGAP-π mainly contributes at high neutron energies, while DJANGOH is significant at low energies. In [7] the contributions of RAPGAP-π and CDM were added using weighting factors 0.65 and 1.2 for the respective predictions. In the present analysis the best description of the neutron energy distribution is obtained by the combination of RAPGAP-π and CDM using the respective weights 0.6 and 1.4, or by the combination of RAPGAP-π and LEPTO using the respective weights 0.6 and 0.7. The difference between the weighting factors for the combination of RAPGAP-π and CDM in this analysis and in [7] is due to the different neutron energy range 9

and the resulting different neutron energy dependence in the two analyses. Compared to [7] the current analysis is extended to much lower neutron energies, at which the contribution from the fragmentation model dominates. The measurements are also compared to predictions of several hadronic interaction models which are commonly used for the simulation of cosmic ray air shower cascades: EPOS LHC [55, 56], QGSJET 01 [57, 58], QGSJET II-04 [59, 60] and SIBYLL 2.1 [61, 62]. These models are based on Regge theory [63], on the Reggeon calculus of Gribov [64] and on perturbative QCD. They use an unitarisation procedure to reconstruct amplitudes for exclusive processes and to determine the total and elastic cross sections. Central elements of these models are the production of mini-jets and the formation of colour strings that fragment into hadrons. Whereas the Regge-Gribov approximation is applied to hadrons as interacting objects in the case of QGSJET and SIBYLL, it is extended to include partonic constituents in EPOS LHC. Compared to the earlier EPOS simulation [55], which was used for comparison with the previous H1 forward photon analysis [8], the new EPOS LHC model [56] includes a modified treatment of central diffraction and the diffractive remnant in order to reproduce rapidity gap measurements at the LHC. The CR models also differ in the treatment of saturation effects at high parton densities at small Bjorken-x and in the treatment of the hadronic remnants in collisions. The programs are interfaced with the PHOJET program [65] for the generation of the ep scattering kinematics. It was pointed out [66] that the hadronic interaction models have been developed for hadron-hadron interactions and therefore the simulation of DIS events might be affected by the superfluous contribution of multi-parton interactions. In order to investigate this assumption, the QGSJET 01 model has been modified [67] to exclude the multi-parton interactions. In the comparison with the measurements this modified model is denoted as ‘QGSJET 01 (no mi)’. The measured distributions may contain background arising from several sources. The background from photoproduction processes is estimated using the PHOJET MC generator. It is found to be about 1% on average and is subtracted from the data distributions bin-by-bin. Background from misidentification of photons or neutrons in the FNC is estimated from the DJANGOH MC simulation to be 2% on average and is subtracted from the data distributions bin-by-bin. Background also arises from a random coincidence of DIS events, causing activity in the central detector, with a beam-related background signal in the FNC, produced from the interaction of another beam proton with a positron or with residual gas in the beampipe. This contribution is estimated by combining DIS events with FNC clusters originating from interactions in the bunch-crossings adjacent to the bunch-crossings of the DIS events. It is found to be smaller than 1% and is neglected. The MC simulations are used to correct the distributions at the level of reconstructed particles back to the hadron level on a bin-by-bin basis. The correction factors are determined from the MC simulations as the ratios of the normalised cross sections obtained from particles at hadron level without QED radiation to the normalised cross sections calculated using reconstructed particles and including QED radiation effects. For the forward photon analysis the average of the correction factors determined from LEPTO and CDM is used. For the forward neutron analysis the correction factors are calculated using the combination of RAPGAP-π and CDM simulations, with the weighting factors 0.6 and 1.4, as described above. The size of the correction factors varies between 2 and 3.5 for the forward photon and between 2 and 6 for the 10

forward neutron xF distributions, and is about 3.2 for the W distribution in both cases. The correction factors are large mainly due to the non-uniform azimuthal acceptance of the FNC, which is about 30% on average. The bin purity, defined as the fraction of events reconstructed in a particular bin that originate from the same bin on hadron level, varies between 75% and 95%.

2.6 Systematic uncertainties Several sources of experimental uncertainties are considered and their effect on the measured cross section is quantified. The systematic uncertainties on the cross section measurements are determined using MC simulations, by propagating the corresponding uncertainty through the full analysis chain. As the cross sections are normalised to the inclusive DIS cross section measured in this analysis, some important systematic uncertainties, such as those involving the trigger efficiency and the integrated luminosity and those related to the reconstruction of the scattered positron and the hadronic final state are largely reduced or cancel. The following sources are considered for both the DIS sample and the FNC samples: • The uncertainty on the measurements of the scattered positron energy (1%). • The uncertainty on the measurements of the scattered positron angle (1 mrad). • The uncertainty on the measurements of the energy of the hadronic final state (4%). • The uncertainty on the trigger efficiency (1%). These uncertainties are strongly correlated between the DIS and the forward photon and neutron samples. The resulting combined uncertainty of the cross section is about 2% on average and is considered as uncorrelated between the measurement points. Several sources of uncertainties related to the reconstruction of the forward photons and neutrons in the FNC are considered: • The acceptance of the FNC calorimeter is defined by the interaction point and the geometry of the HERA magnets and is determined using MC simulations. The uncertainty of the impact position of the particle on the FNC, due to beam inclination and the uncertainty on the FNC position, is estimated to be 5 mm. This results in uncertainties on the FNC acceptance determination of up to 15% for the xF distributions. • The absolute electromagnetic energy scale of the FNC is known to a precision of 5% [8]. This leads to an uncertainty of 1% on the forward photon cross section measurement at low energies, increasing to 11% for the largest xF values. • The uncertainty in the neutron detection efficiency and the 2% uncertainty on the absolute hadronic energy scale of the FNC [7] lead to systematic errors on the cross section of 2% and up to 10%, respectively. 11

These systematic uncertainties related to the FNC are strongly correlated between measurement bins and mainly contribute to the overall normalisation uncertainty. For the normalised cross sections studied as a function of W all above-mentioned FNC related systematic uncertainties contribute to a normalisation uncertainty of approximately 7%. In the procedure of correcting the measurements to the hadron level, using MC simulations, the following sources of systematic uncertainties are considered: • The systematic uncertainty arising from the radiative corrections and the model dependence of the data correction for the forward neutron cross section is estimated by varying the DJANGOH and RAPGAP-π scaling factors within values permitted by the data and by switching between the CDM and LEPTO models within DJANGOH. The resulting uncertainty on the cross section is typically 4%, increasing to 5% at lowest and highest xF values. • For the forward photon cross section the systematic uncertainty, taken as the difference of the acceptance corrections calculated using the LEPTO and CDM models, increases from 1% at low xF to 7% at higher xF . • The use of different parton distribution functions in the MC simulation results in a negligible change of the correction factors. • A 2% uncertainty is attributed to the bin-by-bin subtraction of background arising from the wrong identification of photon and neutron candidates in the FNC and from photoproduction processes. These uncertainties are assumed to be equally shared between correlated and uncorrelated parts. The systematic uncertainties shown in the figures and tables are calculated using the quadratic sum of all contributions, which may vary from point to point. They are significantly larger than the statistical uncertainties in all measurement bins. The total systematic uncertainty for the normalised cross section measurements ranges between 8% and 22% for the xF dependent cross sections and is about 8% for the W dependent cross sections.

3 Results 3.1 Normalised cross sections as a function of W The ratios of the forward photon and forward neutron production cross sections to the inclusive γ,n DIS cross section, σDIS /σDIS , in the kinematic range 6 < Q2 < 100 GeV2 and 0.05 < y < 0.6 and as a function of W are presented in Tables 2 and 3 and are shown in Figure 2. Within uncertainties the measured ratios are consistent with constant values of about 0.027 (forwards photons) and 0.083 (forward neutrons). In other words, within uncertainties the W dependence 12

of the cross section is independent of the presence of a forward neutron or a forward photon, as predicted by the limiting fragmentation hypothesis [9, 10]. In Figure 2 the MC model calculations are compared to the measurements. Both CDM and LEPTO predict a forward photon rate of about 70% higher than observed. A similar excess was observed earlier [8]. The photon production rate as a function of W is rather flat in CDM and shows a slight increase with W in LEPTO. The shape of the W distribution is in both models consistent with the data, within errors. The rate of forward neutron production predicted by LEPTO is consistent with the data, while CDM predicts a much lower rate. However, as was shown in the previous measurement [7], the energy distribution of forward neutrons can be described by MC simulation only if this includes contributions both from standard fragmentation as simulated in DJANGOH, and from a pion exchange mechanism as explicitly simulated in RAPGAP-π but not included in DJANGOH. In Figure 2b the combinations of the RAPGAP-π and DJANGOH simulations, as described in section 2.5, are compared to the measurement. The weighting factors 1.4 for the CDM, 0.7 for the LEPTO and 0.6 for the RAPGAP-π predictions are determined by fitting the observed neutron energy distributions integrated over the full W range. The cross sections for inclusive DIS events, used for the normalisation of the forward neutron cross sections, σDIS , are taken from the CDM and LEPTO simulations without additional weights. The model combination describes the observed W dependence well. It is remarkable that the factors for the CDM and LEPTO contributions differ by a factor two (1.4 and 0.7, respectively). It is also notable that the CDM model, which overestimates the rate of forward photons by about 70%, has to be scaled up in the combination to describe the forward neutron data. In Figure 3 predictions of various cosmic ray hadronic interaction models (EPOS LHC, SIBYLL 2.1 and the two versions of QGSJET) are compared to the measured normalised cross sections as a function of W . The CR model predictions show significant differences in absolute values, for both forward photons and forward neutrons. For photons all models predict too high rates by 30 to 40%, and these rates, with the exception of EPOS LHC, show a slight decrease with increasing W , not confirmed by data. For forward neutrons all CR predictions show a W independent behaviour, in accordance with the measured W dependence. The QGSJET 01 model predicts a much too high and SIBYLL 2.1 a much too low neutron rate, while the EPOS LHC and QGSJET II-04 models are closer to the measurement.

3.2 Normalised cross sections as a function of xF and test of Feynman scaling The measured normalised differential cross sections, 1/σDIS dσ/dxF , of the most energetic photon are presented as a function of xF in Table 4 and in Figures 4 and 5 for the kinematic region defined in Table 1. In order to study the energy dependence of the xF distributions, these cross sections are measured in three W intervals. The normalised differential cross sections as a function of xF are similar for the three W ranges. As shown in Figure 4 and already seen in the comparison of the W dependence, the 13

LEPTO and CDM models predict a rate of forward photons about 70% higher than measured. The shapes of the measured distributions are well described by LEPTO, while the CDM description is very poor by showing a significantly harder spectrum than observed in data. In Figure 5 the predictions of the CR hadronic interaction models are compared to the same measurements. Large differences between the CR models are observed, both in shape and in normalisation. All models tested here overestimate the forward photon rate by 30% to 40% at low xF . The EPOS LHC model describes the shapes of the photon xF distributions well. The SIBYLL 2.1 model predicts a harder xF dependence, while the spectra obtained from the different variants of QGSJET are softer than observed in the data. Forward photon and neutral pion measurements at the LHC also revealed differences with respect to a similar selection of CR models [68–70]. The normalised differential cross sections for forward neutrons are presented in Table 5 and in Figures 6 and 7 for the kinematic region defined in Table 1. The xF distributions are well reproduced by a combination of CDM and RAPGAP-π, using the weighting factors and normalisation as described in section 3.1. The individual contributions of the two models are shown in Figure 6 as well. Fragmentation, as simulated by CDM, dominates the neutron production at lower xF , while the contribution from pion exchange becomes significant at xF > ∼ 0.7. The combination of LEPTO and RAPGAP-π (not shown) also provides a good description of the measurements for the three W ranges. In Figure 7 the predictions of the CR hadronic interaction models are compared to the forward neutron production cross sections. The EPOS LHC model provides a reasonable description of the neutron xF distributions, except at the highest xF values. The SIBYLL 2.1 model describes the shape of the xF spectra but fails in the absolute rate. The QGSJET II-04 model shows a harder xF dependence, and QGSJET 01 predicts a much too high neutron rate. A modified version of the QGSJET 01 model, denoted ’QGSJET 01 (no mi)’ [67], in which the contribution of multi-parton interactions is excluded (see section 2.5), is also compared to the measurements. When multi-parton interactions are switched off, the predicted xF spectra become harder without improving the data description. The W dependence of the xF distributions allows a test of the Feynman scaling hypothesis for particle production. For this test, the ratios of the normalised cross sections for different CM energy intervals are studied as a function of xF . Figures 8 and 9 show the ratios of the second to the first and the third to the first W range for photons. The predictions from CDM, LEPTO and the CR models are also shown. In Figure 10 the same ratios are shown for forward neutrons and the CR models are compared to the data. For all data distributions the values of these ratios are consistent with unity and with being constant within uncertainties, suggesting that Feynman scaling in the target fragmentation region holds for photons and neutrons. The LEPTO and CDM MC models, used for the comparison with forward photon data, show a similar behaviour. All CR models indicate deviations from scaling for the forward photons, such that the production rate decreases with increasing W . In particular, this effect is strong for the SIBYLL 2.1 and QGSJET 01 models. For forward neutrons the CR models are consistent with Feynman scaling, with exception of SIBYLL 2.1.

14

4 Summary The production of high energy forward neutrons and photons has been studied at HERA in deep-inelastic ep scattering in the kinematic region 6 < Q2 < 100 GeV2 and 0.05 < y < 0.6. The normalised DIS cross sections 1/σDIS dσ/dxF for the production of photons and neutrons at pseudorapidities η > 7.9 and in the range of Feynman-x of 0.1 < xF < 0.7 for the photons and 0.1 < xF < 0.94 for neutrons are presented. The measured cross sections as a function of xF at different centre-of-mass energies of the virtual photon-proton system agree within uncertainties, confirming the validity of Feynman scaling in the energy range of the virtual photon-proton system 70 < W < 245 GeV. Different Monte Carlo models are compared to the measurements. All these models overestimate the rate of photons by 30 − 70%. The shapes of the measured forward photon cross sections are well described by the LEPTO MC simulation, while predictions based on the colour dipole model fail, especially at high xF . The cross sections for forward neutrons are well described by a linear combination of the standard fragmentation model, as implemented in DJANGOH, and the one-pion-exchange model RAPGAP-π. Predictions of models, which are commonly used for the simulation of cosmic ray cascades, are also compared to the forward photon and neutron measurements. None of the models describes the photon and neutron data simultaneously well. The best description of the shapes of the photon and the neutron xF distributions is provided by the EPOS LHC model. Within the kinematic range of the measurements, the relative rate of forward photons and neutrons in DIS events is observed to be independent of the energy of the virtual photon-proton CM, and therefore also consistent with the hypothesis of limiting fragmentation. The present measurement provides new information to further improve the understanding of proton fragmentation in collider and cosmic ray experiments.

Acknowledgements We are grateful to the HERA machine group whose outstanding efforts have made this experiment possible. We thank the engineers and technicians for their work in constructing and maintaining the H1 detector, our funding agencies for financial support, the DESY technical staff for continual assistance and the DESY directorate for support and for the hospitality which they extend to the non DESY members of the collaboration. We would like to give credit to all partners contributing to the EGI computing infrastructure for their support for the H1 Collaboration. We also wish to thank Tanguy Pierog, Ralph Engel and Sergey Ostapchenko for providing the predictions of the cosmic ray models and for fruitful discussions.

15

References [1] R. Engel, Nucl. Phys. Proc. Suppl. 75A (1999) 62 [astro-ph/9811225]. [2] A. Buniatyan et al., ”Cosmic Rays, HERA and the LHC: Working group summary.” Proceedings of the Workshop on the Implications of HERA for LHC Physics, Geneva, Switzerland, 26-30 May 2008, p.566, DESY-PROC-2009-02 [arXiv:0903.3861]. [3] C. Adloff et al. [H1 Collaboration], Eur. Phys. J. C6 (1999) 587 [hep-ex/9811013]. [4] S. Chekanov et al. [ZEUS Collaboration], Nucl. Phys. B637 3 (2002) 3 [hep-ex/0205076]. [5] S. Chekanov et al. [ZEUS Collaboration], Nucl. Phys. B776 (2007) 1 [hep-ex/0702028]. [6] S. Chekanov et al. [ZEUS Collaboration], JHEP 06 (2009) 074 [arXiv:0812.2416]. [7] F. Aaron et al. [H1 Collaboration], Eur. Phys. J. C68 (2010) 381 [arXiv:1001.0532]. [8] F. Aaron et al. [H1 Collaboration], Eur. Phys. J. C71 (2011) 1771 [arXiv:1106.5944]. [9] J. Benecke et al., Phys. Rev. 188 (1969) 2159. [10] T. T. Chou and C.-N. Yang, Phys. Rev. D50 (1994) 590. [11] R. P. Feynman, Phys. Rev. Lett. 23 (1969) 1415. [12] G. Arnison et al. [UA1 Collaboration], Phys. Lett. B118 (1982) 167. [13] G. Arnison et al. [UA1 Collaboration], Phys. Lett. B123 (1983) 108. [14] M. Banner et al. [UA2 Collaboration], Z. Phys. C27 (1985) 329. [15] D. Bernard et al. [UA4 Collaboration], Phys. Lett. B166 (1986) 459. [16] G. Alner et al. [UA5 Collaboration], Z. Phys. C33 (1986) 1. [17] G. Alner et al. [UA5 Collaboration], Phys. Rep. 154 (1987) 247. [18] T. Alexopoulos et al., Phys. Rev. Lett. 60 (1988) 1622. [19] F. Abe et al. [CDF Collaboration], Phys. Rev. D41 (1990) 2330. [20] F. Aaron et al. [H1 Collaboration], Phys. Lett. B654 (2007) 148 [arXiv:0706.2456]. [21] H. Abramowicz et al. [ZEUS Collaboration], JHEP 1006 (2010) 009 [arXiv:1001.4026]. [22] E. Pare et al. [UA7 Collaboration], Phys. Lett. B242 (1990) 531. [23] M. Albrow et al. [CHLM Collaboration], Nucl. Phys. B56 (1973) 333. [24] M. Albrow et al. [CHLM Collaboration], Nucl. Phys. B73 (1974) 40. 16

[25] P. Capiluppi et al., Nucl. Phys. B70 (1974) 1. [26] P. Capiluppi et al., Nucl. Phys. B79 (1974) 189. [27] A. Ohsawa, Prog. Theor. Phys. 92 (1994) 1005. [28] I. Abt et al. [H1 Collaboration], Nucl. Instrum. Meth. A386 (1997) 310. [29] I. Abt et al. [H1 Collaboration], Nucl. Instrum. Meth. A386 (1997) 348. [30] R. D. Appuhn et al. [H1 SPACAL Group], Nucl. Instrum. Meth. A386 (1997) 397. [31] D. Pitzl et al., Nucl. Instrum. Meth. A454 (2000) 334 [hep-ex/0002044]. [32] B. Andrieu et al. [H1 Calorimeter Group], Nucl. Instrum. Meth. A336 (1993) 460. [33] P.J. Laycock et al., JINST 7 (2012) T08003 [arXiv:1206:4068]. [34] T. Nicholls et al. [H1 SPACAL Group], Nucl. Instrum. Meth. A374 (1996) 149. [35] F. Aaron et al. [H1 Collaboration], Eur. Phys. J. C72 (2012) 2163, Erratum-ibid. C74 (2014) 2733 [arXiv:1205.2448]. [36] A. Uraev, “A new Forward Neutron Calorimeter for the H1 Experiment at the HERA collider.” (In Russian), Diploma thesis, Moscow Engineering Physics Institute (2001) (available at http://www-h1.desy.de/publications/theses list.html). [37] M. Peez, “Search for deviations from the standard model in high transverse energy processes at the electron proton collider HERA.” (In French), PhD thesis, Univ. Lyon (2003), CPPM-T-2003-04 (available at http://www-h1.desy.de/publications/theses list.html). [38] S. Hellwig, “Investigation of the D∗ − πslow double tagging method in charm analyses.” (In German), Diploma thesis, Univ. Hamburg (2004) (available at http://wwwh1.desy.de/publications/theses list.html). [39] C. Adloff et al. [H1 Collaboration], Z. Phys. C76 (1997) 613 [hep-ex/9708016]. [40] R. Brun et al., GEANT3, CERN-DD/EE/84-1. [41] K. Charchula, G. A. Schuler and H. Spiesberger, DJANGOH 1.4, Comput. Phys. Commun. 81 (1994) 381. [42] G. Ingelman, A. Edin and J. Rathsman, LEPTO 6.5, Comput. Phys. Commun. 101 (1997) 108 [hep-ph/9605286]. [43] L. Lönnblad, ARIADNE 4.10, Comput. Phys. Commun. 71 (1992) 15. [44] B. Andersson et al., Phys. Rep. 97 (1983) 31. [45] T. Sjöstrand, PYTHIA 5.7 and JETSET 7.4, hep-ph/9508391. 17

[46] A. Kwiatkowski, H. Spiesberger and H. J. Möhring, Comp. Phys. Commun. 69 (1992) 155. [47] A. Edin, G. Ingelman and J. Rathsman, Phys. Lett. B366 (1996) 371 [hep-ph/9508386]. [48] F. D. Aaron et al. [H1 Collaboration], Eur. Phys. J. C64 (2009) 561 [arXiv:0904.3513]. [49] H. Jung, RAPGAP 3.1, Comp. Phys. Commun. 86 (1995) 147. [50] M. Bishari, Phys. Lett. B38 (1972) 510. [51] H. Holtmann et al., Phys. Lett. B338 (1994) 363. [52] B. Kopeliovich, B. Povh and I. Potashnikova, Z. Phys. C73 (1996) 125 [hep-ph/9601291]. [53] M. Przybycień, A. Szczurek and G. Ingelman, Z. Phys. C74 (1997) 509 [hep-ph/9606294]. [54] A. Szczurek, N. N. Nikolaev and J. Speth, Phys. Lett. B428 (1998) 383 [hep-ph/9712261]. [55] K. Werner, F.-M. Liu and T. Pierog, Phys. Rev. C74 (2006) 044902 [hep-ph/0506232]. [56] T. Pierog et al., DESY-13-125, arXiv:1306.0121. [57] N. N. Kalmykov and S. S. Ostapchenko, Phys. Atom. Nucl. 56 (1993) 346. [58] N. N. Kalmykov, S. S. Ostapchenko and A. I. Pavlov, Nucl. Phys. Proc. Suppl. 52B (1997) 17. [59] S. S. Ostapchenko, Phys. Rev. D74 (2006) 014026 [hep-ph/0505259]. [60] S. S. Ostapchenko, Phys. Rev. D83 (2011) 014018 [arXiv:1010.1869]. [61] J. Engel et al., Phys. Rev. D46 (1992) 5013. [62] E.-J. Ahn et al., Phys. Rev. D80 (2009) 094003 [arXiv:0906.4113]. [63] P. D. B. Collins, “An Introduction to Regge Theory and High-Energy Physics.” Cambridge 1977, 445p. [64] V. N. Gribov, Sov. Phys. JETP 26 (1968) 414. [65] R. Engel and J. Ranft, Phys. Rev. D54 (1996) 4244 [hep-ph/9509373]. [66] S. S. Ostapchenko, Private communication. [67] S. S. Ostapchenko and T. Pierog, Private communication. [68] O. Adriani et al. [LHCf Collaboration], Phys. Lett. B703 (2011) 128 [arXiv:1104.5294]. [69] O. Adriani et al. [LHCf Collaboration], Phys. Rev. D86 (2012) 092001 [arXiv:1205.4578]. [70] O. Adriani et al. [LHCf Collaboration], Phys. Lett. B715 (2012) 298 [arXiv:1207.7183]. 18

correlated sys. uncertainty

W range [GeV ]

γ σDIS (W ) σDIS (W )

δstat.

δtotal sys.

δuncorrel.sys.

δ EF N C

δXYF N C

δmodel

70. ÷ 115. 115. ÷ 160. 160. ÷ 205. 205. ÷ 245.

0.0269 0.0269 0.0265 0.0265

0.0002 0.0002 0.0002 0.0002

0.0022 0.0022 0.0022 0.0022

0.0006 0.0007 0.0007 0.0007

0.0011 0.0011 0.0011 0.0011

0.0018 0.0018 0.0018 0.0018

0.0003 0.0003 0.0003 0.0003

Table 2: The fraction of DIS events with forward photons in the kinematic region given in Table 1. For each measurement, the statistical (δstat. ), the total systematic (δtotal sys. ), the uncorrelated (δuncorrel.sys. ) systematic uncertainties and the bin-to-bin correlated systematic uncertainties due to the FNC absolute energy scale (δEF N C ), the measurement of the particle impact position in the FNC (δXYF N C ) and the model dependence of the data correction (δmodel ) are given. correlated sys. uncertainty

W range [GeV ]

n (W ) σDIS σDIS (W )

δstat.

δtotal sys.

δuncorrel.sys.

δ EF N C

δXYF N C

δmodel

70. ÷ 115. 115. ÷ 160. 160. ÷ 205. 205. ÷ 245.

0.0843 0.0830 0.0815 0.0826

0.0004 0.0004 0.0005 0.0006

0.0074 0.0074 0.0072 0.0073

0.0020 0.0021 0.0020 0.0022

0.0008 0.0008 0.0008 0.0008

0.0057 0.0056 0.0055 0.0055

0.0042 0.0042 0.0041 0.0041

Table 3: The fraction of DIS events with forward neutrons in the kinematic region given in Table 1. For each measurement, the statistical (δstat. ), the total systematic (δtotal sys. ), the uncorrelated (δuncorrel.sys. ) systematic uncertainties and the bin-to-bin correlated systematic uncertainties due to the FNC absolute energy scale (δEF N C ), the measurement of the particle impact position in the FNC (δXYF N C ) and the model dependence of the data correction (δmodel ) are given.

19

W = 70 − 130 GeV xF range

dσ σDIS dxF 1

0.10 ÷ 0.22 0.22 ÷ 0.34 0.34 ÷ 0.46 0.46 ÷ 0.58 0.58 ÷ 0.70

0.130 0.0542 0.0221 0.00743 0.00202

δstat.

δtotal sys.

δuncorrel.sys.

δ EF N C

δXYF N C

δmodel

0.001 0.0007 0.0005 0.00024 0.00010

0.011 0.0060 0.0031 0.00122 0.00044

0.003 0.0015 0.0007 0.00032 0.00016

0.001 0.0027 0.0018 0.00059 0.00022

0.011 0.0051 0.0024 0.00099 0.00031

0.001 0.0008 0.0003 0.00026 0.00014


1

dσ σDIS dxF

0.10 ÷ 0.22 0.22 ÷ 0.34 0.34 ÷ 0.46 0.46 ÷ 0.58 0.58 ÷ 0.70

0.128 0.0553 0.0222 0.00724 0.00192

xF range 0.10 ÷ 0.22 0.22 ÷ 0.34 0.34 ÷ 0.46 0.46 ÷ 0.58 0.58 ÷ 0.70

dσ σDIS dxF 0.124 0.0568 0.0222 0.00754 0.00190


δstat.

δtotal sys.

δuncorrel.sys.

δ EF N C

δXYF N C

δmodel

0.001 0.0008 0.0005 0.00027 0.00011

0.011 0.0063 0.0031 0.00120 0.00041

0.003 0.0016 0.0007 0.00032 0.00015

0.001 0.0028 0.0018 0.00058 0.00021

0.011 0.0053 0.0024 0.00097 0.00029

0.001 0.0008 0.0003 0.00025 0.00013

W = 190 − 245 GeV

1



δstat.

δtotal sys.

δuncorrel.sys.

δ EF N C

δXYF N C

δmodel

0.001 0.0010 0.0006 0.00034 0.00014

0.011 0.0064 0.0031 0.00125 0.00041

0.003 0.0017 0.0007 0.00033 0.00015

0.001 0.0028 0.0018 0.00060 0.00021

0.011 0.0054 0.0024 0.00101 0.00029

0.001 0.0008 0.0003 0.00026 0.00013

Table 4: Normalised cross sections of forward photon production in DIS as a function of xF . The kinematic phase space of the measurements is given in Table 1. For each measurement, the statistical (δstat. ), the total systematic (δtotal sys. ), the uncorrelated (δuncorrel.sys. ) systematic uncertainties and the bin-to-bin correlated systematic uncertainties due to the FNC absolute energy scale (δEF N C ), the measurement of the particle impact position in the FNC (δXYF N C ) and the model dependence of the data correction (δmodel ) are given.

20

W = 70 − 130 GeV xF range 0.10 ÷ 0.22 0.22 ÷ 0.34 0.34 ÷ 0.46 0.46 ÷ 0.58 0.58 ÷ 0.70 0.70 ÷ 0.82 0.82 ÷ 0.94

1

dσ σDIS dxF 0.0456 0.0823 0.1096 0.1309 0.1407 0.1266 0.0656

δstat.

δtotal sys.

δuncorrel.sys.

δ EF N C

δXYF N C

δmodel

0.0015 0.0016 0.0016 0.0016 0.0015 0.0013 0.0008

0.0042 0.0079 0.0114 0.0151 0.0199 0.0179 0.0096

0.0012 0.0024 0.0033 0.0056 0.0108 0.0069 0.0036

0.0023 0.0049 0.0077 0.0105 0.0127 0.0127 0.0066

0.0023 0.0044 0.0064 0.0076 0.0088 0.0085 0.0050

0.0023 0.0037 0.0044 0.0053 0.0063 0.0063 0.0033

W = 130 − 190 GeV xF range 0.10 ÷ 0.22 0.22 ÷ 0.34 0.34 ÷ 0.46 0.46 ÷ 0.58 0.58 ÷ 0.70 0.70 ÷ 0.82 0.82 ÷ 0.94

dσ σDIS dxF 1

0.0426 0.0801 0.1077 0.1286 0.1359 0.1224 0.0617

0.10 ÷ 0.22 0.22 ÷ 0.34 0.34 ÷ 0.46 0.46 ÷ 0.58 0.58 ÷ 0.70 0.70 ÷ 0.82 0.82 ÷ 0.94


δstat.

δtotal sys.

δuncorrel.sys.

δ EF N C

δXYF N C

δmodel

0.0017 0.0019 0.0019 0.0018 0.0017 0.0014 0.0009

0.0038 0.0077 0.0112 0.0148 0.0192 0.0172 0.0090

0.0010 0.0023 0.0032 0.0055 0.0105 0.0066 0.0033

0.0021 0.0048 0.0075 0.0103 0.0122 0.0122 0.0062

0.0021 0.0043 0.0063 0.0075 0.0085 0.0082 0.0047

0.0021 0.0036 0.0043 0.0051 0.0061 0.0061 0.0031




1 dσ σDIS dxF

qδstat.

δtotal sys.

δuncorrel.sys.

δ EF N C

δXYF N C

δmodel

0.0454 0.0796 0.1093 0.1273 0.1357 0.1250 0.0621

0.0022 0.0024 0.0024 0.0023 0.0021 0.0018 0.0011

0.0042 0.0077 0.0114 0.0146 0.0191 0.0176 0.0090

0.0012 0.0023 0.0033 0.0054 0.0104 0.0067 0.0033

0.0023 0.0048 0.0077 0.0102 0.0122 0.0125 0.0062

0.0023 0.0043 0.0064 0.0074 0.0085 0.0084 0.0047

0.0023 0.0036 0.0044 0.0051 0.0061 0.0062 0.0031

Table 5: Normalised cross sections of forward neutron production in DIS as a function of xF . The kinematic phase space of the measurements is given in Table 1. For each measurement, the statistical (δstat. ), the total systematic (δtotal sys. ), the uncorrelated (δuncorrel.sys. ) systematic uncertainties and the bin-to-bin correlated systematic uncertainties due to the FNC absolute energy scale (δEF N C ), the measurement of the particle impact position in the FNC (δXYF N C ) and the model dependence of the data correction (δmodel ) are given.

21

H1

0.05

n / σ σDIS DIS

γ σDIS / σDIS

Forward Photons

0.3

0.25

(a)

0.04 0.03 0.02

0.2

0.15 H1 Data

Forward Neutrons

H1

H1 Data correlated uncertainty LEPTO

(b)

CDM RAPGAP-π LEPTO × 0.7 + RAPGAP-π × 0.6 CDM × 1.4 + RAPGAP-π × 0.6

0.1

correlated uncertainty

0.01

0.05

LEPTO CDM

0 100

150

0

200

100

150

W [GeV]

200

W [GeV]

Figure 2: The fraction of DIS events with forward photons (a) and forward neutrons (b) as a function of W in the kinematic region given in Table 1. Also shown are the predictions of the LEPTO (solid line) and CDM (dashed line) MC models. In the case of forward neutron production, the predictions of RAPGAP-π and the linear combinations of LEPTO and RAPGAP-π, as well as CDM and RAPGAP-π are also shown.

H1

0.05

(a)

0.04

n / σ σDIS DIS

γ σDIS / σDIS

Forward Photons

Forward Neutrons

0.3

0.2

H1

H1 Data correlated uncertainty SIBYLL 2.1 EPOS LHC QGSJET II-04 QGSJET 01

(b)

0.03 0.02

H1 Data correlated uncertainty

0.01

SIBYLL 2.1 EPOS LHC QGSJET II-04 QGSJET 01

0.1

0

0 100

150

200

100

W [GeV]

150

200

W [GeV]

Figure 3: The fraction of DIS events with forward photons (a) and forward neutrons (b) as a function of W in the kinematic region given in Table 1. Also shown are the predictions of the cosmic ray hadronic interaction models SIBYLL 2.1 (solid line), QGSJET 01 (dashed line), QGSJET II-04 (dotted line) and EPOS LHC (dash-dotted line).

22

Forward Photons

70 < W < 130 GeV

1/σDIS dσ/dx F

1/σDIS dσ/dx F

Forward Photons

H1 -1

10

10-2

H1 -1

10

10-2

H1 Data -3

10

H1 Data -3

10

LEPTO CDM

0.1

0.2

LEPTO CDM

0.3

0.4

0.5

0.6

0.7

0.1

xF Forward Photons

1/σDIS dσ/dx F

130 < W < 190 GeV

0.2

0.3

0.4

0.5

0.6

0.7

xF

190 < W < 245 GeV

H1 -1

10

10-2

H1 Data -3

10

LEPTO CDM

0.1

0.2

0.3

0.4

0.5

0.6

0.7

xF Figure 4: Normalised cross sections of forward photon production in DIS as a function of xF in three W intervals in the kinematic region given in Table 1. The error bars show the total experimental uncertainty, calculated using the quadratic sum of the statistical and systematic uncertainties. Also shown are the predictions of the LEPTO (solid line) and CDM (dashed line) MC models.

23

70 < W < 130 GeV

H1 -1

10

MC/Data

1/σDIS dσ/dx F

Forward Photons

6

Forward Photons

70 < W < 130 GeV

H1

H1 Data SIBYLL 2.1

5

EPOS LHC QGSJET II-04

4

QGSJET 01 H1 Data

10-2

QGSJET 01 (no mi)

3

SIBYLL 2.1 EPOS LHC

2

QGSJET II-04

10-3

QGSJET 01

1

QGSJET 01 (no mi)

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.1

0.2

0.3

0.4

0.5

0.6

xF 130 < W < 190 GeV

H1 10-1

MC/Data

1/σDIS dσ/dx F

Forward Photons

0.7

xF 6

Forward Photons

130 < W < 190 GeV

H1

H1 Data SIBYLL 2.1

5

EPOS LHC QGSJET II-04

4

QGSJET 01 H1 Data

-2

10

QGSJET 01 (no mi)

3

SIBYLL 2.1 EPOS LHC

2

QGSJET II-04 -3

QGSJET 01

10

1

QGSJET 01 (no mi)

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.1

0.2

0.3

0.4

0.5

0.6

xF 190 < W < 245 GeV

H1 -1

10

MC/Data

1/σDIS dσ/dx F

Forward Photons

0.7

xF 6

Forward Photons

190 < W < 245 GeV

H1

H1 Data

5

SIBYLL 2.1 EPOS LHC

4

QGSJET II-04 QGSJET 01

-2

10

H1 Data

3

SIBYLL 2.1

QGSJET 01 (no mi)

EPOS LHC

2

QGSJET II-04 -3

10

QGSJET 01

1

QGSJET 01 (no mi)

0.1

0.2

0.3

0.4

0.5

0.6

0.7

xF

0.1

0.2

0.3

0.4

0.5

0.6

0.7

xF

Figure 5: Normalised cross sections of forward photon production in DIS as a function of xF in three W intervals in the kinematic region given in Table 1. The inner error bars show the statistical uncertainty, while the outer error bars show the total experimental uncertainty, calculated using the quadratic sum of the statistical and systematic uncertainties. Also shown are the predictions of the cosmic ray hadronic interaction models SIBYLL 2.1 (solid line), QGSJET 01 (dashed line), QGSJET 01 (no mi) (dash-double dotted line), QGSJET II-04 (dotted line) and EPOS LHC (dash-dotted line). In the right column the ratios of the CR model predictions to the data are shown. 24

Forward Neutrons

70 < W < 130 GeV

H1 Data CDM × 1.4 + RAPGAP-π × 0.6

0.15

1/σDIS dσ/dx F

1/σDIS dσ/dx F

0.2

H1

CDM

0.1

Forward Neutrons

130 < W < 190 GeV


0.15

RAPGAP-π

0.05

1/σDIS dσ/dx F

0.2

H1

CDM RAPGAP-π

0.1

0.05

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

xF

xF

0.2

Forward Neutrons

190 < W < 245 GeV


0.15

H1

CDM RAPGAP-π

0.1

0.05

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

xF Figure 6: Normalised cross sections of forward neutron production in DIS as a function of xF in three W intervals in the kinematic region given in Table 1. The inner error bars show the statistical uncertainty, while the outer error bars show the total experimental uncertainty, calculated using the quadratic sum of the statistical and systematic uncertainties. Also shown are the predictions of CDM (dotted line), RAPGAP-π (dashed line) and a linear combination of CDM and RAPGAP-π predictions (solid line).

25

0.3

H1

MC/Data

1/σDIS dσ/dx F

0.4

Forward Neutrons 70 < W < 130 GeV H1 Data SIBYLL 2.1 EPOS LHC QGSJET II-04 QGSJET 01 QGSJET 01 (no mi)


H1

3

2 0.2 1

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

xF

xF

0.4

0.3


H1

MC/Data

1/σDIS dσ/dx F

0.1


H1

3

2 0.2 1

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

xF

xF

0.4

0.3


H1

MC/Data

1/σDIS dσ/dx F

0.1


H1

3

2 0.2 1 0.1

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

xF

xF

Figure 7: Normalised cross sections of forward neutron production in DIS as a function of xF in three W intervals in the kinematic region given in Table 1. The inner error bars show the statistical uncertainty, while the outer error bars show the total experimental uncertainty, calculated using the quadratic sum of the statistical and systematic uncertainties. Also shown are the predictions of the cosmic ray hadronic interaction models SIBYLL 2.1 (solid line), QGSJET 01 (dashed line), QGSJET 01 (no mi) (dash-double dotted line), QGSJET II-04 (dotted line) and EPOS LHC (dash-dotted line). In the right column the ratios of the CR model predictions to the data are shown. 26

Forward Photons

1.1

H1

(a)

1

0.9

0.6

1.1

(b)

1

0.8

H1 Data

0.7

LEPTO

0.6

CDM

0.5 0.1

H1

0.9

0.8 0.7

Forward Photons

1.2

1/σDIS dσ/dx F(190