Improving the Simulation of Quark and Gluon Jets with Herwig 7

The properties of quark and gluon jets, and the differences between them, are increasingly important at the LHC. However, Monte Carlo event generators are normally tuned to data from $e^+e^-$ collisions which are primarily sensitive to quark-initiated jets. In order to improve the description of gluon jets we make improvements to the perturbative and the non-perturbative modelling of gluon jets and include data with gluon-initiated jets in the tuning for the first time. The resultant tunes significantly improve the description of gluon jets and are now the default in Herwig 7.1.


Introduction
Monte Carlo generators are essential tools, both for the design of future experiments and the analysis of data from the LHC, and previous collider experiments. Modern event generators [1][2][3] provide a simulation of exclusive events based on the combination of fixed-order perturbative results, resummation of large logarithms of scales using the partonshower approach and non-perturbative models of hadronization and multiple-parton scattering. 1 These simulations rely on universality and factorization in order to construct a simulation of the complex final states observed in hadronic collisions. This allows the simulation of final-state radiation in the parton shower and the nonperturbative hadronization models to be first developed, and the parameters of the the model tuned, using the simpler and cleaner environment of e + e − collisions, and then applied to more complicated hadronic collisions. These models are then combined with the parton-shower simulation of initialstate radiation, a multiple scattering model of the underlying a e-mail: daniel.reichelt@tu-dresden.de b e-mail: peter.richardson@durham.ac.uk c e-mail: andrzej.siodmok@ifj.edu.pl 1 For a recent review of modern Monte Carlo event generators see [4]. event and a non-perturbative colour reconnection model in order to describe hadronic collisions. In principle universality requires that the colour reconnection model is also used to describe leptonic collisions. In practice however colour reconnection has little effect on the distributions which so far have been used to develop and tune the models. These models are therefore usually either not included at all for the simulation of leptonic collisions, or if they are, the parameters are determined by tuning to hadronic data sensitive to multiple partonic scattering.
As the LHC accumulates data at an unprecedented rate there are a number of observables which are not well described by current Monte Carlo event generators, and where the limitations of this approach have started to become obvious, for example: the difference in the properties of jets initiated by quarks and gluons is not well described with generators predicting either a larger or smaller difference between the jets than is observed by the LHC experiments [5]; the transverse momentum spectra of identified baryons and strange hadrons which are not well described by current generators. [6]; long-range correlations in high multiplicity events [7,8].
In this paper we will focus on improvements to the perturbative and non-perturbative modelling to give a better description of both quark-and gluon-initiated jets, as well as the differences between them in Herwig 7. Beyond leading order there is no clear distinction between quark and gluon jets and the definition will depend on the analysis. 2 As e + e − annihilation to hadrons starts with an initial partonic quark-antiquark configuration the data used to develop the final-state parton-shower algorithm, tune its parameters and those of the hadronization model, are dominated by quarkinitiated jets. However at the LHC jets initiated by gluons can often dominate, depending on the production process, rapidity and transverse momentum of the jets. Regrettably while there is great interest in the differences between quark and gluon jets at the LHC most of the experimental studies have concentrated on differentiating between quark and gluon jets using neural network, or similar, techniques which makes a direct comparison with simulated hadron-level events impossible. We will therefore use some recent data from the ATLAS experiment [10] which is sensitive to both quark and gluon jet properties, together with data on gluon jets in e + e − collisions from the OPAL experiment [11,12] which has not previously been used in the development and tuning of the current generation of Monte Carlo event generators to study the properties of gluon jets.
In the next section we will first recap the default partonshower algorithm used in Herwig 7 focusing on recent changes we have made to improve the simulation of both quark and gluon jets. In Section 3 we will briefly review the important parameters in the cluster hadronization model used in Herwig 7 and identify the issues which may lead to different treatments of quark and gluon jets. We will then discuss the tuning strategy used to produce the tunes presented in this paper. We present our results in Section 5 3 followed by our conclusions.

