1 Introduction

In Quantum Chromodynamics (QCD), the final state of an inelastic hadron-hadron collision can be described by contributions from hard and soft scattering occurring between the constituents of the hadrons, initial- and final-state (gluon) radiation and the fragmentation of the initially coloured partonic final state into colour-neutral hadrons. The soft component of a collision is called the underlying event. Its precise theoretical description remains a challenge, while the dynamics of hard scattering processes is well described by perturbative QCD. One source of the underlying event activity is multi-parton interactions (MPI). These arise mainly in the region of a very low parton momentum fraction, where parton densities are high so that the probability of more than a single parton-parton interaction per hadron-hadron collision is large. MPI effects become increasingly important at LHC collision energies, where inelastic interactions between very soft partons are sufficiently energetic to contribute to final state particle production [1].

MPI phenomena can be probed by measuring in the centre-of-mass system the amount of energy created in inelastic hadron-hadron interactions at large values of the pseudorapidity η=−ln[tan(θ/2)], with θ being the polar angle of particles with respect to the beam axis. The energy flow is expected to be directly sensitive to the amount of parton radiation and MPI [2]. For a particular pseudorapidity interval with width Δη, the total energy flow, which is normalised to the number of inelastic pp interactions N int, is defined as

$$ \frac{1}{N_{\mathrm{int}}} \frac{dE_{\mathrm{total}}}{d\eta} = \frac{1}{\Delta\eta} \Biggl( \frac{1}{N_{\mathrm{int}}}\sum_{i=1}^{N_{\mathrm{part},\eta}}E_{i,\eta} \Biggr) , $$
(1)

where N part,η is the total number of stable particles and E i,η is the energy of the individual particles.

In this study, the energy flow is measured in pp collisions at \(\sqrt{s}=7\ \mbox{TeV}\) within the pseudorapidity range 1.9<η<4.9. This extends the previous measurements that have been made in \(p\bar{p}\) [3] and ep collisions [4] to larger pseudorapidity values and higher centre-of-mass energies, and complements the studies performed by the CMS and ATLAS collaborations [5, 6]. Experimental results are compared to predictions given by Pythia-based [7, 8] and cosmic-ray event generators [9, 10], which model the underlying event activity in different ways. In order to probe various aspects of multi-particle production in high-energy hadron-hadron collisions, the measurements are performed for the following four classes of events: inclusive minimum-bias, hard scattering, diffractive, and non-diffractive enriched interactions.

2 The LHCb detector

The LHCb detector [11] is a single-arm forward spectrometer with an angular coverage from 10 mrad to 300 (250) mrad in the bending (non-bending) plane, designed for the study of b- and c-hadrons. The detector includes a high precision tracking system consisting of a silicon-strip vertex detector (VELO) surrounding the pp interaction region, a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes placed downstream. The VELO has a larger angular acceptance than the rest of the spectrometer, including partial coverage of the backward region. It allows reconstruction of charged particle tracks in the pseudorapidity ranges 1.5<η<5.0 and −4<η<−1.5. The combined tracking system has a momentum resolution Δp/p that varies from 0.4 % at 5 GeV/c to 0.6 % at 100 GeV/c, and an impact parameter resolution of 20 μm for tracks with high transverse momentum, p T. Charged hadrons are identified using two ring-imaging Cherenkov detectors. Photon, electron and hadron candidates are distinguished by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter (ECAL) and a hadronic calorimeter (HCAL). The calorimeters have an energy resolution of \(\sigma(E)/E=10\ \%/\sqrt{E}\oplus 1\ \%\) and \(\sigma(E)/E=69\ \%/\sqrt{E}\oplus 9\ \%\) (with E in GeV), respectively. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers.

The trigger consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage which applies a full event reconstruction. For the minimum-bias data used in this analysis, the hardware trigger was accepting all beam–beam crossings, while the presence of at least one reconstructed track was required in the software stage to record an event.

3 Data analysis

3.1 Data and Monte Carlo samples

The analysis is performed using a sample of minimum-bias data collected in pp collisions at \(\sqrt{s}=7\ \mbox{TeV}\) during the initial running period of the LHC with low interaction rate. The fraction of bunch crossings with two or more collisions (“pile-up events”) is estimated to be approximately 5 %. The total number of events available in the sample is 5.8×106, corresponding to an integrated luminosity of about 0.1 nb−1.

