A search for Ξcc++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\Xi}_{\mathrm{cc}}^{++} $$\end{document}→ D+pK−π+ decays

A search for the Ξcc++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\Xi}_{cc}^{++} $$\end{document} baryon through the Ξcc++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\Xi}_{cc}^{++} $$\end{document}→ D+pK−π+ decay is performed with a data sample corresponding to an integrated luminosity of 1.7 fb−1 recorded by the LHCb experiment in pp collisions at a centre-of-mass energy of 13 TeV. No significant signal is observed in the mass range from the kinematic threshold of the decay to 3800 MeV/c2. An upper limit is set on the ratio of branching fractions ℛ=ℬΞcc++→D+pK−π+ℬΞcc++→Λc+K−π+π+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathrm{\mathcal{R}}=\frac{\mathrm{\mathcal{B}}\left({\Xi}_{cc}^{++}\to {D}^{+}{pK}^{-}{\pi}^{+}\right)}{\mathrm{\mathcal{B}}\left({\Xi}_{cc}^{++}\to {\Lambda}_c^{+}{K}^{-}{\pi}^{+}{\pi}^{+}\right)} $$\end{document} with ℛ < 1.7 (2.1) × 10−2 at the 90% (95%) confidence level at the known mass of the Ξcc++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\Xi}_{cc}^{++} $$\end{document} state.


Introduction
The first observed doubly charged and doubly charmed baryon was the Ξ ++ cc (ccu) state found through the Ξ ++ cc → Λ + c K − π + π + and Ξ ++ cc → Ξ + c π + decay modes by the LHCb collaboration [1,2]. With two heavy constituent quarks, this baryon provides a unique system for testing quantum chromodynamics. The average mass of the Ξ ++ cc baryon from the two LHCb measurements now stands at 3621.24 ± 0.65 (stat) ± 0.31 (syst) MeV/c 2 and its lifetime is 0.256 +0.024 −0.022 (stat) ± 0.014 (syst) ps [3], consistent with a weakly decaying state. However, many features of the Ξ ++ cc baryon remain unknown, including its spin and parity. Previously, signals of the singly charged Ξ + cc state were reported in the Λ + c K − π + and pD + K − final states by the SELEX collaboration [4,5]. The masses of the Ξ ++ cc and Ξ + cc ground states are expected to be approximately equal according to isospin symmetry [6]. Searches in different production environments at the FOCUS, BaBar, Belle and LHCb experiments have however not shown evidence for a Ξ + cc state with the properties reported by the SELEX collaboration [7][8][9][10].
To further understand the dynamics of weakly decaying doubly heavy baryons, it is of prime importance to pursue searches for additional decay modes of the Ξ ++ cc baryon. These decays may differ significantly from those of singly heavy hadrons due to interference effects between decay amplitudes of the two heavy quarks. From an experimental viewpoint, the decay Ξ ++ cc → D + (→ K − π + π + )pK − π + is a suitable search channel, since the  Figure 1. The Feynman diagram contributing to the inclusive (left) Ξ ++ cc → D + pK − π + decay with the analogous (right) Ξ ++ cc → Λ + c K − π + π + diagram. D + → K − π + π + trigger is proven to be very efficient at LHCb [11]. 1 The tree-level amplitudes of the inclusive decays of Ξ ++ cc → D + pK − π + and Ξ ++ cc → Λ + c K − π + π + , as shown in figure 1, are comparable, which suggests that the branching fractions of these two modes could be similar. Theoretical calculations have been performed on pseudo-two-body decays of doubly-charmed baryons [12]. The Ξ ++ cc → D + pK − π + decay could proceed as a pseudo-two-body decay if it decays via an excited Σ + * state with a mass greater than 1572 MeV/c 2 , which would then decay to a pK − π + final state. However, the properties of such Σ + * decays are not well known [13]. The Ξ ++ cc → D + pK − π + decay also has a energy release of 180 MeV, compared to 560 MeV for the Ξ ++ cc → Λ + c K − π + π + decay, which means it is expected to have a lower branching fraction because of the smaller available phase space.
The analysis presented in this paper searches for the Ξ ++ cc baryon, at its known mass, through Ξ ++ cc → D + pK − π + decays and also explores a larger mass range to identify the hypothetical isospin partner of the Ξ + cc state that the SELEX collaboration reported. The analysis uses pp collision data corresponding to an integrated luminosity of 1.7 fb −1 recorded by the LHCb experiment in 2016 at a centre-of-mass energy of 13 TeV. The branching fraction of the Ξ ++ cc → D + pK − π + decay is normalised to Ξ ++ cc → Λ + c K − π + π + to reduce systematic uncertainties.
For convenience, the single-event sensitivity, α s , is defined as such that eq. (1.1) reduces to R = α s × N (D + pK − π + ). All aspects of the analysis are fixed before the data in the [3300, 3800] MeV/c 2 mass region are examined.