Herwig 7 Parton-Shower Algorithm
The default Herwig 7 parton-shower algorithm [13] is an improved angular-ordered parton shower. In this approach the momenta of the partons produced in the parton shower are decomposed in terms of the 4-momentum of the parton initiating the jet, p (p 2 = m 2 , the on-shell parton mass-squared), a light-like reference vector, n, in the direction of the colour partner of the parton initiating the jet and the momentum transverse to the direction of p and n. The four momentum of any parton produced in the evolution of the jet can be decomposed as where α i and β i are coefficients and q ⊥i is the transverse four momentum of the parton (q ⊥i · p = q ⊥i · n = 0). If we consider the branching of a final-state parton i to two partons j and k, i.e. i → jk as shown in Fig. 1, the branching is described by the evolution variablẽ where q 2 i is the square of the virtual mass developed by the parton i in the branching, m i is the physical mass of parton i, 3 Additional results on quark and gluon jet discrimination power are included in the Appendix.q iq k q j z i 1 − z i Fig. 1 Branching of the parton i to produce the partons j, k which then undergo subsequent branching. and z i is the momentum fraction of the parton j defined such that The transverse momenta of the partons produced in the branching are where k ⊥i is the transverse momentum generated in the branching. In this case the virtuality of the parton i is where p T is the magnitude of the transverse momentum produced in the branching defined such that k 2 ⊥i = −p 2 Ti . In this case the probability for a single branching to happen is where P i→ jk (z,q) is the quasi-collinear splitting function, and φ i is the azimuthal angle of the transverse momentum k ⊥i generated in the splitting. As the branching probability is singular for massless partons an infrared cut-off is required to regularise the singularity. In HERWIG 6 [14] and early versions of Herwig++ [15] the cut-off was implemented by giving the partons an infrared mass. However while this remains an option in later versions of Herwig++ and Herwig 7 [1] the default cut-off is now on the minimum transverse momentum of the branching [16].
In order to resum the dominant subleading logarithms [17] the transverse momentum of the branching is used as the scale for the strong coupling constant. This also means that the strong coupling used in the parton shower is that defined in the Catani-Marchesini-Webber (CMW) scheme which includes the subleading terms via a redefinition of QCD scale, Λ QCD .
While this specifies both the branching probability and kinematics of the partons for a single emission in the case of subsequent emission from the daughter partons j and/or k we must decide which properties of the originally generated kinematics to preserve once the masses of j and/or k in Eqn. 5 are no longer the infrared cut-off masses but the virtualities generated by any subsequent emissions. While this choice is formally subleading it can have a large effect on physical observables.
In Herwig++ the transverse momentum of the branching was calculated using Eqn. 5 and the infrared cut-off masses when the emission was generated and then preserved during the subsequent evolution of the daughter partons. In Herwig 7.0 the default option was to instead preserve the virtuality of the branching and calculate the transverse momentum of the branching using the virtual masses the daughter partons develop due to subsequent emissions. This means that if the daughter partons develop large virtual masses the transverse momentum of the branching is reduced, and in some cases the branching has to be vetoed if there is no solution of Eqn. 5. However, this choice inhibits further soft emission and significantly changes the evolution by vetoing emissions and leads to incorrect evolution of observables. We therefore consider a further choice in which if it is possible to preserve the virtuality and still have a solution for p 2 T > 0 we do so, however if this is not possible instead of vetoing the emission we set p T = 0 and allow the virtuality to increase.
The most important parameters which affect the behaviour of the parton shower and which we will tune in this paper are: the choice of whether to preserve p T or q 2 during the subsequent evolution; the value of the strong coupling constant AlphaMZ, taken to be α CMW S (M Z ), value of the coupling constant in the CMW scheme at the mass of the Z boson, M Z ; the cut-off in the parton shower 4 . For a cut-off in p T this is the minimum transverse momentum allowed for the branchings in the shower, p min T . For a virtuality cut-off we parameterize the threshold for different flavours as where a and b are parameters chosen to give a threshold which is slightly reduced for heavier quarks. The parameter c = 0.3 GeV is chosen to prevent the cutoff becoming too small, we also keep the default value of b = 2.3. Only the parameters δ (cutoffKinScale) and a (aParameter) are tuned to the data. There is one other major feature of the angular-ordered parton shower which we need to consider. The angular ordering of the parton shower, which is used to implement the phenomenon of colour coherence, leads to regions of phase space in which there is no gluon emission. Consider for example the process e + e − → qqg. In this case there is a deadzone which is not filled by one emission from the parton shower, as shown in Fig. 2. Given this deficit of hard, wideangle emission it is necessary to combine the parton-shower with the fixed-order calculation of e + e − → qqg. There are now a range of techniques which can achieve this including both the next-to-leading order normalization of the total cross section, or including the fixed-order results for multiple emissions. However, for our purposes it is sufficient to consider the simplest matrix-element correction approach where the dead-zone is filled using the leading-order matrix element for e + e − → qqg, as shown in Fig. 3, together with the reweighting of emission probability, Eqn. 6, to the exact leading-order result, for any emission which could have the highest transverse momentum in the parton shower. 5 The choice of whether to preserve the transverse momentum or virtuality of the branching affects the phase-space region which is filled by the shower in the case of multiple emission. In this case we cluster the partons using the Durham jet algorithm [19], using the p-scheme as implemented in FastJet [20], keeping track of the partons emitted by the quark and antiquark and then take the hardest additional jet to be the gluon. The resulting Dalitz plots of e + e − → qq show that while the choice to preserve the transverse momentum of the branching leads to a significant number of events in the dead-zone, Fig. 4, if the virtuality of the branching is preserved, Fig. 5, there is little emission outside the original angular-ordered region.
where E i is the energy of parton i and Q is the centre-of-mass energy of the collision.