Fully simulated minimum-bias pp events at \(\sqrt{s}=7\ \mbox{TeV}\) were generated using the LHCb tune [12] of the Pythia 6.4 event generator [7]. Here, decays of hadronic particles are described by EvtGen [13] in which final state QED radiation is generated using Photos [14]. The interaction of the generated particles with the detector and its response are implemented using the Geant4 toolkit [15, 16] as described in Ref. [17]. Additional Monte Carlo (MC) samples with fully simulated minimum-bias pp interactions at \(\sqrt{s}=7\ \mbox{TeV}\) were generated using the Perugia 0 and Perugia NOCR [18] tunes of Pythia 6.4. These models along with the LHCb tune use different values for the MPI energy scaling parameter and MPI p T cut-off, which entails a sizeable deviation in the amount of MPI predicted by these tunes.

The LHCb tune utilises the CTEQ6L parton density functions (PDFs) [19], while both Perugia tunes use the CTEQ5L PDFs [20]. Colour reconnection effects are not included in the Perugia NOCR tune. In the MC samples generated with the Perugia 0 and Perugia NOCR tunes, diffractive pp interactions are not included, whereas the sample generated with the LHCb tune contains the contributions from both single and double diffractive processes. A sample of fully simulated diffractive events generated with Pythia 8.130 [8], which utilises the CTEQ5L PDFs, is used in addition. This event generator gives a more accurate description of diffractive pp interactions than Pythia 6, especially at high-p T, as it includes the contribution from hard diffractive processes, which is absent in Pythia 6 [21].

In addition to the models above, experimental results are compared to generator level predictions given by the Pythia 8.135 model with default parameters. Furthermore, the measurements are compared with predictions given by the cosmic-ray interaction models Epos 1.99 [22], Qgsjet01, QgsjetII-03 [23], and Sibyll 2.1 [24], which are widely used in extensive air shower simulations and are not tuned to LHC data. These generate inelastic pp interactions taking into account the contributions from both soft and hard scattering processes. Soft contributions are described with Gribov’s Reggeon field theory [25] via exchanges between virtual quasi-particle states (Pomerons), while hard processes are described by perturbative QCD via exchanges of hard or semi-hard Pomerons. The predictions given by these models diverge mainly because of different treatments of non-linear interaction effects related to parton saturation [26] and shadowing [27]. The Qgsjet01 model describes hadronic multiple scattering processes as multiple exchanges of Pomerons without specific treatment of saturation effects. A distinct feature of the QgsjetII model is the treatment of non-linear parton effects via Pomeron interactions taking into account all order re-summation of the corresponding Reggeon field theory diagrams. Based on the dual parton model [28], Sibyll utilises the Lund string model [29] for hadronisation and describes soft and hard processes using the Pomeron formalism and the minijet model [30], correspondingly. The treatment of non-linear effects in this model is based on a simple geometrical approach of parton saturation. The Epos model takes into account energy-momentum correlations between multiple re-scatterings and describes non-linear effects using an effective treatment of lowest order Pomeron–Pomeron interaction graphs. It also accounts for the final state interaction of the produced particles.

3.2 Analysis strategy

The energy flow, as defined in Eq. (1), is the energy-weighted pseudorapidity distribution of particles, normalised to the number of inelastic interactions and the η-bin size. The measurements are performed in ten equidistant pseudorapidity bins of width Δη=0.3 over the range 1.9<η<4.9. The primary measurement is the energy flow carried by charged particles (charged energy flow). It is performed with reconstructed tracks which contain hits in the VELO and downstream tracking stations and have momentum in the range 2<p<1000 GeV/c. Particle identification is not required in this analysis, as the energy is taken from the momentum, neglecting particle masses. In order to be able to compare the results of the measurements with generator level predictions, the reconstructed charged energy flow is corrected for detector effects. The total energy flow is determined by using a data-constrained MC estimate of the neutral component based on information from the ECAL, while the HCAL is not used. Details of the procedure are discussed below.

3.3 Event classes

