Physics Letters B 777 (2018) 369–373 Contents lists available at ScienceDirect Physics Letters B www.elsevier.com/locate/physletb Event-plane dependent di-hadron correlations with harmonic vn subtraction in a hydrodynamic model Wagner M. Castilho a, Wei-Liang Qian b,a,∗, Yogiro Hama c, Takeshi Kodama d,e a Faculdade de Engenharia de Guaratinguetá, Universidade Estadual Paulista, 12516-410, Guaratinguetá, SP, Brazil b Escola de Engenharia de Lorena, Universidade de São Paulo, 12602-810, Lorena, SP, Brazil c Instituto de Física, Universidade de São Paulo, Caixa Postal 66318, 05389-970 São Paulo, SP, Brazil d Instituto de Física, Universidade Federal do Rio de Janeiro, Caixa Postal 68528, 21941-972 Rio de Janeiro, RJ, Brazil e Instituto de Física, Universidade Federal Fluminense, 24210-346, Niterói, RJ, Brazil a r t i c l e i n f o a b s t r a c t Article history: Received 28 July 2017 Received in revised form 18 December 2017 Accepted 22 December 2017 Available online 27 December 2017 Editor: J.-P. Blaizot In this work, a hydrodynamic study of the di-hadron azimuthal correlations for the Au+Au collisions at 200 GeV is carried out. The correlations are evaluated using the ZYAM method for the centrality windows as well as the transverse momentum range in accordance with the existing data. Event-plane dependence of the correlation is obtained after the subtraction of contributions from the most dominant harmonic coefficients. In particular, the contribution from the triangular flow, v3, is removed from the proper corre- lations following the procedure implemented by the STAR collaboration. The resultant structure observed in the correlations was sometimes attributed to the mini-jet dynamics, but the present calculations show that a pure hydrodynamic model gives a reasonable agreement with the main feature of the published data. A brief discussion on the physical content of the present findings is presented. © 2018 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3. 1. Introduction Measurements on the two-particle correlations in the relativis- tic heavy-ion collisions, expressed in terms of the pseudorapidity difference �η and the angular spacing �φ, were carried out by various experimental collaborations [1–7] at both RHIC and LHC. The shape of the two-particle correlations for different collision systems, at various transverse momentum range and its evolution as a function of trigger particle azimuthal angle have been both extensively studied. They are understood to provide relevant in- formation on the jets originating from small momentum transfer scatterings as well as the hot, dense medium created in the col- lisions [8,9]. The observed correlation yields are characterized by an enhancement on the near side around �φ ≈ 0, known as the “ridge”, which possesses a long �η extension in the longitudinal direction. Besides heavy-ion collisions, such ridge structures were also observed in pp [10] and pA [11–13] collisions at LHC. For pA collisions, it is found that the ridge yields vary with centrality. The measured jet-like yields obtained by subtracting long-range pseudorapidity correlations observed in high-multiplicity events, * Corresponding author. E-mail address: wlqian@usp.br (W.-L. Qian). https://doi.org/10.1016/j.physletb.2017.12.057 0370-2693/© 2018 The Authors. Published by Elsevier B.V. This is an open access article SCOAP3. on the other hand, are approximately constant for different cen- tralities [14]. This observed feature shows that the physics behind jet-like yield and ridge yield are indeed distinct, being the latter attributed to the collective flow of the system. The correlation on the away side is found to be more significant in AA collisions than in pp and pA, and it presents a double-peak structure, usually called “shoulders”, which evolves continuously from the double peak for central to one peak for peripheral collisions [2,4]. In order to interpret the “ridge” and “shoulders” in AA collisions by a uniform picture, we proposed the so-called peripheral-tube model [15–19]. In this model, the phenomenon is attributed to the local (nonlinear) behavior of hydrodynamics. The phenomenon can also be explained regarding the triangular flow, as usually done [20,21]. It is understood that hydrodynamical evolution transforms the spatial inhomogeneity of participating nucleons in the initial conditions into the momentum anisotropy of the observed hadrons [22–24]. The triangular flow in the one-particle distribution func- tion generates three peaks in the two-particle correlations: one peak on the near side at �φ = 0 and two others at �φ = 2π/3 and �φ = 4π/3 corresponding to the double peak on the away- side. It is noted that studies by using the AMPT model were carried out [25,26] which showed that the double peak disappears when the contributions due to the elliptic flow and triangular flow are subtracted. The above finding seems to indicate that the triangular under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by https://doi.org/10.1016/j.physletb.2017.12.057 http://www.ScienceDirect.com/ http://www.elsevier.com/locate/physletb http://creativecommons.org/licenses/by/4.0/ mailto:wlqian@usp.br https://doi.org/10.1016/j.physletb.2017.12.057 http://creativecommons.org/licenses/by/4.0/ http://crossmark.crossref.org/dialog/?doi=10.1016/j.physletb.2017.12.057&domain=pdf 370 W.M. Castilho et al. / Physics Letters B 777 (2018) 369–373 flow indeed plays an essential role in the observed structure on the away side. The two-particle correlations are also investigated as a function of the trigger angle φs , known as the event-plane de- pendence of the two-particle correlations [27,28]. It was found that the away side structure evolves from only one peak in the in-plane trigger direction with φs = 0 to double peaks at the out-of-plane trigger direction with φs = π/2. It is worth noting that the above data were extracted by mainly subtracting the background contri- butions of the elliptic flow. The observed features of the data can be understood in terms of a hydrodynamic interpretation known as peripheral one-tube model [29] where geometric fluctuations are manifested as high-energy tubes randomly distributed in the ini- tial conditions. The model [30] is also able to explain the observed centrality dependence of the away side structure in the correla- tions [4]. Based on analysis of the earlier STAR data [28], Luzum pointed out [31] that the observed two-particle correlations are consis- tent with being entirely generated by the collective flow. In the work, the author analyzed the elliptic flow coefficients v(a) 2 , v(t,R) 2 and showed that the flow background used in ZYAM subtraction was probably underestimated, and the non-flow contribution of the second Fourier coefficient is likely insignificant. However, it is not clear that the above reasoning also applies to the third Fourier coefficient, especially because no direct measurement of v3 were made in STAR’s analysis [28]. This is because V 3� was found to be weakly dependent on φs . Even if V 3� , the third Fourier coefficient of the proper correlations, and v(a) 3 v(t,R) 3 have roughly the same dependence on φs provided they are not identical, it is more likely that the subtraction was correct and the triangular non-flow sig- nal is mostly independent of the trigger angle. In other words, one does not know for sure that the non-flow contribution of the third Fourier coefficient is also insignificant. Moreover, in our previous work, it is shown that [29] V 3� is slightly dependent on φs while v3 remains constant. Such feature is consistent with the STAR data (see Fig. 1 of ref. [31]), and is understood as owing to an interplay between the multiplicity fluctuations and the background flow. The above considerations, therefore, strengthen the role played by the triangular flow. Recently, such a study was carried out by the STAR collaboration [32]. In their work, in order to eliminate the con- tributions from the collective flow, the background includes those of the elliptic, triangular, and quadrangular flow. It is somewhat surprising to find that resultant correlation yields still maintain the same feature on the away side: it shows one peak in the in-plane trigger direction and double peaks in the out-of-plane trigger direction. In other words, after the subtraction of the con- tributions from the flow coefficients, some non-trivial structure remains. As the STAR experiment claims, it may be attributed to the pathlength-dependent jet-quenching. What we understand by ZYAM is that it is a method to single out the relevant signal in the resultant correlations, by subtraction of the background evaluated in terms of the average flow harmonics, from the proper correla- tions including the information on event-by-event fluctuating flow. Whether the method is good or not, we don’t know. But the initial fluctuations should remain in the resultant correlations. What can be drawn from the similarity shown in the two STAR papers [28, 32] is the following. The contribution of the average triangular flow 〈v3〉 seems to be much smaller than 〈v2〉. Also, because of the different |�η| restrictions applied in the data analysis, the contri- bution from the jet should be either very small or insensitive to any |�η| cut. In this context, it is worthwhile to carry out an ex- plicit calculation to verify whether the data presented in [32] can be reproduced by a hydrodynamic approach. The above discussions motivated the present study. Here, we present results on event plane dependence of the di-hadron corre- lations by using NeXSPheRIO [33], an ideal hydrodynamic model, with event-by-event fluctuating initial conditions. Since all the pa- rameters of the model are fixed by reproducing the observed mul- tiplicity yields [34,35], there is no free parameter in the present calculation. We then implement the same flow background due to the elliptic, triangular and quadrangular flows as the STAR col- laboration did, and compare the calculated correlations with the data [32]. As a confirmation, we also evaluate and present the correlations by a subtraction only of the elliptic and quadrangu- lar flows, in comparison with the data [28]. The rest of the paper is organized as follows. In Section 2, we show the results of our hy- drodynamical simulations together with discussions, and the con- cluding remarks are given in Section 3. 2. Event plane dependence of di-hadron correlations The event-plane dependence analysis evaluates the correlated di-hadron pairs as a function of the azimuthal angle difference �φ at different trigger angles φs = |φtrigger − �2|, which is measured with respect to the event plane of the second harmonic coefficient. In accordance with the experimental data [32], we study the colli- sions in the 20–60% centrality window. The transverse momentum range of the trigger particles is chosen to be 3 ≤ pt T ≤ 4 GeV, and that of the associated particles is 1 ≤ pa T ≤ 2 GeV. To accommodate the experimental setup of the STAR collaboration, the calculations are carried out in the pseudorapidity interval |η| < 1. Then the results are divided into six equally sized slices in the azimuthal angle of the trigger particles. To calculate the event plane, we make use of all charged particles within the transverse momen- tum range pT < 2 GeV. By carrying out the modified reaction-plane (MRP) method [36] as employed by the STAR collaboration, parti- cle pairs with |�η| < 0.5 are excluded from the construction of the event planes. Correlated pairs with |�η| < 0.7 between the trig- ger particles and the associated particles are also excluded from the analysis because the original purpose was to minimize the near-side jet contribution [32]. The ZYAM method is then used to construct the correlation pattern due to the anisotropic flow. The primary flow correlated background is from the elliptic flow, caused by the average almond shape of the superposition region of the initial energy distribution, the triangular flow, caused by the fluctuations of the initial conditions that occur on an event- by-event basis [20], as well as the quadrangular flow. The above flow harmonics contribute to the two-particle correlations and are to be subtracted from the proper correlation pattern. In [32], the flow correlated background from the elliptic, triangular as well as quadrangular flow are expressed as follows dN d�φ = B ( 1 + 2va 2 vt 2 cos 2�φ + 2va 3vt 3 cos 3�φ + 2va 4{�2}vt 4{�2} cos 4�φ ) (1) where B is the background normalization, va 2, va 4{�2} (vt 2 and vt 4{�2}) are the second and fourth harmonic coefficients of the as- sociated (trigger) particles with respect to the event plane of the second harmonic �2, vt 3 and va 3 are triangular flow of the trigger and the associated particles, calculated with respect to the event plane of the third harmonic �3. The harmonic coefficients of the trigger particle are the average value obtained in the respective slice of the azimuthal angle φs . In the calculations using NeXSPheRIO, the hydrodynamic simu- lations are carried out in smaller centrality bins of 10% each. The reason for such a division is to increase the statistics for more peripheral collisions, as we generate a total of 1300 events in the 20–30% centrality window, 2400 events in the 30–40% cen- trality windows, 4400 events in the of 40–50% and 4600 events W.M. Castilho et al. / Physics Letters B 777 (2018) 369–373 371 Fig. 1. The azimuthal di-hadron �φ correlations for different values of φs = |φt −ψ2|. The proper correlations are obtained by applying a cut on the pseudo-rapidity difference |�η| < 0.7 between the trigger and associated particles. The resultant correlation is obtained by using ZYAM method, with v2 and v3 subtracted. The NeXSPheRIO results are shown by the solid purple curves, and the STAR data [32] are represented by the red stars whereas the gray area between the solid lines indicates the uncertainties. (For interpretation of the references to color in this figure, the reader is referred to the web version of this article.) Fig. 2. The azimuthal di-hadron �φ correlations for different values of φs = |φt − ψ2|. The correlations are obtained by using ZYAM method, with subtraction of v2 and v4, but not v3. The results are shown in solid purple curves. The STAR data [28] are represented by the red stars whereas the gray area between the solid lines indicates the uncertainties. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) in the 50–60% centrality windows. To further increase the statis- tics, for each event with given initial conditions, the Monte Carlo hadron generator is invoked 200 times. The obtained correlations of individual bins are then averaged to obtain the desired proper correlation. To evaluate the background correlation, the harmonic coefficients are obtained by the event-plane method [37–39]. Sub- sequently, the ZYAM method is made use of and is implemented according to Eq. (1). The resultant correlations are shown in Fig. 1. In the plots, the solid purple curves are those obtained by NeX- SPheRIO and the data are represented by the red stars whereas the gray areas between the solid lines indicate the uncertainties. Notice that, even though the contribution of the triangular flow is explicitly subtracted from the proper correlation as shown in (1), the resultant correlation is still featured by one peak in the away- side for the in-plane direction, with its maximum at �φ ≈ π , which evolves to double peaks for the out-of-plane direction. Re- mark that similar feature had been reported by STAR collaboration in [27,28] where only the elliptic and quadrangular flows were subtracted. As a comparison, we proceed to evaluate the two-particle cor- relation with the same specifications of [28]. We note that the experimental configurations, such as the pseudo-rapidity cut of |�η| > 0.7 for the counting of particle pairs as well as the mo- mentum ranges in consideration, are exactly the same in this case, and the differences come from the pseudo-rapidity filters in the construction of event planes. Therefore the same data sets generated previously by the hydrodynamic simulations are used, while the event plane �2 is evaluated by taking into considera- tion all the particles while implementing a cut in pseudo-rapidity |�η| < 0.35. Both v2 and v4 are calculated with respected to the event plane �2. Now, the correlated flow background is estimated as dN d�φ = B ( 1 + 2va 2vt 2 cos 2�φ + 2va 4{�2}vt 4{�2} cos 4�φ ) (2) The resulting two-particle correlations are shown in Fig. 2, com- pared with the STAR data from Ref. [28]. It is found that the dou- ble peak observed on the away side in the out-of-plane direction is also reasonably reproduced by the hydrodynamic calculations, which is consistent with our previous results [29] obtained by us- ing the cumulant method. It is somewhat surprising to find that the subtracted two-particle correlations in Fig. 1 and 2 are not very different. This is because the triangular flow is understood to play a significant role in the explanation of the double-peak on the away side. Consequently, it was somehow expected that the subtraction of v3 might lead to the disappearance of the same on the away side. In order to quantitatively understand the contribu- tion of v3 in our hydrodynamical results, we present in Tables 1 and 2 the calculated values of flow harmonics used in the ZYAM subtraction to produce Fig. 1. We note that the data analysis carried out in this work is not for the most central window, and therefore, the average values of v3 for the associated particles are smaller in comparison with those of v2, as shown in Table 2. This is also the case for the trig- ger particles at the in-plane directions, as can be readily verified in the first few columns of Table 1. One also observes that the mag- nitude of v4{�} for the associated particles follows the same trend of that of v2. However, since its size is much smaller than that of v2, it does not play a significant role in the resultant correla- tions. On the other hand, for the trigger particles, as one goes to the out-of-plane directions, the elliptic flow vt 2 starts to decrease and eventually becomes negative, as confirmed by the three right- most columns in Table 1. This is expected, since vt 2 is evaluated with respect to the event plane �2 of the associated particles, so this is implied by definition. Owing to the change of the sign of v2, the contribution of v3 in the background Eq. (1) is relatively small in the in-plane direction in comparison to that in the out-of-plane direction. On the other hand, for the calculated correlations, the appearance of the double peak in the out-of-plane direction indi- cated that the observed modulation of the resultant correlation is dominated by the second order harmonic coefficient V 2 ≡ va 2 vt 2 in the in-plane direction and the remnant of V 3 ≡ va 3 vt 3 in the out- of-plane direction. By Eq. (1), this is consistent with the previous observations on the magnitudes of v2 and v3. In Ref. [32], it is discussed that the trends of the away- side correlation might underscore the importance of path-length- dependent jet-medium interactions. As the present hydrodynam- 372 W.M. Castilho et al. / Physics Letters B 777 (2018) 369–373 Table 1 The calculated flow harmonics vn of trigger particles for transverse momentum range of 3 < pt T < 4 GeV/c for different azimuthal angles and centrality windows. φs ≡ |φt − ψ2| 0 − π/12 π/12 − π/6 π/6 − π/4 π/4 − π/3 π/3 − 5π/12 5π/12 − π/2 20–30% v2 0.1899 0.1222 0.0345 −0.0235 −0.0525 −0.0649 v3 0.0304 0.0065 0.0065 0.0314 0.0288 0.0059 v4{�2} 0.1649 0.0036 −0.1043 −0.0804 −0.0011 −0.0559 30–40% v2 0.2030 0.1239 0.0344 −0.0226 −0.0488 −0.0594 v3 0.0311 0.0073 0.0066 0.0298 0.0297 0.0068 v4{�2} 0.1767 0.0038 −0.1039 −0.0779 −0.0016 −0.0513 40–50% v2 0.2038 0.1277 0.0345 −0.0219 −0.0497 −0.0616 v3 0.0323 0.0075 0.0063 0.0325 0.0306 0.0063 v4{�2} 0.1773 0.0054 −0.1028 −0.0774 −0.0008 −0.0532 50–60% v2 0.2035 0.1260 0.0368 −0.0236 −0.0505 −0.0664 v3 0.0313 0.0052 0.0069 0.0295 0.0282 0.0079 v4{�2} 0.1767 0.0028 −0.1075 −0.0806 −0.0014 −0.0574 Table 2 The calculated flow harmonics vn of associated particles for transverse momentum range of 1 < pa T < 2 GeV/c for different azimuthal angles and centrality windows. 20–30% 30–40% 40–50% 50–60% v2 0.1028 0.1146 0.1199 0.1189 v3 0.0399 0.0421 0.0445 0.0449 v4{�2} 0.0106 0.0138 0.0144 0.0135 ical simulations are able to reproduce observed feature of two- particle correlation, it strongly indicates that the observed cor- relations are likely to be a collective-flow effect of the system. Moreover, it seems that the ZYAM procedure, devised to essen- tially subtract the contribution of collective flow from the proper two-particle correlation, has somehow failed in its purpose. In par- ticular, it is found that the subtraction of v3 does not affect the essential feature of the resultant correlations, namely, the relative magnitudes between V 2 and V 3. To us, this might be related to the event-by-event fluctuating initial conditions and their impact on flow harmonics. This is because the subtracted vn in Eq. (1) is, in fact, the event average value, 〈vn〉. Due to the event-by-event fluctuations, the event average value of a product of harmonic co- efficients can be significantly different from the product of the corresponding average values. In other words, not only the mag- nitude of the triangular flow, v3, is understood to be related to the event-by-event fluctuations, its fluctuations might also play a non- trivial role in the particle correlations. Moreover, the event planes between different harmonics might be correlated for a given event but uncorrelated among the various events, which further compli- cates the problem. As a result, the average of the Fourier expansion in azimuthal angle cannot be simply approximated by a Fourier ex- pansion in terms of the products of average harmonic coefficients, even rescaled by the ZYAM scheme. The present calculations em- ploying NeXSPheRIO give reasonable results for two different sets of data obtained by different procedures. It implies that the fluc- tuating initial conditions generated by NeXuS are mostly realistic. We also note that this topic is closely related to the correlation between different flow harmonics, as recently explored by several authors [40–42]. It is noted that in the above calculations, the event-plane method is employed for the background subtraction, as a part of our strategy to closely follow the same steps in the data analysis. However, the event-plane method carried out in a hydrodynamic study might not be equivalent to that in experiment analysis. Ac- cording to the calculations by the Monte Carlo Glauber model [43], the obtained v2 depends on the resolution R of the event-plane method. As long as the resolution is high enough, the calculated result gives the desired mean value of vn . However, as the res- olution decreases, the obtained result gradually approaches the RMS value of corresponding harmonic coefficients. Since in a hy- drodynamic calculation, the resolution is typically better than that of experimental measurements, the current calculations likely un- derestimate the flow background. In practice, the resolution of a hydrodynamic approach can be controlled by the number of Monte Carlo simulations during the hadronization phase. One may there- fore roughly estimate the discrepancy owing to the resolution, by extrapolating the present results to the corresponding resolution of the data by varying the number of Monte Carlo simulations. We have carried out the analysis, and it shows that the flow harmon- ics used in the flow background subtraction is around 10% smaller than those employed in STAR analysis. Since the magnitudes of flow harmonics shown in Tables 1 and 2 are small, the deviations of their contributions in Eqs. (1) and (2) are even less. As a result, the values of B in Eqs. (1) and (2) should be slightly smaller but mostly unchanged. We have verified that the resultant correlation after the flow background subtraction will be slightly smaller, but the qualitative results shown in Fig. 1 and 2 remain unchanged. 3. Concluding remarks As discussed above, we understand that ZYAM is a method aimed to single out relevant signal in the resultant correlations, by subtraction of the flow background evaluated in terms of the average harmonic coefficients, from the proper correlations includ- ing the event-by-event fluctuating flow. To us, it seems possible that the initial fluctuations still present, or even become domi- nant in the resultant correlations. The similarity presented in the two STAR papers [28,32] indicates that the flow is fluctuating and 〈v3〉 � 〈v2〉. Also, owing to the different |�η| restrictions applied in the data analysis, the jet contribution is likely to be very in- significant. In this work, we show by explicit calculations that hydrodynamics, with fluctuating initial conditions, is able to re- produce the observed double-peak structure of two-particle corre- lation on the away-side for the out-of-plane triggers, even when the triangular flow is subtracted by using ZYAM method. There- fore, the present study further strengthens the idea that observed correlations are mostly of hydrodynamic origin. W.M. Castilho et al. / Physics Letters B 777 (2018) 369–373 373 From a hydrodynamic viewpoint, it is understood that the role of v3 is closely associated with the event-by-event geometrical fluctuations in the initial conditions. The existence of two peaks, in the away side correlation, clearly shows the need of v3. Our results indicate that the observed event plane dependence of two-particle correlations, for the most part, can be reproduced by a hydrody- namic approach. The physical mechanism behind the findings may be attributed to the correlation and/or fluctuation of flow harmon- ics, which is closely associated with the event-by-event fluctua- tions in the initial conditions. This is because of the assumption of the background correlation module, Eq. (1), assumes the picture of average flow with no correlation neither fluctuation of signifi- cant consequence. Ongoing efforts on the event-plane correlations, symmetric cumulant [40–42] may likely provide a better insight of the problem, and a clear criterion to distinguish between different mechanisms. Some alternative methods may also be meaningful, as shown in our recent work concerning centrality dependence of di- hadron correlations [30]. There, the cumulant method gives quite a similar result as the ZYAM one. We note that the 2+1 correla- tion [44,45] employed by STAR Collaboration may also serve as a good tool to disentangle the jet signal. On the other hand, the peripheral tube model [15], in which the interplay between random hot tubes in the initial conditions and collective background flow plays an important role, gives a unified description of the “ridge” structure, both for the near-side and the away-side ones. As a result, in this interpretation, differ- ent flow harmonics of an individual event are naturally “born” together, and consequently correlated, their appearance does not depend on any global structure of the initial condition. From the viewpoint of the peripheral tube model, it is possible that the correlation between v2 and v3 owing to the flow deflected by high energy tube may provide an intuitive explanation of the ob- served data. In particular, the lack of the near-side peak for the out-of-plane triggers, which would be produced by this two-peak structure in the single-particle distribution, might be studied ana- lytically in this framework. Recently, a Reaction Plane Fit method was proposed [46] and employed to estimate the correlation func- tions in the background dominated region on the near-side [47]. The resulting correlation does not show the double peak on the away side, neither any dramatic shape modification as a function of centrality. Therefore, the authors conclude that the Mach cone is an artifact of the background subtraction and the jets do not fully equilibrate with the medium. These results further indicate that the effect of the jet in the di-hadron correlation is indeed a subtle subject. We plan to explore these topics further and look for a possible criterion to distinguish between different approaches in the near future. Acknowledgments We are thankful for valuable discussions with Jiangyong Jia. We gratefully acknowledge the financial support from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES). A part of the work was developed under the project INCTFNA Proc. No. 464898/2014-5. This research is also supported by the Center for Scientific Computing (NCC/GridUNESP) of the São Paulo State University (UNESP). References [1] STAR Collaboration, B. Abelev, et al., Phys. Rev. C 80 (2009) 064912, arXiv: 0909.0191. [2] STAR Collaboration, B. Abelev, et al., Phys. Rev. Lett. 102 (2009) 052302, arXiv: 0805.0622. [3] PHENIX Collaboration, A. Adare, et al., Phys. Rev. C 77 (2008) 011901, arXiv: 0705.3238. [4] PHENIX Collaboration, A. Adare, et al., Phys. Rev. C 78 (2008) 014901, arXiv: 0801.4545. [5] PHOBOS Collaboration, B. Alver, et al., Phys. Rev. Lett. 104 (2010) 062301, arXiv: 0903.2811. [6] CMS Collaboration, S. Chatrchyan, et al., Eur. Phys. J. C 72 (2012) 2012, arXiv: 1201.3158. [7] CMS Collaboration, S. Chatrchyan, et al., J. High Energy Phys. 1107 (2011) 076, arXiv:1105.2438. [8] E.V. Shuryak, Phys. Rev. C 76 (2007) 047901, arXiv:0706.3531. [9] C.B. Chiu, R.C. Hwa, Phys. Rev. C 79 (2009) 034901, arXiv:0809.3018. [10] CMS, V. Khachatryan, et al., J. High Energy Phys. 09 (2010) 091, arXiv: 1009.4122. [11] ATLAS, G. Aad, et al., Phys. Rev. Lett. 110 (2013) 182302, arXiv:1212.5198. [12] ALICE, B. Abelev, et al., Phys. Lett. B 719 (2013) 29, arXiv:1212.2001. [13] CMS, S. Chatrchyan, et al., Phys. Lett. B 718 (2013) 795, arXiv:1210.5482. [14] ALICE, B.B. Abelev, et al., Phys. Lett. B 741 (2015) 38, arXiv:1406.5463. [15] R. Andrade, F. Grassi, Y. Hama, W.-L. Qian, J. Phys. G 37 (2010) 094043, arXiv: 0912.0703. [16] Y. Hama, R.P.G. Andrade, F. Grassi, W.-L. Qian, Nonlinear Phenom. Complex Syst. 12 (2009) 466, arXiv:0911.0811. [17] R.P.G. Andrade, F. Grassi, Y. Hama, W.-L. Qian, Phys. Lett. B 712 (2012) 226, arXiv:1008.4612. [18] Y. Hama, R.P. Andrade, F. Grassi, J. Noronha, W.-L. Qian, Acta Phys. Pol. Suppl. 6 (2013) 513, arXiv:1212.6554. [19] W.-L. Qian, et al., J. Phys. G 41 (2014) 015103, arXiv:1305.4673. [20] B. Alver, G. Roland, Phys. Rev. C 81 (2010) 054905, arXiv:1003.0194. [21] D. Teaney, L. Yan, Phys. Rev. C 83 (2011) 064904, arXiv:1010.1876. [22] H. Petersen, G.-Y. Qin, S.A. Bass, B. Muller, Phys. Rev. C 82 (2010) 041901, arXiv: 1008.0625. [23] B. Schenke, S. Jeon, C. Gale, Phys. Rev. Lett. 106 (2011) 042301, arXiv: 1009.3244. [24] G.-Y. Qin, H. Petersen, S.A. Bass, B. Muller, Phys. Rev. C 82 (2010) 064903, arXiv: 1009.1847. [25] J. Xu, C.M. Ko, Phys. Rev. C 83 (2011) 021903, arXiv:1011.3750. [26] G.-L. Ma, X.-N. Wang, Phys. Rev. Lett. 106 (2011) 162301, arXiv:1011.5249. [27] STAR Collaboration, A. Feng, J. Phys. G 35 (2008) 104082, arXiv:0807.4606. [28] STAR Collaboration, H. Agakishiev, et al., arXiv:1010.0690, 2010. [29] W.-L. Qian, R. Andrade, F. Gardim, F. Grassi, Y. Hama, Phys. Rev. C 87 (2013) 014904, arXiv:1207.6415. [30] W.M. Castilho, W.-L. Qian, F.G. Gardim, Y. Hama, T. Kodama, Phys. Rev. C 95 (2017) 064908, arXiv:1610.04108. [31] M. Luzum, Phys. Lett. B 696 (2011) 499, arXiv:1011.5773. [32] STAR, H. Agakishiev, et al., Phys. Rev. C 89 (2014) 041901, arXiv:1404.1070. [33] Y. Hama, T. Kodama, O. Socolowski Jr., Braz. J. Phys. 35 (2005) 24, arXiv:hep- ph/0407264. [34] W.-L. Qian, et al., Braz. J. Phys. 37 (2007) 767, arXiv:nucl-th/0612061. [35] D.M. Dudek, et al., arXiv:1409.0278, 2014. [36] STAR, J. Adams, et al., Phys. Rev. C 72 (2005) 014904, arXiv:nucl-ex/0409033. [37] S. Voloshin, Y. Zhang, Z. Phys. C 70 (1996) 665, arXiv:hep-ph/9407282. [38] A.M. Poskanzer, S.A. Voloshin, Phys. Rev. C 58 (1998) 1671, arXiv:nucl-ex/ 9805001. [39] S.A. Voloshin, A.M. Poskanzer, R. Snellings, arXiv:0809.2949, 2008. [40] ATLAS, G. Aad, et al., Phys. Rev. C 90 (2014) 024905, arXiv:1403.0489. [41] J. Jia, J. Phys. G 41 (2014) 124003, arXiv:1407.6057. [42] ALICE, J. Adam, et al., Phys. Rev. Lett. 117 (2016) 182301, arXiv:1604.07663. [43] B. Alver, et al., Phys. Rev. C 77 (2008) 014906, arXiv:0711.3724. [44] STAR, H. Pei, J. Phys. Conf. Ser. 316 (2011) 012016. [45] STAR, K. Kauder, Nucl. Phys. A 830 (2009) 685C, arXiv:0907.4673. [46] N. Sharma, J. Mazer, M. Stuart, C. Nattrass, Phys. Rev. C 93 (2016) 044915, arXiv:1509.04732. [47] C. Nattrass, N. Sharma, J. Mazer, M. Stuart, A. Bejnood, Phys. Rev. C 94 (2016) 011901, arXiv:1606.00677. http://refhub.elsevier.com/S0370-2693(17)31043-2/bib737461722D726964676532s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib737461722D726964676532s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib737461722D726964676533s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib737461722D726964676533s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7068656E69782D726964676534s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7068656E69782D726964676534s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7068656E69782D726964676535s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7068656E69782D726964676535s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib70686F626F732D726964676536s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib70686F626F732D726964676536s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib636D732D726964676533s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib636D732D726964676533s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib636D732D726964676534s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib636D732D726964676534s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7163642D636F72722D3031s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7163642D636F72722D3032s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib636D732D726964676537s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib636D732D726964676537s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib61746C61732D72696467652D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib616C6963652D72696467652D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib636D732D726964676539s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib616C6963652D72696467652D32s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D32s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D32s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D33s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D33s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D34s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D34s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D37s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D37s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D766E2D34s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib687964726F2D76332D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib687964726F2D76332D32s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib687964726F2D76332D38s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib687964726F2D76332D38s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib687964726F2D76332D35s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib687964726F2D76332D35s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib687964726F2D76332D34s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib687964726F2D76332D34s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib687964726F2D76332D36s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib687964726F2D76332D37s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib737461722D706C616E6531s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib737461722D706C616E6532s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D65762D34s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D65762D34s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D65762D36s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D636F72722D65762D36s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib70682D636F72722D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib737461722D706C616E6533s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D7265766965772D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D7265766965772D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D656F732D32s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7370682D656F732D33s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib737461722D726964676536s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib6576656E742D706C616E652D6D6574686F642D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib6576656E742D706C616E652D6D6574686F642D32s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib6576656E742D706C616E652D6D6574686F642D32s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib6576656E742D706C616E652D6D6574686F642D33s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib61746C61732D766E33s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib61746C61732D766E35s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib616C6963652D766E35s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib70686F626F732D6D63676C61756265722D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib737461722D636F72722D32706C7573312D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib737461722D636F72722D32706C7573312D32s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7A79616D2D7270662D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7A79616D2D7270662D31s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7A79616D2D7270662D32s1 http://refhub.elsevier.com/S0370-2693(17)31043-2/bib7A79616D2D7270662D32s1 Event-plane dependent di-hadron correlations with harmonic vn subtraction in a hydrodynamic model 1 Introduction 2 Event plane dependence of di-hadron correlations 3 Concluding remarks Acknowledgments References