Detector and software
The LHCb detector [14,15] is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5, designed for the study of particles containing b or c quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region [16], a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes [17] placed downstream of the magnet. The tracking system provides a measurement of the momentum, p, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200 GeV/c. The minimum distance of a track to a primary vertex (PV), the impact parameter (IP), is measured with a resolution of (15 + 29/p T ) µm, where p T is the component of the momentum transverse to the beam, in GeV/c. Charged hadrons are identified using two ring-imaging Cherenkov detectors [18]. Photon, electron and hadron candidates are identified by a calorimeter system consisting of scintillating-pad, pre-shower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers [19]. The trigger consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction. The online reconstruction incorporates near-real-time alignment and calibration of the detector [11]. The same alignment and calibration information is propagated to the offline reconstruction, ensuring consistent and high-quality information between the trigger and offline software. The identical performance of the online and offline reconstruction offers the opportunity to perform physics analyses directly using candidates reconstructed in the trigger [20]. The analysis described in this paper makes use of these features. Simulated Ξ ++ cc → D + pK − π + decays are used to design the candidate selection and to calculate the efficiency of such a selection. The proton-proton interactions are generated using Pythia [21] with a specific LHCb configuration [23]. Genxicc v2.0 [24], the dedicated generator for doubly-heavy-baryon production at LHCb, is used to produce the signal. Decays of hadronic particles are described by EvtGen [25], in which final-state radiation is generated using Photos [26]. The interaction of the generated particles with the detector and their response are implemented using the Geant4 toolkit [27] as described in ref. [29]. The Ξ ++ cc → D + pK − π + decays are generated with a Ξ ++ cc mass of 3621.40 MeV/c 2 and the decay products of Ξ ++ cc and D + hadrons are distributed uniformly in phase space. The procedure to trigger, reconstruct and select candidates is designed to retain Ξ ++ cc signal and to suppress three primary sources of background: combinatorial background, which arises from random combination of tracks; misreconstructed charm or beauty hadron decays, which typically have displaced decay vertices; and combinations of a real D + meson with other tracks to form a fake Ξ ++ cc candidate. To better control systematic uncertainties, the selection of Ξ ++ cc → D + pK − π + decays is also designed to be as similar as possible to that of the Ξ ++ cc → Λ + c K − π + π + normalisation channel, described in ref. [1]. The D + candidates are reconstructed in the final state K − π + π + . At least one of the three tracks used to reconstruct the D + candidate must be selected by the inclusive software trigger, which requires that the track has p T > 250 MeV/c and χ 2 IP > 4 with respect to any PV, where χ 2 IP is defined as the difference in χ 2 of a given PV reconstructed with and without the considered track. The D + candidate then must be reconstructed and accepted by a dedicated D + → K − π + π + selection algorithm in the software trigger. This algorithm applies several geometric and kinematic requirements; at least one of the three tracks must have p T > 1 GeV/c and χ 2 IP > 50, at least two of the tracks must have p T > 0.4 GeV/c and χ 2 IP > 10 and the scalar sum of the p T of the three tracks must be larger than 3 GeV/c. Furthermore, the D + candidate must have a good vertex-fit quality with χ 2 /ndf < 6. The candidate must also point back to its associated PV, where the angle between its flight path and momentum vector should be less than 0.01 radians. The associated PV is that which best fits the flight direction of the reconstructed candidate. The D + vertex must also be displaced from this PV such that the estimated D + decay time is longer than 0.4 ps. Only candidates whose invariant mass is within ±80 MeV/c 2 of the known mass of the D + meson (1869.65 MeV/c 2 [13]) are retained. Finally, candidates are required to pass a MatrixNet classifier [11] within the software trigger, which has been trained on p T and vertex χ 2 information prior to data taking. For events that pass the online trigger, the offline selection of D + candidates proceeds in a similar fashion to that used in the software trigger: three tracks are required to form a common vertex that is significantly displaced from the associated PV of the candidate and its combined invariant mass must be in the range [1847, 1891] MeV/c 2 . Particle identification (PID) requirements are imposed on all three tracks to suppress combinatorial background and misidentified charm decays. The Ξ ++ cc candidates are formed by combining a D + candidate with three more charged tracks, each with p T > 500 MeV/c and separately identified as a proton, kaon and pion with good track quality. The three tracks and the D + candidate are required to form a vertex in which each pairwise combination of the four particles is required to have a distance of closest approach of less than 10 mm and the fitted Ξ ++ cc vertex must have χ 2 /ndf < 10. The Ξ ++ cc candidate is also required to point back to the PV, and to have p T > 4.5 GeV/c. Only events that passed the hardware trigger based on information from the muon and calorimeter systems that are not part of the reconstructed Ξ ++ cc event are used in the analysis [11]. Hence, the event is triggered independently of the reconstructed Ξ ++ cc candidate, which reduces the systematic uncertainty on the efficiency ratios between the Ξ ++ cc → D + pK − π + and Ξ ++ cc → Λ + c K − π + π + decay modes.

