Tagging the initial-state gluon

We study the production of an electroweak boson in association with jets, in processes where the jet with the highest transverse momentum is identified as quark-initiated. The quark/gluon tagging procedure is realised by a cut on a jet angularity and it is therefore theoretically well-defined and exhibits infrared and collinear safety. In this context, exploiting resummed perturbation theory, we are able to provide theoretical predictions for transverse momentum distributions at a well-defined and, in principle, systematically improvable accuracy. In particular, tagging the leading jet as quark-initiated allows us to enhance the initial-state gluon contribution. Thus, these novel transverse momentum distributions are potentially interesting observables to probe the gluonic degrees of freedom of the colliding protons.


Introduction
Over the past decade our understanding of the internal structure of jets has increased tremendously. Thanks to the applications of the methods of perturbative QCD, the field has become mature and jet substructure algorithms that are, at the same time, performant and robust have been developed. Observables originally designed for searching new physics are now the target of unfolded measurements that can be compared to theoretical predictions, at the precision level, see e.g. [1][2][3][4][5][6][7][8]. Furthermore, ideas developed by the jet substructure community have found applications in other contexts of particle physics. For instance, after it was realised that the momentum fraction that characterises the splitting identified by the SoftDrop algorithm [9] follows a distribution dictated by the QCD splitting functions [10][11][12][13], this a e-mail: simone.caletti@ge.infn.it b e-mail: oleh.fedkevych@ge.infn.it c e-mail: simone.marzani@ge.infn.it d e-mail: daniel.reichelt@uni-goettingen.de observable has become a standard way of probing interactions with the quark/gluon plasma, see e.g. [14] and references therein.
This work is part of an incipient effort to find innovative ways of applying jet substructure techniques to the broader LHC physics program and to provide the necessary tools so that this cross-pollination can bear its fruits. In this context, measurements of the internal structure of highly energetic jets produced in association with a boosted electroweak boson offer an ideal playground for such studies, because the leptonic decay of the Z boson offers a valuable trigger. These events are not only relevant for background studies 1 but, as we shall detail in the following, they open up novel possibilities to employ jet substructure for Standard Model measurements.
Substructure observables, such as jet angularities [15], measure the pattern of the hadronic activity within a jet. For this reason, they are often employed as tagging variables that aim to distinguish jets that have been originated by elementary particles carrying different colour degrees of freedom, e.g. colour singlets versus QCD partons or, even, quarks versus gluons. Although there exist more powerful quark/gluon (q/g) discriminants than a cut on a jet angularity, this procedure is theoretically well-defined, infra-red and collinear (IRC) safe and, hence, its behaviour can be understood with perturbative methods.
In the following, we study Z+jet production, requiring that the jet with the highest transverse momentum has been identified as quark-initiated. As preliminarily explored in [16], by tagging the final-state, we indirectly bias particular partonic sub-processes. If the leading jet is tagged as quarkinitiated, then, at leading order, the subprocess that features a quark and gluon in the initial state is enhanced. Thus, this procedure could potentially provide us with a new handle on the determination of the gluon parton distribution function (PDF). The main result of this study is the calculation of the transverse momentum distribution of the Z boson in events where the leading jet is identified as quark-initiated. Because of IRC safety of the tagging procedure, we are able to compute this distribution at a well-defined and, in principle improvable, accuracy. Therefore, this observable could be directly included in standard fits of PDFs. The results presented here account for the resummation of the tagging parameter at the next-to-leading logarithmic (NLL) accuracy, matched to fixed order predictions at O α 2 s with respect to the Born contribution, henceforth denoted as next-to-leading order (NLO) 2 .

