Elliptic flow of electrons from beauty-hadron decays extracted from Pb--Pb collision data at $\sqrt{s_{\rm NN}}$ = 2.76 TeV

We present a calculation of the elliptic flow of electrons from beauty-hadron decays in semi-central Pb--Pb collisions at centre-of-mass energy per colliding nucleon pair, represented as $\sqrt{s_{\rm NN}}$, of 2.76 TeV. The result is obtained by the subtraction of the charm quark contribution in the elliptic flow of electrons from heavy-flavour hadron decays in semi-central Pb--Pb collisions at \snn = 2.76 TeV recently made publicly available by the ALICE collaboration.


Introduction
The nuclear matter exposed to conditions of high temperature and energy density is expected to undergo a phase transition to a colour deconfined state of matter [1,2], the Quark-Gluon Plasma (QGP). The conditions for the phase transition can be achieved in the laboratory with collisions of heavy ions at high energies [1]. The properties of the medium formed in the laboratory can be probed with a unique degree of control by particles from decays of heavy flavours (charm and beauty), since heavy quarks are mainly produced in hard parton scattering processes [3][4][5] at the initial stage of heavy-ion collisions [6,7] and participate in the entire evolution of the created system. The partons traversing the medium lose energy via collisional and radiative processes [8][9][10][11][12] in the interaction with the medium constituents. The energy loss of partons is predicted to be dependent on their colour charge and mass, resulting in a hierarchy where beauty quarks lose less energy than charm quarks, charm quarks lose less energy than light quarks, and quarks lose less energy than gluons [10,13,14]. The heavy-flavour energy loss can be investigated experimentally with the nuclear modification factor (R AA ) of heavy-flavour particles, defined as where dN AA /dp T is the transverse momentum (p T ) differential yield in nucleus-nucleus (AA) collisions; T AA is the average nuclear overlap function in nucleus-nucleus collisions, given by the ratio of the average number of binary collisions and the inelastic cross section; and dσ pp /dp T is the p T -differential cross section in proton-proton (pp) collisions. A suppression of the yield of D mesons and leptons from heavy-flavour hadron decays (R AA < 1) for p T > 3 GeV/c was observed in gold-gold (Au-Au) collisions at √ s NN = 200 GeV at the Relativistic Heavy Ion Collider (RHIC) [15][16][17][18] and in lead-lead (Pb-Pb) collisions at √ s NN = 2.76 and 5.02 TeV at the Large Hadron Collider (LHC) [19][20][21][22][23][24][25], indicating energy loss of heavy flavours in the medium. An experimental hint to the quark mass dependence of the heavyflavour energy loss has been found in the comparison of the R AA of D mesons and non-prompt J/ψ from B-hadron decays in central Pb-Pb collisions at √ s NN = 2.76 TeV at the LHC [20,26,27]. The observed difference of these measurements is described by model calculations [28,29] as predominantly due to the quark mass dependence of the parton energy loss.
The interaction of heavy quarks with the medium can be further investigated with the azimuthal anisotropy of heavy-flavour particles, extracted from the coefficients v n of the Fourier decomposition of the particle azimuthal distribution in the transverse plane [30] where ϕ is the azimuthal angle of the heavy-flavour particles and Ψ n is the symmetry-plane angle of the n th -order harmonic. The second Fourier coefficient v 2 , called elliptic flow, quantifies the elliptic azimuthal anisotropy of the emitted particles. The origin of the elliptic azimuthal anisotropy of heavy-flavour particles in non-central heavy-ion collisions depends on the transverse-momentum interval. While the v 2 at low p T is sensitive to the collective motion of the medium constituents caused by pressure gradients, the v 2 at high p T can constrain the path-length dependence of the in-medium energy loss of heavy quarks, resulting from the direction of the particles that traverse the ellipsoidal nuclear overlap region. The elliptic flow of prompt D mesons at mid-rapidity is observed to be positive in 30-50% Pb-Pb collisions at √ s NN = 2.76 TeV at the LHC [31,32] with 5.7σ significance in the interval 2 < p T < 6 GeV/c, indicating that charm quarks participate in the collective motion of the system. Measurements of the prompt D-meson v 2 in Pb-Pb collisions at √ s NN = 5.02 TeV [33,34] have smaller uncertainties compared to the ones in Pb-Pb collisions at √ s NN = 2.76 TeV. The results at the two collision energies are compatible within uncertainties. The prompt D + s v 2 in semi-central Pb-Pb collisions at √ s NN = 5.02 TeV is compatible within uncertainties with the average of prompt D 0 , D + , and D * + v 2 in the same collision system [34]. A positive v 2 is also observed for leptons from heavy-flavour hadron decays at low and intermediate p T in semi-central Au-Au collisions at √ s NN = 200 GeV at RHIC [35,36] and in semi-central Pb-Pb collisions at √ s NN = 2.76 TeV at the LHC [25,37,38]. In particular, the v 2 of electrons from heavy-flavour hadron decays is observed to be positive with 5.9σ significance in the range 2 < p T < 2.5 GeV/c in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV. In view of the experimental results on the elliptic flow of heavy-flavour particles, an important question that remains open is whether beauty quarks take part in the collective motion in the medium. The first measurement of the v 2 of non-prompt J/ψ mesons from B-hadron decays is compatible with zero within uncertainties in two kinematic regions, 6.5 < p T < 30 GeV/c and |y| < 2.4, and 3 < p T < 6.5 GeV/c and 1.6 < |y| < 2.4, in 10-60% Pb-Pb collisions at √ s NN = 2.76 TeV at the LHC [27]. In this paper, we present a method to subtract the contribution of charm quarks in the published measurement of the elliptic flow of electrons from heavy-flavour hadron decays in semi-central Pb-Pb collisions at √ s NN = 2.76 TeV performed by the ALICE collaboration. The calculation uses as input the v 2 coefficients of prompt D mesons and electrons from heavy-flavour hadron decays measured by the ALICE collaboration [32,37] and three different results for the relative contribution of electrons from beauty-hadron decays to the yield of electrons from heavy-flavour hadron decays [39][40][41][42]].

