Identification of High-Momentum Top Quarks, Higgs Bosons, and W and Z Bosons Using Boosted Event Shapes

Identification of High-Momentum Top Quarks, Higgs Bosons, and W and Z Bosons Using Boosted Event Shapes


At the Large Hadron Collider, numerous physics processes expected within the standard model and theories beyond it give rise to very high momentum particles decaying to multihadronic final states. Development of algorithms for efficient identification of such “boosted” particles while rejecting the background from multihadron jets from light quarks and gluons can greatly aid in the sensitivity of measurements and new particle searches. This paper presents a new method for identifying boosted high-mass particles using event shapes in Lorentz-boosted reference frames. Variables calculated in these frames for multihadronic jets can then be used as input to a large artificial neural network to discriminate their origin.

I Motivation

The center-of-mass energy of collisions at the Large Hadron Collider (LHC) has now reached 13 TeV in 2015 and 2016 running, up from 7 TeV and 8 TeV in the 2011 and 2012 running periods, respectively. The large LHC general-purpose experimental collaborations ATLAS and CMS, working with theoretical colleagues, have developed powerful new analysis techniques to efficiently identify high-momentum (“highly boosted”) hadronically decaying top quarks, Higgs bosons, and and bosons, while rejecting the background from multijet final states produced in QCD processes. As energies increase, the rate of production of highly boosted heavy particles increases dramatically due to the increased parton-parton luminosity at large invariant masses. Numerous production mechanisms for highly boosted particles exist within the standard model (SM) and are predicted in theories beyond the standard model (BSM).

Figure 1: Momentum vectors of the final state particles from top quark decay to three jets in the laboratory frame and in a boosted frame corresponding to a top quark with the original jet laboratory momentum.

Up to now, the primary method for identifying, for example, decays of boosted top quarks decaying to a b quark and boson, where the W boson decays to a quark pair, has been to use the pattern of energy (or momentum) flow in the two-dimensional space of pseudorapidity () and azimuthal angle () within a “fat” jet reconstructed from observed final state particles (1); (2); (3); (4); (5). By “fat” we mean here a jet reconstructed using a large value of the cone size parameter governing the jet finding algorithm, for example the anti- algorithm (6) with a cone size paramter of 0.8. The algorithms look for substructure within such jets, where there are distinct and separable deposits of energy from the individual jets formed from the initial-state quark and light quark jets from the top decay. Similar algorithms can be used to identify boosted and bosons decaying to quark pairs, and Higgs bosons decaying to .

We present here a new approach in the identification of boosted massive particles. The goal is a generic algorithm which can be used to determine the origin of individual “fat” jets with large transverse energy () produced in collisions at the LHC. In the ideal case such an algorithm will yield estimators allowing analysts to set operating points in a multidimensional output space and optimize the efficiency for identifying hadronically decaying quarks, bosons, and and bosons, discriminating against jets from light quarks and gluons.

Ii Boosted Event Shapes

In this new algorithm, one begins with a list of the fat jet constituents, which are reconstructed from charged particle tracks, calorimeter energy deposits, and information from the outer muon system. For example the CMS experiment utilizes a sophisticated algorithm known as “particle flow” (PF) (7) to arrive at a list of particle candidates from which jets are found using jet clustering algorithms. Typically the PF algorithm can identify 80-90% of the particles in a jet, leaving relatively small corrections to determine the overall jet energy and momentum.

Given the list of fat jet constituents, then, we successively Lorentz-boost all of them into hypothetical reference frames corresponding to our heavy particle () masses, along the jet momentum direction, and assuming the momentum of the heavy particle is equal to the observed jet momentum. In the correct reference frame one would expect to see back-to-back jets for a , , or event, and three jets in a plane for a quark decay. Figure 1 shows the momentum vectors of the final state particles from a top quark decay in the laboratory frame and in a boosted frame corresponding to that of a top quark with the momentum of the jet. As the figure shows, one clearly sees a three-jet structure in the boosted frame, which in the laboratory appears to have a two-jet substructure. This illustrates that this method has the potential to resolve jets which strongly overlap in the laboratory frame momentum space, and thus the detector itself.

