abstract: To truly eliminate Cartesian ghosts from the science of consciousness, we must describe consciousness as an aspect of the physical. Integrated Information Theory states that consciousness arises from intrinsic information generated by dynamical systems; however existing formulations of this theory are not applicable to standard models of fundamental physical entities. Modern physics has shown that fields are fundamental entities, and in particular that the electromagnetic field is fundamental. Here I hypothesize that consciousness arises from information intrinsic to fundamental fields. This hypothesis unites fundamental physics with what we know empirically about the neuroscience underlying consciousness, and it bypasses the need to consider quantum effects.
author: Adam B. Barrett[^1]

*Sackler Centre for Consciousness Science* and *Department of Informatics*
University of Sussex, Brighton BN1 9QJ, UK
date: \[Published Feb. 4, 2014 in the *Consciousness Research* specialty section of *Frontiers in Psychology*, article no. 5(63).\]
title: An Integration of Integrated Information Theory with Fundamental Physics

# Introduction

The key question in consciousness science is: "Given that consciousness (i.e., subjective experience) exists, what are the physical and biological mechanisms underlying the generation of consciousness?". From a basic property of our phenomenology, namely that conscious experiences are integrated representations of large amounts of information, Integrated Information Theory (IIT) hypothesizes that, at the most fundamental level of description, consciousness is integrated information, defined as information generated by a whole system, over and above its parts (Tononi, 2008). Further, given the private, non-externally observable nature of consciousness, IIT considers consciousness to be an intrinsic property of matter, as fundamental as mass, charge or energy. Thus, more precisely, IIT posits that consciousness is intrinsic integrated information, where by intrinsic information it is meant that which is independent of the frame of reference imposed by outside observers of the system. The quantity of consciousness generated by a system is the amount of intrinsic integrated information generated (Balduzzi and Tononi, 2008), whilst the qualities of that consciousness arise from the precise nature of informational relationships between the parts of the system (Balduzzi and Tononi, 2009).

IIT has garnered substantial attention amongst consciousness researchers. However, it has been criticized for its proposed measures of integrated information not successfully being based on an intrinsic perspective (Gamez, 2011; Beaton and Aleksander, 2012; Searle, 2013). The proposed "$\Phi$" measures are applicable only to networks of discrete nodes, and thus for a complex system depend on the observer choosing a particular graining. More broadly, information can only be intrinsic to fundamental physical entities, and descriptions of information in systems modeled at a non-fundamental level necessarily rely on an extrinsic observer's choice of level (Floridi, 2009, 2010; Gamez, 2011). Here I propose a potential solution to this problem, what might be called the field integrated information hypothesis (FIIH). Modern theoretical physics describes the universe as being fundamentally composed of continuous fields. Electrical signals are the predominant substrate of information processing in brains, and the electromagnetic field that these produce is considered fundamental in physics, i.e., it is not a composite of other fields. Thus, I hypothesize that consciousness arises from information intrinsic to fundamental fields, and propose that, to move IIT forward, what is needed is a measure of intrinsic information applicable to the configuration of a continuous field.

The remainder of this article is laid out as follows. First I discuss the concept of fundamental fields in physics, and how if one takes the view that consciousness is an intrinsic property of matter, then it must be a property arising from configurations of fields. In the following section, I discuss the hypothesis that consciousness arises from integrated information intrinsic to fundamental fields, the shortcomings of existing approaches to integrated information, and the possibility of constructing a measure that can successfully measure this quantity for field configurations. I then explain how IIT and the FIIH imply a limited form of panpsychism, and why this should not be considered a problem, before contrasting the FIIH with previously proposed field theories of consciousness, such as that of Pockett (2000). Finally, the summary includes some justification for this theoretical approach to consciousness.

# Fundamental fields and consciousness

|  | **Mass (GeV/$c^2$)** | **Electric charge** | **Strong charge** | **Weak charge** |
|:---|:---|:---|:---|:---|
| **LEPTONIC MATTER** |  |  |  |  |
| electron neutrino ($\nu_e$) | $<1.3\times10^{-10}$ | 0 | No | Yes |
| electron (e) | 0.0005 | -1 | No | Yes |
| muon neutrino ($\nu_\mu$) | $<1.3\times10^{-10}$ | 0 | No | Yes |
| muon ($\mu$) | 0.106 | -1 | No | Yes |
| tau neutrino ($\nu_\tau$) | $<1.4\times10^{-10}$ | 0 | No | Yes |
| tau ($\tau$) | 1.78 | -1 | No | Yes |
| **QUARK MATTER** |  |  |  |  |
| up (u) | 0.002 | 2/3 | Yes | Yes |
| down (d) | 0.005 | -1/3 | Yes | Yes |
| charm (c) | 1.3 | 2/3 | Yes | Yes |
| strange (s) | 0.1 | -1/3 | Yes | Yes |
| top (t) | 173 | 2/3 | Yes | Yes |
| bottom (b) | 4.2 | -1/3 | Yes | Yes |
| **BOSONS** |  |  |  |  |
| **Electromagnetic force:** |  |  |  |  |
| photon ($\gamma$) | 0 | 0 | No | No |
| **Strong force:** |  |  |  |  |
| gluon (g) | 0 | 0 | Yes | No |
| **Weak force:** |  |  |  |  |
| $W^-$ | 80 | -1 | No | No |
| $W^+$ | 80 | 1 | No | No |
| Z | 91 | 0 | No | No |
| **Gravity:** |  |  |  |  |
| graviton$^*$ | 0 | 0 | No | No |
| **Higgs mechanism:** |  |  |  |  |
| Higgs (H) | 126 | 0 | No | Yes |
|  |  |  |  |  |