Enhancing the gluon contribution
We start by considering the channel fractions f i j that measure, to a given order in perturbation theory, the contribution to the Z (or jet) transverse momentum distributions from the subprocess initiated by partons i and j, which, for brevity, we indicate as σ a i j , where a could be the Z boson or the leading jet J: where qq and qg include any combination of quarks and anti-quarks. Our aim is to study how tagging a quark-jet in the final-state changes the relative contributions of the various partonic subprocesses. To this purpose, we define where the tilde indicates "after tagging". Furthermore, in what follows we will be interested in the fractions of the event (before and after tagging) which feature at least one gluon in the initial state. To this purpose, we define f a g = f a qg + f a gg , and f a g = f a qg + f a gg . (3) We will refer to this fraction as gluon channel purity. We begin our discussion in a simplified setting, which is nevertheless enough to capture the essential physics points. We will then validate our conclusions using Monte Carlo (MC) parton shower simulations. If we consider the Born approximation, then we have no gg contribution and, because the Z boson and the jet are back-to-back, we obtain 2 The fixed-order counting is slightly different than what is usually employed in standard transverse-momentum distribution. This has to do with the fact that, in order to obtain a non-vanishing value of any angularity, the jet must have at least two constituents. Thus, in what follows we will refer to the lowest-order (2 → 2) scattering as Born approximation, while (N)LO will be reserved for the contributions with one (two) additional emission(s).
Fragmentation of the hard partons can lead to transverse momentum imbalance, however, we can limit ourselves to the leading logarithmic (LL) regime, where all parton splittings happen in the soft and collinear limit, with no recoil. 3 Therefore, dropping the superscript a, we have We note that, within our approximation, the fraction of events with a (properly defined) final-state quark can be considered as a proxy for the gluon channel purity. We find that, for p t J ≥ 100 GeV, f g 0.85 and it exhibits a rather mild dependence on the transverse momentum. Next, we note that in our approximation we simply have where ε q is the efficiency of the tagger to correctly label quark jets and ε g the false-positive rate, which measures how often gluon jets are wrongly labelled as quarks. We note that a perfect quark-tagger with ε q = 1 and ε g = 0 returns tagged events that have been originated by one gluon in the initial state, i.e. f g = 1. On the other hand, with an efficiency of 50%, which corresponds to tossing a coin, we recover Eq. (4). A rather common class of q/g taggers exhibits at LL a property known as Casimir scaling (CS) [15], namely the tagging efficiencies are related by ε g = (ε q ) C A /C F , where the exponent is given by the ratio of the Casimir operators in the appropriate colour representation, C F for quarks and C A for gluons. This property emerges because the LL distributions of both quarks and gluons exhibit the same Sudakov-like behaviour, with a coefficient determined by the appropriate colour factor. Thus, for a CS tagger, we find For instance, a CS tagger with ε q = 0.65, yields f a g 0.9. This can be pushed to 0.95 if the tighter working point ε q = 0.35 is considered. Taggers that obey CS are not the most performant, but they are interesting for us because they are under good theoretical control. The efficiencies ε i can be computed in QCD using resummed perturbation theory to a well-defined and, in principle improvable, theoretical accuracy. Furthermore, the inclusion of higher logarithmic corrections generally leads to an improvement with respect to strict CS, so that, depending on the specifics of the tagging procedure, the target gluon channel purity f a g 0.95 can be achieved at a reasonable working point. In the following we construct a CS tagger that is based on a particular class of jet substructure observables known as jet angularities. The dashed portion of these lines indicate the region sensitive to splittings with relative transverse momentum below 1 GeV, where NP effects are expected to be sizeable. This region of phase-space is sensitive to shower cut-off effects, which cause the observed kinks.
Jet angularities as a quark/gluon tagger Jet angularities [15] are defined as where the sum runs over the constituents of the hardest jet in the event and ∆ i = (y i − y J ) 2 + (φ i − φ J ) 2 is the distance in the azimuth-rapidity plane of particle i from the jet axis. We define jets with the anti-k t clustering algorithm [17] with radius R 0 and standard E-scheme for recombination. IRC safety requires α > 0, while angularities with α ≤ 1 are sensitive to recoil against soft emissions [18]. In order to circumvent this issue, when α ≤ 1, the jet axis is obtained using the Winner-Take-All (WTA) recombination scheme [19]. We also consider groomed jets. In this case, we recluster the jet with the Cambridge-Aachen algorithm [20,21] and apply the SoftDrop grooming algorithm with parameters z cut and β [9]. The angularity is then computed on the constituents of the groomed jet, with the WTA prescription adopted for angularities with α ≤ 1.
Because of the different colour factor, angularity distributions for quark-and gluon-initiated jets peak at different values. We can exploit this separation and define our q/g tagger through a cut on the jet angularity. In particular, a jet with λ α < λ cut will be labelled as a quark jet: where we have introduced the differential angularity distribution. The implicit dependence of the tagged distribution on the transverse momentum, on the angularity exponent α and, optionally, on the SoftDrop parameters z cut and β , is understood. From the simple CS analysis above, we have concluded that we should work with quark efficiencies ε q 0.35, in order to reach a purity of initial-state gluons around 0.95. The value of the angularity cut that is necessary to achieve this working point for the tagger clearly depends on the angular exponent α in Eq. (7) as well as on the parameters z cut and β of the SoftDrop algorithm, should we wish to employ groomed jets. Different theoretical considerations can guide us with this choice. First of all, we would like to preserve calculability, i.e. we want to cement our findings in perturbative field theory. Thus, we would like our tagger to be as insensitive as possible to non-perturbative (NP) contributions such as hadronisation corrections and the Underlying Event (UE). Secondly, although we can calculate transverse momentum spectra in resummed perturbation theory [22][23][24], we aim for perturbative stability. Thus, we favour working points for the tagger for which λ cut is not too small. In order to turn the above considerations into a quantitative study we use simulated data obtained with the MC event generator PYTHIA 8.303 [25]. The UE is simulated according to the model presented in [26][27][28] and hadronisation effects according to the Lund string model [29,30]. Throughout the paper, we use the NNPDF 3.0 NLO set of PDFs [31]. We consider the inclusive production of a pair of oppositely charged muons in proton-proton collisions at 13 TeV centre-of-mass energy, requiring that the invariant mass of the muon pair to be within 70 and 110 GeV. Jets are clustered with the anti-k t algorithm with R 0 = 0.4 and ordered in transverse momentum. Henceforth, the jet will be implicitly considered to be the hardest one and we will refer to the muon-antimuon pair as the Z boson. With this in mind, the fiducial volume is defined following Ref.
[8]: p t µ > 26 GeV, p t Z > 30 GeV, p t J > 15 GeV, |η µ | < 2.4, and |y jet | < 1.7. Furthermore, in order to enforce back-to- For event selection and analysis we employ RIVET [32,33]. Jet reconstruction is done with FASTJET [34], and the SoftDrop implementation in the FASTJET contrib is used. We distinguish two different stages of the simulation: "partonlevel", i.e. with parton shower effects only, and "hadronlevel", i.e. with hadronisation and UE included. By keeping the two partonic processes of interest separate, we compute Receiver Operating Characteristic (ROC) curves that show the mis-tag rate (gluon efficiency) ε g as a function of the signal (quark) efficiency. They are computed for different values of the angularity exponent α = 0.5, 1, 2 in the ungroomed case and for SoftDrop jets with z cut = 0.1 and β = 0, 1. We show hadron-level results in Fig. 1, as well as the ratios to their parton-level counterparts, which we take as a measure of NP contributions. The dotted portion of the latter indicates that the efficiency ε q is dominated by splittings in the non-perturbative region, as determined, for instance, in Ref. [24]. We also show the target line, which fixes the gluon efficiency ε g as a function of ε q , for given f g and f g , which is easily derived from Eq. (5). The slope shown in Fig. 1 is determined by the original f g = 0.85 and target gluon purity f g = 0.95. The intersections of each ROC curve with the target line set our tagger working points. Corresponding values of λ cut are reported in the figure. The choice of the tagger working points is clearly not unique. For instance, we could have optimised the signal-to-background ratio by choosing on each ROC curve the point that is closest to the (1, 0) corner. The analysis of Ref. [24] tells us that larger values of α are under better theoretical control. However, as it is clear from Fig. 1, lower values of α have better performance, essentially because of their increased sensitivity to the collinear region. Thus, the choice α = 1 appears to be a good compromise between performance and robustness. We note that the use of SoftDrop does not always provide us with improvements on the size of the NP contributions. This might be related to the fact that in order to obtain the same efficiency ε q we need to cut groomed jets at lower values of λ cut , where NP physics may be more prominent. Furthermore, we note that we are working with a rather small jet radius, which prevents large contributions from the UE. We expect (light) grooming to be beneficial, should one consider larger jet radii. Finally, configurations with aggressive SoftDrop, i.e. β = 0, typically result in worse performance, because important information is groomed away. Thus, in what follows, we shall focus on the α = 1 case either with no grooming or with β = 1, z cut = 0.1.
We now quantify the gain in the initial-state gluon purity that we obtain after tagging. Fig. 2 shows the fractionsf Z g as a function of the Z boson transverse momentum for the selected taggers, with f Z g also shown for comparison. We first note that the performance of both taggers is very good, leading to gluon channel purities that exceed our 95% target. Analogous conclusions can also be drawn if we plot our results as a function of p t J . However, as we will argue shortly, in this context, the transverse momentum distribution of the Z boson is a more robust observable. We also notice that the gluon channel purities decrease with p t Z . This is due to the fact that our taggers are defined looking at their efficiencies with p t J > 100 GeV and, therefore, we expect them to work better at the lower end of the transverse momentum spectrum. This loss in performance could be cured by adjusting the cut on the angularity as a function of p t J . However, in this first study, we prefer to keep our framework simple.