Methodology
The particle azimuthal distribution of electrons from heavy-flavour hadron decays (e ← c + b) can be separated into the contributions of electrons from charm-hadron decays (e ← c) and from beauty-hadron decays (e ← b). Consequently, the elliptic flow of electrons from beauty-hadron decays can be expressed as v e←b where R represents the relative contribution of electrons from beauty-hadron decays to the yield of electrons from heavy-flavour hadron decays.
In the following, we present the currently published measurements and, in case there is no available measurement, our calculations of the three observables required to obtain the elliptic flow of electrons from beauty-hadron decays. Based on available results on open heavy flavours at RHIC and LHC, the most suitable system for this analysis is the Pb-Pb collision system at √ s NN = 2.76 TeV in the 20-40% centrality class, which corresponds to the centrality range where the measured v 2 of electrons from heavyflavour hadron decays is observed to be positive with a maximum significance [37] and thus a possible elliptic flow of electrons from beauty-hadron decays is expected to be more significant. In this analysis, the v 2 and R AA of heavy-flavour particles are assumed to be the same at slightly different mid-rapidity ranges (|y| < 0.5, 0.7 and 0.8) in which the measurements needed in the calculation are available. Indeed, no dependence on rapidity was observed in recent ALICE results on those observables for electrons from heavy-flavour hadron decays at mid-rapidity (|y| < 0.7 for v 2 and |y| < 0.6 for R AA measurements) and muons from heavy-flavour hadron decays at forward rapidity (2.5< y <4) [24,37].

Elliptic flow of electrons from heavy-flavour hadron decays
The result on the elliptic flow of electrons from heavy-flavour hadron decays (v e←c+b 2 ) at mid-rapidity (|y| < 0.7) in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV published by the ALICE collaboration [37] is used in this analysis. The v e←c+b 2 is measured in the interval 0.5 < p T < 13 GeV/c with the event plane method [30]. A positive value is observed in the interval 2 < p T < 2.5 GeV/c with significance of 5.9σ [37].

