arxiv: v1 [nucl-th] 5 Dec 2018

Transcription

1 EPJ manuscrit No. (will be inserted by te editor NLO Productions of ω and KS wit a Global Extraction of te Jet ransort Parameter in Heavy Ion collisions Guo-Yang Ma, Wei Dai a, Ben-Wei Zang b, and Enke Wang arxiv:8.33v [nucl-t] 5 Dec 8 Key Laboratory of Quark & Leton Pysics (MOE and Institute of Particle Pysics, Central Cina Normal University, Wuan 4379, Cina Scool of Matematics and Pysics, Cina University of Geosciences, Wuan 4379, Cina Introduction Received: date / Revised version: date Abstract. In tis work, we ave te way to calculate te roductions of ω and K S mesons at large in + and A+A collisions at te RHIC and te LHC. e ω meson fragmentation functions (FFs in vacuum at next-to-leading order (NLO are obtained by evolving NLO DGLAP evolution equations wit rescaled ω FFs at initial scale Q =.5 GeV from a broken SU(3 model, and te K S FFs in vacuum are taken from AKK8 arametrization directly. Witin te framework of te NLO QCD imroved arton model, we make good descritions of te exerimental data on ω and K S in + bot at te RHIC and te LHC. Wit te iger-twist aroac to take into account te jet quencing effect by medium modified FFs, te nuclear modification factors for ω meson and K S meson at te RHIC and te LHC are resented wit different sets of jet transort coefficient ˆq. en we make a global extraction of ˆq at te RHIC and te LHC by confronting our model calculations wit all available data on 6 identified mesons: π, η, ρ, φ, ω, and K S. e minimum value of te total χ /d.o.f for roductions of tese mesons gives te best value of ˆq =.5 GeV /fm for Au+Au collisions wit s NN = GeV at te RHIC, and ˆq =. GeV /fm for Pb+Pb collisions wit s NN =.76 ev at te LHC resectively, wit te QGP sacetime evolution given by an event-by-event viscous ydrodynamics model IEBE-VISHNU. Wit tese global extracted values of ˆq, te nuclear modification factors of π, η, ρ, φ, ω, and K S in A+A collisions are resented, and redictions of yield ratios suc as ω/π and K S/π at large in eavy-ion collisions at te RHIC and te LHC are rovided. PACS..38.M Quark- lasma q Relativistic eavy-ion collisions 3.85.Ni Inclusive roduction wit identified adrons Jet quencing effect describes te energy dissiation of an energetic arton wen it traverses troug te ot and dense QCD medium, wic is roduced sortly after te ig energy nuclear collisions []. e single adron roduction suression at ig region wen comared to scaled + data is one of te rimary erturbative robes to study te roerties of tis de-coule and QCD matter []. π as te most well measured final state adron, its nuclear modification factor as a function of transverse momentum is interreted as te consequence of jet quencing effect and it did el us to constrain te strengt of te jet-medium interaction and also te roerties of te QCD medium [3, 4,5,6,7,8]. In 3, Jet collaboration summarized different jet quencing teoretical frameworks and comared te jet transort arameters extracted by different energy a weidai@cug.edu.cn b bwzang@mail.ccnu.edu.cn loss models using te same ydro descrition of te QCD medium [9]. Meanwile, te for different final state identified adrons ave been measured, and teir roduction suressions as well as atterns of teir yield ratios ave been observed [4,5,,,,3,4,5]. It is of interest and callenge to describe te cross sections of leading adron of different tyes and teir yield ratios wit eac oter in eavy-ion collisions (HIC at te RHIC and te LHC wit a unified model of jet quencing, wic sould sed ligt on flavor deendence of arton energy loss and te intrinsic roerties of identified adron roductions in + and A+A reactions [6,7,8,9,,,,3]. Previous efforts on te leading adron roductions in A+A collisions mainly focused on te roduction suression of ions. However we ave acieved more understanding of te suression atterns of different mesons by conducting calculations and analysis of te and article ratios for oter identified adron roductions [4, 5, 6]. In a revious study [4] we ave firstly calculated te η roduction in A+A collisions bot at te RHIC and te LHC, teoretically exlained te coincidence of te η/π

2 Guo-Yang Ma et al.: ω and K S Productions wit a Global Extraction of te Jet ransort Parameter ω Q =.5 GeV - - ω z =.4-6 K s Q = 4. GeV - K s z = ω Q =. GeV z - - ω z = Q K s Q =. GeV z K s z = Q Fig.. Left: e NLO DGLAP evolved Fragmentation Functions of ω meson as functions of z at fixed Q =.5 GeV and Q = GeV. Rigt: e NLO DGLAP evolved Fragmentation Functions of ω meson as functions of Q at fixed z =.4 and z =.6. Fig.. Left: e NLO DGLAP evolved Fragmentation Functions of K S meson as functions of z at fixed Q = 4. GeV and Q = GeV. Rigt: e NLO DGLAP evolved Fragmentation Functions of K s meson as functions of Q at fixed z =.4 and z =.6. ratios in + and Au+Au collision at larger region, and furter redict te similar overlaing of η/π ratios in + and Pb+Pb wit resect to at te LHC, wic as been very recently confirmed by ALICE measurement on π and η mesons at te LHC [3]. is teoretical aroac as later been extended to investigate te roductions of ρ meson and φ meson at large in eavy-ioncollisions[5,6],andsowstatρ /π ratiosin A+A sould aroac to tat in + at larger region, wereas te coincidence of yield ratios φ/π in te + and A+A article will not sow u over a wide region of [6]. We emasize tat inclusiveleadingadronyields in eavy-ion collisions sould be determined by tree factors: te initial ard arton sectrum, te arton energy loss mecanism, and te secific sae of arton fragmentation functions in vacuum; and te yield ratios of identified adrons at large are te combined effects of tese tree factors, not determined solely by one factor [4]. In tis article, we investigate te roduction of oter twomesonsω and KS atlarge in A+Acollisions,wic as never been comuted so far, to te best of our knowledge. e ω meson is constituted of similar valence quark of π wit larger mass MeV and sin. Kaons are a grou of ligtest mesons, carrying strangeness comonents, KS (d s+s d is one tye of te Kaons, wic is consisted by s quark, d quark and teir corresonding antiquarks. e roductions of tese two mesons ave been measured in + and A+A collisions at te RHIC and also at te LHC, but lacking te teoretical descrition. Witin te NLO QCD imroved arton model, we calculate ω and KS yields at large in eavy-ion collisions,by emloying te medium-modification fragmentation functions (FFs due to radiation in te ot/dense QCD in te iger-twist aroac of jet quencing [,, 7, 8,9], te same aroac as in te calculations on roductions of η, ρ and φ mesons in HIC [4,5,6]. It is noted in te revious studies, for te consistency, we imlemented Hirano ydro descrition [3, 3] to de- scribe te sace-time evolution of te QGP fireball. In tis work, we may utilize a state-of-art, event-by-event (+- D viscous ydrodynamics model (IEBE-VISHNU [3] to give te sace-time evolution information of te ot and dense medium. Due to te cange of te medium descrition, te strengt of te jet-medium interaction caracterized by te jet transort arameter ˆq sould be reextracted. aking te advantage of te systematical study of te nuclear modification factors wit resect to and a large amount of exerimental data of π, η, ρ and φ bot at te RHIC and te LHC, wit two more mesons ω and KS calculated in tis article, we can make a global extraction of te jet transort arameter ˆq wit all te available exerimental data of tese 6 identified mesons. e article is organized as follows. We firstly resent te teoretical framework of comuting single adron cross sections in + collision in Sec., and give te + baseline in investigating te in-medium modification of tese roductions. In Sec. 3, we discuss te inclusive adron roduction in A+A collisions wit te medium-modified FFs, basing on te iger-twist aroac of arton energy loss, and resent te nuclear modification factors for ω and KS at te RHIC and te LHC. Sec. 4sowsaglobal extraction of jet transort coefficient ˆq at te RHIC and te LHC by confronting our model calculations wit all available data on 6 identified mesons: π, η, ρ, φ, ω, and KS. In Sec.5,wemakeredictionsonte identified adron yield ratios ω/π and KS /π in + and A+A collisions, and comare tem wit exerimental data if alicable. We give a brief summary in Sec. 6. Large yield of ω and K S meson in + We start wit te roductions of ω and KS mesons at NLO in + collisions. In te framework of te QCD imroved arton model, te inclusive adron roduction

3 Guo-Yang Ma et al.: ω and K S Productions wit a Global Extraction of te Jet ransort Parameter 3 c 3 - (mb GeV 3 σ/d 3 Ed PHENIX ω RHIC ALICE (7GeV ω LHC + -> ω+x Fig. 3. e teoretical results of te ω meson roduction in + collisions at s = GeV confronted wit te PHENIX data [4] and confronted wit te ALICE data [33] at s = 7 ev. c 3 - (mb GeV 3 σ/d 3 Ed K s RHIC SAR LHC ALICE (76GeV K s + -> K +X s Fig. 4. e teoretical results of te K s meson roduction in + collisions at s = GeV RHIC comared wit SAR data [34] and at s =.76 ev LHC comared wit ALICE data [5]. in + collision can be given by te convolution of tree arts, arton distribution functions (PDFs in a roton, ard artonic scattering cross section denoted as dσ/dˆt (u to te order of α 3 s, and te arton fragmentation functions (FFs to te final state adron D q(g (z. One may convolute PDFs and artonic cross section dσ/dˆt into te initial ard (arton- jet sectrum F q,g (, and we ave: dσ ω,k S = d + F q ( z D q ω,k S (z, dz z F g ( z D g ω,k S (z, dz z. ( in wic te quark and fragmenting contributions are written searately to facilitate future discussions. A next-to-leading order (NLO Monte Carlo code as been emloyed to calculate leading adron roductions in + collision [35], and te C4 arametrization of PDFs for free roton [36] as been imlemented. In te framework, as long as te arametrization of te secific final state adron FFs in vacuum is available, one can redict its roduction yield at large in elementary + collision in rincile. In order to make te NLO calculation of ω and KS mesons in +, te NLO arton FFs of tese two mesons are needed. e arton FFs of KS meson at NLO can be found in AKK8 arametrization [37], wile tere is not suc kind of global arametrization for ω FFs in vacuum. So we aveto rely on teoretical models and utilize ω arton FFs at a starting scale Q =.5 GeV rovided by a broken SU(3 model [38,39]. And we note tat our receding investigations on leading ρ and φ roduction [5, 6] ave benefited from tis broken SU(3 model. In tis model, te indeendent arton FFs of different flavor are reduced into 3 indeendent functions named as valence V(x, sea γ(x and D g (x by considering te SU(3 flavor symmetry wit a symmetry breaking arameter also te isosin and carge conjugation invariance of te vector mesons. We can write tese functions into a standard olynomial at te starting low energy scale of Q =.5 GeV as: F i (x = Ka i x bi ( x ci (+d i x+e i x ( Were te full set of arametersdefined in V, γ, D g and a few additional arameters defined for eac vector meson suc as strangeness suression factor λ, vector mixing angle θ, sea suression factor f ω sea, fu (ω and fω g ave been determined in te broken SU(3 model in Ref. [38, 39], and K is a rescaling factor to make te best fit to leading adron sectrum in + at te RHIC, and we fix K = 3 for ω meson in our NLO numerical simulations. o obtain a NLO arton FFs for ω at any energy scale Q, we emloy a numerical NLO DGLAP evolution rogram rovided in Ref. [4] wit te initial arton FFs starting scale Q =.5 GeV as inut. We demonstrate in Fig. te NLO DGLAP evolved FFs of ω meson as functions of z at fixed Q =.5 GeV and Q = GeV on te left and as functions of Q at fixed z =.4 and z =.6 on te rigt. We find in te tyical fraction region z =.4.7 D ω g > Dω s Dω u(d, unlike te case Ds φ > Dg φ D φ u(d in te φ FFs [6]. Due to te fact tat in arton FFs, exceed contribution, an overwelming advantage of fragmenting contribution of te final state ω roduction in + collision is exected in te cometition wit quark fragmenting contribution. In Fig., we lot te NLO FFs of KS from AKK8 arametrization in te same manner as ω meson in Fig.. We find te strange quark FF D K S s > D K S g in te tyical z =.4.7 region, terefore a cometition of te an fragmenting contributions is exected in te KS roduction in +

4 4 Guo-Yang Ma et al.: ω and K S Productions wit a Global Extraction of te Jet ransort Parameter.5.5 S=GeV, y < P ω: osu(vis=.8,q =.3 GeV osu(vis=.8,q =.4 GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =.6 GeV.5 osu(vis=.8,q =. GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =. GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =3.5 GeV.5 S=76GeV, y < P Fig. 5. Left: Comarison between te PHENIX Data [4] of ω nuclear modification factor in Au + Au collisions at GeV and numerical simulations at NLO. Rigt: e numerical rediction of ω nuclear modification factor in Pb + Pb collisions at 76 GeV at NLO..5.5 K sort : S=GeV, y <.5 osu(vis=.8,q =.3 GeV osu(vis=.8,q =.4 GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =.6 GeV P.5.5 K sort : S=76GeV, y <.8 osu(vis=.8,q =. GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =. GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =3.5 GeV P Fig. 6. Left:Comarison betweentesardata[5]ofk S nuclear modification factor in Au + Au collisions at GeV and numerical simulations at NLO. Rigt: Comarison between te ALICE Data [5] of K S nuclear modification factor in Pb + Ob collisions at 76 GeV and numerical simulations at NLO. collision based on te attern discovered in te revious study on π, η, ρ and φ meson. Wit te NLO FFs in vacuum of bot ω and KS meson, we are able to confront te numerical results of te ω and KS meson roductions in + collisions u to te NLO wit te available exerimental data bot at te RHIC and at te LHC sown in Fig. 3 and Fig. 4. Note tat in our simulations, te ard scales suc as factorization scale, renormalization scale and fragmentation scale are cosen to be te same and roortional to te of te final state adron. We see te numerical results can matc well wit te exerimental data on ω meson sectra bot at te RHIC s NN = GeV and te LHC snn = 7 ev wit te rescaling factor K = 3 and te ardscalesµ =.5, assownin Fig. 3. We alsodemonstrate in Fig. 4 te good agreements of te teoretical calculations of leading KS meson, wit te SAR data and te ALICE data, were te ard scales are fixed to be µ =. to give te best fit. 3 Large yield of ω and K S Meson in A+A ω: o calculate te ω and KS roductions in A+A collisions, jet quencing effect sould be included. In te article we utilize te iger-twist aroac of arton energy loss, wic relates arton energy loss due to its multile scattering in QCD medium and medium-induced radiation to twist-four rocesses, and sows tat tese rocesses may give rise to additional terms in QCD evolution equations and lead to te effectively medium-modified fragmentation functions; tus te artonic energy loss effect can be taken into accountbyrelacingte FFs in vacuum wit te effectively medium modified FFs [7, 8, 9, 4, 4]. By assuming a termal ensemble of quasi-article states in te ot/dense medium and also neglecting te multile article correlations inside te medium, te effective medium modified FFs are given as [,,4,5, 6]: Q D q (z = Dq (z + α s(q dl π l dz [ γ q qg (z,x,x L,l z z Dq( z z ] + γ q gq (z,x,x L,l D g (z z. (3 e contribution of te medium-induced radiation is attributed to te medium modified slitting functions reresented by γ q qg and γ q gq. In tis formula, we assume te energy loss of te jet roagating troug te medium is totally carried away by te radiative s, ten te convolution of tese energy loss kernels γ q qg ( γ q gq wit te (DGLAP evolved vacuum FFs at any scale D q,g( z z imlies te assumtion tat te fast arton first loses its energy in te medium and ten fragments into final state adrons in vacuum. e medium modified slitting functions deend on te roerties of te local medium wic can not be determined directly by te teoretical calculation itself. e jet transort arameter ˆq wic defined as te average squared transverse momentum broadening er unit lengt is terefore introduced to rofile te deendency of te local medium roerties in tese energy loss kernels [9]. e jet transort arameter ˆq can be enomenological assumed to be roortional to te local arton density in te ot/dense medium to adot furter medium descrition: ˆq(τ,r = ˆq ρ med (τ,r ρ med (τ, µ u µ, (4 were ρ med is te arton (quarks and density in te medium at a given temerature, q is te initial jet transort arameter at te center of te bulk medium at te initial time τ, µ is te four-momentum of te jet and u µ is te four flow velocity of te robed local medium in te collision frame. In tis article we emloy te stateof-art, event-by-event viscous ydrodynamics descrition IEBE-VISHNU [3] to give te sace-time evolution information of te ot/dense medium suc as temerature, arton density, four velocity of te local medium at any evolution time and osition, as well as te formation time of te QGP τ =.6 fm. erefore te only undetermined arameteristecalculationis ˆq reresentingtestrengt of te jet-medium interaction, wic may be fixed by fitting te data of leading adrons (usually π meson suression in A+A collisions [9, 4, 5, 6]. Wen averaging te medium modified FFs over te initial roduction osition and jet roagation direction,

5 Guo-Yang Ma et al.: ω and K S Productions wit a Global Extraction of te Jet ransort Parameter 5 we are able to directly relace te vacuum FFs in te + formalism wit tese averaged medium modified FFs D c(z,e,b, terefore te formalism of te calculated cross section of te single adron roductions in HIC would take te form as: π (ALICE: S=76GeV, y <.8 osu(vis=.8,q =. GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =. GeV osu(vis=.8,q =3. GeV π osu(vis=.8,q =. GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =. GeV osu(vis=.8,q =3. GeV ρ N AB coll dσ AB (b dyd = abcd dx a dx b f a/a (x a,µ f b/b (x b,µ dσ dˆt (ab cd D c (z,e,b πz c +O(α 3 s. (5 in wic N AB coll (b = d rt A (rt B ( b r is te number of binary nucleon-nucleon collisions at certain imact arameter b in A+B collisions and it can be calculated by using Glauber Model [43]. f a/a (x a,µ reresents te effective PDFs inside a nucleus. In our calculations, we emloyed EPPS6 NLO nuclear PDFs [44] to include initialstate cold nuclear matter effects on single adron roductions [45] π : S=GeV,PHENIX(8 π : S=GeV,PHENIX( osu(vis=.8,q =.4 GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =.6 GeV osu(vis=.8,q =.7 GeV η: S=GeV,PHENIX osu(vis=.8,q osu(vis=.8,q =.5 GeV osu(vis=.8,q =.6 GeV osu(vis=.8,q =.7 GeV =.4 GeV ω: S=GeV,PHENIX osu(vis=.8,q =.4 GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =.6 GeV osu(vis=.8,q =.7 GeV π η ω ρ : S=GeV,SAR osu(vis=.8,q =.4 GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =.6 GeV osu(vis=.8,q =.7 GeV φ: S=GeV,PHENIX osu(vis=.8,q osu(vis=.8,q =.5 GeV osu(vis=.8,q =.6 GeV osu(vis=.8,q =.7 GeV =.4 GeV K s: S=GeV,SAR osu(vis=.8,q =.4 GeV osu(vis=.8,q =.5 GeV osu(vis=.8,q =.6 GeV osu(vis=.8,q =.7 GeV Fig. 7. eoretical calculation results of nuclear modification factors as functions of at ˆq =.4.7 GeV /fm confronted wit te RHICexerimental data of π [4,5], ρ [5], K S [5], η [] and φ [], ω [4]. e nuclear modification factor as a function of is calculated as cross sections in A + A collisions divided by te ones in + collision, scaled by te averaged number of binary nucleon-nucleon collisions wit a cosen imact arameter b [3]: R b AB (,y = dσ AB /dyd N AB coll (b dσ /dyd. (6 e teoretical results of at various values of ˆq =.4.7 GeV /fm for bot ω and K s mesons ave been resented in Fig. 5 and Fig. 6. ρ φ K s η: S=76GeV,ALICE osu(vis=.8,q =. GeV.6 osu(vis=.8,q =.5 GeV osu(vis=.8,q =. GeV η.4 osu(vis=.8,q =3. GeV osu(vis=.8,q =. GeV.6 ω osu(vis=.8,q =.5 GeV.4 osu(vis=.8,q =. GeV. osu(vis=.8,q =3. GeV φ: S=76GeV,ALICE osu(vis=.8,q osu(vis=.8,q =.5 GeV osu(vis=.8,q =. GeV osu(vis=.8,q =3. GeV =. GeV K s: S=76GeV,ALICE osu(vis=.8,q =. GeV osu(vis=.8,q =.5 GeV K s osu(vis=.8,q =. GeV osu(vis=.8,q =3. GeV φ Fig. 8. eoretical calculation results of nuclear modification factors as functions of at ˆq =. 3. GeV /fm confronted wit te LHC exerimental data of π [], η [3], φ [5], K S [5], ρ and ω. Wen obtaining suitable value of ˆq by comaring te teoretical calculations of for ω and K s wit te corresonding data, tere exist two cautions. Firstly, since te data of ω and KS meson at large in A+A collisions are rater limited and wit large uncertainty, it is difficult to make a good constrain on ˆq wit tese data. Secondly, because we emly te IEBE-VISHNU ydrodynamics model to describe te sace-time evolution of te fireball, wic may give different information of ysics quantities suc as temerature and density from tose rovided by oter ydro models suc as Hirano ydro descrition [3, 3]. erefore, we could not take advantage of te extracted value of ˆq in Ref. [9,,,4,5,6], were Hirano ydro descrition as been utilized. Wit tese two cautions in mind, and realizing tat wit our model we are now ready to make a systematic studyof6tyesofidentifiedmesonssucasπ,η,ρ,φ,ω, and KS in eavy-ion collisions, it will be of great interest to make a global extraction of jet transort coefficient ˆq at te RHIC and te LHC by confronting our model calculations (wit te IEBE-VISHNU ydro model against all available data on 6 identified mesons: π, η, ρ, φ, ω, and KS, and ten make recise calculations on nuclear modification modification factors of tese mesons including ω, and KS, as well as yield ratios of tese 6 mesons in HIC. 4 Global Extraction of ˆq wit for Six Identified Mesons We erform a systematic calculation of te of 6 identifiedmesons(π,ρ,η,φ,ω,k S tocomarewitallteir

6 6 Guo-Yang Ma et al.: ω and K S Productions wit a Global Extraction of te Jet ransort Parameter χ /d.o.f χ /d.o.f π ρη φ ω K s RHIC-GeV q (GeV S=GeV, y <.35 S=76GeV, y <.8 error bars. We also derive te best value of ˆq in te Pb+Pb collisions at te LHC s NN =.76 ev to be ˆq =. GeV /fm, toug te curve χ /d.o.f at te LHC is muc flatter tat tat at te RHIC. So teoretical results wit ˆq =..4 GeV /fm sould all give decent descritions on data at te LHC. It is noted tat in our current model te extracted values of jet transort coefficient ˆq at te RHIC and te LHC are smaller tan tat by JE Collaboration in Ref. [9] as well as te ones in our receding calculations [4, 5, 6]. ese differences come mainly from te different ydro models utilized in te studies. We ave cecked[46] tat if we emloy Hirano ydro descrition in te calculation, te extracted values of ˆq from our global fitting will be consistent wit tose in Ref. [9,4,5,6], toug te model in tis article as te otential to give more recise extraction of jet transort coefficients wen more data of identified adrons in A+A collisions become available in te near future. 5 5 Particle Ratios of ω/π and K s /π in A+A q (GeV Fig. 9. o: demonstration of te χ /d.o.f between te teoretical rediction for of π, ρ, η, φ, ω, K S at different ˆq and te current ublic exerimental data at te RHIC GeV; Bottom: te global χ /d.o.f of taking into account all tese 6 identified adrons at te RHIC GeV and LHC.76 ev. available exerimental data at te RHIC and te LHC in Fig. 7 and Fig. 8. In order to extract te best value of te initial jet transort arameter ˆq, we erform a χ fit to comare te teoretical results wit different ˆq and te available exerimental data of all te different final adrons. χ (a = i [D i i (a] σ i (7 In te aboveequations,d i reresentste exerimental grids and i is our teoretical rediction at inut arameter a. σ i means te i t systematic and statistical exerimental errors. We sow in te to anel of Fig. 9, te derived χ averaged by te number of te comared data oints for different final state mesonsat various ˆq =.4.7 GeV /fm at te RHIC wit s NN = GeV. In te bottomaneloffig.9welottecurveχ /d.o.f asafunction of ˆq at te RHIC and te LHC, were te minimum of te curve χ /d.o.f wit resect to ˆq resents te best fit of teory wit data. We ten observe tat at te RHIC te minimum oint of χ /d.o.f of all 6 identified mesons gives te best value of ˆq =.5 GeV /fm. It is noted tat te roduction of π will carry te largest weigt due to its more abundant data wit relatively smaller Wit te global extracted value of ˆq discussed in Sec. 4, we are able to furter investigate te article ratio of ω and KS bot at te RHIC and te LHC. We first calculate ω/π ratio as a function of and sow te results in + and Au+Au at te RHIC in te left anel of Fig., were te PHENIX exerimental data on ω/π ratio in + are also illustrated. An enancement of te ratio in A+A relative to tat in + is found in small region, wereas a small suression in large region. We also redict te ω/π ratio as a function of in + and Pb+Pb collisions wit s NN =.76 ev at te LHC. In Fig. we don t see te overlaing of te curves ω/π in + and A+A at large, as we ave observed for te ratios η/π [4] and ρ /π [5] in + and A+A. o understand tis feature deeer, we lot te an(fragmentation contribution fractions to ω (and π in + and Au+Au at te RHIC in Fig.. One see in te + collision at te RHIC, te roduction of ω is dominated by fragmentation. In te A+A collisions, tesuffersmoreenergylosstante ligtquarkdue to its larger color factor, wic decreases te contribution fraction and increase te ligt quark contribution fraction, but te dominant contribution fraction of ω 6% is still observed u to = GeV in Au+Au at te RHIC. On te oter and, π meson roduction is ligt quark fragmentation dominant in +, and jet quencing effect furter enance tis dominance of ligt quark fragmentation to π in A+A. erefore, we see te yield ratio of ω/π in A+A sould be suressed relative to tat in + due to energy loss effect, and it searates wit te one in + even at very ig. We also comute te KS /π ratio as a function of bot at te RHIC and te LHC in Fig.. We find te curves in A+A and in + are aroacing to eac oter wit increasing, and an obvious coincidence of tese two curves is seen at te LHC. We sow te

7 Guo-Yang Ma et al.: ω and K S Productions wit a Global Extraction of te Jet ransort Parameter 7 ω/π 3 - PHENIX ω->π π + π PHENIX ω->π γ GeV q =.5 GeV /fm ω/π GeV q =. GeV /fm q =. GeV /fm q =.4 GeV /fm Fig.. Left: te yield ratios of ω/π as functions of in + and Au+Au collisions wit GeV at te RHIC, and PHENIX data in + [4]; Rigt: redictions of te yield ratios of ω/π as functions of in + and Pb+Pb collisions wit.76 ev at te LHC.. + ω-g ω-q π -g π -q. ω-g Au+Au ω-q π-g π -q Fraction.8.6 Fraction Fig.. Gluon an contribution fractions of te total yields of ω and π mesons in + and Au+Au collisions wit snn = GeV at te RHIC. K s/π.5 GeV q =.5 GeV /fm K s/π GeV q =. GeV /fm q =. GeV /fm q =.4 GeV /fm Fig.. Left: te roduction ratios of K S/π as functions of in + and Au+Au collisions wit GeV at te RHIC; Rigt: redictions of te roduction ratios of K S/π as functions of in + and Pb+Pb collisions wit.76 ev at te LHC.

8 8 Guo-Yang Ma et al.: ω and K S Productions wit a Global Extraction of te Jet ransort Parameter Fraction. K s -g + K s -q π -g π -q Fraction. K s -g Au+Au K s -q π -g π -q Fig. 3. Gluon an contribution fractions of K S and π yields in + and Au+Au collisions wit s NN = GeV at te RHIC. an contribution fractions to K S ( and π meson as functions of in Fig. 3. We find in + collision, te roductions of bot K S and π at very large are dominated by quark fragmentation. In A+A collisions, te contribution may be furter suressed because generally loses more energy. us, bot in + and A+A collisions, te ratio K S /π sould be largely determined by te ratio of quark FFs for KS (DK S q (z to quark FFs for π (Dq π(z at very ig, were tese FFs vary slowly wit te momentum fraction z, very similar to te case of η /π at ig [4]. Even toug in A+A collisions, jet quencing effect may sift z of quark FFs, if quark FFs ave a rater weak deendence on z and, we may see at veryig te curves KS /π in A+A and + curves may come close to eac oter, and even coincide at te LHC. 6 Summary In summary, we obtain te NLO FFs of ω meson in vacuum by evolving te rescaled ω FFs from a broken SU(3 model at a starting scale Q =.5 GeV, and directly emloy NLO KS FFs in vacuum from te AKK8 arameterizations, te numerical simulation of roductions of bot ω and KS matces well wit te exerimental data in + reactions. Wit te IEBE-VISHNU ydro rofile of te QCD medium, we calculate te nuclear modification factors of ω and KS meson as well as π, η, φ, ρ in A+A collisions at te RHIC and te LHC, including jet quencing effect in iger-twist aroac. e global extraction of jet transort arameter ˆq are made, wit comarison of te teoretical calculation and te exerimental data of all six identified mesons: π, ρ, η, φ, ω, KS. Furtermore,weredict teyield ratiosofω/π bot at te RHIC and te LHC, and a fairly good agreement of te teoretical results and exerimental data is found at te RHIC. eoretical redictions of KS /π ratios as functions of at te RHIC and te LHC are resented as well. Acknowledgments: e researc is suorted by te NSFC of Cina wit Project Nos and 8567, and artly suorted by te Fundamental Researc Funds for te Central Universities, Cina University of Geosciences (Wuan (No References. X. N. Wang and M. Gyulassy, Pys. Rev. Lett. 68, 48 (99.. M. Gyulassy, I. Vitev, X. N. Wang and B. W. Zang, In *Hwa, R.C. (ed. et al.: Quark lasma* 3-9 [nucl-t/377]. 3. S. S. Adler et al. [PHENIX Collaboration], Pys. Rev. Lett. 96, 3 (6 [nucl-ex/637]. 4. A. Adare et al. [PHENIX Collaboration], Pys. Rev. Lett., 33 (8 [arxiv:8.4 [nucl-ex]]. 5. G. Agakisiev et al. [SAR Collaboration], Pys. Rev. Lett. 8, 73 ( [arxiv:.579 [nucl-ex]]. 6. A. Adare et al. [PHENIX Collaboration], Pys. Rev. C 88, no., 496 (3 [arxiv:34.34 [nucl-ex]]. 7. B. B. Abelev et al. [ALICE Collaboration], Pys. Lett. B 736, 96 (4 [arxiv:4.5 [nucl-ex]]. 8. S. Acarya et al. [ALICE Collaboration], Eur. Pys. J. C 77, no. 5, 339 (7 [Eur. Pys. J. C 77, no. 9, 586 (7] [arxiv:7.97 [e-ex]]. 9. K. M. Burke et al. [JE Collaboration], Pys. Rev. C 9, no., 499 (4 [arxiv:3.53 [nucl-t]].. A. Adare et al. [PHENIX Collaboration], Pys. Rev. C 8, 9 ( [arxiv:5.496 [nucl-ex]].. A. Adare et al. [PHENIX Collaboration], Pys. Rev. C 83, 499 ( [arxiv:4.353 [nucl-ex]].. J. Adam et al. [ALICE Collaboration], Pys. Rev. C 93, no. 3, 3493 (6 [arxiv: [nucl-ex]]. 3. S. Acarya et al. [ALICE Collaboration], Pys. Rev. C 98, no. 4, 449 (8 [arxiv: [nucl-ex]]. 4. A. Adare et al. [PHENIX Collaboration], Pys. Rev. C 84, 449 ( [arxiv: [nucl-ex]].

9 Guo-Yang Ma et al.: ω and K S Productions wit a Global Extraction of te Jet ransort Parameter 9 5. J. Adam et al. [ALICE Collaboration], Pys. Rev. C 95, no. 6, 6466 (7 [arxiv:7.555 [nucl-ex]]. 6. B. I. Abelev et al. [SAR Collaboration], Pys. Rev. Lett. 97, 53 (6 [nucl-ex/663]. 7. Y. Xu, J. Pys. G 37, 9459 ( [arxiv:.38 [nucl-ex]]. 8. W. Liu, C. M. Ko and B. W. Zang, Pys. Rev. C 75, 59 (7. 9. S.J. BrodskyandS.Gardner, Pys.Lett. B643, (6 [e-/689].. X. Cen, H. Zang, B. W. Zang and E. Wang, J. Pys. 37, 54 (.. X. Cen, C. Greiner, E. Wang, X. N. Wang and Z. Xu, Pys. Rev. C 8, 6498 (.. X. Cen,. Hirano, E. Wang, X. N. Wang and H. Zang, Pys. Rev. C 84, 349 (. 3. Z.Q. Liu,H.Zang, B. W. ZangandE. Wang, Eur.Pys. J. C 76, no., (6 [arxiv:56.84 [nucl-t]]. 4. W. Dai, X. F. Cen, B. W. Zang and E. Wang, Pys. Lett. B 75, 39 (5 [arxiv: [nucl-t]]. 5. W. Dai, B. W. Zang and E. Wang, Pys. Rev. C 98, 49 (8. 6. W. Dai, B. W. Zang, H. Z. Zang, E. Wang and X. F. Cen, Eur. Pys. J. C 77, no. 8, 57 (7 [arxiv:7.64 [nucl-t]]. 7. B. W. Zang and X. N. Wang, Nucl. Pys. A 7, 49 (3 [e-/395]. 8. B. W. Zang, E. Wang and X. N. Wang, Pys. Rev. Lett. 93, 73 (4 [nucl-t/394]. 9. A. Scafer, X. N. Wang and B. W. Zang, Nucl. Pys. A 793, 8 (7 [arxiv:74.6 [e-]]. 3.. Hirano, Pys. Rev. C 65, 9 ( [nucl-t/84]. 3.. Hirano and K. suda, Pys. Rev. C 66, 5495 ( [nucl-t/543]. 3. C. Sen, Z. Qiu, H. Song, J. Bernard, S. Bass and U. Heinz, Comut. Pys. Commun. 99, 6 ( D. Peresunko [ALICE Collaboration], Nucl. Pys. A 94-95, 755c (3 [arxiv:.5749 [nucl-ex]]. 34. B. I. Abelev et al. [SAR Collaboration], Pys. Rev. C 75, 649 (7 [nucl-ex/6733]. 35. N. Kidonakis and J. F. Owens, Pys. Rev. D 63, 549 ( [e-/768] J. Hou et al., JHEP 73, 99 (7 [arxiv: [e-]]. 37. S. Albino, B. A. Kniel and G. Kramer, Nucl. Pys. B 83, 4 (8 [arxiv: [e-]]. 38. D. Indumati and H. Saveeta, Int. J. Mod. Pys. A 7, 53 ( [arxiv:.5594 [e-]]. 39. H. Saveeta, D. Indumati and S. Mitra, Int. J. Mod. Pys. A 9, no. 7, 4549 (4 [arxiv:39.34 [e]]. 4. M. Hirai and S. Kumano, Comut. Pys. Commun. 83, ( [arxiv:6.553 [e-]]. 4. X. F. Guo and X. N. Wang, Pys. Rev. Lett. 85, 359 ( [e-/544]. 4. X. N. Wang and X. f. Guo, Nucl. Pys. A 696, 788 ( [e-/3]. 43. D. G. d Enterria, nucl-ex/ K. J. Eskola, P. Paakkinen, H. Paukkunen and C. A. Salgado, Eur. Pys. J. C 77, no. 3, 63 (7 [arxiv:6.574 [e-]]. 45. S. Y.Cen, K.M. Sen,W.Dai, B. W. Zang,H. Z. Zang and E. K. Wang, Commun. eor. Pys. 64, no., 95 ( B. W. Zang, G. Y. Ma, W. Dai, S. Wang and S. L. Zang, arxiv:8.93 [nucl-t].

10