Measurement of the production of charm jets tagged with D$^{0}$ mesons in pp collisions at $\sqrt{s}$= 7 TeV

The production of charm jets in proton-proton collisions at a center-of-mass energy of $\sqrt{s}=7$ TeV was measured with the ALICE detector at the CERN Large Hadron Collider. The measurement is based on a data sample corresponding to a total integrated luminosity of $6.23$ ${\rm nb}^{-1}$, collected using a minimum-bias trigger. Charm jets are identified by the presence of a D$^0$ meson among their constituents. The D$^0$ mesons are reconstructed from their hadronic decay D$^0\rightarrow$K$^{-}\pi^{+}$. The D$^0$-meson tagged jets are reconstructed using tracks of charged particles (track-based jets) with the anti-$k_{\mathrm{T}}$ algorithm in the jet transverse momentum range $5<p_{\rm{T,jet}}^{\mathrm{ch}}<30$ ${\rm GeV/}c$ and pseudorapidity $|\eta_{\rm jet}|<0.5$. The fraction of charged jets containing a D$^0$-meson increases with $p_{\rm{T,jet}}^{\rm{ch}}$ from $0.042 \pm 0.004\, \mathrm{(stat)} \pm 0.006\, \mathrm{(syst)}$ to $0.080 \pm 0.009\, \rm{(stat)} \pm 0.008\, \rm{(syst)}$. The distribution of D$^0$-meson tagged jets as a function of the jet momentum fraction carried by the D$^0$ meson in the direction of the jet axis ($z_{||}^{\mathrm{ch}}$) is reported for two ranges of jet transverse momenta, $5<p_{\rm{T,jet}}^{\rm{ch}}<15$ ${\rm GeV/}c$ and $15<p_{\rm{T,jet}}^{\rm{ch}}<30$ ${\rm GeV/}c$ in the intervals $0.2<z_{||}^{\rm{ch}}<1.0$ and $0.4<z_{||}^{\rm{ch}}<1.0$, respectively. The data are compared with results from Monte Carlo event generators (PYTHIA 6, PYTHIA 8 and Herwig 7) and with a Next-to-Leading-Order perturbative Quantum Chromodynamics calculation, obtained with the POWHEG method and interfaced with PYTHIA 6 for the generation of the parton shower, fragmentation, hadronisation and underlying event.


te for Nuclea
Research (JINR)
DubnaRussia

Korea Institute of Science and Technology Information
DaejeonRepublic of Korea

KTO Karatay University
KonyaTurkey

Laboratoire de Physique Subatomique et de Cosmologie
Université Grenoble-Alpes
CNRS-IN2P3
GrenobleFrance

Lawrence Berkeley National Laboratory
BerkeleyCaliforniaUnited States

Department of Physics
Division of Particle Physics
Lund University
LundSweden

Nagasaki Institute of Applied Science
NagasakiJapan

Nara Women's University (NWU)
NaraJapan

School of Science
Department of Physics
National and Kapodistrian University of Athens
AthensGreece

National Centre for Nuclear Research
WarsawPoland

National Institute of Science Education and Research
Homi Bhabha National Institute
JatniIndia

National Nuclear Research Center
BakuAzerbaijan

National Research Centre Kurchatov Institute
MoscowRussia

Niels Bohr Institute
University of Copenhagen
CopenhagenDenmark

National institute for subatomic physics
Nikhef, AmsterdamNetherlands

NRC Kurchatov Institute IHEP
ProtvinoRussia

NRNU Moscow Engineering Physics Institute
Mo y Laboratory
DaresburyUnited Kingdom Nuclear Physics Institute of the Czech Academy of Sciences
Řež uPrahyCzech Republic

Oak Ridge National aboratory
Oak RidgeTennesseeUnited States

Ohio State University
ColumbusOhioUnited States

Petersburg Nuclear Physics Institute
GatchinaRussia

Physics department
Faculty of science
University of Zagreb
ZagrebCroatia

Physics Department
Panjab University
ChandigarhIndia

Physics Department
University of Jammu
JammuIndia

Physics Department
University of Rajasthan
101 Physikalisches InstitutJaipurIndia

Eberhard-Karls-Universität Tübingen
102 Physikalisches InstitutTübingenGermany

Ruprecht-Karls-Universität Heidelberg
HeidelbergGermany

103 Physik Department
Technische Universität München
104 Politecnico di BariMunich, BariGermany, Italy

105 Research Division and ExtreMe Matter Institute EMMI
GSI Helmholtzzentrum für Schwerionenforschung GmbH
DarmstadtGermany

106 Rudjer Bošković Institute
ZagrebCroatia

Russian Federal Nuclear Center (VNIIEF)
SarovRussia

Saha Institute of Nuclear Physics
Homi Bhabha National Institute
KolkataIndia

University of Cape Town ape TownSouth Africa

University of Houston
HoustonTexasUnited States

University of Jyväskylä
JyväskyläFinland

University of Liverpool
LiverpoolUnited Kingdom

University of Science and Techonology of China
HefeiChina

University of Tennessee
KnoxvilleTennesseeUnited States

University of the Witwatersrand
JohannesburgSouth Africa

University of Tokyo
TokyoJapan

University of Tsukuba
TsukubaJapan

Université Clermont Auvergne
CNRS/IN2P3
LPC
Clermont-FerrandFrance

Université de Lyon
Université Lyon 1
CNRS/IN2P3
IPN-Lyon
Villeurbanne, LyonFrance

UMR 7178
Université de Strasbourg
CNRS
IPHC
F-67000Strasbourg, StrasbourgFrance, France

Départment de Physique Nucléaire (DPhN)
Université Paris-Saclay Centre d'Etudes de Saclay (CEA)
IRFU
SaclayFrance

Università degli Studi di Foggia
FoggiaItaly

Università degli Studi di Pavia
PaviaItaly

Università di Brescia
BresciaItaly

141 Variable Energy Cyclotron Centre
Homi Bhabha National Institute
KolkataIndia

Warsaw University of Technology
WarsawPoland