Hadronization and Colour Reconnection
All the Herwig family of event generator generators use the cluster hadronization model [21]. This model is based on the phenomena of colour pre-confinement, i.e. if we nonperturbatively split the gluons left at the end of the parton shower into quark-antiquark pairs and cluster quarks and antiquarks into colour-singlet clusters the mass spectrum of these clusters is peaked at masses close to the cut-off in the parton shower, falls rapidly as the cluster mass increases, and is universal, i.e. the mass distribution of these clusters is independent of the hard scattering process and its centre-ofmass energy. The cluster model assumes that these clusters are a superposition of heavy hadronic states and uses a simple phase-space model for their decay into two hadrons. The main parameters of the model are therefore: the non-perturbative gluon mass, which is not very sensitive and we do not tune; the parameters which control the probability of producing baryons and strange quarks during cluster decay; the parameter which controls the Gaussian smearing of the direction of the hadrons produced which contain a parton from the perturbative evolution about the direction of that parton, with separate values for light, charm and bottom quarks.
There are however a small fraction of large mass clusters for which the two hadron decay ansatz is not reasonable and these must first be fissioned into lighter clusters. While only a small fraction of clusters undergo fission due to the larger masses of these clusters they produce a significant fraction of the hadrons.
A cluster is split into two clusters if the mass, M, is such that where Cl max and Cl pow are parameters of the model, and m 1,2 are the masses of the constituent partons of the cluster. For clusters that need to be split, a qq pair is selected to be popped from the vacuum. The mass distribution of the new clusters is given by where m q is the mass of the parton popped from the vacuum, M 1,2 are the masses of the clusters formed by the splitting and R 1,2 are pseudo-random numbers uniformly distributed between 0 and 1. The distribution of the masses of the clusters is controlled by the parameter P split .
In order to improve the description of charm and bottom hadron production these parameters for cluster fission all depend on the flavour of the partons in the cluster so that there are separate parameters for light, charm and bottom quarks.
In practice there is always a small fraction of clusters that are too light to decay into two hadrons. Before Herwig 7.1 these clusters were decayed to the lightest hadron, with the appropriate flavours. However in some cases, for example for clusters containing a charm or bottom quarkantiquark pair, or a bottom quark and a light antiquark, there can be a number of hadrons of the appropriate flavour below the threshold. In these cases the lightest meson with the appropriate flavours is the pseudoscalar 1 S 0 state and the vector 3 S 1 state is also below the threshold 6 which leads to a lower production rate for the vector state with respect to the pseudoscalar state than expected. For the mesons composed of a bottom quark and a light quark the rate is significantly less than that expected from the counting of spin states, or indeed observed experimentally [22][23][24][25]. For charmonium and bottomonium states as this mechanism is the only way the vector states can be produced via hadronization it leads to a complete absence of direct J/ψ and ϒ production. In Herwig 7.1 we therefore include the possibility that instead of just producing the lightest state all states below the threshold are produced with a probability proportional to 2S + 1, where S is the spin of the particle.
In order to improve the behaviour at the threshold for charm and bottom clusters the option exists of allowing clusters above the threshold mass, M threshold , for the production of two hadrons to decay into a single hadron such that a single hadron can be formed for masses where SingleHadronLimit is a free parameter of the model. The probability of such a single-meson cluster decay is assumed to decrease linearly for M threshold < M < M limit and there are separate parameters for charm and bottom clusters. In order to explain the rising trend of p t vs N ch (average transverse momentum as a function of the number of charged particles in the event) observed already by UA1 [26] and describe Underlying Event [27][28][29][30] and the Minimum Bias data [31][32][33][34], the hadronization model is supplemented with a model of colour reconnections (CR) [35]. The default version of the model implemented in Herwig 7.0 is not very sophisticated. The colour reconnection model defines the distance between two partons based on their invariant mass, i.e. the distance is small when their invariant mass (cluster mass) is small. The aim of the CR model is to reduce the colour length λ ≡ ∑ N cl i=1 m 2 i , where N cl is the number of clusters in an event and m i is the invariant mass of cluster i. The colour reconnection of the clusters leading to a reduction of λ is accepted with a given probability which is a parameter of the model. Although the default model is quite simple it should be stressed that its results resemble the more sophisticated statistical colour reconection model [35] which implements the minimization of λ as Metropolis-like algorithm and requires a quick "cooling" of the random walk.
In this model the only possible reconnections which are not allowed are connecting the quark and antiquark produced in the non-perturbative splitting of the gluon. It is therefore possible that the colour lines of a gluon produced at any other stage of the shower can be reconnected leading to the production of a colour-singlet object. While this is physically possible we would expect that it occurs at a rate which is suppressed in the number of colours, N C , as ∼ 1 , not the much higher reconnection rate ∼ 2/3 7 which is necessary to describe the underlying event data. This can lead to the production of a colour-singlet gluon jet at a much higher rate than expected. This is particularly problematic in the theoretically clean, but experimentally inaccessible, colour-singlet gluon pair production processes often used to study gluon jets [9].
Consider, for example, the simple process of coloursinglet gluon pair production followed by the branching of all the gluons via g → gg, shown in Fig. 6a. After the nonperturbative splitting of the gluons into quark-antiquark pairs, as shown in Fig. 6b, without colour reconnection the quarks and antiquarks will be formed into colour-singlet clusters as (q 1 ,q 3 ), (q 3 ,q 4 ), (q 4 ,q 2 ) and (q 2 ,q 1 ). Given the configuration it is likely that the clusters containing partons from the 7 The value from the tune of Herwig 7.1 with a new soft and diffractive model [36]. and (q 4 ,q 2 ), will have large masses and the rearrangement to give the clusters (q 1 ,q 2 ) and (q 4 ,q 3 ) will be kinematically favoured, although it means the original gluons will effectively become colour singlets rather than octets. In Herwig 7.1 we have therefore included the possibility to forbid the colour reconnection model making any reconnection which would lead to a gluon produced in any stage of the parton-shower evolution becoming a colour-singlet after hadronization. We will investigate the effect of this change on the simulation of quark and gluon jets.

