Constraints on jet quenching from a multi-stage energy-loss approach
C. Park, A. Angerami, S. A. Bass, S. Cao, Y. Chen, J. Coleman, L. Cunqueiro, T. Dai, L. Du, R. Ehlers, H. Elfner, D. Everett, W. Fan, R. Fries, C. Gale, F. Garza, Y. He, M. Heffernan, U. Heinz, B. V. Jacak, P. M. Jacobs, S. Jeon, W. Ke, E. Khalaj, B. Kim, M. Kordell II, A. Kumar, D. Liyanage, T. Luo, M. Luzum, A. Majumder, M. McNelis, J. Mulligan, C. Nattrass, D. Oliinychenko, L. G. Pang, J.-F. Paquet, J. H. Putschke, G. Roland, B. Schenke, L. Schwiebert, C. Shen, A. Silva, C. Sirimanna, R. A. Soltz, Y. Tachibana, G. Vujanovic, X. -N. Wang, R. L. Wolpert, Y. Xu
CConstraints on jet quenching from a multi-stageenergy-loss approach
Chanwook Park ∗ for the JETSCAPE collaboration Department of physics, McGill University,3600 University Street, Montréal, QC, H3A 2T8, Canada
E-mail: [email protected]
We present a multi-stage model for jet evolution through a quark-gluon plasma within theJETSCAPE framework. The multi-stage approach in JETSCAPE provides a unified descrip-tion of distinct phases in jet shower contingent on the virtuality. We demonstrate a simultaneousdescription of leading hadron and integrated jet observables as well as jet v n using tuned param-eters. Medium response to the jet quenching is implemented based on a weakly-coupled recoilprescription. We also explore the cone-size dependence of jet energy loss inside the plasma. HardProbes20201-6 June 2020Austin, Texas ∗ Speaker © Copyright owned by the author(s) under the terms of the Creative CommonsAttribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND 4.0). https://pos.sissa.it/ a r X i v : . [ nu c l - t h ] S e p onstraints on jet quenching from a multi-stage energy-loss approach Chanwook Park
1. Introduction
Jet evolution through the QGP is characterized by several distinct phases depending on jetvirtualities, and different energy loss mechanisms are essential to describe each stage. A multi-stage approach within the JETSCAPE framework provides a unified description of the jet shower,including a high-virtuality gluon-splitting phase and a low-virtuality scattering-dominated phase.In these proceedings, we report a comprehensive study of multi-stage jet evolution by performinga model-to-data comparison to constrain the jet quenching parameter in heavy-ion collisions.
2. Unified approach in JETSCAPE
Throughout this study, a dynamically evolving QGP created in Pb-Pb collisions at √ s N N = . + Q . The simulation of p+p collisions isperformed by MATTER vacuum showers using the JETSCAPE PP19 tune [1].MATTER is a Monte-Carlo event generator for partons with virtuality Q > Q . Parton splittingsare described by a generalized Sudakov form factor, which includes vacuum and medium-modifiedparton splitting functions. The in-medium contribution, which induces transverse momentumbroadening of jets, ˆ q , in a QGP, is estimated based on the Higher-Twist energy loss model [9–11].We have used a hard thermal loop technique [12] to formulate ˆ q .The time-ordered in-medium shower in LBT for low virtuality partons relies on solving alinearized Boltzmann equation with in-medium kernels. The model contains leading order 2 → → + n inelastic scatterings, where n indicates multiple gluon radiation. The Higher-Twist formalism evaluates the average number of emitted gluons from a hard parton, which followsthe Poisson distribution.The switching virtuality Q is set to 1, 2, and 3 GeV, and a value of α s = .
25 is used forthe strong coupling to determine the quenching parameter ˆ q . Our previous analysis of the singlehadron and jet nuclear modification factor R AA constrained these model parameters [13]. Both theMATTER and the LBT in-medium showers implemented recoil partons based on a weakly-coupledpicture to reproduce the medium response to jet quenching. The energies and momenta originatingfrom incoming thermal partons during jet-medium scattering (holes) are subtracted from the jetsignals in the final state.
3. Results
The left panel of Fig. 1 shows the jet cross section in p+p collisions at √ s N N = .
02 TeVwith two rapidity cuts, normalized by the PYTHIA predictions. The p T dependence of the jetcross-section is consistent with data for jet p T >
200 GeV at mid-rapidity. The ratio of the jet2 onstraints on jet quenching from a multi-stage energy-loss approach
Chanwook Park (GeV)
Tjet p T j e t p d j e t y d σ d R a t i o t o P Y T H I A . ATLAS [PLB 790, 108 (2019)]JETSCAPEPYTHIA 8.230 =5.02 TeV s pp, =0.4 R T k anti-|<0.3 jet y | (GeV) Tjet p T j e t p d j e t y d σ d R a t i o t o P Y T H I A . ATLAS [PLB 790 108 (2019)]JETSCAPEPYTHIA 8.230 =5.02 TeV s pp, =0.4 R T k anti-|<2.8 jet y | − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary (GeV)
Tjet p T j e t p d j e t y d σ d R a t i o t o P Y T H I A . ATLAS [PLB 790, 108 (2019)]JETSCAPEPYTHIA 8.230 =5.02 TeV s pp, =0.4 R T k anti-|<0.3 jet y | (GeV) Tjet p T j e t p d j e t y d σ d R a t i o t o P Y T H I A . ATLAS [PLB 790 108 (2019)]JETSCAPEPYTHIA 8.230 =5.02 TeV s pp, =0.4 R T k anti-|<2.8 jet y | − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary (a) (b) (c) (c) (GeV)
Tjet p T j e t p d j e t y d σ d R a t i o t o P Y T H I A . ATLAS [PLB 790, 108 (2019)]JETSCAPEPYTHIA 8.230 =5.02 TeV s pp, =0.4 R T k anti-|<0.3 jet y | (GeV) Tjet p T j e t p d j e t y d σ d R a t i o t o P Y T H I A . ATLAS [PLB 790 108 (2019)]JETSCAPEPYTHIA 8.230 =5.02 TeV s pp, =0.4 R T k anti-|<2.8 jet y | − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary (GeV)
Tjet p T j e t p d j e t y d σ d R a t i o t o P Y T H I A . ATLAS [PLB 790, 108 (2019)]JETSCAPEPYTHIA 8.230 =5.02 TeV s pp, =0.4 R T k anti-|<0.3 jet y | (GeV) Tjet p T j e t p d j e t y d σ d R a t i o t o P Y T H I A . ATLAS [PLB 790 108 (2019)]JETSCAPEPYTHIA 8.230 =5.02 TeV s pp, =0.4 R T k anti-|<2.8 jet y | − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary (a) (b) (c) (c)
Figure 1:
Comparison between the results obtained from the JETSCAPE PP19 tune at √ s N N = | y jet | < . R = . . R = . − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary (a) (b) (c)(d) (e) (f)
Figure 2:
Inclusive jet R AA in central (top panel) and peripheral (bottom panel) Pb-Pb collisions at √ s N N = .
02 TeV with various R and Q values [15]. cross-section with various R with respect to R = . R AA is well reproduced by the MATTER vacuum shower inJETSCAPE with the PP19 tune.We present the jet R AA with various R and switching virtualities Q in central and peripheralPb+Pb collisions in Fig. 2. We consistently observe stronger jet quenching with larger values of Q .The parton shower in the high-virtuality phase (MATTER) is dominated by virtuality splitting, butthe low-virtuality phase (LBT) is largely affected by scatterings, which induce jet p T broadening.This accounts for the jet R AA being more suppressed when the LBT phase starts at higher virtuality Q . The jet R AA independent to R leads to the R AA ratio with respect to R = . p T ,implying that the jet energy contained within R < . R AA value. Thesteeply falling jet shape function shown in the left panel of Fig. 4 supports this interpretation.3 onstraints on jet quenching from a multi-stage energy-loss approach Chanwook Park (a) (b) (c) (d) − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary
Figure 3:
Ratio of jet R AA as a function of R with respect to R = . p T intervals. The data is calculated from the jet R AA results shown in Fig. 7in [15]. − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary . . . . . . r P ( r )( G e V ) PbPb 5.02 TeV, 0-10% p jet T >
120 GeV , | ¥ jet | < . p trk T > . − −
10 10.811.21.41.61.8
ATLAS (0-10%) =2 GeV) Q MATTER+AdS/CFT ( =1 GeV) Q MATTER+AdS/CFT ( − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+MARTINIMATTER(vacuum)+MARTINI − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT D P b P b ( z ) / D pp ( z ) z z z JETSCAPE,2.76TeV,PbPb : 0 - - k T R = 0.4,100 < p jetT < 398GeV,0 < | Y jet | < 2.1, p trkT > 1GeV − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY − −
10 10.811.21.41.61.8
ATLAS (0-10%)MATTER+LBTMATTER(vacuum)+LBT
JETSCAPE PRELIMINARY (a) (b) (c)
Preliminary Preliminary Preliminary . . . . . . r P ( r )( G e V ) PbPb 5.02 TeV, 0-10% p jet T >
120 GeV , | ¥ jet | < . p trk T > . . . . . . . r P ( r )( G e V ) PbPb 5.02 TeV, 0-10% p jet T >
120 GeV , | ¥ jet | < . p trk T > . (a) (b) anti- k T , R = 0 . Figure 4: (Left) Jet shape function for R = . v and v for jets at peripheral Pb-Pb collisions [ ? ]. However, a rigorous investigation of recoils would be necessary as their influence on jet shape isexpected to be significant at larger R .The right panel in Fig. 4 shows the jet v and v in peripheral Pb-Pb collisions. The observednon-zero v for high-energy jets originates from the path-length dependent jet quenching in analmond-shaped QGP. The vanishing jet v within the statistical uncertainties is consistent with thedata.
4. Conclusion
We have studied jet modification using a unified approach within the JETSCAPE framework.The results for the jet cross-section in pp collisions using the JETSCAPE PP19 tune show goodagreement with data. The multi-stage model with a combination of MATTER and LBT providesa simultaneous description of the integrated and differential jet observables. Our future work willinvestigate recoils for the detailed jet quenching mechanism at large R . References [1] A. Kumar et al. [JETSCAPE], [arXiv:1910.05481 [nucl-th]].4 onstraints on jet quenching from a multi-stage energy-loss approach
Chanwook Park[2] C. Shen, Z. Qiu, H. Song, J. Bernhard, S. Bass and U. Heinz, Comput. Phys. Commun. ,61-85 (2016) doi:10.1016/j.cpc.2015.08.039 [arXiv:1409.8164 [nucl-th]].[3] J. S. Moreland, J. E. Bernhard and S. A. Bass, Phys. Rev. C , no.1, 011901 (2015)doi:10.1103/PhysRevC.92.011901 [arXiv:1412.4708 [nucl-th]].[4] T. Sjöstrand, S. Ask, J. R. Christiansen, R. Corke, N. Desai, P. Ilten, S. Mrenna, S. Pres-tel, C. O. Rasmussen and P. Z. Skands, Comput. Phys. Commun. , 159-177 (2015)doi:10.1016/j.cpc.2015.01.024 [arXiv:1410.3012 [hep-ph]].[5] A. Majumder, Phys. Rev. C , 014909 (2013) doi:10.1103/PhysRevC.88.014909[arXiv:1301.5323 [nucl-th]].[6] R. Abir and A. Majumder, Phys. Rev. C , no.5, 054902 (2016)doi:10.1103/PhysRevC.94.054902 [arXiv:1506.08648 [nucl-th]].[7] X. N. Wang and Y. Zhu, Phys. Rev. Lett. , no.6, 062301 (2013)doi:10.1103/PhysRevLett.111.062301 [arXiv:1302.5874 [hep-ph]].[8] Y. He, T. Luo, X. N. Wang and Y. Zhu, Phys. Rev. C , 054908 (2015)doi:10.1103/PhysRevC.91.054908 [arXiv:1503.03313 [nucl-th]].[9] X. N. Wang and X. f. Guo, Nucl. Phys. A , 788-832 (2001) doi:10.1016/S0375-9474(01)01130-7 [arXiv:hep-ph/0102230 [hep-ph]].[10] A. Majumder, Phys. Rev. D , 014023 (2012) doi:10.1103/PhysRevD.85.014023[arXiv:0912.2987 [nucl-th]].[11] G. Y. Qin and A. Majumder, Phys. Rev. Lett. , 262301 (2010)doi:10.1103/PhysRevLett.105.262301 [arXiv:0910.3016 [hep-ph]].[12] S. Caron-Huot and C. Gale, Phys. Rev. C , 064902 (2010) doi:10.1103/PhysRevC.82.064902[arXiv:1006.2379 [hep-ph]].[13] A. Kumar et al. [JETSCAPE], [arXiv:2002.07124 [nucl-th]].[14] M. Aaboud et al. [ATLAS], Phys. Lett. B , 108-128 (2019)doi:10.1016/j.physletb.2018.10.076 [arXiv:1805.05635 [nucl-ex]].[15] [CMS], CMS-PAS-HIN-18-014.[16] A. M. Sirunyan et al. [CMS], JHEP05