Relative contribution of electrons from beauty-hadron decays to the yield of electrons from heavy-flavour hadron decays
The measurement of the relative contribution of electrons from beauty-hadron decays to the yield of electrons from heavy-flavour hadron decays (R) has been published by the ALICE collaboration only in pp collisions at √ s = 2.76 TeV [39,40]. The coefficient R measured in pp collisions is used in the analysis with the caveat that initial-and final-state effects modify the yield of electrons from heavyflavour hadron decays in heavy-ion collisions. In particular, the coefficient R at high p T is expected to be higher in Pb-Pb collisions compared to pp collisions, since the in-medium energy loss of charm quarks is predicted to be larger than the one of beauty quarks [23]. Therefore, the factor R at high p T in Pb-Pb collisions at √ s NN = 2.76 TeV is expected to have an exclusive value between the measured factor R in pp collisions at √ s = 2.76 TeV and unity. Consequently, according to Eq. 3, the minimum value of the v 2 of electrons from beauty-hadron decays can be computed with the R measured in pp collisions.
In addition to the available measurement, the coefficient R is obtained with a Monte Carlo (MC) simulation based on POWHEG [43], which provides the calculation of the heavy-flavour production in hadronic collisions at Next-to-Leading Order (NLO) accuracy. The POWHEG results are interfaced to PYTHIA [44,45] in order to generate the shower, hadronisation and decay. In agreement with other heavy-flavour production tools, e.g. pQCD calculation at Fixed Order plus Next-to-Leading Logarithms (FONLL) [4,46] and earlier pp calculations [5], the square root of the quadratic sum of the quark mass (m Q ) and p T are used as renormalization and factorization scales, i.e. µ = µ f = µ r = m 2 Q + p 2 T . The charm-and beauty-quark masses are set as 1.5 GeV and 4.75 GeV, respectively. Even though the calculated coefficient R is sensitive to the choice of heavy-quark masses and scales, only the central value is shown in this analysis. Admittedly, the described framework is designed for pp collisions, but by making use of the EPS09 [47] NLO nuclear Parton Distribution Functions (nPDFs) the framework is able to account for initial-state cold nuclear effects. The nPDF gluon shadowing results in reduced p T -differential cross sections for electrons from heavy-flavour hadron decays for p T < 6 GeV/c and affects contributions from charm quarks stronger than those from beauty quarks. Thus, it provides a lower baseline for  [43][44][45] at NLO accuracy using EPS09 NLO nPDFs, with PYTHIA at LO accuracy using EPS09 LO nPDFs, and with FONLL calculation [4,46] using CTEQ6.6 PDFs. Only the central values are shown. the factor R, which is suggested to be further enhanced by medium interactions as it will be discussed in this paper. In addition, the R coefficient is also obtained with a leading order (LO) calculation based on PYTHIA using EPS09 LO nPDFs to study the impact of NLO corrections [5]. The comparison of the calculations with LO and NLO approaches, shown in Fig. 1, reveals that the factor R is reduced with the NLO corrections, stressing the importance of NLO calculations. In fact, the additional processes of heavy-flavour production at NLO give rise to large logarithmic corrections to the charm-and beautyquark cross sections depending on the heavy-flavour mass. The corresponding FONLL calculation of the factor R using CTEQ6.6 PDFs, which is also shown in Fig. 1, is similar to the POWHEG+PYTHIA (EPS09NLO) result at high p T .
The result on the factor R from the BAMPS heavy-flavour transport model [41,42], which includes collisional and radiative in-medium energy loss of heavy quarks, is also employed in the analysis to obtain the v 2 of electrons from beauty-hadron decays. The choice of the BAMPS model is justified by the good agreement of the predictions for the R AA of electrons from beauty-and heavy-flavour hadron decays for p T > 3 GeV/c in central Pb-Pb collisions at √ s NN = 2.76 TeV with what measured by the ALICE collaboration [23,24].
In this analysis, the R coefficient measured in pp collisions using the electron-hadron azimuthal technique [39,40] and the ones obtained with POWHEG+PYTHIA (EPS09NLO) and with the BAMPS model are used to estimate the elliptic flow of electrons from beauty-hadron decays. The R coefficient measured in pp collisions is also obtained with the track impact parameter method [39,40]. Results obtained with both techniques are compatible within uncertainties.

