Measurement of azimuthal correlations of D mesons with charged particles in pp collisions at \(\sqrt{s}=7\) TeV and p–Pb collisions at \({\sqrt{{{s}}_{\scriptscriptstyle {\mathrm{NN}}}}}=5.02\) TeV

The European Physical Journal C, Apr 2017

The azimuthal correlations of D mesons with charged particles were measured with the ALICE apparatus in pp collisions at \({\sqrt{s}}=7~\mathrm {TeV}\) and p–Pb collisions at \({\sqrt{{{s}}_{\scriptscriptstyle {\mathrm{NN}}}}}=5.02~\mathrm {TeV}\) at the Large Hadron Collider. \(\mathrm{D}^{0}\), \({\mathrm{D}^{+}}\), and \({\mathrm{D}^{*+}}\) mesons and their charge conjugates with transverse momentum \(3<{p}_{\mathrm{T}}<16~\mathrm {GeV}/c\) and rapidity in the nucleon-nucleon centre-of-mass system \(|y_\mathrm{cms}|<0.5\) (pp collisions) and \(-0.96<y_\mathrm{cms}<0.04\) (p–Pb collisions) were correlated to charged particles with \({p}_{\mathrm{T}}>0.3~\mathrm {GeV}/c\). The yield of charged particles in the correlation peak induced by the jet containing the D meson and the peak width are compatible within uncertainties in the two collision systems. The data are described within uncertainties by Monte-Carlo simulations based on PYTHIA, POWHEG, and EPOS 3 event generators.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-017-4779-8.pdf

Measurement of azimuthal correlations of D mesons with charged particles in pp collisions at \(\sqrt{s}=7\) TeV and p–Pb collisions at \({\sqrt{{{s}}_{\scriptscriptstyle {\mathrm{NN}}}}}=5.02\) TeV