Wayne State University
DetroitMichiganUnited States

Institut für Kernphysik
Wigner Research Centre for Physics
Westfälische Wilhelms-Universität Münster
145MünsterGermany

Hungarian Academy of Sciences
BudapestHungary

Yale University
New Haven, ConnecticutUnited St tes

Yonsei University
SeoulRepublic of Korea

EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH
25 April 20198968AFBBC1E08F361CF2767DB67D6A69arXiv:1905.02510v2[nucl-ex]
The production of charm jets in proton-proton collisions at a center-of-mass energy of √ s = 7 TeV was measured with the ALICE detector at the CERN Large Hadron Collider.The measurement is based on a data sample corresponding to a total integrated luminosity of 6.23 nb −1 , collected using a minimum-bias trigger.Charm jets are identified by the presence of a D 0 meson among their constituents.The D 0 mesons are reconstructed from their hadronic decay D 0 →K − π + .The D 0 -meson tagged jets are reconstructed using tracks of charged particles (track-based jets) with the anti-k T algori hm in the jet transverse momentum range 5 < p ch T,jet < 30 GeV/c and pseudorapidity |η jet | < 0.5.The fraction of charged jets containing a D 0 -meson increases with p ch T,jet from 0.042 ± 0.004 (stat) ± 0.006 (syst) to 0.080 ± 0.009 (stat) ± 0.008 (syst).The distribution of D 0 -meson tagged jets as a function of the jet momentum fraction carried by the D 0 meson in the direction of the jet axis (z ch | | ) is reported for two ranges of jet transverse momenta, 5 < p ch T,jet < 15 GeV/c and 15 < p ch T,jet < 30 GeV/c in the intervals 0.2 < z ch | | < 1.0 and 0.4 < z ch | | < 1.0, respectively.The data are compared with results from Monte Carlo event generators (PYTHIA 6, PYTHIA 8 and Herwig 7) and with a Next-to-Leading-Order perturbative Quantum Chromodynamics calculation, obta ned with the POWHEG method and interfaced with PYTHIA 6 for the generation of the parton shower, fragmentation, hadronisation and underlying event.

Introduction

The study of heavy-flavour production in high-energy interactions provides important tests for Quantum-Chromodynamics (QCD) calculations [1][2][3].The transverse-momentum (p T )-differential production cross section of D mesons from charm-quark fragmentation (referred to as "prompt" D mesons) was measured in proton-proton (pp) and pp collisions at several center-of-mass energies, from √ s = 0.2 TeV at RHIC up to the energies of Tevatron ( √ s = 1.96TeV) and the LHC ( √ s = 13 TeV) [4][5][6][7][8][9][10][11][12].The data are described reasonably well by calculatio

based on perturbative QCD
(pQCD) that rely either on the collinear-factorisation approach, like FONLL [13][14][15] and GM-VFNS [16], or on the k T -factorisation approach [17].In comparison to single-particle measurements, the reconstruction of jets containing charm hadrons allows for more differential studies to characterise the heavy quark production and fragmentation.A relevant observable is the fraction (z | | ) of the jet momentum ( ì p jet ) carried by the D meson along the jet axis direction:
z | | = ì p jet • ì p D ì p jet • ì p jet ,(1)
where ì p D is the D-meson momentum.

Pioneering measurements of charm jets were performed at the CERN SPS [18] and at the Tevatron [19,20].

The STAR experiment at RHIC mea ured the D * ± -meson production in jets in pp collisions at √ s = 200 GeV [21].The jets were measured in the interval 8 < p T,jet < 20 GeV/c.The yield at low z | | values is higher than that obtained with a Monte Carlo simulation performed with PYTHIA 6 [22] using only the direct charm flavour creation processes, gg → cc and qq → cc.This suggests that higher order processes (gluon splitting, flavour excitation) are not negligible in the charm production at RHIC energies.In a more recent analysis, the PHENIX collaboration measured azimut al correlations of charm and bottom hadrons in their semi-leptonic decays using unlike-and like-sign muon pairs [23].Overall they found good agreement with a PYTHIA 6 [22] simulation.Through a Bayesian analysis based on PYTHIA 6 templates, the PHENIX collaboration found that while leading order pair creation is dominant for bottom production, higher order processes dominate for charm one.

At the LHC, the analysis of the angular correlations of b-hadron decay vertices, measured by CMS [24], indicated that the collinear region, where the contributions of gluon splitting processes are expected to be large, is not adequately described by PYTHIA 6 nor by predictions based on Next-To-Leading (NLO) order QCD calculations.The ATLAS experiment measured the D * ± -meson production in jets in pp collisions at √ s = 7 TeV [25], finding that the z | | distribution differs from expectations of PYTHIA 6, HERWIG 6 [26,27] and POWHEG [28][29][30][31] event generators, both in overall normalization and shape, with data displaying a higher probability for low z | | values and a steeper decrease towards z | | = 1.The discrepancy between data and generator expectations is maximum in the lowest jet p T interval, 25 < p T < 30 GeV/c.The ATLAS data are well described in a recent global QCD analysis of fragmentation functions based on the ZM-VFNS [32] scheme, in which the in-jet fragmentation data were combined with previous D- eson measurements in a global fit [33].This global QCD analysis evidences the importance of in-jet fragmentation data in order to pin down the otherwise largely unconstrained momentum fraction dependence of the gluon fragmentation function.

In this paper, we report the first ALICE measurements of the D 0 -meson tagged track-based jet p Tdifferential cross section in pp collisions t √ s = 7 TeV and of the D 0 -meson z ch | | distribution.The z ch | | is defined as in Eq. 1 but using the momenta of the track-based jet ì p ch jet .With track-based jets we indicate jets reconstructed with only their charged-particle constituents [34].As described in Section 2, the excellent low-and intermediate-momentum tracking capabilities of the ALICE apparatus allow the measurement of jets at very low p T , particularly in the charged jet transverse momentum range 5 < p ch T,jet < 30 GeV/c considered in this paper.This kinematic region is still larg ly unconstrained by previous measurements.

