University of Groningen
Observation of the Lambda(0)(b) -> Lambda+cK+K-pi(-) decay
Aaij, R.; Adeva, B.; Adinolfi, M.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; De Bruyn, K.;
Dufour, L.; Mulder, M
Published in:
Physics Letters B
DOI:
10.1016/j.physletb.2021.136172
IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from
it. Please check the document version below.
Document Version
Publisher's PDF, also known as Version of record
Publication date:
2021
Link to publication in University of Groningen/UMCG research database
Citation for published version (APA):
Aaij, R., Adeva, B., Adinolfi, M., Ajaltouni, Z., Akar, S., Albrecht, J., Alessio, F., De Bruyn, K., Dufour, L.,
Mulder, M., Onderwater, C. J. G., Pellegrino, A., Tolk, S., van Veghel, M., & LHCb Collaboration (2021).
Observation of the Lambda(0)(b) -> Lambda+cK+K-pi(-) decay. Physics Letters B, 815, [136172].
https://doi.org/10.1016/j.physletb.2021.136172
Copyright
Other than for strictly personal use, it is not permitted to download or to forward/distribute the text or part of it without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license (like Creative Commons).
Take-down policy
If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim.
Downloaded from the University of Groningen/UMCG research database (Pure): http://www.rug.nl/research/portal. For technical reasons the number of authors shown on this cover page is limited to 10 maximum.
Physics Letters B 815 (2021) 136172
Contents lists available atScienceDirect
Physics
Letters
B
www.elsevier.com/locate/physletb
Observation
of
the
Λ
b
0
→ Λ
+
c
K
+
K
−
π
−
decay
.
LHCb
Collaboration
a
r
t
i
c
l
e
i
n
f
o
a
b
s
t
r
a
c
t
Articlehistory:
Received4December2020
Receivedinrevisedform18February2021 Accepted18February2021
Availableonline23February2021 Editor:L.Rolandi
The Λ0b→ Λ+cK+K−
π
− decay is observed for the first time using a data sample of proton-proton collisionsatcentre-of-massenergiesof√s=7 and8TeV collectedbytheLHCb detector,corresponding toanintegratedluminosityof3fb−1.TheratioofbranchingfractionsbetweentheΛ0b→ Λ+cK+K−π
− andtheΛ0b→ Λ+cD−s decaysismeasuredtobeB(Λ0b→ Λ+cK+K−
π
−) B(Λ0b→ Λ+cD−s)= (9.26±0.29±0.46±0.26)×10−2,
wherethefirstuncertaintyisstatistical,thesecondsystematicandthethirdisduetotheknowledgeof the D−s →K+K−
π
− branchingfraction.Nostructureonthe invariantmassdistributionoftheΛ+cK+ systemisfound,consistentwithnoopen-charmpentaquarksignature.©2021TheAuthor(s).PublishedbyElsevierB.V.ThisisanopenaccessarticleundertheCCBYlicense (http://creativecommons.org/licenses/by/4.0/).FundedbySCOAP3.
1. Introduction
Over the last two decades, a wealth of information has been accumulated on the decays of hadrons containing b quarks [1]. Measurementsoftheirdecayratesandpropertieshavebeenused to testthe Cabibbo-Kobayashi-Maskawamechanism [2,3] describ-ing weakinteractionswithintheStandardModel,andtoexamine various theoretical approaches, such asthe heavy quark effective theory [4] and thefactorizationhypothesis [5–8]. Althoughmany
b-hadron decays have been observed with their branching
frac-tions measured, a large number of them remains either unob-served or poorly measured, most notably decays of
Λ
0b,Ξ
b andΩ
b− baryons.Inthelast years,theLHCbexperimenthasobserved many newΛ
0b decays to final states such asΛ
+cπ
−π
+π
− [9],Λ
+cπ
−p p [10],Λ
+c D−s [11],χ
c1p K−,χ
c2p K−[12],ψ(
2S)
p K−andJ
/
ψ
π
+π
−p K−[13].1In thisLetter, the first observationof the
Λ
b0→ Λ
+c K+K−π
− decay (referred to hereafter assignal channel)is reported, along with a measurement of its branching fraction relative to that of theΛ
0b→ Λ
+c D−s decay(normalisationchannel).Theanalysisuses a data sample of proton-proton(pp) collisions atcentre-of-mass energies of√
s=
7 and 8 TeV collectedby the LHCb experiment, corresponding to an integrated luminosity of 3 fb−1. The obser-vation oftheΛ
0b→ Λ
+cK+K−π
− decay provides a laboratory to search for open-charm pentaquarks with valence quark content1 Thecharge-conjugateprocessisimpliedthroughoutthisLetter.
Fig. 1. FeynmandiagramoftheleadingcontributiontotheΛb0→ Λ+cK+K−π− sig-naldecay.
c
¯
suud that could decay strongly to theΛ
+cK+ final state. Thesestatesareanaturalextensionofthethreenarrowpentaquark can-didateswithquark contentccuud observed
¯
inΛ
0b→
J/
ψ
p K− de-cays [14], with the c quark¯
replaced by an¯
s quark. The recent discoveryofa D+K− structureinB−→
D−D+K−decays [15,16], consistentwithopen-charmtetraquarks,alsomotivatesthesearch foropen-charmpentaquarks.Fig.1showstheleadingdiagramcontributingtothesignal de-cay. Contributions to the companion K+K−
π
− system could be through intermediate a−1 mesons, such as the a1(
1260)
− state,whichis found to dominatein B
→
D(∗)K∗0K− decays [17]. De-caysofΣ
c0→ Λ
+cπ
−orevenΞ
c0→ Λ
+cK− could alsocontribute tothesignal.https://doi.org/10.1016/j.physletb.2021.136172
0370-2693/©2021TheAuthor(s).PublishedbyElsevierB.V.ThisisanopenaccessarticleundertheCCBYlicense(http://creativecommons.org/licenses/by/4.0/).Fundedby SCOAP3.
2. Detectorandsimulation
The LHCb detector [18,19] is a single-armforward spectrom-eter covering the pseudorapidity range 2
<
η
<
5, designed for the study ofparticles containing b or c quarks.The detector in-cludesahigh-precisiontrackingsystemconsistingofasilicon-strip vertex (VELO) detector surrounding the pp interaction region, a large-areasilicon-stripdetectorlocatedupstreamofadipole mag-net with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes placed downstream of the magnet. The tracking system provides a measurement of themomentum,p,ofchargedparticleswitharelativeuncertainty that varies from 0.5% at low momentum to 1.0% at 200 GeV/c. The minimum distance of a trackto a primary pp collision ver-tex (PV), the impact parameter (IP), is measured with a reso-lution of(
15+
29/
pT)
μ
m, where pT is the component of themomentum transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished using information from two ring-imagingCherenkovdetectors.Photons,electrons andhadrons are identified by a calorimeter systemconsisting of scintillating-pad and preshowerdetectors,an electromagneticanda hadronic calorimeter(HCAL).Muonsareidentifiedbyasystemcomposedof alternatinglayersofironandmultiwireproportionalchambers.
The onlineeventselectionis performedbya triggerbased on signal information only.The trigger consistsofa hardware stage, basedon informationfromthecalorimetersystem, followedbya software stage, whichapplies a full event reconstruction [20]. At the hardware trigger stage, referred to as L0 trigger in the fol-lowing, the
Λ
0b→ Λ
c+K+K−π
− andΛ
0b→ Λ
+c D−s candidates are required to include a hadron having high transverse energy de-posited in the calorimeters. The transverse energy threshold is 3.5 GeV. Thesoftwaretrigger,alsonamedhigh-level trigger(HLT), requiresa two-,three- orfour-track vertexwithasignificant dis-placement fromany PV. At leastone charged particlemust have a large transverse momentum andbe inconsistent with originat-ing from any PV. A multivariate algorithm [21] is used for the identificationofdisplaced verticesconsistent withthedecayofab-hadron.
Simulation is used to model the effects of the detector ac-ceptanceandthe selection requirements,tovalidate thefit mod-els and to evaluate efficiencies. In the simulation, pp collisions
are generated using Pythia 8 [22] with a specific LHCb con-figuration [23]. Decays of unstable particles are described by EvtGen [24], in which final-state radiation is generated using Photos [25]. The interaction of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [26] asdescribedinRef. [27].
3. Eventselection
Candidate
Λ
+c andD−s hadronsarereconstructedthroughtheir decaystothe p K−π
+andK+K−π
− finalstates,respectively.The offline candidateselection isperformedbyapplyingaloose pres-election,followedbyamultivariateanalysis(MVA)tofurther sup-presscombinatorialbackgroundoriginatingfromrandom combina-tions.Toreducesystematicuncertaintiesontheratioofefficiencies betweenthe signal andthenormalisation channels, theselection criteriaofΛ
+c candidatesareidenticalbetweenthetwochannels.A good-quality track with pT
>
100 MeV/c and p>
1 GeV/c isrequired foreach final-stateparticle. Protonsandantiprotonsare required to have a momentum greater than 10 GeV/c to improve their identification.Allfinal-stateparticlesare alsorequiredtobe inconsistent with originating from any PV by requiring a large
χ
2IP,where
χ
IP2 isdefined asthe difference in theχ
2 of a givenPV fit withandwithout the trackunder consideration. Each
Λ
+c baryon candidateis required to have at leastone decayproductwith pT
>
500 MeV/c and p>
5 GeV/c, a good-quality vertex (i.e.small
χ
2vtx), andinvariant mass within
±
15 MeV/c2 of theknownΛ
+c mass [1]. FortheΛ
+c candidates, thesumof transverse mo-mentaoftheirdecayproductsmustexceed1.8 GeV/c.Theselection criteriafor D−s candidates aresimilar to thoseofΛ
+c candidates. The K+K−π
−invariantmassisrequiredtobewithin±
35 MeV/c2fromtheknown D−s mesonmass.
Thesignalchannelisreconstructedbycombining
Λ
+c,K+, K−and
π
− candidates, while the normalisation channel is recon-structed by combining aΛ
+c with a D−s candidate. The combi-nationsabove formΛ
0b candidates, whichare required to havea smallχ
2vtx and
χ
IP2, and a decay time with respect to itsasso-ciated PV greater than 0
.
2 ps. The associated PV is the one that givesthe smallestχ
2IP, wherethe
χ
IP2 denotes the IP significanceofcandidate’s trajectory returned by the kinematicalfit. The an-gle between the
Λ
b0 momentum and the vector pointing from the associated PV to theΛ
0b decay vertex,θ
p, is required to be smallerthan 11mrad.TheΛ
0b candidateis alsorequiredtohave atleastonefinal-state particlewith pT>
1.
7 GeV/c, andits decayvertexsignificantly displaced fromanyPV. The latteris achieved byrequiringthesignificanceoftheflightdistancebetweenthe
Λ
0b decayvertexandanyPVtobe largerthan4.Final-statetracksof signal andnormalisationcandidates mustpass stringent particle-identification requirements based on the information from RICH detectors,calorimetersystemand muonstations. Toreject tracks thatsharethesamesegmentintheVELOdetector,anytwotracks withthesamechargeusedtoformtheΛ
0b candidatearerequired tohaveanopeninganglelargerthan0.5mrad.Akinematicfit [28] ofthedecaychain constrains theΛ
0b candidateto originate from the associated PV and theΛ
+c candidate invariant mass to its knownvalue [1].The
Λ
0b candidatecould originate from B0→
D+K+K−π
− orB0s
→
D+sK+K−π
− decays, where a pion or kaon inD+
→
K+π
−π
+ or D+s(
D+)
→
K+K−π
+ decays is misidenti-fiedasaproton.Thesebackgroundcontributionsarevetoedifthe invariant massesof theΛ
+c andΛ
b0 candidates,evaluated by re-placingtheprotonbyeitherthepionorkaonmasshypothesis,are within±
15 MeV/c2oftheknownD+(
D+s
)
massand±
25 MeV/c2 of theknown B0(
B0s
)
mass [1].Thesevetoesare appliedtoboththe signal and the normalisation channels. For the signal decay, ad-ditionalvetoesare appliedifthe invariant massofthe K+π
− orK+K−
π
− companiontracksfallswithin±
30 MeV/c2 ofthe D0 orD−s knownmass,respectively [1].
Reconstructedcandidates arefurther requiredtopassan MVA output threshold based upon a multilayer perceptron (MLP) fil-ter [29], designed to reject the combinatorial background. The MLP classifier is trained using a signal sample of simulated
Λ
0b→ Λ
+cK+K−π
− decays tuned on datato reproduce correctly theΛ
0b production kinematics based on the pT and ydistribu-tionsandabackgroundsample takenfromtheuppersidebandof the
Λ
b0 invariant mass spectrum in the range of 5.75 – 7 GeV/c2.Afour-bodyphase-spacesimulationisusedforthesignal sample tokeeptheMLPefficiencyasuniformaspossible,asincluding in-termediate resonancesinthe simulationcould potentially leadto small MLP efficiencies for less represented phase-space regions. The lower sideband is not used to avoid potential background contributions frompartially reconstructed decays.The MLP input includes thefollowing variables: pT sumof the
Λ
+c decay prod-ucts,minimalχ
2IPamongthe
Λ
c+decayproducts,minimal pT andminimal
χ
2IPamongthekaonsoriginatingdirectlyfromthe
Λ
b0 de-cay,pTandχ
IP2 oftheπ
−fromtheΛ
0bdecay,pT sumofallπ
and K originating directlyfrom theΛ
0b decay,χ
2vtx of the
Λ
+c candi-date,χ
2 ofthe flightdistancebetweentheΛ
0b decayvertexand theassociatedPV,cos
θ
p,χ
2 probabilityoftheΛ
0ver-LHCb Collaboration Physics Letters B 815 (2021) 136172
Fig. 2. Invariantmassdistributionof(a)Λ0
b→ Λ+cK+K−π− and(b)Λ0b→ Λ+cD−s candidates.Fitprojectionsareoverlaidasabluesolidline.For(a),theredsolidline representsthesignalcomponent,thebluedashedlineisthebackgroundduetorandomcombinations,andthevioletdottedlineisthecontributionfromΛ0
b→ Σc+K+K−π− decays.For(b),theredsolidlineisthenormalizationchannelcomponent,thevioletdottedlineisthe Λ0b→ Λ+cD∗−s background,thegreendashed-dottedlineisthe contributionfromΛ0b→ Λ+cK+K−π−decays,andthebluedashedlinerepresentscombinatorialbackground.
texfit,andthedifferenceoflongitudinalpositionbetweenthe
Λ
+c andtheΛ
0b decayvertices.The MLP response obtained from the training is also applied to the normalisation channel sample. The optimal thresholds on the MLP response are obtained for the signal and normalisation channels separately by maximising a figure-of-merit, defined as
S
/
√
S+
B,whereS andB aretheexpectedsignalandbackground yieldsforΛ
0bcandidateswithina±
2.
5σ
masswindowaroundthe knownΛ
0bmass [1],whereσ
isthemassresolutioncorresponding to about12 MeV/c2.Both S and B aredetermined bymultiplyingtheinitialyieldsofsignal andbackgroundwiththecorresponding MLPselectionefficienciesestimatedfromsimulationandsideband data,respectively.Theinitialsignalandbackgroundyieldsare ob-tainedfromapreliminaryfittothepreselecteddatasamplebefore theMLPrequirementapplied,wherethesignal
Λ
b0peakisalready seenintheΛ
+cK+K−π
− invariantmassdistribution.Theoptimal point corresponds toa signalefficiency of90%anda background rejectionof 85%.About 0.6%eventsinthe signal channelcontain multiple candidates, onlyone candidate is retainedby a random selection.4. Signalyieldsandsearchforintermediatestates
The yields in both the signal andnormalisation channels are determined from an unbinned extended maximum-likelihood fit to the corresponding invariant mass spectraof the
Λ
+c K+K−π
− system. The signal component is modelled by a sum of two Crystal Ball functions [30] with a common mean of the Gaus-sian cores, with tail parameters fixed to the values obtained from simulation. For both the signal and normalisation chan-nels, the combinatorial background is described by an expo-nential function, whose parameters are varied freely and al-lowed to be different between the signal and normalisation channels. For the signal channel, a significant contribution fromΛ
0b→ Σ
c+[→ Λ
+cπ
0]
K+K−π
− decays ispresentin thelower in-variant mass region, which has the same final state as theπ
0is not reconstructed. The shape of this background is obtained from a simulation of
Λ
b0→ Σ
c(
2455)
+K+K−π
− decays. For the normalisation channel, theΛ
0b→ Λ
+c D∗−s decay may be recon-structed asΛ
0b→ Λ
+cDs− due to photon emission in the D∗−sdecay. The shape of thisbackground is obtained from simulated
Λ
0b→ Λ
+cDs∗− decays. The signal decay can also contribute to the normalisation channel forming a background under the D−smass peak. This background contribution is estimated from the
D−s sidebandsofthenormalisationdatasample,wherethewidth of the sideband is chosen to be the same as that of the D−s
mass window used in the normalisation channel selection. The
Fig. 3. InvariantmassdistributionsofΛ+cK+candidatesintheΛ0b→ Λ+cK+K−π− signalchannelforthe simulation(red line)andthe background-subtracteddata (bluepointswitherrorbars).
invariantmassdistributionsforthesignalandnormalisation chan-nelsare showninFig.2 withthefitprojects overlaid.The signal yieldsareobtainedtobeN
(Λ
0b→ Λ
+c K+K−π
−)
=
3400±
80 andN
(Λ
0b→ Λ
+cD−s[
K+K−π
−])
=
2550±
60, respectively,wherethe uncertaintiesarestatisticalonly.Anopen-charmpentaquarkstate couldbe revealedasa struc-ture in the invariant mass distribution of the
Λ
+c K+ system, shownin Fig.3 fordata andsimulation. The datadistribution is backgroundsubtractedthroughthesPlot weightingtechnique [31], usingtheΛ
+c K+K−π
− invariantmassasdiscriminating variable. Nostructureisobserved.Afullamplitudeanalysisisneededto es-timatethe limitof thepentaquark contribution,whichis beyond thescopeofthisLetter.Instead,arich structureofknownhadroncontributions is vis-ible in thebackground-subtracted invariant mass distributions of the
Λ
+cπ
−, K+π
− and K+K−π
− systems, shown inFig. 4. TheΣ
c(
2455)
0andΣ
c(
2520)
0resonancesarevisibleintheΛ
+cπ
− dis-tribution. A large K∗(
892)
0 resonance is observed in the K+π
−projection.In the K+K−
π
− system,a broadpeakingstructure at about1.
5 GeV/c2 isalsoobserved.Asimilar structureis alsoseen inB→
D(∗)K∗0K−decaysbytheBelleexperiment [17],andis ex-plainedasthetailcontributionofthea1(
1260)
−resonance.5. Branchingfractionratioandefficiencies
Theratioofthebranchingfractionsofthe
Λ
b0→ Λ
+cK+K−π
− decayincludingresonance contributionswith respectto the nor-malisationchannelisdeterminedbyFig. 4. Invariantmassdistributionsof(a)theΛ+cπ−,(b)K+π−,and(c)K+K−π− systemsintheΛ0b→ Λ+cK+K−π−signalchannel,forthebackground-subtracteddata. ThereddashedverticallinesindicatethevetomassintervalsforD0mesonsintheK+π−distribution,andD−
s intheK+K−π−distribution.
B
(Λ
0b→ Λ
+c K+K−π
−)
B
(Λ
b0→ Λ
+c D−s)
=
N(Λ
0b→ Λ
+c K+K−π
−)
N(Λ
b0→ Λ
+c Ds−[
K+K−π
−])
×
(1)tot
(Λ
0b→ Λ
+c D−s[
K+K−π
−])
tot
(Λ
0b→ Λ
+c K+K−π
−)
×
B
(
D−s→
K+K−π
−),
where
B
standsforthebranchingfractionofthecorresponding de-cay.ThesignalandnormalisationyieldsarereportedinSec.4.The totalefficienciestot ofthesignalandthenormalisationchannels
aredeterminedbytheproduct
tot
=
acc
×
sel
×
L0
×
HLT
×
PID
,
(2)where
acc accountsfor the LHCb geometrical acceptance,
sel is
the efficiency of reconstructing and selecting a candidate within the acceptance,
L0 is the L0 trigger efficiency for the selected
candidates,
HLT is the HLT efficiency for the selected
candi-dates passingtheL0triggerrequirement,and
PID isthe
particle-identification(PID)efficiencyfortheselectedcandidatesthat sur-vive all trigger requirements. All efficiencies except for
L0 and
PIDaredeterminedfromsimulation,andthe
(
pT,
y)
distributionsofthe simulated
Λ
0b baryonsareweighted tomatchthat ofdata, where y istherapidityofthecandidate.Theweightsareobtained usingthenormalisationchannelandappliedtothesignaldecay.To take into account the resonance contributions to the sig-nal decay channel, the simulation uses a mixture of three de-cay modes:
Λ
b0→ Λ
+ca1(
1260)
−(
→
K∗0K−)
,Λ
+c K∗0K− and non-resonant four-body phase space. The fractions are determined by fitting the two-dimensional data distribution of K+π
− andK+K−
π
− invariantmasses.The L0 efficiency of each hadron is computed using samples of well identified pions and kaons from D0
→
K−π
+ decays and protons fromΛ
→
pπ− decays [32]. The efficiency is cal-culated in bins of transverse energy for the particles incident on the HCAL surface, separately for its inner and outer regions. The PID efficiency is determined by the calibration samples ofD∗+
→
D0(
→
K−π
+)
π
+ andΛ
→
pπ− decays andis evaluated asafunctionoftrackmomentum,trackpseudorapidityandevent multiplicity,wherethelatterisrepresentedbythenumberofthe reconstructedtracksintheevent.Theratiobetweenthetotalefficiencies forthesignal and nor-malisation channels in Eq. (1), is determined to be 0
.
78±
0.
02, where the uncertainty accounts only for the size of the simula-tionsample.Thevaluediffersfromunityprimarilyduetodifferent selectionefficienciesontheMVAresponsesforthesignaland nor-malisationchannels.External inputs are used for the branching fractions
B(
D−s→
K+K−π
−)
= (
5.
39±
0.
15)
×
10−2 [1] andB(Λ
0b
→ Λ
+c D−s)
= (
1.
10±
0.
10)
×
10−2 [11]. In the latter case, whilethevalueismeasuredbytheLHCbcollaboration [11],its un-certainty isdominated by thebranching fractionof B0→
D+D−s decays,andisessentially uncorrelatedwiththepresent measure-ment.
6. Systematicuncertainties
Allsystematicuncertaintiesonthemeasurementoftheratioof branching fractions are listed inTable 1. The total uncertaintyis determined fromthe sumofall contributions in quadrature.The dominantuncertaintyisrelatedtotheresonancestructurethatis notperfectlymodelled bythesimulation.
Uncertaintiesduetothefitmodelareconsidered.Forthe back-groundduetorandomcombinationsoffinal-stateparticlesinboth
LHCb Collaboration Physics Letters B 815 (2021) 136172
Table 1
Summaryofsystematicuncertaintiesontheratioofbranchingfractions.
Source Uncertainty (%)
Combinatorial background 0.9
Shape ofΛ0
b→ Σc+K+K−π−contribution 0.3
Λ0
b→ Λ+cK+K−π−background in normalisation channel 0.8
Signal fit model 0.5
Simulation sample size 2.5
PID efficiency 0.4
Trigger efficiency 0.1
(pT, y)weight 0.8
Track multiplicity weight 0.8
Λ+c Dalitz structure 1.4
Mixture fraction in simulation 0.2
Resonance structure 3.6
Multiple candidates 0.3
MVA selection 0.5
Total 4.9
thesignalandnormalisationchannels, theexponentialfunctionis replaced by a second-order polynomial function. From the com-parison tothedefaultresult, therelativeuncertaintyon theratio of branching fractions is 0.9%. In the signal channel, the uncer-tainty duetothe
Λ
b0→ Σ
c+K+K−π
− background contributionis assessedbyperformingthefitwithawidenedmassregion, result-inginarelativeuncertaintyof0.3%.Forthenormalisationchannel, changingtheyieldoftheΛ
0b→ Λ
+cK+K−π
− contributionwithin itsuncertaintyresultsinarelative0.
8% variation.Thesystematicuncertaintyduetothemodelforbothsignaland normalisationchannels,isstudiedbychangingtoasingleHypatia function [33],wherethemeanandwidthparametersareleftfree whileallother parametersaretakenfromsimulation.Thisresults inarelativeuncertaintyof0
.
5%.Theuncertaintiesontheratioofefficienciesareevaluated.The uncertainty dueto the finite simulationsample size isevaluated fromtheexpectedefficiencyvariation inbins of pT and y of the
Λ
0bcandidateasσ
=
ii
(
1−
i
)
Niwi i Niwi,
(3)for each bin i, where Ni isthe number of generated events, wi isacorrectionweight,and
i isthecandidateefficiency.The nor-malisationof theweightsis chosensuch that thedenominator is equaltototalnumberofgeneratedeventswithouttheweighting. Therelativeuncertaintyisfoundtobe2
.
5%.Pseudoexperiments areused toevaluate thesystematiceffects due to uncertainties on the weights or efficiencies in different bins. Foragiven source,manypseudoexperimentsare generated, in which each produces a new set ofweights or efficiencies ac-cordingtothecentralvaluesanduncertaintiesfollowingGaussian distributions.Theefficiencyratiobetweenthesignaland normali-sationchannelsisrecomputed.Theresultingefficiencyratiosfrom manypseudoexperimentsofthissourceproduceaGaussian distri-bution centering atthe baseline value. The standard deviationof the Gaussian distributionis takenasabsolute uncertaintyon the efficiency ratiofor the givensource. The procedure is applied to obtainthesystematicuncertaintyrelatedtothePIDandtrigger ef-ficienciesandto
(
pT,
y)
andtrackmultiplicityweighting.Thetrackingefficiencyreturnedby thesimulationiscalibrated using adata-driven method [34].The uncertainty onthe calibra-tionsamplesizeispropagatedtotheefficiencyratiousing pseudo-experiments,resultinginasystematicuncertaintyof0.8%.Because the final states forsignal and normalisation modesare identical,
possible data-simulation differences in hadron interactions with thedetectormaterialareassumedtobenegligible.
The agreement between data and simulation for the
Λ
+c→
p K−
π
+ channel is tested by comparing the Dalitz structure. The signal simulation sample is weighted in the m(
p K−)
versusm
(
K−π
+)
plane to match the distribution of the background-subtracteddata.Theuncertaintyrelatedtothelimitedsamplesize used for obtaining these weights is 1.1%, obtained from pseudo-experiments. The uncertainty related to the choice of binning is 0.8%,determinedbyusinganalternativebinning.Atotalof1.4%is assignedassystematicuncertainty.Thecontributionsofthe
Λ
0b decaysthroughthemixtureofthe threedecaymodesareconsideredwhengeneratingthesimulated eventsofthesignalchannel,andtheirfractionsareobtainedby fit-tingthetwo-dimensionaldistributionoftheK+K−π
−andK+π
−systemsinthebackground-subtractedsignaldata.Thefractionsare changed according to the statistical uncertainty of the fit result, yielding0.2%ofrelativeuncertainty.
The simulation does not fully model the resonance structure,
e.g. the contribution of
Σ
c0 resonances, which is clearly seen inthe
Λ
+cπ
− invariant mass distribution, as illustrated in Fig. 4. Byweighting thesimulation tomatchthem(Λ
+cπ
−)
distribution in the data, a 1.3% variation of the ratio of branching fractions is found and assigned as systematic uncertainty. Besides, differ-ences between background-subtracted data andsimulated signal eventsarealsoobservedintheinvariantmassdistributionsoftheΛ
+c K+K− and K+K− systems. To account for this discrepancy, the simulatedsample is weighted according to theΛ
+c K+K− orK+K− mass distribution of background-subtracted data, andthe ratioofbranchingfractionsisreevaluated.Thetwoprocedures re-turnchanges of 2.6%and 2.0%,respectively. The three valuesare added inquadrature to account for the uncertainty dueto reso-nancestructure.
Simulation does not account well for multiple candidates, whichisfoundto beabout0.6%ofthe datasampleinthe signal channel.Halfofthisfractionisassignedassystematicuncertainty duetotherandomchoicetoretainonlyonecandidate.
TheMVAselectioncriteriaareoptimizedseparatelyforthe sig-nalandnormalisationchannels.Asanalternativechoice,theMVA selectionofthe normalisationchannel isfixed tobe thesameas thatofthesignalchanneltotesttherobustnessoftheMVA selec-tion.Therelative variationofthe branchingfractionratiois0.5%, whichisassignedassystematicuncertainty.
7. Resultsandsummary
The first observation of the
Λ
0b→ Λ
c+K+K−π
− decay is presented, and the branching fraction is determined using theΛ
0b→ Λ
+cD−s decay as a normalisation channel. The relative branchingfractionismeasuredtobeB
(Λ
0b→ Λ
+c K+K−π
−)
B
(Λ
0b→ Λ
+c D−s)
= (
9.
26±
0.
29±
0.
46±
0.
26)
×
10−2,
wherethefirstuncertaintyisstatistical,thesecondsystematic,and thethirdisduetotheknowledgeofthe D−s
→
K+K−π
− branch-ingfraction [1].Usingthisratio,theΛ
b0→ Λ
+c K+K−π
−branching fractionisdeterminedtobeB
(Λ
b0→ Λ
+c K+K−π
−)
= (
1.
02±
0.
03±
0.
05±
0.
10)
×
10−3,
where the third term includes the uncertainty on the branching fractionoftheΛ
0b→ Λ
+c D−s decay [1]. Theinvariant mass distri-butionoftheΛ
+cK+systemisinspectedforpossiblestructuredue toopen-charmpentaquarks,andnocontributionisobserved.Declarationofcompetinginterest
Theauthorsdeclarethattheyhavenoknowncompeting finan-cialinterestsorpersonalrelationshipsthatcouldhaveappearedto influencetheworkreportedinthispaper.
Acknowledgements
We express our gratitude to our colleagues in the CERN ac-celerator departments for the excellent performance of the LHC. We thankthe technicaland administrative staff at the LHCb in-stitutes. We acknowledge support from CERN and from the na-tional agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); MOST and NSFC (China); CNRS/IN2P3 (France); BMBF, DFG and MPG (Germany); INFN (Italy); NWO (Netherlands); MNiSW and NCN (Poland); MEN/IFA (Romania); MSHE (Russia); MICINN (Spain); SNSF andSER (Switzerland); NASU(Ukraine); STFC (United King-dom); DOE NP and NSF (USA). We acknowledge the computing resourcesthatareprovidedbyCERN,IN2P3(France),KITandDESY (Germany), INFN (Italy), SURF (Netherlands), PIC (Spain), GridPP (United Kingdom),RRCKIandYandexLLC(Russia), CSCS (Switzer-land),IFIN-HH(Romania),CBPF(Brazil),PL-GRID(Poland)andOSC (USA). We are indebtedto the communities behind the multiple open-source software packages on which we depend. Individual groups or members havereceived support fromAvH Foundation (Germany); EPLANET, Marie Skłodowska-Curie Actions and ERC (European Union); A*MIDEX, ANR, Labex P2IO and OCEVU, and Région Auvergne-Rhône-Alpes (France); Key Research Program of FrontierSciencesofCAS,CASPIFI,ThousandTalentsProgram,and Sci. & Tech. Program of Guangzhou (China); RFBR, RSF and Yan-dex LLC (Russia); GVA,XuntaGal and GENCAT (Spain); The Royal Society andtheLeverhulmeTrust(UnitedKingdom).
Appendix A. Supplementarymaterial
Supplementarymaterialrelatedtothisarticlecanbefound on-lineathttps://doi.org/10.1016/j.physletb.2021.136172.
References
[1]ParticleDataGroup,P.A.Zyla,etal.,Reviewofparticlephysics,tobepublished inProg.Theor.Exp.Phys.6(2020)083C01.
[2]N.Cabibbo,Unitarysymmetryandleptonicdecays,Phys.Rev.Lett.10(1963) 531.
[3]M.Kobayashi,T.Maskawa,CP-violationintherenormalizabletheoryofweak interaction,Prog.Theor.Phys.49(1973)652.
[4]E.Eichten,B.Hill,Aneffectivefieldtheoryforthecalculationofmatrix ele-mentsinvolvingheavyquarks,Phys.Lett.B234(1990)511.
[5]D.Fakirov,B.Stech,FandDdecays,Nucl.Phys.B133(1978)315.
[6]N.Cabibbo,L.Maiani,Two-bodydecaysofcharmedmesons,Phys.Lett.B73 (1978)418,Erratum:Phys.Lett.B76(1978)663.
[7]T.Mannel,W.Roberts,Z.Ryzak,Factorizationhypothesisandthenon-leptonic decaysofheavyhadrons,Phys.Lett.B259(1991)359.
[8]M.Beneke,G.Buchalla,M.Neubert,C.T.Sachrajda,QCDfactorizationfor exclu-sivenon-leptonicb-mesondecays:generalargumentsandthecaseofheavy– lightfinalstates,Nucl.Phys.B591(2000)313.
[9]LHCbcollaboration,R.Aaij,etal.,Measurementsofthebranchingfractionsfor B0
(s)→D(s)π π πandΛb0→ Λ+cπ π π,Phys.Rev.D84(2011)092001,Erratum:
Phys.Rev.D85(2012)039904,arXiv:1109.6831.
[10]LHCbcollaboration,R.Aaij,etal.,ObservationofthedecayΛ0
b→ Λ+cp pπ−,
Phys.Lett.B784(2018)101,arXiv:1804.09617.
[11]LHCbcollaboration,R.Aaij,etal.,Studyofbeautyhadrondecaysintopairsof charmhadrons,Phys.Rev.Lett.112(2014)202001,arXiv:1403.3606. [12]LHCbcollaboration,R.Aaij,etal.,ObservationofthedecaysΛ0b→χc1p K−and
Λ0
b→χc2p K−,Phys.Rev.Lett.119(2017)062001,arXiv:1704.07900.
[13]LHCbcollaboration,R.Aaij,etal.,ObservationofΛ0
b→ ψ(2S)p K− andΛ0b→
J/ψπ+π−p K− decays andameasurementoftheΛb0baryonmass,J.High
EnergyPhys.05(2016)132,arXiv:1603.06961.
[14]LHCbcollaboration,R.Aaij,etal.,Observationofanarrowpentaquarkstate, Pc(4312)+,andoftwo-peakstructureofthe Pc(4450)+,Phys.Rev.Lett.122
(2019)222001,arXiv:1904.03947.
[15]LHCb collaboration,R. Aaij,et al., Model-independentstudyofstructure in B+p
→
D+D−K+ decays, Phys. Rev.Lett. 125 (2020) 242001, arXiv:2009. 00025.[16]LHCbcollaboration,R.Aaij,etal.,AmplitudeanalysisoftheB+→D+D−K+ decay,Phys.Rev.D102(2020)112003,arXiv:2009.00026.
[17]Bellecollaboration,A.Drutskoy,etal.,ObservationofB→D(∗)K−K0(∗)decays,
Phys.Lett.B542(2002)171,arXiv:hep-ex/0207041.
[18]LHCbcollaboration,A.A.AlvesJr.,etal.,TheLHCbdetectorattheLHC,J. In-strum.3(2008)S08005.
[19]LHCbcollaboration,R.Aaij,etal.,LHCbdetectorperformance,Int.J.Mod.Phys. A30(2015)1530022,arXiv:1412.6352.
[20]R.Aaij,etal.,TheLHCbtriggeranditsperformancein2011,J.Instrum.8(2013) P04022,arXiv:1211.3055.
[21]V.V.Gligorov,M.Williams,Efficient,reliableandfasthigh-leveltriggeringusing abonsaiboosteddecisiontree,J.Instrum.8(2013)P02013,arXiv:1210.6861. [22]T.Sjöstrand,S.Mrenna,P.Skands,AbriefintroductiontoPYTHIA8.1,Comput.
Phys.Commun.178(2008)852,arXiv:0710.3820.
[23]I.Belyaev,etal.,HandlingofthegenerationofprimaryeventsinGauss,the LHCbsimulationframework,J.Phys.Conf.Ser.331(2011)032047.
[24]D.J.Lange,TheEvtGenparticledecaysimulationpackage,Nucl.Instrum. Meth-odsA462(2001)152.
[25]P.Golonka,Z.Was,PHOTOSMonteCarlo:aprecisiontoolforQEDcorrections inZ andW decays,Eur.Phys.J.C45(2006)97,arXiv:hep-ph/0506026. [26]Geant4collaboration,J.Allison,etal.,Geant4developmentsandapplications,
IEEETrans.Nucl.Sci.53(2006)270;
Geant4collaboration,S.Agostinelli,etal.,Asimulationtoolkit,Nucl.Instrum. MethodsA506(2003)250.
[27]M.Clemencic,etal.,TheLHCbsimulationapplication,Gauss:design,evolution andexperience,J.Phys.Conf.Ser.331(2011)032023.
[28]W.D.Hulsbergen,DecaychainfittingwithaKalmanfilter,Nucl.Instrum. Meth-odsA552(2005)566,arXiv:physics/0503191.
[29]A.Hoecker,etal.,TMVA4—toolkitformultivariatedataanalysiswithROOT. Usersguide,arXiv:physics/0703039.
[30]T.Skwarnicki,Astudyoftheradiativecascadetransitionsbetweenthe Upsilon-primeandUpsilonresonances,PhDthesisInstituteofNuclearPhysics,Krakow, 1986,DESY-F31-86-02.
[31]M.Pivk,F.R.LeDiberder,sPlot:astatisticaltooltounfolddatadistributions, Nucl.Instrum.MethodsA555(2005)356,arXiv:physics/0402083.
[32]A. Martin Sanchez, P. Robbe, M.-H. Schune, Performances ofthe LHCb L0 CalorimeterTrigger,LHCb-PUB-2011-026,CERN-LHCb-PUB-2011-026,2012. [33]D. MartínezSantos,F.Dupertuis,Mass distributionsmarginalizedover
per-eventerrors,Nucl.Instrum.MethodsA764(2014)150,arXiv:1312.5000. [34]LHCbcollaboration,R.Aaij,etal.,Measurementofthetrackreconstruction
ef-ficiencyatLHCb,J.Instrum.10(2015)P02007,arXiv:1408.1251.
LHCbCollaboration
R. Aaij
31,
C. Abellán Beteta
49,
T. Ackernley
59,
B. Adeva
45,
M. Adinolfi
53,
H. Afsharnia
9,
C.A. Aidala
84,
S. Aiola
25,
Z. Ajaltouni
9,
S. Akar
64,
J. Albrecht
14,
F. Alessio
47,
M. Alexander
58,
A. Alfonso Albero
44,
Z. Aliouche
61,
G. Alkhazov
37,
P. Alvarez Cartelle
47,
S. Amato
2,
Y. Amhis
11,
L. An
21,
L. Anderlini
21,
A. Andreianov
37,
M. Andreotti
20,
F. Archilli
16,
A. Artamonov
43,
M. Artuso
67,
K. Arzymatov
41,
E. Aslanides
10,
M. Atzeni
49,
B. Audurier
11,
S. Bachmann
16,
M. Bachmayer
48,
J.J. Back
55,
S. Baker
60,
P. Baladron Rodriguez
45,
V. Balagura
11,
W. Baldini
20,
J. Baptista Leite
1,
R.J. Barlow
61,
S. Barsuk
11,
W. Barter
60,
M. Bartolini
23,
i,
F. Baryshnikov
80,
J.M. Basels
13,
G. Bassi
28,
B. Batsukh
67,
A. Battig
14,
A. Bay
48,
M. Becker
14,
F. Bedeschi
28,
I. Bediaga
1,
A. Beiter
67,
V. Belavin
41,
S. Belin
26,
V. Bellee
48,
LHCb Collaboration Physics Letters B 815 (2021) 136172