Within each reference frame we then calculate event shape parameters and mometum balance estimators. These include Fox-Wolfram moments (8), the eigenvalues of the sphericity tensor (9), and thrust (10). We also run the anti- jet finding algorithm in the boosted frame, relative to the boost axis and using a cone size parameter of 0.5. We then calculate the of the four leading jets, the masses of pairs of the three leading jets, and the invariant mass of the four leading jets. If we have boosted into the correct frame, corresponding to the true origin of the fat jet, then we would expect that the overall momentum of the boosted constituents should be near zero, either in terms of the total momentum or the momentum along the boost direction. Hence we also calculate the magnitude of the total momentum of the four leading jets and the ratio of their longitudinal momentum sum to the total. Generically we call these quantities “boosted event shapes” (BES). This approach is inspired by earlier work at LEP aimed at discriminating decays of the boson using the “boosted sphericity product” (11). It differs from the approach described in (12), and employed by the ATLAS experiment in their measurement of high- vector boson pair prodiction cross section (13), in that we do not assume a rest frame based on the invariant mass of the sum of the observed jet constituents.

To demonstrate the separation that these BES variables can provide, we generate simulated events with the PYTHIA Monte Carlo generator (14) in which we produce a new hypothetical particle which decays to , , or . The mass of the is taken to be 3.0 TeV. This results in a distribution of the transverse energy of the decay prodicts which is relatively uniform up to about 1.5 TeV.

We then simulate the detector response with the PGS program (15), and find jets reconstructed from calorimeter energy deposits using the anti- algorithm with a cone size of 0.8. A crude simulation of PF constituents is then achieved by selecting all generator-level final state particles with GeV matching a jet, and imposing a uniform 90% efficiency. Note that this results in a list of jet constituents which does not reflect detector resolution effects or the effect of additional (“pileup”) interactions.

We take the list of jet constituent four-momenta and boost each one into a reference frame with velocity , where to calculate we use the masses of the top quark, Higgs boson, , and bosons. In each frame we then calculate the above-mentioned event shape quantities. Since we use the four Fox-Wolfram moments, the three sphericity tensor eigenvalues, the thrust, and ten quantities characterizing the jets found in the boosted frame, for a given jet this results in 418=72 different quantities. Space does not permit displaying the distributions of all these quantites, but we show two of them in Figure 2. These distributions are based on fat jets found in events in which a 3.0 TeV is produced and decays to , , or , as implemnented in the PYTHIA Monte Carlo generator. In these figures the histograms show how the distributions evolve as one successively increases the mass used to calculate the boosted frame.

Figure 2: Distributions of the first Fox-Wolfram moment (left) and the ratio of the longitudinal to toal momentum of anti- jets (right) in various boosted reference frames, for top quarks, bosons, and quark jets from 3.0 TeV decays.

Iii Neural Network BES Discriminator

With the large number of BES quantities characterizing a given jet, it is natural to deploy a multivariate technique such as artificial neural networks (ANN) or boosted decision trees (BDT) to arrive at a classification of jets as to their source ( or light quark or gluon jet). As an initial test we created training and test samples using the 3.0 TeV events. From jets matched to generator level , , and , required to decay all-hadron final states, training and test sets are made. For each event we have 50 input variables. From the 72 variables enumerated above we eliminate the boosted jet and invariant masses for all but the frame, since the quantities in the different frames are so highly correlated from one frame to the next as to be redundant, which is understandable given their kinematic nature. This leaves 48 quantities, to which we add two more: the jet and . We arrive at training samples with 20k events each of , , and jets, and 15k events in each of the three test sets. To avoid the multivariate method using the as a discriminating variable, we draw training and test events from the and samples according to the distribution of for jets from top decays.