Table of the fields/particles that are considered fundamental. Familiar matter arises from leptons and quarks, while the forces of nature arise from interactions of matter with "carrier" bosons. Mass is given in giga electron volts per speed of light squared (Gev/$c^2\approx 2\times10^{-27}$kg). Electric charge is in standard units relative to minus the charge of the electron, i.e., one unit equals $1.6\times10^{-19}$ Coulombs. A description of the group theoretic strong and weak charges is beyond the scope of this article, but the table shows which fields have strong and weak charges. \*The gravity field is considered fundamental and is well-studied, but the gravity particle (graviton) has not to date explicitly been observed; at quantum (i.e., very microscopic) spatial scales, a consistent set of field equations for gravity have yet to be constructed.

Contemporary physics postulates that "fields" are the fundamental physical ingredients of the universe, with the more familiar quantum particles arising as the result of microscopic fluctuations propagating across fields, see e.g., Oerter (2006) for a lay person's account, or Coughlan et al. (2006) for an introduction for scientists. In theoretical terms, a field is an abstract mathematical entity, which assigns a mathematical object (e.g., scalar, vector) to every point in space and time. (Formally a field is a mapping $F$ from the set $S$ of points in spacetime to a scalar or vector field $X$, $F: S \to X$.) So, in the simplest case, the field has a number associated with it at all points in space. At a very microscopic scale, ripples, i.e., small perturbations, move through this field of numbers, and obey the laws of quantum mechanics. These ripples correspond to the particles that we are composed of, and there is precisely one fundamental field for each species of fundamental particle. At the more macroscopic level, gradients in field values across space give rise to forces acting on particles. The Earth's gravitational field, or the electromagnetic field around a statically charged object, are examples of this, and the classical (as opposed to quantum) description is a good approximation at this spatial scale. However, both levels of description can be considered equally fundamental if the field is fundamental, i.e., not some combination of other simpler fields. Note that the electromagnetic and gravitational fields are both examples of fundamental fields, with the corresponding fundamental particles being the photon and the graviton. Particles are divided up into matter particles and force-carrying particles, but all types of particle have associated fields; all the forces of nature can be described by field theories which model interactions, i.e., exchanges of energy, between fields. See Table 1 for a list of fields/particles that are considered fundamental according to this so-called "Standard Model" of particle physics.

To be consistent with modern theoretical physics, a theory of consciousness that considers consciousness to be a fundamental attribute of matter must describe how consciousness manifests itself in the behavior of either fundamental fields or quantum particles. Since we know that the brain generates electric fields with a rich spatiotemporal structure, and that, for the main part, information processing in the brain is carried out by electrical signaling between neurons operating mostly in the classical (as opposed to quantum) regime (Koch and Hepp, 2006), empirical evidence favors the former. Thus, on the view that consciousness is a fundamental attribute of matter, it must be the structure and/or dynamics of the electromagnetic field (which is an example of a fundamental field) that is fundamentally the generator of brain-based consciousness.