The measurements reported in this paper are also important to define a pp reference baseline for future measurements in Pb-Pb and p-Pb collisions at the LHC.Charm quarks, interacting with the constituents of the Quark-Gluon Plasma formed in these collisions, lose energy via both radiative and collisional processes, as evidenced by the strong suppression of high-p T D-meson production measured by AL-ICE [35][3

S [38].Contrary to si
gle particles, jets allow one to capture more details of the parton shower dynamics in the medium.In particular, the study of jet substructure, pioneered for QCD studies and beyond standard model searches [39], can be important to investigate the microscopic properties of hadronic matter at high densities and temperatures [40][41][42][43].

The paper is structured as follows.Section 2 describes the components of the ALICE apparatus, the data sample and Monte Carlo simulations used in the analysis.In Section 3, the analysis procedure to obtain the raw spectrum of D 0 -meson tagged jets and the z ch | distribution is outlined.Section 4 describes several corrections that are required to account for the D 0 -meson and jet reconstruction efficiency, the jet momentum scale and the contribution from D 0 mesons coming from b-hadron decays.The systematic uncertainties affecting the measurements are reported in Section 5.The results and physics implications are discussed in Section 6.Finally, Section 7 closes the paper with conclusions and future perspectives.


Apparatus and data sample

The measurements presented in this paper were carried out using data recorded by the ALICE apparatus [ 4,45] in 2010.ALICE is composed of a central barrel embedded in a 0.5 T magnetic field parallel to the beam direction (z axis in the ALICE reference frame) and a set of forward-and backward-rapidity detectors.The Inner Tracking System (ITS) and the Time-Projection Chamber (TPC) were used for charged-particle track reconstruction and the combined information from the TPC and the Time-Of-Flight (TOF) detectors was used to provide particle identification (PID).These detectors are located in the central barrel, which has a full azimuthal coverage and a pseudorapidity interval of |η| < 0.9.

The ITS is the closest detector to the interaction point and consists of six cylindrical layers of silicon detectors, using three dif erent technologies: Silicon Pixel Detectors (SPD), whose radius of the first layer is 3.9 cm, Silicon Drift Detectors (SDD) and Silicon Strip Detectors (SSD).The proximity of the SPD to the interaction point, combined with its high spatial resolution, provides a resolution on the track impact parameter with respect to the primary vertex better than 75 μm for tracks with transverse momentum p T > 1 GeV/c.The TPC consists of a 510 cm long cylinder with an inner radius of 85 cm and an outer radius of 250 cm.The detector is divided into two halves at the center by a high voltage electrode that generates a uniform electric field in the longitudinal direction pointing from the endplates to the center.The TPC is filled with a mixture of Ne (90%) and CO 2 (10%) gases.The trajectories of charged particles traversing the TPC volume are reconstructed from the ionisation produced in the gas.The ALICE apparatus is capable of reconstructing charged-particle tracks down to p T = 0.15 GeV/c with a p T -resolution better than 2% up to p T = 20 GeV/c.The PID information from the TPC is based on the par

cle specific ionisation energy loss dE/dx
n the gas.The TOF provides particle identification based on the time-of-flight of the particle from the interaction point to the hit in the Multi-Gap Resistive Plate Chambers (MRPCs) that compose the detector.For events with sufficiently large multiplicity, the best estimate of the collision time is obtained from the particle arrival times at the TOF [46]; for lower-multiplicity events the collision time is measured by the T0 detector, which consists of two arrays of Cherenkov counters located at +350 cm and −70 cm along the beam line.The combined PID information from both detectors provides up to 3σ separation power for pions/kaons in the range 0.5 < p T < 2 GeV/c [47].

The V0 detector was used for triggering minimum-bias events.The detector consists of two scintillator arrays located around the beam pipe on each side of the interaction point covering the pseudo-rapidity interval −3.7 < η < −1.7 and 2.8 < η < 5.1, respectively.The minimum-bias condition is defined by the presence of at least one hit in one of the V0 scintillators or in the SPD.

I the work presented in this paper, pp collisions at √ s = 7 TeV were analysed.The sample consists of about 388 × 10 6 minimum-bias events, corresponding to an integrated luminosity of L int = 6.23 nb −1 [48].Events were selected offline by using the timing information from the V0 and the correlation between the number of hits and track segments in the SPD detector to remove background due to beam-gas interactions.Only events with the primary vertex reconstructed within |z| < 10 cm with respect to the center of the detector were used for this analysis.

Monte Carlo (MC) simulations were employed to calculate corrections as described in Section 4. The simulations were performed using PYTHIA 6.4.24 [22] with the Perugia 2011 tune [49].The generated particles were transported through the ALICE apparatus using the GEANT3 transport model [50].The luminous region distribution, the geometry of the apparatus, as well as the conditions of all the ALICE detectors were reproduced in detail in the simulations.


Analysis


D 0 -meson selection

The D 0 mesons were reconstructed via their hadronic decay D 0 → K − π + (and charge conjugate) which has a branching ratio of (3.89 ± 0.04)% [51].In each event, D 0 -meson candidates and their decay vertices were constructed from pairs of tracks with opposite charge.The tracks were required to have |η| < 0.8, p T > 0.3 GeV/c, at least 70 associated TPC space points (out of a maximum of 159), χ 2 /ndf < 4 in the TPC (where ndf is the number of degrees of freedom involved in the tracking procedure), at least one hit in either of the two layers of the SPD and a minimum of 3 hits in the entire ITS.

The D 0 -meson selection criteria were established in previously published works by the ALICE Collaboration [7,9].D 0 mesons were required to be within the rapidity interval comprised by a fiducial detector acceptance, |y| < y fid (p T,D ), with y fid (p T,D ) increasing from 0.5 to 0.8 in the D 0 -meson transverse momentum interval 2 < p T,D < 5 GeV/c and y fid (p T,D ) = 0.8 for p T,D > 5 GeV/c.Outside of this selection the D 0 -meson reconstruction efficiency drops rapidly as a consequence of the detector pseudorapidity acceptance and the kinematic selections applied on the tracks.

In order to suppress the combinatorial background, we exploited the specific decay topology of the D 0 mesons.D 0 mesons have a mean proper decay length cτ = 123 μm [51].Their d

ay vertices are therefore typically dis
laced by a few hundred μm from the primary vertex of the interaction.The selection requirements were tuned to maximise the statistical significance of the signal along with good reconstruction efficiency.The geometrical selections were based on the displacement of the tracks from the interaction vertex, the distance between the D-meson decay vertex and the primary vertex (decay length, L) and the pointing of the reconstructed D-meson momentum to the primary vertex in the laboratory reference frame.

Further reduction of the combinatorial background was achieved by applying PID to the decay track candidates.The PID selection is based on the dE/dx and the time-of-flight signals measured with the TPC and TOF detectors, respectively.The selection was applied by requiring that the difference between the measured and expected PID signals was below 3σ, where σ is the experimental uncertainty associated with the measured signals.Based on the PID information, D 0 -meson candidates were accepted (as D 0 , D 0 , or both) or rejected, according to the compatibility with the K ∓ π ± final state.In the cases where both decay track candidates are found to be compatible with both the kaon and pion hypotheses, the D 0 -meson candidate was considered twice in either mass combinations corresponding to one of the two possible final states K − π + and K + π − .The candidates corresponding to a real D 0 meson but with the wrong daughter particle mass assignment are referred to as reflections.This component of the background was subtracted using Monte Carlo templates as described in Section 3.3.


Jet reconstruction and D 0 -meson tagging

For jet reconstruction, looser track selection criteria were employed as ompared to those used to identify D 0 -meson candidates.The pseudorapidity and momentum acceptance windows were extended to |η| < 0.9 and p T > 0.15 GeV/c, respectively.The requirement on the SPD hits was lifted to increase the trackreconstruction efficiency and improve its uniformity as a function of η and azimuthal angle ϕ.Tracks without SPD hits were required to contain at least 3 hits in the ITS and were constrained to the primary vertex of the interaction.Tracks without SPD hits comprise about 19% of the track sample used for jet reconstruction.The track reconstruction efficiency obtained with these selection criteria is uniform as a function of η and ϕ.As a function of the track transverse momentum, the efficiency is about 70% for p T t 85%.The relative track transverse momentum resolution is better than 3% in the range 0.15 < p T < 40 GeV/c.Jet reconstruction was performed with the anti-k T algorithm [52], as implemented in the FastJet [53] software package, with a resolution parameter R = 0.4 and the p T recombination scheme.From simple kinematic considerations we evaluated that more than 50% of the D 0 mesons with p T,D = 3 GeV/c have their decay products emitted at an angle larger than 0.4 rad with respect to the D 0 -meson momentum direction.This fraction approaches zero only for p T,D > 7 GeV/c.As a consequence, the decay products of low-momentum D 0 mesons are often found outside of the reconstructed jet cone that is physically correlated with the D 0 meson.It follows that the decay products of a single D 0 -meson candidate may be wrongly associated to two different j

candidates
n the jet finding phase.In order to avoid ambiguities in the charm jet tagging and to im ring that pairs of kaons and pions identified as the decay products of the same D 0meson candidate are part of the same jet.This constraint was implemented by removing the 4-momenta of the decay products of identified D 0 -meson candidates from the pool of particle tracks used in the jet finding, and replacing them with the 4-momenta of their respectively associated D 0 -meson candidates.Events containing more than one D 0 -meson candidate passing all the selection riteria are very rare and amount to approximately 0.9% of the events that contain at least one accepted candidate.In these cases, the jet reconstruction procedure was repea uded in the jet reconstruction as single tracks.This ensures that the combinatorial background of K − π + track pairs, which dominates the D 0 -meson candidates at low p T , does not influence the reconstruction of signal jets.Jets containing a D 0meson candidate among their constituents were tagged and retained for the next steps of the analysis.Jets with p ch T,jet > 5 GeV/c and |η jet | < 0.5 were accepted.The requirement on the jet pseudorapidity ensures that jets are fully contained in the detector acceptance.No correction to the reconstructed jet p T was performed to account for the background coming from the underlying event (UE), e.g.via multi-parton interactions (MPI).


D 0 -meson tagged jet yield extraction

The jet raw yields were extracted with an invariant mass analysis of

he D 0 -meson candidates u
ed to tag the charm jet candidates.These candidates were first divided in bins of p T,D .For each interval of p T,D the invariant mass distribution was fit with a function composed of a Gaussian function for the signal and an exponential term for the background.The position m fit and width σ fit of the D 0 -meson invariant mass peak were extracted from the corresponding parameters of the Gaussian component of the fit function.The top panels in Fig. 1  The contribution from residual D 0 -meson reflections not rejected by PID was accounted for by including in the fit a template consisting of the sum of two Gauss an functions with centroids and widths fixed to values obtained in the simulation.The amplitudes were normalized using the signal observed in data, keeping the ratio of the reflection component over the D 0 -meson signal fixed to the value obtained in the Monte Carlo simulation.In the wide invariant mass interval 1.715 < m inv < 2.015 GeV/c 2 used in the fitting procedure, the reflections over signal ratio varies in the range 0.15 − 0.30 as a function of p T,D .

The peak region contains a mixture of signal and combinatorial background, while the side-band region is far enough from the D 0 -meson peak to be signal-free.The total background N bkg (p T,D ) under the peak was extracted from the exponential and reflection components of the invariant-mass fit function by integrating them in the interval |m inv − m fit | < 2 σ fi tributions as a function of these observables are constructed for both the peak and side-band regions.The side-band distribution is scaled such that its total integral is equal to N bkg (p T,D ) and then subtracted from the peak-region distribution to obtain the raw yield as a function of p ch T,jet :
N raw (p T,D , p ch T,jet ) = N PR (p T,D , p ch T,jet ) − N bkg (p T,D ) N tot,SB (p T,D ) N SB (p T,D , p ch T,jet ),(2)
where
N raw (p T,D , p ch T,jet ), N PR (p T,D , p ch T,jet

Corrections

The p T -differential cross section of

arm jets tagged with D 0 mesons is def
ned as:
d 2 σ dp ch T,jet dη jet (p ch T,jet ) = 1 L int 1 BR N(p ch T,jet ) ∆η jet ∆p ch T,jet ,(3)
where N(p ch T,jet ) is the measured yield in each bin of p ch T,jet , corrected for the reconstruction efficiency, acceptance and b-hadron feed-down fraction, and unfolded for the detector jet momentum resolution; ∆p ch T,jet is the width of the histogram bin; ∆η jet = 1 is the jet reconstruction acceptance.Details on the corrections are discussed in the following sections.The reported yield includes jets containing either D 0 or D 0 mesons with p T,D > 3 GeV/c.The distribution of the jet momentum fraction carried by the D 0 meson in the direction of the jet axis (z ch | | ) is reported as a differential cr ss section defined as:
d 3 σ dz ch | | dp ch T,jet dη jet (p ch T,jet , z ch | | ) = 1 L int 1 BR N(p ch T,jet , z ch | | ) ∆η jet ∆p ch T,jet ∆z ch | | . (4)
The distribution was measured in the range 0.     ) are normalized such that one count corresponds to a single D 0 meson.It follows that a jet containing two D 0 mesons will enter this definition twice.While this choice may seem unnatural for the definition of a jet cross section, it has the advantage o having a model-independent tagging efficiency.In fact, if we were to count only once those jets containing two D 0 mesons, then their tagging efficiency would be twice as large.Then, the overall tagging efficiency would depend on the model-dependent fraction of jets with two D 0 mesons. 1


Reconstruction efficiency

The reconstruction efficiency of D 0 -meson tagged jets depends mainly on the track- econstruction efficiency and on the topological selections applied to find the D 0 -meson candidates.The efficiency was estimated with a Monte Carlo simulation using the PYTHIA 6 (Perugia 2011) [22,49] event generator and the GEANT3 [50] transport code.As shown in Fig. 2   In order to minimize the dependence of the efficiency correction on the fragmentation model and on the p T -spectrum shape of the simulated D 0 -meson tagged jet sample, the p ch T,jet distributions were multiplied by the inverse of the efficiency before summi g over p T,D :
N corr (p ch T,det jet ) = p T,D N raw (p T,D , p ch T,det jet ) P p T,D ,(5)
where N corr is the efficiency-corrected jet raw yield as a function of reconstructed jet transverse momentum p ch T,det jet , N raw (p T,D , p ch T,det jet ) was defined in Eq. 2, P p T,D is the prompt D 0 -meson reconstruction efficiency as a function of p T,D .The sum p T,D is intended over all p T,D ranges used in the invariant mass analysis (3 < p T,D < 30 GeV/c).The same procedure is applied to obtain the yields as a function of z ch | | .The corresponding equation is obtained by replacing p ch T,jet with z ch | | in Eq. 5.


Subtraction of the b-jet contribution

The efficiency of prompt D 0 -meson tagged jets is lower compared to the efficiency of those coming from the fragmentation of a beauty quark for whi

the non-p
ompt D 0 -meson is produced by the decay of a beauty hadron.The prompt and non-prompt acceptance and reconstruction efficiency correction factors are compared in Fig. 2 (right panel).Due to the longer decay length of beauty hadrons (cτ ≈ 500 μm [51]), some topological selections are more efficient for non-prompt D 0 mesons.The non-prompt efficiency is higher by about a factor 2 for p T,D = 3 GeV/c compared to the prompt efficiency.The separation between the two efficiencies decreases with p T,D , until they almost converge for p T,D > 15 GeV/c.

Due to the higher reconstruction efficiency of the non-prompt D 0 mesons, the natural admixture of the prompt and the non-prompt components is biased towards the non-prompt in a detector-and analysisspecific way.In order to simplify comparisons with other experimental results and theoretical calculations, the fraction of D 0 -meson tagged jets coming from the fragmentation of b quarks (via the decay of a beauty hadron into a D 0 ) was subtracted as follows.

The non-prompt fraction was estimated with POWHEG interfaced with the PYTHIA 6 (Perugia 2011) Monte Carlo parton shower.The decays of beauty adrons were turned off in PYTHIA 6, to allow EvtGen [54]  Th parton distribution function (PDF) was obtained using the LHAPDF 6 [55] interpolator with the PDF set CT10nlo [56].

where R NP is the matrix representing the p ch T,jet detector response for non-prompt D 0 -meson ta olding

The reconstructed jet momentum is affected by the finite detector resolution.The main factor impacting the jet momentum resolution is the track-reconstruction efficiency, which causes an average negative shift and a smearing of the reconstructed jet momentum compared to the true jet momentum.The detector resolution was quantified with the same Monte Carlo simulation used to estimate the efficiency.It was verified that the simulation is able to reproduce at the detector level the main features of the data, such as jet and D 0 -meson p T distributions, and the average number of jet constituents.D 0 -meson tagged jets at the detector level were uniquely matched with the corresponding jets at the generator level.The matching criteria are based on the presence of the same D 0 meson, which was followed from the generator level throughout its decay and transport in the detector volume.The jet transverse momentum resolution can be quantified from the probability density distribution of the relative difference (∆ p T ) between the reconstructed jet transverse momentum p ch T,det jet and the generated jet transverse momentum p ch T,gen jet :
∆ p T = (p ch T,det jet − p ch T,gen jet )/p ch T,gen jet .(7)
A similar quantity is defined for the jet momentum fraction carried by the D 0 :
∆ z = (z ch | |,det − z ch | |,gen )/z ch | |,gen .(8)
Figure 4 shows the probability density distributions of ∆ p T (left) and ∆ z (right) for a selection of p ch T,jet and z ch | | ranges.

The mean relative shift of the reconstructed jet momentum varies monotonically from −2% for p ch T,gen jet = 5 GeV/c to −7% for p ch T,gen jet = 30 GeV/c.The resolution, defined as the standard deviation fr m the mean of the probability density distribution, also varies monotonically as a function of p ch T,gen jet from 10% to 15%.The resolution is slightly better compared to the inclusive jet measurem nt performed on the same dataset with similar techniques [57].This difference can be ascribe to the requirement of the presence of a D 0 meson with p T,D > 3 GeV/c in the jet.T e finite detector resolution modifies the dependence of the measured yields as a function of p ch T,jet and z ch | | .The relationship between the raw and the generated yields can be written as:
N det (p ch T,det jet ) = R P (p ch ,det jet , p ch T,gen jet ) • N gen (p ch T,gen jet ),(9)N det (z ch | |,det ) = R P (z ch | |,det , z ch | |,gen ) • N gen (z ch | |,gen ),(10)
where R P is the matrix representing the p ch T,jet detector response for prompt D 0 -meson tagged jets; N det and N gen are the vectors corresponding to the measured and generated yields in bins of either p ch T,jet or z ch | | .

The effects of the limited detector resolution discussed above were corrected through an unfolding procedure.The measured distributions N sub were unfolded using an iterative approach based on Bayes' theorem [58].The iterative unfolding algorithm successfully converged after three iterations.The N sub (p ch T,jet ) distribution and the two N sub (z ch | | ) distributions for 5 < p ch T,jet < 15 GeV/c and for 15 < p ch T,jet < 30 GeV/c were each unfolded separately with their corresponding detector response matrices.For the z ch | | distributions an additional correction, based on the same PYTHIA 6 + GEANT3 simulation, was applied to account for the effect of the detector resolution on p ch T,jet , which causes jet candidates to fall in or out of the considered p ch T,jet intervals.This correction is about +15% for jets tagged with D 0 mesons with p T,D < 5 GeV/c and negligible for p T,D > 5 GeV/c.The same histogram binning was used for the measured and the unfolded distributions.Underflow (p ch T,jet < 5 GeV/c) and overflow (p ch T,jet > 30 GeV/c) bins were excluded from the p ch T,det jet axis of the response matrix, but were kept as degrees of freedom in the p ch T,gen jet axis that could be populated according to the probabilities mapped by the response matrix.The same applies, only for the underflow bin, for the z ch | | distributions (z ch | | ≤ 1 by constr

tion).The overall unfoldi
g corrections on the yields are: between +2% and +14% for the p ch T,jet distribution; between −6% and +5% for the z ch | | distribution with 5 < p ch T,jet < 15 GeV/c; between −30% and +10% for the z ch | | distribution with 15 < p ch T,jet < 30 GeV/c.In all cases, the unfolding correction is smaller than the sta istical uncertainties.