The event classes studied in this analysis are defined as follows. Inclusive minimum-bias events are selected by requesting the presence of at least one track originating from the luminous region in order to suppress pollution from beam–gas interactions and beam halo related background. Events with two or more reconstructed primary vertices are rejected to suppress pile-up contamination. To minimise biases on the track multiplicity of the event, the information on the primary vertex is not used. The selected inclusive minimum-bias interactions are further classified as hard scattering, diffractive and non-diffractive enriched events using the following criteria:

  • Hard scattering events: at least one track with p T>3 GeV/c and 1.9<η<4.9.

  • Diffractive enriched events: no tracks reconstructed with −3.5<η<−1.5.

  • Non-diffractive enriched events: at least one track reconstructed with −3.5<η<−1.5.

The selection requirements applied for the last two event classes are motivated by the fact that a sizeable rapidity gap is an experimental signature of diffractive processes [31]. The level of enrichment of the diffractive and non-diffractive samples was studied in simulation, by retrieving the Pythia process type of the pp interaction for every selected diffractive and non-diffractive candidate. In the case of the LHCb tune of Pythia 6.4, the purities of the selected diffractive and non-diffractive enriched samples are found to be about 70 % and 90 %, respectively. Although the actual percentages are only meaningful within the specific model, the study shows that the applied selection criteria indeed lead to sizeable enhancement of the respective event classes.

To minimise the experimental corrections, the definition of the event classes at generator level is similar to that at detector level. Inclusive minimum-bias events at generator level are selected by requiring the presence of at least one outgoing final-state charged particle (lifetime τ>10−8 s) in the pseudorapidity range 1.9<η<4.9, but without imposing any condition on its energy. The sample of hard scattering events is selected by requesting at least one final-state charged particle with p T>3 GeV/c and 1.9<η<4.9. The absence or presence of at least one final-state charged particle in −3.5<η<−1.5 is used as criterion to select diffractive and non-diffractive enriched events among inclusive minimum-bias interactions, respectively. For the selected events, the energy flow at generator level is determined using the outgoing final-state charged and neutral particlesFootnote 1 which are either prompt, originating directly from the fragmentation, or the decay products of unstable particles. Since neutrinos are not reconstructed by the LHCb spectrometer the energy carried by these particles is not taken into account. Only MC events simulated with exactly one inelastic pp interaction are considered in this study.

3.4 Corrections

The reconstructed charged energy flow measured with data is corrected for detector effects using bin-by-bin correction factors, which are estimated as the ratio of the charged energy flow at generator and detector level in simulation for each η region and event class under consideration. The overall correction factor for each bin is taken as the average of the correction factors obtained with different MC models used in this analysis. For inclusive, hard scattering and non-diffractive enriched events, the average and standard deviation of the correction factors, which is included in the model-dependent systematic uncertainty, are determined from the LHCb, Perugia 0 and Perugia NOCR tunes of Pythia 6.4. In the case of the diffractive enriched event class, only the LHCb tune and the Pythia 8 diffractive simulation are used. Except for the lowest η bin, which suffers from reduced acceptance for low-p T particles and thus exhibits large corrections and a sizeable model dependence, the correction factors are found to be stable among the models with a slight rise towards the edges of the detector acceptance. The majority of the factors are well below two, indicating that most of the energy is measured by the detector. In the case of diffractive enriched events, the correction factors obtained with the LHCb tune are slightly smaller than unity for some of the bins, i.e. the energy flow at detector level is found to be larger than at generator level. This is due to detection inefficiency for charged particles over the pseudorapidity range −3.5<η<−1.5. As a result, some of the events containing backward going charged particles migrate into the diffractive sample, which leads to enhanced energy flow at detector level.

For the measurement of the total energy flow, the neutral component F neut,η is estimated in the following way. To first order F neut,η is assumed to be proportional to the corrected charged energy flow F char,η with a factor R gen,η , which is the average ratio of the neutral energy flow to the charged energy flow obtained at generator level for each η bin and event class with different Pythia tunes. This ratio is found to be rather stable over the entire pseudorapidity range of the measurements with only small variations between the Pythia tunes. This reflects the usage of the same hadronisation mechanism governed in the Pythia generator by the Lund string model [7, 8]. The latter successfully describes the hadronisation of quarks and gluons emerging from high energy interactions and was rigorously tested for high-p T processes [3234]. The R gen,η ratio is found to be around 0.6 for all event types except the hard scattering interactions. For the latter, it is about 15 % smaller for all η bins. This feature is found to be a consequence of the requirement of a high-p T charged particle in the definition of this event class.

