University of Birmingham Measurement of electrons from beauty-hadron decays in pPb collisions at sNN = 5 . 02 TeV and Pb-Pb collisions at sNN = 2 . 76 TeV

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 1 < pT < 8 GeV/c in minimum-bias p-Pb collisions at √ sNN = 5.02 TeV and in 1.3 < pT < 8 GeV/c in the 20% most central Pb-Pb collisions at √ sNN = 2.76 TeV. The pp reference spectra at √ s = 5.02 TeV and √ s = 2.76 TeV, needed for the calculation of the nuclear modification factors RpPb and RPbPb, were obtained by a pQCD-driven scaling of the cross section of electrons from beauty-hadron decays measured at √ s = 7 TeV. In the pT interval 3 < pT < 8 GeV/c, a suppression of the yield of electrons from beauty-hadron decays is observed in Pb-Pb compared to pp collisions. Towards lower pT, the RPbPb values increase with large systematic uncertainties. The RpPb 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 RpPb 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 RPbPb below unity at high pT may be ascribed to an effect of the hot and dense medium formed in Pb-Pb collisions.


Introduction
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 (m Q Λ QCD ), 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/c for charm and beauty quarks, respectively [3]. 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/c [4] and its decoupling time is about 10 fm/c for collisions at LHC energies [5]. 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 massdependent 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, p T 10 GeV/c, where the ratio m c /p T approaches zero [13]. Therefore, due to the larger mass, beauty quarks can be sensitive -1 -

JHEP07(2017)052
probes for testing the mass dependence of the parton energy loss up to transverse momenta well above 10 GeV/c [13]. 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, k T 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, R AA (R pA ), defined as the ratio of the p T distributions measured in AA (p-A) collisions with respect to the one in pp collisions: where dN AA /dp T and dσ pp /dp T are the p T -differential yield and production cross section of a given particle species in AA and pp collisions, respectively, and T AA is the average of the nuclear overlap function for the centrality range under study [19]. Previous beauty-hadron production measurements in pp collisions at various energies at RHIC [20,21], the Tevatron [22] 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-p T 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 p T [32][33][34][35][36]. The D meson and pion R PbPb 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 p T spectrum and the fragmentation into hadrons [13,39]. Furthermore, the nuclear modification factors R PbPb of prompt D mesons and of J/ψ from B meson decays were compared in the p T interval 8 < p T < 16 GeV/c for D mesons and 6.5 < p T < 30 GeV/c for J/ψ mesons in order to have a similar average p T (≈10 GeV/c) for the heavy hadrons [34,40,41]. This comparison with models indicates that charm quarks lose more energy than beauty quarks in this p T 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-p T interval 70 < p T < 250 GeV/c [42]. 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 [43]. There is a hint that in the momentum interval 3 < p T < 4 GeV/c the R AuAu of electrons from beauty-hadron decays is larger than that of electrons from charm-hadron decays.

JHEP07(2017)052
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][46][47][48]. Therefore, the observed suppression of charm and beauty yields at high p T 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 √ s NN = 200 GeV at RHIC, an enhancement was measured at backward rapidity by means of R dAu of muons from heavy-flavour hadron decays [50]. 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 v 2 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 R dAu [57]. In this paper, the invariant cross section in p-Pb and yield in Pb-Pb collisions are presented together with the nuclear modification factors, R pPb and R PbPb , of electrons from beauty-hadron decays in p-Pb and Pb-Pb collisions at √ s NN = 5.02 TeV and √ s NN = 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 |y lab | ≤ 0.6 (|y lab | ≤ 0.8) and the p T interval 1.0 < p T < 8.0 GeV/c (1.3 < p T < 8.0 GeV/c). 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 ∆y = 0.465 in the proton beam direction, resulting in the rapidity coverage −1.06 < y cms < 0.14 for electrons. Given the cms energies and the rapidity coverages in the p-Pb and Pb-Pb collisions, both measurements probe, at the lowest p T , similar values of Bjorken-x of about 10 −3 for electrons from beauty-hadron decays [58]. 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 R pPb and R PbPb is reported in section 6. The results are presented and discussed in section 7. Section 8 summarises the results.

Experimental apparatus and data samples
A comprehensive description of the ALICE apparatus and its performance can be found in [59] and [60], 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.