Systematic uncertainties

The relative systematic uncertainties on the p T -differential D 0 -meson tagged jet cross section and on the z ch | | distributions for 5 < p ch T,jet < 15 GeV/c and 15 < p ch T,jet < 30 GeV/c are summarized in Tables 1, 2 and 3, respectively.In the following, each source of systematic uncertainty is discussed.

The uncertainty on the track-reconstruction efficiency affects our measurement via an uncertainty on the reconstruction efficiency of the D 0 meson and an uncertainty on the jet momentum resolution.For the D 0 -meson reconstruction efficiency, a p T -independent systematic uncertainty of 4% was assigned based on the D 0 -meson studies in [9].The relative systematic uncertainty on the track-reconstruction efficiency for the set of tracks used for jet reconstruction was estimated to be 5% in [57].Therefore, the detector response matrix was modified by randomly rejecting 5% of the tracks reconstructed in the detector simulation.The jet p T distribution and z ch | | distributions were unfolded using this modified matrix and compared with the distributions unfolded with the nominal matrix.The relative differences were found to be less than 7% in most cases.The uncertainty on the track momentum resolution was determined to have a negligible effect on the jet momentum resolution.Tracks of charged particles produced in the decays of neutral strange hadrons or in secondary interactions with the detector material (including photon conversions) are largely suppressed by the track selection criteria used in the jet finding.The residual contamination is reproduced fairly well by the Monte Carlo simulation used to estimate the detector response: this residual contamination is corrected for in the unfolding procedure.However, PYTHIA 6 does not adequately reproduce the strange particle production [59].An uncertainty of about 0.5% on the jet momentum arises from this [60], which causes an uncertainty of about 2% on the p Tand z ch | | -differential yields.A possible influence of the si ulated p T -spectrum shape of charm jets on the D 0 -meson reconstruction efficiency was investigated by re-calculating the corrections using an alternative p T -spectrum shape obtained from an independent simulation in which POWHEG replaced PYTHIA 6 for the generation of the hard scattering.The effect on the final results was found to be negligible.

