Top anti-top pairs at the LHC heavy ion collision: a new interesting probe of quark gluon plasma
We investigate top anti-top quark pair production in lead-lead collisions at the Large Hadron Collider with nucleon-nucleon center of mass energy of 5.5 TeV. Due to the very high temperature and energy density created in heavy ion collision, a new state of QCD matter known as Quark-Gluon Plasma (QGP) is expected to be produced. Top decay products loose energy inside the QGP medium. Therefore, we also study the medium modifications of different kinematic distributions. We observe significant modification in the dijets and trijets invariant mass distributions.We also found that the peak position and shape of the distributions could be used to characterize the nature of jet energy loss in the QGP.
The primary goal of heavy ion collisions (HIC) at Relativistic Heavy Ion Collider (RHIC) at BNL, and the Large Hadron Collider (LHC) at CERN is to produce and study the properties of a hot/dense state of QCD matter known as Quark-Gluon Plasma (QGP) dk1 (). QGP is a deconfined state of matter where quarks and gluons are the effective degrees of freedom rather than nucleons or hadrons dk2 (). Fast partons propagating in a hot/dense nuclear medium are expected to loose a large fraction of their energy jq (). The observation of the suppression of energetic partons in the QGP, that is jet quenching jq_RHIC (), and centrality-dependent dijet asymmetry da_LHC () are the most important results from the HIC at RHIC and LHC experiment, respectively.
-bosons are massive weakly interacting Standard Model (SM) particles. The larger LHC HIC energies open the possibility to probe the nucleus-nucleus collisions via the -bosons. The vector bosons are produced early ( fm/c) and their decay time is small ( fm/c and fm/c pdg ()). Whereas, in the most accepted picture of the QGP formation and evolution, at the LHC, QGP is expected to form after fm/c of the initial hard scattering, thermalize quickly and it might last fm/c. Therefore, weak bosons are produced and decay before the formation of QGP and the decay products of pass through the QGP. -boson dominantly decays to a pair of quarks. Quarks loose energy in the QGP and thus the hadronic decays of -boson could be an interesting probe to characterize QGP. However, in presence of huge QCD dijet background, it is extremely challenging to study the hadronic decays of -bosons. The leptonic BF of -boson is small however, due to small background, the signature could be easily detected at the LHC HIC. Unlike jets, leptons interact electromagnetically with the QGP and loose experimentally insignificant amount of energy within the QGP ConesadelValle:2009vp (). Therefore, the leptonic decays of weak bosons behave as a medium blind reference. ATLAS collaboration has already measured the -boson yield in the leptonic decay channels for TeV ATLAS_WZ ().
A new regime of heavy ion physics will be reached at the LHC with TeV where hard and semi-hard particle production can dominate over the underlying soft events. The higher LHC energies open the possibility to study top quarks for the first time at the HIC. The top quark was discovered at the Tevatron experiment top (). After the discovery, properties of top quark have been extensively studied at the Tevatron and LHC. Presently, different properties of top quark is known with good precision. As an example, top quark mass and full decay width are determined to be GeV and GeV pdg (). Different decay channels and branching ratios of top quark have also been observed and measured. As a result, top quark can now be considered as standard benchmark for other experimental observations. In view of this fact, it is important to investigate top quarks at the LHC HIC and study the medium (QGP) influence on different kinematic distributions which are precisely known from the previous and collider experiments. In this letter, we have for the first time studied production at the LHC HIC and proposed few kinematic distributions for the study of the QGP created in the LHC HIC.
The main source of top quarks at the LHC HIC is the top anti-top () pair production. At leading order in perturbation theory there are two processes that contribute to production: quark-antiquark annihilation, and gluon-gluon fusion, . With TeV center-of-mass energy per nucleon, the NLO+NNLL production cross-section per nucleon-nucleon collision in Pb-Pb reaction is given by pb top_cross (). Therefore, the total production cross-section in minimum bias Pb-Pb scattering is estimated to be 3.5 in the frame work of the Glauber model GM (). The instantaneous luminosity of the LHC Pb-Pb collision at TeV is expected to be . Therefore, nb integrated luminosity data will be accumulated within one month ( second) of Pb-Pb collision at TeV. These estimations indicate that about 3500 events are expected to be produced in one month of LHC HIC running at TeV.
Due to the large top mass, pairs are produced early in the HIC. The decay width of top quark is large and thus the decay time is small. As a result, top quarks are produced and decay before formation of QGP. Top quark decays to a bottom quark () and -boson with almost 100% branching fraction: . -boson subsequently decays hadronically () with 67.7% BF or leptonically () with % BF pdg (). Therefore, for production, there are only three possible final state topologies: (i) Hadronically decaying pairs: Both the top quarks decay hadronically () and give rise to 2- jets and four light quarks jets in the final state with 46% effective branching ratio. (ii) Semi-leptonically decaying pairs: One top quark decays hadronically and the other decays leptonically (). The final state is characterized by - jets+-light quark jets+one charged lepton + one neutrino. The effective branching fraction of semi-leptonic decay mode of pairs is %. (iii) Leptonically decaying pairs: In this case, both the top quarks decay leptonically giving rise to - jets + -lepton + -neutrino final state with 10.4% effective branching fraction.
The jets (-jets as well as light jets) and leptons can be observed at the LHC. The neutrinos remain invisible at the detectors and give rise to a imbalance in the transverse momentum known as missing transverse momentum. However, in the HIC environment, faithful measurement of missing transverse momentum will be challenging in the presence of a continuum energy deposit from the produced QGP.
Hadronically decaying events suffer from huge QCD background. Moreover, due to a large combinatorial background it is difficult to reconstruct top quark and -boson from hadronically decaying pairs. Leptonically decaying channel is a clean signal channel due to less background. However, the rate of this channel is suppressed by the top quark leptonic branching fraction. Moreover, in the absence of proper knowledge about missing transverse momentum, the reconstruction of top quark and -boson mass will be difficult from leptonic decay channel. As a result, in this letter, we have investigated the semi-leptonic decay channel of pairs due to the following advantages: (i) The rate of semi-leptonic final state is relatively large. If we consider only electron and muon decay modes, Pb-Pb collision with TeV gives rise to about 1000 semi-leptonically decaying events for nb integrated luminosity. (ii) Due to the presence of a lepton, the semi-leptonic final state of pairs suffers less from the QCD background. (iii) At the parton level, the semi-leptonic final state contains only -light quark jets arising from the -decay. Therefore, invariant mass can be constructed with out any ambiguity. However, there is a two fold ambiguity in the reconstruction of top quark mass.
After the decay of pairs, the decay products pass through the hot and dense QGP medium and thus, loose energy. The energy loss of energetic partons, so called ”jet quenching”, leads to a number of phenomena which are already seen at RHIC and LHC. In this work, we have investigated the quenching of decay products and proposed some new phenomena which could be seen at the LHC HIC. We have used PYTHIA pythia () to simulate the production and decay of pairs. Subsequently, the PYTHIA generated events are passed in to a fast Monte-Carlo simulator PYQUEN pyquen () for simulating the energy loss (quenching) of decay products. Finally, quarks and gluons are hadronized according to the Lund string mode ṖYQUEN simulates the radiative and collisional (Coll.) energy loss coll_loss () of hard partons in longitudinally expanding QGP taking into account the realistic nuclear geometry. The radiative energy loss is calculated in the framework of BDMS model bdms () with the simple generalization to a massive quark case using the âdead-coneâ approximation dead-cone (). Measuring jet energy as a sum of the energies of final hadrons moving inside an angular cone with a given finite size allow some of the radiated gluons to belong to the jet and thus some part of the radiated energy to be reconstructed. Therefore, the knowledge of angular structure of medium-induced radiation is very important for any phenomenological prediction using jets. In this analysis, we have used the simple parametrizations of the gluon distribution over the emission angle available in PYQUEN: (i) Small-angular radiation (SAR): , where is the typical angle of the coherent gluon radiation as estimated in Ref. SA (); (ii) Wide-angular radiation (WAR): . The strength of the energy loss in PYQUEN is determined mainly by the initial maximal temperature of hot matter in Pb-Pb collisions. The energy loss also depends on the number of active flavors in the medium. In our analysis, we have used GeV and .
In this letter, we have investigated -jets out of which -jets are -tagged plus one lepton signature as signal of production in HIC. Therefore, before going into the details of our analysis, it is important to discuss the status of jet reconstruction and -tagging in the context of HIC. The main obstacle to studying jets in HIC is the presence of the huge background given by the underlying event (UE) This UE needs to be properly subtracted from the momentum of a given jet in order to reconstruct its âtrueâ momentum. It was shown in Ref. Cacciari:2010te () that in presence of QGP, faithful reconstruction of jets are possible using anti- algorithm. Jets have already been successfully used as an observable in HIC by the ATLAS collaboration da_LHC (). Moreover, in Ref. vitev (), different jet shape variables have been studied for the jets passing through QGP. In Ref White:2005au (), -tagging in the environment of HIC has been studied by examining the reconstruction efficiency and rejection power (against light quark jets) using secondary vertex finding. Their study suggests that -tagging efficiency can be achieved at the HIC for a rejection power of 50.
|Energy loss||For distribution||For distribution|
In our analysis, we have introduced a set of basic selection criteria to identify electrons, muons, jets etc. The object selection is described in brief in the following: (i) Jets are constructed using anti- algorithm with and only jets with GeV and are considered for further analysis. To take into account the effects of finite detector resolution, we have smeared the jet energies with Gaussian functions in Ref. smear (). (ii) We demand that lepton candidates (both electron and muon) have GeV and are separated from jets by at least . After reconstructing different objects, we consider events with one lepton (electron or muon) and 4-jets for further analysis. We also demand that out of the 4-jets, two jets are -tagged. The dominant background for semileptonic signal arises from +jets production followed by the leptonic decay of -boson. Here, -jets results from the mistagging of light jets. Since the mistagging efficiency of light quark jets to be tagged as -jets is small, -tagging significantly reduce this background. Production of +jets also contributes to the background. However, we have estimated that these cross-sections are very small compared to the cross-section. Therefore, semileptonic decay products of could be easily detected over the background.
To study the influence of the hot/dense QGP, we have constructed the following kinematic distributions.
Di-jets invariant mass distribution: We have ordered the non jets according to their hardness () and constructed the invariant mass of the hardest () and the second hardest () jet, . In the semileptonic decay of pairs, non -jets arise from the decay of one -boson. Therefore, in absence of QGP, distribution should be peaked at the -mass, GeV. However, in presence of QGP, -decay product suffers energy loss. In Fig. 1, we have presented the distributions for different energy loss scenarios. To determine the peak position of the distributions, we have fitted the distributions with asymmetric Gaussian functions:
where, , and are the parameters of fitting. In Fig. 1 (left panel), we have presented distributions for collisional and radiative (SAR as well as WAR) energy loss scenario separately. Corresponding fitting parameters are presented in Table 1. Fig. 1 (left panel) and Table 1 show that collisional energy loss significantly (about 5 GeV) shifts the peak position of distribution from -mass. Whereas, radiative energy loss (both SAR and WAR) gives rise to small change in the peak position of distribution. Due to the high boost of hard partons, the radiated gluons shift towards the parent parton direction and thus, resulting jets include large part of radiated gluon energies. As a result, radiative energy loss of hard partons has negligible impact on the peak position of distribution. However, Fig. 1 (left panel) shows that WAR significantly changes the shape (which could be quantified by the fitting parameters in Table 1) of distribution. In the passage of a fast parton through QGP, collisional and radiative energy loss occur simultaneously. Therefore, in Fig. 1 (right panel), we have presented distribution in presence of both collisional and radiative energy loss. Corresponding fitting parameters are presented in the last two rows of Table 1. Table 1 shows that the peak position of distribution is determined by the collisional energy loss. Whereas, the shape parameters are governed by the nature of radiative loss. In Fig. 1 (right panel), we have also presented the background. The background contributions are substantially small compared to the contribution.
Tri-jets () invariant mass distribution: Hadronic decay of one top quark gives rise to one -jet and two light quark jets. Therefore, it is important to study invariant mass distribution of jet and di-jets system arising from the hadronic top quark decay. However, it is difficult to identify the jet and di-jets arising from the same top decay. There are several algorithm available in the literature for the reconstruction of top quark. However, most of the algorithms rely on the knowledge of peak position (top mass: GeV). In presence of QGP, energy loss of top decay products shift the peak position of distribution. Therefore, most of the top reconstruction algorithm are not applicable for Pb-Pb collision in their present form. In our analysis, we have used the following simplified algorithm for reconstructing peak position. We first order jets according to their hardness () and constructed two invariant masses: . Out of these two invariant masses, we consider the invariant mass which is close to the GeV for plotting distribution. In Fig. 2, we have presented distributions. In Fig. 2 (left panel), we have shown the effect of collisional and radiative energy loss on distribution separately. Whereas, right panel of Fig. 2 shows the resulting distribution if we consider both collisional and radiative energy loss simultaneously. We have also fitted these distributions with asymmetric Gaussian functions and the fitting parameters are presented in Table 1. Due to the energy loss, the position of peak shifts about 12 GeV from . The shift could be easily observed at the LHC as a signature of QGP.
To summarize, we have investigated semileptonic signature at the LHC HIC. Semileptonic signature could be easily observed at the LHC HIC with TeV. However, due to the presence of QGP decay products suffer energy loss and thus shape of different kinematic distributions are modified significantly. As for example, we have studied dijets and trijets invariant mass distributions and predicted significant change and shift in the shape and peak position of these distributions, respectively.
- (1) B. Muller, The Physics of Quark Gluon Plasma, Springer, Heidelberg, 1985.
- (2) C. Y. Wong, Introduction of High Energy Heavy Ion Collisions, World Scientific, Singapore, 1994; R. C. Hwa (ed.), Quark Gluon Plasma, Vol. I II, World Scientific, Singapore, 1990, 1995.
- (3) J. D. Bjorken, FERMILAB-PUB-82-059-THY (1982).
- (4) J. Adams et al. [STAR Collaboration], Nucl. Phys. A 757, 102 (2005); K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 757, 184 (2005).
- (5) G. Aad et al. [Atlas Collaboration], Phys. Rev. Lett. 105, 252303 (2010); S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 712, 176 (2012).
- (6) C. Amsler et al. (Particle Data Group), Phys.Lett.B 667, (2008).
- (7) Z. Conesa del Valle, Eur. Phys. J. C 61, 729 (2009).
- (8) The ATLAS Collaboration, ATLAS-CONF-2011-078; The ATLAS Collaboration, ATLAS-CONF-2012-052.
- (9) F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 80, 2779 (1998); F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 82, 271 (1999) [Erratum-ibid. 82, 2808 (1999)]; B. Abbott et al. [D0 Collaboration], Phys. Rev. Lett. 80, 2063 (1998).
- (10) http://www.lpthe.jussieu.fr/ cacciari/ttbar/
- (11) M. L. Miller, K. Reygers, S. J. Sanders and P. Steinberg, Ann. Rev. Nucl. Part. Sci. 57, 205 (2007); B. Alver, M. Baker, C. Loizides and P. Steinberg, arXiv:0805.4411 [nucl-ex].
- (12) T. Sjostrandet al., JHEP 0605, 026 (2006).
- (13) I. P. Lokhtin, A.M. Snigirev, Eur. Phys. J. C 45, (2006) 211.
- (14) J. D. Bjorken, Fermilab publication Pub-82/29-THY, 1982; E. Braaten, M. Thoma, Phys. Rev. D 44, (1991) 1298; I.P. Lokhtin, A.M. Snigirev, Eur. Phys. J. C 16, (2000).
- (15) R. Baier, Yu. L. Dokshitzer, A.H. Mueller, D. Schiff, Phys. Rev. C 60, (1999) 064902; R. Baier, Yu. L. Dokshitzer, A.H. Mueller, D. Schiff, Phys. Rev. C 64, (2001) 057902.
- (16) Yu.L. Dokshitzer, D. Kharzeev, Phys. Lett. B 519, (2001) 199
- (17) I.P. Lokhtin, A.M. Snigirev, Phys. Lett. B 440, (1998) 163.
- (18) M. Cacciari, J. Rojo, G. P. Salam and G. Soyez, Eur. Phys. J. C 71, 1539 (2011).
- (19) I. Vitev, S. Wicks and B. -W. Zhang, JHEP 0811, 093 (2008); B. -W. Zhang, Y. He, R. B. Neufeld, I. Vitev and E. Wang, arXiv:1207.6558 [nucl-th]; I. Vitev, J. Phys. G 38, 124087 (2011); Y. He, I. Vitev and B. -W. Zhang, Phys. Lett. B 713, 224 (2012).
- (20) S. N. White, Acta Phys. Hung. A 25, 531 (2006).
- (21) G. Aad et al., [ATLAS Collaboration], arXiv:0901.0512.