JHEP07(2017)052
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 [61]. In this analysis, the ITS was also used to reconstruct the primary (interaction) vertex and the track impact parameter d 0 , 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 d 0 is better than 65 µm and 70 µm for charged particles with momenta larger than 1 GeV/c in Pb-Pb and p-Pb collisions [60], 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 dE/dx 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 dE/dx 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) [63] surrounding the TPC provides hadron and electron identification via the measurement of the specific energy loss dE/dx 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) [64], based on Multi-gap Resistive Plate Chambers (MRPCs), provides hadron rejection at low transverse momentum via the timeof-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 [65] 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 -4 -JHEP07(2017)052 extent [66]. 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 [19], as described in [66]. The Glauber-model simulation yields a value of 18.93 ± 0.74 mb −1 for the average nuclear overlap function T AA 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 L pPb int = 47.8 ± 1.6 µb −1 and L PbPb int = 2.2 ± 0.2 µb −1 , respectively.

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 (b → e, branching ratio: ≈11%) as well as cascade decays (b → c → e, branching ratio: ≈10%) of hadrons that contain a beauty (or anti-beauty) quark [67]. 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 |y lab | ≤ 0.8 and in p-Pb to |y lab | ≤ 0.6. 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 (dE/dx) and the time-offlight. The selection variable (hereafter n TPC σ or n TOF σ ) 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 (σ TPC or σ TOF ). The expected signal and the resolution originate from parametrisations of the TPC and TOF detector signals, described in detail in [60]. For both analyses, particles were accepted with the TPC as electron candidates if they satisfied the condition −0.5 < n TPC σ < 3. This asymmetric selection was chosen to remove hadrons, that are mainly found at negative n TPC σ 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 |n TOF σ | ≤ 3 was applied for the whole p T range in the Pb-Pb analysis and for p T -5 -JHEP07(2017)052 ≤ 2.5 GeV/c 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.
• (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 (cτ ≈ 500 µm [67]) 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 d 0 < 0. 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 p T 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 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 √ s NN = 2.76 TeV was generated with HIJING v1.36 [68] 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 [69] with Perugia-0 tune [70]. Each added pp event contains one cc or bb pair. For the p-Pb analysis, the same procedure was used. The generated particles were propagated through the ALICE apparatus using GEANT3 [71] 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 π 0 from strange-hadron decays (K 0 S , K 0 L , K ± , Λ) and three prong decays of strange hadrons (K 0 L ,K ± ), 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 -7 -

JHEP07(2017)052
'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 |d 0 |, 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.

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 √ s = 2.76 TeV and √ s = 7 TeV [23, 24]. 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 p T distributions of the background sources in data, respectively, and then subtracted.

Extraction of electrons from beauty-hadron decays
Electron candidates with an impact parameter |d 0 | > 0.0054 + 0.078 × exp(−0.56 × p T ) (with d 0 in cm and p T in GeV/c) were selected. This selection criterion was determined from Monte Carlo simulations to maximise the significance for electrons from beautyhadron decays. The selection of the minimum impact parameter is p T 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 [23] are p T dependent.
The number of hadrons passing the track selection, eID, and the minimum impact parameter requirement was estimated at high transverse momentum (p T ≥ 4 GeV/c) by parametrising the TPC n TPC σ distribution in momentum slices, and it was subtracted [72]. Above a p T of 4 GeV/c, the hadron contamination increases with transverse momentum and reaches 10% at 8 GeV/c, see figure 2 (left). At low transverse momentum (p T ≤ 4 GeV/c), the hadron contamination is negligible except in the transverse momentum interval 1 < p T < 1.2 GeV/c, see figure 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 -8 -

JHEP07(2017)052
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 p T -bins by selecting electrons and protons with tight selection criteria in TOF and TPC. The estimated proton contribution, which is ≈10% (4%) in the p T interval 1 < p T < 1.1 GeV/c (1.1 < p T < 1.2 GeV/c), was subtracted statistically from the measured electron candidate p T 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 p T , 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 p T 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 π 0 production cross section in p-Pb collisions available, the π 0 input was based on the measured charged-pion spectra [73,74] assuming N π 0 = (N π + + N π − )/2. Due to the requirement of a minimum impact parameter, the contribution of electrons from decays of secondary π 0 from strange-hadron decays is comparable with the one from primary decays. Therefore the measured p T spectra of K ± , K 0 S and Λ [73] 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 [72]. The contribution of electrons originating from secondary pions from strange-hadron decays or three-body decays of strange hadrons is shown in figure 2 (left). The other light mesons (η, ρ, ω, η and φ), which contribute little, via Dalitz decays and photon conversions compared to primary π 0 decays, were re-weighted via m T -scaling of the π 0 spectrum [72]. The electron background from charmhadron decays was estimated based on the D 0 , D + and D + s meson production cross section measurements with ALICE [48] in the transverse momentum intervals 1 < p T < 16 GeV/c, 2 < p T < 24 GeV/c and 2 < p T < 12 GeV/c, respectively. In a first step the measurements were extrapolated to the p T interval 1 < p T < 24 GeV/c by assuming constant ratios D 0 /D + and D + s /D 0 from the measured D meson production cross sections. Next the p T differential production cross sections were extrapolated to p T = 50 GeV/c via FONLL pQCD calculations. About 10% of the electrons with p T ≤ 8 GeV/c originate from the extrapolated D meson high-p T region (p T ≥ 24 GeV/c). The electron contribution from Λ + c decays was estimated using the ratio σ (Λ + c )/σ D 0 + D + measured by the ZEUS Collaboration [75]. Analogous to the light mesons, the measured D meson p T spectra were also used to re-weight the p T 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 p T spectrum is shown in figure 2 (left). At p T = 1 GeV/c, 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 1/3. With increasing p T 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 p T interval 4.5 < p T < 5 GeV/c, the S/B ratio reaches its maximum of 3. Here the electron background mostly originates from charm-hadron decays. At higher p T , 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 p T range [72].
The raw yield of electrons from beauty-hadron decays N e,raw was then corrected for the geometrical acceptance and for the efficiency ( rec ) 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 ( TPCeID ) was determined to be 69% via a data-driven approach based on the n TPC σ distributions [72]. The transverse-momentum dependence of the efficiencies is shown in figure 2 (right). The total efficiency shows a significant p T dependence, mainly due to the d 0 cut. The effects of the finite momentum resolution and the energy loss due to Bremsstrahlung were taken into account in a bin-by-bin p T resolution correction step based on a Monte Carlo simulation [23, 76].

JHEP07(2017)052
The p T -differential invariant cross section of electrons from beauty-hadron decays, (e + + e − )/2, is thus given as: where p centre T is the centre of the p T bin with width ∆p T and ∆y denotes the geometrical acceptance in |y lab | to which the analysis was restricted. N mb 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 [77], is σ V0 mb = 2.09 ± 0.07 b [78].

Systematic uncertainties estimation
An overview of the relative systematic uncertainties is shown in table 1. The systematic uncertainties were estimated as a function of p T 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 [79], 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 p T 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 p T .
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 [23]. At 1 GeV/c (8 GeV/c) this corresponds to a ≈10% (≈ 25%) variation of the cut value.
The number of electrons from photon conversions increases quickly with decreasing transverse momentum (see figure 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 K 0 S decays identified via topological and invariant mass cuts [80] 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 K 0 S 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. Re--11 -

JHEP07(2017)052
peating 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 p T distributions in Monte Carlo simulations. Uncertainties due to the m Tscaling of the background yields, estimated as 30% [72], and the extrapolation of the D meson p T 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 p T reach to 50 GeV/c. The uncertainty of the contribution of electrons from Λ c decays was estimated by varying the ratio σ (Λ c )/σ D 0 + D + by ± 50% of its original value. The resulting uncertainty is negligible compared to the overall systematics, because the Λ c contribution is small, less than 10%.
Over the whole p T 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 p T , but decreases quickly with increasing p T as does the overall yield of this background source, as shown in figure 2 (left). At high p T , 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 p T interval and decreases to 12% at p T = 8 GeV/c. The systematic uncertainty due to the determination of the nucleon-nucleon cross section for the minimum-bias trigger condition is 3.7% [78].

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 p T 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 figure 1. The template fits were performed using the method proposed in [81]. 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 p T 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.

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 figure 1). The number of counts from a particular electron source j in a particular bin i is called a ji . Its unknown expectation value is called A ji 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 p j : f i = j p j A ji . 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 (d i ) within each impact parameter bin leading to the likelihood distribution [81] log This gives one free amplitude parameter for each electron source (p j ) and one free expectation value parameter for each electron source and impact parameter bin (A ji ). The main parameters of interest are the p j , in particular p beauty , while the nuisance parameters A ji 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 A ji for fixed p j is suggested in [81]. Numerical minimisation is then performed only for the p j . Equations for the iterative procedure can be found by setting the differentials dL/dA ji to zero. Solving these equations for A ji , yields an iterative rule for each bin.
For a bin i with a finite number of entries from data, but zero counts in any of the templates, the likelihood distribution of the A ji is not well represented by its maximum. This happens mostly in the tails of the distributions (see figure 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 p T interval, the number of electrons in the template was scaled by the amplitude parameter p beauty . 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 ( TPCeID ) was determined via a data-driven approach using electrons from photon conversions identified via topological cuts and the invariant mass [82]. The corresponding n TPC σ distributions were fitted with the function Landau · Exp ⊗ Gauss [72], 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 p T spectrum was unfolded.

JHEP07(2017)052
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 [76]. Due to the restricted p T 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 (1.1 < p T < 1.3 GeV/c and 8 < p T < 12 GeV/c) 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.

Systematic uncertainties estimation
The systematic uncertainties are summarized in table 2. They were estimated using datadriven methods where possible. An overview of the efficiencies of the different track selection steps may be found in figure 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 2%. The uncertainty for non-geometric effects was estimated to be smaller than 3%. 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 (−1.5 < n TPC σ < 4), in more peripheral collisions (20−40%, 40−80%), 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 3% for the requirement of charged tracks with good quality and about 10% 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 n TPC σ .
The uncertainty due to the modelling of the n TPC σ 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 5% for the TPC eID is the quadratic sum of the following contributions: 2% from the extrapolation, 2% from the pseudorapidity dependence, 3% from a possible p T dependence and 2% from the tail of the n TPC σ distribution.

JHEP07(2017)052
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 p T intervals with downward fluctuations of the measured charm yield. The systematic uncertainty varies between 19% and 6% between the different p T 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 (d 0 ) 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 p T 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 10%. The yield using a correction of 12% instead, differs by about 10% at p T = 1.3 GeV/c with the difference decreasing quickly towards higher p T . 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 n TPC σ distribution, suggesting a contamination of the electron candidate sample of about 15% 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 p T ≥ 5 GeV/c, 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 p T = 1.3 GeV/c decreasing towards higher transverse momentum. The proton contamination is significant only below p T = 1.3 GeV/c.
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 -16 -

JHEP07(2017)052
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 crosschecked using charged pions from K 0 S 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 p T = 1.3 GeV/c 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 [37]. The resulting difference in the yield of electrons from beauty-hadron decays was used as the estimate for the systematic uncertainty. It is 1.3% for low p T , decreasing towards higher transverse momentum.
Electrons at a fixed transverse momentum have mother particles in a range of p T 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 p T distribution [83]. 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 p T distribution of the templates and observing the resulting change in the measured electron p T distribution. The beauty-hadron p T 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 p T distribution was studied by introducing a momentum-dependent nuclear modification factor R AA . An R AA based on a theoretical calculation was used for the central points [84]. It has values near unity for low transverse momenta and drops to about 0.5 from a hadron p T of 5 to 10 GeV/c. This was varied to half its effect (R AA → (1 + R AA )/2) in order to estimate the associated uncertainty. For the charm case, the variation was done according to the measurement uncertainties [83]. The difference in the resulting measured yield of electrons from beauty-hadron decays is about 8%, with no visible p T 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 calcu- Associated uncertainty Tracking and matching 4.7% TOF matching and eID 10% TPC eID 5% Signal extraction 17% to 12% d 0 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 p T distribution 8% Baryon/meson ratio 5% Total 26% to 17% Table 2. Systematic uncertainties in the Pb-Pb analysis. Individual sources of systematic uncertainties are p T dependent, which is reported using intervals. The lower and upper value of the interval, respectively, lists the uncertainty at p T = 1.3 GeV/c and p T = 8.0 GeV/c. The second group of entries in the table is related to the method used to extract the electrons from beautyhadron decays.
lations [85]. This led to a change in the measured yield of electrons from beauty-hadron decays of about 5% with no clear momentum dependence. Decreasing the baryon ratio even to 0 has a smaller effect.
6 Reference pp cross sections at √ s = 2.76 TeV and √ s = 5.02 TeV For the calculations of the nuclear modification factors R pPb and R PbPb , corresponding pp reference spectra at √ s = 5.02 TeV and √ s = 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 √ s = 5.02 TeV exists. Therefore, the cross section of electrons from beauty-hadron decays measured in the momentum interval 1 < p T < 8 GeV/c at √ s = 7 TeV [23] was scaled to √ s = 5.02 TeV by applying a pQCD-driven √ sscaling [86]. The p T -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 √ s = 5.02 TeV and √ s = 7 TeV. Both the direct (b → e) and the cascade decay (b → c → e) were considered. For the calculations at both energies the same parameters were used for the beauty-quark mass (m b = 4.75 GeV/c 2 ), the PDFs (CTEQ6.6 [87]) as well as the factorisation µ F and renormalisation µ R scales with µ R = µ F = µ 0 = m 2 b + p 2 T,b , where p T,b denotes the transverse momentum of the beauty quark. The uncertainties of the p T -dependent scaling function were estimated by varying the parameters. The beauty-quark mass was set to m b = 4.5 and 5 GeV/c 2 . The uncertainties for the PDFs were obtained by using the CTEQ6.6 PDF uncertainties [87]. The contribution from the scale uncertainties was estimated by using six different sets: (µ R /µ 0 , µ F /µ 0 ) = (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 µ R = µ F = µ 0 . The resulting √ s-scaling uncertainty is almost independent of p T . It ranges from +4 −2 % at 1 GeV/c to about +2 −2 % at 8 GeV/c. The total systematic uncertainty of the pp reference spectrum at √ s = 5.02 TeV is then given as the bin-by-bin quadratic sum of the √ s-scaling uncertainty and the relative systematic uncertainty of the measured spectrum at √ s = 7 TeV. For the statistical uncertainties the relative uncertainties of the spectrum measured at √ s = 7 TeV were taken.
For the R PbPb analysis, the measured spectrum at √ s = 7 TeV was scaled to √ s = 2.76 TeV using FONLL pQCD calculations at the respective energies. The systematic scaling uncertainty is about +11 − 7 % at 1 GeV/c and about +7 −5 % at 8 GeV/c. The resulting pp reference spectrum was found to be consistent with the measurement of electrons from beauty-hadron decays in pp collisions at √ s = 2.76 TeV [88], shown in figure 4. The measured spectrum at √ s = 2.76 TeV was not taken as a reference for the R PbPb , because of larger statistical and systematic uncertainties than the reference obtained via the √ s-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 √ s = 7 TeV is 3.5% and also holds for the obtained pp reference spectra at √ s = 5.02 TeV and √ s = 2.76 TeV.
The systematic uncertainties of the input p T -differential cross section of electrons from beauty-hadron decays measured at √ s = 7 TeV, the normalisation uncertainty, as well as the scaling uncertainties for the reference spectra are summarised in table 3.
-19 -JHEP07(2017)052 pp spectrum 7 TeV 45% to 35% for 1 < p T < 1.5 GeV/c 35% to 20% for 1.5 < p T < 2.5 GeV/c ≤ 20% for p T ≥ 2.5 GeV/c Normalisation uncertainty 3.5% scaling uncertainty for p-Pb ( √ s = 5.02 TeV) Pb-Pb ( √ s = 2.76 TeV) Table 3. Systematic uncertainties of the p T -differential cross section of electrons from beautyhadron decays measured at √ s = 7 TeV [23], the normalisation uncertainty, as well as the scaling uncertainties for the reference spectra at √ s = 5.02 TeV and √ s = 2.76 TeV. The scaling uncertainties for the reference spectra are slightly p T dependent; the uncertainties are given for the two extreme p T intervals. Details are described in the text.

Results
The p T -differential cross section and invariant yield of electrons from beauty-hadron decays at mid-rapidity in minimum-bias p-Pb collisions at The nuclear modification factors R PbPb and R pPb are shown in figure 6 (left). The R PbPb was obtained using equation (1.1). The R pPb 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 p T bins. The normalisation uncertainty of the pp spectrum and the uncertainty of the nuclear overlap function T AA 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 figure 6.
The R pPb is consistent with unity within uncertainties (of about 20% for p T > 2 GeV/c) 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 R PbPb for the 20% most  [17, 18, 29-31, 57, 84, 90-97], see text for details. The vertical bars represent the statistical uncertainties, while the boxes indicate the systematic uncertainties. The normalisation uncertainty, common to all points, is shown as a filled box at high p T for both collision systems.
central Pb-Pb collisions increase, for p T ≤ 3 GeV/c, with sizeable uncertainties of 30-45%. In the interval 3 < p T < 6 GeV/c, the R PbPb is about 0.7 with a systematic uncertainty of about 30%; in 6 < p T < 8 GeV/c the ratio is 0.48 with an uncertainty of about 25%. In the latter transverse momentum range the suppression with respect to R PbPb = 1 is a 3.3σ effect taking into account the statistical and systematic uncertainties.
A comparison of the R PbPb of electrons from beauty-hadron decays with the one from charm-and beauty-hadron decays is shown in figure 6 (right) for the 20% most central Pb-Pb collisions. For the latter R PbPb , the p T -differential invariant yields of electrons from charm-and beauty-hadron decays published in [89] for the centrality classes 0-10% and 10-20% were combined. For the pp reference in the momentum range up to p T ≤ 12 GeV/c, the corresponding invariant cross section measurement at √ s = 2.76 TeV [24], which has uncertainties of about 20%, was used. For p T ≥ 12 GeV/c, the ATLAS measurement [72] at √ s = 7 TeV was extrapolated to √ s = 2.76 TeV applying a FONLL pQCD-driven √ sscaling 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 p T , where the beauty contribution is larger than the charm contribution [24]. In the p T interval 3 < p T < 6 GeV/c, the suppression of the R PbPb 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 R pPb is described by pQCD calculations including modifications of the parton distribution functions (FONLL [29][30][31] + EPS09NLO [90] nuclear PDFs) as shown in figure 7 (left). The data and the calculation suggest that cold nuclear matter effects are small at high transverse momentum. Recent measurements of long-range -22 -

JHEP07(2017)052
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 [73] 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. [57], in which the p T 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 p T -spectra of light hadrons [73] 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 [17]. An enhancement at intermediate p T is predicted by the calculations based on incoherent multiple scattering [18]. 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 √ s = 2.76 TeV (FONLL [29-31] + EPS09NLO [90] nuclear PDFs) cannot describe the R PbPb at high transverse momentum (see figure 7, right), indicating that the suppression, particularly evident in the interval 6 < p T < 8 GeV/c, 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 R PbPb 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 [91], the partonic transport description BAMPS [96,97], and in WHDG [93][94][95]. The non-perturbative transport model TAMU [84] includes only collisional processes, while the POWLANG [92] 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 [95].
The right-hand side of figure 7 shows the comparison of the various models with the measured R PbPb . The MC@sHQ+EPOS2 calculation with EPOS initial conditions [98,99], including the Landau-Pomeranchuk-Migdal (LPM) effect [100], is consistent with the data at high p T . 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 p T 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 p T spectra and elliptic flow measurements of light hadrons. The calculations are consistent with the data at high p T , indicating a limited sensitivity of the current data to radiative energy loss effects. The POWLANG [92] transport calculation takes into account initial-state nuclear effects via -23 -

JHEP07(2017)052
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 [101] 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 T dec = 155 MeV; the results with a temperature of T dec = 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 R PbPb .
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 R PbPb values around unity or slightly above. The data show a larger increase with decreasing transverse momentum, however exhibit large systematic and statistical uncertainties.

Summary
The p T -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 √ s-scaling of the cross section of electrons from beauty-hadron decays measured at √ s = 7 TeV. The R pPb is consistent with unity within uncertainties of about 20% at high transverse momentum p T , which increase towards low p T . The R pPb 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 R PbPb is about 0.7 with an uncertainty of about 30% in the interval 3 < p T < 6 GeV/c and 0.48 with an uncertainty of about 25% for 6 < p T < 8 GeV/c. 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 R PbPb increase for p T ≤ 3 GeV/c with uncertainties of about 30-45%. The measured R PbPb is described within uncertainties by pQCD-inspired models of beauty-quark energy loss in the QGP. In the interval 3 < p T < 6 GeV/c, we observe that the suppression of the R PbPb 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.

Acknowledgments
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 Collab- [28] ATLAS collaboration, Measurement of the differential cross-sections of inclusive, prompt and non-prompt J/ψ production in proton-proton collisions at -27 -