Tuning
The Rivet [37] program was used to analyse the simulated events and compare the results with the experimental measurements. The Professor program [38] was then used to interpolate the shower response and tune the parameters by minimising the chi-squared. 8 In general we use a heuristic chi-squared function where p is the set of parameters being tuned, O are the observables used each with weight w O , b are the different bins in each observable distribution with associated experimental measurement R b , error ∆ b and Monte Carlo prediction f b (p). Weighting of those observables for which a good description of the experimental result is important is used in most cases. The parameterisation of the event generator response, f (p), is used to minimize χ 2 and find the optimum parameter values. We take w O = 1 in most cases except for the particle multiplicities where we use w O = 10 and total charged particle multiplicities where we use w O = 50. This ensures that particle multiplicities influence the result of the fit and are required due to the much higher quantity of event 8 While tuning the parameters sensitive to bottom quarks it proved impossible to get a reliable interpolation of the generator response with Professor and therefore a random scan of the bottom parameters was performed and the values adjusted by hand about the minimum to minimise the χ 2 .
shape and spectrum data used in the tuning. Given the aim of this paper is to improve the description of gluon jets this data was also included with w O = 10 in order to avoid the fit being dominated by the large quantity of data sensitive to quark jets. In addition as we do not except a Monte Carlo event generator to give a perfect description of all the data and in order to avoid the fit being dominated by a few observables with very small experimental errors we use rather than the true experimental error, ∆ b , in the fit. The standard procedure which was adopted to tune the shower and hadronization parameters of the Herwig++ and Herwig 7 event generators to data is first the shower and those hadronization parameters which are primarily sensitive to light quark-initiated processes are tuned to LEP1 and SLD measurements of event shapes, the average charged multiplicity and charged multiplicity distribution, and identified particle spectra and rates which only involve light quark mesons and baryons; the hadronization parameters for bottom quarks are tuned to the bottom quark fragmentation function measured by LEP1 and SLD together with LEP1 and SLD measurements of event shapes and identified particle spectra from bottom events; the hadronization parameters involving charm quarks are then tuned to identified particle spectra, from both the Bfactories and LEP1, and LEP1 and SLD measurements of event shapes and identified particle spectra from charm events; the light quark parameters are then retuned using the new values of the bottom and charm parameters together with different weights for the charged multiplicity distributions in e + e − collisions at energies between 12 GeV and 209 GeV due to the difficulty in fitting the charged multiplicity.
Only e + e − annihilation data from the continuum region near the ϒ (4s) meson, for charm meson spectra, and at the Z-pole from LEP1 and SLD were used in the tune.
In this paper we have extended this approach in order to better constrain the energy evolution to include data from a wider range of centre-of-mass energies both below the Zpole, from the JADE and TASSO experiments, and above the Z-pole, from LEP2.
In order to study the various effects we have discussed we have produced tunes for the shower and hadronization parameters in the case that either the transverse momentum or virtuality in the shower is preserved. In each case we first tuned the shower and light quark parameters without the data on charged particle multiplicities as centre-ofmass energies below the mass of the Z 0 boson. In the final stage of the process where we retune these parameters three tunes were produced for each choice of cut-off and preserved quantity, one (labelled A) without the low-energy charged multiplicity data, one (labelled B) where all the charged multiplicity data was included with in the tune with weight w O = 100 and a final tune (labelled C) where this data had weight w O = 1000.
Unfortunately due to the CPU time required it is impossible to include the ATLAS data [10] directly in the tune, therefore we compare the results of the different tunes to this data.

