PLB31416S0370-2693(15)00794-710.1016/j.physletb.2015.10.037The AuthorsExperimentsFig. 1(Color online.) Two-dimensional Δϕ vs. Δη correlation functions for charged hadron (left), pion (middle), and non-pion (right) triggers from 0–10% most-central Au+Au (top row) and minimum-bias d+Au (bottom row) data at 200 GeV. All trigger and associated charged hadrons are selected in the respective pT ranges 4<pTtrig<5  GeV/c and 1.5<pTassoc<4  GeV/c.Fig. 2(Color online.) The Δϕ and Δη projections of the pure-cone correlations for |Δη|<0.78 and |Δϕ|<π/4, respectively, for pion triggers (left two panels) and non-pion triggers (right two panels). Filled symbols show data from the 0–10% most-central Au+Au collisions at 200 GeV; open symbols show data from minimum-bias d+Au data at the same energy. Shaded boxes centered at zero show the uncertainty in background level determination; colored bands show the remaining systematic uncertainties.Fig. 3(Color online.) (a) Δϕ projections over |Δη|<1.5 in 0–10% Au+Au and minimum-bias d+Au data at 200 GeV after subtracting the jet-like fit components. Overlapping dashed and dotted lines illustrate the results of the two-dimensional fits for the flow- and mini-jet based models, respectively. Only statistical errors are shown. (b) Solid symbols show extracted Fourier coefficients from fitting in (π2<Δϕ<3π2)×(|Δη|<1.5) (1.52 for d+Au), for pion, non-pion, and charged hadron triggers in the flow-based model. Open symbols show d+Au data, scaled by the ratio of background levels. (c) V3/V2 ratio for pion and non-pion triggers in Au+Au, and the extrapolated value for “pure protons”, as described in the text. In all panels, statistical errors are shown as lines (smaller than symbol size for some points), and systematic uncertainties as colored boxes in panels (b) and (c).Table 1Fiducial ([|Δη|<0.78]×[|Δϕ|<π/4]) and extrapolated pure-cone yields for pion, non-pion and charged hadron (unidentified) triggers (see text), and the associated yield ratios.TriggerAu+Au 0–10%d+Au MinBiasFid.Ext.Stat.Sys.Fid.Ext.Stat.Sys.π0.2110.2143%7%0.1710.1714%6%non-π0.1360.1425%6%0.1420.1487%8%All0.1760.1802%5%0.1610.1682%5%Y(non-π)Y(π)0.6430.6626%5%0.8350.8668%8%Table 2First and second harmonic extracted using the mini-jet model in 0–10% most central Au+Au data at 200 GeV. Note that the amplitudes are multiplied by 100.Trigger100(V1±stat.±sys.)100(V2±stat.±sys.)π−0.86±0.09±0.07−0.017±0.045±0.034non-π−1.53±0.13±0.08−0.343±0.059±0.041all−1.19±0.05±0.01−0.173±0.025±0.004Di-hadron correlations with identified leading hadrons in 200 GeV Au+Au and d+Au collisions at STARSTAR CollaborationL.AdamczykaJ.K.AdkinssG.AgakishievqM.M.AggarwalacZ.AhammedauI.AlekseevoA.AparinqD.ArkhipkincE.C.AschenauercG.S.AverichevqX.BaihV.BairathizA.BanerjeeauR.BellwiedaqA.BhasinazA.K.BhatiacP.BhattaraiapJ.BielcikjJ.BielcikovakL.C.BlandcI.G.BordyuzhinoJ.BouchetrD.BrandenburgaiA.V.BrandinyI.BunzarovqJ.ButterworthaiH.CainesayM.Calderón de la Barca SánchezeJ.M.CampbellaaD.CebraeM.C.CervantesaoI.ChakaberiacP.ChaloupkajZ.ChangaoS.ChattopadhyayauX.ChenuJ.H.ChenalJ.ChengarM.CherneyiW.ChristiecG.ContinvH.J.CrawforddS.DaslL.C.De SilvaiR.R.DebbecT.G.DedovichqJ.DengakA.A.DerevschikovaeB.di RuzzacL.DidenkocC.DilksadX.DongvJ.L.DrachenbergatJ.E.DrapereC.M.DuuL.E.DunkelbergerfJ.C.DunlopcL.G.EfimovqJ.EngelagedG.EppleyaiR.EshafO.EvdokimovhO.EysercR.FatemisS.FaziocP.FederickJ.FedorisinqZ.FenggP.FilipqY.FisyakcC.E.FloreseL.FulekaC.A.GagliardiaoD.GarandafF.GeurtsaiA.GibsonatM.GirardavL.GreinervD.GrosnickatD.S.GunarathneanY.GuoajS.GuptaazA.GuptaazW.GuryncA.HamadrA.HamedaoR.HaquezJ.W.HarrisayL.HeafS.HeppelmanncS.HeppelmannadA.HirschafG.W.HoffmannapD.J.HofmanhS.HorvatayT.HuangcB.HuanghH.Z.HuangfX.HuangarP.HuckgT.J.HumanicaaG.IgofW.W.JacobsnH.JangtJ.JiacK.JiangajE.G.JudddS.KabanarD.KalinkinoK.KangarK.Kauderaw⁎H.W.KecD.KeanerA.KechechyanqZ.H.KhanhD.P.KikołaavA.KisielavL.KochendayD.D.KoetkeatL.K.KosarzewskiavA.F.KraishananP.KravtsovyK.KruegerbL.KumaracM.A.C.LamontcJ.M.LandgrafcK.D.LandryfJ.LauretcA.LebedevcR.LednickyqJ.H.LeecX.LianW.LialZ.M.LigY.LiarC.LiajX.LiajM.A.LisaaaF.LiugT.LjubiciccW.J.LlopeawM.LomnitzrR.S.LongacrecX.LuogG.L.MaalY.G.MaalR.MacL.MaalN.MagdyamR.MajkaayA.ManionvS.MargetisrC.MarkertapH.MasuivH.S.MatisvD.McDonaldaqK.MeehaneJ.C.MeiakN.G.MinaevaeS.MioduszewskiaoD.MishrazB.MohantyzM.M.MondalaoD.A.MorozovaeM.K.MustafavB.K.NandimMd.NasimfT.K.NayakauG.NigmatkulovyT.NiidaawL.V.NogachaeS.Y.NohtJ.NovakxS.B.NurushevaeG.OdyniecvA.OgawacK.OhagV.OkorokovyD.OlvittJr.anB.S.PagecR.PakcY.X.PanfY.PandithY.PanebratsevqB.PawlikabH.PeigC.PerkinsdA.PetersonaaP.PilecJ.PlutaavK.PoniatowskaavJ.PortervM.PosikanA.M.PoskanzervN.K.PruthiacJ.PutschkeawH.QiuvA.QuinterorS.RamachandransR.RaniwalaahS.RaniwalaahR.L.RayapH.G.RittervJ.B.RobertsaiO.V.RogachevskiyqJ.L.RomeroeA.RoyauL.RuancJ.RusnakkO.RusnakovajN.R.SahooaoP.K.SahulI.SakrejdavS.SalurvJ.SandweissayA.SarkarmJ.SchambachapR.P.ScharenbergafA.M.SchmahvW.B.SchmidkecN.SchmitzwJ.SegeriP.SeybothwN.ShahalE.ShahalievqP.V.ShanmuganathanrM.ShaoajB.SharmaacM.K.SharmaazW.Q.ShenalS.S.ShigQ.Y.ShoualE.P.SichtermannvR.SikoraaM.SimkokS.SingharM.J.SkobynN.SmirnovayD.SmirnovcL.SongaqP.SorensencH.M.SpinkabB.SrivastavaafT.D.S.StanislausatM.StepanovafM.StrikhanovyB.StringfellowafM.SumberakB.SummaadY.SunajZ.SunuX.M.SungX.SunvB.SurrowanD.N.SviridaoM.A.SzelezniakvA.H.TangcZ.TangajT.TarnowskyxA.TawfikaxJ.H.ThomasvA.R.TimminsaqD.TlustyaiT.TodorokicM.TokarevqS.TrentalangefR.E.TribbleaoP.TribedycS.K.TripathylO.D.TsaifT.UllrichcD.G.UnderwoodbI.UpsalaaG.Van BurencG.van NieuwenhuizencM.VandenbrouckeanR.VarmamA.N.VasilievaeR.VertesikF.VidebækcY.P.ViyogiauS.VokalqS.A.VoloshinawA.VossennJ.S.WanguF.WangafH.WangcG.WangfY.WangarY.WanggG.WebbcJ.C.WebbcL.WenfG.D.WestfallxH.WiemanvS.W.WissinknR.WittasY.F.WugWurZ.G.XiaoarW.XieafK.XinaiH.XuuZ.XucQ.H.XuakY.F.XualN.XuvS.YangajY.YangcQ.YangajY.YanguC.YangajY.YanggZ.YehZ.YehP.YepesaiL.YiayK.YipcI.-K.YooagN.YugH.ZbroszczykavW.ZhaajY.ZhangajZ.ZhangalJ.B.ZhanggJ.ZhangakX.P.ZhangarS.ZhangalJ.ZhanguJ.ZhaogC.ZhongalL.ZhouajX.ZhuarY.ZoulkarneevaqaAGH University of Science and Technology, Cracow 30-059, PolandAGH University of Science and TechnologyCracow30-059PolandbArgonne National Laboratory, Argonne, IL 60439, USAArgonne National LaboratoryArgonneIL60439USAcBrookhaven National Laboratory, Upton, NY 11973, USABrookhaven National LaboratoryUptonNY11973USAdUniversity of California, Berkeley, CA 94720, USAUniversity of CaliforniaBerkeleyCA94720USAeUniversity of California, Davis, CA 95616, USAUniversity of CaliforniaDavisCA95616USAfUniversity of California, Los Angeles, CA 90095, USAUniversity of CaliforniaLos AngelesCA90095USAgCentral China Normal University (HZNU), Wuhan 430079, ChinaCentral China Normal University (HZNU)Wuhan430079ChinahUniversity of Illinois at Chicago, Chicago, IL 60607, USAUniversity of Illinois at ChicagoChicagoIL60607USAiCreighton University, Omaha, NE 68178, USACreighton UniversityOmahaNE68178USAjCzech Technical University in Prague, FNSPE, Prague, 115 19, Czech RepublicCzech Technical University in PragueFNSPEPrague115 19Czech RepublickNuclear Physics Institute AS CR, 250 68 Řež/Prague, Czech RepublicNuclear Physics Institute AS CRŘež/Prague250 68Czech RepubliclInstitute of Physics, Bhubaneswar 751005, IndiaInstitute of PhysicsBhubaneswar751005IndiamIndian Institute of Technology, Mumbai 400076, IndiaIndian Institute of TechnologyMumbai400076IndianIndiana University, Bloomington, IN 47408, USAIndiana UniversityBloomingtonIN47408USAoAlikhanov Institute for Theoretical and Experimental Physics, Moscow 117218, RussiaAlikhanov Institute for Theoretical and Experimental PhysicsMoscow117218RussiapUnivenrsity of Jammu, Jammu 180001, IndiaUnivenrsity of JammuJammu180001IndiaqJoint Institute for Nuclear Research, Dubna, 141 980, RussiaJoint Institute for Nuclear Research141 980DubnaRussiarKent State University, Kent, OH 44242, USAKent State UniversityKentOH44242USAsUniversity of Kentucky, Lexington, KY, 40506-0055, USAUniversity of KentuckyLexingtonKY40506-0055USAtKorea Institute of Science and Technology Information, Daejeon 305-701, Republic of KoreaKorea Institute of Science and Technology InformationDaejeon305-701Republic of KoreauInstitute of Modern Physics, Lanzhou 730000, ChinaInstitute of Modern PhysicsLanzhou730000ChinavLawrence Berkeley National Laboratory, Berkeley, CA 94720, USALawrence Berkeley National LaboratoryBerkeleyCA94720USAwMax-Planck-Institut fur Physik, Munich 80805, GermanyMax-Planck-Institut fur PhysikMunich80805GermanyxMichigan State University, East Lansing, MI 48824, USAMichigan State UniversityEast LansingMI48824USAyMoscow Engineering Physics Institute, Moscow 115409, RussiaMoscow Engineering Physics InstituteMoscow115409RussiazNational Institute of Science Education and Research, Bhubaneswar 751005, IndiaNational Institute of Science Education and ResearchBhubaneswar751005IndiaaaOhio State University, Columbus, OH 43210, USAOhio State UniversityColumbusOH43210USAabInstitute of Nuclear Physics PAN, Cracow 31-342, PolandInstitute of Nuclear Physics PANCracow31-342PolandacPanjab University, Chandigarh 160014, IndiaPanjab UniversityChandigarh160014IndiaadPennsylvania State University, University Park, PA 16802, USAPennsylvania State UniversityUniversity ParkPA16802USAaeInstitute of High Energy Physics, Protvino 142281, RussiaInstitute of High Energy PhysicsProtvino142281RussiaafPurdue University, West Lafayette, IN 47907, USAPurdue UniversityWest LafayetteIN47907USAagPusan National University, Pusan 609735, Republic of KoreaPusan National UniversityPusan609735Republic of KoreaahUniversity of Rajasthan, Jaipur 302004, IndiaUniversity of RajasthanJaipur302004IndiaaiRice University, Houston, TX 77251, USARice UniversityHoustonTX77251USAajUniversity of Science and Technology of China, Hefei 230026, ChinaUniversity of Science and Technology of ChinaHefei230026ChinaakShandong University, Jinan, Shandong 250100, ChinaShandong UniversityJinanShandong250100ChinaalShanghai Institute of Applied Physics, Shanghai 201800, ChinaShanghai Institute of Applied PhysicsShanghai201800ChinaamState University Of New York, Stony Brook, NY 11794, USAState University Of New YorkStony BrookNY11794USAanTemple University, Philadelphia, PA 19122, USATemple UniversityPhiladelphiaPA19122USAaoTexas A&M University, College Station, TX 77843, USATexas A&M UniversityCollege StationTX77843USAapUniversity of Texas, Austin, TX 78712, USAUniversity of TexasAustinTX78712USAaqUniversity of Houston, Houston, TX 77204, USAUniversity of HoustonHoustonTX77204USAarTsinghua University, Beijing 100084, ChinaTsinghua UniversityBeijing100084ChinaasUnited States Naval Academy, Annapolis, MD, 21402, USAUnited States Naval AcademyAnnapolisMD21402USAatValparaiso University, Valparaiso, IN 46383, USAValparaiso UniversityValparaisoIN46383USAauVariable Energy Cyclotron Centre, Kolkata 700064, IndiaVariable Energy Cyclotron CentreKolkata700064IndiaavWarsaw University of Technology, Warsaw 00-661, PolandWarsaw University of TechnologyWarsaw00-661PolandawWayne State University, Detroit, MI 48201, USAWayne State UniversityDetroitMI48201USAaxWorld Laboratory for Cosmology and Particle Physics (WLCAPP), Cairo 11571, EgyptWorld Laboratory for Cosmology and Particle Physics (WLCAPP)Cairo11571EgyptayYale University, New Haven, CT 06520, USAYale UniversityNew HavenCT06520USAazUniversity of Jammu, Jammu 180001, IndiaUniversity of JammuJammu180001India⁎Corresponding author.Editor: D.F. GeesamanAbstractThe STAR Collaboration presents for the first time two-dimensional di-hadron correlations with identified leading hadrons in 200 GeV central Au+Au and minimum-bias d+Au collisions to explore hadronization mechanisms in the quark gluon plasma. The enhancement of the jet-like yield for leading pions in Au+Au data with respect to the d+Au reference and the absence of such an enhancement for leading non-pions (protons and kaons) are discussed within the context of a quark recombination scenario. The correlated yield at large angles, specifically in the ridge region, is found to be significantly higher for leading non-pions than pions. The consistencies of the constituent quark scaling, azimuthal harmonic model and a mini-jet modification model description of the data are tested, providing further constraints on hadronization.Experimental data from heavy-ion collisions at ultra-relativistic energies achieved at the Relativistic Heavy Ion Collider (RHIC), and more recently at the Large Hadron Collider (LHC), are conventionally interpreted in terms of a unique form of matter, the strongly-interacting Quark Gluon Plasma (sQGP). It is estimated that temperatures reached in those collisions [1,2] are well above the critical values predicted by lattice quantum chromodynamics calculations for the phase transition between hadronic and de-confined (partonic) matter [3]. The RHIC experiments concluded that the formed medium displays the properties of a nearly perfect liquid [4]. A distinct feature of the sQGP is jet quenching, which describes the large energy loss of “hard” (high transverse momentum, pT) probes observed for example in measurements of inclusive hadron distributions [5].Jet quenching is also evident in modifications of back-to-back di-hadron correlations with a leading (high-pT) “trigger” hadron in Au+Au collisions in comparison to p+p and d+Au data [6–10]. One of the striking features found in di-hadron correlations from heavy-ion collisions is the emergence of a long-range plateau in relative pseudorapidity (Δη) on the near-side of a trigger hadron (small relative azimuth Δϕ), referred to as the “ridge” [6–8]. The majority of recent theoretical descriptions of this phenomenon invoke a transport of initial-state to final-state anisotropy via hydrodynamic expansion, thus connecting measured observables to transport coefficients and properties of the medium [11–13]. Most of the proposed alternative explanations also require hydrodynamic evolution of the medium to reproduce the ridge [14]. The latest observations of ridge-like correlations in high-multiplicity p+p and p+Pb collisions at the LHC provide new tests of theoretical explanations of the ridge [15,16].Another anomaly, the enhancement of the relative baryon-to-meson production, was discovered at RHIC in the intermediate-pT range between 2 and 5 GeV/c, where the ridge happens to be most prominent [17–20]. The ratio of proton to pion yields in central Au+Au collisions exceeds by more than a factor of two that in d+Au and p+p events. Similar baryon enhancements were reported in the strange-hadron sector [21,22]. In the same kinematic region, baryons and mesons exhibit different trends in azimuthal anisotropy, which at RHIC appear to scale with the number of constituent quarks [23]. Recombination models, which incorporate the coalescence of two or three constituent quarks as a formation mechanism for mesons and baryons, are able to reproduce the observed enhancements in inclusive measurements [24]. Description of hadronization processes remains challenging for theoretical calculations (see for example the unexpected measurements reported in Ref. [25]); we expect these new measurements will facilitate further developments in this area.In this Letter, we use angular correlations of intermediate-pT particles with identified leading hadrons to further explore possible hadronization mechanisms in the quark gluon plasma, including changes to parton fragmentation patterns, dilution effects (reduction in per-trigger yields) due to recombination contributions, and quark number scaling behavior in correlations at large relative angles.Two-dimensional di-hadron correlations in Δϕ and Δη, with statistically separated pion and non-pion triggers, are studied for the 0–10% most-central Au+Au and minimum-bias d+Au collisions at the center-of-mass energy per nucleon pair sNN=200 GeV. No charge separation is considered; in this paper, the terms pion, proton, and kaon will be used to refer to the sum of particles and their respective anti-particles. We separately study the short-range (jet-like peak) and long-range (ridge) correlations in central Au+Au data, comparing to reference measurements performed in d+Au collisions. Details of the short-range correlations can shed light on the interplay between parton fragmentation, energy loss, and recombination processes in the quark gluon plasma. The long-range correlations are studied with two approaches: (1) Fourier decomposition where we extract the azimuthal harmonic amplitudes which in some approaches are interpreted as hydrodynamic “flows” [12], and (2) a mini-jet (defined in [26]) modification model [26,27].The analysis was conducted using 1.52×108 central-triggered Au+Au events at sNN=200 GeV from STAR's 2010 data run, and 4.6×107 events from the 2008 minimum-bias 200 GeV d+Au data set. Particle densities as well as Glauber Monte Carlo results for these centrality selections can be found in Ref. [28]. The STAR Time Projection Chamber (TPC) [29] was used for tracking, momentum reconstruction and particle identification. Contamination by tracks from another collision (“pileup”), which can distort the shape of di-hadron correlations [26], was removed by rejecting events with an abnormally large (over three standard deviations above the average) number of tracks not originating from the primary vertex.Trigger particles are defined as the highest-pT charged hadron in a given event with pTtrig between 4 and 5 GeV/c; charged hadrons between 1.5 and 4 GeV/c are associated with each trigger particle. This kinematic range focuses on a particularly interesting region where trigger particle production is thought to be dominated by fragmentation, at least in pp and d+Au collisions. For the same range in Au+Au events, the baryon-to-meson enhancement is large, suggesting significant recombination contributions [30]. Since the medium induced jet quenching affects the correlations in essentially an opposite way from thermal parton recombination contributions, comparing the correlations for proton and pion triggers provides an additional handle for separating these effects. Also, the ridge and away-side modifications in this pT range are significant, and elliptic flow, and its pT dependence, are minimal, thus facilitating the present tests of constituent quark number scaling. All trigger particles are required to have at least 30 TPC points per track (for optimal identification), otherwise, standard quality cuts and corrections are applied as described in Ref. [8]. After quality cuts, 3.5×106 Au+Au and 1×105 d+Au events with a trigger particle were used for the analysis. The statistical hadron identification procedure relies on the measured ionization energy loss (dE/dx) in the TPC gas. The dE/dx calibration was carried out individually for five pseudorapidity and two trigger pT bins. Details of the particle identification (PID) technique are identical to those in Refs. [17,28,31].We construct a two-dimensional correlation with each trigger and all associated hadrons in an event, following the procedure outlined in Ref. [8]. Pion identification is straightforward: selecting triggers with dE/dx above the central (expected) pion value provides a sample with 98% pion purity and, by construction, 50% selection efficiency. The “pure-pion” correlation is constructed with those triggers. The remaining triggers are comprised of all protons, about 97% of all kaons, and the remaining 50% of pions.We remove the pion contribution from the correlation with those remaining triggers by direct subtraction of the pure-pion-triggered correlation. The resulting “non-pion” correlation is then associated with a mixture of proton and kaon triggers (about three protons for every two kaons [32]). Separating kaons from protons is complicated by the small dE/dx difference between the two and was not attempted in this Letter. The systematic uncertainty due to the pion subtraction procedure is included in the PID uncertainty. The evaluation procedure is similar to previous identified particle analyses with the STAR Time Projection Chamber, where sensitivities of the final observables to systematic variations in the dE/dx cut parameters were determined. The feed-down contribution from weak decay daughters to the trigger particles cannot be disentangled directly. Due to decay kinematics, the dominant feed-down contribution originates from Λ→pπ and is estimated to constitute about 5% of the non-pion triggers. Resonance contributions are greatly suppressed at high pT and shown to give only minor contributions to correlation structures [26]. All raw correlation functions are corrected for detector inefficiency derived from Monte Carlo tracks embedded into real data as in Refs. [8,9,28]. Pair-acceptance effects are corrected using the mixed-event technique as in Ref. [9]. The resulting correlations are shown in Fig. 1, with visible differences between the two trigger types in both jet-like peak and large Δη region in Au+Au. A significantly larger ridge amplitude is seen for non-pion triggers, while the jet-like peak is more pronounced for the pion triggers. By comparison, the correlations in d+Au show no discernible ridge on the near-side, while differences between trigger types in the jet-like region are qualitatively similar to Au+Au, suggesting that these may be partly due to kinematic effects. In the following, we analyze these modifications individually.Initially, we study the small-angle jet-like correlated signals. Assuming that all background contributions are Δη-independent, as shown in Refs. [8,33], we subtract those contributions averaged over large |Δη|=0.9–1.5 from the full correlations, resulting in “pure-cone” distributions. This procedure is supported by the two-dimensional fits to the data described below. We then calculate the fiducial jet-like yield in |Δη|<0.78, |Δϕ|<π/4 as in Ref. [7]. To isolate medium effects from initial-state nuclear effects, the Au+Au results are then compared to the correlation function constructed in an identical way for d+Au data (see Fig. 2). We report significant differences in the jet-like yield per trigger between the two systems for pion triggers. At the same time, correlations with non-pion triggers show, within uncertainties, similar yields for the two systems. For quantitative comparisons, the integrated yields are presented in Table 1. The yield extrapolation outside the fiducial range is performed using cone-shape modeling described below. The systematic errors are dominated by the tracking efficiency uncertainty (5%); other sources include uncertainties from pT resolution (3%), PID uncertainty (2–3%), background level determination (2% for Au+Au, 2–5% for d+Au; found by varying the range for the Δη-independent ridge structure between |Δη|=0.8–1.4 and |Δη|=1.0–1.6), track splitting/merging correction (1%), and pair acceptance (<1%). The effect of feed-down protons on the jet-like yield is estimated to be less than 1%.The jet-like yield in the pT range 1.5-4 GeV/c associated with pion triggers in central Au+Au collisions is enhanced by 24±6(stat.)±11(sys.)% with respect to the reference measurement in d+Au. The yields for non-pion triggers are found to be similar between the two systems. A previous work found similar trends in near-side associated yields [34]; however, in that one-dimensional analysis, no separation between jet-like peak and ridge contributions was possible. We find that the jet-like yield for unidentified charged hadron triggers is also enhanced, consistent with our identified trigger results. The enhancement of the jet-like yield of soft hadrons associated with pion triggers could be caused by the jet-quenching effect and/or medium-induced modification of fragmentation functions, and is qualitatively consistent with other observations from non-identified correlations [35] and direct jet measurements [36–38] for low pT hadrons. It is expected that a larger fraction of non-pion triggers are produced from gluon-jets rather than quark-jets compared to pion triggers [32,39,40]. A predicted higher energy loss for in-medium gluons should then result in even larger jet-like yields for non-pion triggers [41]. On the other hand, particle production from recombination should produce smaller yields than particle production from hard processes (fragmentation) [30], thus diluting (reducing) per-trigger associated yields. This dilution effect would be stronger for baryons, as more intermediate-pT baryons than mesons are expected to be formed through such a mechanism. The associated yields for non-pion triggers combine both of these competing effects. Thus the observed reduction could be due to a larger recombination effect relative to that from the increased energy loss expected for non-pion leading particles.The ratio of associated yields for non-pion and pion leading hadrons is shown in Table 1 for Au+Au and d+Au systems. In these ratios, dominant contributions to systematic uncertainties from the tracking efficiency estimate cancel out. The double-ratio constructed from these two results quantifies the relative decrease in associated jet-like yields for non-pion triggers with respect to leading pion results in Au+Au compared to d+Au. This double-ratio, 0.76±0.08(stat.)±0.07(sys.), can measure the net effect of the competition between higher energy loss/higher associated yields at lower pT for gluon jets versus reduced yields due to recombination in central Au+Au collisions. Currently, no quantitative predictions for either of these two mechanisms are available for direct comparison with data.Outside of the jet-like cone region, we find no Δη-dependence in the correlated yields within our fiducial range. To characterize the long-range contributions in Au+Au, we perform two-dimensional fits to the full correlation with two different models. One model attributes the ridge to modified fragmentation of produced mini-jets, and the other explains it in terms of higher-order hydrodynamic flows. In both models, the near-side jet-like peak is mathematically characterized by a two-dimensional generalized Gaussian e−(|Δϕ|/αϕ)βϕe−(|Δη|/αη)βη. The resulting fit parameters for the jet cone are found to be identical between the two models and were used for extrapolation of the jet-like cone yields presented in Table 1.The Δϕ projections of the pseudorapidity-independent parts of the two-dimensional correlations (after subtracting the jet-like peak), are shown in Fig. 3, panel (a), together with both fit functions discussed below. Also shown is the corresponding projection for d+Au data, shifted by an arbitrary offset. We performed the Fourier analysis on this data as well; without an appreciable near-side ridge, the above mini-jet model is not applicable here.In the flow-based approach, based on a hydrodynamic expansion of an anisotropic medium, all Δη-independent parts of the correlations are described via Fourier expansion: A(1+∑n=1N2Vncos⁡nΔϕ), where A describes the magnitude of the uncorrelated background, V2 is conventionally associated with “elliptic flow”, and V3 with “triangular flow”. In this work, the first five terms (N=5) exhaust all features of the correlation to the level of statistical uncertainty, and Vn represents the combined trigger and associated hadron anisotropy parameters. We note that in this approach, the fragmentation contributions to the away-side correlations in central Au+Au data are expected to be strongly suppressed relative to flow effects by quenching, and they are therefore neglected [11,12]. In the d+Au data, on the contrary, the away-side jet contributions dominate and no appreciable near-side correlated yield at large Δη is present.The Fourier fit results are shown in Fig. 3 (b). In Au+Au, the second harmonic is dominant in all long-range correlations for the central data, followed by the triangular (V3) term. Higher-order harmonic amplitudes rapidly decrease. All harmonic amplitudes for non-pion triggers are found to be larger than those for pion triggers, which is qualitatively consistent with recombination expectations.The corresponding Fourier harmonics in d+Au are found to be consistent with expectations for a decomposition of a Gaussian peak at π, i.e. rapidly falling and with alternating signs. For n=3, the harmonics are already consistent with zero within errors, confirming that the Fourier and mini-jet approaches are indistinguishable in this case.In central Au+Au collisions at RHIC, elliptic flow parameters of identified hadrons have been shown to scale with the number of constituent quarks nq, suggesting collective behavior at the partonic level [23]. The estimated baryon/meson ratio for V2 in this analysis is also consistent with 3/2, see below. We note that in our trigger pT range, azimuthal anisotropy is approximately independent of pT [42], eliminating the need to address quark momentum dependence. To test whether this scaling extends to the triangular term, we examine the V3/V2 ratios. This test assumes that the measured Fourier coefficients factorize into Vn=〈vntrig〉〈vnassc〉, where vntrig and vnassc measure azimuthal anisotropies of trigger and associated hadrons, respectively [12]. The factorization has been demonstrated experimentally for V2 in [43]. The extracted V2 coefficients are found consistent with the product of previously measured identified and unidentified v2 values. Since the selection of associated particles is identical for all correlations in this analysis, the anisotropy contributions from associated hadrons should cancel in the ratios of Vn coefficients. Fig. 3 (c) shows V3/V2 ratios extracted from long-range correlations versus average nq per particle for pion and non-pion triggers. The systematic uncertainty, determined by varying the fitting range and the dE/dx cut position for pion/non-pion separation, was found to be similar to, or smaller than, the statistical uncertainty. We find that the ratio of triangular and elliptic flow is 0.546±0.025(stat.)±0.018(sys.) for pion triggers and 0.681±0.025(stat.)±0.015(sys.) for non-pions. If the measured final-state azimuthal anisotropies are indeed of collective partonic origin which transform into final-state hadronic observables through the coalescence/recombination of constituent quarks, then we would expect the same dependence of all vntrig on constituent quark number. Even with the significant meson contribution to non-pion triggers, the ratios give a strong indication of a breaking of the simple nq scaling behavior between the second and third Fourier harmonics. Assuming that kaons, as mesons, adhere to the pion scaling trend, and using the known p/π ratio reported in Ref. [32], we construct an estimate of the V3/V2 ratio for pure protons in Fig. 3 (c). The systematic uncertainty in the estimated “pure-proton” V3/V2 value of 0.736±0.038(stat.)±0.032(sys.) includes an additional 1% uncertainty from PID. Feed-down protons from Λ closely preserve the original parent direction, and we expect no measurable effect on the Δη-independent terms, as Λ and protons have very similar azimuthal anisotropy in our kinematic range.The observed violation of constituent quark number scaling for V3, based on the “pure proton” extrapolated value for V3, V3(baryon)/V3(meson)=2.03±0.12(stat.)±0.20(sys.) compared to V2(baryon)/V2(meson)=1.50±0.06(stat.)±0.07(sys.), is intriguing because recombination/coalescence models are the only ones presently capable of describing constituent quark scaling behavior among the V2 parameters for many identified hadrons.The difference between the V3 and V2 scaling behavior demonstrated in Fig. 3(c) therefore suggests the need for other contributions to long-range correlations to explain the data. We note that deviations from nq scaling of elliptic flow have been observed at the LHC [44], and at RHIC for non-central collisions [45]. The vn scaling proposed in Ref. [46] better describes our measured V3/V2 ratios, but still under-predicts the enhancement for non-pion triggers.In the mini-jet model, in which the major component is in-medium modification of fragmentation, only the first two terms (N=2) of the Fourier expansion are kept and the near-side ridge in this analysis is modeled by a one-dimensional Gaussian, resulting in A(1+2V1cos⁡Δϕ+2V2cos⁡2Δϕ)+Be−Δϕ2/2σ2. Here A is the uncorrelated yield, B is the ridge amplitude, and σ is the ridge width parameter. The dipole V1 is designated to describe the away-side jet and/or momentum conservation effects, and V2 describes a non-jet quadrupole (potentially of flow origin). The addition of the 1D near-side Gaussian, which differs from the original model elements in Ref. [26], was necessary to reproduce the data, as noted in Ref. [27]. The mini-jet model fit describes the measured Au+Au correlations for all three trigger types as well as the flow-based approach (Fig. 3(a)), yielding identical uniformly distributed residuals and χ2 values. The extracted harmonic amplitudes are shown in Table. 2. As the away-side structure is for the most part described by the dipole term, the magnitude of the V1 amplitude is found to be significantly larger for leading non-pions than for pions. For back-to-back jets, this V1 increase is supposed to balance the near-side (leading) jet contributions, which would have to consist of both the jet-like peak and the ridge because the jet-like peak alone decreases for non-pion trigger particles. Understanding the behavior of the V2 term in the mini-jet model fits is challenging: the V2 amplitude, while consistent with zero for pion triggers, is significantly negative for non-pion triggers. This negative value for V2, which is conventionally associated with elliptic flow, is not expected from any known source and calls into question the applicability of the assumed parameterization for the centrality and pT range studied here, the validity of the “mini-jets + quadrupole only” physics scenario, or both.In summary, a statistical separation of pion and non-pion triggers was performed to study the systematic behavior of di-hadron correlations from central Au+Au and minimum-bias d+Au collisions at 200 GeV with the STAR experiment. The correlations, decomposed into short- and long-range parts in Δη, are analyzed for different identified trigger types to test the consistency of two models in order to improve our understanding of hadronization mechanisms in the quark gluon plasma. We find significant enhancement of intermediate-pT charged-hadron jet-like yields associated with pion triggers relative to a d+Au reference measurement. The enhancement is qualitatively consistent with observed modifications of jet fragmentation functions measured at the LHC, suggesting it results from the energy loss process. For the non-pion trigger sample, a larger contribution from gluon fragmentation is expected compared to pion triggers [32,39,40]. Due to the color-charge factor, a larger energy loss for gluons is expected relative to that of quarks. No enhancement is observed for non-pion triggers in contrast to pQCD-based expectations for color charge dependence of energy loss. This lack of enhancement may indicate a competition between parton–medium interaction effects and dilution of jet triggers by quark recombination contributions.No statistically significant ridge is found associated with either trigger type in minimum bias d+Au data. In Au+Au data, we find a significantly larger ridge-like yield and away-side correlation strength for non-pion than for pion triggers. Two fitting models which are mathematically similar but which are based on distinct physical assumptions were applied to the Au+Au data. Both models, while describing the correlations well, attain parameter values which are problematic within the assumed physical scenarios. In the flow model, the observed differences of V3/V2 ratios imply that the explanation of the ridge and away-side modifications as resulting only from hydrodynamic flow of a partonic medium with constituent quark recombination at hadronization is incomplete. On the other hand, the negative V2 result for the mini-jet based model for leading non-pions indicates that for the data reported here, either the assumed scenario or the mathematical parameterization for jets and dijets is inadequate, or both. These results may have significant implications for understanding the origin of the ridge and hadronization in the QGP.AcknowledgementsWe thank the RHIC Operations Group and RCF at BNL, the NERSC Center at LBNL, the KISTI Center in Korea, and the Open Science Grid consortium for providing resources and support. This work was supported in part by the Offices of Nuclear Physics within the U.S. DOE Office of Science, the U.S. NSF, the Ministry of Education and Science of the Russian Federation, NSFC, CAS, MOST and MOE of China, the National Research Foundation of Korea, GACR and MSMT of the Czech Republic, FIAS of Germany, DAE, DST, and UGC of India, the National Science Centre of Poland, National Research Foundation, the Ministry of Science, Education and Sports of the Republic of Croatia, and RosAtom of Russia.References[1]A.AdarePHENIX CollaborationPhys. Rev. Lett.1042010132301[2]M.WildeALICE CollaborationNucl. Phys. A904–9052013573c[3]F.KarschNucl. Phys. A6982002199[4]J.AdamsSTAR CollaborationNucl. Phys. A7572005102K.AdcoxPHENIX CollaborationNucl. Phys. A7572005184B.BackPHOBOS CollaborationNucl. Phys. A757200528I.ArseneBRAHMS CollaborationNucl. Phys. A75720051[5]J.AdamsSTAR CollaborationPhys. Rev. Lett.912003172302S.S.AdlerPHENIX CollaborationPhys. Rev. Lett.912003072303[6]J.AdamsSTAR CollaborationPhys. Rev. Lett.952005152301[7]B.I.AbelevSTAR CollaborationPhys. Rev. C802009064912[8]G.AgakishievSTAR CollaborationPhys. Rev. C852012014903[9]M.AggarwalSTAR CollaborationPhys. Rev. C822010024912[10]J.AdamsSTAR CollaborationPhys. Rev. Lett.912003072304A.AdarePHENIX CollaborationPhys. Rev. Lett.982007232302A.AdarePHENIX CollaborationPhys. Rev. C782008014901H.AgakishievSTAR CollaborationPhys. Rev. C832011061901B.I.AbelevSTAR CollaborationPhys. Rev. Lett.1052010022301[11]B.AlverG.RolandPhys. Rev. C812010054905B.AlverG.RolandPhys. Rev. C822010039903(Erratum)[12]M.LuzumPhys. Lett. B6962011499[13]C.GaleS.JeonB.SchenkeInt. J. Mod. Phys. A2820131340011[14]C.-Y.WongPhys. Rev. C842011024901R.C.HwaL.ZhuarXiv:1101.1334C.AndrésA.MoscosoC.PajaresarXiv:1405.3632G.MoschelliS.GavinNucl. Phys. A836201043H.PetersenV.BhattacharyaS.A.BassC.GreinerPhys. Rev. C842011054908X.ZhangJ.LiaoarXiv:1311.5463A.DumitruK.DuslingF.GelisJ.Jalilian-MarianT.LappiR.VenugopalanPhys. Lett. B697201121[15]V.KhachatryanCMS CollaborationJ. High Energy Phys.10092010091S.ChatrchyanCMS CollaborationPhys. Lett. B7182013G.AadATLAS CollaborationPhys. Lett. B725201360ATLAS Collaboration, ATLAS-CONF-2015-027.ATLAS Collaboration, ATLAS-COM-CONF-2015-033.LHCb Collaboration, LHCb-CONF-2015-004.LHCb Collaboration, CERN-LHCb-CONF-2015-004.ALICE CollaborationarXiv:1506.08032 [nucl-ex][16]K.WernerB.GuiotI.KarpenkoT.PierogProceedings of 24th International Conference on Ultra-Relativistic Nucleus–Nucleus Collisions (Quark Matter 2014)Nucl. Phys. A93120148391arXiv:1411.1048 [nucl-th]K.DuslingR.VenugopalanProceedings of 24th International Conference on Ultra-Relativistic Nucleus–Nucleus Collisions (Quark Matter 2014)Nucl. Phys. A9312014283I.BautistaA.F.TellezP.GhosharXiv:1505.00924 [hep-ph]T.AltinolukN.ArmestoG.BeufA.KovnerM.LublinskyarXiv:1503.07126 [hep-ph][17]B.I.AbelevSTAR CollaborationPhys. Rev. Lett.972006152301[18]K.AdcoxPHENIX CollaborationPhys. Rev. Lett.882002242301[19]B.AbelevSTAR CollaborationPhys. Lett. B6552007104[20]B.B.AbelevALICE CollaborationPhys. Lett. B7362014196[21]M.AggarwalSTAR CollaborationPhys. Rev. C832011024901[22]B.B.AbelevALICE CollaborationPhys. Rev. Lett.1112013222301[23]B.AbelevSTAR CollaborationPhys. Rev. C752007054906A.AdarePHENIX CollaborationPhys. Rev. Lett.982007162301[24]R.C.HwaC.B.YangPhys. Rev. C702004024905R.J.FriesB.MüllerC.NonakaS.A.BassPhys. Rev. C682003044902[25]A.AirapetianHERMES CollaborationNucl. Phys. B78020071A.AirapetianHERMES CollaborationEur. Phys. J. A472011113[26]G.AgakishievSTAR CollaborationPhys. Rev. C862012064902[27]R.L.RayD.J.PrindleT.A.TrainorarXiv:1308.4367[28]B.AbelevSTAR CollaborationPhys. Rev. C792009034909[29]M.AndersonNucl. Instrum. Methods Phys. Res., Sect. A, Accel. Spectrom. Detect. Assoc. Equip.4992003659[30]R.J.FriesNucl. Phys. A7832007125R.C.HwaarXiv:0904.2159 [nucl-th][31]M.ShaoNucl. Instrum. Methods Phys. Res., Sect. A, Accel. Spectrom. Detect. Assoc. Equip.5582006419[32]J.AdamsSTAR CollaborationarXiv:nucl-ex/0601042J.AdamsSTAR CollaborationPhys. Lett. B6372006161[33]B.AbelevPhys. Rev. C772008054901[34]S.AdlerPHENIX CollaborationPhys. Rev. C712005051902[35]K.AamodtALICE CollaborationPhys. Rev. Lett.1082012092301[36]L.AdamczykSTAR CollaborationPhys. Rev. Lett.1122014122301[37]S.ChatrchyanCMS CollaborationPhys. Rev. C902014024908arXiv:1406.0932 [nucl-ex][38]G.AadATLAS CollaborationPhys. Lett. B7392014320arXiv:1406.2979 [hep-ex][39]S.AlbinoB.KniehlG.KramerNucl. Phys. B7252005181[40]B.MohantySTAR CollaborationProceedings of the 23rd Winter Workshop on Nuclear Dynamics2007arXiv:0705.0953[41]X.-N.WangPhys. Rev. C5819982321Q.WangX.-N.WangPhys. Rev. C712005014903S.WicksW.HorowitzM.DjordjevicM.GyulassyNucl. Phys. A7842007426N.ArmestoA.DaineseC.A.SalgadoU.A.WiedemannPhys. Rev. D712005054027[42]J.AdamsSTAR CollaborationPhys. Rev. C722005014904[43]K.AamodtALICE CollaborationPhys. Lett. B7082012249[44]B.B.AbelevALICE CollaborationarXiv:1405.4632 [nucl-ex][45]A.AdarePHENIX CollaborationPhys. Rev. C852012064914[46]R.LaceyJ. Phys. G382011124048