Jet flavour tagging for future colliders with fast simulation

Jet flavour identification algorithms are of paramount importance to maximise the physics potential of future collider experiments. This work describes a novel set of tools allowing for a realistic simulation and reconstruction of particle level observables that are necessary ingredients to jet flavour identification. An algorithm for reconstructing the track parameters and covariance matrix of charged particles for an arbitrary tracking sub-detector geometries has been developed. Additional modules allowing for particle identification using time-of-flight and ionizing energy loss information have been implemented. A jet flavour identification algorithm based on a graph neural network architecture and exploiting all available particle level information has been developed. The impact of different detector design assumptions on the flavour tagging performance is assessed using the FCC-ee IDEA detector prototype.


Introduction
Precision measurements of standard model (SM) parameters are key objectives of the physics program of future lepton and hadron machines [1][2][3][4][5][6].In particular, the measurement of the Higgs couplings to bottom (b) and charm (c) quarks, and gluons (g) [7][8][9][10][11][12][13], the Higgs self-coupling [14] and the precise characterisation of top quark properties, such as the top quark mass [15] and its electroweak couplings [16,17] require an efficient reconstruction and identification of hadronic final states.Being able to efficiently identify the flavour of the parton that initiated the formation of a jet, known as jet flavour tagging, is therefore critical for the success of the physics program of future electroweak factories [18].The large statistics of hadronic Z boson decays at both lepton circular machines and future hadron machines would provide copious control samples to calibrate jet tagging algorithms in data.
Jets originating from b and c quark decays contain a b or c hadron that typically travels a macroscopic distance before decaying into lighter hadrons.Compared to b, c, up or down jets (collectively referred to as ud or light jets in what follows), strange quark (s) jets contain a larger fraction of s mesons and baryons.Gluons (g) carry a larger colour charge than quarks and thus tend to produce jets with a large particle multiplicity.Quarks have a harder fragmentation function compared to g, which results in a larger fraction of the jet momentum carried by a smaller fraction of the constituents.Jet flavour tagging algorithms aim at identifying these characteristic end-products of the fragmentation and hadronization of the initial parton.
The first b and c quark tagging algorithms were developed at LEP [19,20] and the Tevatron [21,22].These algorithms typically rely on the detector capability to identify and measure charged tracks with a significant displacement (cτ ∼ 500 (150) µm) from the beam axis originated from long-lived B (D) meson weak decays.On the other hand, tracks from the charged hadrons produced in ud quark decays feature a small distance at closest approach to the interaction point.Therefore, in the case of b and c tagging, tracks are typically clustered to reconstruct possible secondary vertices (SVs).However, c tagging is more challenging than b tagging, due to its properties laying between the b and ud or g jets.The track multiplicity and the mass of the SV (expected to be large for heavy flavour jets), together with the presence of a non-isolated electron or muon indicating a semi-leptonic heavy flavour decay, are also used as discriminating variables in traditional heavy quark tagging algorithms.Such taggers, widely used in the early days of current LHC experiments [23][24][25][26] are implemented by directly applying a selection on a combination of the tracks and SV properties, by constructing a likelihood ratio or a multi-variate discriminant based on a set of jet-level properties.
Recently, a new generation of advanced machine learning based jet tagging algorithms has been developed [27][28][29][30], bringing more than an order of magnitude improvement in background rejection compared to the traditional approaches in heavy flavour and g tagging.Three are the primary reasons for this success.First, significant advancements in the architecture of the neural networks used, as well as new jet representations that allow to better capture the jet properties have been achieved.Second, these algorithms exploit directly low-level information, e.g., from reconstructed particles (as in the Particle-Flow algorithm [31]) or even reconstructed hits, compared to traditional methods.This allows to explore in much more depth the true potential of the detectors and the event reconstruction, and also better capture the jet properties compared to algorithms relying on jet-level observables.Moreover, the nature of each of the jet constituents, via particle identification techniques (PID), is expected to provide an additional useful handle in discriminating between different jet species.Powerful particle identification capabilities based on ionisation energy loss (via dE/dx or cluster counting), or via precise time-of-flight measurements, are expected to be highly beneficial for jet flavour tagging, in particular for s tagging where the identification of charged kaons is crucial.Finally, the developments in computing, e.g., graphics processing units, and the availability of very large Monte Carlo simulated and collision data samples, were critical for the development of these advanced methods.
In this paper we present a general framework for building a jet flavour tagging at future colliders using fast detector simulation and state-of-the art machine learning techniques.A major goal of the present work has been to allow for the evaluation of the impact of specific detector design options on the jet flavour tagging performance (and in turn on the physics potential) in an efficient yet precise way.To this end, we have implemented two key additions to the official Delphes fast simulation framework.The TrackCovariance [32], described in Section 2, which allows for a simple definition of a tracker geometry, and the fast simulation and the reconstruction of the parameters and covariance matrix of charged particles tracks.The TimeOfFlight [33] and ClusterCounting [34], described in Section 3, open up the possibility to model particle identification in Delphes.Section 4 describes the input observables and the implementation of the jet flavour identification algorithm.The tagging algorithm is based on ParticleNet [35], using state-of-the-art jet representation and a graph neural network (GNN) architecture.The performance of the algorithm is evaluated using one of the FCC-ee/CEPC baseline detector concepts, the IDEA [1,36,37] detector.Variations around the baseline using Higgs decays taken from a Higgsstrahlung sample at √ s = 240 GeV are discussed.Finally, a discussion of the results, together with limitations of the current approach and perspectives for future work are presented in Section 5.