Results
We have produced 12 tunes for different choices of the cutoff variable in the shower, the choice of which quantity to preserve in the parton shower, and different weightings of the charged particle multiplicities. The parameters obtained in the fits are given in Table. 1 while the χ 2 values are given in Table. 2.
The effects of changing the colour reconnection model can be seen in Fig. 7. In the results of Herwig++ 2.7.1 or Herwig 7.0 there is an unphysical tendency of the gluon jets to contain an even number of charged particles due to the production of colour-singlet gluons by the reconnection model, this feature is not present in any of the new tunes which provide a much better description of the distribution of charged particles in the gluon jets, see also the Appendix.
The choice of which tune and choice of cut-off variable and preserved quantity has to be a balance between how well we wish to describe the various different data sets, as unfortunately no choice provides a good description of all the data sets.
If we first consider the choice of cut-off it is clear that using a virtual mass provides a larger χ 2 for all sets of observables used in the tuning apart from those sensitive to bottom quarks. In addition it displays an unphysical energy dependence in the difference in charged particle multiplicities between bottom (or charm) quark and light quark events, as shown in Fig. 8 where the results which use a cutoff on the virtual mass, Herwig++ 2.7.1 and the new tune q 2 -q 2 -B, show a strong dependence on the centre-of-mass energy while those which use a p ⊥ cut-off, Herwig 7.0 and the new tune p ⊥ -q 2 -B, are relatively independent of energy. We therefore prefer a cut-off on the minimum transverse momentum of the branching.
In order to obtain a reasonable evolution of the number of charged particles with centre-of-mass energy in e + e − collisions, see Fig. 9, without ruining the description of particle spectra and event shape observables we choose to use the B tune as our default.
The choice of whether to preserve the p ⊥ or q 2 of the branching is more complicated. While the data on light quark jets, in particular event shapes measured at LEP (for example the thrust Fig. 10), favour preserving q 2 the data on the charged particle multiplicity in gluon jets at LEP Fig. 11, and in jets at the LHC Figs. 12,13 favours preserving the p ⊥ of the branching.
Our preferred choice, in particular in the presence of higher-order matching, is to preserve the q 2 of the branching in order to ensure that the parton shower does not overpopulate the dead-zone. This also ensures a more reasonable value of strong coupling, α CMW S (M Z ) = 0.126 which gives α MS S (M Z ) = 0.118. However given the better description of gluon jets it is reasonable to also consider the alternative of  Table 1 The Monte Carlo parameters obtained for different choices of the cut-off option, the preserved quantity in the shower and weight of the charged particle multiplicity data.  Table 2 The values of χ 2 per degree of freedom obtained in the fit for different choices of the cut-off option, the preserved quantity in the shower and weight of the charged particle multiplicity data. The values are χ 2 as described in the text for the tuning observables, normalised to the sum of the weights for the different bins, and the true χ 2 using the experimental error for the charged particle multiplicities. The number of degrees of freedom for each set of observables is given together with the sum including weights in brackets, where this is different.
preserving the p ⊥ , see for example Fig. 16 from the Appendix.