JHEP10(2019)124
To improve the mass resolution, the following mass estimator is used in the analysis where M (D + pK − π + ) is the measured invariant mass of the Ξ ++ cc candidate, M ([K − π + π + ] D + ) is the measured invariant mass of the K − π + π + combination corresponding to the intermediate D + candidate and M PDG (D + ) is the known mass of the D + meson. By using the mass definition in eq. (3.1), a mild correlation between decay time and mass is reduced and the mass resolution is improved by 0.15 MeV/c 2 . The Ξ ++ cc candidates are accepted if they have a reconstructed mass in the range 3300 ≤ m(D + pK − π + ) ≤ 3800 MeV/c 2 .
Following a comparison study of different multivariate methods, a classifier based on the multilayer perceptron (MLP) algorithm [30] is used to further suppress combinatorial background. Simulated Ξ ++ cc decays are used to train the MLP classifier to recognise signal. Dedicated software triggers reconstruct an unphysical combination of D + pK + π + (wrong-sign-plus, WSP) and D + pK − π − (wrong-sign-minus, WSM) data. The WSP and WSM samples are expected to be good proxies for combinatorial background in the Ξ ++ cc → D + pK − π + (right-sign, RS) channel. For this analysis, WSP data in the 3550 ≤ m(D + pK − π + ) ≤ 3700 MeV/c 2 mass region is used to train the MLP classifier to identify background, while the WSM data is used to cross-check the results. Fifteen input variables are used in the MLP training. The variables with the best discriminating power between signal and background are: the Ξ ++ cc vertex fit with a kinematic refit [31] of the Ξ ++ cc decay chain requiring it to originate from its PV; the smallest p T of the four decay products of the Ξ ++ cc candidate; the angle between the Ξ ++ cc momentum vector and the direction from the PV to the Ξ ++ cc decay vertex; the χ 2 IP of the Ξ ++ cc candidate with respect to its PV; the maximum distance of the closest approach between all pairs of Ξ ++ cc tracks forming the Ξ ++ cc candidate; and the maximum distance of the closest approach between all pairs of tracks from the decay of the D + candidate. To maintain a sizeable number of signal events, the hardware-trigger requirements are not applied to the signal and background samples. In addition to the training samples, disjoint testing samples are acquired from the same source. After training, the response of the MLP is compared between the training and testing samples. No signs of the MLP classifier being overtrained are found based on the Kolmogorov-Smirnov test statistic. Candidates are retained only if the MLP response output exceeds a certain threshold. The threshold is chosen by maximising the Punzi figure of merit [32], with a target significance of five sigma. To test for potential misreconstruction effects, the same selection criteria are applied to the WSP and WSM data; no peaking structures are visible in either control sample, as expected.
After the multivariate selection, events may contain multiple Ξ ++ cc candidates. This can arise from mistakes in the reconstruction of Ξ ++ cc → D + (→ K − π + π + )pK − π + decays. For instance, there can be cases when Ξ ++ cc candidates in the same event have used the same track more than once. To deal with this, the angle between any two tracks of the same charge is required to be greater than 0.5 mrad. If a Ξ ++ cc candidate has been formed from at least one pair of these cloned tracks, then the candidate is removed. This requirement removes around 6% of Ξ ++ cc candidates in RS data following the multivariate selection. In a separate scenario, the same six final-state tracks may be used to reconstruct more π + , wrong-signplus (red) D + pK + π + and wrong-sign-minus (blue) D + pK − π − data combinations are shown. The control samples have been normalised to the right-sign sample. than one Ξ ++ cc candidate in the same event but with the tracks wrongly interchanged (e.g., the K − track originating from the Ξ ++ cc decay vertex and the K − track coming from the D + decay vertex). In this situation, only one of the Ξ ++ cc candidate from such an event, chosen at random, is retained. This requirement discards less than 1% of candidates at this stage of the selection.