Fast Tracking simulation
The tracking system is a major part of modern detectors for high energy physics experiments and arguably the most relevant for jet flavour tagging since it is responsible for reconstructing and identifying charged particles.The design of this system, its optimization, and the evaluation of its performance on many specific physics benchmarks is a fundamental step in the planning of future experiments.To this end, we have developed and included in Delphes, a versatile and modular framework to easily study different detector configurations, and provide for each of those a fast simulation of the tracking performance.The corresponding module is named TrackCovariance [32].In this section we present the general implementation of the algorithm, while technical details on the speed optimisation and randomisation can be found in Appendix A.
While various attempts to calculate the track resolution analytically have been made (see for instance [38]), they usually make highly simplifying assumptions such as equal spacing and equal detector resolution, that make them unsuitable to use for a realistic combined tracking system.The tracking system geometry is described in terms of layers.Only two types of layer geometries are considered: cylinders coaxial with the beam axis and planar disks orthogonal to the beam axis (z-direction).Each layer can be either associated to a measurement with a given resolution or else be just included to describe passive material in the system.An accurate description of the material inside the tracking volume is important to estimate appropriately the contribution of multiple coulomb scattering to the track resolution.Several measurement geometries are allowed: axial or stereo strips and wires, and pixels.
The tracking system is located inside the solenoid magnet generating a constant field, B, directed parallel to the z-direction.With these assumptions, charged tracks follow a helix trajectory that is described with a set of five parameters: α = (D, ϕ 0 , C, z 0 , λ).These parameters are defined in the point of closest approach (PCA) of the track to the z-axis; D is the signed transverse distance of the PCA from the z-axis, ϕ 0 is the track azimutal angle, C is the signed half curvature, z 0 the z-coordinate of the PCA and λ the cotangent of the track polar angle.Given a charged particle originating at x with momentum p and charge Q, the parameters α are uniquely defined, as is the associated trajectory.
During its motion the charged particle will cross some of the layers described in the geometry.At each crossing the particle will undergo small random changes of its direction due to multiple scattering and, in the case of measurement layers, a generalized coordinate, d * , will be measured with an uncertainty given by the specific detector resolution.The track parameters are reconstructed from the measured coordinates by minimizing the following χ 2 with respect to the track parameters α.The χ 2 is defined as: where d * is the array of measured coordinates and d that of predicted coordinates, that can be computed from the track parameters α and the geometry of each measurement layer.S is the covariance matrix of all the measurements and includes contributions from the detector resolution and from the multiple scattering.The superscript t indicates the transpose of a vector or a matrix.Assuming d 0 to be the array of predicted coordinates for which the χ 2 is minimized, then for small variations of the track parameters relative to this minimum, δ α, we have: where A is the derivative matrix.Including Equation .(2.2) in Equation .(2.1) we obtain: Differentiating with respect to the track parameters we obtain the track parameter covariance matrix, C: This equation highlights the key ingredients to estimate the track covariance matrix; the derivative matrix and the covariance matrix of the measurements.The former is straightforward and can be derived for every type of measurement given the track equation.The latter requires the combination of two elements: the intrinsic detector resolution and the multiple scattering contribution, as shown in the following equation: where the indices i and j identify the measurement layers, σ i is the detector resolution for layer i and M ij is the multiple scattering contribution.The M ij includes contributions from all scattering layers below the smallest of the two indices, as shown in the following equation: where L i is the distance traveled by the track to the layer i and θ k (i, j) the standard deviation of the multiple scattering angle generated by layer k after correcting for projection factors specific for layers i and j.
Once A and S have been determined for a given track, the parameter covariance matrix can be computed analytically using Equation.(2.4).The obtained resolution as a function of momentum of the track parameters (p T , D, z 0 , θ) for two different reference detector configurations proposed for a future e + e − collider is shown in Figure 1.