Discrepancies between simulation and data that affect the D 0 -meson reconstruction and selection efficiency introduce a systematic uncertainty.For example, the selections based on the displacement of the decay vertex from the collision point are sensitive to the resolution on the track impact parameter with respect to the primary vertex; residual misalignment of the silicon pixel detector can also intro-duce irreducible differences between data and simulation.The systematic uncertainty arising from these discrepancies was determined by repeating the analysis with different sets of selection criteria.In the D-meson cross-section analysis [9], this uncertainty was estimated to be 5%.Since the uncertainty was also found to depend weakly on the p T,D , for the range considered in this analysis, and the D 0 -meson tagged jet reconstruction efficiency does not depend significantly on the p ch T,jet (see Fig. 2), the same value of 5% was assigned as uncertainty on the yield for this measurement.

The systematic uncertainties on the raw yield extraction were estimated by repeating the fitting procedure of the invariant mass distributions several times with different fit conditions.The tests included the following: (i) variations of the upper and lower limits of the fit range; (ii) variations of the invariant-mass distribution bin width; (iii) background fit function (default exponential replaced by first and second order polynomial functions); (iv) mean and/or σ parameters of the Gaussian function fixed to the values expected from Monte Carlo simulations.The root-mean-square of the differences between signal-yield distributions obtained from the various trials was taken as the systematic uncertainty.An additional systematic uncertainty was assigned by varying the assumed relative contribution of the D 0 -meson reflections to the signal by ± 50%.The total uncertainty on the p T -differential jet cross section varies between 4-15% and rises with p ch T,jet .The uncertainty on the z ch | | distributions is 2-23% with higher values for the two lowest z ch | | intervals.In order to estimate the systematic uncertainty on the simulation used to subtract the b-hadron feed-down, the following parameters of the POWHEG simulation were varied: the b-quark mass, the perturbative scales (µ R,F ) and the PDF.The systematic uncertainties were obtained by taking the largest upward and downward variations in the final yields.In addition, another source of uncertainty as taken into consideration by using the PYTHIA 6 decayer instead of EvtGen to decay the beauty hadrons.The b-hadron feed-down fractions with their systematic uncertainties are shown in Fig. 3  It was verified that after the 3 rd Bayesian iteration the unfolding procedure converges and subsequent iterations do not differ significantly from the previous one.In addition, the prior spectrum used as initial guess was varied in a wide range, by using power-law functions with exponents differing by up to 4 units from each other.Finally, a different unfolding technique based on the Singular Value Decomposition (SVD) method [61] was used, as well as a simple bin-by-bin correction technique.In all of these tests, the deviations from the nominal result were found to be smaller than the statistical uncertainties of the measurement.Therefore, the systematic uncertainties were assigned using a Monte Carlo closure test.