Mass distributions
To determine the yield of Ξ ++ cc and D + particles following the selection of Ξ ++ cc → D + pK − π + candidates, the m(D + pK − π + ) and M ([K − π + π + ] D + ) mass distributions are fitted using models that are developed using simulation.
The invariant-mass distribution M ([K − π + π + ] D + ) of the D + candidates after the candidate selection is shown in figure 2 (left). A Crystal Ball function with exponential tails on both sides [33] is used to model the signal component and a linear function is used to fit the background contribution. The parameters of the signal model are fixed to values obtained from simulation, while all parameters in the background model are free. The selection retains 2697 D + candidates with a purity of 80% according to the results of the fit to the mass spectrum.
The invariant-mass distributions in the RS, WSP and WSM data samples after the candidate selection are shown in figure 2 (right). All the samples have similar smoothly shaped distributions across the entire mass range studied.
The invariant-mass distribution of the Ξ ++ cc candidates, m(D + pK − π + ), for the signal decay mode after applying all requirements of the analysis, is shown in figure 3 (left). The mass distribution is fitted with an unbinned extended maximum-likelihood method, assuming only a background contribution, described by a second-order Chebyshev polynomial.  No signal peak is visible in the spectrum and the local p-value is calculated as a function of mass and shown in figure 3 (right). The local p-value is defined as the probability of observing data that is less compatible with the background-only hypothesis than the data set. The test statistic used is based on q 0 in ref. [34], but instead of assigning it the value zero when observing fewer than expected candidates, it is assigned the value −q 0 to achieve a more intuitive behaviour of the p-value for downward fluctuations. The likelihoods are evaluated with Poisson statistics using the predicted number of background candidates and observed number of signal candidates in regions of ±3σ m around each hypothetical mass, where σ m = 2.8 MeV/c 2 is the Ξ ++ cc mass resolution determined from simulated Ξ ++ cc → D + pK − π + decays. There is no visible signal near the mass of 3620 MeV/c 2 where a Ξ ++ cc signal would be expected, nor is there any excess of candidates near the mass of 3520 MeV/c 2 where the hypothetical isospin partner was observed by the SELEX collaboration [4,5]. The global p-value, including the look-elsewhere effect in the mass range 3500-3800 MeV/c 2 , is 26% and only one signal candidate is observed in the mass range from the kinematic threshold of 3441 MeV/c 2 to 3500 MeV/c 2 . Hence, no significant signal is observed in the mass range from the kinematic threshold to 3800 MeV/c 2 and we proceed to set a limit on the relative branching fraction R.
The invariant-mass distribution of the Ξ ++ cc candidates, m(Λ + c K − π + π + ), for the normalization decay mode, Ξ ++ cc → Λ + c K − π + π + , is shown in figure 4. In this case a signal peak is clearly visible. Both the candidate selection and the modelling of the mass spectrum are identical to that in ref. [1], except for the additional requirements on the hardware trigger. An extended unbinned maximum-likelihood fit to this invariant-mass distribution returns a signal yield of 184 ± 29.