As a demonstration that this technique can potentially achieve a high degree of accurate classification, we train a set of feedforward multilayer perceptrons with various numbers of hidden layers and nodes in each layer, and two output nodes. The output nodes are trained to be (0,0) for jets, (0,1) for jets, and (1,0) for jets. The ANN is trained using backpropagation with a learning strength parameter of 0.001, a momentum parameter of 0.5, and updating the weights and thresholds after every 50 training events. The learning strength parameter decays exponentially during training with a time constant of 500 epochs (complete presentations of the training sets). Performance of the networks is monitored during training, and we observe stable performance asymptotically, well before the maximum number of 1000 training epochs is reached. We see little variation in performance as a function of the number of hidden layers and number of nodes in each layer - in fact we see more variation from one network to another having the same number of nodes and layers, but starting at a different random weights and thresholds. The results presented here are attained with nets with three hidden layers of 40 nodes each.

Figure 3 illustrates the neural network output distributions for the three jet types. The figure shows that there is a high degree of separation of the three sources. One can also characterize the accuracy of classification using a matrix showing how many true , , and jets are classified as such. Using a simple classification in which the category is taken to be the one nearest the neural network outputs, we see from Table 1 that approximately 85% of jets are classified correctly. One can also see the performance in Figure 4, which show the efficiency of and classification, respectively, versus the efficiency for the other two categories. From these curves one can see, for example, that at an efficiency of 60% for identification, one has a background efficiency for jets of just over 2%.

Figure 3: Distributions of the neural network outputs for test samples of jets originating from all hadronic top quark decay, decay, and quarks for networks trained with charged and neutral constituents (top) and charged conatituents only (bottom).
Classified as 84.2% 3.9% 10.2%
Classified as 5.6% 86.8% 5.9%
Classified as 10.4% 9.3% 83.9%
Table 1: Probabilities for classification of , , and jets as such using the neural network outputs, training on charged and neutral jet constituents.
Figure 4: Efficiency for identification of jets from top decay (horozintal axis) versus efficiency for jets (green) and jets (blue), showing the results using charged and neutral jet constituents (solid curves) and using charged constituents only (dashed curves).

While the above results are encouraging, it is clear that in an actual LHC analysis the performance will likely not attain these levels. Perhaps the main reason is that in real analyses the identification of individual jet constituents will not be as efficient or accurate as assumed in this prototype study. At very high momentum, in particular, when there is more confusion between neutral particles, it may be necessary to rely solely on charged particles. For this reason we tested the performance of the BES algorithm using charged tracks only, selecting only those generator-level final state particles which were charged, having GeV. Figures 3 and 4 compare the results using charged constituents only to the results from using both charged and neutral constituents. Remarkably, the performance for the separation of top quarks from bosons remains nearly the same, but the performance for separation of top from jets is degraded, with background efficiency a factor of two larger at a top quark efficiency of 60%.

It is also of interest to examine whether this new algorithm is capable of distinguishing all-hadronic top quark decays, Higgs boson decays to , , and when the decaying particles are higly boosted. To test this we also used samples of a 3 TeV mass decaying to , , , and . As before the events are simulated with PYTHIA and passed through the PGS simulation. Fat jets (anti-, ) are then matched to the decaying boosted objects and 90% of the jet constituents with GeV and lying within of the jet direction are recorded for the BES analysis. We consider here the performance using both charged and neutral objects.

We train a 6-layer neural network with the same 50 inputs, 4 layers of 40 hidden nodes each, and an output layer with four outputs. The network is trained to give values of (1,0,0,0) for jets, (0,1,0,0) for jets, (0,,0,1,0) for , (0,0,0,1) for jets, and (0,0,0,0) for quark jets. The resulting output, which are points in a 4-dimensional space, can nevertheless be usefully displayed on the two-dimensional plane of the difference between the first and third outputs ( axis) and the difference between the second and fourth outputs ( axis). Figure 5 shows these distributions for the five jet categories, and illustrates the great potential to discriminate them well, even from decays. As an example, Figure 6 shows the efficiency for identifying top quark decays as a function of the efficiency for identifying the other particle types.

Figure 5: Distributions of the neural network output differences for a 4-output neural network trained to discriminate jets originating from all-hadronic top quark decay, decay, decay, decay, and quarks for networks trained with charged and neutral constituents.
Figure 6: Efficiency for identification of jets from top decay (horozintal axis) versus efficency (vertical axis) for jets originating from all-hadronic decay of , , and bosons, and quarks, for a 4-output neural network based on BES variables from charged and neutral jet constituents.