Elliptic flow of electrons from charm-hadron decays
The elliptic flow of electrons from charm-hadron decays (v e←c 2 ) is estimated using a MC simulation of decays of D 0 mesons into electrons with PYTHIA. The MC simulation is based on two observables measured for D 0 mesons in Pb-Pb collisions at √ s NN = 2.76 TeV: -the p T -differential yield, which is used as a probability distribution for finding a D 0 meson with a certain p T ; -the p T -differential v 2 , which is used to obtain the ϕ D 0 − Ψ 2 probability distribution with Eq. 2.
In fact, the p T -differential yield of D 0 mesons at mid-rapidity (|y| < 0.8) in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV is estimated from the ALICE results on the p T -differential yield and R AA of prompt D 0 mesons at mid-rapidity (|y| < 0.5) in 0-20% Pb-Pb collisions at the same collision energy [19] as dN AA dp T 20−40% where the coefficient C ∆y = 1.6 corresponds to the scaling factor of the yield from |y| < 0.5 to |y| < 0.8 in pp collisions, assuming a uniform distribution of the D 0 -meson yield within the rapidity range. The terms C T AA = 0.362 ± 0.020 [19] and C R AA are the ratios of the average nuclear overlap function and the D 0 -meson R AA , respectively, in Pb-Pb collisions at √ s NN = 2.76 TeV in the 20-40% centrality class to the ones in the 0-20% centrality class. Note that the terms C ∆y and C T AA are constant, so they do not play a role in the determination of the D 0 -meson p T probability distribution. The non-measured R AA of D 0 mesons in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV is estimated by the average of the ALICE results on the D 0 -meson R AA in Pb-Pb collisions at the same collision energy in the 0-20% and 40-80% centrality classes [19] weighted by the corresponding yield of D 0 mesons in each centrality class. The resulting p T distribution of D 0 mesons in 20-40% Pb-Pb collisions obtained from Eq. 4 is then fitted by a power-law function (left panel of Fig. 2), considering the statistical uncertainty of the experimental results. The fit function is used as the D 0 -meson p T probability distribution.
The v 2 of prompt D 0 mesons in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV (right panel of Fig. 2) is obtained by the arithmetic average of the measured prompt D 0 -meson v 2 in Pb-Pb collisions at the same collision energy in the 10-30% and 30-50% centrality classes [32]. Indeed, experimental results show that the v 2 of heavy-flavour particles increases with the centrality class [32,35,37,38], which is consistent with the qualitative expectation of increasing of the elliptic anisotropy from central to peripheral nucleusnucleus collisions. The statistical and systematic uncertainties are propagated considering the prompt D 0meson v 2 in the 10-30% and 30-50% centrality classes as uncorrelated as a conservative estimation. In the D 0 -meson v 2 measurement by the ALICE collaboration, the central value was obtained by assuming that the v 2 coefficients of prompt D mesons and D mesons from B-meson decays are the same [32]. However, the systematic uncertainty related to this assumption, referred to as systematic uncertainty from the B feed-down subtraction, was evaluated by the ALICE collaboration. It was assumed that the v 2 of prompt D mesons from B-meson decays should be between zero and v 2 of prompt D mesons, resulting in the upper and lower limits of the systematic uncertainty, respectively. Therefore, the B feeddown contribution decreases the absolute value of the D 0 -meson v 2 and thus the systematic uncertainty is restricted to the upper (lower) limit when the v 2 is positive (negative). Since the measured D 0 -meson v 2 coefficients are negative in the 8 < p T < 12 GeV/c and 12 < p T < 16 GeV/c intervals in the 10-30% and 30-50% centrality classes, respectively, the resulting propagated systematic uncertainty from the B feed-down subtraction contains lower and upper limits.
Finally, the estimated p T and v 2 distributions of D 0 mesons in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV are used to obtain the v e←c 2 = cos [2 (ϕ e − Ψ 2 )] in the same collision system using the PYTHIA event generator. The azimuthal angle of electrons (ϕ e ) takes into account the angular separation between electrons and their parent D 0 mesons.

Statistical uncertainty
The statistical uncertainty of the D 0 -meson v 2 is used as input for the MC simulation to obtain the statistical uncertainty of the v e←c 2 . The statistical uncertainties of the measurements used to obtain the p T -differential yield of D 0 mesons are considered in the fit of the D 0 -meson probability distribution. Further variations are considered as systematic uncertainties.