Efficiency determination
To set an upper limit on the ratio R, it is necessary to evaluate the ratio of efficiencies between the Ξ ++ cc → D + pK − π + and Ξ ++ cc → Λ + c K − π + π + decay modes. The efficiency ratio may be factorised as where efficiencies are evaluated for the geometric acceptance (acc), the reconstruction and selection excluding particle identification requirements (sel), the particle identification requirements (PID) and the trigger (trig). Each factor is the efficiency relative to all previous steps in the order given above. The individual ratios are evaluated with simulated Ξ ++ cc decays, assuming a uniform phase space model, except for PID which is derived from data [18,35]. The efficiencies are corrected for known differences between simulation and data, apart from the geometric acceptance.
The individual efficiency components, shown in eq. (5.1), are found to be similar between the two Ξ ++ cc decay modes, except for the reconstruction and selection efficiency, ε sel|acc , where in the Ξ ++ cc → D + pK − π + channel it is found to be approximately twice as large as that of the Ξ ++ cc → Λ + c K − π + π + decay. This leads to a total efficiency ratio of ε(Λ + c K − π + π + )/ε(D + pK − π + ) = 0.46 ± 0.01, where the uncertainty is statistical only. Combining this total relative efficiency with the value for N (Λ + c K − π + π + ) obtained in section 4 and the known values for the branching fractions B(D + → K − π + π + ) and B(Λ + c → pK − π + ), then according to eq. (1.2), the single-event sensitivity is α s = (1.74 ± 0.29) × 10 −3 . The uncertainty on α s includes the total uncertainty -8 -

Systematic uncertainties
The statistical uncertainty on the measured signal yield in the Ξ ++ cc → Λ + c K − π + π + channel is the dominant uncertainty on α s and the systematic uncertainties on α s have small effect on the upper limits on the ratio R.
The largest systematic uncertainty arises from the evaluation of the efficiency of the hardware-trigger requirement. Only candidates that are triggered independently of the Ξ ++ cc candidate's final-state tracks are used in the branching fraction ratio limit to minimise this systematic uncertainty. The ratio of these efficiencies is equal to one if the kinematic distributions of the Ξ ++ cc candidate in the Ξ ++ cc → D + pK − π + and Ξ ++ cc → Λ + c K − π + π + decay modes are identical. However, the efficiencies can be different if the respective selection requirements of the Ξ ++ cc → D + pK − π + and Ξ ++ cc → Λ + c K − π + π + analyses select different kinematic regions of the Ξ ++ cc candidate. This effect is studied by weighting the p T distributions in simulated samples. The change in efficiency of the hardware trigger after the weighting is evaluated and results in a systematic uncertainty of 3.5%. The impact of the model used to fit the m(Λ + c K − π + π + ) invariant-mass distribution on the yield of Ξ ++ cc candidates, N (Λ + c K − π + π + ), is investigated by using alternative signal and background models and performing the fit over different mass ranges. The largest variation in the yield of Ξ ++ cc candidates is 3.1% and this is taken as a systematic uncertainty on α s . The effect of the uncertainty associated with the Ξ ++ cc baryon's lifetime on the relative reconstruction and selection efficiency between the Ξ ++ cc → D + pK − π + and Ξ ++ cc → Λ + c K − π + π + channels is investigated by varying the lifetime within its uncertainty and a systematic uncertainty of 2.9% is assigned to the α s parameter. The PID efficiency is determined in bins of particle momentum and pseudorapidity using calibration samples taken from data [35]. The size of the bins is increased or decreased by a factor of two and the largest deviation on α s of 1.5% is assigned as systematic uncertainty. Finally, since the simulation may not describe the signal perfectly, simulated Ξ ++ cc → D + pK − π + decays are weighted to make their p T distribution match that observed in the Ξ ++ cc → Λ + c K − π + π + data. The selection and software-trigger efficiencies are similarly calculated using p T -corrected simulated Ξ ++ cc decays. The number of p T bins used is increased or decreased by a factor of two and the efficiencies are recalculated for both decay channels. This results in a change in α s of 1.2%. All efficiencies calculated from simulation are averaged over the entire phase space assuming a uniform distribution for both the Ξ ++ cc → D + pK − π + and Ξ ++ cc → Λ + c K − π + π + decays. The phase-space distributions of the selected candidates are uniform and show agreement in data and simulation. Therefore, no systematic uncertainty is assigned to the relative selection and reconstruction efficiencies for the effect of intermediate resonances in their decay. Table 1 summarises the systematic and statistical uncertainties on α s . The statistical uncertainty is dominated by the uncertainty of the yield of the normalisation mode but includes a small contribution from the finite size of the simulated samples. The ratio of  the branching fractions B(D + → K − π + π + ) and B(Λ + c → pK − π + ) have a combined uncertainty of 5.7%. The systematic uncertainties from the different sources discussed above are considered uncorrelated and are added in quadrature to give a total systematic uncertainty of 5.8%. Adding all sources of uncertainty in quadrature gives a total uncertainty of 17.7% on the α s parameter.