Transverse momentum distributions
We now provide theoretical predictions for our observables of interest, namely transverse-momentum distributions, in the presence of tagging. Our calculation includes the resummation of logarithms of λ cut at NLL accuracy and is matched to NLO. Thanks to a flavour-dependent matching procedure, cf. also [35,36], we are able to achieve NLO+NLL accuracy. We also include a bin-by-bin NP correction factor obtained with MC simulations. 4 Our results are reported in Fig. 3, where we show the p t Z distribution for events where the highest-p t jet is quark initiated, for the taggers selected for this study. We show 4 See Ref. [24] for details about the calculation and its numerical implementation in the resummation plugin [37,38] to the SHERPA [39,40] framework, including perturbative uncertainties, obtained by varying the perturbative (renormalisation, factorisation and resummation) scales and NP corrections. We use COMIX [41] in conjunction with OpenLoops [42] and Recola [43,44] for the fixed order calculation. The NP corrections are based on SHERPA parton shower simulations at MC@NLO accuracy [45,46] hadronised with SHERPA's cluster fragmentation model [47]. the p t Z distribution with no-tagging, with tagging on standard jets and with tagging on SoftDrop jets. We note that NP corrections are rather sizeable at low p t Z , making this observable most reliable in the high transverse momentum region. The latter is actually per se interesting because it allows us to probe the proton dynamics described by the PDFs at large values of the momentum fraction x, i.e. in a kinematic region where they are less constrained. To qualitatively assess the values of x which would be accessible, we show, on the upper horizontal axis, the Born-level momentum fraction computed at central (zero) rapidities:x = .
We have also studied p t J distributions. However, as anticipated, these distributions turn out to be less robust, essentially because jet dynamics can be significantly altered by the cut on the angularity, as well as by the grooming procedure 5 . In contrast, the p t Z spectrum is inclusive with respect to the jet activity and thus, for a given working point of the tagger, less dependent on the details of the tagging procedure itself. We can therefore take modifications in such distribution as more directly related to the bias that the tagger induces on the composition of the partonic initial state, which is what we want to achieve.