Systematic uncertainty
The systematic uncertainties from data and from the B feed-down subtraction of the D 0 -meson v 2 (right panel of Fig. 2) are used as input for the MC simulation to obtain the systematic uncertainty of the v e←c 2 . The following is a discussion on other sources of systematic uncertainty that can influence the v e←c 2 estimation.
In order to validate the Eq. 4, the p T -differential yield and R AA of prompt D 0 mesons in 40-80% Pb-Pb collisions at √ s NN = 2.76 TeV [19] are also used as reference to obtain the p T -differential yield of D 0 mesons in the 20-40% centrality class. The result is the same as the one obtained with the 0-20% centrality class (left panel of Fig. 2).
The ALICE result on the D 0 -meson R AA in the 30-50% centrality class [32] is used as an alternative for the D 0 -meson R AA estimation in the 20-40% centrality class. No significant difference is observed in the resulting v e←c 2 with respect to the one obtained with the R AA estimated by the average of the D 0meson R AA measurements in the 0-20% and 40-80% centrality classes weighted by the corresponding yield of D 0 mesons in each centrality class.
The systematic uncertainties of the measurements of the p T -differential yield and R AA of prompt D 0 mesons are considered in the fit of the D 0 -meson p T distribution in 20-40% Pb-Pb collisions. No significant difference is observed in the resulting v e←c 2 with respect to the one considering only the statistical uncertainty in the fit. For further investigation, The BAMPS result on the p T distribution of D mesons at |y| < 0.8 in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV [41,42] is also used to compute the v e←c 2 . The relative difference of the obtained v e←c 2 using the estimated p T distribution of D 0 mesons and the BAMPS result, which increases from 1% to 20% in the interval 2 < p T < 8 GeV/c, is included in the systematic uncertainty.
The effect of the D 0 -meson v 2 estimation in 20-40% Pb-Pb collisions using the arithmetic average of the D 0 -meson v 2 measurements in the 10-30% and 30-50% centrality classes is investigated in this analysis. For this purpose, the trend of the unidentified charged particle v 2 as a function of the average number of binary collisions ( N coll ) [48] is assumed to be the same as the one for D 0 mesons. The v 2 as a function of N coll is obtained from a parametrisation of the centrality-dependent v 2 measurement of unidentified charged particles integrated over the interval 0.2 < p T < 5 GeV/c [49]. The corresponding result exhibits a linear dependence between v 2 and N coll with a negative slope for N coll > 220.
For comparison, the parametrisation is also obtained from the centrality-dependent v 2 measurement of unidentified charged particles integrated over the interval 10 < p T < 20 GeV/c [50]. The linear dependence between v 2 and N coll is the same as the one obtained for particles in a lower p T interval. The D 0 -meson v 2 is then obtained by the average of the D 0 -meson v 2 in the 10-30% and 30-50% centrality classes weighted by the v 2 coefficients of the corresponding N coll values [48]. The relative difference of the obtained D 0 -meson v 2 with respect to the one obtained with the arithmetic average is negligible for p T < 8 GeV/c and its average is 19% for p T > 8 GeV/c, which is still compatible within uncertainties. The v e←c 2 coefficients obtained with the two approaches show a relative difference of 2% in the range 2 < p T < 8 GeV/c. This deviation is considered as a consequence of statistical fluctuations in the D 0 -meson v 2 measurement for p T > 8 GeV/c and thus no systematic uncertainty is assigned for this effect.
In order to investigate the impact of the assumption of the particle mass ordering of the elliptic flow [51] used to determine the systematic uncertainty from the B feed-down subtraction in the D 0meson v 2 measurement, one can assume that the v 2 of prompt D mesons from B-meson decays should be between zero and the unidentified charged particle v 2 . The unidentified charged particle v 2 in 20-40% Pb-Pb collisions is obtained by the average of the v 2 measurements in the 20-30% and 30-40% centrality classes [50] weighted by the corresponding N coll values. The v 2 coefficients of prompt D 0 mesons and unidentified charged particles are compatible within uncertainties as well as the v e←c 2 obtained with these two results. Therefore, the lower limit of the systematic uncertainty from the B feed-down subtraction can be positioned at the central values of the prompt D-meson v 2 and v e←c 2 without strictly considering that the B-meson v 2 is expected to be lower than the D-meson v 2 .
As a consequence of the p T interval (2 < p T < 16 GeV/c) of the D 0 -meson v 2 and p T -differential yield measurements, the v e←c 2 is obtained in the range 2 < p T < 8 GeV/c. The fraction of electrons with p T > 2 GeV/c that come from D 0 mesons with p T < 2 GeV/c is negligible according to PYTHIA simulations. The effect of the p T upper limit of the D 0 -meson measurements is studied by evaluating the v e←c 2 with extrapolation of the p T and v 2 distributions of D 0 mesons up to 26 GeV/c. The transverse momentum extrapolation is obtained from the power-law fit function shown in the left panel of Fig. 2, while the impact of the v 2 of D 0 mesons is estimated by explicitly setting its value, in the interval 16 < p T < 26 GeV/c, to either zero, or constant at high p T , or maximum value of the prompt D 0 -meson v 2 (shown in Fig. 2). The highest relative difference in these three scenarios, which increases from 0.3% to 40% in the interval 2 < p T < 8 GeV/c, is assigned as a conservative systematic uncertainty.
The effect of the mid-rapidity range of D 0 mesons is investigated by obtaining the v e←c 2 using the D 0 -meson p T distribution in the rapidity range |y| < 1.6 as input for the simulation. The D 0 -meson v 2 is considered to be the same in this rapidity range, because no dependence on rapidity was observed in ALICE results on leptons from heavy-flavour hadron decays [24,37] as discussed previously. The relative difference of the obtained v e←c 2 with respect to the one using the D 0 -meson p T distribution in the rapidity range |y| < 0.8 is negligible and thus no additional systematic uncertainty is considered due to the rapidity effect.
In this analysis, the v 2 and shape of the p T -differential yields of charm hadrons are assumed to be the same as the ones measured for D 0 mesons. This is justified by the fact that the v 2 coefficients of D 0 , D + and D * + mesons are compatible within uncertainties in 30-50% Pb-Pb collisions at √ s NN = 2.76 TeV [31], also the prompt D + s v 2 is compatible within uncertainties with the prompt non-strange D meson v 2 in 30-50% Pb-Pb collisions at √ s NN = 5.02 TeV [34]. In addition, the ratios of the yields of D + /D 0 and D * + /D 0 were observed to be constant within uncertainties in pp collisions at √ s = 7 TeV and no modification of the ratios was observed within uncertainties in central and semi-central Pb-Pb collisions at √ s NN = 2.76 TeV [52]. A possible hint for an enhancement of the D + s /D 0 ratio is observed in 0-10% Pb-Pb collisions at √ s NN = 2.76 TeV [53], but the current uncertainties do not allow for a conclusion. The effect of different decay kinematics of charm particles is estimated by simulating the v 2 of electrons from combined D, D * , D s , and Λ c particle decays taking into account the fraction of charm quarks that hadronise into these particles [54] and using the same simulation input as used in the analysis (p T -differential yield and v 2 of D 0 mesons). The obtained v e←c 2 is compatible with the one using D 0 -meson decay and thus no systematic uncertainty is considered due to this effect. In order to exemplify the impact of a possible production enhancement of D + s and Λ c particles in Pb-Pb collisions with respect to pp collisions, their fragmentation fractions are increased by a factor 2 and 5, respectively, in the simulation of the combined charm meson v 2 . The relative difference of the obtained v e←c 2 and the one using D 0 -meson decay is negligible for p T < 3 GeV/c and its average is 5% for p T > 3 GeV/c.
Finally, the D 0 -meson v 2 systematic uncertainty from data is summed in quadrature with other sources of systematic uncertainty that affect significantly the v e←c 2 estimation, which are the p T distribution of D 0 mesons and the limited p T interval of the D 0 -meson measurements. They are considered as uncorrelated since the effect from the p T distribution of D 0 mesons is obtained with the BAMPS result and the effect from the limited p T interval of the D 0 -meson measurements is obtained by extrapolations. The term "from data" is maintained later in this paper to distinguish all sources of systematic uncertainty from the systematic uncertainty related to the B feed-down subtraction of the D 0 -meson v 2 measurement, which is shown separately.