Once one ascribes electromagnetic fields with the potential to generate consciousness, it is natural to ask whether other fields might also have the potential to generate consciousness. According to modern physics, there was a symmetry between all fields at the origin of the universe, although these symmetries were broken as the universe began to cool (Georgi and Glashow, 1974; see Hawking, 2011 for a lay-person's account). It could be argued by Occam's razor that it makes more sense to posit that potential for consciousness existed at the outset, and hence potential for consciousness is a property of all fields, than that it emerged only during symmetry breaking. However, in practice, it is unlikely that any complex consciousness could exist in any field other than the electromagnetic field, for reasons to do with the physics and chemistry of the electromagnetic field compared with other fields. Considering the four forces: strong, weak, electromagnetic and gravitational, the strong and weak forces don't propagate over distances much larger than the width of the nucleus of an atom, and gravity alone cannot generate complex structures by virtue of being solely attractive; in contrast, the electromagnetic field can propagate over macroscopic scales, is both repulsive and attractive, and is fundamentally what enables non-trivial chemistry and biology. Considering fields associated with matter, these in general do not have any undulations at spatial scales larger than the quantum scale; the non-trivial structures in these fields are essentially just the ripples associated with the familiar quantum matter particles, i.e., electrons and quarks, and various "exotic" particles detectable in particle physics experiments (see Table 1). Finally, the recently discovered Higgs field has essentially a uniform structure; quantum interactions exist between the Higgs field and many of the other fields, and this is fundamentally the origin of mass in the universe (see e.g., Coughlan et al., 2006; Oerter, 2006). Thus, the physics of the electromagnetic field uniquely lends itself to the generation of complex structures.

# The Field Integrated Information Hypothesis

Given the above, I propose that the principal conceptual postulates of IIT should be restated as follows. Consciousness arises from information intrinsic to the configuration of a fundamental field. The amount of consciousness generated by a patch of field is the amount of integrated information intrinsic to it. When a patch of field generates a large quantity of intrinsic integrated information, mathematically there is a high-dimensional informational structure associated with it (Tononi, 2008; Balduzzi and Tononi, 2009). The geometrical and topological details of this structure determine the contents of consciousness. The task now is to correctly mathematically characterize intrinsic integrated information, and construct equations to measure it.

A true measure of intrinsic integrated information must be frame invariant, just like any fundamental quantity in physics. That is, it must be independent of the point of view of the observer: independent of the units used to quantify distance or time, independent of which direction is up, and independent of the position of the origin of the coordinate system; and also independent of the scale used for quantifying charge, or field strength.

The "$\Phi$" measures put forth by existing formulations of IIT (Balduzzi and Tononi, 2008; Barrett and Seth, 2011) are not applicable to fields because they require a system with discrete elements, and fields are continuous in space. One could ask, however, whether a perspective on a system in terms of discrete elements could actually be equivalent to an intrinsic field-based perspective, thus obviating the need for a field-based measure. To see explicitly that this is not the case, let us revisit the photodiode, which, according to the existing theory (Tononi, 2008), has 1 bit of intrinsic information by virtue of having two states, on or off. There is a wire inside the photodiode, and the electrons inside the wire are all individually fluctuating amongst many different states. The electromagnetic field generated by the diode, and the circuit to which it is connected has two stable configurations for as long as the circuit is connected. But other more general configurations for an electromagnetic field are ruled out by each of these states. Considering the system at this level of description yields a distinct perspective, and would lead one to deduce that the amount of information generated by the system's states is some quantity other than 1 bit. Thus the field-based perspective is not equivalent to the observer-dependent discrete perspective.

The idea here is that a formula should be obtained that could in theory be applied universally to explore the intrinsic information in any patch of spacetime, without requiring an observer to do any modeling, i.e., one would just measure field values in as fine a graining as possible to get the best possible approximations to the intrinsic informational structure. Only a formula in continuous space and time would allow this. If a discrete formula were to be applied, there would always be the possibility of encountering an informational structure on a finer scale than that of the formula. (Unless the graining required by the formula were the Planck scale, i.e., the scale of the hypothesized superstring, on which continuous models of physics break down; however there do not exist complex structures at that scale.) In practice however, observations of systems are necessarily discrete, so discrete approximations to a continuous formula could be useful for empirical application. See Balduzzi (2012) for some recent work on the information-theoretic structure of distributed measurements.

We don't yet know how to properly calculate intrinsic information, so must remain agnostic on the precise amount of intrinsic integrated information generated by photodiodes, or of anything. However, the failure of existing approaches does not rule out the construction in the future of a successful formula. While it is beyond the scope of this present paper to make a serious attempt at solving this problem, I speculate that a formula in terms of thermodynamic entropy as opposed to Shannon entropy might be more likely to succeed, as the former is inherently an intrinsic property, whereas the latter was constructed for the purpose of describing an external observer's knowledge of a system (Floridi, 2009, 2010; Gamez, 2011; Beaton and Aleksander, 2012).

# Integrated Information Theory and panpsychism

Searle (2013) criticizes IIT for its stance that integrated information always produces consciousness, stating that this ludicrously ascribes consciousness to all kinds of everyday objects and would mean that consciousness is "spread thinly like a jam across the universe". Koch and Tononi (2013) counter that only "local maxima" of integrated information exist (over spatial and temporal scales): "my consciousness, your consciousness, but nothing in between". If local maxima of intrinsic integrated information in field configurations always generate consciousness, then there must be minute amounts, say "germs", of consciousness all over the universe, even though there would be no superordinate consciousness amongst groups of people. Thus, IIT and the FIIH do imply a form of panpsychism. However, the phenomenology assigned to an isolated electron in a vacuum, or even a tree, which has no complex electromagnetic field, would be very minimal. Since the only consciousness we can be certain of is our own, the positing by integrated information theories of germs of consciousness everywhere is no reason to dismiss them. A theory should stand or fall on whether or not it can elegantly and empirically describe human consciousness.

For those uncomfortable with subscribing to a panpsychist theory, a possible way round the problem is to assign an attribute "potential consciousness" to matter at the most fundamental level. Then, the quantity of potential consciousness is simply the quantity of integrated intrinsic information. But only when there is a large amount of intrinsic integrated information with a sufficiently rich structure to be worthy of being compared to a typical healthy adult human waking conscious moment, should we say that the integrated information has "actual consciousness" associated with it. A line could thus be drawn somewhere between the potential consciousness of an isolated electron in a vacuum and the actual consciousness generated by my brain as I write this article. The problem with such a distinction however is that potential consciousness would still be assigned phenomenal content, so it is perhaps more elegant to just use a single term "consciousness" for the whole spectrum of integrated information. On the other hand, since consciousness is defined by some as any mental content, but by others as only self-reflective mental content, there is no single terminology that appeals to everybody. The key point, irrespective of the precise definition of consciousness, is that on the theory discussed here, intrinsic integrated information is what underlies subjective experience at the most fundamental level of description. Alternatively, one could further imagine different lines being drawn for different purposes. For example, a threshold of conscious awareness above which surgery cannot be performed; or thresholds at which various people are comfortable eating animals.

# Relation to previous electromagnetic field theories of consciousness

There have been several other theories of consciousness put forward that identify consciousness with various types or configurations of fields, see Pockett (2013) for a review. Notably, Pockett's electromagnetic field theory (EMT) of consciousness (Pockett, 2000, 2011, 2012) posits that "conscious perceptions (and sensations, inasmuch as they can be said to have independent existence) are identical with certain spatiotemporal electromagnetic patterns generated by the normal functioning of waking mammalian brains" (Pockett, 2013). In the most recent formulation of this theory, the key feature of field patterns underlying consciousness is the presence of a neutral region in the middle of a radial pattern. This hypothesis was motivated by the observation that such field patterns appear during recurrent cortical activity, (with the neutral region in layer 4), and the empirical association of consciousness with recurrent processing (Pockett, 2012).

A problem common to previous field theories of consciousness (Libet, 1994; Pockett, 2000, 2013; McFadden, 2002) is that they claim that cutting outgoing neural connections from a slab of cortex that generates a conscious experience will not affect the ability to report that conscious experience. EMT argues that the electromagnetic field within such an isolated hypothetical slab would still propagate through space and enable communication between the conscious field generated by the slab and the spatially contiguous larger conscious mental field. This is not however compatible with the laws of physics. Any cutting of synapses to or from regions of cortex that are generating consciousness will alter the field, and will therefore alter the conscious experience. There is no electromagnetic field residing in the brain other than that generated specifically by all of the neural and chemical activity. And it does not make sense to talk of the brain's electromagnetic field and its firing neurons and synapses as being able to exist independently of each other. On the theory put forward here, neurons can be considered the scaffolding that enable very complex electromagnetic field configurations to be sustained. As far as describing the mechanisms of perception and cognition that generate the specific contents of consciousness in any given scenario, the current paradigm of associating it with neural activity is of course the only valid and useful level of description. However, in terms of explaining more fundamentally how matter gives rise to consciousness, a description in terms of fields would be much more elegant than a description in terms of the complex entities that are neurons.

Another shortcoming of previous field theories of consciousness is that none of them relate physical properties of proposed correlates of consciousness to properties of phenomenology, i.e., they do not posit "explanatory correlates of consciousness" (Seth, 2009). The FIIH raises for the first time the possibility of constructing a field theory of consciousness that can account for a fundamental aspect of phenomenology, namely that conscious experiences are integrated representations of large amounts of information.

# Discussion

In this paper I have hypothesized that, at the most fundamental level of description, human consciousness arises from information intrinsic to the complex electromagnetic fields generated by the brain. This "FIIH" builds on the axioms of IIT, namely that consciousness is integrated information, and that consciousness is an intrinsic and fundamental property of matter analogous to mass or charge. However, it also implies that a new mathematical formalism is required to properly quantify intrinsic integrated information, since electromagnetic fields are continuous in space, and existing "$\Phi$"-type measures of integrated information are applicable only to discrete systems (which require an observer dependent perspective). The idea that consciousness can be identified with certain spatiotemporal electromagnetic patterns has been previously put forward in other electromagnetic field theories of consciousness. But by suggesting that integrated information is the key factor, the theory here connects, for the first time, such electromagnetic field theories of consciousness to basic aspects of phenomenology.

The hypothesis is admittedly rather speculative, and any proposed mathematical formula for conscious level in terms of information intrinsic to an electromagnetic field will be difficult to test directly, simply because we do not have the technological tools or the computational resources to record in full detail the three-dimensional electromagnetic field structure generated by the brain. Rather, this can only be sampled at a spatial scale that is sparse compared to the finest scale of its undulations. However, there is a strong case to be made that the theoretical development of the ideas presented here has substantial value. Theories in physics have been vigorously pursued for their logic and beauty, in the absence of imminent direct experimental tests. For example, there is a vast amount of work being conducted on string theory; there, rather than experimental verification, the goal is an elegant explanation of our existing empirical knowledge of particle physics and gravity. If there already existed several analogous theories of consciousness, then one could argue that it would not be useful to add to the speculation. However, there is as yet no compellingly believable set of equations for describing, fundamentally, how consciousness is generated. IIT has potential in this direction, but a major step forward for the theory would be a truly plausible formula for intrinsic information applicable to fundamental physical entities. The FIIH provides a conceptual starting point for achieving this. All this is not to say that such a theory will aid understanding of all aspects of consciousness; indeed the multi-faceted nature of consciousness requires descriptions at many different levels. Non-reductionist frameworks are required to understand the complexity of the biological machinery that enables the brain to do any kind of information processing, conscious or unconscious, and to understand the differences between conscious and unconscious cognitive processes neural dynamics and behavior must necessarily be modeled at multiple levels of description.

Finally, any theory can potentially indirectly make predictions. Indeed IIT has already inspired heuristic measures of information integration/complexity that have been successfully applied to recorded electrophysiological data and are able to distinguish the waking state from diverse unconscious states, i.e., sleep and anaesthesia under various anaesthetics (Massimini et al., 2005; Ferrarelli et al., 2010; Casali et al., 2013). The results are in broad agreement with the predictions of IIT and provide encouragement for further theoretical work on the relationship between information integration and consciousness. Theories built from the FIIH could make new and distinct predictions about the types of structural and/or functional neuronal architectures that are capable of generating consciousness; and new theory can only further inform the quest for ever more reliable measures of consciousness that can be applied to observable brain variables.

# Acknowledgements

I thank Emily Lydgate and Anil Seth for invaluable discussions during the writing of this paper, and Daniel Bor and David Gamez for very useful comments on draft manuscripts. ABB is funded by EPSRC grant EP/L005131/1.

[^1]: email@example.com


author: [^1]

The rejection of the contamination, or background, from low-energy strong interactions at hadron collider experiments is a topic that has received significant attention in the field of particle physics. This article builds on a particle-level view of collision events, in line with recently-proposed subtraction methods. While conventional techniques in the field usually concentrate on probability distributions, our study is, to our knowledge, the first attempt at estimating the frequency distribution of background particles across the kinematic space inside individual collision events. In fact, while the probability distribution can generally be estimated given a model of low-energy strong interactions, the corresponding frequency distribution inside a single event typically deviates from the average and cannot be predicted a priori. We present preliminary results in this direction, and establish a connection between our technique and the particle weighting methods that have been the subject of recent investigation at the Large Hadron Collider.

**Keywords:** 29.85.Fj; High Energy Physics; Particle Physics; Large Hadron Collider; LHC; soft QCD; pile-up; mixture models; Gibbs sampler; Markov Chain Monte Carlo; Expectation Maximisation.

# Nomenclature and general remarks

-   Collisions: proton-proton collisions at the Large Hadron Collider.

-   Events: triggered proton-proton collisions.

-   Bunch crossings: intersections between colliding proton beam bunches.

-   Physics processes: either the high-energy parton scattering of interest or low-energy strong interactions.

-   Missing transverse energy: event-level energy imbalance measured on a plane perpendicular to the direction of the colliding particle beams.

-   Particle transverse momentum, $p_T$: absolute value of the component of the particle momentum vector on a plane perpendicular to the direction of the colliding beams.

-   Particle pseudorapidity, $\eta$: kinematic quantity expressed in terms of the particle polar angle in the laboratory frame, $\theta$, by $\eta=-\mbox{log}\left[\mbox{tan}(\theta/2)\right]$.

-   Whenever neutral particles are referred to in the text, neutrinos are not considered.

# Introduction

The subtraction of the contamination, or background, from soft, i.e. low-energy, physics processes described by Quantum Chromodynamics (QCD) that take place in proton-proton collisions is a critical task at the Large Hadron Collider (LHC). The impact of the correction is going to become even more significant in the upcoming scenarios whereby the hard, i.e. high-energy, parton scattering of interest will be superimposed with a higher number of low-energy interactions from collisions between other protons, the so-called pile-up events. This is an important aspect at the LHC, and one that is going to have an even more significant impact at the High-Luminosity LHC (HL-LHC), i.e. at the accelerator that is going to be built following the LHC upgrade project. In fact, the contribution of pile-up particles to the events of interest often makes the study of rare processes particularly challenging.

Subtraction techniques are well established, and typically combine tracking information for charged particles with estimates of the energy flow associated with neutral particles that originate from low-energy QCD interactions . In particular, pile-up subtraction is a key component in the data processing pipelines responsible for the reconstruction and calibration of jets, i.e. of collections of particles interpreted as originating from the same scattered parton.

In this context, an important role has been played by correction procedures based on jet area , which provides a measure of the susceptibility of reconstructed jets to the soft QCD energy flow. Such methods work by subtracting from the total momentum of the high-energy jets a quantity proportional to an event-level estimate of the background momentum density as well as to the area of the jet of interest. Therefore, this takes into account event-to-event background variability, and, since the correction can be calculated in a kinematics-depedent way, also the presence of different levels of pile-up activity in different kinematic regions inside events. However, due to the quantum nature of the underlying physics, the density of pile-up particles can be different even inside jets with very similar kinematics in the same collision event. While techniques based on jet area cannot correct for this, more recent methods exploiting information encoded in the substructure of jets can effectively take this into account.

In this article, we explore a different perspective in order to estimate the frequency distribution of soft QCD particles inside events. We build on a view of collision events as collections of particles whereby soft QCD particles are rejected upstream of jet reconstruction, in line with the particle-level pile-up subtraction algorithms that are currently being evaluated at the LHC .

Due to the quantum nature of the underlying physics processes and the limited number of final-state particles inside individual collisions, the particle multiplicity across the kinematic space inside events will generally vary across collisions even when the physics processes involved are exactly the same. More precisely, the soft QCD particle-level frequency distribution will normally deviate from the corresponding probability distribution, and will be different in different events. What is discussed in this article is a data-driven method of estimating the soft QCD particle multiplicity across the kinematic space inside each event, using the following:

-   The kinematic probability distributions of soft QCD particles and of particles originating from the signal hard scattering, e.g. obtained using simulated data;

-   The average fraction of soft QCD particles in the events;

-   The observed particle multiplicity, i.e. the observed number of particles in different kinematic regions in the event.

Our approach relies on particles from high-energy scattering processes having pseudorapidity distribution more peaked at $\eta=0$, as well as higher values of $p_T$, on average, than those originating from low-energy strong interactions. This is essentially due to the higher transverse momentum transferred between the colliding protons, and results in different kinematic distributions of the final-state particles on the ($\eta, p_T$) plane, as illustrated in Fig. with reference to the control sample. Although different signal processes will generally be associated with different kinematic signatures, the dissimilarity between background soft QCD and hard scattering signal particles in terms of their $\eta$ and $p_T$ distributions typically outweighs the variability associated with the choice of signal process.

Moreover, as a filtering stage upstream of physics analysis, reconstructed events at the experiments are usually subdivided, or "skimmed", into multiple data streams enriched in different signal processes. For the purpose of this discussion, the signal model can be thought of as describing the particle-level kinematics corresponding to the high-energy processes that the events analysed are enriched in.

The kinematic variables used in this study are those that the relevant signal and background probability distributions can be written as functions of, namely particle pseudorapidity, $\eta$, and transverse momentum, $p_T$. To our knowledge, this is the first method of estimating how the frequency distribution of soft QCD particles inside individual events deviates from the expectation due to the non-deterministic nature of the underlying processes.

This article reports preliminary results on simulated collision events at the LHC, showing that the algorithm produces reasonable estimates of the number of soft QCD particles in different $(\eta, p_T)$ regions inside events regardless of the presence in those regions of particles from the hard scattering. Given that assessing the performance of this method requires knowledge of the true frequency distribution of soft QCD particles in each event, which is not available at the experiments, the validation was performed on simulated data, using an event generator commonly employed in the field . Specifically, background and signal particles were generated using soft QCD processes and $gg\rightarrow t\bar{t}$, respectively.

Our interest in the estimation of the multiplicity of soft QCD particles across the kinematic space inside individual collision events relates to the development of further-improved methods of rejecting pile-up in high-luminosity regimes. Since our approach is based on a different principle and works in a different way as compared to established techniques, we expect its combined use with existing methods to result in enhanced pile-up subtraction in the upcoming higher-luminosity regimes at the LHC. We speculate that this can also lead to improved missing transverse energy resolution and to higher-quality estimates of particle isolation as the pile-up rates increase.

It should be noted that, in addition to pile-up, another source of soft QCD particles at the LHC is the so called Underlying Event (UE), which consists of particles from low-energy parton interactions taking place in the same proton-proton collision that contains the particles produced by the hard parton scattering. Pile-up and UE particles originate from similar processes: for this reason, with regard to estimating the frequency distribution of soft QCD particles inside events, it is expected that UE particles will contribute to the background category, i.e. that they will count towards the number of soft QCD particles together with those that originate from pile-up. In any case, although the distinction between pile-up and UE particles is conceptually important, pile-up is by far the primary source of soft QCD contamination at the ATLAS and CMS experiments at the LHC.

The algorithm that we describe in this article is a simplified deterministic variant of the Markov Chain Monte Carlo technique used in , where we discussed the idea of filtering collision events particle by particle upstream of jet reconstruction. Both our previous contributions and the present article relate to the development of new subtraction methods, with a view to improving further on the rejection of contamination from low-energy strong interactions in high-luminosity hadron collider environments. In particular, it is our opinion that the simplicity and parallelisation potential of this technique make it a promising candidate for inclusion in particle-by-particle event filtering procedures at the reconstruction level at future high-luminosity hadron collider experiments.

# The algorithm

By construction, the probability density function (PDF) describing the kinematics of particles originating from a given process, e.g. with reference to soft QCD interactions, can be estimated as the limit of the corresponding frequency distribution, averaged over a large enough number of events. On the other hand, the corresponding frequency distribution inside a single event normally deviates from the PDF due to the limited number of particles. In fact, even when the processes involved are exactly the same, different collisions contain independent, and therefore different, realisations of the underlying quantum processes. For this reason, the shapes of the corresponding particle-level frequency distributions, e.g. that of soft QCD particles, generally vary across collisions.

Let $f_0$ and $f_1$ denote the kinematic PDFs of background and signal particles, respectively, normalised such that $\int\int f_i(\eta, p_T) d\eta dp_T = 1,~i=0,1$. For the purpose of this study, we describe collision events as statistical populations of particles originating from soft QCD interactions and from the hard scattering, using a mixture model of the form $\alpha_0 f_0(\eta, p_T) + \alpha_1 f_1(\eta, p_T)$, where $\alpha_0$ is the fraction of soft QCD particles in the events, and $\alpha_1 = 1-\alpha_0$.

In this context, the probability for a given particle to originate from soft QCD interactions can be expressed using the following quantity:

$$w_0(\eta, p_T) \equiv \displaystyle \frac{\alpha_0 f_0(\eta, p_T)}{\alpha_0 f_0(\eta, p_T) + \alpha_1 f_1(\eta, p_T)}.
\label{eq:w0}$$

Inside each collision event, although the actual numbers of background and signal particles in the different $(\eta, p_T)$ bins are not known, it is possible to estimate the corresponding expected numbers, $\nu_b(\eta, p_T)$ and $\nu_s(\eta, p_T)$, given a background and a signal model, respectively.

For the purpose of this discussion, we express $\nu_b$ in terms of $\nu_b(\eta, p_T) = N \alpha_0 f_0(\eta, p_T)\Delta\eta \Delta p_T$, where $N$ is the total number of particles in the event, and $\Delta\eta$ and $\Delta p_T$ are the bin widths along the $\eta$ and $p_T$ axes, respectively. The corresponding expected number of signal particles in the bin, $\nu_s(\eta, p_T)$, can be calculated in a similar way using $f_1$.

If one denotes the unknown true numbers of signal and soft QCD particles as functions of particle $\eta$ and $p_T$ in each event by $n_s^*(\eta, p_T)$ and $n_b^*(\eta, p_T)$, respectively, then $n(\eta, p_T) = n_s^*(\eta, p_T) + n_b^*(\eta, p_T)$, where $n(\eta, p_T)$ is the corresponding number of particles in the data. When one considers LHC events with a number of proton-proton interactions in line with what is expected in the upcoming higher-luminosity regimes, the final-state particle multiplicities associated with soft QCD interactions and with the signal hard scattering are such that the expected number of signal particles in each bin, $\nu_s(\eta, p_T)$, is on average much lower than the corresponding number of soft QCD particles, i.e. $\left<\nu_s(\eta, p_T)\right> \ll \left<n_b^*(\eta, p_T)\right>$, where the average is taken over the $(\eta, p_T)$ space.

One therefore expects that the statistical fluctuations on the observed number of particles in each $(\eta, p_T)$ bin will also be dominated by those on the number of soft QCD particles, i.e. $\left<\sigma_{n_s}(\eta, p_T)\right> \ll \left<\sigma_{n_b}(\eta, p_T)\right>$.

It should be noted that, if the variability of the number of signal particles can be neglected, the quantity

$$\hat{n}_b(\eta, p_T) = w_0(\eta, p_T) n(\eta, p_T)
\label{eq:nbhat}$$

is expected to provide a more reliable estimate of the unknown number of soft QCD particles, $n^*_b(\eta, p_T)$, than $\nu_b(\eta, p_T)$ does. In fact, if the true number of soft QCD particles in each $(\eta, p_T)$ bin inside a given event deviates from the expected value by an amount $\Delta n_b$, i.e. if $n^*_b = \nu_b + \Delta n_b$, then a fraction $w_0$ of $\Delta n_b$ will contribute to $\hat{n}_b$ in that bin. On the other hand, $\nu_b$ reflects an expectation and therefore does not contain any information about statistical fluctuations inside events.

Given the estimated number of soft QCD particles in each ($\eta, p_T$) bin, $\hat{n}_b(\eta, p_T)$, the corresponding unknown actual number can be treated as a random variable following a binomial distribution with mean given by expression () and standard deviation

$$\sigma_{n_b} = \sqrt{n w_0 (1-w_0)}.
\label{eq:sigma_nbhat}$$

Based on expression (), for the purpose of estimating the background particle multiplicity inside each event, the discrimination between soft QCD interactions and the hard parton scattering exploits the different shapes of the corresponding PDFs as functions of $\eta$ and $p_T$. Specifically, as anticipated, the discrimination relies on particles from the hard scattering having a pseudorapidity distribution more peaked at $\eta=0$, as well as having on average higher values of $p_T$.

The use of expression () for $w_0(\eta, p_T)$ in order to estimate the multiplicity of soft QCD particles across the kinematic space inside the events is essentially equivalent to weighting particles against PDFs that summarise our knowledge of the kinematics of the underlying processes, thereby taking into account the average fraction of soft QCD particles in the events. As far as the hard scattering is concerned, in addition to $gg\rightarrow t\bar{t}$, which is used to illustrate our method in the present article, the algorithm has also been run on particles originating from decays of the Standard Model Higgs boson produced via vector boson fusion . In fact, such a process does not involve colour exchange between the colliding protons, and is therefore expected to lead to a lower degree of particle activity around $\eta=0$ when compared to $gg\rightarrow t\bar{t}$.

The following discussion refers to neutral particles, since the identification of charged pile-up particles is made significantly easier by the availability of information from the tracking detectors at the experiments.

It is worth pointing out that, although the signal and background PDFs were derived using simulated collision events in the context of this study, similar information can in principle be obtained using control samples from the data. As for the estimation of the soft QCD particle fraction among the neutral particles in the events, $\alpha_0$, it was decided to use the corresponding fraction of charged particles averaged over the events generated. In fact, the investigation of a more sophisticated approach including the use of a kinematic correction factor obtained from Monte Carlo showed no significant performance improvement . It should also be emphasised that the present estimate of $\alpha_0$ based on charged particles was exclusively obtained for the purpose of this investigation, and that more information will typically be available at the experiments, e.g. in the form of data on the neutral energy flow provided by the calorimeters.

Whereas $\alpha_0$ was defined as a global event-level quantity, the possibility of introducing a dependence on $\eta$ and $p_T$ is worth investigating in the future, as this could lead to further-improved results. Nonetheless, it was decided to rely on as simple an approach as possible for the purpose of this feasibility study.

An overview of this method is given in Fig. , which highlights what information is extracted from the models and what comes from the data. We will show that this approach produces a more accurate estimate of the background particle multiplicity inside the events than would be obtained using exclusively the expected number, $\nu_b$, as long as the statistical fluctuations in the data are dominated by those on the number of soft QCD particles.

# Results

We discuss the results of a proof-of-concept study on Monte Carlo data at the generator level. We used Pythia 8.176 to generate 1,000 events, each consisting of a $gg\rightarrow t\bar{t}$ hard parton scattering at $\sqrt{s} =$ 14 TeV, superimposed with 50 soft QCD interactions to simulate the presence of pile-up. Soft QCD interactions were generated with "SoftQCD:all", "PartonLevel:ISR", "PartonLevel:FSR" and "PartonLevel:MI" set to "on". The performance of this method on a reference event will be discussed for the sake of illustration, and distributions on all events generated will then be shown in order to confirm the consistency of the results.

As a prerequisite for the execution of this method, the particle-level $(\eta, p_T)$ space in each event was subdivided into bins of widths $\Delta\eta = 0.5$ and $\Delta p_T = 0.05~\mbox{GeV}/c$. Whereas our $\Delta p_T$ binning will have to be revised in the context of a full detector simulation study, which is outside the scope of this article, we are using this choice of bins as a starting point to illustrate our method. Our analysis focusses on particles with $0 < p_T < 1~\mbox{GeV}/c$, which are the majority of those produced by soft QCD interactions.

Signal and background PDF templates as functions of particle $\eta$ and $p_T$ were obtained using a control sample dataset containing ${\sim}300,000$ particles from $gg\rightarrow t\bar{t}$ and ${\sim}300,000$ from soft QCD interactions. The corresponding Monte Carlo truth ($\eta, p_T$) distributions of neutral soft QCD particles and of neutral particles from the hard scattering are shown in Fig. (a) and Fig. (b), respectively, each normalised to unit volume.

The above-mentioned collections of ${\sim}300,000$ particles, although significantly-lower statistics than the large datasets normally used in the field, are considered adequate for the purpose of estimating the signal and background probability distributions in the context of this study. In fact, our emphasis is on how the soft QCD frequency distributions inside individual events deviate from the corresponding probability distribution. For this purpose, ${\sim}300,000$ particles are high-enough statistics for the local features in the frequency distributions due to the presence of statistical fluctuations to be averaged out.

The distributions shown in Fig. were used together with the previously-mentioned estimate of the average fraction of soft QCD particles over all neutrals in the events, $\alpha_0$, in order to calculate $w_0$ according to expression (). The distribution of $w_0(\eta, p_T)$ is shown in Fig. (a) in relation to the Monte Carlo event chosen to illustrate our results.

Fig. (a) displays the true multiplicity of neutral soft QCD particles as a function of particle $\eta$ and $p_T$ in the reference event. As expected, the frequency distribution deviates from the corresponding higher-statistics distribution, shown in Fig. (a), due to the presence of local features that are typically washed out when multiple events are lumped together.

The estimate of the multiplicity of neutral soft QCD particles across the $(\eta, p_T)$ space obtained using this method is shown in Fig. (b) with reference to the same event. A comparison with the corresponding true distribution in Fig. (a) suggests that the local features of the distribution are reasonably well described, e.g. the excess at $\eta\simeq 2.5$ and $p_T\simeq 0.2~\mbox{GeV}/c$. The performance is discussed further below with reference to all the events generated.

The absolute and relative statistical uncertainties on $\hat{n}_b$ in the reference event are displayed in Fig. (b) and Fig. (a), respectively, where the absolute uncertainty, $\sigma_{n_b}$, was estimated using expression (). In particular, Fig. (b) suggests that the precision of this method is better than 1 particle, although, in order for this claim to be made, the precision over all events generated also needs to be assessed, as discussed in the following.

It is worth emphasising that, in general, it is not known which $(\eta, p_T)$ bins in the event contain signal particles and which do not. Therefore, once it is observed that $\hat{n}_b$ is a better estimate of $n^*_b$ than $\nu_b$ is, and that it is sufficiently precise, it is also necessary to verify that $\hat{n}_b$ is more accurate than the estimate that would be obtained if the possible presence of signal particles in the bins was simply neglected. For this purpose, the absolute deviation of the estimated number of neutral soft QCD particles from the true value, normalised to the true number of signal particles, $\left|\hat{n}_b-n^*_b\right|/n^*_s$, is displayed in Fig. (b) above $0.4~\mbox{GeV}/c$ particle $p_T$ in those ($\eta, p_T$) bins that contain at least 1 signal particle. As it can be seen, $\left|\hat{n}_b-n^*_b\right|/n^*_s \lesssim 1$, i.e. the absolute deviation of the estimated number of background particles from the true number is lower than the number of signal particles across the kinematic space in the event, corresponding to the pile-up rate considered.

As anticipated, Fig. (b) and Fig. relate to one single event. In order to verify the precision and the accuracy of the algorithm in more detail, the corresponding heat maps were obtained over all events generated. Fig. (a) and Fig. (b) display $\left<\sigma_{n_b}\right>$ and $\left<\left|\hat{n}_b-n^*_b\right|/n^*_s\right>$ across the $(\eta, p_T)$ kinematic space, respectively, where the average is taken over events. The plots confirm that the algorithm produces consistent results, with $\left<\sigma_{n_b}\right>$ below 1 particle and $\left<\left|\hat{n}_b-n^*_b\right|/n^*_s\right>$ significantly lower than 1 across the ($\eta, p_T$) space.

It is worth noticing that the present investigation relies on a particle-level kinematic comparison between soft QCD interactions and a specific hard scattering process, namely $gg\rightarrow t\bar{t}$. Choosing a different signal process will normally change the final-state particle kinematics, although the difference between particles originating from a hard scattering and soft QCD particles is generally expected to be more pronounced than differences across signal processes. As discussed in section , this is supported by the study documented in , where this method was applied to vector boson fusion Standard Model Higgs production. In any case, the potential dependence of the performance of this technique on the choice of signal process deserves further investigation, in order for the results presented in this article to be generalised.

# Conclusion

The contamination, or background, from particles produced by low-energy strong interactions is a major issue at the Large Hadron Collider. Although well-established correction procedures are in use at the experiments, new techniques are also being investigated in the field. The primary objective of this new line of development is to meet the requirements that will be posed by the upcoming higher-luminosity operational regimes of the accelerator.

We have investigated a different perspective to mainstream methods, thereby estimating the kinematics of background particles inside individual collision events in terms of their frequency distribution. Whereas the use of probability distributions is traditional, our emphasis on the frequency distributions is, to the best of our knowledge, a distinctive and unique feature of our approach.

Our hope and expectation is that the ability to describe the kinematic frequency distribution of the contaminating particles collision by collision will help towards the development of improved subtraction methods at higher luminosity. In particular, we expect that our method will become useful as a complement to existing correction algorithms applied to jets, i.e. to collections of final-state particles interpreted as originating from the same scattered parton.

The preliminary results discussed in this article suggest that our method can produce more accurate estimates of the number of contaminating particles in different kinematic regions inside collision events than would be possible by relying exclusively on the expected numbers. Although the possible impact of mismodelling remains to be investigated in more detail, this encourages further studies in this direction.

However, it should be stressed that a proper quantitative test of the proposed method will require a detailed assessment on specific observables, which is beyond the scope of the present feasibility study.

It should also be emphasised that the algorithm is inherently parallel. In fact, different kinematic regions inside events can be processed independently, and the calculation of the only global variable, i.e. of the average fraction of contaminating particles per event, can be performed in advance on a control sample. For this reason, this method is potentially suitable for inclusion in future particle-level event filtering procedures upstream of jet reconstruction.

The possible complementarity between our method and the particle weighting techniques that have recently been proposed at the Large Hadron Collider will also be worth exploring. In fact, our estimate of the probability for individual particles to originate from low-energy QCD processes as opposed to the high-energy signal scattering is based on information that is currently not employed by any particle weighting algorithms. In this context, we anticipate that multivariate combinations of different weighting schemes can prove beneficial with a view to improving further on the rejection of contaminating particles.

Finally, it will be useful to study the impact of this method, when used in conjunction with existing techniques, on the resolution of the missing transverse energy as well as on estimates of particle isolation, with a view to quantifying the associated benefit at the level of physics analysis. The source code used for the purpose of this study is made available upon request.

# Competing interests

The authors declare that there is no conflict of interest regarding the publication of this paper.

[^1]: Email: firstname.lastname@example.com