Under the assumption outlined above, the total energy flow for a particular event class and pseudorapidity bin F total,η can be written as

$$ F_{\mathrm{total},\eta} = F_{\mathrm{char},\eta}+F_{\mathrm{neut},\eta} = F_{\mathrm{char},\eta} \times (1+R_{\mathrm{gen},\eta} ). $$
(2)

In order to constrain this initially purely model-based estimate of the neutral energy flow to data, the total energy flow is further multiplied by an additional correction factor k η . It accounts for differences between simulation and data being defined for every η bin as

$$ k_{\eta} = \frac{1+R_{\mathrm{data},\eta}}{1+R_{\mathrm{mc},\eta}} . $$
(3)

Here, R data,η and R mc,η are ratios of the uncorrected neutral to charged energy flow measured in data and simulation, respectively. The R mc,η ratio is obtained with different Pythia tunes and its average is taken for the estimation of the k η factors. The neutral component of R data,η and R mc,η is measured using reconstructed photon candidates which are selected from neutral clusters in the ECAL with an energy greater than 2 GeV and p T>0.2 GeV/c. Since the polar angular coverage of the ECAL begins at about 30 mrad, there are no measurements of the neutral energy for the last two η bins (η>4.3). The k η factors for this pseudorapidity region are estimated using a linear extrapolation of the k η factors obtained for the pseudorapidity interval 3.1<η<4.3. The bins with η<3.1 are not considered for the extrapolation, as these are affected by the detection inefficiency for low-p T charged particles. The latter have a low average momentum in this η region and thus are unlikely to reach downstream tracking stations. Except for the lowest η bin, which suffers most from the detection inefficiency especially in the case of the diffractive enriched event class, the k η factors are found to be rather close to unity, reflecting the fact that the ratio of the neutral to charged energy flow is well simulated at detector level.

4 Systematic uncertainties

The total uncertainties on the results are dominated by systematic effects, as the statistical uncertainties are found to be negligible for all η bins and event classes. The various contributions to the systematic uncertainties are summarised in Table 1.

Table 1 Relative systematic uncertainties (in percent) affecting the energy flow measurements for all event classes. The total uncertainties are obtained by adding the individual sources in quadrature. The ranges indicate the variation of the uncertainty as a function of η

For all event types except hard scattering interactions, the largest uncertainty on the charged energy flow arises from the model dependence of the bin-by-bin correction factors, which is estimated as the standard deviation of the correction factors obtained with different Pythia tunes. Here, the largest impact is at low η, reaching 9 % for inclusive and non-diffractive enriched events, 4 % for hard scattering interactions and up to 43 % for diffractive enriched events. At large η this effect generally drops to about 1–2 % for all event classes except diffractive enriched interactions for which it stays above 15 %.

Systematic uncertainties related to the track selection requirements are estimated by comparing the fraction of the energy flow from tracks which are rejected by the selection cuts in data and simulation. For the majority of the bins the resulting systematic uncertainty is found to be less than 4 %. Only for hard scattering events this uncertainty approaches 9 % at low η.

To account for differences between the true tracking efficiency and that estimated using simulation, a global 3 % systematic uncertainty is assigned across the entire η range following the analysis presented in Ref. [35]. This applies for all event classes under consideration.

The other tracking related factors having an influence on the charged energy flow measurements are contaminations from multiply reconstructed tracks and tracks created from random combinations of hits (spurious tracks). The impact of the former is estimated by removing from the measurement all tracks found within the same event with similar momentum vectors. It is observed that the charged energy flow drops by less than 1 % for all η bins and event classes in case of both data and simulation. For the final results, a global 1 % systematic uncertainty for multiply reconstructed tracks is conservatively assigned. The effect of spurious tracks is estimated in simulation by determining the energy flow carried by reconstructed tracks which cannot be associated with particles at generator level and accounting for the difference between the rate of spurious tracks in data and simulation. The corresponding systematic uncertainty is found to vary between 0.2 % and 2 %.