Elliptic flow of electrons from beauty-hadron decays
The v 2 of electrons from beauty-hadron decays (v e←b 2 ) is obtained from Eq. 3 using the R, v e←c+b 2 and v e←c 2 results presented in their respective sections. The three results are considered as statistically independent. First, the factor R was measured in a different collision system (pp collisions) or obtained with calculations. Second, the v e←c 2 is obtained with a simulation using measurements of D 0 mesons reconstructed via the hadronic decay channel D 0 → K − π + in a different centrality class than in the v e←c+b 2 measurement. Even though the systematic uncertainties of the R, v e←c+b 2 and v e←c 2 results might be partially correlated, especially concerning the particle identification selection criteria, the limited public information prevents a more accurate treatment of these uncertainties. Therefore, they are assumed to be uncorrelated as a conservative estimation. As an example of the effect of a possible overestimation, if the systematic uncertainties of the v e←c 2 and v e←c+b 2 results decrease by 30% in the interval 2 < p T < 8 GeV/c, the systematic uncertainty from data of the v e←b 2 result is expected to decrease by approximately 24%. Therefore, the statistical and systematic uncertainties of the R, v e←c+b 2 and v e←c 2 results are propagated as independent variables. The v e←b 2 systematic uncertainties from data and from the B feed-down subtraction are asymmetric as a consequence of the systematic uncertainty asymmetry of the measurements used in this analysis. The systematic uncertainty from data is evaluated according to the method described in [55], where the positive and negative deviations are obtained separately and their average is added in quadrature. For verification, the alternative approach presented in [56] is also applied in this analysis. No significant difference between these methods is observed. Since the asymmetry of the systematic uncertainty from the B feed-down subtraction only comes from the v e←c