The performance of the algorithm may benefit from the various jet “pruning” techniques such as the soft-drop method (19) and PUPPI (16), but this is best studied in the context of a full simulation including all detector effects and the effets of multiple interactions.

The theoretical uncertainties on the modeling of rest-frame quantities were studied in detail elsewhere (17). We note that these can be at least partly overcome by training the neural network using background (light quark and gluon) samples derived from observed dijet events, where the other jet is anti-tagged as coming from .

Finally, we emphasize that the variables arising from the boosted event shapes can be readily combined with variables used in other approaches, including the presently used  (18) and soft-drop mass, and other novel approaches such as those based on two-dimensional imaging of jets (20), to arrive at an even more selective algorithm. One can expect, hoever, that the the BES variables will retain their discriminating power up to ultrahigh boosts at which 2D variables become less discriminating.

Iv Summary

We have described here a new approach to identifying very high momentum hadronically decaying top quarks, Higgs bosons, and vector bosons, based on boosting the individual constituents of jets sequentially into rest frames corresponding to the hypothetical particle masses and calcualting in each of those frames a set of event shape quantities. These quantities include the Fox-Wolfram moments, sphericity tensor and its eigenvalues, thrust, and quantities derived from anti- jets calculated using the boost axis as the “beam”. The algorithm is capable of separating jets which strongly overlap in laboratory momentum space. In an initial, idealized scenario using generator-level particles the algorithm shows excellent separation of hadronic jets arising from the decay of top quarks, Higgs bosons, and bosons and from quarks. This level of separation is largely maintained when using charged jet constituents only, as might be necessary for ultra-high-energy jets. The variables arising in this algorithm can be readily combined with those from other boosted particle tagging techniques.


  1. (CMS Collaboration) Top Tagging with New Approaches, CMS-PAS-JME-15-002 (2016).
  2. G. Aad, et al. Identification of high transverse momentum top quarks in collisions at = 8 TeV with the ATLAS detector (ATLAS Colalboration), CERN-EP-2016-010 (2016).
  3. G. Aad, et al. (ATLAS Collaboration), JHEP 1309:076 (2013).
  4. D. E. Kaplan, et al., Phys. Rev. Lett. 101, 142001 (2008).
  5. S. D. Ellis, C. K. Vermilion, and J. R. Walsh, Phys. Rev. D81 094023 (2010).
  6. M. Cacciari, G. P. Salam, and G. Soyez, JHEP 0804:063 (2008).
  7. (CMS Collaboration) Particle-Flow Event Reconstruction in CMS and Performance for Jets, Taus and Missing Transverse Energy. (CMS-PAS-PFT-09-001), 2009.
  8. G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978).
  9. J. D. Bjorken and S. J. Brodsky, Phys. Rev. D1, 1416 (1970).
  10. S. Brandt, Ch. Peyrou, R. Sosnowski and A. Wroblewski, Phys. Lett. 12, 57 (1964); E. Fahri, Phys. Rev. Lett. 39, 1587 (1977).
  11. L. Bellantoni et al., Nucl. Instrum. Meth. A310, 618 (1991).
  12. C. Chen, Physical Review D 85, 034007 (2012).
  13. G. Aad, et al. (ATLAS Collaboration), New J.Phys., 16(11):113013 (2014).
  14. T. Sjostrand, S. Mrenna, and P. Skands, JHEP 0605:026 (2006).
  15. M. Carena, et al., FERMILAB-CONF-00-279-T (2000).
  16. D. Bertolini, et al., JHEP 1410:59 (2014).
  17. I. Feige, et al., Phys. Rev. Lett. 109, 092001 (2012).
  18. J. Thaler and K. Van Tilburg, JHEP 1103:015 (2011).
  19. A. J. Larkowski, et al., JHEP 1405:146 (2014).
  20. P. Baldi, et al., Phys. Rev. D93, 094034 (2016).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description