It has been checked that reversing the LHCb magnet polarity has only an influence on the measurements of the charged energy flow for the diffractive enriched event class, which mainly consists of low-multiplicity events. Here, the corresponding effect is assigned as a systematic uncertainty.

Events with more than one reconstructed primary vertex are vetoed in the analysis in order to suppress pile-up contamination. Its residual effect is estimated to be 1.7 % by taking the efficiencies to accept pile-up events from simulation. This factor is included in the normalisation of the energy flow and conservatively taken as the systematic uncertainty.

The total energy flow acquires an additional uncertainty from the variation of the R gen,η and k η factors between the Pythia tunes and the extrapolation procedure used for the k η factors in two highest η bins. No systematic uncertainty is assigned to account for inaccuracies of the Lund string model in describing the ratio of the neutral energy flow to the charged energy flow. The uncertainties associated with the ECAL energy calibration and photon reconstruction efficiency also affect the accuracy of the total energy flow. To account for a possible difference in the photon reconstruction efficiency between simulation and data, a global 3.7 % systematic uncertainty is assigned following the analysis presented in Ref. [36]. The energy calibration of the ECAL has been studied by measuring the invariant masses of diphoton resonances (π 0γγ and ηγγ) and has been assigned a global systematic uncertainty of 1.5 %. Both uncertainties are scaled with a factor F neut,η /F total,η and are listed in Table 1.

Other potential sources of systematic uncertainties such as momentum- and η-smearing, effect of the beam crossing angle, neglecting the masses of charged particles, pollution from elastic scattering and beam-gas interactions have been studied as well. Their impacts on the accuracy of the measurements are found to be negligible.

The total systematic uncertainties on the corrected charged and total energy flow are listed in Tables 2 and 3 for all event classes and η bins. It should be noted that the uncertainties are strongly correlated between the bins.

Table 2 Charged energy flow for all event classes and η bins with the corresponding systematic uncertainties. The statistical uncertainties are insignificant and not listed. All values are in GeV per unit pseudorapidity interval
Table 3 Total energy flow for all event classes and η bins with the corresponding systematic uncertainties. The statistical uncertainties are insignificant and not listed. All values are in GeV per unit pseudorapidity interval

5 Results

The fully-corrected measurements for the charged energy flow are shown in Fig. 1 for each event class together with the generator level predictions given by the Pythia tunes and the corresponding systematic uncertainties. By comparing experimental results obtained for different event classes one can clearly see that the amount of energy flow strongly correlates with the momentum transfer in an underlying pp inelastic interaction. The charged energy flow rises more steeply with pseudorapidity in data than predicted by the majority of the Pythia tunes. As a consequence, the discrepancy between the measurements and generator level predictions increases towards large η rising to 20 % in the last η bin. At lower η the data are reasonably well described by the Pythia tunes. This is the case for all event classes except the diffractive enriched one. For the latter, the measurements are well described by the Pythia 8.135 generator with default parameters. However, this model overestimates the charged energy flow in the case of hard scattering events over the entire pseudorapidity range of the measurements.

Fig. 1
figure 1

Charged energy flow as a function of η for all event classes as indicated in the figures. The corrected measurements are given by points with error bars, while the predictions by the Pythia tunes are shown as histograms. The error bars represent the systematic uncertainties, which are highly correlated between the bins. The statistical uncertainties are negligible. The ratios of MC predictions to data are shown in addition

Figure 2 illustrates the charged energy flow along with the predictions given by the cosmic-ray interaction models. It is interesting to note that the measurements performed with inclusive minimum-bias and non-diffractive enriched events are well described by the Epos 1.99 and Sibyll 2.1 models, while the Qgsjet01 and QgsjetII-03 generators overestimate the charged energy flow for these event classes. The latter also occurs at large η in the case of hard scattering interactions for all cosmic-ray interaction models except the QgsjetII-03. The diffractive enriched charged energy flow is underestimated by all cosmic-ray interaction models.

Fig. 2
figure 2

Charged energy flow as a function of η for all event classes as indicated in the figures. The corrected measurements are given by points with error bars, while the predictions by the cosmic-ray interaction models are shown as histograms. The error bars represent the systematic uncertainties, which are highly correlated between the bins. The statistical uncertainties are negligible. The ratios of MC predictions to data are shown in addition