Particle identification
Particle identification techniques can play a major role in the identification of the jet flavour.In particular, as will be discussed in Section 4, s jets contain a significant fraction of charged kaons (K ± ) compared to u or d jets that are mostly composed of charged pions (π ± ).Given that the performance of such algorithms heavily depends on explicit detector design choices, it is crucial to be able to first simulate appropriately the detector response and then to implement such particle identification algorithms.
Two complementary particle identification techniques have been included in the Delphes fast-simulation.The precise measurement of the time of arrival of tracks in the outermost part of the tracking volume, together with the momentum and the path length, provide an indirect measurement of the particle mass via the well known time-of-flight method.This method has been implemented in the TimeOfFlight module [33].The cluster counting method, dN/dx, implemented in the ClusterCounting module [34], consists in counting the multiplicity of the primary ionization clusters produced along the track in gaseous detectors, which together with the particle momentum, can also be used to infer the particle mass.In this section we discuss the implementation of these two methods within the simulation framework.

Time-of-flight
The time-of-flight (t flight ) of a particle can be expressed as: where t F is the measured time after propagation, t V is the particle time of production at vertex, L is total path length, and p, E and m are the momentum, energy and mass of the particle, respectively.Provided that the quantities L and p (or E) and t V can be measured, the measurement of t flight provides an estimate for the particle mass and thus a powerful handle for particle identification.
For charged particles the reconstructed mass is given by: The initial position (and therefore L) and the particle momentum p are reconstructed by means of the inner/outer tracking system, and simulated with the procedure described in Section 2. The time of a particle production at vertex t V can be estimated indirectly, with the following procedure.Assuming that the beamspot has a small time (σ B,t ) and longitudinal (σ B,z ) spread compared to the precision of the timing measurement device, the time of the primary vertex can be simply taken as t PV = 0.However, if the particle originates from a highly displaced vertex (e.g. from K S or Λ.), assuming t V = 0 can lead to a severe over-estimate of t flight .A more accurate estimate for the vertex time corresponds to t V = r V β V , where r V is the distance of the vertex to the origin and β V is the vertex velocity, computed from its outgoing particles.In the current study we assume we are able to reconstruct the initial time at vertex perfectly and therefore we take the initial time from Monte Carlo simulation.The time of flight distribution of charged Kaons and Pions emitted at 90 • is shown in Figure 2 (a), assuming a 30 ps timing resolution, which allows for an efficient K/π discrimination for momenta p < 3.5 GeV.
For neutral particles the mass can be reconstructed from the energy measurement provided by the calorimeters: At low momenta, where the time-of-flight method is expected to provide good identification capabilities, the calorimetric energy measurement is sub-optimal and leads to poor m t.o.f.resolution for neutral particles compared to charged particles.Moreover, the vertex time determination is inaccessible for neutral particles.The assumption t PV = 0 for all neutral particles leads to an additional uncertainty on the m t.o.f.estimate.The reconstructed m t.o.f. for K L and neutrons is shown in Figure 2 (b).

Cluster counting
The cluster counting technique is expected to provide improved particle identification relative to the more commonly used dE/dx methods in large drift chambers or TPCs [39,40].In addition it does not require the tuning of truncated mean algorithms to suppress the large Landau tails present in the dE/dx distribution.The number of ionization clusters per unit length is obtained from a very detailed simulation program, Heed++ [41], now fully integrated into Garfield++ [42].An array of number of ionization clusters per unit length for several values of βγ is obtained from Garfield and used to interpolate the average cluster density.The total mean number of clusters is found by multiplying for the track length in the chamber.Finally the observed cluster number is obtained by extraction over a Poisson distribution with that mean.Four common gas options are available: pure Helium or Argon, He 90% + Isobutane 10%, Argon 50% + Ethane 50%.This library can be easily extended if needed to a larger collection of gas mixtures.
In Figure 3 (a) the potential for K/π separation is shown for a He 90% + Isobutane 10% mixture over a wide range of momenta (2 < p < 30 GeV) .The combination of the cluster counting and time-of-flight techniques is displayed in Figure 3 (b) and shows an efficient separation of K ± / π ± separation (≥ 3σ) for momenta p < 30 GeV.