Conclusions and future developments
We have shown how q/g tagging can be successfully applied to Z+jets events in order to significantly enhance the gluoninitiated contributions. In particular, our tagger is realised through a simple cut on a jet angularity, which is an IRC safe observable and therefore can be studied in perturbation theory. Exploiting MC simulations, we have performed a study of their efficiencies and their dependence on NP effects, exploring different angularities and different levels of grooming. We have explicitly studied the transverse momentum of the Z boson, providing theoretical predictions that included both all-order resummation and matching to NLO. We have shown that we can achieve initial-state gluon purities close to 95%.
We see several possible directions for future work in this context. First, we would like to assess the impact of this type of observables on PDF fits. In particular, while our study does show an increase in the gluon purity, this comes at the cost of a noticeable reduction of the available dataset. Thus, despite the relative large cross-section of the process we are considering, we may need to optimise the tagger's working point, taking this further aspect into consideration. For instance, as mentioned above, we could choose λ cut so that the signal-to-background ratio is maximised. Second, the results presented in this study are based on a calculation of the angularity spectra, and hence of closely related efficiencies, at NLO+NLL . The resummed calculation can be promoted to higher accuracy [48][49][50]. The fixed-order can also be improved by including the two-loop correction to Z+1 jet (see [51] and references therein). An improvement in the description of the angularity distribution away from λ α = 0 is much more challenging because it requires Z+2 partons at NNLO, which may become available in the near future. Finally, with the aim of enhancing performance while maintaining calculability, it would be interesting to consider more powerful, albeit more sophisticated, taggers. In this context, the Les Houches multiplicity [16] is very promising and, although its theoretical understanding is only in its infancy, we believe that achieving NLL accuracy is within reach.