The total energy flow is shown for each event class in Fig. 3 along with the generator level predictions given by the Pythia tunes and the corresponding systematic uncertainties. It can be clearly seen that all Pythia 6 tunes underestimate the amount of energy flow at large pseudorapidity for all event classes. The Pythia 8.135 generator gives the best description of the measurements performed with inclusive minimum-bias, diffractive and non-diffractive enriched events among the Pythia tunes, except for the pseudorapidity range 1.9<η<2.5. None of these models provide an accurate description of the energy flow measured with hard scattering events. The predictions given by the LHCb and Perugia NOCR tunes for non-diffractive enriched and hard scattering events are rather similar, while the Perugia 0 tune significantly underestimates the energy flow for all event classes. For diffractive enriched events, the inconsistency between the data and the prediction given by the LHCb tune is found to be rather large throughout the entire pseudorapidity range 1.9<η<4.9, while the Pythia 8.135 generator with default parameters gives a good description of the corresponding energy flow at large η.

Fig. 3
figure 3

Total energy flow as a function of η for all event classes as indicated in the figures. The corrected measurements are given by points with error bars, while the predictions by the Pythia tunes are shown as histograms. The data obtained with extrapolated k η factors are shown in grey. The error bars represent the systematic uncertainties, which are highly correlated between the bins. The statistical uncertainties are negligible. The ratios of MC predictions to data are shown in addition

Figure 4 illustrates the total energy flow together with the predictions given by the cosmic-ray interaction models. It is observed that the Sibyll 2.1 generator gives the best description of the energy flow measured with inclusive minimum-bias and non-diffractive enriched events at large η. The predictions given by the Epos 1.99 generator for these event classes also describe the measurements reasonably well. In the case of hard scattering interactions, the best description of the data at large η is given by the QgsjetII-03 generator. The total energy flow measured with diffractive enriched events is underestimated at large η by all cosmic-ray interaction models used in this study. The measurements of the charged and total energy flow are summarised in Tables 2 and 3 for all event classes and η bins.

Fig. 4
figure 4

Total energy flow as a function of η for all event classes as indicated in the figures. The corrected measurements are given by points with error bars, while the predictions by the cosmic-ray interaction models are shown as histograms. The data obtained with extrapolated k η factors are shown in grey. The error bars represent the systematic uncertainties, which are highly correlated between the bins. The statistical uncertainties are negligible. The ratios of MC predictions to data are shown in addition

The results obtained in this study cannot be directly compared with the measurements performed by the CMS collaboration [5], since different event selection criteria are applied in the analyses. Nevertheless, both measurements demonstrate that the energy flow is underestimated by Pythia 6 tunes at large pseudorapidity, while the results of the ATLAS collaboration indicate that the amount of transverse energy is also underestimated by various Pythia tunes at large η [6].

6 Conclusions

The energy flow is measured in the pseudorapidity range 1.9<η<4.9 with data collected by the LHCb experiment in pp collisions at \(\sqrt{s}=7\ \mbox{TeV}\) for inclusive minimum-bias interactions, hard scattering processes and events with enhanced or suppressed diffractive contribution. The primary measurement is the energy flow carried by charged particles. For the measurement of the total energy flow, a data-constrained MC estimate of the neutral component is used. The energy flow is found to increase with the momentum transfer in an underlying pp inelastic interaction. The evolution of the energy flow as a function of pseudorapidity is reasonably well reproduced by the MC models. Nevertheless, the majority of the Pythia tunes underestimate the measurements at large pseudorapidity, while most of the cosmic-ray interaction models overestimate them, except for diffractive enriched interactions. For inclusive and non-diffractive enriched events, the best description of the data at large η is given by the Sibyll 2.1 and Pythia 8.135 generators. The latter also provides a good description of the energy flow measured with diffractive enriched events, especially at large η. The comparison shows that the absence of hard diffractive processes moderates the amount of the forward energy flow meaning that their inclusion is vital for a more precise description of partonic interactions. It also demonstrates that higher-order QCD effects as contained in the Pomeron phenomenology play an important role in the forward region. None of the event generators used in this analysis are able to describe the energy flow measurements for all event classes that have been studied.