Conclusions
We have performed a tuning the the Herwig 7 event generator using data on gluon jets from LEP for the first time. To-gether with changes to the non-perturbative modelling this gives a significantly better description of gluon jets, in particular their charge particle multiplicity. It is however impossible to get a good description of the LEP particle spectra and the charged particle multiplicities, particularly in gluon jets, at the same time. We therefore choose the tune p ⊥ -q 2 -B as the default for Herwig 7.1. However for jets at the LHC the tune p ⊥ -p ⊥ -B gives a better description of jet properties.   The data is from [59, 61-63, 73, 79-88] as compiled in [89] While the tunes presented in this paper are an improvement on their predecessors there is a tension between the data on charged particle multiplicities, for both quark and gluon initiated jets, and the data on event shapes and particle spectra from LEP. The cluster hadronization model also continues to have problems describing final states in events with bottom quarks. Any further improvement in the description of this data will require improvements to the nonperturbative modelling.
Acknowledgements This work was supported in part by the European Union as part of the FP7 and H2020 Marie Skłodowska-Curie Initial Training Networks MCnetITN and MCnetITN3 (PITN-GA-2012-315877 and 722104). Daniel Reichelt thanks CERN for the award of a summer studentship during which this work was initiated and acknowl- Appendix A: Generalized angularities and quark and gluon jet discrimination power In this appendix we investigate how the improvements of the simulation of quark and gluon proposed in the manuscript affect the quark and gluon jet discrimination power recently  studied in [9] 9 . For this purpose, we present results for five generalized angularities λ κ β [90]:

