ALICE Collaboration††thanks: See Appendix A for the list of collaboration members
The production of beauty hadrons was measured via semi-leptonic decays at mid-rapidity with the ALICE detector at the LHC in the transverse momentum interval in minimum-bias p–Pb collisions at = 5.02 TeV and in in the 20% most central Pb–Pb collisions at = 2.76 TeV. The pp reference spectra at = 5.02 TeV and = 2.76 TeV, needed for the calculation of the nuclear modification factors and , were obtained by a pQCD-driven scaling of the cross section of electrons from beauty-hadron decays measured at = 7 TeV. In the interval a suppression of the yield of electrons from beauty-hadron decays is observed in Pb–Pb compared to pp collisions. Towards lower , the values increase with large systematic uncertainties. The is consistent with unity within systematic uncertainties and is well described by theoretical calculations that include cold nuclear matter effects in p–Pb collisions. The measured and these calculations indicate that cold nuclear matter effects are small at high transverse momentum also in Pb–Pb collisions. Therefore, the observed reduction of below unity at high may be ascribed to an effect of the hot and dense medium formed in Pb–Pb collisions.
In collisions of heavy nuclei at ultra-relativistic energies, a high-density colour-deconfined state of strongly-interacting matter, called Quark–Gluon Plasma (QGP), is expected to be produced [1, 2]. Due to their large masses (), heavy quarks (charm and beauty) are almost exclusively produced in the early stage of the collision via hard parton scatterings characterised by production-time scales of less than 0.1 and 0.01 fm/ for charm and beauty quarks, respectively . They can, therefore, serve as probes to test the mechanisms of medium-induced parton energy loss, because the formation time of the QGP medium is expected to be about 0.3 fm/  and its decoupling time is about 10 fm/ for collisions at LHC energies . Due to their stronger colour coupling to the medium gluons are argued to lose more energy than quarks [6, 7, 8]. Furthermore, the radiative energy loss of heavy quarks is predicted to be reduced with respect to light quarks due to the mass-dependent restriction of the phase space into which medium-induced gluon radiation can take place (dead-cone effect) [9, 10, 11, 12]. The effect of the charm-quark mass on energy loss becomes negligible at high transverse momentum, 10 , where the ratio approaches zero . Therefore, due to the larger mass, beauty quarks can be sensitive probes for testing the mass dependence of the parton energy loss up to transverse momenta well above 10 . Final-state effects, such as colour-charge and mass dependence of parton energy loss, can be studied experimentally through the spectra of hadrons containing heavy quarks in comparison with light-flavour hadrons in heavy-ion (AA) collisions.
The understanding of final-state effects requires measurements of initial-state effects in Cold Nuclear Matter (CNM), which are inherent to nuclei in the collision system and thus present in AA collisions. Measurements in proton–nucleus (p–A) collisions are used to investigate cold nuclear matter effects such as the modification of the Parton Distribution Functions (PDF) inside the nucleus with respect to those in the proton, broadening via parton collisions inside the nucleus prior to the hard scattering and energy loss in cold nuclear matter [14, 15, 16, 17, 18]. The effects of hot (cold) nuclear matter can be studied using the nuclear modification factor, (), defined as the ratio of the distributions measured in AA (p–A) collisions with respect to the one in pp collisions:
where and are the -differential yield and production cross section of a given particle species in AA and pp collisions, respectively, and is the average of the nuclear overlap function for the centrality range under study .
Previous beauty-hadron production measurements in pp collisions at various energies at RHIC [20, 21], the Tevatron  and the LHC [3, 23, 24, 25, 26, 27, 28] are described by Fixed Order plus Next-to-Leading-Log perturbative Quantum Chromodynamics (FONLL pQCD) calculations [29, 30, 31] within uncertainties.
At both RHIC and the LHC, a suppression of the yield of D mesons and high- electrons and muons from heavy-flavour hadron decays was observed in AA collisions. The suppression is nearly as large as that of light-flavour hadrons at high [32, 33, 34, 35, 36]. The D meson and pion were found to be consistent within uncertainties and described by model calculations that include a colour-charge dependent energy loss [34, 37, 38]. However, in addition to energy loss, the nuclear modification factor is also influenced by e.g. the parton spectrum and the fragmentation into hadrons [39, 13]. Furthermore, the nuclear modification factors of prompt D mesons and of J/ from B meson decays were compared in the interval for D mesons and for J/ mesons in order to have a similar average (10 ) for the heavy hadrons [34, 40, 41]. This comparison with models indicates that charm quarks lose more energy than beauty quarks in this interval in central Pb–Pb collisions. The b-jet yield as measured in Pb–Pb collisions also shows a suppression compared with the yield expected from pp collisions in the jet- interval . Recently, the relative contributions of electrons from charm- and beauty-hadron decays were measured as a function of transverse momentum in Au–Au collisions at RHIC . There is a hint that in the momentum interval the of electrons from beauty-hadron decays is larger than that of electrons from charm-hadron decays.
In p–Pb collisions at the LHC, the nuclear modification factors of B mesons , b-jets , J/ from beauty-hadron decays [46, 47], leptons from heavy-flavour hadron decays and D mesons [48, 49] were investigated extensively. The results are consistent with unity within uncertainties and compatible with theoretical calculations including cold nuclear matter effects [45, 47, 46, 48]. Therefore, the observed suppression of charm and beauty yields at high in Pb–Pb collisions is not explained in terms of initial-state effects but is due to strong final-state effects induced by hot partonic matter.
In central d–Au collisions at = 200 GeV at RHIC, an enhancement was measured at backward rapidity by means of of muons from heavy-flavour hadron decays . Theoretical calculations including modified PDFs cannot describe the data, implying that models incorporating only initial-state effects are not sufficient and suggesting the possible importance of final-state effects in the d–Au collision system. Recently, a potential signature of collective behaviour in small systems was observed via the anisotropic flow parameter of charged hadrons in p–Pb collisions [51, 52, 53, 54] and in d–Au collisions[55, 56], suggesting radial flow as a possible explanation of the enhancement of the .
In this paper, the invariant cross section in p–Pb and yield in Pb–Pb collisions are presented together with the nuclear modification factors, and , of electrons from beauty-hadron decays in p–Pb and Pb–Pb collisions at = 5.02 TeV and = 2.76 TeV, respectively. The identification of electrons from beauty-hadron decays is based on their separation from the interaction vertex, induced by the sizable lifetime of beauty hadrons. The p–Pb (Pb–Pb) measurement covers the rapidity range () and the interval (). In the p–Pb collisions, due to the different energy per nucleon of the proton and lead beam, the centre-of-mass system (cms) is shifted by in the proton beam direction, resulting in the rapidity coverage for electrons. Given the cms energies and the rapidity coverages in the p–Pb and Pb–Pb collisions, both measurements probe, at the lowest , similar values of Bjorken- of about for electrons from beauty-hadron decays . The Pb–Pb measurement is restricted to the 20% most central Pb–Pb collisions, where the largest effect of energy loss on heavy-flavour production is expected.
The paper is organised as follows: Section 2 describes the experimental apparatus and the data samples used in both analyses, which are outlined in Section 3. Details of the analysis in p–Pb and Pb–Pb collisions are given in Sections 4 and 5, respectively. The determination of the pp reference spectra for the calculations of the and is reported in Section 6. The results are presented and discussed in Section 7. Section 8 summarises the results.
2 Experimental apparatus and data samples
A comprehensive description of the ALICE apparatus and its performance can be found in  and , respectively. Electron tracks were reconstructed and identified using detectors located inside the solenoid magnet that generates a field of 0.5 T parallel to the beam direction. Forward and backward detectors inside and outside the magnet were employed for triggering, background rejection and event characterisation.
Charged particles are tracked with the Inner Tracking System (ITS) [59, 61] and the Time Projection Chamber (TPC)  in the pseudorapidity range 0.9. The ITS consists of six cylindrical layers of silicon detectors. The two innermost layers are made of Silicon Pixel Detectors (SPD), the two middle layers of Silicon Drift Detectors (SDD) and the two outermost layers of Silicon Strip Detectors (SSD). In the direction perpendicular to the detector surface, the total material budget of the ITS corresponds on average to 7.7% of a radiation length . In this analysis, the ITS was also used to reconstruct the primary (interaction) vertex and the track impact parameter , defined as the distance of closest approach of the track to the interaction vertex in the plane transverse to the beam direction. The resolution on is better than 65 m and 70 m for charged particles with momenta larger than 1 in Pb–Pb and p–Pb collisions , respectively, including the resolution of the primary vertex determination. The particle identification capability of the four outer layers of the ITS via the measurement of the ionisation energy loss was used at low transverse momentum in the p–Pb analysis. The TPC, which provides up to 159 space points per track, is used for particle identification via the measurement of the specific energy loss in the detector gas. The tracks reconstructed in the ITS and the TPC are matched to hits in the other detectors inside the magnet located at larger radii. The Transition Radiation Detector (TRD)  surrounding the TPC provides hadron and electron identification via the measurement of the specific energy loss and transition radiation. During the Pb–Pb (p–Pb) data taking period it covered 7/18 (13/18) of the full azimuth. Therefore, only in the Pb–Pb analysis it was used to verify the amount of hadron contamination within the electron identification strategy at low transverse momentum (see Section 5). The Time-Of-Flight array (TOF) , based on Multi-gap Resistive Plate Chambers (MRPCs), provides hadron rejection at low transverse momentum via the time-of-flight measurement, within the electron identification strategy applied in both analyses. The T0 detectors, arrays of Cherenkov counters, located at 350 cm and 70 cm from the interaction point along the beam direction  provided, together with the TOF detector, the precise start time for the time-of-flight measurement in the p–Pb analysis. For central Pb–Pb events the start time was estimated only using the particle arrival times at the TOF detector.
The SPD, the T0 detectors as well as the V0 scintillator arrays, placed on both sides of the interaction point at 2.8 5.1 (V0-A) and 3.7 1.7 (V0-C), respectively, can be employed to define a minimum-bias trigger. The two Zero Degree Calorimeters (ZDC), that are symmetrically located 112.5 m from the interaction point on either side, were used in the offline event selection to reject beam-gas interactions by correlating the time information with the one from the V0 detectors.
The Pb–Pb and p–Pb data presented here were recorded in 2010 and 2013, respectively. Minimum-bias p–Pb collisions were selected by requiring coincident signals in V0-A and V0-C (V0AND condition). Beam-gas interactions were rejected offline by the aforementioned correlation of the ZDC and V0 time information. The Pb–Pb collisions were collected with two different minimum-bias interaction triggers. The first trigger condition required signals in two of the following three detectors: SPD (two hits in the outer SPD layer), V0-A and V0-C. The second trigger condition required a coincidence between V0-A and V0-C. Both minimum-bias trigger conditions had efficiencies larger than 95% for hadronic interactions, whereas the second rejected electromagnetic processes to a large extent . Only events with a primary vertex within 10 cm from the centre of the detector along the beam direction were considered in the p–Pb and Pb–Pb analyses. The Pb–Pb events were categorised into centrality classes by fitting the sum of the two V0 signal amplitudes with a geometrical Glauber-model simulation , as described in . The Glauber-model simulation yields a value of 18.93 0.74 for the average nuclear overlap function for the 20% most central Pb–Pb collisions considered in the analysis. About 100 and 3 million p–Pb and 20% most central Pb–Pb events passed the offline selection criteria corresponding to an integrated luminosity of b and b, respectively.
3 Analysis overview and electrons from background sources
The identification of electrons from beauty-hadron decays is divided into the following steps:
selection of tracks with good quality,
electron identification (eID),
determination of the electron yield from beauty-hadron decays.
The signal contains both electrons from direct decays (, branching ratio: 11%) as well as cascade decays (, branching ratio: 10%) of hadrons that contain a beauty (or anti-beauty) quark . Throughout the paper the term ‘electron’ denotes both electron and positron. The track selection procedure is identical to previous analyses on the production of electrons from beauty-hadron decays [23, 24]. The selection criteria are the same in the p–Pb and Pb–Pb analyses, except for the restriction of the geometrical acceptance in rapidity, which was adjusted in each collision system to the region where the TPC could provide optimal electron identification, taking into account the detector and running conditions during each data-taking period. In Pb–Pb collisions this corresponds to the rapidity range and in p–Pb to . The tracks were required to have associated hits in both SPD layers, in order to minimise the contribution of electrons from photon conversions in the ITS detector material and the fraction of tracks with misassociated hits (see below).
The electrons were identified with the TPC and the TOF detectors via the measurement of their respective signal, specific energy loss in the gas (d/d) and the time-of-flight. The selection variable (hereafter or ) is defined as the deviation of the measured signal of a track with respect to the expected signal for an electron in units of the corresponding detector resolution ( or ). The expected signal and the resolution originate from parametrisations of the TPC and TOF detector signals, described in detail in . For both analyses, particles were accepted with the TPC as electron candidates if they satisfied the condition . This asymmetric selection was chosen to remove hadrons, that are mainly found at negative values. However, at low and high transverse momentum, the eID strategy based on TPC is subject to contamination from pions, kaons, protons and deuterons. To resolve these ambiguities, a selection cut of was applied for the whole range in the Pb–Pb analysis and for 2.5 in the p–Pb analysis. The remaining hadron contamination was determined via data-driven methods in the p–Pb analysis and subtracted statistically (see Section 4). The technique used for the Pb–Pb analysis is described in Section 5.
The electrons passing the track and eID selection criteria originate, besides from beauty-hadron decays, from the following background sources. In what follows, prompt and non-prompt contributions are marked in parentheses as ‘P’ and ‘NP’, respectively:
(P) Dalitz and di-electron decays of prompt light neutral mesons (),
(P) di-electron decays of prompt heavy quarkonia (J/, etc.).
(NP) decay chains of hadrons carrying a strange (or anti-strange) quark,
(NP) photon conversions in the detector material,
(NP) semi-leptonic decays of prompt hadrons carrying a charm (or anti-charm) quark.
The measurement of the production of electrons from beauty-hadron decays exploits their larger mean proper decay length 500 m ) compared to that of charm hadrons and most other background sources, resulting in a larger average impact parameter. The sign of the impact parameter value is attributed based on the relative position of the track and the primary vertex, i.e. if the primary vertex lies on the left- or right-hand side of the track with respect to the particle momentum direction in the transverse plane.
For the presented analyses, the impact parameter was multiplied with the sign of the particle charge and of the magnetic field component along the beam axis (plus or minus for the two field orientations). With this definition, the sign of the impact parameter depends on whether the primary vertex lies inside or outside of the circle defined by the track projection in the transverse plane. Electrons from the conversion of photons in the detector material have an initial momentum with a very small angle to the direction of the photon. The magnetic field bends the track away from the primary vertex. Thus, they typically have an impact parameter . The asymmetric shape helps to differentiate this background source. It is important to include the field configuration, because the magnetic field direction was reversed during the Pb–Pb data taking period, which motivated this redefinition.
Figure 1 shows for two intervals the resulting distribution of the measured impact parameter value multiplied by the sign of the charge of each track and the sign of the magnetic field for the interval in the 20% most central Pb–Pb collisions. The impact parameter distributions for electrons from beauty- and charm-hadron decays, from Dalitz decays of light mesons, and from photon conversions are also drawn for comparison. The distributions were obtained from Monte Carlo simulations and normalised to the data using the fit values described in Section 5. The distribution for electrons from photon conversions is, as explained before, visible as an asymmetric and shifted distribution. The impact parameter distribution of electrons from prompt sources, such as Dalitz and quarkonium decays, is determined by the impact parameter resolution. The electrons from these sources are thus categorised as Dalitz decays within both analyses.
The Monte Carlo simulations were produced as follows. A sample of minimum-bias Pb–Pb collisions at = 2.76 TeV was generated with HIJING v1.36  for efficiency and acceptance corrections as well as to obtain the impact parameter distributions for photon conversions and Dalitz decays. To increase the statistics of electrons from charm- and beauty-hadron decays, a signal enhanced sample was generated using pp events produced by the generator PYTHIA v6.4.21  with Perugia-0 tune . Each added pp event contains one or pair. For the p–Pb analysis, the same procedure was used. The generated particles were propagated through the ALICE apparatus using GEANT3  and a realistic detector response was applied to reproduce the performance of the detector system during data taking.
The inclusive yield of electrons originating from strange-hadron decays is small compared to the other background sources. However, as these electrons originate from secondary from strange-hadron decays (, , , ) and three prong decays of strange hadrons (,), the impact parameter distribution is broader than that of electrons from Dalitz and di-electron decays of other light neutral mesons. Sections 4 and 5 describe how the analyses handle this background contribution.
Although requiring hits in both SPD layers, electrons from photon conversions in detector material with production radii outside the SPD layers were observed to have passed the track selection. These electron tracks are wrongly associated with signals of other particles in the inner detector layers. Within this paper these electrons are called ‘mismatched conversions’. The amount of mismatched conversions depends on the track multiplicity within the event and thus has a larger impact for the Pb–Pb analysis. Sections 4 and 5 outline how the analyses deal with the mismatched conversions.
The impact parameter distributions of electrons from most background sources are narrow compared to the one of electrons from beauty-hadron decays. By applying a minimum cut on the absolute value of the impact parameter , the fraction of electrons from beauty-hadron decays can thus be enhanced. The remaining background can be described using a cocktail method and subtracted statistically to obtain electrons from beauty-hadron decays [23, 24]. This method was applied in the p–Pb analysis and is described in detail in Section 4. Another technique, used in the Pb–Pb analysis (see Section 5), is to make use of the whole impact parameter distribution, i.e. to compare the impact parameter distributions of the various electron sources from simulation (templates) with the impact parameter distribution of all measured electron candidates to estimate the individual contributions.
4 Data analysis in p–Pb collisions
The identification of electrons from beauty-hadron decays in the p–Pb analysis is based on the selection of electrons with large impact parameters. This method was already applied in pp collisions at = 2.76 TeV and = 7 TeV [24, 23]. Since the impact parameter distribution of electrons from beauty-hadron decays is broader compared to the one of electrons from most background sources (see Section 3), the requirement of a minimum absolute impact parameter enhances the signal-to-background (S/B) ratio of electrons from beauty-hadron decays. The remaining background due to hadron contamination and electrons from background sources was obtained via a data-driven method and from Monte Carlo simulations re-weighted to match the distributions of the background sources in data, respectively, and then subtracted.
4.1 Extraction of electrons from beauty-hadron decays
Electron candidates with an impact parameter (with in cm and in ) were selected. This selection criterion was determined from Monte Carlo simulations to maximise the significance for electrons from beauty-hadron decays. The selection of the minimum impact parameter is dependent, because the width of the impact parameter distribution, the S/B ratio as well as the true impact parameter distribution of the various electron sources  are dependent.
The number of hadrons passing the track selection, eID, and the minimum impact parameter requirement was estimated at high transverse momentum ( ) by parametrising the TPC distribution in momentum slices, and it was subtracted . Above a of , the hadron contamination increases with transverse momentum and reaches 10% at 8 , see Fig. 2 (left). At low transverse momentum ( ), the hadron contamination is negligible except in the transverse momentum interval , see Fig. 2 (left), where electrons cannot be distinguished from protons via the measurement of specific energy loss in the TPC gas. In addition, the requirement of a minimum impact parameter increases the relative contribution of secondary protons originating from e.g. and decays, which have larger impact parameter values compared to electrons from beauty-hadron decays. The relative abundance of protons in the electron candidate sample was determined by using the ITS particle identification capabilities, because electrons and protons can be separated with ITS in this momentum interval. The ITS energy loss signal was fitted with data-derived templates for electrons and protons. The templates were obtained in -bins by selecting electrons and protons with tight selection criteria in TOF and TPC. The estimated proton contribution, which is 10% (4%) in the interval (), was subtracted statistically from the measured electron candidate distribution.
Figure 2 (left) shows the transverse momentum distribution of electrons passing the track, eID, and impact parameter selection, before efficiency corrections. The contributions due to the proton and hadron contamination at low and high , respectively, determined via the aforementioned methods are shown. Also shown are the distributions of electrons originating from the various background sources, which were obtained using the Monte Carlo simulations described in Section 3. To match the measured shapes, the differential yields of the background sources were re-weighted in the Monte Carlo simulations prior to the propagation through the ALICE apparatus with GEANT3. As there is no measurement of the production cross section in p–Pb collisions available, the input was based on the measured charged-pion spectra [73, 74] assuming . Due to the requirement of a minimum impact parameter, the contribution of electrons from decays of secondary from strange-hadron decays is comparable with the one from primary decays. Therefore the measured spectra of , and  were used to compute the corresponding weights. To obtain the weights, the pion and strange-hadron spectra were parameterised with a Tsallis function as described in . The contribution of electrons originating from secondary pions from strange-hadron decays or three-body decays of strange hadrons is shown in Fig. 2 (left). The other light mesons (, , , and ), which contribute little, via Dalitz decays and photon conversions compared to primary decays, were re-weighted via -scaling of the spectrum . The electron background from charm-hadron decays was estimated based on the , and meson production cross section measurements with ALICE  in the transverse momentum intervals , and , respectively. In a first step the measurements were extrapolated to the interval by assuming constant ratios / and / from the measured D meson production cross sections. Next the differential production cross sections were extrapolated to = 50 via FONLL pQCD calculations. About 10% of the electrons with 8 originate from the extrapolated D meson high- region ( ). The electron contribution from decays was estimated using the ratio / measured by the ZEUS Collaboration . Analogous to the light mesons, the measured D meson spectra were also used to re-weight the distributions in the Monte Carlo simulations.
The signal of electrons from beauty-hadron decays was obtained after subtraction of the aforementioned background contributions from the measured electron candidate sample after track selection, eID and impact parameter requirement. The resulting spectrum is shown in Fig. 2 (left). At , the number of electrons from beauty-hadron decays is approximately equal to the one from charm-hadron decays, from Dalitz decays of light mesons, from strange-hadron decays and from photon conversions, resulting in a S/B ratio of approximately . With increasing the background electron yield from Dalitz decays of light mesons, from strange-hadron decays and from photon conversions quickly decreases compared to the contribution of electrons from charm-hadron decays. In the interval , the S/B ratio reaches its maximum of 3. Here the electron background mostly originates from charm-hadron decays. At higher , the S/B ratio decreases again due to the increasing hadron contamination. Other background sources, such as di-lepton decays of J/ mesons are negligible due to the minimum impact parameter selection. The yield of electrons from Drell-Yan processes is negligible over the whole range .
The raw yield of electrons from beauty-hadron decays was then corrected for the geometrical acceptance and for the efficiency () of the track reconstruction, matching and selection criteria, TOF electron identification and minimum impact parameter requirement using the Monte Carlo simulations. The efficiency of the TPC electron identification selection () was determined to be 69% via a data-driven approach based on the distributions . The transverse-momentum dependence of the efficiencies is shown in Fig. 2 (right). The total efficiency shows a significant dependence, mainly due to the cut. The effects of the finite momentum resolution and the energy loss due to Bremsstrahlung were taken into account in a bin-by-bin resolution correction step based on a Monte Carlo simulation [76, 23].
The -differential invariant cross section of electrons from beauty-hadron decays, , is thus given as:
where is the centre of the bin with width and denotes the geometrical acceptance in to which the analysis was restricted. is the total number of analysed minimum-bias events. The p–Pb cross section for the minimum-bias V0 trigger condition, which has an efficiency of more than 99% for non-single-diffractive (NSD) p–Pb collisions , is b .
4.2 Systematic uncertainties estimation
An overview of the relative systematic uncertainties is shown in Table 1. The systematic uncertainties were estimated as a function of by repeating the analysis with modified track selection and eID criteria and by varying the background yields within their estimated uncertainties.
The uncertainty of the tracking results from differences in data and Monte Carlo simulations for the track reconstruction with the ITS and the TPC, which includes the uncertainty of finding a hit in the ITS for a track reconstructed in the TPC. The latter uncertainty (3%) was taken from , where the effect was studied for charged particles. The TOF-TPC matching uncertainty (5%) was obtained by comparing the matching efficiency of electrons from photon conversions identified via topological selections in data and Monte Carlo simulations. The TOF eID uncertainty was derived by repeating the analysis with different eID selection criteria. At high the TOF was not used in the analysis and thus the corresponding uncertainty does not apply in this region. The uncertainty of the TPC eID was estimated in the same way as for the TOF eID. The systematic uncertainty of the determination of the hadron contamination ranges from 1% to 6%, i.e. increasing as the contamination itself with increasing .
The systematic uncertainty of the minimum impact parameter requirement was evaluated by varying this selection criterion by 1 , where corresponds to the measured impact parameter resolution . At 1 (8 ) this corresponds to a 10% (%) variation of the cut value.
The number of electrons from photon conversions increases quickly with decreasing transverse momentum (see Fig. 2, left). The difference in yield of mismatched conversions in data and Monte Carlo simulations was estimated and assigned as a systematic uncertainty. For this purpose pions from decays identified via topological and invariant mass cuts  can be used, because their decay vertex can be reconstructed, in contrast to electrons from photon conversions, for which it is more difficult due to their small opening angle. The yield of pions from decays was studied as a function of the production vertex with and without requiring a signal in both SPD layers and compared with the corresponding results from Monte Carlo simulations. The difference in yield was propagated into the simulation by renormalising the number of electrons from photon conversions. Repeating the full analysis with the varied conversion yield results in the uncertainties listed in Table 1.
The systematic uncertainty arising from the subtraction of electrons from the various background sources was evaluated by propagating the statistical and systematic uncertainties of the light-meson, strange- and charm-hadron measurements used as input to re-weight the distributions in Monte Carlo simulations. Uncertainties due to the -scaling of the background yields, estimated as 30% , and the extrapolation of the D meson distributions to the unmeasured transverse momentum regions were included. The latter was obtained by using the uncertainties of the various D meson ratios and by using a power-law fit instead of FONLL pQCD calculations for the extrapolation of the reach to 50 . The uncertainty of the contribution of electrons from decays was estimated by varying the ratio / by 50% of its original value. The resulting uncertainty is negligible compared to the overall systematics, because the contribution is small, less than 10%.
Over the whole range, the systematic uncertainty due to the subtraction of electrons from charm-hadron decays dominates. The uncertainty due to the subtraction of electrons from light-hadron decays is large at very low , but decreases quickly with increasing as does the overall yield of this background source, as shown in Fig. 2 (left). At high , the uncertainty of the hadron contamination increases.
The influence due to the form factor of electrons from charm and beauty hadron decays as well as light neutral mesons was studied using different decayers and estimated to be negligible.
As the individual sources of systematic uncertainties are uncorrelated, they were added in quadrature to obtain the total systematic uncertainty. The total uncertainty amounts to 38% for the lowest interval and decreases to 12% at = 8 .
The systematic uncertainty due to the determination of the nucleon–nucleon cross section for the minimum-bias trigger condition is 3.7% .
|Tracking and matching||5.6%||5.2%|
|TOF matching and eID||5.4%||n/a|
|Hadron contamination||n/a||1% to 6%|
|Mismatched conversions||4% to 0.3%||negligible|
|Light- and strange-hadron decay bkg.||17% to 1.5%||1.3% to 0%|
|Charm-hadron decay bkg.||32% to 9.6%||8.9% to 6.2%|
|Total||38% to 14%||12%|
5 Data analysis in Pb–Pb collisions
In the Pb–Pb analysis, the yield of electrons from beauty-hadron decays was extracted using the full information contained within the impact parameter distribution of all electron candidates. From the shape of the impact parameter distribution within one interval, it is possible to infer the contributions from the different electron sources (see Section 3). Templates for these distributions were obtained from Monte Carlo simulations including effects such as particle lifetime and the detector response. The templates were then added with appropriate weights to reproduce the measured impact parameter distribution for all electron candidates. Examples are shown in Fig. 1. The template fits were performed using the method proposed in . The approach relies on the accurate description of the impact parameter distributions in Monte Carlo simulations. Thus, detailed studies of differences between the impact parameter distributions in data and Monte Carlo simulations were performed. Differences were corrected for, while the related uncertainties were propagated as detailed below. For the template fit method, four classes of electron sources were distinguished. Their impact parameter distributions, as provided by the Monte Carlo simulations for each interval, will be referred to as fit templates in the following. The four categories correspond to electrons from beauty-hadron decays, from charm-hadron decays, from photon conversions and electrons from other processes, which will be referred to as ‘Dalitz electrons’. The latter is dominated by electrons from Dalitz decays of light neutral mesons. Given that these electrons essentially originate from the interaction point with respect to the detector resolution, the measured impact parameter distribution depends only on the transverse momentum of the electron. Similarly, the remaining hadron contamination consists of particles mostly produced close to the interaction point making its impact parameter distribution similar to that of the Dalitz electrons.
5.1 Extraction of electrons from beauty-hadron decays
The fit templates from the Monte Carlo simulations can be considered as random samples of the unknown true distributions. For each of the four electron sources considered in the previous section, there is a number of counts in the template for each impact parameter bin (see Fig. 1). The number of counts from a particular electron source in a particular bin is called . Its unknown expectation value is called and is considered as a free parameter of the fit. The fit function is the sum of the expectation values, each weighted with the appropriate amplitude parameter : . The bin counts of the impact parameter templates are connected to their expectation values via Poisson statistics. The same relation holds between the fit function and the data () within each impact parameter bin leading to the likelihood distribution 
This gives one free amplitude parameter for each electron source () and one free expectation value parameter for each electron source and impact parameter bin (). The main parameters of interest are the , in particular , while the nuisance parameters arise due to the finite statistics of the templates. Evaluating the full likelihood distribution in several hundred dimensions is challenging. Therefore a simpler approach is to use the maximum likelihood as an estimator for the amplitudes of the electron sources.
An iterative procedure to find the maximum likelihood with respect to the for fixed is suggested in . Numerical minimisation is then performed only for the . Equations for the iterative procedure can be found by setting the differentials to zero. Solving these equations for , yields an iterative rule for each bin.
For a bin with a finite number of entries from data, but zero counts in any of the templates, the likelihood distribution of the is not well represented by its maximum. This happens mostly in the tails of the distributions (see Fig. 1), where the contribution of electrons from beauty-hadron decays dominates. Thus, for this case only the contribution from this source was considered.
To obtain the raw yield of the signal, i.e. electrons from beauty-hadron decays, in a given interval, the number of electrons in the template was scaled by the amplitude parameter . As in the p–Pb analysis, the raw yield was then corrected for the geometrical acceptance, the track reconstruction and selection criteria, the TOF acceptance, and the TOF eID using Monte Carlo simulations. The TPC eID efficiency () was determined via a data-driven approach using electrons from photon conversions identified via topological cuts and the invariant mass . The corresponding distributions were fitted with the function , which describes the distributions including fluctuations, and the efficiency determined as the ratio of electrons before and after the TPC eID selection criterion (see Section 3). Next the spectrum was unfolded. The off-diagonal elements of the response matrix are small. For this reason no regularisation was used in the unfolding procedure to avoid additional systematic uncertainties. The unfolding was done using a matrix inversion of the response matrix . Due to the restricted range of the measurement there is some dependence of the unfolded values on bins that have not been measured, mainly the adjacent bins. To solve this, the yield was measured in two further bins ( and ) and used only in the unfolding calculations. The statistical uncertainties were propagated accordingly.
To validate this signal extraction method, the template fit method was also applied to the p–Pb data, where results were found to be consistent with the cut method described in Section 4.
5.2 Systematic uncertainties estimation
The systematic uncertainties are summarized in Table 2. They were estimated using data-driven methods where possible. An overview of the efficiencies of the different track selection steps may be found in Fig. 3.
The efficiency due to the ITS track selection criteria (hits in both SPD layers) does not depend strongly on the particle species. Thus, charged tracks could be used as a representative sample with respect to the geometric effects, such as inactive areas of the detector. The normalisation for the efficiency was performed by making use of phase space (pseudorapidity and azimuthal angle) regions where the efficiency was close to unity. Averaging over the phase space yields a proxy for the total efficiency which was compared between data and Monte Carlo simulations and yielded a difference of . The uncertainty for non-geometric effects was estimated to be smaller than . The efficiencies of the requirements on charged tracks with good quality, the TOF matching and TOF eID depend more strongly on the particle type. Therefore, only an electron sample could be representative. It was obtained by selecting electrons from photon conversions. Due to the large particle multiplicity in central Pb–Pb collisions (resulting in a sizeable hadron contamination), the comparison was done using weak additional particle identification (), in more peripheral collisions (), and with different ITS track selection criteria (excluding signals in the innermost layer). To account for biases due to these additional criteria, they were varied and the results were checked for consistency. The estimated systematic uncertainties are about for the requirement of charged tracks with good quality and about for the TOF matching and eID. The systematic uncertainty of the TPC eID includes differences in the eID efficiency for electrons from beauty-hadron decays and for electrons from photon conversion (due to the different pseudorapidity distributions) in the sample as well as the uncertainty of the extrapolation towards lower . The uncertainty due to the modelling of the distribution was checked by comparing different model descriptions with the standard one and by comparing with a sample of pions selected with the TRD and TOF. The total uncertainty of for the TPC eID is the quadratic sum of the following contributions: from the extrapolation, from the pseudorapidity dependence, from a possible dependence and from the tail of the distribution.
To estimate the statistical and systematic uncertainty on the extracted signal yield due to the maximum-likelihood fit, a Monte Carlo closure test was used. For this purpose, the templates were slightly smoothed and the result sampled with the statistics present in the measurement. The pseudo-data was created by using the measured contributions as input. The application of the template fit allowed for a comparison of the measured and true value. Repetitions of this process gave an estimation of the statistical and systematic contribution to the uncertainty. The charm yield of the test was varied to avoid underestimating the uncertainty in intervals with downward fluctuations of the measured charm yield. The systematic uncertainty varies between 19% and 6% between the different intervals.
There is an uncertainty in how well the impact parameter distributions of the different electron sources are described by the Monte Carlo simulations. Where possible, any differences were corrected for. The remaining uncertainty was propagated to the measured spectrum of electrons from beauty-hadron decays by changing the fit templates within their uncertainties.
The different resolution of the impact parameter () with the given track and event selection criteria in Monte Carlo simulations and data was corrected for. The size of the correction was estimated by comparing the impact parameter distributions of primary pions, yielding a 10–12% worse resolution in data compared to the Monte Carlo simulations in the range of the measurement. To correct for this effect, a Gaussian distributed random number was added to each impact parameter value such that the resolution in the Monte Carlo simulations matched that of the data. The central values of the yield of electrons from beauty-hadron decays were estimated using a resolution correction of . The yield using a correction of instead, differs by about at = 1.3 with the difference decreasing quickly towards higher . The effect of the correction was found to be negligible for the p–Pb analysis.
Despite the strong eID requirements, there is a significant contamination of the electron sample by hadrons (mostly charged pions). The contribution was estimated using a clean TPC energy loss signal of pions identified with the TRD, which was fitted to the distribution, suggesting a contamination of the electron candidate sample of about even for low transverse momentum. The contamination was not explicitly subtracted. The impact parameter distribution of charged hadrons is similar to that of the Dalitz template. This means that the contribution of the hadron contamination to the impact parameter distribution was absorbed into the Dalitz template by the fit method. To account for slight differences between the distributions, the result was compared with a fit using the hadron impact parameter template instead. A hypothetical template with the same mixture of Dalitz electrons and hadrons as in data would yield a result between these two extreme cases. For 5 , the fit using the hadron template was used for the central points as the contribution from hadrons dominates compared to that of the Dalitz electrons. The difference in the measured yield of electrons from beauty-hadron decays after exchanging the Dalitz template for the hadron template is 7% at = 1.3 decreasing towards higher transverse momentum. The proton contamination is significant only below = 1.3 .
Like for the p–Pb analysis, the influence of the difference in yield of mismatched conversions in data and Monte Carlo simulations had to be considered, especially as it increases with the multiplicity of the event. By making use of the multiplicity dependence, it was possible to create templates that either over- or underestimate this effect. This was cross-checked using charged pions from decays as done in the p–Pb analysis (see Section 4.2). The change of the resulting measured spectra of electrons from beauty-hadron decays was used as an estimate for the systematic uncertainty, which is 14% at = 1.3 and decreases quickly towards higher transverse momentum.
As for the p–Pb analysis, electrons from secondary pion and three-body decays of hadrons carrying a strange (or anti-strange) quark had to be considered, especially as these have broader impact parameter distributions than Dalitz electrons (see Section 3). Due to the different final states, both the template for electrons from photon conversions and the template for Dalitz electrons are affected. These were split into a contribution from the decay of strange particles and the rest. For the fit they were considered as separate templates, but the amplitude parameters were coupled to have a fixed ratio. This was necessary because the contribution from strangeness is very small and could not be constrained by the information from the impact parameter distribution alone. The relative strength of the strangeness content was varied by a factor of two which includes the variation expected from the measured kaon/pion ratio . The resulting difference in the yield of electrons from beauty-hadron decays was used as the estimate for the systematic uncertainty. It is for low , decreasing towards higher transverse momentum.
Electrons at a fixed transverse momentum have mother particles in a range of values. The impact parameter distributions of electrons depend on the momentum distributions of the mother particles. For the charm case this can be disentangled by making use of the measured charm distribution . For the beauty case this means that the result of the measurement depends on the input beauty-hadron spectrum in the Monte Carlo simulation. The effect was estimated by varying the beauty-hadron distribution of the templates and observing the resulting change in the measured electron distribution. The beauty-hadron distribution was obtained according to PYTHIA simulations with a Perugia-0 tune which describes the measured p–Pb data well. Therefore, an effect of the variation of the distribution was studied by introducing a momentum-dependent nuclear modification factor . An based on a theoretical calculation was used for the central points . It has values near unity for low transverse momenta and drops to about from a hadron of to . This was varied to half its effect () in order to estimate the associated uncertainty. For the charm case, the variation was done according to the measurement uncertainties . The difference in the resulting measured yield of electrons from beauty-hadron decays is about , with no visible dependence.
For the template fit, all species of charmed hadrons were combined into one template. The same holds for the beauty case. The baryon fraction of heavy hadrons is currently not known for Pb–Pb collisions and might be different than for pp collisions. Because of the different masses and decay channels, the various heavy-flavour hadron decays produce electrons with different impact parameter distributions. The templates were split into their contributions from only mesons or only baryons, with fixed ratios of the fit amplitudes. To estimate the uncertainty, the baryon fraction was increased by a factor of three for both charm and beauty simultaneously, motivated by the results of thermal model calculations . This led to a change in the measured yield of electrons from beauty-hadron decays of about with no clear momentum dependence. Decreasing the baryon ratio even to 0 has a smaller effect.
|Tracking and matching||4.7%|
|TOF matching and eID||10%|
|Signal extraction||17% to 12%|
|resolution correction||10% to 0.4%|
|Hadron contamination||7% to 1.4%|
|Mismatched conversions||14% to 0.02%|
|Strangeness||1.3% to 0.3%|
|Mother particle distribution||8%|
|Total||26% to 17%|
6 Reference pp cross sections at = 2.76 TeV and = 5.02 TeV
For the calculations of the nuclear modification factors and , corresponding pp reference spectra at = 5.02 TeV and = 2.76 TeV are needed. To obtain these, the same method is used in both analyses. It is described in more detail in the following for the p–Pb analysis.
At present no pp measurement at = 5.02 TeV exists. Therefore, the cross section of electrons from beauty-hadron decays measured in the momentum interval at = 7 TeV  was scaled to = 5.02 TeV by applying a pQCD-driven -scaling . The -dependent scaling function was obtained by calculating the ratio of the production cross sections of electrons from beauty-hadron decays from FONLL pQCD calculations [29, 30, 31] at = 5.02 TeV and = 7 TeV. Both the direct () and the cascade decay () were considered. For the calculations at both energies the same parameters were used for the beauty-quark mass ( ), the PDFs (CTEQ6.6 ) as well as the factorisation and renormalisation scales with = = = , where denotes the transverse momentum of the beauty quark. The uncertainties of the -dependent scaling function were estimated by varying the parameters. The beauty-quark mass was set to 4.5 and 5 . The uncertainties for the PDFs were obtained by using the CTEQ6.6 PDF uncertainties . The contribution from the scale uncertainties was estimated by using six different sets: () = (0.5,0.5),(1,0.5),(0.5,1),(2,1),(1,2),(2,2). The uncertainties originating from the mass and PDF variations are negligible. The uncertainty stemming from the variation of the scales was defined as the largest deviation from the scaling factor obtained with = = . The resulting -scaling uncertainty is almost independent of . It ranges from % at 1 to about % at 8 . The total systematic uncertainty of the pp reference spectrum at = 5.02 TeV is then given as the bin-by-bin quadratic sum of the -scaling uncertainty and the relative systematic uncertainty of the measured spectrum at = 7 TeV. For the statistical uncertainties the relative uncertainties of the spectrum measured at = 7 TeV were taken.
For the analysis, the measured spectrum at = 7 TeV was scaled to = 2.76 TeV using FONLL pQCD calculations at the respective energies. The systematic scaling uncertainty is about % at 1 and about % at 8 . The resulting pp reference spectrum was found to be consistent with the measurement of electrons from beauty-hadron decays in pp collisions at = 2.76 TeV , shown in Fig. 4. The measured spectrum at = 2.76 TeV was not taken as a reference for the , because of larger statistical and systematic uncertainties than the reference obtained via the -scaling.
The systematic uncertainty of the normalisation related to the determination of the cross section of the minimum-bias trigger used for the measurement at = 7 TeV is 3.5% and also holds for the obtained pp reference spectra at = 5.02 TeV and = 2.76 TeV.
The systematic uncertainties of the input -differential cross section of electrons from beauty-hadron decays measured at = 7 TeV, the normalisation uncertainty, as well as the scaling uncertainties for the reference spectra are summarised in Table 3.
|pp spectrum 7 TeV||
|scaling uncertainty for||p–Pb ( = 5.02 TeV)||Pb–Pb ( = 2.76 TeV)|
|at = 1||%||%|
|at = 8||%||%|
The -differential cross section and invariant yield of electrons from beauty-hadron decays at mid-rapidity in minimum-bias p–Pb collisions at = 5.02 TeV and in the 20% most central Pb–Pb collisions at = 2.76 TeV, respectively, are shown in Fig. 5. The markers are plotted at the centre of the bin. The vertical bars indicate the statistical uncertainties, the boxes represent the systematic uncertainties. The pp reference spectra, obtained via the pQCD-driven -scaling from the measurement in pp collisions at = 7 TeV as described in Section 6, are shown for comparison. The pp reference spectra were multiplied by the number of nucleons in the Pb nucleus (A = 208) for the p–Pb and with the nuclear overlap function () for the Pb–Pb comparison. The Pb–Pb result shows a suppression of electrons from beauty-hadron decays at high compared with the yield in pp collisions. Such a suppression is not seen in the comparison of the p–Pb spectrum with the corresponding pp reference.
The nuclear modification factors and are shown in Fig. 6 (left). The was obtained using Equation 1. The was calculated as the ratio of the cross section of electrons from beauty-hadron decays in p–Pb and pp collisions scaled by the number of nucleons in the Pb nucleus (A = 208). The statistical and systematic uncertainties of the Pb–Pb or p–Pb and the pp spectra were propagated as independent uncertainties. The systematic uncertainties of the nuclear modification factors are partially correlated between the bins. The normalisation uncertainty of the pp spectrum and the uncertainty of the nuclear overlap function or the normalisation uncertainties of the p–Pb spectrum, respectively, were added in quadrature. The normalisation uncertainties are shown as filled boxes at high transverse momentum in Fig. 6.
The is consistent with unity within uncertainties (of about 20% for 2 ) for all shown transverse momenta. The production of electrons from beauty-hadron decays is thus consistent with binary-collision scaling of the corresponding measurement in pp collisions at the same centre-of-mass energy. The values of the for the 20% most central Pb–Pb collisions increase, for 3 , with sizeable uncertainties of 30–45%. In the interval , the is about 0.7 with a systematic uncertainty of about 30%; in the ratio is 0.48 with an uncertainty of about 25%. In the latter transverse momentum range the suppression with respect to = 1 is a 3.3 effect taking into account the statistical and systematic uncertainties.
A comparison of the of electrons from beauty-hadron decays with the one from charm- and beauty-hadron decays is shown in Fig. 6 (right) for the 20% most central Pb–Pb collisions. For the latter , the -differential invariant yields of electrons from charm- and beauty-hadron decays published in  for the centrality classes 0–10% and 10–20% were combined. For the pp reference in the momentum range up to 12 , the corresponding invariant cross section measurement at = 2.76 TeV , which has uncertainties of about 20%, was used. For 12 , the ATLAS measurement  at = 7 TeV was extrapolated to = 2.76 TeV applying a FONLL pQCD-driven -scaling analogous to the method described in Section 6. The uncertainty of the pp reference in this momentum range is about 15%. As expected, the results agree within uncertainties at high , where the beauty contribution is larger than the charm contribution . In the interval , the suppression of the for electrons from beauty-hadron decays is about 1.2 less. This difference is consistent with the ordering of charm and beauty suppression seen in the prompt D meson and J/ from B meson comparison [34, 40, 41].
Within uncertainties, the is described by pQCD calculations including modifications of the parton distribution functions (FONLL [29, 30, 31] + EPS09NLO  nuclear PDFs) as shown in Fig. 7 (left). The data and the calculation suggest that cold nuclear matter effects are small at high transverse momentum. Recent measurements of long-range correlations for charged hadrons [51, 53, 54] and studies of the mean transverse momentum as a function of the charged-particle multiplicity in the event  suggest that there might be collective effects in p–Pb collisions. The figure also reports the result of a calculation based on the idea proposed in Ref. , in which the distribution of beauty hadrons from a hydrodynamically expanding medium is obtained from a blast-wave model. The blast-wave parameters were extracted from fits to the -spectra of light hadrons  in p–Pb collisions. The uncertainties of the measurement do not allow for a conclusion on possible flow effects. The data are also described by calculations which include CNM energy loss, nuclear shadowing and coherent multiple scattering at the partonic level . An enhancement at intermediate is predicted by the calculations based on incoherent multiple scattering . Presently, the large systematic uncertainties of the measurement do not allow one to discriminate between the aforementioned theoretical approaches.
Perturbative QCD calculations including initial-state effects for Pb–Pb collisions at = 2.76 TeV (FONLL [29, 30, 31] + EPS09NLO  nuclear PDFs) cannot describe the at high transverse momentum (see Fig. 7, right), indicating that the suppression, particularly evident in the interval , is induced by the presence of a hot and dense medium in the final state. At lower transverse momentum, the large uncertainties do not allow one to conclude whether the measured is larger than that obtained from this calculation.
In order to gain further insight into the energy loss mechanisms, particularly the relative importance of radiative and collisional energy loss, the data are compared with several models of heavy-quark transport and energy loss in the QGP. Both radiative and collisional energy loss are included in the pQCD model MC@sHQ+EPOS2 , the partonic transport description BAMPS [96, 97], and in WHDG [93, 94, 95]. The non-perturbative transport model TAMU  includes only collisional processes, while the POWLANG  transport calculation simulates the production of heavy quarks using POWHEG and their propagation in the plasma via a relativistic Langevin equation. Heavy-quark energy loss can also be calculated using the AdS/CFT heavy-quark drag model .
The right-hand side of Fig. 7 shows the comparison of the various models with the measured . The MC@sHQ+EPOS2 calculation with EPOS initial conditions [98, 99], including the Landau-Pomeranchuk-Migdal (LPM) effect , is consistent with the data at high . The BAMPS [96, 97] model is based on pQCD cross sections including the running of the coupling and scaled by a constant factor . The two shown values of cannot be distinguished given the uncertainties in the data. In the WHDG calculation, the medium density is assumed to be proportional to the charged particle multiplicity and a 1-D Bjorken-expansion is included. The WHDG model describes the measurement well within the restricted range shown.
The TAMU model includes collisional processes and incorporates resonance formation close to the critical temperature as well as diffusion of heavy-flavour mesons in the hadronic phase. The hydrodynamic expansion is constrained by spectra and elliptic flow measurements of light hadrons. The calculations are consistent with the data at high , indicating a limited sensitivity of the current data to radiative energy loss effects. The POWLANG  transport calculation takes into account initial-state nuclear effects via EPS09 modifications of the PDFs and describes the medium using an underlying hydrodynamical model. The transport coefficients used for the evolution of the heavy quark in the medium are either extracted from lattice-QCD calculations or Hard-Thermal-Loop (HTL) resummation  of medium effects. The hadronisation via in-vacuum fragmentation functions or via in-medium string-fragmentation routines occurs once the decoupling temperature is reached. The calculations are shown for different transport coefficients with a decoupling temperature = 155 MeV; the results with a temperature of = 170 MeV look similar. No scenario is clearly favoured by the current data set. The AdS/CFT model, which includes energy loss fluctuations in a realistic strong-coupling energy loss mode, clearly shows a stronger suppression than the measured .
The MC@sHQ+EPOS2, the BAMPS as well as the TAMU calculation describe the suppression seen in data at high transverse momentum. They also show an increase towards lower momentum reaching values around unity or slightly above. The data show a larger increase with decreasing transverse momentum, however exhibit large systematic and statistical uncertainties.
The -differential cross section and invariant yield of electrons from beauty-hadron decays in minimum-bias p–Pb collisions and in the 20% most central Pb–Pb collisions, respectively, were measured at mid-rapidity. The measurements are compared via the nuclear modification factors with pp reference spectra, obtained by a pQCD-driven -scaling of the cross section of electrons from beauty-hadron decays measured at = 7 TeV. The is consistent with unity within uncertainties of about 20% at high transverse momentum , which increase towards low . The is described by pQCD calculations including initial-state effects, energy loss approaches as well as by a blast wave model calculation that parametrises possible hydrodynamic effects. The is about 0.7 with an uncertainty of about 30% in the interval and 0.48 with an uncertainty of about 25% for . The suppression seen in the higher transverse momentum interval is not described by pQCD calculations including only initial-state effects, indicating a final-state effect as the origin. The values of the increase for with uncertainties of about 30–45%. The measured is described within uncertainties by pQCD-inspired models of beauty-quark energy loss in the QGP. In the interval , we observe that the suppression of the for electrons from beauty-hadron decays is about 1.2 less than that from charm- and beauty hadron decays. This difference is consistent with the ordering of charm and beauty suppression seen in the prompt D meson and J/ from B meson comparison.
The ALICE Collaboration would like to thank all its engineers and technicians for their invaluable contributions to the construction of the experiment and the CERN accelerator teams for the outstanding performance of the LHC complex. The ALICE Collaboration gratefully acknowledges the resources and support provided by all Grid centres and the Worldwide LHC Computing Grid (WLCG) collaboration. The ALICE Collaboration acknowledges the following funding agencies for their support in building and running the ALICE detector: A. I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation (ANSL), State Committee of Science and World Federation of Scientists (WFS), Armenia; Austrian Academy of Sciences and Nationalstiftung für Forschung, Technologie und Entwicklung, Austria; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Financiadora de Estudos e Projetos (Finep) and Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Brazil; Ministry of Education of China (MOE of China), Ministry of Science & Technology of China (MOST of China) and National Natural Science Foundation of China (NSFC), China; Ministry of Science, Education and Sport and Croatian Science Foundation, Croatia; Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Cuba; Ministry of Education, Youth and Sports of the Czech Republic, Czech Republic; Danish National Research Foundation (DNRF), The Carlsberg Foundation and The Danish Council for Independent Research | Natural Sciences, Denmark; Helsinki Institute of Physics (HIP), Finland; Commissariat à l’Energie Atomique (CEA) and Institut National de Physique Nucléaire et de Physique des Particules (IN2P3) and Centre National de la Recherche Scientifique (CNRS), France; Bundesministerium für Bildung, Wissenschaft, Forschung und Technologie (BMBF) and GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany; Ministry of Education, Research and Religious Affairs, Greece; National Research, Development and Innovation Office, Hungary; Department of Atomic Energy Government of India (DAE), India; Indonesian Institute of Science, Indonesia; Centro Fermi - Museo Storico della Fisica e Centro Studi e Ricerche Enrico Fermi and Istituto Nazionale di Fisica Nucleare (INFN), Italy; Institute for Innovative Science and Technology , Nagasaki Institute of Applied Science (IIST), Japan Society for the Promotion of Science (JSPS) KAKENHI and Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan; Consejo Nacional de Ciencia (CONACYT) y Tecnología, through Fondo de Cooperación Internacional en Ciencia y Tecnología (FONCICYT) and Dirección General de Asuntos del Personal Academico (DGAPA), Mexico; Nationaal instituut voor subatomaire fysica (Nikhef), Netherlands; The Research Council of Norway, Norway; Commission on Science and Technology for Sustainable Development in the South (COMSATS), Pakistan; Pontificia Universidad Católica del Perú, Peru; Ministry of Science and Higher Education and National Science Centre, Poland; Ministry of Education and Scientific Research, Institute of Atomic Physics and Romanian National Agency for Science, Technology and Innovation, Romania; Joint Institute for Nuclear Research (JINR), Ministry of Education and Science of the Russian Federation and National Research Centre Kurchatov Institute, Russia; Ministry of Education, Science, Research and Sport of the Slovak Republic, Slovakia; National Research Foundation of South Africa, South Africa; Korea Institute of Science and Technology Information and National Research Foundation of Korea (NRF), South Korea; Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT) and Ministerio de Ciencia e Innovacion, Spain; Knut & Alice Wallenberg Foundation (KAW) and Swedish Research Council (VR), Sweden; European Organization for Nuclear Research, Switzerland; National Science and Technology Development Agency (NSDTA), Office of the Higher Education Commission under NRU project of Thailand and Suranaree University of Technology (SUT), Thailand; Turkish Atomic Energy Agency (TAEK), Turkey; National Academy of Sciences of Ukraine, Ukraine; Science and Technology Facilities Council (STFC), United Kingdom; National Science Foundation of the United States of America (NSF) and United States Department of Energy, Office of Nuclear Physics (DOE NP), United States.
-  Wuppertal-Budapest Collaboration, S. Borsanyi, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg, C. Ratti, and K. K. Szabo, “Is there still any mystery in lattice QCD? Results with physical masses in the continuum limit III,” JHEP 09 (2010) 073, arXiv:1005.3508 [hep-lat].
-  HotQCD Collaboration, A. Bazavov et al., “Equation of state in ( 2+1 )-flavor QCD,” Phys. Rev. D90 (2014) 094503, arXiv:1407.6387 [hep-lat].
-  A. Andronic et al., “Heavy-flavour and quarkonium production in the LHC era: from proton-proton to heavy-ion collisions,” Eur. Phys. J. C76 (2016) 107, arXiv:1506.03981 [nucl-ex].
-  F.-M. Liu and S.-X. Liu, “Quark-gluon plasma formation time and direct photons from heavy ion collisions,” Phys. Rev. C89 (2014) 034906, arXiv:1212.6587 [nucl-th].
-  ALICE Collaboration, K. Aamodt et al., “Two-pion Bose-Einstein correlations in central Pb–Pb collisions at 2.76 TeV,” Phys. Lett. B696 (2013) 328–337, arXiv:1012.4035 [nucl-ex].
-  M. Gyulassy and X.-N. Wang, “Multiple collisions and induced gluon bremsstrahlung in QCD,” Nucl. Phys. B420 (1994) 583–614, arXiv:nucl-th/9306003 [nucl-th].
-  R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne, and D. Schiff, “Radiative energy loss and p(T) broadening of high-energy partons in nuclei,” Nucl. Phys. B484 (1997) 265–282, arXiv:hep-ph/9608322 [hep-ph].
-  U. A. Wiedemann, “Gluon radiation off hard quarks in a nuclear environment: opacity expansion,” Nucl. Phys. B588 (2000) 303–344, arXiv:hep-ph/0005129 [hep-ph].
-  Y. L. Dokshitzer and D. E. Kharzeev, “Heavy quark colorimetry of QCD matter,” Phys. Lett. B519 (2001) 199–206, arXiv:hep-ph/0106202 [hep-ph].
-  N. Armesto, C. A. Salgado, and U. A. Wiedemann, “Medium induced gluon radiation off massive quarks fills the dead cone,” Phys. Rev. D69 (2004) 114003, arXiv:hep-ph/0312106 [hep-ph].
-  B.-W. Zhang, E. Wang, and X.-N. Wang, “Heavy Quark Energy Loss in Nuclear Medium,” Phys. Rev. Lett. 93 (2004) 072301, arXiv:nucl-th/0309040 [nucl-th].
-  M. Djordjevic and M. Gyulassy, “Heavy quark radiative energy loss in QCD matter,” Nucl. Phys. A733 (2004) 265–298, arXiv:nucl-th/0310076 [nucl-th].
-  N. Armesto, A. Dainese, C. A. Salgado, and U. A. Wiedemann, “Testing the color charge and mass dependence of parton energy loss with heavy-to-light ratios at RHIC and CERN LHC,” Phys. Rev. D71 (2005) 054027, arXiv:hep-ph/0501225 [hep-ph].
-  N. Armesto, “Nuclear shadowing,” J. Phys. G32 (2006) R367–R394, arXiv:hep-ph/0604108 [hep-ph].
-  H. Fujii, F. Gelis, and R. Venugopalan, “Quark pair production in high energy pA collisions: General features,” Nucl. Phys. A780 (2006) 146–174, arXiv:hep-ph/0603099 [hep-ph].
-  J. W. Cronin, H. J. Frisch, M. J. Shochet, J. P. Boymond, R. Mermod, P. A. Piroue, and R. L. Sumner, “Production of hadrons with large transverse momentum at 200, 300, and 400 GeV,” Phys. Rev. D11 (1975) 3105.
-  R. Sharma, I. Vitev, and B.-W. Zhang, “Light-cone wave function approach to open heavy flavor dynamics in QCD matter,” Phys. Rev. C80 (2009) 054902, arXiv:0904.0032 [hep-ph].
-  Z.-B. Kang, I. Vitev, E. Wang, H. Xing, and C. Zhang, “Multiple scattering effects on heavy meson production in p+A collisions at backward rapidity,” Phys. Lett. B740 (2015) 23–29, arXiv:1409.2494 [hep-ph].
-  M. L. Miller, K. Reygers, S. J. Sanders, and P. Steinberg, “Glauber Modeling in High Energy Nuclear Collisions,” Ann. Rev. Nucl. Part. Sci. 57 (2007) 205–243, arXiv:nucl-ex/0701025 [nucl-ex].
-  STAR Collaboration, M. M. Aggarwal et al., “Measurement of the Bottom Quark Contribution to Nonphotonic Electron Production in Collisions at =200 GeV,” Phys. Rev. Lett. 105 (2010) 202301, arXiv:1007.1200 [nucl-ex].
-  STAR Collaboration, H. Agakishiev et al., “High non-photonic electron production in collisions at GeV,” Phys. Rev. D83 (2011) 052006, arXiv:1102.2611 [nucl-ex].
-  CDF Collaboration, D. Acosta et al., “Measurement of the meson and -hadron production cross sections in collisions at TeV,” Phys. Rev. D71 (2005) 032001, arXiv:0412071 [hep-ex].
-  ALICE Collaboration, B. Abelev et al., “Measurement of electrons from beauty hadron decays in collisions at TeV,” Phys. Lett. B721 (2013) 13–23, arXiv:1208.1902 [hep-ex]. [Erratum: Phys. Lett.B763,507(2016)].
-  ALICE Collaboration, B. Abelev et al., “Measurement of electrons from semileptonic heavy-flavor hadron decays in collisions at TeV,” Phys. Rev. D91 (2015) 012001, arXiv:1405.4117 [nucl-ex].
-  ALICE Collaboration, B. Abelev et al., “Measurement of prompt and beauty hadron production cross sections at mid-rapidity in collisions at TeV,” JHEP 11 (2012) 065, arXiv:1205.5880 [hep-ex].
-  LHCb Collaboration, R. Aaij et al., “Measurement of the production cross-section in collisions at TeV,” JHEP 04 (2012) 093, arXiv:1202.4812 [hep-ex].
-  CMS Collaboration, S. Chatrchyan et al., “Measurement of the Strange Meson Production Cross Section with J/Psi Decays in Collisions at TeV,” Phys. Rev. D84 (2011) 052008, arXiv:1106.4048 [hep-ex]. and references therein.
-  ATLAS Collaboration, G. Aad et al., “Measurement of the differential cross-sections of inclusive, prompt and non-prompt production in proton-proton collisions at TeV,” Nucl. Phys. B850 (2011) 387–444, arXiv:1104.3038 [hep-ex].
-  M. Cacciari, M. Greco, and P. Nason, “The spectrum in heavy flavor hadroproduction,” JHEP 05 (1998) 007, arXiv:hep-ph/9803400 [hep-ph].
-  M. Cacciari, S. Frixione, and P. Nason, “The spectrum in heavy flavor photoproduction,” JHEP 03 (2001) 006, arXiv:hep-ph/0102134 [hep-ph].
-  M. Cacciari, S. Frixione, N. Houdeau, M. L. Mangano, P. Nason, and G. Ridolfi, “Theoretical predictions for charm and bottom production at the LHC,” JHEP 10 (2012) 137, arXiv:1205.6344 [hep-ph].
-  PHENIX Collaboration, A. Adare et al., “Energy Loss and Flow of Heavy Quarks in Au+Au Collisions at GeV,” Phys. Rev. Lett. 98 (2007) 172301, arXiv:nucl-ex/0611018 [nucl-ex].
-  PHENIX Collaboration, A. Adare et al., “Heavy quark production in and energy loss and flow of heavy quarks in Au+Au collisions at GeV,” Phys. Rev. C84 (2011) 044905, arXiv:1005.1627 [nucl-ex].
-  ALICE Collaboration, J. Adam et al., “Centrality dependence of high- D meson suppression in Pb–Pb collisions at = 2.76 TeV,” JHEP 11 (2015) 205, arXiv:1506.06604 [nucl-ex].
-  ALICE Collaboration, J. Adam et al., “Transverse momentum dependence of D-meson production in Pb–Pb collisions at 2.76 TeV,” JHEP 03 (2016) 081, arXiv:1509.06888 [nucl-ex].
-  ALICE Collaboration, B. Abelev et al., “Production of Muons from Heavy Flavor Decays at Forward Rapidity in pp and Pb–Pb Collisions at = 2.76 TeV,” Phys. Rev. Lett. 109 (2012) 112301, arXiv:1205.6443 [hep-ex].
-  ALICE Collaboration, B. Abelev et al., “Production of charged pions, kaons and protons at large transverse momenta in pp and Pb–Pb collisions at =2.76 TeV,” Phys. Lett. B736 (2014) 196–207, arXiv:1401.1250 [nucl-ex].
-  M. Djordjevic, M. Djordjevic, and B. Blagojevic, “RHIC and LHC jet suppression in non-central collisions,” Phys. Lett. B737 (2014) 298–302, arXiv:1405.4250 [nucl-th].
-  M. Djordjevic, “Heavy flavor puzzle at LHC: a serendipitous interplay of jet suppression and fragmentation,” Phys. Rev. Lett. 112 (2014) 042302, arXiv:1307.4702 [nucl-th].
-  ALICE Collaboration, J. Adam et al., “Inclusive, prompt and non-prompt J/ production at mid-rapidity in Pb–Pb collisions at = 2.76 TeV,” JHEP 07 (2015) 051, arXiv:1504.07151 [nucl-ex].
-  CMS Collaboration, S. Chatrchyan et al., “Suppression of non-prompt , prompt , and Y(1S) in PbPb collisions at TeV,” JHEP 05 (2012) 063, arXiv:1201.5069 [nucl-ex].
-  CMS Collaboration, S. Chatrchyan et al., “Evidence of b-Jet Quenching in PbPb Collisions at TeV,” Phys. Rev. Lett. 113 (2014) 132301, arXiv:1312.4198 [nucl-ex]. [Erratum: Phys. Rev. Lett.115,029903(2015)].
-  PHENIX Collaboration, A. Adare et al., “Single electron yields from semileptonic charm and bottom hadron decays in AuAu collisions at GeV,” Phys. Rev. C93 (2016) 034904, arXiv:1509.04662 [nucl-ex].
-  CMS Collaboration, V. Khachatryan et al., “Study of Meson Production in Collisions at Using Exclusive Hadronic Decays,” Phys. Rev. Lett. 116 (2016) 032301, arXiv:1508.06678 [nucl-ex].
-  CMS Collaboration, V. Khachatryan et al., “Transverse momentum spectra of inclusive b jets in pPb collisions at 5.02 TeV,” Phys. Lett. B754 (2016) 59, arXiv:1510.03373 [nucl-ex].
-  CMS Collaboration, “ production in p Pb collisions,”. CMS-PAS-HIN-14-009.
-  LHCb Collaboration, R. Aaij et al., “Study of production and cold nuclear matter effects in collisions at TeV,” JHEP 02 (2014) 072, arXiv:1308.6729 [nucl-ex].
-  ALICE Collaboration, B. Abelev et al., “Measurement of Prompt -Meson Production in Collisions at = 5.02 TeV,” Phys. Rev. Lett. 113 (2014) 232301, arXiv:1405.3452 [nucl-ex].
-  ALICE Collaboration, J. Adam et al., “Measurement of electrons from heavy-flavour hadron decays in p–Pb collisions at TeV,” Phys. Lett. B754 (2016) 81–93, arXiv:1509.07491 [nucl-ex].
-  PHENIX Collaboration, A. Adare et al., “Cold-Nuclear-Matter Effects on Heavy-Quark Production at Forward and Backward Rapidity in d+Au Collisions at =200 GeV,” Phys. Rev. Lett. 112 (2014) 252301, arXiv:1310.1005 [nucl-ex].
-  ALICE Collaboration, B. Abelev et al., “Long-range angular correlations on the near and away side in -Pb collisions at TeV,” Phys. Lett. B719 (2013) 29–41, arXiv:1212.2001 [nucl-ex].
-  ALICE Collaboration, B. Abelev et al., “Long-range angular correlations of , K and p in p–Pb collisions at = 5.02 TeV,” Phys. Lett. B726 (2013) 164–177, arXiv:1307.3237 [nucl-ex].
-  ATLAS Collaboration, G. Aad et al., “Measurement with the ATLAS detector of multi-particle azimuthal correlations in p+Pb collisions at =5.02 TeV,” Phys. Lett. B725 (2013) 60–78, arXiv:1303.2084 [hep-ex].
-  CMS Collaboration, S. Chatrchyan et al., “Observation of long-range near-side angular correlations in proton-lead collisions at the LHC,” Phys. Lett. B718 (2013) 795–814, arXiv:1210.5482 [nucl-ex].
-  PHENIX Collaboration, A. Adare et al., “Quadrupole Anisotropy in Dihadron Azimuthal Correlations in Central Au Collisions at =200 GeV,” Phys. Rev. Lett. 111 (2013) 212301, arXiv:1303.1794 [nucl-ex].
-  STAR Collaboration, L. Adamczyk et al., “Long-range pseudorapidity dihadron correlations in +Au collisions at GeV,” Phys. Lett. B747 (2015) 265–271, arXiv:1502.07652 [nucl-ex].
-  A. M. Sickles, “Possible evidence for radial flow of heavy mesons in d+Au collisions,” Phys. Lett. B731 (2014) 51–56, arXiv:1309.6924 [nucl-th].
-  A. Dainese, Charm production and in-medium QCD energy loss in nucleus nucleus collisions with ALICE: A Performance study. PhD thesis, Padua U., 2003. arXiv:nucl-ex/0311004 [nucl-ex].
-  ALICE Collaboration, K. Aamodt et al., “The ALICE experiment at the CERN LHC,” JINST 3 (2008) S08002.
-  ALICE Collaboration, B. Abelev et al., “Performance of the ALICE experiment at the CERN LHC,” Int. J. Mod. Phys. A29 (2014) 1430044, arXiv:1402.4476 [nucl-ex].
-  ALICE Collaboration, K. Aamodt et al., “Alignment of the ALICE Inner Tracking System with cosmic-ray tracks,” JINST 5 (2010) P03003, arXiv:1001.0502 [physics.ins-det].
-  J. Alme et al., “The ALICE TPC, a large 3-dimensional tracking device with fast readout for ultra-high multiplicity events,” Nucl. Instrum. Meth. A622 (2010) 316–367, arXiv:1001.1950 [physics.ins-det].
-  ALICE Collaboration, P. Cortese et al., “ALICE Transition-Radiation Detector : Technical Design Report,”. https://cds.cern.ch/record/519145. CERN-LHCC-2001-021.
-  ALICE Collaboration, P. Cortese et al., “ALICE Time-Of Flight system (TOF) : addendum to the Technical Design Report,”. https://cds.cern.ch/record/545834. CERN-LHCC-2002-016.
-  ALICE Collaboration, P. Cortese et al., “ALICE Forward Detectors: FMD, TO and VO : Technical Design Report,”. https://cds.cern.ch/record/781854. CERN-LHCC-2004-025.
-  ALICE Collaboration, B. Abelev et al., “Centrality determination of Pb–Pb collisions at = 2.76 TeV with ALICE,” Phys. Rev. C88 (2013) 044909, arXiv:1301.4361 [nucl-ex].
-  Particle Data Group Collaboration, K. A. Olive et al., “Review of Particle Physics,” Chin. Phys. C38 (2014) 090001.
-  M. Gyulassy and X.-N. Wang, “HIJING 1.0: A Monte Carlo program for parton and particle production in high-energy hadronic and nuclear collisions,” Comput. Phys. Commun. 83 (1994) 307, arXiv:nucl-th/9502021 [nucl-th].
-  T. Sjostrand, S. Mrenna, and P. Z. Skands, “PYTHIA 6.4 physics and manual,” JHEP 05 (2006) 026, arXiv:hep-ph/0603175 [hep-ph].
-  P. Z. Skands, “Tuning Monte Carlo generators: The Perugia tunes,” Phys. Rev. D82 (2010) 074018, arXiv:1005.3457 [hep-ph].
-  R. Brun et al., “GEANT: Detector Description and Simulation Tool,”. https://cds.cern.ch/record/1082634. CERN-W-5013.
-  ALICE Collaboration, B. Abelev et al., “Measurement of electrons from semileptonic heavy-flavour hadron decays in pp collisions at = 7 TeV,” Phys. Rev. D86 (2012) 112007, arXiv:1205.5423 [hep-ex].
-  ALICE Collaboration, B. Abelev et al., “Multiplicity dependence of pion, kaon, proton and lambda production in p–Pb collisions at = 5.02 TeV,” Phys. Lett. B728 (2014) 25–38, arXiv:1307.6796 [nucl-ex].
-  ALICE Collaboration, J. Adam et al., “Multiplicity dependence of charged pion, kaon, and (anti)proton production at large transverse momentum in p–Pb collisions at = 5.02 TeV,” arXiv:1601.03658 [nucl-ex].
-  ZEUS Collaboration, H. Abramowicz et al., “Measurement of charm fragmentation fractions in photoproduction at HERA,” JHEP 09 (2013) 058, arXiv:1306.4862 [hep-ex].
-  G. Cowan, Statistical data analysis. Oxford science publications. Clarendon Press, Oxford, Repr. ed., 2004.
-  ALICE Collaboration, B. Abelev et al., “Pseudorapidity Density of Charged Particles in p–Pb Collisions at =5.02 TeV,” Phys. Rev. Lett. 110 (2013) 032301, arXiv:1210.3615 [nucl-ex].
-  ALICE Collaboration, B. Abelev et al., “Measurement of visible cross sections in proton-lead collisions at = 5.02 TeV in van der Meer scans with the ALICE detector,” JINST 9 (2014) P11003, arXiv:1405.1849 [nucl-ex].
-  ALICE Collaboration, B. Abelev et al., “Transverse momentum dependence of inclusive primary charged-particle production in p–Pb collisions at ,” Eur. Phys. J. C74 (2014) 3054, arXiv:1405.2737 [nucl-ex].
-  ALICE Collaboration, B. Abelev et al., “ and Production in Pb–Pb Collisions at = 2.76 TeV,” Phys. Rev. Lett. 111 (2013) 222301, arXiv:1307.5530 [nucl-ex].
-  R. Barlow and C. Beeston, “Fitting using finite Monte Carlo samples,” Comp. Phys. Comm. 77 no. 2, (1993) 219.
-  ALICE Collaboration, B. Abelev et al., “Neutral pion production at midrapidity in pp and Pb–Pb collisions at ,” Eur. Phys. J. C74 (2014) 3108, arXiv:1405.3794 [nucl-ex].
-  ALICE Collaboration, B. Abelev et al., “Suppression of high transverse momentum D mesons in central Pb–Pb collisions at TeV,” JHEP 09 (2012) 112, arXiv:1203.2160 [nucl-ex].
-  M. He, R. J. Fries, and R. Rapp, “Heavy flavor at the large hadron collider in a strong coupling approach,” Phys. Lett. B735 (2014) 445–450, arXiv:1401.3817 [nucl-th].
-  Y. Oh, C. M. Ko, S. H. Lee, and S. Yasui, “Heavy baryon/meson ratios in relativistic heavy ion collisions,” Phys. Rev. C79 (2009) 044905, arXiv:0901.1382 [nucl-th].
-  R. Averbeck, N. Bastid, Z. C. del Valle, P. Crochet, A. Dainese, and X. Zhang, “Reference heavy flavour cross sections in pp collisions at = 2.76 TeV, using a pQCD-driven -scaling of ALICE measurements at = 7 TeV,” arXiv:1107.3243 [hep-ph].
-  P. M. Nadolsky, H.-L. Lai, Q.-H. Cao, J. Huston, J. Pumplin, D. Stump, W.-K. Tung, and C. P. Yuan, “Implications of CTEQ global analysis for collider observables,” Phys. Rev. D78 (2008) 013004, arXiv:0802.0007 [hep-ph].
-  ALICE Collaboration, B. Abelev et al., “Beauty production in pp collisions at = 2.76 TeV measured via semi-electronic decays,” Phys. Lett. B738 (2014) 97–108, arXiv:1405.4144 [nucl-ex]. [Erratum: Phys. Lett.B763,507(2016)].
-  ALICE Collaboration, J. Adam et al., “Measurement of the production of high- electrons from heavy-flavour hadron decays in Pb-Pb collisions at = 2.76 TeV,” arXiv:1609.07104 [nucl-ex].
-  K. J. Eskola, H. Paukkunen, and C. A. Salgado, “EPS09 – A new generation of NLO and LO nuclear parton distribution functions,” JHEP 04 (2009) 065, arXiv:0902.4154 [hep-ph].
-  M. Nahrgang, J. Aichelin, P. B. Gossiaux, and K. Werner, “Influence of hadronic bound states above on heavy-quark observables in Pb + Pb collisions at at the CERN Large Hadron Collider,” Phys. Rev. C89 (2014) 014905, arXiv:1305.6544 [hep-ph].
-  A. Beraudo, A. De Pace, M. Monteno, M. Nardi, and F. Prino, “Heavy flavors in heavy-ion collisions: quenching, flow and correlations,” Eur. Phys. J. C75 (2015) 121, arXiv:1410.6082 [hep-ph].
-  W. A. Horowitz and M. Gyulassy, “The surprisingly transparent sQGP at LHC,” Nucl. Phys. A872 (2011) 265–285, arXiv:1104.4958 [hep-ph].
-  S. Wicks, W. Horowitz, M. Djordjevic, and M. Gyulassy, “Elastic, inelastic, and path length fluctuations in jet tomography,” Nucl. Phys. A784 (2007) 426–442, arXiv:nucl-th/0512076 [nucl-th].
-  W. A. Horowitz, “Fluctuating heavy quark energy loss in a strongly coupled quark-gluon plasma,” Phys. Rev. D91 (2015) 085019, arXiv:1501.04693 [hep-ph].
-  J. Uphoff, O. Fochler, Z. Xu, and C. Greiner, “Elastic and radiative heavy quark interactions in ultra-relativistic heavy-ion collisions,” J. Phys. G42 (2015) 115106, arXiv:1408.2964 [hep-ph].
-  J. Uphoff, F. Senzel, O. Fochler, C. Wesp, Z. Xu, and C. Greiner, “Elliptic Flow and Nuclear Modification Factor in Ultrarelativistic Heavy-Ion Collisions within a Partonic Transport Model,” Phys. Rev. Lett. 114 (2015) 112301, arXiv:1401.1364 [hep-ph].
-  K. Werner, I. Karpenko, T. Pierog, M. Bleicher, and K. Mikhailov, “Event-by-event simulation of the three-dimensional hydrodynamic evolution from flux tube initial conditions in ultrarelativistic heavy ion collisions,” Phys. Rev. C82 (2010) 044904, arXiv:1004.0805 [nucl-th].
-  K. Werner, I. Karpenko, M. Bleicher, T. Pierog, and S. Porteboeuf-Houssais, “Jets, bulk matter, and their interaction in heavy ion collisions at several TeV,” Phys. Rev. C85 (2012) 064907, arXiv:1203.5704 [nucl-th].
-  R. Baier, Y. L. Dokshitzer, S. Peigne, and D. Schiff, “Induced gluon radiation in a QCD medium,” Phys. Lett. B345 (1995) 277–286, arXiv:hep-ph/9411409 [hep-ph].
-  E. Braaten and R. D. Pisarski, “Soft amplitudes in hot gauge theories: A general analysis,” Nucl. Phys. B337 (1990) 569–634.
Appendix A The ALICE Collaboration
G. Aglieri Rinella,
R. Alfaro Molina,
C. Alves Garcia Prado,
A. Batista Camejo,
H. Bello Martinez,
E. Calvo Villar,
J. Castillo Castellanos,
C. Ceballos Sanchez,
V. Chibante Barroso,
G. Conesa Balbastre,
Z. Conesa del Valle,
Y. Corrales Morales,
I. Cortés Maldonado,
R. Cruz Albino,
A. De Caro,
G. de Cataldo,
C. de Conti,
J. de Cuveland,
A. De Falco,
D. De Gruttola,
N. De Marco,
S. De Pasquale,
R.D. De Souza,
D. Di Bari,
A. Di Mauro,
P. Di Nezza,
B. Di Ruzza,
M.A. Diaz Corchero,
D. Domenicis Gimenez,