In this test, a detector-level simulation of D 0 -meson tagged jets with a statistical precision comparable to our data was unfolded using a detector response matrix obtained from a different and larger Monte Carlo sample.The unfolded result w s compared with the generator-level spectrum.A p T -independent 5% systematic uncertainty was assigned based on the maximum deviations observed between the unfolded results and the truth.

Finally, he normalization of the p T -differential cross section is affected by uncertainties on the D 0 → K − π + branching ratio (1% [51]) and on the minimum-bias trigger efficiency (3.5% [48]).

The total systematic uncertainties on the cross section were obtained by summing in quadrature the uncertainties estimated for each of the sources outlined above.They rise slightly with increasing p ch T,jet and are comparable to the statistical uncertainties, except for p ch T,jet > 20 GeV/c, where the statistical uncertainty dominates.Similarly, the z ch | | di

ribution
for 5 < p ch T,jet < 15 GeV/c is affected by statistical and systematic uncertainties at a comparable level, while for 15 < p ch T,jet < 30 GeV/c statistical ones dominate.


Results

Figure 5 (left) shows the p T -differential cross section of charm jets containing a D 0 meson in pp collisions at √ s = 7 TeV.The cross section is shown for 5 < p ch T,jet < 30 GeV/c.The D 0 mesons used to tag the j ght) shows the rate of the D 0 -meson tagged jets over the inclusive jet production as a function of p ch T,jet :
R(p ch T,jet ) = N D 0 jet (p ch T,jet ) N jet (p ch T,jet ) .(11)
The inclusive jet production cross section in pp collisions at √ s = 7 TeV was reported by ALICE in [57] and more recently in [60], where the kinematic reach was extended down to p ch T,jet = 5 GeV/c.The rate increases from about 0.04 to about 0.08 in the range 5 < p ch T,jet < 10 GeV/c; it then tends to flatten at a value around 0.08 in the range 8 < p ch T,jet < 20 GeV/c.According to Monte Carlo simulations based on POWHEG + PYTHIA 6, the increase of the charm-jet fraction in the interval 5 < p T < 8 GeV/c is not  The measurements are compared with PYTHIA 6.4.28 [62] (Perugia-2011 tune [49]), PYTHIA 8.2.1 [62] (Monash-2013 tune [63]) and Herwig 7 [64,65] (MEPP2QQ and MEMinBias).Both versions of PYTHIA overestimate the yield by a factor ≈ 1.5 which appears to be approximately constant in the measured p ch T,jet range as shown in Fig. 5 (left).However, since they also overestimate the inclusive jet cross section by a similar amount [60], they provide a good description of the ratio of D 0 -meson tagged jets over the inclusive jet production, as shown in Fig. 5 (right).For the purpose of this comparison the most prominent difference between the Herwig processes MEPP2QQ and MEMinBias is that the former implements massive quarks in the matrix element calculations, whereas the latter treats all quarks as massless partons [64].When calculating the ratio, the inclusive jet cross section from the Herwig MEMinBias process was used for both the MEMinBias and the MEPP2QQ processes in the numerator.Both Herwig processes end to overestimate the measured cross section, with MEPP2QQ describing the data better.The MEMinBias process reproduces well the ratio to the inclusive jet cross section.