Results
The relative contribution of electrons from beauty-hadron decays to the yield of electrons from heavyflavour hadron decays at √ s = 2.76 TeV obtained with POWHEG+ PYTHIA at NLO accuracy using EPS09 NLO nPDFs is shown in the left panel of Fig. 3. The result is compared with the R in pp collisions at √ s = 2.76 TeV measured by the ALICE collaboration using the electron-hadron azimuthal correlation technique [39,40] and with the BAMPS result in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV [41]. The comparison shows that R is higher when in-medium effects are present, which is consistent with the expectation of the mass hierarchy of the energy loss of charm and beauty quarks in the medium. The right panel of Fig. 3 shows the v 2 of electrons from charm-hadron decays at mid-rapidity in 20-40% 76 TeV obtained with POWHEG+PYTHIA at NLO accuracy using EPS09 NLO nPDFs. The result is compared with the R in pp collisions at √ s = 2.76 TeV measured by the ALICE collaboration using the electron-hadron azimuthal correlation technique [39,40] and with the BAMPS result in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV [41]. The statistical and systematic uncertainties of the R coefficient obtained with POWHEG+PYTHIA and from the BAMPS model are zero. Right: Elliptic flow of electrons from charm-hadron decays at mid-rapidity in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV estimated using a MC simulation with PYTHIA based on ALICE results [19,32]. The vertical error bars represent the statistical uncertainties and the horizontal error bars indicate the bin widths. The empty and filled boxes represent the systematic uncertainties from data and from the B feed-down subtraction, respectively, in the D 0 -meson v 2 measurement [32].
Pb-Pb collisions at √ s NN = 2.76 TeV obtained with a MC simulation with PYTHIA using as input the p T -differential yield and v 2 distributions of D 0 mesons in the same collision system. A positive v 2 of electrons from charm-hadron decays is found in all p T intervals, with a maximum significance of 3.2σ − , where σ − is the combined statistical and systematic uncertainties of the lower limit, in the interval 2 < p T < 3 GeV/c.
The v 2 coefficients of electrons from beauty-hadron decays in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV obtained with different approaches of the factor R (left panel of Fig. 3) are shown in Fig. 4. The result computed with the coefficient R in pp collisions is an estimation of the minimum value, as discussed previously. The v 2 of electrons from beauty-hadron decays in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV is compatible with zero within approximately 1σ of the total uncertainty, obtained by summing in quadrature the different uncertainty contributions, in all p T intervals and different R coefficients. However, the large statistical and systematic uncertainties prevent a definite conclusion. The result is consistent with the measured v 2 of non-prompt J/ψ mesons from B-hadron decays in 10-60% Pb-Pb collisions at √ s NN = 2.76 TeV [27], which is also compatible with zero within uncertainties. Figure 5 shows the v 2 of electrons from charm-and beauty-hadron decays, inclusive [37] and separated, in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV. The v 2 of electrons from beauty-hadron decays is lower than the v 2 of electrons from charm-hadron decays, although they are compatible within uncertainties. The average of the v 2 coefficients of electrons from charm-and beauty-hadron decays obtained in the interval 2 < p T < 8 GeV/c are listed in Table 1. Because of the asymmetric uncertainties, the average is obtained numerically with an iterative sum of the likelihood functions parametrised by variable-width Gaussians [56,57]. The standard deviation, which is the combination of statistical and systematic uncertainties, is assumed to vary linearly. The maximum value of the summed likelihood function corresponds to the average v 2 , while the points at which the function is -0.5 correspond to the lower and upper limits of the total uncertainty.