Eur. Phys. J. C Measurement of azimuthal correlations of D mesons with charged √ particles in pp collisions at s = 7 TeV and p-Pb collisions at √sNN = 5.02 TeV ALICE Collaboration 0 0 CERN , 1211 Geneva 23 , Switzerland The azimuthal correlations of D mesons with charged particles were measured with the ALICE apparatus in pp collisions at √s = 7 TeV and p-Pb collisions at √sNN = 5.02 TeV at the Large Hadron Collider. D0, D+, and D∗+ mesons and their charge conjugates with transverse momentum 3 < pT < 16 GeV/c and rapidity in the nucleonnucleon centre-of-mass system |ycms| < 0.5 (pp collisions) and −0.96 < ycms < 0.04 (p-Pb collisions) were correlated to charged particles with pT > 0.3 GeV/c. The yield of charged particles in the correlation peak induced by the jet containing the D meson and the peak width are compatible within uncertainties in the two collision systems. The data are described within uncertainties by Monte-Carlo simulations based on PYTHIA, POWHEG, and EPOS 3 event generators. 1 Introduction The study of the angular correlation of D mesons with charged particles, i.e. the distribution of the differences in azimuthal angles, ϕ = ϕch − ϕD, and pseudorapidities, η = ηch − ηD, allows for the characterisation of charm production and fragmentation processes in proton–proton (pp) collisions and of their possible modifications due to nuclear effects in proton–Pb and Pb–Pb collisions [1,2]. For leadingorder (LO) Quantum-ChromoDynamic (QCD) processes, charm quark–antiquark pairs are produced back-to-back in azimuth: the angular correlation of D mesons with charged particles features a “near-side” peak at ( ϕ, η) = (0, 0) and an “away-side” peak at ϕ = π . The former originates from the jet containing the “trigger” D meson, the latter is induced by the recoil jet, which can also include the decay products of the other charmed hadron produced in the collision. The away-side peak extends over a wide range in η. The two peaks lie on top of an approximately flat distribution arising from the correlation of D mesons with charged particles from the underlying event. Next-to-leading order (NLO) production processes can give rise to significantly different correlation patterns [3,4]. For example, the radiation of a hard gluon from a charm quark smears the back-to-back topology of LO production and broadens both the near- and the away-side peak. In addition, quark–antiquark charm pairs originating from the splitting of a gluon can be rather collimated and, especially at high transverse momentum ( pT), can generate sprays of hadrons contributing to a unique and broader “near-side” peak. In such cases, the away-side peak stems from the particles coming from the fragmentation of the recoil parton (typically a gluon or a light quark), which is not aligned with the trigger D meson. Finally, in hardscattering topologies classified as “flavour excitation” (see e.g. [4]) a charm quark (antiquark) from an initial splitting g → cc¯ undergoes a hard interaction. The hadrons originating from the antiquark (quark) can be significantly separated in rapidity with respect to the trigger D meson and contribute with a rather flat term to the ϕ distribution. Since the first measurement performed by STAR in Au–Au collisions at √sNN = 200 GeV [5], two-particle azimuthal correlations have been exploited at both RHIC and the LHC [6–8] to investigate the possible modifications of jet and dijet properties that can be caused by the interaction of high-energy partons with the constituents of the Quark Gluon Plasma (QGP) formed in ultra-relativistic heavy-ion collisions. The most evident effect is the suppression of the awayside correlation peak, commonly attributed to in-medium partonic energy loss. The results allow one to constrain the dependence of the energy loss on the distance covered by partons in the QGP as well as the initial gluon density [9,10]. The correlation pattern of hadron–hadron pairs primarily arises from the back-to-back production of gluons or light-quarks produced in hard-scattering processes, and their subsequent fragmentation. PHENIX measured the azimuthal correlation of electrons from heavy-flavour hadron decays with charged particles in Au–Au collisions at √sNN = 200 GeV [11]. The near- and away-side peaks are suppressed by factors compatible, within uncertainties, to those observed for hadron– hadron correlations, if a similar pT is considered for the trigger hadron and the electron parent hadron. The proper interpretation of nucleus–nucleus results and the connection of the modifications of the correlation peak properties to the parton dynamics in the QGP requires the comparison of data with model predictions. It is crucial that the models reproduce the correlation pattern measured in pp collisions, where nuclear effects are absent, as well as the production spectra in both pp and nucleus–nucleus collisions. Therefore, the measurement of azimuthal correlations of D mesons with charged particles in pp and p–Pb collisions serves not only as a reference for future measurements in Pb–Pb collisions but it also allows for validation of Monte-Carlo generator expectations, which is fundamental for the understanding of the results in all collision systems. Perturbative QCD calculations relying on the collinearfactorisation approach, like FONLL [12] and GM-VFNS [13], or based on the kT-factorisation approach [14] describe reasonably well the pT-differential production cross sections of D mesons from charm-quark fragmentation measured at central rapidity (|y| < 0.5) in pp collisions at √s = 7 and 2.76 TeV using the ALICE detector [15,16]. These calculations represent the state of the art for the computation of ( pT, y)-differential cross sections of charm quarks and charmed hadrons. However, the kinematic relationship between D mesons and particles from charm fragmentation and the underlying event is accessible only with event generators coupled with parton-shower Monte-Carlo programs like those provided by PYTHIA [17] and HERWIG [18]. The order of hard-scattering matrix elements used, the specific implementation of parton shower and hadronisation, as well as the modeling of the underlying event have an influence on the angular correlations of D mesons with charged particles produced in the event. For heavy quarks with mass M and energy EQ, the suppression of gluon radiation off the quark inside the forward cone with opening angle = M/EQ (the so-called “dead-cone” effect) reduces the phase space for primary gluon radiation [19]. This implies a harder fragmentation of the quarks into the heavy hadrons and leads to essential differences in the profiles of gluon-, light-quark- and heavy-quark-initiated jets resulting in shape differences of pT-spectra and multiplicity distributions of primary hadrons in the jets [20,21]. Correlations between D mesons were measured at the LHC in pp collisions at √s = 7 TeV with the LHCb experiment [22], providing information on charm production mechanisms and on the properties of events containing heavy flavours. ATLAS measured the production of D∗+ mesons in jets in pp collisions at √s = 7 TeV for jets with 25 < pT < 70 GeV/c and D∗+ carrying a jet momentum fraction (z) in the range 0.3 < z < 1. The results indicate that the production of charm-quark jets or charm-quark fragmentation into D∗+ mesons is not properly modeled in state-of-the-art Monte-Carlo generators [23]. Azimuthal correlations of electrons from heavy-flavour hadron decays with charged particles were also exploited to study the relative beauty contribution to the population of electrons from heavy-flavour hadron decays in pp collisions at RHIC and at the LHC [24,25]. The angular distribution of particles produced in an event is sensitive to collective effects that correlate particle production over wide phase-space regions. This is particularly relevant in Pb–Pb collisions with non-zero collision impact parameter, where the azimuthal asymmetry of the overlapping region of the colliding nuclei gives rise to anisotropic pressure gradients inducing an anisotropy in the azimuthal distribution of particle momenta [26,27]. The main component of the Fourier decomposition used to describe the resulting ϕ distribution of two particle correlations is the 2nd order term, proportional to cos(2 ϕ), called elliptic flow or v2. Given that correlations induced by the collective motion of the system extend over large pseudorapidity ranges, the elliptic-flow term manifests itself with the presence of two long-range ridge-like structures in the near and away sides of two-particle angular correlations. Unexpectedly, similar long-range correlation structures were observed in high-multiplicity pp and p–Pb collisions at the LHC [28– 33]. Also in central d–Au collisions at RHIC [34,35] similar results were obtained, although contributions from jet-like correlations due to biases on the event selection could not be excluded [36]. The origin of such v2-like structures is still debated. Positive v2 values in high-multiplicity pp collisions and p–Pb (d–Au) collisions at LHC (RHIC) are expected in models that include final-state effects [37–41], as well as initial-state effects related to the Color Glass Condensate [42] or to gluon bremsstrahlung by a quark–antiquark string [43]. A modification of the azimuthal correlations of D mesons with charged particles in p–Pb with respect to pp collisions could be a signal of the presence of longrange v2-like correlations for particles originating from hardscattering processes. This would yield complementary information to that obtained from correlations of light-flavour particles, which at low pT are primarily produced in soft processes. The D-meson pT-differential production cross section in p–Pb collisions at √sNN = 5.02 TeV was measured with ALICE in the interval of rapidity in the nucleon-nucleon centre-of-mass system −0.96 < ycms < 0.04 [44]. The data are compatible, within uncertainties, with a Glauber-modelbased geometrical scaling of a pp collision reference obtained f√rosm=t2h.e76crToesVs .sTehcitsiosnusggmesetasstuhraetdnuatcle√asr ef=fect7s aTreeVratahnedr small for D mesons in the range 1 < pT < 24 GeV/c. However, they could still affect angular correlations as observed at RHIC for azimuthally-correlated pairs of electrons and muons from decays of heavy-flavour hadrons in d–Au collisions at √sNN = 200 GeV [45]. A modification of the azimuthal correlation of heavy-flavour particles in p–Pb collisions could occur at the LHC due to gluon saturation in the heavy nucleus [46]. Moreover, transport models based on the Langevin equation [2,47] describe, within uncertainties, the nuclear modification factor of D mesons measured in p–Pb collisions at the LHC and that of electrons from heavy-flavour hadron decays measured in d–Au collisions at RHIC [48]. These models assume the formation of a smallsize QGP in p–Pb and d–Au collisions and include the possibility of heavy-flavour hadron formation via coalescence of heavy quarks with thermalised light quarks from the medium. These transport calculations predict a positive D-meson v2 in central p–Pb collisions. As an example, in the case of the POWLANG model [2] the maximum expectation for the 20% most central p–Pb collisions is v2 ∼ 5% at pT = 4 GeV/c. A finite v2 of muons from heavy-flavour hadron decays in high-multiplicity p–Pb collisions was also suggested in [31] as one of the possibilities for reconciling the measured values of v2 of inclusive muons with the expectations based on the multi-phase transport model AMPT [49]. In this paper we report the first measurements of azimuthal correlations of D mesons with charged primary particles in pp and p–Pb collisions at √s = 7 TeV and √sNN = 5.02 TeV, respectively. Unless differently specified we always refer to “prompt” D mesons from charm-quark fragmentation. In what follows, primary particles are defined as particles originated at the collision point, including those deriving from strong and electromagnetic decays of unstable particles, and those from decays of hadrons with charm or beauty. The paper is organised as follows. In Sect. 2 the data samples used and the details of the ALICE experimental apparatus relevant for this analysis are described. The analysis strategy, the Dmeson signal extraction, the associated-track selection criteria, and the corrections applied to measure the correlations between D mesons and charged primary particles are reported in Sect. 3. In the same section, the fit procedure adopted to quantify the correlation peak properties is described. Section 4 reports the systematic uncertainties affecting the measurement. The results are discussed in Sect. 5. The paper is then summarised in Sect. 6. 2 Experimental apparatus and data samples 2.1 The ALICE detector and event selection The ALICE apparatus [50,51] consists of a central barrel embedded in a 0.5 T solenoidal magnetic field, a forward muon spectrometer, and a set of detectors located in the forward- and backward-rapidity regions dedicated to trigger and event characterisation. The analysis reported in this paper is performed using the central barrel detectors. Charged particle tracks are reconstructed using the Inner Tracking System (ITS), consisting of six layers of silicon detectors, and the Time Projection Chamber (TPC). Particle identification (PID) is based on the specific energy loss dE /dx in the TPC gas and on the time of flight from the interaction vertex to the Time-Of-Flight (TOF) detector. The ITS, TPC and TOF detectors provide full azimuthal coverage in the pseudorapidity interval |η| < 0.9. The pp data sample consists of about 3 × 108 minimumbias events, corresponding to an integrated luminosity of Lint = 5 nb−1. These collisions are triggered by the presence of at least one hit in one of the V0 scintillator arrays, covering the ranges −3.7 < η < −1.7 and 2.8 < η < 5.1, or in the Silicon Pixel Detector (SPD), constituting the two innermost layers of the ITS, with an acceptance of |η| < 2 (inner layer) and |η| < 1.4 (outer layer). The p–Pb data sample consists of about 108 minimum-bias events, corresponding to an integrated luminosity of about Lint = 50 µb−1. In this case the minimum-bias trigger requires signals in both the V0 detectors. Only events with a reconstructed primary vertex within ±10 cm from the centre of the detector along the beam line are considered for both pp and p–Pb collisions. This choice maximises the detector coverage of the selected events, considering the longitudinal size of the interaction region, and the detector pseudorapidity acceptances (for more details see [51]). For p–Pb collisions, the center-of-mass reference frame of the nucleon-nucleon collision is shifted in rapidity by yNN = 0.465 in the proton direction with respect to the laboratory frame, due to the different per-nucleon energies of the proton and the lead beams. Beam-gas events are removed by offline selections based on the timing information provided by the V0 and the Zero Degree Calorimeters (two sets of neutron and proton calorimeters located around 110 m from the interaction point along the beam direction), and the correlation between the number of hits and track segments in the SPD detector. The minimum-bias trigger efficiency is 100% for events with D mesons with pT > 1 GeV/c for both pp and p–Pb data sets. For the analyzed data samples, the probability of pile-up from collisions in the same bunch crossing is below 4% per triggered pp event and below the percent level per triggered p–Pb event. Events in which more than one primary interaction vertex is reconstructed with the SPD detector are rejected, which effectively removes the impact of in-bunch pile-up events on the analysis. The contribution of particles from pile-up of pp collisions in different bunch crossings is also negligible due to the selections applied to the tracks used in this analysis and the large interval between subsequent bunch crossings in the data samples used. 2.2 Monte-Carlo simulations Monte-Carlo simulations including a complete description of the ALICE detector are used to calculate the corrections for the azimuthal-correlation distributions evaluated from data. The distribution of the collision vertex along the beam line, the conditions of all the ALICE detectors, and their evolution with time during the pp and p–Pb collision runs are taken into account in the simulations. Proton-proton collisions are simulated with the PYTHIA 6.4.21 event generator [17] with the Perugia-0 tune (tune number 320) [52] while p–Pb collisions are simulated using the HIJING v1.36 event generator [53]. For the calculation of D-meson reconstruction efficiencies PYTHIA simulations of pp collisions are used, requiring that in each event a cc¯ or b b¯ pair is present. In the simulation used for the analysis of p–Pb data, a p–Pb collision simulated with HIJING is added on top of the PYTHIA event. The generated particles are transported through the ALICE apparatus using the GEANT3 package [54]. The measured angular-correlation distributions are compared to simulation results obtained with the event generators PYTHIA 6.4.25 [17] (tunes number 320, 327, and 350, corresponding to the reference versions of the Perugia-0, Perugia2010, and Perugia-2011 sets [52], respectively), PYTHIA 8.1 (tune 4C) [55], POWHEG [56, 57] coupled to PYTHIA (Perugia-2011 tune), and EPOS 3.117 [58–60] (referred to as EPOS 3 hereafter). PYTHIA simulations utilise LOpQCD matrix elements for 2 → 2 processes, along with a leading-logarithmic pT-ordered parton shower, the Lund string model for hadronisation, and an underlying-event simulation including Multiple-Parton Interactions (MPI). With respect to older tunes, the Perugia tunes use different initialstate radiation and final-state radiation models. One of the main differences is that the parton shower algorithm is based on a pT-ordered evolution rather than a virtuality-ordered one. Significant differences in the treatment of colour reconnection, MPI, and the underlying event were also introduced. Perugia 0 is the first of the series. The Perugia-2010 tunes differ from those of Perugia-0 in the amount of final-state radiation and by a modification of the high-z fragmentation (inducing a slight hardening of the spectra). They are expected to better reproduce observables related to the jet shape. The first LHC data, mainly from multiplicity and underlying-event related measurements, were considered for the Perugia-2011 tunes. PYTHIA 8.1 also includes several improvements in the treatment of MPI and colour reconnection [55]. In the simulations done with √s = 5.02 TeV, the centre-of-mass frame is boosted in rapidity by yNN = 0.465 in order to reproduce the rapidity shift of the reference frame of the nucleon-nucleon collision in the p–Pb collision system. POWHEG is a NLO-pQCD generator [56, 57] that, coupled to parton shower programs (e.g. from PYTHIA or HERWIG [18]), can provide exclusive final-state particles, maintaining the next-to-leading order accuracy for inclusive observables. The charm-production cross sections obtained with POWHEG+PYTHIA are consistent with FONLL [12] and GM-VFNS [13] calculations within the respective uncertainties, and are in agreement with measured D-meson production cross sections within the model and experimental uncertainties [61, 62]. The POWHEG+PYTHIA simulations presented in this paper are obtained with the POWHEG BOX framework [63, 64] and the tune Perugia 2011 of PYTHIA 6.4.25. For the comparison with the measured p–Pb collision data, parton distribution functions (PDFs) corrected for nuclear effects (CT10nlo [65] with EPS09 [66]) are used. In addition, a boost in rapidity by yNN = 0.465 is applied to the partons generated with POWHEG before the PYTHIA parton shower process. EPOS 3 [58–60] is a Monte-Carlo event generator based on a 3+1D viscous hydrodynamical evolution starting from flux tube initial conditions, which are generated in the Gribov-Regge multiple-scattering framework. Individual scatterings are referred to as Pomerons, and are identified with parton ladders. Each parton ladder is composed of a pQCD hard process with initial and final state radiation. Non-linear effects are considered by means of a saturation scale. The hadronisation is performed with a string fragmentation procedure. Based on these initial conditions, the hydrodynamical evolution can be applied on the dense core of the collision. An evaluation within the EPOS 3 model shows that the energy density reached in pp collisions at √s = 7 TeV is high enough to apply such hydrodynamic evolution [60]. 3 Data analysis The analysis procedure consists of three main parts, which are described in the following subsections: D-meson reconstruction and selection of primary particles to be used in the correlation analysis (Sect. 3.1), construction of azimuthalcorrelation distribution and corrections, including the subtraction of combinatorial background and beauty feed-down contributions (Sect. 3.2), extraction of correlation properties via fits to the azimuthal distributions (Sect. 3.3). 3.1 D-meson and associated-particle reconstruction The correlation analysis is performed by associating D mesons (D0, D+, D∗+ mesons and their antiparticles), defined as “trigger” particles, with charged primary particles in the same event, and excluding those coming from the decay of the trigger D mesons themselves. The D0, D+, D∗+ mesons and their charge conjugates are reconstructed via their hadronic decay channels D0 → K−π +, with Branching Ratio (BR) of (3.88±0.05)%, D+ → K−π +π +, BR of (9.13±0.19)%, and D∗+ → D0π +, BR of (67.7±0.5)% [67]. The D-meson signal extraction is based on the reconstruction of decay vertices displaced from the primary vertex by a few hundred microns and on the identification of the decayparticle species. The same selection procedures used for the measurements of D-meson production in pp and p–Pb collisions at √s = 7 TeV and √sNN = 5.02 TeV, respectively, are adopted [15,44]. For both the pp and p–Pb data sets, D0 and D+ candidates are formed by combining two or three tracks, respectively, with each track satisfying |η| < 0.8 and pT > 0.3 GeV/c. Additionally, D0 and D+ daughter tracks are required to have at least 70 out of a maximum of 159 possible associated space points in the TPC, a χ 2/NDF of the momentum fit in the TPC smaller than 2, and at least 2 out of 6 associated hits in the ITS. D∗+ candidates are formed combining D0 candidates with tracks with one point in the SPD, |η| < 0.8 and pT > 0.1 GeV/c. The main variables used to reject the combinatorial background are the separation between primary and secondary vertices, the distance of closest approach (DCA) of the decay tracks to the primary vertex, and the angle between the reconstructed D-meson momentum and the flight line defined by the primary and secondary vertices. A tighter selection is applied for p–Pb collisions with respect to pp collisions to reduce the larger combinatorial background. Charged kaons and pions are identified using the TPC and TOF detectors. A ±3σ cut around the expected value for pions and kaons is applied on both TPC and TOF signals. The D mesons are selected in a fiducial rapidity range varying from |ylab| < 0.5 at low pT to |ylab| < 0.8 for D mesons with pT > 5 GeV/c in order to avoid cases in which the decay tracks are close to the edge of the detector, where the acceptance decreases steeply. The D0 and D+ raw yields are extracted using fits to the distributions of invariant mass M (K−π +) and M (K−π +π +), respectively, with a function composed of a Gaussian term for the signal and an exponential term that models the combinatorial background. In the case of the D∗+, the raw yield is obtained by fitting the invariant-mass difference M = M (K−π +π +) − M (K−π +), using a Gaussian function for the signal and a threshold function multiplied by an exponential (a√ M − Mπ ·eb( M−Mπ )) to describe the background. Relatively wide D-meson pT intervals (3 < pT < 5 GeV/c, 5 < pT < 8 GeV/c, 8 < pT < 16 GeV/c for pp collisions and 5 < pT < 8 GeV/c, 8 < pT < 16 GeV/c for p–Pb collisions) are chosen to reduce the statistical fluctuations in the azimuthal-correlation distributions. Figure 1 shows the D0 and D+ invariant mass, and D∗+ invariant-mass difference distributions in the 3 < pT < 5 GeV/c interval for pp collisions and in the 5 < pT < 8 GeV/c, 8 < pT < 16 GeV/c intervals for p–Pb collisions. The fits used to evaluate the raw yields are also shown. The statistical uncertainty of the D-meson raw yields in the pT intervals analyzed varies from about 5 to 8% (3 to 5%) in pp (p–Pb) collisions for the D0 and D+ mesons and from about 5 to 6% (5 to 10%) for the D∗+ mesons, depending on pT. For both collision systems, the signal over background ratio of the signal peaks is between 0.2 and 1 for the D0 and D+ mesons, and up to 2.6 for the D∗+ meson. In the interval 3 < pT < 5 GeV/c the D-meson yield can be extracted from the invariant mass distribution with statistical uncertainty smaller than 3% in both pp and p–Pb collisions. However, in the latter case, the near- and away-side peaks of the azimuthal-correlation distribution, that have a small amplitude at low D-meson pT, cannot be disentangled from the statistical fluctuations of the baseline, which is related to the multiplicity of the event and thus higher in p–Pb than in pp collisions. Therefore, for this pT interval, the results are shown only for pp collisions. Associated particles are defined as all charged primary particles with passoc > 0.3 GeV/c and with pseudorapid T ity |η| < 0.8, except for the decay products of the trigger D meson. Particles coming from other weak decays or originating from interactions with the detector material are defined as secondary particles and are discarded. Reconstructed tracks with at least 70 points in the TPC and 3 in the ITS, and a χ 2/NDF of the momentum fit in the TPC smaller than 2 are associated to D-meson candidates. Using Monte Carlo simulations (see Sect. 2.2), these selection criteria yield an average track reconstruction efficiency for charged primary particles of about 85% in the pseudorapidity range |η| < 0.8 and in the interval 0.3 < pT < 24 GeV/c, with variations contained within ≈5% for pT < 1.5 GeV/c. Negligible variations are observed at higher pT. The contamination of secondary particles is removed by requiring the DCA of the associated tracks to the primary vertex to be less than 2.5 mm in the transverse (x , y) plane and less than 1 cm along the beam line (z direction). This selection identifies primary particles with a purity ( pprim) of approximately 96% and an efficiency higher than 99%, also for particles originating from decays of charm or beauty hadrons, which can be displaced by several hundred micrometers from the primary vertex. The purity is independent of pT in the measured pT range. For the D0-meson case, the low- pT pion produced from the D∗+ → D0π + decay is removed from the sample of associated particles by rejecting tracks that yield a M compatible within 3σ with the value expected for D∗+ mesons. It was verified with Monte Carlo simulations that this selection rejects more than 99% of the pions from D∗+ decays in all D-meson pT intervals considered and has an efficiency larger than 99% for primary particles with pT > 0.3 GeV/c. 0 1.75 1.8 1.85 1.9 1.95 2 2.05 Invariant Mass (Kπ) (GeV/c2) 100 p-Pb, sNN = 5.02 TeV 0 1.75 1.8 1.85 1.9 1.95 2 2.05 Invariant Mass (Kπ) (GeV/c2) p-Pb, sNN = 5.02 TeV 0 1.75 1.8 1.85 1.9 1.95 2 2.05 Invariant Mass (Kπ) (GeV/c2) /c 500 V e /8M400 s e i trn 300 E Fig. 1 Distributions of D0 (left column) and D+ (middle column) candidate invariant mass and of the D∗+ candidate invariant-mass difference (right column). The distributions are shown for pp collisions in the 3 < pT < 5 GeV/c range (top row) and for p–Pb collisions in the 3.2 Azimuthal-correlation distributions and corrections D-meson candidates with invariant mass (M ) in the range |M − μ| < 2σ (peak region), where μ and σ denote the mean and width of the Gaussian term of the invariant-mass fit function, are correlated to tracks selected with the criteria described above, and the difference in the azimuthal angle ( ϕ) and in pseudorapidity ( η) of each pair is computed. In order to correct for the acceptance and reconstruction efficiency (Acc × ε) of the associated tracks and for the variation of (Acc × ε) of prompt D mesons inside a given pT interval, a weight equal to the inverse of the product of both (Acc × ε) is assigned to each pair. The dependence of the 5 < pT < 8 GeV/c (middle row) and 8 < pT < 16 GeV/c (bottom row) ranges. The fits to the invariant mass distributions and the Gaussian mean and sigma values are also shown associated-track efficiency on transverse momentum, pseudorapidity, and position of the primary vertex along the beam axis is taken into account. The dependence of the track reconstruction efficiency on the event multiplicity is negligible and therefore neglected. The reconstruction efficiency of prompt D mesons is calculated as a function of pT and event multiplicity. It is on the order of few percent in the lowest D-meson pT interval, about 20% at high pT [15,44], and varies within each pT interval by up to a factor 2–3 (1.5–2) at low (high) pT, depending on the D-meson species and collision system. The D-meson (Acc×ε) factor also accounts for the pT-dependent fiducial rapidity range of the selected D mesons (Sect. 3.1) in order to normalise the results to one unit of rapidity. The obtained distribution, C ( ϕ, η)peak, also includes the angular correlation of combinatorial D-meson candidates in the peak range, which is a source of background and needs to be subtracted. This contribution is estimated via the per-trigger correlation distribution of background candidates in the sideband invariant-mass range, 1/Bsidebands × C ( ϕ, η)sidebands, where Bsidebands is the amount of background in the sideband region 4σ < |M − μ| < 8σ (right side only, 4σ < M − μ < 15σ , in the case of D∗+ mesons). The term C ( ϕ, η)sidebands represents the correlation distribution obtained as described above, but selecting trigger D-meson candidates with invariant mass in the sidebands. The background contribution is then subtracted from C ( ϕ, η)peak after being normalised to the amount of combinatorial background in the peak region, Bpeak. The latter is obtained from the counts in the invariant-mass distribution in the peak region, after subtracting the signal, Speak, estimated from the invariant-mass fit. Note that Speak, Bpeak and Bsidebands are calculated from the invariant-mass distributions weighted by the inverse of the prompt D-meson reconstruction efficiency. The correlation distributions C ( ϕ, η)peak and C ( ϕ, η)sidebands are corrected for the limited detector acceptance and spatial inhomogeneities using the event mixing technique. In this approach, D-meson candidates found in a given event are correlated with charged tracks from other events with similar multiplicity and primary-vertex position along the beam axis. The distribution obtained from the mixed events, ME( ϕ, η), shows a typical triangular shape as a function of η, due to the limited η coverage of the detector, and is approximately flat as a function of ϕ. The event-mixing distribution is rescaled by its average value in the range (−0.2 < ϕ < 0.2,−0.2 < η < 0.2) and its inverse is used as a map to weight the distributions C ( ϕ, η)peak and C ( ϕ, η)sidebands. A correction for the purity of the primary-particle sample ( pprim, see Sect. 3.1) is applied and the per-trigger normalisation is obtained dividing by Speak. The above procedure is summarised in Eq. 1, where the notation C˜ refers to angular-correlation distributions normalised by the number of trigger particles: C( ϕ, η) Bpeak C( ϕ, η) ME( ϕ, η) peak − Bsidebands ME( ϕ, η) sidebands Finally, the per-trigger azimuthal distribution C˜ inclusive ( ϕ) is obtained by integrating C˜ inclusive( ϕ, η) in the range | η| < 1. It was verified using Monte-Carlo simulations based on PYTHIA (Perugia-2011 tune) that the per-trigger azimuthal correlation of D mesons and secondary particles not rejected by the track selection has a ϕ-dependent modulation with a maximum variation of 7% with respect to the azimuthal correlation of D mesons and primary particles. This ϕdependent contamination has a negligible impact on the final results, considering the 4% level of contamination of secondary particles in the sample of associated tracks, hence, it was neglected. A fraction of the reconstructed D mesons consists of secondary D mesons coming from B-meson decays. The topological cuts, applied to reject combinatorial background, preferentially select displaced vertices, yielding a larger (by about a factor 2 for D0 mesons in the measured pT range) efficiency for secondary D mesons than for prompt D mesons. Therefore, the fraction fprompt of reconstructed prompt D mesons does not coincide with the natural fraction and depends on the analysis details. The different fragmentation, as well as the contribution of B-meson decay particles and a possible different contribution of gluon splitting to charm- and beauty-quark production, imply a different angular-correlation distribution of prompt and secondary D mesons with charged particles, as it was verified with the Monte-Carlo simulations described in Sect. 2.2. The contribution of feed-down D mesons to the measured angular correlation is subtracted as follows: −(1 − fprompt)C˜ fMeeCd-tdeomwpnl( ϕ)). In Eq. 2, C˜ prompt( ϕ) is the per-trigger azimuthal-correlation distribution after the subtraction of the feed-down contribution, fprompt is the fraction of prompt D mesons and C˜ fMeeCd −tedmopwln( ϕ) is a template for the azimuthal-correlation distribution of the feed-down component. Using the same method described in [15], fprompt was evaluated on the basis of FONLL calculations of charm and beauty pT-differential production cross sections [12] and of the reconstruction efficiencies of prompt and secondary D mesons, calculated using Monte-Carlo simulations. The value of fprompt, which depends on the D-meson species and varies as a function of pT, is estimated to be larger than 75%. The azimuthal correlation of feed-down D mesons, C˜ fMeeCd −tedmopwln, was obtained from PYTHIA (tune Perugia 2011 [52]) simulations of pp collisions at √s = 7 TeV and √s = 5.02 TeV for the analysis of pp and p–Pb data, respectively. In order to avoid biases related to the different event multiplicity in real and simulated events, the correlation distribution was shifted to have its minimum coinciding with the baseline of the data azimuthalcorrelation distribution before feed-down subtraction. A difference smaller than 8% was observed in the simulation between the baseline values of the azimuthal-correlation distributions for prompt and feed-down D mesons. Considering the typical values of fprompt, this difference results in a shift of the baseline of C˜ prompt( ϕ) smaller than 2%, negligible with respect to the other uncertainties affecting the measurement. 3.3 Characterization of azimuthal-correlation distributions In order to quantify the properties of the measured azimuthal correlations, the following fit function is used: f ( ϕ) = b + √ It is composed of two Gaussian terms describing the nearand away-side peaks and a constant term describing the baseline. A periodicity condition is also imposed to the function, requiring f (0) = f (2π ). The integrals of the Gaussian terms, ANS and AAS, correspond to the associated-particle yields for the near (NS)- and away (AS)-side peaks, respectively, while σfit,NS and σfit,AS quantify the widths of the correlation peaks. By symmetry considerations, the mean of the Gaussian functions are fixed to ϕ = 0 and ϕ = π . The baseline b represents the physical minimum of the ϕ distribution. To limit the effect of statistical fluctuations on the estimate of the associated yields, b is fixed to the weighted average of the points in the transverse region, defined as π/4 < | ϕ| < π/2, using the inverse of the square of the point statistical uncertainty as weights. Given the symmetry of the correlation distributions around ϕ = 0 and ϕ = π , the azimuthal distributions are reported in the range 0 < ϕ < π to reduce statistical fluctuations. The effect of a v2-like modulation in the ϕ distribution, which could be present in p–Pb collisions, was estimated and assessed in Sect. 5. In the case of the simulations, for which statistical fluctuations are negligible, the baseline is estimated as the minimum of the azimuthal-correlation distribution. An alternative fitting procedure based on a convolution of two Gaussian functions for the description of the NS peak was performed for Monte Carlo simulations. The resulting NS yields were found to be compatible with those obtained with the standard procedure, with a maximum variation of 7% (10%) in pp (p–Pb) collisions in case of EPOS 3 simulations. 4 Systematic uncertainties The fit of the D-meson invariant-mass distribution introduces systematic uncertainties on Speak and Bpeak (Sect. 3.2, Eq. 1). The uncertainty on the correlation distribution was estimated by calculating Bpeak from the integral of the background term of the invariant-mass fit function in the range |M − μ| < 2σ and by varying the fit procedure. In particular, the fit was repeated modeling the background distribution with a linear function and a parabola instead of an exponential function (for D0 and D+ mesons only), considering a different histogram binning, and varying the fit range. A 10% systematic uncertainty was estimated from the corresponding variation of the azimuthal-correlation distribution. No significant trend was observed as a function of ϕ and the same uncertainty was estimated for all D-meson species in all pT-intervals and in both pp and p–Pb collision systems. A 5% uncertainty (10% for D+ mesons in p–Pb collisions) arises from the possible dependence of the shape of C˜ ( ϕ, η)sidebands on the sideband range. This source of uncertainty was determined by restricting the invariant-mass sideband window to the intervals 4σ < |M − μ| < 6σ or to 6σ < |M − μ| < 8σ for all the D mesons, and also by considering, for D0 and D+ mesons, only the left or only the right sideband. The uncertainty on the correction for the associatedparticle reconstruction efficiency was assessed by varying the selection criteria applied to the reconstructed tracks, removing the request of at least three associated clusters in the ITS, or demanding a hit on at least one of the two SPD layers. A ±4% uncertainty was estimated for p–Pb collisions, while a +10% contribution was obtained for the pp analysis, with −5% the +10% contribution arising from the request of hits in the SPD. No significant trend in ϕ was observed. The uncertainty on the residual contamination from secondary tracks was evaluated by repeating the analysis varying the cut on the DCA in the (x , y) plane from 0.1 cm to 1 cm, and re-evaluating the purity of charged primary particles for each variation. This resulted in a 5% (3.5%) systematic uncertainty in pp (p–Pb) collisions, independent of ϕ and pTassoc. A 5% systematic effect originating from the correction of the D-meson reconstruction efficiency was evaluated by applying tighter and looser topological selections on the Dmeson candidates. No significant dependence on ϕ was observed and the same uncertainty was estimated for the three D-meson pT intervals, apart from D+ meson in p–Pb collisions, for which a 10% uncertainty was assigned. The uncertainty on the subtraction of the beauty feeddown contribution was quantified by generating the templates MC templ of feed-down azimuthal-correlation distributions, C˜ feed−down ( ϕ) in Eq. 2, with different PYTHIA 6 tunes (Perugia 0, Perugia 2010, see Sect. 2.2), and by considering the range of fprompt values obtained by varying the prompt and feed-down D-meson pT-differential production cross sections within FONLL uncertainty band, as described in [15]. The effect on the azimuthal-correlation distributions is ϕ dependent and contained within 8% and is more pronounced in the near side, in particular in the low and mid D-meson pT intervals. Table 1 List of systematic uncertainties for the ϕ-correlation distributions in pp and p–Pb collisions. See text for details +10%, −5% The consistency of the whole correction procedure, prior to the feed-down subtraction, was verified by performing the analysis on simulated events (“Monte-Carlo closure test”) separately for prompt and feed-down D mesons. For prompt D mesons, no effect was found for both pp and p–Pb collision systems. Conversely, for feed-down D mesons, an overestimate by about 20% in the near side was found for both collision systems. It was verified that the source of this excess is related to a bias induced by the topological selection applied to D mesons, that tends to favour cases with a small angular opening between the products of the beauty-hadron decay, thus between the D meson and the other decay particles. This effect results in a ϕ-dependent overestimate of the feed-down subtracted correlation distribution in the near side, contained within 2%. The systematic uncertainties affecting the ϕ-correlation distributions are summarised in Table 1 for both pp and p–Pb collision systems. The ϕ-dependent parts of the uncertainties arising from the feed-down subtraction and the MonteCarlo closure test define the ϕ-uncorrelated systematic uncertainties. All the other contributions, correlated in ϕ, act as a scale uncertainty. No significant dependence on the transverse momentum of D mesons and associated particles was observed for both ϕ-correlated and uncorrelated uncertainties, except for the feed-down systematic uncertainty. Different approaches were applied to estimate the systematic uncertainty on the near-side peak associated yield and peak width and on the baseline, obtained from the ANS, σfit,NS, and b parameters of the fit of the azimuthalcorrelation distribution, as described in Sect. 3.3. The main source of uncertainty originates from the definition of the baseline itself, which is connected to the assumption that the observed variation of the azimuthal-correlation distribution in the transverse region is determined mainly by statistical fluctuations rather than by the true physical trend. The variation of ANS, σfit,NS, and b values obtained when considering a ±π/4 variation of the ϕ range defining the transverse region is interpreted as the systematic uncertainty due to the baseline definition. In addition, the fits were repeated by moving upwards and downwards the data points by the corresponding value of the ϕ-uncorrelated systematic uncertainty. The final systematic uncertainty was calculated by summing in quadrature the aforementioned contributions and, for the associated yields and baseline, also the systematic uncertainty correlated in ϕ. The values of the total systematic uncertainties on the near-side peak yield, width, and baseline are reported in Table 2, for two intervals of transverse momentum of D mesons and associated particles. Considering all the measured kinematic ranges, the uncertainties vary from ±12 to ±25% for the near-side peak yield, from ±2 to ±13% for the near-side peak width and from ±11 to ±16% for the baseline. Typically, lower uncertainties are obtained for p–Pb collisions, where the larger available statistics of the correlation distributions allow for a more precise estimate of the baseline height, which constitutes the main source of uncertainty also on the evaluation of the near-side peak associated yield and width. 5 Results The azimuthal-correlation distributions of D0, D+, D∗+ mesons with charged particles with passoc > 1 GeV/c are T compared in Fig. 2 for 5 < pTD < 8 GeV/c in pp collisions Table 2 List of systematic uncertainties for near-side (NS) peak associated yield, near-side peak width, and baseline in pp and p–Pb collisions, for two different kinematic ranges of D mesons and associated particles. See text for details ALICE 1 ND2 Fig. 2 Comparison of the azimuthal-correlation distributions of D mesons with charged particles obtained for D0, D+ and D∗+ mesons for 5 < pTD < 8 GeV/c, pTassoc > 1 GeV/c in pp collisions at √s = 7 TeV (left panel) and for 8 < pTD < 16 GeV/c, passoc > 1 GeV/c in p–Pb collisions at √sNN = 5.02 TeV (right T (left panel) and for 8 < pTD < 16 GeV/c in p–Pb collisions (right panel). The distributions obtained with the three Dmeson species are compatible within the quadratic sum (wi , i = D0, D+, D∗+) of the statistical uncertainty and of the systematic uncertainties on the signal, background normalisation, and on the background shape (see Table 1), that are uncorrelated among the three meson species. The D0-, D+, D∗+-meson data are averaged using 1/wi2 as weights. The averages of the distributions are shown, for all the considered kinematic ranges, in Fig. 3 for pp and p–Pb collisions. A rising trend of the height of the near-side peak with increasing D-meson pT is observed for both collision systems. A similar trend is present for hadron–hadron correlations measured at Tevatron and LHC energies [68–71]: an increase of hadron multiplicity in jets with increasing jet energy is expected from the evolution of parton cascade with the parton energy for both light and heavy quarks [19]. A decrease of the baseline level with increasing pT of the associated particles can also be noticed. Figure 4 shows the ϕ distributions after the subtraction of the baseline, calculated as described in Sect. 3.3. The distributions show a near-side peak along with a wider and lower peak in the away-side region. The results obtained for the two collision systems are compatible within the total uncertainties. According to simulations of pp collisions performed using PYTHIA 6 (Perugia-0, -2010, and -2011 tunes), the different centre-of-mass energy and the slightly different D-meson rapidity range of the two measurements should induce variations in the baseline-subtracted azimuthal-correlation distributions smaller than 7% in the near- and away-side regions. The same estimate is obtained with POWHEG+PYTHIA simulations including the EPS09 parametrisation of nuclear PDFs (see Sect. 2.2). Such differences are well below the current level of uncertainties. panel). The statistical uncertainties are shown as error bars, the ϕuncorrelated systematic uncertainties as boxes, while the part of systematic uncertainty correlated in ϕ is reported as text (scale uncertainty). The latter is largely uncorrelated among the D-meson species A further comparison of the results from pp and p–Pb collisions is done by quantifying the integrals and the widths of the near-side correlation peaks by fitting the measured distributions as described in Sect. 3.3. The fit results are reported only for the near-side peak parameters and the baseline because of the poor statistical precision on the fit parameters of the awayside peaks. Figure 5 shows an exemplary fit to the azimuthalcorrelation distributions of D mesons with charged particles with passoc > 1 GeV/c, for 5 < pTD < 8 GeV/c in pp col T lisions (left panel) and for 8 < pD < 16 GeV/c in p–Pb T collisions (right panel). The curves superimposed to the data represent the three terms of the function defined in Eq. 3. Within the uncertainties, the fit function describes the measured distributions in all kinematic cases considered, yielding χ 2/NDF values close to unity. The evolution of the near-side peak associated yield as a function of the D-meson pT is reported in Fig. 6 (top row), for pp and p–Pb collisions, for passoc > 0.3 GeV/c (left panel) and T for the two sub-intervals 0.3 < passoc < 1 GeV/c (midT dle panel) and passoc > 1 GeV/c (right panel). The near T side peak associated yield exhibits an increasing trend with D-meson pT and has similar values, within uncertainties, for the softer (0.3 < passoc < 1 GeV/c) and the harder T ( pTassoc > 1 GeV/c) sub-ranges of passoc used, in each DT meson pT interval considered. The values obtained for pp and p–Pb collision data are compatible within statistical uncertainties. In the bottom row of the same figure the width of the near-side Gaussian term (σfit,NS) is shown. Although the case with passoc > 0.3 GeV/c seems to suggest that σfit,NS T does not strongly depend on D-meson pT in the range of the measurement, the current level of uncertainty does not allow for quantification of the dependence of σfit,NS on Dmeson and associated charged particle pT, as well as any potential difference between the values extracted using pp Average D0, D+, D*+ −+1103%% scale uncertainty (pp) −+1100%% scale uncertainty (p-Pb) −+1100%% scale uncertainty (p-Pb) −+1100%% scale uncertainty (p-Pb) −+1100%% scale uncertainty (p-Pb) 1 ND4 Fig. 3 Average of the azimuthal-correlation distributions of D0, D+ and D∗+ mesons with 3 < pTD < 5 GeV/c (left column), 5 < pD < 8 GeV/c (middle column), and 8 < pTD < 16 GeV/c (right T column), with charged particles with passoc > 0.3 GeV/c (top row), T 0.3 < passoc < 1 GeV/c (middle row), and passoc > 1 GeV/c (bottom T T row), measured in pp collisions at √s = 7 TeV and in p–Pb collisions at √sNN = 5.02 TeV. The statistical uncertainties are shown as error bars, the ϕ-uncorrelated systematic uncertainties as boxes, while the part of systematic uncertainty correlated in ϕ is reported as text (scale uncertainty) and p–Pb data. In particular, our approach for baseline calculation (Sect. 3.3) guarantees a robust estimate of the minimum, but the baseline uncertainty and its impact on the associated-yield uncertainty are rather large (Sect. 4). This systematic uncertainty is expected to be significantly reduced in future measurements with larger data samples, where a smaller ϕ range for the baseline calculation could be used. A v2-like modulation of the baseline would introduce a bias in the measurement of the associated yield and peak width and that needs to be taken into account while interpreting the measured quantities in terms of charm-jet properties. In order to get an estimate of this possible effect, for the p–Pb case the fit was repeated by subtracting from the correlation distribution a v2-like modulation assuming v2 = 0.05 for D mesons and v2 = 0.05 (0.1) for associated charged particles with pT > 0.3 (1) GeV/c. These values were chosen on the basis of charged-particle measurements in highmultiplicity p–Pb collisions [30] and assuming for D mesons the maximum value predicted in [2] for the 20% most central p–Pb collisions as a test case. With such assumptions, rather extreme also considering that this measurement is performed without any selection on event multiplicity, ANS varies by −10% (−6%) for D mesons with 5 < pT < 8 GeV/c and for 0.3 < passoc < 1 GeV/c ( pTassoc > 1 GeV/c). The vari T ations on σfit,NS and on the baseline are below 4 and 1%, respectively. Significantly smaller modifications result for D mesons with 8 < pT < 16 GeV/c. With the available statis 1 ND 0 Average D0, D+, D*+ −+1103%% scale uncertainty (pp) −+1100%% scale uncertainty (p-Pb) −+1103%% scale uncertainty (pp) −+1100%% scale uncertainty (p-Pb) −+1114%% scale uncertainty (pp) −+1100%% scale uncertainty (p-Pb) −+1100%% scale uncertainty (p-Pb) −+1100%% scale uncertainty (p-Pb) tics, the precision of the measurement is not sufficient to observe or exclude these modifications. Figure 7 shows the comparison of the averaged azimuthalcorrelation distributions measured in pp collisions with expectations from simulations performed with PYTHIA, POWHEG+PYTHIA, and EPOS 3 (see Sect. 2.2), after the baseline subtraction. The average of the two lowest values of the azimuthal-correlation distribution is used to define the uncertainty related to the baseline definition in Monte-Carlo simulations (see Sect. 3.3). This uncertainty is negligible and not displayed in the figures. The distributions obtained with the different generators and tunes do not show significant differences in the near side, except from EPOS 3 which tends to have higher and wider distributions. In the away side, the statistical uncertainties are shown as error bars, the ϕ-uncorrelated systematic uncertainties as boxes around the data points, the part of systematic uncertainty correlated in ϕ is reported as text (scale uncertainty), the uncertainties deriving from the subtraction of the baselines are represented by the boxes at ϕ > π PYTHIA 6 tunes Perugia 0 and Perugia 2010 tend to have higher correlation values, especially for passoc > 1 GeV/c, T compared to the other simulation results. Similar considerations hold for EPOS 3 in the case of D mesons with 8 < pT < 16 GeV/c. The considered Monte-Carlo simulations describe, within the uncertainties, the data in the whole ϕ range. The comparison of the associated yield in the near-side peak in data and in simulations is displayed in the top row of Figs. 8 and 9, for pp and p–Pb collisions, respectively. The simulations obtained with EPOS 3 provide a better description of the near-side yields for D mesons with 8 < pT < 16 GeV/c in both pp and p–Pb collisions. At lower D-meson pT a better agreement is obtained with PYTHIA and POWHEG+PYTHIA simulations. The width −+1103%% scale uncertainty −+1100%% scale uncertainty systematic uncertainties as boxes, while the part of systematic uncertainty correlated in ϕ is reported as text (scale uncertainty). The terms of the fit function described in Sect. 3.3 are also shown separately: nearside Gaussian function (blue dashed line), away-side Gaussian function (green dashed–dotted line) and baseline constant term (magenta dotted line) <7% variation expected from energy and rapidity difference (PYTHIA6, Perugia 2011) )d 0.4 a r ( ,SN0.3 ift σ p-Pb points and error boxes shifted by ΔpT = +0.3 GeV/c Fig. 6 Comparison of the near-side peak associated yield (top row) and peak width (bottom row) in pp and p–Pb collisions as a function of pD, for passoc > 0.3 GeV/c (left column), 0.3 < passoc < 1 GeV/c T T T (middle column), and passoc > 1 GeV/c (right column). The points and T error boxes for p–Pb collisions are shifted by pT = +0.3 GeV/c. Statistical and systematic uncertainties are shown as error bars and boxes, respectively of the near-side peaks, shown in the second row of the same figures, is better reproduced by the simulations in the case of p–Pb than of pp results. The evolution of the baseline value as a function of the D-meson pT is compared for pp-collision data to expectations from PYTHIA simulations in the bottom row of Fig. 8 for the three ranges of passoc considered in the T Average D0,D+,D*+ ALICE baseline-subtraction uncertainty −+1103%% scale uncertainty −+1103%% scale uncertainty −+1114%% scale uncertainty −+1103%% scale uncertainty −+1103%% scale uncertainty −+1114%% scale uncertainty −+1103%% scale uncertainty −+1103%% scale uncertainty −+1114%% scale uncertainty Fig. 7 Comparison of ϕ-correlation distributions of D mesons with charged particles measured in pp collisions at √s = 7 TeV and MonteCarlo simulations performed with different event generators, after the analysis. The value of the baseline, mainly determined by the event multiplicity, does not show substantial variations as a function of D-meson pT, as expected also from PYTHIA and EPOS 3 simulations, which reproduce the observed values within the uncertainties. 6 Summary The first measurements of the azimuthal correlations between D mesons with charged particles in pp and p–Pb collisions at √s = 7 TeV and √sNN = 5.02 TeV, respectively, performed with the ALICE apparatus at the LHC were presented. The ϕ distributions were studied in pp collisions in three different D-meson transverse-momentum intervals, 3 < pTD < subtraction of the baseline. The statistical and systematic uncertainties of the measured distributions are displayed as in Fig. 4 5 GeV/c, 5 < pTD < 8 GeV/c, and 8 < pTD < 16 GeV/c, for associated charged particles with passoc > 0.3 GeV/c, T and in the two sub-ranges 0.3 < passoc < 1 GeV/c and T passoc > 1 GeV/c. For p–Pb collisions, the results were T reported in two D-meson pT ranges, 5 < pTD < 8 GeV/c, and 8 < pD < 16 GeV/c. The baseline-subtracted azimuthal T correlation distributions observed in the two collision systems are compatible within uncertainties. The variations expected from the lower nucleon-nucleon centre-of-mass energy of p–Pb collisions and from the slightly different Dmeson rapidity ranges used for the p–Pb analysis were studied with simulated pp collisions at the two centre-of-mass energies and are well below the sensitivity of the measurements. )d 0.4 a r ( i,tfSN 0.3 σ 0.2 pp, s = 7 TeV Fig. 8 Comparison of near-side peak associated yield (top row), nearside peak width (middle row), and baseline (bottom row) values measured in pp collisions at √s = 7 TeV with the expectations from simulations performed with different Monte-Carlo event generators. Statistical and systematic uncertainties are shown as error bars and boxes, respectively The properties of the near-side correlation peak, sensitive to the characteristics of the jet containing the D meson, were described in terms of the yield of associated charged particles and peak width, obtained by fitting the ϕ distributions with a function composed of a constant term, representing the physical minimum of the distribution, and two Gaussian terms modeling the near- and away-side peaks. The values measured in the two collision systems are compatible within uncertainties. The measured azimuthal distributions, as well as the properties of the correlation peaks, were compared to expectations from simulations performed with different Monte-Carlo generators. The simulations reproduce the correlation distributions within uncertainties. Considering that the overall uncertainty is dominated by the statistical component, the data collected from pp collisions at √s = 13 TeV in the ongoing Run 2 at the LHC will allow for a more precise measurement. In particular, the predicted increase of the cross section for charm production by more than a factor 2 at pT = 10 GeV/c at the higher collision energy [12], along with the foreseen larger integrated luminosity, will allow for a significant reduction of the statistical uncertainty, providing a more quantitative and constraining comparison of the data with expectations from Monte-Carlo generators. As mentioned in Sect. 5, with larger data samples a different determination of the baseline of the azimuthal-correlation distribution will become possible, bringing to a significant reduction of the system)d 0.4 a r ( ,SN0.3 itf σ p-Pb, sNN = 5.02 TeV Fig. 9 Comparison of near-side peak associated yield (top row) and near-side peak width (bottom row) values measured in p–Pb collisions at √sNN = 5.02 TeV with the expectations from simulations performed atic uncertainty on the measurement of the associated yields. The data that will be collected in next p–Pb collision runs at the LHC may also allow for a study of the evolution of the azimuthal-correlation distribution as a function of the event multiplicity, searching for possible long-range ridgelike structures already observed with angular correlation of light particles. The results reported in this paper represent a first step towards the measurement of possible modifications concerning the azimuthal correlation of D mesons with charged particles in Pb–Pb collisions, which has the potential to provide important information on the charm-quark energy-loss mechanisms in the presence of the medium formed in heavyion collisions at LHC energies. Given the same collision energy, the p–Pb results presented in this paper could serve as a reference to study medium effects in Pb–Pb collisions at √sNN = 5.02 TeV collected during the LHC Run 2. Acknowledgements 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: State Committee of Science, World Federation of Sciwith different Monte-Carlo event generators. Statistical and systematic uncertainties are shown as error bars and boxes, respectively entists (WFS) and Swiss Fonds Kidagan, Armenia; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Financiadora de Estudos e Projetos (FINEP), Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP); Ministry of Science & Technology of China (MSTC), National Natural Science Foundation of China (NSFC) and Ministry of Education of China (MOEC)”; Ministry of Science, Education and Sports of Croatia and Unity through Knowledge Fund, Croatia; Ministry of Education and Youth of the Czech Republic; Danish Natural Science Research Council, the Carlsberg Foundation and the Danish National Research Foundation; The European Research Council under the European Community’s Seventh Framework Programme; Helsinki Institute of Physics and the Academy of Finland; French CNRS-IN2P3, the ‘Region Pays de Loire’, ‘Region Alsace’, ‘Region Auvergne’ and CEA, France; German Bundesministerium fur Bildung, Wissenschaft, Forschung und Technologie (BMBF) and the Helmholtz Association; General Secretariat for Research and Technology, Ministry of Development, Greece; National Research, Development and Innovation Office (NKFIH), Hungary; Council of Scientific and Industrial Research (CSIR), New Delhi; Department of Atomic Energy and Department of Science and Technology of the Government of India; Istituto Nazionale di Fisica Nucleare (INFN) and Centro Fermi – Museo Storico della Fisica e Centro Studi e Ricerche “Enrico Fermi”, Italy; Japan Society for the Promotion of Science (JSPS) KAKENHI and MEXT, Japan; National Research Foundation of Korea (NRF); Consejo Nacional de Cienca y Tecnologia (CONACYT), Direccion General de Asuntos del Personal Academico(DGAPA), México, Amerique Latine Formation academique – European Commission (ALFA-EC) and the EPLANET Program (European Particle Physics Latin American Network); Stichting voor Fundamenteel Onderzoek der Materie (FOM) and the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO), Netherlands; Research Council of Norway (NFR); Pontificia Universidad Católica del Perú; National Science Centre, Poland; Ministry of National Education/Institute for Atomic Physics and National Council of Scientific Research in Higher Education (CNCSI-UEFISCDI), Romania; Joint Institute for Nuclear Research, Dubna; Ministry of Education and Science of Russian Federation, Russian Academy of Sciences, Russian Federal Agency of Atomic Energy, Russian Federal Agency for Science and Innovations and The Russian Foundation for Basic Research; Ministry of Education of Slovakia; Department of Science and Technology, South Africa; Centro de Investigaciones Energeticas, Medioambientales y Tecnologicas (CIEMAT), E-Infrastructure shared between Europe and Latin America (EELA), Ministerio de Economía y Competitividad (MINECO) of Spain, Xunta de Galicia (Consellería de Educación), Centro de Aplicaciones Tecnolgicas y Desarrollo Nuclear (CEADEN), Cubaenergía, Cuba, and IAEA (International Atomic Energy Agency); Swedish Research Council (VR) and Knut & Alice Wallenberg Foundation (KAW); National Science and Technology Development Agency (NSDTA), Suranaree University of Technology (SUT) and Office of the Higher Education Commission under NRU project of Thailand; Ukraine Ministry of Education and Science; United Kingdom Science and Technology Facilities Council (STFC); The United States Department of Energy, the United States National Science Foundation, the State of Texas, and the State of Ohio. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3. ALICE Collaboration 4 Department of Physics, Centre for Astroparticle Physics and Space Science (CAPSS), Bose Institute, Kolkata, India 5 Budker Institute for Nuclear Physics, Novosibirsk, Russia 6 California Polytechnic State University, San Luis Obispo, CA, USA 7 Central China Normal University, Wuhan, China 8 Centre de Calcul de l’IN2P3, Villeurbanne, Lyon, France 9 Centro de Aplicaciones Tecnológicas y Desarrollo Nuclear (CEADEN), Havana, Cuba 10 Centro de Investigaciones Energéticas Medioambientales y Tecnológicas (CIEMAT), Madrid, Spain 11 Centro de Investigación y de Estudios Avanzados (CINVESTAV), Mexico City, Mérida, Mexico 12 Centro Fermi-Museo Storico della Fisica e Centro Studi e Ricerche “Enrico Fermi’, Rome, Italy 13 Chicago State University, Chicago, IL, USA 14 China Institute of Atomic Energy, Beijing, China 15 Commissariat à l’Energie Atomique, IRFU, Saclay, France 16 COMSATS Institute of Information Technology (CIIT), Islamabad, Pakistan 17 Departamento de Física de Partículas and IGFAE, Universidad de Santiago de Compostela, Santiago de Compostela, Spain 18 Department of Physics, Aligarh Muslim University, Aligarh, India 19 Department of Physics, Ohio State University, Columbus, OH, USA 20 Department of Physics, Sejong University, Seoul, South Korea 21 Department of Physics, University of Oslo, Oslo, Norway 22 Department of Physics and Technology, University of Bergen, Bergen, Norway 23 Dipartimento di Fisica dell’Università ‘La Sapienza’ and Sezione INFN, Rome, Italy 24 Dipartimento di Fisica dell’Università and Sezione INFN, Cagliari, Italy 25 Dipartimento di Fisica dell’Università and Sezione INFN, Trieste, Italy 26 Dipartimento di Fisica dell’Università and Sezione INFN, Turin, Italy 27 Dipartimento di Fisica e Astronomia dell’Università and Sezione INFN, Bologna, Italy 28 Dipartimento di Fisica e Astronomia dell’Università and Sezione INFN, Catania, Italy 29 Dipartimento di Fisica e Astronomia dell’Università and Sezione INFN, Padua, Italy 30 Dipartimento di Fisica ‘E.R. Caianiello’ dell’Università and Gruppo Collegato INFN, Salerno, Italy 31 Dipartimento DISAT del Politecnico and Sezione INFN, Turin, Italy 32 Dipartimento di Scienze e Innovazione Tecnologica dell’Università del Piemonte Orientale and INFN Sezione di Torino, Alessandria, Italy 33 Dipartimento Interateneo di Fisica ‘M. Merlin’ and Sezione INFN, Bari, Italy 34 Division of Experimental High Energy Physics, University of Lund, Lund, Sweden 35 European Organization for Nuclear Research (CERN), Geneva, Switzerland 36 Excellence Cluster Universe, Technische Universität München, Munich, Germany 37 Faculty of Engineering, Bergen University College, Bergen, Norway 38 Faculty of Mathematics, Physics and Informatics, Comenius University, Bratislava, Slovakia 39 Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, Prague, Czech Republic 40 Faculty of Science, P.J. Šafárik University, Kosice, Slovakia 41 Faculty of Technology, Buskerud and Vestfold University College, Tonsberg, Norway 42 Frankfurt Institute for Advanced Studies, Johann Wolfgang Goethe-Universität Frankfurt, Frankfurt, Germany 43 Gangneung-Wonju National University, Gangneung, South Korea 44 Department of Physics, Gauhati University, Guwahati, India 45 Helmholtz-Institut für Strahlen- und Kernphysik, Rheinische Friedrich-Wilhelms-Universität Bonn, Bonn, Germany 46 Helsinki Institute of Physics (HIP), Helsinki, Finland 47 Hiroshima University, Hiroshima, Japan 48 Indian Institute of Technology Bombay (IIT), Mumbai, India 49 Indian Institute of Technology Indore, Indore, India 50 Indonesian Institute of Sciences, Jakarta, Indonesia 51 Inha University, Incheon, South Korea 52 Institut de Physique Nucléaire d’Orsay (IPNO), Université Paris-Sud, CNRS-IN2P3, Orsay, France 53 Institute for Nuclear Research, Academy of Sciences, Moscow, Russia 54 Institute for Subatomic Physics of Utrecht University, Utrecht, The Netherlands 55 Institute for Theoretical and Experimental Physics, Moscow, Russia 56 Institute of Experimental Physics, Slovak Academy of Sciences, Kosice, Slovakia 57 Institute of Physics, Academy of Sciences of the Czech Republic, Prague, Czech Republic 58 Institute of Physics, Bhubaneswar, India 59 Institute of Space Science (ISS), Bucharest, Romania 60 Institut für Informatik, Johann Wolfgang Goethe-Universität Frankfurt, Frankfurt, Germany 61 Institut für Kernphysik, Johann Wolfgang Goethe-Universität Frankfurt, Frankfurt, Germany 62 Institut für Kernphysik, Westfälische Wilhelms-Universität Münster, Münster, Germany 63 Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Mexico City, Mexico 64 Instituto de Física, Universidad Nacional Autónoma de México, Mexico City, Mexico 65 Institut Pluridisciplinaire Hubert Curien (IPHC), Université de Strasbourg, CNRS-IN2P3, Strasbourg, France 66 iThemba LABS, National Research Foundation, Somerset West, South Africa 67 Joint Institute for Nuclear Research (JINR), Dubna, Russia 68 Konkuk University, Seoul, South Korea 69 Korea Institute of Science and Technology Information, Daejeon, South Korea 70 KTO Karatay University, Konya, Turkey 71 Laboratoire de Physique Corpusculaire (LPC), Clermont Université, Université Blaise Pascal, CNRS-IN2P3, Clermont-Ferrand, France 72 Laboratoire de Physique Subatomique et de Cosmologie, Université Grenoble-Alpes, CNRS-IN2P3, Grenoble, France 73 Laboratori Nazionali di Frascati, INFN, Frascati, Italy 74 Laboratori Nazionali di Legnaro, INFN, Legnaro, Italy 75 Lawrence Berkeley National Laboratory, Berkeley, CA, USA 76 Moscow Engineering Physics Institute, Moscow, Russia 77 Nagasaki Institute of Applied Science, Nagasaki, Japan 78 National Centre for Nuclear Studies, Warsaw, Poland 79 National Institute for Physics and Nuclear Engineering, Bucharest, Romania 80 National Institute of Science Education and Research, Bhubaneswar, India 81 National Research Centre Kurchatov Institute, Moscow, Russia 82 Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark 83 Nikhef, Nationaal instituut voor subatomaire fysica, Amsterdam, The Netherlands 84 Nuclear Physics Group, STFC Daresbury Laboratory, Daresbury, UK 85 Nuclear Physics Institute, Academy of Sciences of the Czech Republic, Rˇež u Prahy, Czech Republic 86 Oak Ridge National Laboratory, Oak Ridge, TN, USA 87 Petersburg Nuclear Physics Institute, Gatchina, Russia 88 Physics Department, Creighton University, Omaha, NE, USA 89 Physics Department, Panjab University, Chandigarh, India 90 Physics Department, University of Athens, Athens, Greece 91 Physics Department, University of Cape Town, Cape Town, South Africa 92 Physics Department, University of Jammu, Jammu, India 93 Physics Department, University of Rajasthan, Jaipur, India 94 Physikalisches Institut, Universität München, Munich, Germany 95 Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany 96 Physik Department, Technische Universität München, Munich, Germany 97 Purdue University, West Lafayette, IN, USA 98 Pusan National University, Pusan, South Korea 99 Research Division and ExtreMe Matter Institute EMMI, GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany 100 Rudjer Boškovic´ Institute, Zagreb, Croatia 101 Russian Federal Nuclear Center (VNIIEF), Sarov, Russia 102 Saha Institute of Nuclear Physics, Kolkata, India 103 School of Physics and Astronomy, University of Birmingham, Birmingham, UK 104 Sección Física, Departamento de Ciencias, Pontificia Universidad Católica del Perú, Lima, Peru 105 Sezione INFN, Bari, Italy 106 Sezione INFN, Bologna, Italy 107 Sezione INFN, Cagliari, Italy 108 Sezione INFN, Catania, Italy 109 Sezione INFN, Padua, Italy 110 Sezione INFN, Rome, Italy 111 Sezione INFN, Trieste, Italy 112 Sezione INFN, Turin, Italy 113 SSC IHEP of NRC Kurchatov institute, Protvino, Russia 114 Stefan Meyer Institut für Subatomare Physik (SMI), Vienna, Austria 115 SUBATECH, Ecole des Mines de Nantes, Université de Nantes, CNRS-IN2P3, Nantes, France 116 Suranaree University of Technology, Nakhon Ratchasima, Thailand 117 Technical University of Košice, Kosice, Slovakia 118 Technical University of Split FESB, Split, Croatia 119 The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Kraków, Poland 120 Physics Department, The University of Texas at Austin, Austin, TX, USA 121 Universidad Autónoma de Sinaloa, Culiacán, Mexico 122 Universidade de São Paulo (USP), São Paulo, Brazil 123 Universidade Estadual de Campinas (UNICAMP), Campinas, Brazil 124 Universidade Federal do ABC, Santo Andre, Brazil 125 University of Houston, Houston, TX, USA 126 University of Jyväskylä, Jyväskylä, Finland 127 University of Liverpool, Liverpool, UK 128 University of Tennessee, Knoxville, TN, USA 129 University of the Witwatersrand, Johannesburg, South Africa 130 University of Tokyo, Tokyo, Japan 131 University of Tsukuba, Tsukuba, Japan 132 University of Zagreb, Zagreb, Croatia 133 Université de Lyon, Université Lyon 1, CNRS/IN2P3, IPN-Lyon, Villeurbanne, Lyon, France 134 Università di Brescia, Brescia, Italy 135 V. Fock Institute for Physics, St. Petersburg State University, St. Petersburg, Russia 136 Variable Energy Cyclotron Centre, Kolkata, India 137 Warsaw University of Technology, Warsaw, Poland 138 Wayne State University, Detroit, MI, USA 139 Wigner Research Centre for Physics, Hungarian Academy of Sciences, Budapest, Hungary 140 Yale University, New Haven, CT, USA 141 Yonsei University, Seoul, South Korea 142 Zentrum für Technologietransfer und Telekommunikation (ZTT), Fachhochschule Worms, Worms, Germany a Deceased b Also at: Georgia State University, Atlanta, Georgia, USA c Also at: Also at Department of Applied Physics, Aligarh Muslim University, Aligarh, India d Also at: M.V. Lomonosov Moscow State University, D.V. Skobeltsyn Institute of Nuclear, Physics, Moscow, Russia 1. A. Beraudo , A. De Pace , M. Monteno , M. Nardi , F. Prino , Heavy flavors in heavy-ion collisions: quenching, flow and correlations . Eur. Phys. J. C 75 ( 3 ), 121 ( 2015 ). doi:10.1140/epjc/ s10052- 015 - 3336 -6. arXiv: 1410 .6082 [hep-ph] 2. A. Beraudo , A. De Pace , M. Monteno , M. Nardi , F. Prino , Heavy-flavour production in high-energy d-Au and p-Pb collisions . JHEP 1603 , 123 ( 2016 ). doi:10.1007/ JHEP03(2016)123 . arXiv: 1512 .05186 [hep-ph] 3. M.L. Mangano , P. Nason , G. Ridolfi , Heavy quark correlations in hadron collisions at next-to-leading order . Nucl. Phys B 373 , 295 - 345 ( 1992 ). doi:10.1016/ 0550 -3213(92) 90435 - E 4. E. Norrbin , T. Sjöstrand , Production and hadronization of heavy quarks . Eur. Phys. J C 17 , 137 - 161 ( 2000 ). doi:10.1007/ s100520000460. arXiv:hep-ph/ 0005110 [hep-ph] 5. STAR Collaboration , C. Adler et al., Disappearance of back-toback high pT hadron correlations in central Au + Au collisions at √ sN N = 200 GeV Phys. Rev. Lett . 90 , 082302 ( 2003 ). doi:10. 1103/PhysRevLett.90.082302. arXiv:nucl-ex/0210033 6. STAR Collaboration , J. Adams et al., Evidence from d + Au measurements for final state suppression of high pT hadrons in Au + Au collisions at RHIC . Phys. Rev. Lett . 91 , 072304 ( 2003 ). doi:10. 1103/PhysRevLett.91.072304. arXiv:nucl-ex/0306024 7. PHENIX Collaboration, A. Adare et al., Trends in yield and azimuthal shape modification in dihadron correlations in relativistic heavy ion collisions . Phys. Rev. Lett . 104 ( 2010 ) 252301,doi:10. 1103/PhysRevLett.104.252301. arXiv: 1002 .1077 [nucl-ex] 8. ALICE Collaboration , K. Aamodt et al., Particle-yield modification in jet-like azimuthal di-hadron correlations in Pb-Pb collisions at √ sN N = 2.76 TeV. Phys. Rev. Lett . 108 , 092301 ( 2012 ). doi:10. 1103/PhysRevLett.108.092301. arXiv: 1110 .0121 [nucl-ex] 9. H. Zhang , J.F. Owens , E. Wang , X.-N. Wang , Dihadron tomography of high-energy nuclear collisions in NLO pQCD . Phys. Rev. Lett . 98 , 212301 ( 2007 ). doi:10.1103/PhysRevLett.98.212301. arXiv:nucl-th/0701045 10. T. Renk , The phenomenology of elastic energy loss . Phys. Rev C . 76 , 064905 ( 2007 ). doi:10.1103/PhysRevC.76.064905. arXiv: 0708 .4319 [hep-ph] 11. PHENIX Collaboration, A. Adare et al., Azimuthal correlations of electrons from heavy-flavor decay with hadrons in p+ p and Au + Au collisions at √sN N = 200 GeV . Phys. Rev C . 83 , 044912 ( 2011 ). doi:10.1103/PhysRevC.83.044912. arXiv: 1011 .1477 [nucl-ex] 12. M. Cacciari , M. Greco , P. Nason , The pT spectrum in heavy flavor hadroproduction . JHEP 05 , 007 ( 1998 ). doi:10.1088/ 1126 - 6708 / 1998/05/007. arXiv:hep-ph/ 9803400 13. B.A. Kniehl , G. Kramer , I. Schienbein , H. Spiesberger , Inclusive charmed-meson production at the CERN LHC . Eur . Phys. J C 72 , 2082 ( 2012 ). doi:10.1140/epjc/s10052- 012 - 2082 -2. arXiv: 1202 .0439 [hep-ph] 14. R. Maciula , A. Szczurek , Open charm production at the LHC-kt - factorization approach . Phys. Rev. D 87 , 094022 ( 2013 ). doi:10. 1103/PhysRevD.87.094022. arXiv: 1301 .3033 [hep-ph] 15. ALICE Collaboration, B. Abelev et al., Measurement of charm production at central rapidity in proton-proton collisions at √s = 7 TeV . JHEP 01 , 128 ( 2012 ). doi:10.1007/ JHEP01(2012)128. arXiv:1111.1553 [hep-ex] 16. ALICE Collaboration, B. Abelev et al., Measurement of charm production at central rapidity in proton-proton collisions at √ s = 2.76 TeV. JHEP 07 , 191 ( 2012 ). doi:10.1007/ JHEP07(2012)191. arXiv:1205.4007 [hep-ex] 17. T. Sjostrand , S. Mrenna , P.Z. Skands , PYTHIA 6.4 physics and manual.' JHEP 05 , 026 ( 2006 ). doi:10.1088/ 1126 - 6708 /2006/05/ 026. arXiv:hep-ph/ 0603175 18. G. Corcella , I.G. Knowles , G. Marchesini , S. Moretti , K. Odagiri , P. Richardson , M.H. Seymour , B.R. Webber , HERWIG 6: an event generator for hadron emission reactions with interfering gluons (including supersymmetric processes) . JHEP 0101 , 010 ( 2001 ) 19. Y.L. Dokshitzer , V.A. Khoze , S.I. Troian , Particle spectra in light and heavy quark jets . J. Phys. G17 , 1481 - 1492 ( 1991 ). doi:10. 1088/ 0954 - 3899 /17/10/003 20. Y.L. Dokshitzer , F. Fabbri , V.A. Khoze , W. Ochs , Multiplicity difference between heavy and light quark jets revisited . Eur. Phys. J. C 45 , 387 - 400 ( 2006 ). doi:10.1140/epjc/s2005- 02424 -5. arXiv:hep-ph/ 0508074 21. R. Perez Ramos , V. Mathieu , M.-A. Sanchis-Lozano , Heavy quark flavour dependence of multiparticle production in QCD jets . JHEP 08 , 047 ( 2010 ). doi:10.1007/ JHEP08(2010)047 . arXiv: 1005 .1582 [hep-ph] 22. LHCb Collaboration , R. Aaij et al., Observation of double charm production involving open charm in pp collisions at √s = 7 TeV . JHEP 06 , 141 ( 2012 ). doi:10.1007/ JHEP06(2012)141. arXiv:1205.0975 [hep-ex] . [Addendum: JHEP 03 ( 2014 ) 108] 23. ATLAS Collaboration, G. Aad et al., Measurement of D∗+ /− meson production in jets from pp collisions at √s = 7 TeV with the ATLAS detector . Phys. Rev. D 85 , 052005 ( 2012 ). doi:10.1103/ PhysRevD.85.052005. arXiv: 1112 .4432 [hep-ex] 24. STAR Collaboration, M.M. Aggarwal et al., Measurement of the bottom contribution to non-photonic electron production in p + p collisions at √s = 200 GeV. Phys. Rev. Lett . 105 , 202301 ( 2010 ). doi:10.1103/PhysRevLett.105.202301. arXiv: 1007 .1200 [nucl-ex] 25. ALICE Collaboration, B. Abelev et al., Beauty production in pp collisions at √s = 2.76 TeV measured via semi-electronic decays . Phys. Lett B . 738 , 97 - 108 ( 2014 ). doi:10.1016/j.physletb. 2014 .09. 026. arXiv: 1405 .4144 [nucl-ex] 26. ALICE Collaboration, K. Aamodt et al., Elliptic flow of charged particles in Pb-Pb collisions at √sNN = 2.76 TeV . Phys. Rev. Lett . 105 , 252302 ( 2010 ). doi:10.1103/PhysRevLett.105.252302. arXiv: 1011 .3914 [nucl-ex] 27. ALICE Collaboration, K. Aamodt et al., Higher harmonic anisotropic flow measurements of charged particles in Pb-Pb collisions at √ sN N = 2.76 TeV. Phys. Rev. Lett . 107 , 032301 ( 2011 ). doi:10.1103/PhysRevLett.107.032301. arXiv: 1105 .3865 [nucl-ex] 28. CMS Collaboration, V. Khachatryan et al., Observation of longrange near-side angular correlations in proton-proton collisions at the LHC . JHEP 09 , 091 ( 2010 ). doi:10.1007/ JHEP09(2010)091. arXiv:1009.4122 [hep-ex] 29. ATLAS Collaboration, G. Aad et al., Observation of long-range elliptic anisotropies in √s =13 and 2.76 TeV pp collisions with the ATLAS detector . Phys. Rev. Lett . 116 , 172301 ( 2016 ). doi:10. 1103/PhysRevLett.116.172301. arXiv: 1509 .04776 [hep-ex] 30. ALICE Collaboration, B. Abelev et al., Long-range angular correlations on the near and away side in p-Pb collisions at √sNN = 5.02 TeV . Phys. Lett . B 719 , 29 - 41 ( 2013 ). doi:10.1016/j.physletb. 2013 . 01.012. arXiv: 1212 . 2001 [nucl-ex] 31. ALICE Collaboration, J. Adam et al., Forward-central two-particle correlations in p-Pb collisions at √sNN = 5.02 TeV . Phys. Lett . B 753 , 126 - 139 ( 2016 ). doi:10.1016/j.physletb. 2015 .12.010. arXiv: 1506 .08032 [nucl-ex] 32. ATLAS Collaboration, G. Aad et al., Observation of associated near-side and away-side long-range correlations in √sN N = 5.02 TeV proton-lead collisions with the ATLAS detector . Phys. Rev. Lett . 110 , 182302 ( 2013 ). doi:10.1103/PhysRevLett.110.182302. arXiv: 1212 .5198 [hep-ex] 33. CMS Collaboration, V. Khachatryan et al., Evidence for collective multiparticle correlations in p-Pb collisions . Phys. Rev. Lett . 115 ( 1 ), 012301 ( 2015 ). doi:10.1103/PhysRevLett.115.012301. arXiv: 1502 .05382 [nucl-ex] 34. PHENIX Collaboration, A. Adare et al., Quadrupole anisotropy in dihadron azimuthal correlations in central d+Au collisions at √sN N = 200 GeV . Phys. Rev. Lett . 111 (21), 212301 ( 2013 ). doi:10. 1103/PhysRevLett.111.212301. arXiv: 1303 .1794 [nucl-ex] 35. STAR Collaboration, L. Adamczyk et al., Long-range pseudorapidity dihadron correlations in d+Au collisions at √sNN = 200 GeV . Phys. Lett . B 747 , 265 - 271 ( 2015 ). doi:10.1016/j.physletb. 2015 .05.075. arXiv: 1502 .07652 [nucl-ex] 36. STAR Collaboration, L. Adamczyk et al., Effect of event selection on jetlike correlation measurement in d+Au collisions at √sNN = 200 GeV . Phys. Lett . B 743 , 333 - 339 ( 2015 ). doi:10. 1016/j.physletb. 2015 .02.068. arXiv: 1412 .8437 [nucl-ex] 37. K. Werner , I. Karpenko , T. Pierog , Ridge in proton-proton scattering at 7 TeV . Phys. Rev. Lett . 106 , 122004 ( 2011 ). doi:10.1103/ PhysRevLett.106.122004. arXiv: 1011 .0375 [hep-ph] 38. S. Alderweireldt , P. Van Mechelen , Obtaining the CMS ridge effect with multiparton interactions . In Proceedings, 3rd International Workshop on Multiple Partonic Interactions at the LHC (MPI@LHC 2011 ), pp. 33 - 40 ( 2012 ). arXiv:1203.2048 [hep-ph] . http://inspirehep.net/record/1093441/files/arXiv:1203.2048.pdf 39. P. Bozek , W. Broniowski , Collective dynamics in high-energy proton-nucleus collisions . Phys. Rev. C 88 , 014903 ( 2013 ). doi:10. 1103/PhysRevC.88.014903. arXiv: 1304 .3044 [nucl-th] 40. P. Bozek , W. Broniowski , Correlations from hydrodynamic flow in p-Pb collisions . Phys. Lett. B 718 , 1557 - 1561 ( 2013 ). doi:10. 1016/j.physletb. 2012 .12.051. arXiv: 1211 .0845 [nucl-th] 41. L. He , T. Edmonds , Z.-W. Lin , F. Liu , D. Molnar , F. Wang , Anisotropic parton escape is the dominant source of azimuthal anisotropy in transport models . Phys. Lett. B 753 , 506 - 510 ( 2016 ). doi:10.1016/j.physletb. 2015 .12.051. arXiv: 1502 .05572 [nucl-th] 42. K. Dusling , R. Venugopalan , Comparison of the color glass condensate to dihadron correlations in proton-proton and protonnucleus collisions . Phys. Rev. D 87 , 094034 ( 2013 ). doi:10.1103/ PhysRevD.87.094034. arXiv: 1302 .7018 [hep-ph] 43. B.A. Arbuzov , E.E. Boos , V.I. Savrin , CMS ridge effect at LHC as a manifestation of bremsstrahlung of gluons due to the quarkantiquark string formation . Eur. Phys. J. C 71 , 1730 ( 2011 ). doi:10. 1140/epjc/s10052- 011 - 1730 -2. arXiv: 1104 .1283 [hep-ph] 44. ALICE Collaboration, B. Abelev et al., Measurement of prompt Dmeson production in p − Pb collisions at √ sN N = 5.02 TeV. Phys. Rev. Lett . 113 (23), 232301 ( 2014 ). doi:10.1103/PhysRevLett.113. 232301. arXiv: 1405 .3452 [nucl-ex] 45. PHENIX Collaboration, A. Adare et al., Heavy-flavor electronmuon correlations in p + p and d+Au collisions at √sN N = 200 GeV . Phys. Rev. C 89 , 034915 ( 2014 ). doi:10.1103/PhysRevC.89. 034915. arXiv: 1311 .1427 [nucl-ex] 46. H. Fujii , K. Watanabe , Heavy quark pair production in high energy pA collisions: open heavy flavors . Nucl. Phys. A 920 , 78 - 93 ( 2013 ). doi:10.1016/j.nuclphysa. 2013 .10.006. arXiv: 1308 .1258 [hep-ph] 47. Y. Xu , S. Cao , G.-Y. Qin , W. Ke , M. Nahrgang , J. Auvinen , S.A. Bass , Heavy-flavor dynamics in relativistic p-Pb collisions at √sNN = 5.02 TeV . arXiv: 1510 .07520 [nucl-th] 48. PHENIX Collaboration, A. Adare et al., Cold-nuclear-matter effects on heavy-quark production in d+Au collisions at √sNN = 200 GeV . Phys. Rev. Lett . 109 , 242301 ( 2012 ). doi:10.1103/ PhysRevLett.109.242301. arXiv: 1208 .1293 [nucl-ex] 49. P. Bozek , A. Bzdak , G.-L. Ma , Rapidity dependence of elliptic and triangular flow in proton-nucleus collisions from collective dynamics . Phys. Lett. B 748 , 301 - 305 ( 2015 ). doi:10.1016/j.physletb. 2015 .06.007. arXiv: 1503 .03655 [hep-ph] 50. ALICE Collaboration, K. Aamodt et al., The ALICE experiment at the CERN LHC . JINST 3, S08002 ( 2008 ). doi:10.1088/ 1748 - 0221 /3/08/S08002 51. ALICE Collaboration, B. Abelev et al., Performance of the ALICE experiment at the CERN LHC. Int. J. Mod. Phys A . 29 , 1430044 ( 2014 ). doi:10.1142/S0217751X14300440. arXiv: 1402 .4476 [nucl-ex] 52. P.Z. Skands , Tuning Monte Carlo generators: the Perugia tunes . Phys. Rev. D 82 , 074018 ( 2010 ). doi:10.1103/PhysRevD.82. 074018. arXiv: 1005 .3457 [hep-ph] 53. X.-N. Wang , M. Gyulassy , HIJING: a Monte Carlo model for multiple jet production in pp, pA and AA collisions . Phys. Rev. D 44 , 3501 - 3516 ( 1991 ). doi:10.1103/PhysRevD.44.3501 54. R. Brun , F. Carminati , S. Giani , GEANT detector description and simulation tool . CERN Program Library Long Write-up W5013 (1994) 55. T. Sjöstrand , S. Mrenna , P.Z. Skands , A brief introduction to PYTHIA 8.1. Comput . Phys. Commun . 178 , 852 - 867 ( 2008 ). doi:10.1016/j.cpc. 2008 .01.036. arXiv: 0710 .3820 [hep-ph] 56. P. Nason , A new method for combining NLO QCD with shower Monte Carlo algorithms . JHEP 11 , 040 ( 2004 ). doi:10.1088/ 1126 - 6708 /2004/11/040. arXiv:hep-ph/ 0409146 57. S. Frixione , P. Nason , C. Oleari , Matching NLO QCD computations with parton shower simulations: the POWHEG method . JHEP 11 , 070 ( 2007 ). doi:10.1088/ 1126 - 6708 /2007/11/070. arXiv: 0709 . 2092 [hep-ph] 58. H.J. Drescher , M. Hladik , S. Ostapchenko , T. Pierog , K. Werner , Parton based Gribov-Regge theory . Phys. Rep . 350 , 93 - 289 ( 2001 ). doi:10.1016/ S0370-1573(00)00122-8. arXiv:hep-ph/0007198 59. K. Werner , I. Karpenko , T. Pierog , M. Bleicher , K. Mikhailov , Event-by-event simulation of the three-dimensional hydrodynamic evolution from flux tube initial conditions in ultrarelativistic heavy ion collisions . Phys. Rev. C 82 , 044904 ( 2010 ). doi:10.1103/ PhysRevC.82.044904. arXiv: 1004 .0805 [nucl-th] 60. K. Werner , B. Guiot , I. Karpenko , T. Pierog , Analysing radial flow features in p-Pb and p-p collisions at several TeV by studying identified particle production in EPOS3 . Phys. Rev. C89(6) , 064903 ( 2014 ). doi:10.1103/PhysRevC.89.064903. arXiv: 1312 .1233 [nucl-th] 61. M. Cacciari , S. Frixione , N. Houdeau , M.L. Mangano , P. Nason , G. Ridolfi , Theoretical predictions for charm and bottom production at the LHC . JHEP 10 , 137 ( 2012 ). doi:10.1007/ JHEP10(2012)137 . arXiv: 1205 .6344 [hep-ph] 62. M. Klasen , C. Klein-Bösing , K. Kovaric , G. Kramer , M. Topp , J.P. Wessels , NLO Monte Carlo predictions for heavy-quark production at the LHC: pp collisions in ALICE . JHEP 08 , 109 ( 2014 ). doi:10.1007/ JHEP08(2014)109 . arXiv: 1405 .3083 [hep-ph] 63. S. Alioli , P. Nason , C. Oleari , E. Re , A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX . JHEP 06 , 043 ( 2010 ). doi:10.1007/ JHEP06(2010)043 . arXiv: 1002 .2581 [hep-ph] 64. S. Frixione , P. Nason , G. Ridolfi , A positive-weight nextto-leading-order Monte Carlo for heavy flavour hadroproduction . JHEP 09 , 126 ( 2007 ). doi:10.1088/ 1126 - 6708 /2007/09/126. arXiv: 0707 .3088 [hep-ph] 65. H.-L. Lai , M. Guzzi , J. Huston , Z. Li , P.M. Nadolsky , J. Pumplin , C.P. Yuan , New parton distributions for collider physics . Phys. Rev. D 82 , 074024 ( 2010 ). doi:10.1103/PhysRevD.82.074024. arXiv: 1007 .2241 [hep-ph] 66. K.J. Eskola , H. Paukkunen , C.A. Salgado , EPS09: a new generation of NLO and LO nuclear parton distribution functions . JHEP 04 , 065 ( 2009 ). doi:10.1088/ 1126 - 6708 /2009/04/065. arXiv: 0902 .4154 [hep-ph] 67. Particle Data Group Collaboration , K.A. Olive et al., Review of particle physics (RPP). Chin. Phys. C 38 , 090001 ( 2014 ). doi:10. 1088/ 1674 - 1137 /38/9/090001 68. CDF Collaboration, T. Affolder et al., Charged jet evolution and the underlying event in p p¯ collisions at 1.8 TeV . Phys. Rev. D 65 , 092002 ( 2002 ). doi:10.1103/PhysRevD.65.092002 69. CMS Collaboration, V. Khachatryan et al., First measurement of the underlying event activity at the LHC with √s = 0.9 TeV. Eur. Phys. J. C 70 , 555 - 572 ( 2010 ). doi:10.1140/epjc/s10052- 010 - 1453 -9. arXiv: 1006 . 2083 [hep-ex] 70. ATLAS Collaboration, G. Aad et al., Measurement of underlying event characteristics using charged particles in pp collisions at √s = 900 GeV and 7 TeV with the ATLAS detector . Phys. Rev. D 83 , 112001 ( 2011 ). doi:10.1103/PhysRevD.83.112001. arXiv: 1012 .0791 [hep-ex] 71. ALICE Collaboration, B. Abelev et al., Underlying event measurements in pp collisions at √s = 0.9 and 7 TeV with the ALICE experiment at the LHC. JHEP 07, 116 (2012). doi:10.1007/ JHEP07(2012)116. arXiv:1112.2082 [hep-ex] J. Adam39, D. Adamová85, M. M. Aggarwal89, G. Aglieri Rinella35, M. Agnello31,112, N. Agrawal48, Z. Ahammed136, S. Ahmad18, S. U. Ahn69, S. Aiola140, A. Akindinov55, S. N. Alam136, D. S. D. Albuquerque123, D. Aleksandrov81, B. Alessandro112, D. Alexandre103, R. Alfaro Molina64, A. Alici12,106, A. Alkin3, J. R. M. Almaraz121, J. Alme22,37, T. Alt42, S. Altinpinar22, I. Altsybeev135, C. Alves Garcia Prado122, C. Andrei79, A. Andronic99, V. Anguelov95, T. Anticˇic´100, F. Antinori109, P. Antonioli106, L. Aphecetche115, H. Appelshäuser61, S. Arcelli27, R. Arnaldi112, O. W. Arnold36,96, I. C. Arsene21, M. Arslandok61, B. Audurier115, A. Augustinus35, R. Averbeck99, M. D. Azmi18, A. Badalà108, Y. W. Baek68, S. Bagnasco112, R. Bailhache61, R. Bala92, S. Balasubramanian140, A. Baldisseri15, R. C. Baral58, A. M. Barbano26, R. Barbera28, F. Barile33, G. G. Barnaföldi139, L. S. Barnby35,103, V. Barret71, P. Bartalini7, K. Barth35, J. Bartke119,a, E. Bartsch61, M. Basile27, N. Bastid71, S. Basu136, B. Bathen62, G. Batigne115, A. Batista Camejo71, B. Batyunya67, P. C. Batzing21, I. G. Bearden82, H. Beck61,95, C. Bedda112, N. K. Behera51, I. Belikov65, F. Bellini27, H. Bello Martinez2, R. Bellwied125, R. Belmont138, E. Belmont-Moreno64, L. G. E. Beltran121, V. Belyaev76, G. Bencedi139, S. Beole26, I. Berceanu79, A. Bercuci79, Y. Berdnikov87, D. Berenyi139, R. A. Bertens54, D. Berzano35, L. Betev35, A. Bhasin92, I. R. Bhat92, A. K. Bhati89, B. Bhattacharjee44, J. Bhom119, L. Bianchi125, N. Bianchi73, C. Bianchin138, J. Bielcˇík39, J. Bielcˇíková85, A. Bilandzic36,82,96, G. Biro139, R. Biswas4, S. Biswas4,80, S. Bjelogrlic54, J. T. Blair120, D. Blau81, C. Blume61, F. Bock75,95, A. Bogdanov76, H. Bøggild82, L. Boldizsár139, M. Bombara40, M. Bonora35, J. Book61, H. Borel15, A. Borissov98, M. Borri84,127, F. Bossú66, E. Botta26, C. Bourjau82, P. Braun-Munzinger99, M. Bregant122, T. Breitner60, T. A. Broker61, T. A. Browning97, M. Broz39, E. J. Brucken46, E. Bruna112, G. E. Bruno33, D. Budnikov101, H. Buesching61, S. Bufalino31,35, S. A. I. Buitron63, P. Buncic35, O. Busch131, Z. Buthelezi66, J. B. Butt16, J. T. Buxton19, J. Cabala117, D. Caffarri35, X. Cai7, H. Caines140, L. Calero Diaz73, A. Caliva54, E. Calvo Villar104, P. Camerini25, F. Carena35, W. Carena35, F. Carnesecchi12,27, J. Castillo Castellanos15, A. J. Castro128, E. A. R. Casula24, C. Ceballos Sanchez9, J. Cepila39, P. Cerello112, J. Cerkala117, B. Chang126, S. Chapeland35, M. Chartier127, J. L. Charvet15, S. Chattopadhyay136, S. Chattopadhyay102, A. Chauvin36,96, V. Chelnokov3, M. Cherney88, C. Cheshkov133, B. Cheynis133, V. Chibante Barroso35, D. D. Chinellato123, S. Cho51, P. Chochula35, K. Choi98, M. Chojnacki82, S. Choudhury136, P. Christakoglou83, C. H. Christensen82, P. Christiansen34, T. Chujo131, S. U. Chung98, C. Cicalo107, L. Cifarelli12,27, F. Cindolo106, J. Cleymans91, F. Colamaria33, D. Colella35,56, A. Collu75, M. Colocci27, G. Conesa Balbastre72, Z. Conesa del Valle52, M. E. Connors140,b, J. G. Contreras39, T. M. Cormier86, Y. Corrales Morales26,112, I. Cortés Maldonado2, P. Cortese32, M. R. Cosentino122,124, F. Costa35, J. Crkovská52, P. Crochet71, R. Cruz Albino11, E. Cuautle63, L. Cunqueiro35,62, T. Dahms36,96, A. Dainese109, M. C. Danisch95, A. Danu59, D. Das102, I. Das102, S. Das4, A. Dash80, S. Dash48, S. De122, A. De Caro12,30, G. de Cataldo105, C. de Conti122, J. de Cuveland42, A. De Falco24, D. De Gruttola12,30, N. De Marco112, S. De Pasquale30, R. D. De Souza123, A. Deisting95,99, A. Deloff78, E. Dénes139,a, C. Deplano83, P. Dhankher48, D. Di Bari33, A. Di Mauro35, P. Di Nezza73, 1 A.I. Alikhanyan National Science Laboratory (Yerevan Physics Institute) Foundation , Yerevan, Armenia 2 Benemérita Universidad Autónoma de Puebla, Puebla, Mexico 3 Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-017-4779-8.pdf

Measurement of azimuthal correlations of D mesons with charged particles in pp collisions at \(\sqrt{s}=7\) TeV and p–Pb collisions at \({\sqrt{{{s}}_{\scriptscriptstyle {\mathrm{NN}}}}}=5.02\) TeV, The European Physical Journal C, 2017, DOI: 10.1140/epjc/s10052-017-4779-8