The measurement is also compared with two NLO pQCD calculations obtained with the POWHEG-BOX V2 framework [28][29][30], matched with PYTHIA 6 (Perugia-2011 tune) for the generation of the parton shower and of the non-perturbative aspects of the simulation, such as hadronization of colored partons and generation of the underlying event.The theoretical uncertainties were estimated by varying the renormalization and factorization scales (0.5µ 0 ≤ µ F,R ≤ 2.0µ 0 with 0.5 ≤ µ R /µ F ≤ 2.0), the mass of the charm quark (m c = 1.3, 1.7 GeV/c 2 with m c,0 = 1.5 GeV/c 2 ) and the parton distribution function (central points: CT10nlo; variation: MSTW2008nlo68cl [66]).Two process implementations of the POWHEG framework were employed: the heavy-quark [67] and the di-jet implementation [31].As shown in Fig. 6 (left), good agreement is found within the theoretical and experimental uncertainties between the measured p T -differential cross section and the cross section obtained with the POWHEG heavy-quark implementation.The POWHEG di-jet implementation systematically overestimates the production yield by a constant factor of ≈ 1.5.The right panel of Fig. 6 shows a comparison of the ratio of the D 0 -meson tagged jet yield over the inclusive jet yield, for the data and POWHEG.The inclusive jet yield was obtained with the POWHEG di-jet implementation for both presented POWHEG ratio cases.In the calculation of the theoretical systematic uncertainties on the cross-section ratios, the variations of the perturbative scales and of the PDF were applied consistently in the numerator and in the denominator.The measured ratio is found to be in agreement with the di-jet implementation, while the heavy-quark implementation systematically underestimates the ratio.It is worth remarking that the excellent agreement between the data and the POWHEG di-jet implementation for R(p ch T,jet ) means that it overestimates both the D 0 -meson tagged and the inclusive jet cross sections by a similar factor.This peak is populated by jets in which the D 0 meson is the only constituent.In the higher p ch T,jet interval single-constituent jets are much rarer and the peak at z ch | | ≈ 1 disappears.In general, as p ch T,jet increases the fragmentation becomes softer, a feature that has been observed also in inclusive jet measurements [68].