LHA width mass
where λ κ β = ∑ i∈jet z κ i θ β i , i runs over the jet constituents, z i ∈ [0, 1] is a momentum fraction, and θ i ∈ [0, 1] is an angle to the jet axis. To quantify discrimination performance, we use 9 The results and the analysis code used for this study is available as a RIVET routine [37], which can be downloaded from https://github. com/gsoyez/lh2015-qg.   Fig. 13 The difference between the average number of particles in central and forward jets compared to data from the ATLAS experiment [10]. classifier separation: where p q (p g ) is the probability distribution for λ in a generated quark jet (gluon jet) sample. ∆ = 0 corresponds to no discrimination power and ∆ = 1 corresponds to perfect discrimination power.
We start with an idealized case of e + e − collisions (see Section 5 of [9] for details). In Fig. 14 we show the discrimination power as a function of an angularity predicted by PYTHIA 8.215 [2], HERWIG++ 2.7.1 [16], SHERPA 2.2.1 [3] , the NNL analytical calculation from [9] and the both p ⊥ -q 2 -B and p ⊥ -p ⊥ -B tunes of Herwig 7.1. Firstly, we see that the both Herwig 7.1 tunes give significantly different results compared to HERWIG++ 2.7.1. In order to understand the source of the difference, in Fig 15 we investigate, for p ⊥ -q 2 -B tune, the following settings variations: -HERWIG: NO g → qq. Turning off g → qq splittings in the parton shower. -HERWIG: NO CR. The variation turns off color reconnections.
We can see that the results are not very sensitive to the change of the settings. This was not the case for HERWIG++ 2.7.1 where the colour reconnection had a huge effect on the discrimination power, see [9]. Therefore, we can conclude that the difference is due to the improvements of the CR model described in Section 3, which as expected reduce effects of CR in the case of e + e − collisions. Secondly, the results of the both Herwig 7.1 tunes are quite similar and closer to the other predictions giving more constrained prediction on the quark/gluon jet discrimination power in e + e − collisions.
In fact just before finishing this paper the new tune was used in [91] confirming that indeed that improvements introduced in the manuscript reduced the tension between Pythia and Herwig and bring Herwig results closer to NNLL' results from [91]. Next, in Fig. 16 we show the results for ∆ in the case of quark/gluon tagging at the LHC (see Section 6 of [9] for details). Here we can see that the differences between HER-WIG++ 2.7.1 and the both Herwig 7.1 tunes are more modest when compared to the previous case of e + e − collisions. However, as expected the largest differences between generators appear for IRC-unsafe observables like multiplicity (0,0) and p D T (2,0), where nonperturbative hadronization plays an important role. It is also worth to notice that the p ⊥ -p ⊥ -B tune which is preferred by the data on the charged particle multiplicity in gluon jets at LEP Fig. 11, and in jets at the LHC Figs. 12,13 gives slightly better discrimination power reducing the gap between predictions of Pythia and the other generators. Finally, it would be interesting to estimate the parton-shower uncertainties [92][93][94][95] in the context of the quark and gluon jet discrimination observables to see whether the remaining discrepancy in the predictions is covered by the uncertainty band.   Fig. 16 Classifier separation ∆ for the five angularities, determined from the various generators at hadron level in the case of quark/gluon tagging at the LHC (see Section 6 of [9] for details). The first two columns correspond to IRC-unsafe distributions (multiplicity and p D T ), while the last three columns are the IRC-safe angularities.