Jet flavour identification
In this Section a novel jet tagging algorithm is presented.The jet flavour discrimination uses reconstructed observables at the level of the jet constituents.For simplicity, the jet flavour discriminant is built and evaluated using e + e − collisions reconstructed with the IDEA detector concept and will thus be referred as ParticleNetIdea.While the obtained performance is specific to the clean e + e − environment and the explicit detector specifications, the inputs and the construction of the discriminant itself are general.We first discuss the event generation and reconstruction details, then introduce the particle- level input observables and the architecture of the neural network discriminant.Finally we address the tagger performance and its robustness with respect to different detector choices.

Simulated data
The simulated sample consists of e + e − → ZH events produced at a center of mass energy √ s = 240 GeV.The Higgs bosons decay to H → gg or H → q q, where q = (u, d), s, c, b with relative fraction as expected for a SM Higgs boson with m = 125 GeV, whereas the Z bosons always decay to a pair of neutrinos.The hard scattering process is generated with MadGraph5_aMC@NLO [43], while Pythia8 [44] is used for modeling the decay, parton-shower and hadronisation processes.Five different samples, corresponding to each jet flavour category ( ud, s, c, b, g ) containing 10 6 events each (or equivalently 2 × 10 6 jets) are used for the training.Final state particles are reconstructed with the Delphes PF algorithm.In particular, charged particles are reconstructed using the latest TrackCovariance module described in Section 2, and the time-of-flight and number of ionisation clusters per unit length (dN/dx), are reconstructed using the TimeOfFlight and ClusterCounting modules, described in Sections 3.1 and 3.2, respectively.Neutral particles (photons and neutral hadrons) are reconstructed by the PF algorithm implemented in the DualReadoutCalorimeter module [45].The time-of-flight (and corresponding reconstructed mass m t.o.f. ) of neutral hadrons is also included and assumes a 100 ps resolution, as opposed to 30 ps assumed for charged particles.The baseline simulation setup assumes the nominal IDEA detector concept [36,37].Jets are clustered with the FastJet-3.3.4 [46] package using the e + e − generalized k T algorithm [47,48] with parameter p = −1.

Input features
The jet constituents in the form of PF candidates are used as inputs to the ParticleNetIdea algorithm.For each PF candidate we define a set of input observables (features) that are summarized in Table 1.The first set of inputs, denoted as kinematics, uses features derived from the 4-momentum of each jet constituent.These include the energy measurement of the constituent relative to the jet energy and the direction of the jet constituents relative to the jet momentum.The second set of features, labelled as displacement, includes observables related to the longitudinal and transverse displacement of the jet constituents which are more relevant to identify jets originating from the hadronization of the b and c quarks.Finally, the third set of inputs, labelled as identification, refers to the nature of each particle using the PF reconstruction and the particle identification (PID) algorithms presented in Section 3. mass calculated from time-of-flight dN/dx number of primary ionisation clusters along track isMuon if the particle is identified as a muon isElectron if the particle is identified as an electron isPhoton if the particle is identified as a photon isChargedHadron if the particle is identified as a charged hadron isNeutralHadron if the particle is identified as a neutral hadron The total number of reconstructed jet constituents, shown in Figure 4(a), is typically larger for g jets compared to quark jets due to their different color factor.We note that the particle multiplicity is shown here for illustrative purposes only as it is not used directly as input to ParticleNetIdea since it is a jet-based variable, while only particle-level observables are used.The remaining distributions of Figure 4 correspond to particle-level observables and are calculated using the charged constituent with the largest displacement.Figure 4 (b) displays the relative energy of the jet constituent with respect to the jet energy.Gluon jets populate lower values of this observables, indicating that the jet energy is more democratically distributed among the constituents.Figures 4 (b) display observables relevant for b and c quark identification, such as SIP 2D (left) and its significance SIP 2D /σ 2D (right) as defined in Table 1.As expected, in b jets, and to smaller extent in c jets, a significantly larger displacement is observed compared to the other jet flavours.Displaced particles can also be present in other jet flavours, e.g. from long-lived K 0 S or Λ hadrons decays, but represent a much smaller fraction.