Conclusions
We presented a method to subtract the charm quark contribution in the elliptic flow of electrons from heavy-flavour hadron decays. The v 2 of electrons from charm-hadron decays was estimated using a MC simulation of D 0 -meson decays into electrons with PYTHIA, based on measurements of the p Tdifferential yield and v 2 of D 0 mesons in Pb-Pb collisions at √ s NN = 2.76 TeV by ALICE. A positive v 2 of electrons from charm-hadron decays is found with a maximum significance of 3.2σ − in the interval 2 < p T < 3 GeV/c. The computed v 2 of electrons from charm-hadron decays was then subtracted from the v 2 of electrons from heavy-flavour hadron decays in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV measured by the ALICE collaboration. The subtraction was weighted by the relative contribution of electrons from beauty-hadron decays to the yield of electrons from heavy-flavour hadron decays. Since this observable is not measured in Pb-Pb collisions, three different approaches were used as estimations in the analysis. The resulting v 2 of electrons from beauty-hadron decays in 20-40% Pb-Pb collisions at √ s NN = 2.76 TeV from the subtraction is compatible with zero within approximately 1σ of the total uncertainty in all p T intervals and different approaches of the relative contribution of electrons from beauty-hadron decays to the yield of electrons from heavy-flavour hadron decays. However, the large statistical and systematic uncertainties prevent a definite conclusion. The v 2 of electrons from beautyhadron decays is found to be lower than the v 2 of electrons from charm-hadron decays.

Outlook
In the presented method, the elliptic flow of electrons from beauty-hadron decays can be determined by using three observables that have largely been measured at the LHC and RHIC. Based on available results of these observables, the procedure was applied using measurements performed by the ALICE collaboration. The method demonstrated to be effective; however, the current statistical and systematic uncertainties of the ALICE results prevent a definite conclusion whether the collective motion of the medium constituents influences beauty quarks. A better accuracy of the results on heavy-flavour particles has been achieved in measurements in Pb-Pb collisions at √ s NN = 5.02 TeV [34] and it is expected to be further improved with the ALICE upgrade, which is foreseen to start in 2019.
In particular, the upgrade of the Inner Tracking System (ITS) detector will improve the determination of the distance of closest approach to the primary vertex, momentum resolution and readout rate capabilities [58]. These improvements will allow for more precise measurements of D mesons down to low transverse momenta and for reducing the systematic uncertainties from data and from the B feeddown subtraction. The latter will be possible with the direct measurement of the fraction of prompt D mesons and D mesons from B-meson decays, which is expected to be accessible with relative statistical and systematic uncertainties smaller than 1% and 5% [58], respectively, for prompt D 0 mesons. In addition, the ITS upgrade will enable the tracking of electrons down to approximately 0.05 GeV/c and enhance the capability to separate prompt from displaced electrons [58], improving the reconstruction of electrons that do not originate from heavy-flavour hadron decays needed for the background subtraction. Moreover, the systematic uncertainty of the elliptic flow of electrons from beauty-hadron decays can be further improved by taking into account correlations among different contributions.
The capability of the heavy-flavour measurements will also enhance with the increase of luminosity. For instance, the current relative statistical uncertainty of the D-meson v 2 measurement in Pb-Pb collisions is 10% for an integrated luminosity of 0.1 nb −1 , while it is expected to be 0.2% for a scenario with an integrated luminosity of 10 nb −1 [58]. Also the elliptic flow coefficients of D s and Λ c particles are expected to be achievable with a relative statistical uncertainty of 8% and 20% [58], respectively. Therefore, the presented method can be used to extract the elliptic flow of electrons from beautyhadron decays with better precision with future measurements of the three needed observables.