Results
In this analysis no significant Ξ ++ cc → D + pK − π + signal is observed so an upper limit is set on the ratio of branching fractions, R. The CLs method [36] is used to determine the ratio of confidence levels (CL) between the signal-plus-background and background-only hypotheses. The upper limit is obtained from the total number of candidates, N obs , observed in the expected signal mass region. This value is calculated by counting the number of candidates within the mass region, 3612 < m(D + pK − π + ) < 3630 MeV/c 2 (indicated by two dashed blue lines in the left-hand plot of figure 3). This mass region corresponds to approximately ±3σ m around the average mass of the Ξ ++ cc state. The CLs score for a possible value of ratio R is calculated as where N s is sampled from the distribution of the expected number of signal candidates for a given ratio R, N b is sampled from the distribution of the expected number of background candidates predicted by the background-only fit (figure 3, left) and P indicates the probability that these statistical quantities are smaller than N obs . The data points in the mass region 3612 < m(D + pK − π + ) < 3630 MeV/c 2 are removed for the fit and N b is determined by performing an integral extrapolation. The probability requirements in the numerator and denominator of eq.

Conclusions
Following observations of the Ξ ++ cc → Λ + c K − π + π + and Ξ ++ cc → Ξ + c π + decay modes, a search for the decay Ξ ++ cc → D + pK − π + is performed using pp collision data recorded by the LHCb experiment in 2016 at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 1.7 fb −1 . No significant signal is found in the mass range from the kinematic threshold of the decay of 3441 MeV/c 2 to 3800 MeV/c 2 . Considering the statistical and systematic uncertainties, an upper limit on the ratio of branching fractions between the Ξ ++ cc → D + pK − π + and Ξ ++ cc → Λ + c K − π + π + decay is set to be R < 1.7 (2.1) × 10 −2 at the 90% (95%) confidence level at the known mass of the Ξ ++ cc baryon. The upper limit on the ratio of branching fractions between the two Ξ ++ cc decay modes is derived assuming a uniform phase space model in the efficiency determinations. A better theoretical understanding of the resonant and nonresonant contributions underpinning the Ξ ++ cc → D + pK − π + and Ξ ++ cc → Λ + c K − π + π + decay processes is required to understand the at least two orders of magnitude difference between the branching fractions of the two Ξ ++ cc decay modes. Dynamical effects or spin constraints in the resonance structures could be suppressing the Ξ ++ cc → D + pK − π + decay. The full dataset from LHCb, or future data taking with the upgraded detector, may reveal evidence of this decay and then shed more light on the production and decay dynamics of the Ξ ++ cc baryon.