The flavour tagging algorithm
The ParticleNetIdea algorithm is based on the ParticleNet jet tagging algorithm [35].
ParticleNet uses an advanced network architecture, based on Graph Neural Networks (GNN) that first developed in the context of proton-proton collisions at the LHC.A novel jet representation was utilized in ParticleNet, where jets are represented as an un-ordered set of particles.As shown in Refs.[35,[49][50][51][52][53][54][55][56][57][58], this provides a more natural jet representation compared to alternative approaches based on jet images [59][60][61][62][63][64][65][66][67][68] or ordered lists of jet constituents [69][70][71][72][73][74][75][76][77] and translates to an improved tagging performance.A hierarchical learning approach using convolution operations [78] is adopted.Different convolutional layers are used to learn features at different scales: the shallower layers explore local neighborhood information, whereas more global structures are learned by deeper layers.The jet constituents are represented as a graph, where each node of the graph is a jet constituent, and relationships between the particles are the edges of the graph.Each node has a set of features related to constituent properties.However, the graph is not static, rather it is updated after each convolutional operation.The ultimate goal is to group jet constituents according to their proximity in the multi-dimensional space defined by the learned features.
The current ParticleNetIdea implementation uses up to 75 constituents for each jet, sorted by the highest momentum, which typically correspond to more than 99% of the total jet momentum.The algorithm is designed to discriminate between five orthogonal jet classes: ud, s, c, b, and g jets.The training is performed using the weaver package [79] on 10M jets (2M per category) over 30 epochs on a NVIDIA GTX 1080Ti GPUs.The network outputs 5 real numbers D i (i = g, (ud), s, c, b) between 0 and 1 (discriminants), one for each jet category.Approximately 1M jets are used to evaluate the ParticleNetIdea performance.For every jet flavour pair (i, j), the binary discriminant is constructed as: where D i(j) are the output scores of the classes i and j.For example, D b,c represents the binary discriminant for tagging b quark jets against c quark jets.The efficiency of tagging flavour i as function of the probability of mis-identifiying the jet as flavour j (mistag rate) can be constructed by computing the probability of selecting jets that satisfy D i,j > α, for α ∈ [0, 1].The receiver operating characteristic (ROC) curve, i.e. the mistag rate as a function of the tagging efficiency (for every α), is used as a figure of merit for evaluating the tagger performance for every jet flavour.

Results
The nominal ParticleNetIdea flavour tagging performance is shown in Figure 5 for different jet flavours.The b tagging performance is shown in Figure 5(a).The most effective discrimination is observed against ud jets since these contain mostly tracks with no displacement.For high b tagging efficiency, g jet rejection is more effective than c jet (with both being less effective than u, d or s jet rejection).Conversely at small tagging efficiencies (i.e. for high tagging purity), c jet rejection becomes more effective than g jet rejection due to a sizeable probability for g to produce b b splittings.For c tagging, at high efficiencies, b jet discrimination is the most effective (due to a large difference of lifetime between B and D mesons), followed by ud and g jet rejection.For large c tagging purity (i.e. at low efficiency and high background rejection), we observe that b jet rejection becomes more challenging than ud jet rejection, which is expected since a fraction of B meson have inherently a comparable decay length to D mesons.We also observe that in this regime g → cc splittings result into more challenging g rejection.In Figure 5(c) the s tagging performance is shown.The most effective discrimination is observed against b jets followed by c jets due to the large displacement of their tracks.The mistag rate against g and ud jets is substantially larger since displacement observables are not discriminating, and the algorithm relies mainly on PID-related variables.Rejection of g jet is more effective than ud jet one since s and ud jets have similar particle multiplicities.Finally, the g tagging performance is displayed in Figure 5(d).Rejection of ud jets is the most challenging, due to similar particle displacement and nature, followed by s, c and b jet rejection.
The modularity of the framework enables the study of the algorithm performance for different detector design choices.In this work, we report two representative examples of possible detector design variations.Figure 6 (a) shows the importance of particle identification information in discriminating s jets from other jet flavours.Exploiting PID information with the nominal dN/dx and t flight resolutions yields to approximately an order of magnitude reduced ud jet mistag rate for the same s tagging efficiency.A timing detector providing an improved t flight resolution of 3 ps for charged particles, yields a small, but detectable improvement compared to the more realistic scenario of 30 ps.The performance obtained using MC truth information for PID ("ideal PID") is also shown for reference.In that case only a marginal improvement in performance is observed, suggesting that the existing detector configurations and PID algorithms are very close to optimal.The usage of PID information brings only modest improvement in other jet flavour tasks.
The distance of the first vertex detector layer to the interaction point is the most important parameter for achieving optimal transverse impact parameter resolution and hence b and c tagging performance.While the nominal IDEA vertex detector provides already an excellent resolution (three layers, with the innermost layer located at 1.5 cm), we study the impact of introducing an additional fourth layer in the pixel detector, closer to the beam pipe, located at 1 cm from the interaction point, on c jet identification.The corresponding performance is displayed in Figure 6 (b).The largest improvement is observed in the discrimination against ud jets, where for the same S , B is reduced by almost a factor of two.Smaller, yet important improvement, is seen in the discrimination against other jet flavours.The impact of an additional pixel layer was studied for other jet flavours treated as signal without significant improvement in the performance.

Conclusion and perspectives
Jet flavour tagging will be a crucial tool for maximising the physics potential at future colliders.This work builds on the design of a fast detector simulation framework, and provides an efficient way to study the impact of different detector design options to the jet flavour tagging problem.A fast tracking module was developed, which allows to easily configure a full tracking geometry including material effects and compute both the charge particle track parameters and the track covariance matrix.Two algorithms that allow for particle identification, the time-of-flight and cluster counting with respectively configurable time resolution and gas composition have also been added.The framework is designed to provide flexibility for further studies, such as the exploration of alternative clustering algorithms, beam energies and final states.
Deep learning techniques based on GNNs have proven very effective for classification problems such as jet flavour tagging and boosted jet tagging at the LHC, and have not been explored yet in the context of future experiments.This paper presents the first algorithm for jet flavour tagging at future e + e − colliders using state-of-the-art jet representation and a GNN architecture.At such future machines where statistics for Higgs processes are moderate, flavour taggers will be required to perform well in the high tagging efficiency regime while still providing excellent background rejection.In this study we have investigated the impact of MIP timing resolution and of an additional inner tracking layer on the tagging performance.More studies are possible and should be pursued: the interplay of MIP and calorimeter timing on PID performance, the impact of the tracker design on displaced tracks performance, K S and Λ reconstruction and hence on s tagging, secondary vertex reconstruction on b, c and s tagging.Another area for future studies is the calibration of the algorithm.The algorithm is designed to have very little dependence on the jet kinematics and therefore a calibration strategy relying on a Z boson sample of unprecedented statistical power expected to be obtained in e + e − experiments seems a promising avenue.
We stress that this study has possible limitations given the inherent optimistic nature of fast simulation.In particular, this tracking simulation include a simplistic particlematter description where multiple scattering is taken into account to derive track parameter resolutions but no secondary emissions are simulated (i.e.electron bremsstrahlung and hadronic interactions are neglected).A natural next step is to assess the limitation of the fast detector simulation framework by validating the results with events produced using Full Simulation.Nevertheless, the set of tools presented in this article should provide robust means for assessing an upper limit of the achievable tagging performance and the relative performance of alternative detector design choices at future e + e − colliders.We also point out that the presented framework should allow for similar optimisations at any future machine, including high energy proton-proton or Muon colliders, acknowledging however that further caution is required due to the lack of simulation of larger background levels.

Figure 1 .
Figure 1.Track parameter resolution for the IDEA and CLD detector concepts for FCC-ee [1].The dashed lines in the top left plot show the multiple scattering contribution.

Figure 2 .
Figure 2. (a) Time-of-flight for K ± and π ± track at θ = 90 • as a function of momentum in the IDEA detector drift chamber.(b) Reconstructed m t.o.f. for K ± , π ± , K L , protons and neutrons.

Figure 3 .
Figure 3. (a) Number of cluster distribution of charged pions and kaons for 90 • tracks in the IDEA detector drift chamber as function of momentum; (b) K/π separation in number of σ as a function of the particle momentum using the dN/dx and time-of-flight methods.

Figure 4 .
Figure 4. Shape comparison of a set of representative observables relevant for jet flavour identification.The different colors correspond to different jet flavours.The FCC-ee IDEA detector concept is used.

Figure 5 .
Figure 5. Evaluation of ParticleNetIdea performance in terms of a receiver operating characteristic (ROC) curve for the identification of different jet flavours i.e., b quarks (upper left), c quarks (upper right), s (lower left), and g (lower right).The different jet flavours considered background are indicated on the labels.The IDEA detector configuration is used.

Figure 6 .
Figure 6.Evaluation of ParticleNetIdea performance of the jet flavour identification for various detector assumptions.(a) Impact of particle identification on s tagging performance.(b) Impact of inner track geometry 3 vs 4 layers on c tagging.

Table 1 .
Set of input variables