In Fig. 7, the data are compared with simulations obtained with the POWHEG heavy-quark implementation and the erwig 7 MEPP2QQ process, both of which showed the best agreement with the D 0 -meson tagged jet p T -differential cross section in Figs. 5 and 6.

The same data are shown in Fig. 8 with a different normalization choice.The z ch | | -differential cross secti n was divided by the inclusive jet cross section integrated in the corresponding p ch T,jet interval:
R(p ch T,jet , z ch | | ) = N D 0 jet (p ch T,jet , z ch | | ) N jet (p ch T,jet ) . (12)
In this case the data 7 MEMinBias process, which showed the best agreement with the ratio of the D 0 -meson tagged jet cross section over the inclusive jet cross section in Figs. 5 and 6.The choice of these two normalization approaches facilitate the comparison between data and simulation of the shapes of the z ch
| | distributions.
All models show an overall good agreement with the z ch | | -differential data for jets with 5 < p ch T,jet < 15 GeV/c, with the only exception of Herwig 7 MEPP2QQ, which features a substantially harder fragmen-   tation of D 0 mesons in jets.For jets with 15 < p ch T,jet < 30 GeV/c the models can describe the data quite well within large uncertainties.A depletion is observed in the last z ch | | bin in data compared to all models; however the discrepancy is only slightly larger than 1 σ.


The measurement of the z ch

| | distribution for 15 < p ch T,jet < 30 GeV/c partially overlaps with the D * ± in-jet fragmentation data reported by ATLAS in [25] for 25 < p T,jet < 30 GeV/c.While the jet measurement reported here includes only charged tracks, the ATLAS measurement also includes neutral constituents of the jets, and this difference should be taken into account while omparing our measurement with that of ATLAS.The mean transverse momenta of the track-based jets considered in our analysis are 7.53 ± 0.07 GeV/c for 5 < p ch T,jet < 15 GeV/c and 19.5 ± 0.1 GeV/c for 15 < p ch T,jet < 30 GeV/c.Using a POWHEG + PYTHIA 6 simulation, it was estimated the transverse momentum of D 0 -meson tagged jets increases on average by 12% for 5 < p ch T,jet < 15 GeV/c and by 14% for 15 < p ch T,jet < 30 GeV/c when neutral particles are included.Furthermore, ATLAS reported jets reconstructed with a resolution parameter R = 0.6 instead of R = 0.4 used throughout this work.Finally, in the case of the ATLAS measurement the contribution from the b-hadron feed-down was not subtracted.ATLAS observed a large disagreement between data and various Monte Carlo event generators, including PYTHIA 6 and POWHEG di-jet.Our data indicate a much better agreement with the simulations, however experimental uncertainties are large.


Conclusions

The measurement of charm jet production and fragmentation in pp collisions at √ s = 7 TeV, in which charm jets are tagged using fully reconstructed D 0 mesons, was presented in this paper.The D 0meson tagged jet p T -differential cross section was reported in the range 5 < p ch T,jet < 30 GeV/c.The fraction of charged jets containing a D 0 -meson increases with p ch T,jet fro

0.042 ± 0.00
(stat) ± 0.006 (syst) to 0.080 ± 0.009 (stat) ± 0.008 (syst).The cross section of D 0 -meson tagged jets was reported also differentially as a function of the jet momentum fraction carried by the D 0 meson in the direction of the jet axis (z ch | | ) for two ranges of jet transverse momenta, 5 < p ch T,jet < 15 GeV/c and 15 < p ch T,jet < 30 GeV/c in the ranges 0.2 < z ch | | < 1 and 0.4 < z ch | | < 1.0, respectively.The data were compared with PYTHIA 6, PYTHIA 8 and Herwig 7.Both versions of PYTHIA are able to describe reasonably well the ratio to inclusive jets, but not the cross section.The Herwig 7 implementation of the heavy-quark production process can describe both the p T -differential cross section of D 0 -meson tagged jets and its ratio to the inclusive jet cross section.The measurement was also compared with two NLO pQCD calculations obtained with the POWHEG heavy-quark and di-jet implementation.The POWHEG heavy-quark implementation can reproduce the absolute cross section; when comparing with the ratio to inclusive jets (using the POWHEG di-jet implementation for the inclusive jets), it significantly underestimates the data.When the POWHEG di-jet implementation is used both for the charm and the inclusive jet production the agreement to the ratio is restored.All reported models can describe the measured D 0 -meson tagged jet fragmentation within the uncertainties.A small tension between the data and simulations is observed for 15 < p ch T,jet < 30 GeV/c, with the data favouring a slightly softer fragmentation.

The experimental uncertainties are dominated by the limited statistics: the analysis of larger data samples, like those collected by ALICE in pp collisions at √ s = 5 and 13 TeV, may allow for more differential measurements and a more conclusive comparison between data and theoretical expectations.Although the uncertainties are still sizeable, agreement of the measurements with calculations provided y PYTHIA, Herwig and POWHEG indicates that the observables studied in this work are well described by pQCD and they can therefore be exploited to address possible modifications to the charm jet production and internal structure induced by the Quark-Gluon Plasma medium formed in heavy-ion collisions.

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.



show the invariant mass distributions of D 0 -meson candidates in tagged jets in different intervals of p T,D and 5 < p ch T,jet < 30 GeV/c.The D 0 -meson tagged jet candida