Observation of the decay $\Lambda_c^+\rightarrow \Sigma^- \pi^+\pi^+\pi^0$

TL;DR

This paper reports the first observation and measurement of the decay $oldsymbol{ ext{Lambda}_c^+ ightarrow ext{Sigma}^- ext{pi}^+ ext{pi}^+ ext{pi}^0}$, providing new data on charmed baryon decay modes using $e^+e^-$ collision data at BESIII.

Contribution

The paper presents the first observation of a specific Lambda_c^+ decay mode and provides precise measurements of its branching fraction, improving understanding of charmed baryon decays.

Findings

First observation of $ ext{Lambda}_c^+ ightarrow ext{Sigma}^- ext{pi}^+ ext{pi}^+ ext{pi}^0$ decay.

Branching fraction measured as $(2.11 \pm0.33 ext{(stat.)} ext{±}0.14 ext{(syst.)}) ext{ extperthousand}$.

Improved measurement of $ ext{Lambda}_c^+ ightarrow ext{Sigma}^- ext{pi}^+ ext{pi}^+$ decay.

Abstract

We report the first observation of the decay $\Lambda^+_{c}\rightarrow \Sigma^- \pi^+\pi^+\pi^0$, based on data obtained in $e^+e^-$ annihilations with an integrated luminosity of 567~pb$^{-1}$ at $\sqrt{s}=4.6$~GeV. The data were collected with the BESIII detector at the BEPCII storage rings. The absolute branching fraction $\mathcal{B}(\Lambda^+_{c}\rightarrow\Sigma^-\pi^+\pi^+\pi^0)$ is determined to be $(2.11\pm0.33({\rm stat.})\pm0.14({\rm syst.}))\%$. In addition, an improved measurement of $\mathcal{B}(\Lambda^+_{c}\rightarrow\Sigma^-\pi^+\pi^+)$ is determined as $(1.81\pm0.17({\rm stat.})\pm0.09({\rm syst.}))\%$.

Tables

Table 1. Table 1: Requirements on Δ ​ E Δ 𝐸 \Delta E and ST yields N Λ ¯ c − subscript 𝑁 superscript subscript ¯ Λ 𝑐 N_{\bar{\Lambda}_{c}^{-}} for the eleven ST modes. The uncertainties are statistical only.

Mode	$\Delta E$ (GeV)	$N_{{\overline{\Lambda }}_c^{-}}$
$\overline{p}{K}_S^0$	$[-0.025,0.028]$	$1066\pm 33$
$\overline{p}{K}^{+}{\pi }^{-}$	$[-0.019,0.023]$	$5692\pm 88$
$\overline{p}{K}_S^0{\pi }^0$	$[-0.035,0.049]$	$593\pm 41$
$\overline{p}{K}^{+}{\pi }^{-}{\pi }^0$	$[-0.044,0.052]$	$1547\pm 61$
$\overline{p}{K}_S^0{\pi }^{+}{\pi }^{-}$	$[-0.029,0.032]$	$516\pm 34$
$\overline{\Lambda }{\pi }^{-}$	$[-0.033,0.035]$	$593\pm 25$
$\overline{\Lambda }{\pi }^{-}{\pi }^0$	$[-0.037,0.052]$	$1864\pm 56$
$\overline{\Lambda }{\pi }^{-}{\pi }^{+}{\pi }^{-}$	$[-0.028,0.030]$	$674\pm 36$
${\overline{\Sigma }}^0{\pi }^{-}$	$[-0.029,0.032]$	$532\pm 30$
${\overline{\Sigma }}^{-}{\pi }^0$	$[-0.038,0.062]$	$329\pm 28$
${\overline{\Sigma }}^{-}{\pi }^{+}{\pi }^{-}$	$[-0.049,0.054]$	$1009\pm 57$

Table 2. Table 2: Summary of the relative systematic uncertainties Δ Σ − ​ π + ​ π + syst subscript superscript Δ syst superscript Σ superscript 𝜋 superscript 𝜋 \Delta^{\rm syst}_{\Sigma^{-}\pi^{+}\pi^{+}} and Δ Σ − ​ π + ​ π + ​ π 0 syst subscript superscript Δ syst superscript Σ superscript 𝜋 superscript 𝜋 superscript 𝜋 0 \Delta^{\rm syst}_{\Sigma^{-}\pi^{+}\pi^{+}\pi^{0}} in ℬ ​ ( Λ c + → Σ − ​ π + ​ π + ) ℬ → subscript superscript Λ 𝑐 superscript Σ superscript 𝜋 superscript 𝜋 \mathcal{B}(\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}) and ℬ ​ ( Λ c + → Σ − ​ π + ​ π + ​ π 0 ) ℬ → subscript superscript Λ 𝑐 superscript Σ superscript 𝜋 superscript 𝜋 superscript 𝜋 0 \mathcal{B}(\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}\pi^{0}) , respectively.

Source	${\Delta }_{{\Sigma }^{-}{\pi }^{+}{\pi }^{+}}^{\mathrm{syst}}$ [%]	${\Delta }_{{\Sigma }^{-}{\pi }^{+}{\pi }^{+}{\pi }^0}^{\mathrm{syst}}$ [%]
${\pi }^{\pm }$ tracking	3.0	3.0
${\pi }^{\pm }$ identification	3.0	3.0
${\pi }^0$ reconstruction	$\cdots$	2.0
Fit to $M_n-M_{n{\pi }^{-}}$	2.0	3.6
Signal modelling	2.0	2.0
MC statistics	0.6	0.7
$N_{{\overline{\Lambda }}_c^{-}}^{\mathrm{tot}}$	1.0	1.0
Total	5.2	6.4

Equations

$$M_n=\sqrt{(E_{\mathrm{beam}}-E_{{\pi }^{+}{\pi }^{+}{\pi }^{-}({\pi }^0)}{)}^2-|{\overrightarrow{p}}_{{\Lambda }_c^{+}}-{\overrightarrow{p}}_{{\pi }^{+}{\pi }^{+}{\pi }^{-}({\pi }^0)}{|}^2}$$

$$M_{n{\pi }^{-}}=\sqrt{(E_{\mathrm{beam}}-E_{{\pi }^{+}{\pi }^{+}({\pi }^0)}{)}^2-|{\overrightarrow{p}}_{{\Lambda }_c^{+}}-{\overrightarrow{p}}_{{\pi }^{+}{\pi }^{+}({\pi }^0)}{|}^2}$$

$$\mathcal{B}({\Lambda }_c^{+}\to {\Sigma }^{-}{\pi }^{+}{\pi }^{+}({\pi }^0))$$

$$\displaystyle=\frac{N^{\rm obs}_{\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})}}{N^{\rm tot}_{\bar{\Lambda}_{c}^{-}}\cdot\varepsilon_{\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})}\cdot\mathcal{B}(\Sigma^{-}\rightarrow n\pi^{-})},\leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\ \leavevmode\nobreak\$$

Full text

Keywords:  branching fraction, charmed baryon, weak decays, $e^{+}e^{-}$ annihilation, BESIII

††journal: Physics Letters B

1 Introduction

The study of hadronic decays of charmed baryons provides important information to understand both the strong and the weak interactions [1]. It also provides essential input to understand background contributions in the study of $b$-baryon physics, as $\Lambda_{b}$ decays dominantly to $\Lambda_{c}^{+}$. More than 30 years have passed since the $\Lambda^{+}_{c}$ baryon was first observed in $e^{+}e^{-}$ annihilations by the Mark II experiment [2] and the knowledge of $\Lambda^{+}_{c}$ decays remains very poor compared to that for charmed mesons. So far, measured decay modes account for only about 60% [3] of all $\Lambda_{c}^{+}$ decays, primarily consisting of modes with a $\Lambda(\Sigma)$ hyperon or a proton in the final state. Decays to the $\Sigma^{-}$ hyperon are Cabibbo-allowed and are expected to have large rates. However, no experimental measurements exist except for $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}$ [3]. Therefore, searching for additional decay modes with $\Sigma^{-}$ in the final state is important to build up knowledge on $\Lambda_{c}^{+}$ decays. In this paper, we report the first observation of the so-far undetermined, but expected to be large, decay of $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}\pi^{0}$ [4]. In addition, we perform the first absolute measurement of the branching fraction for $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}$.

The data analyzed in this work corresponds to an integrated luminosity of $567$ pb*-1* [5] of $e^{+}e^{-}$ annihilations at center-of-mass energy (c.m.) $\sqrt{s}=4.6$ GeV by the BEPCII collider and collected with the BESIII detector [6]. The c.m. energy is slightly above the threshold for the production of $\Lambda_{c}^{+}\bar{\Lambda}_{c}^{-}$, so $\Lambda_{c}^{+}\bar{\Lambda}_{c}^{-}$ pairs are produced with no additional hadrons. The analysis technique in this work, which was first applied in the Mark III experiment [7], is optimized for measuring charm hadron pairs produced near threshold. First, we select the subset of our events in which a $\bar{\Lambda}_{c}^{-}$ is reconstructed in an exclusive hadronic decay mode, designated as the single-tag (ST) sample. Events in this ST sample are then searched for the signal channel $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})$ in the system recoiling against the ST to select double tag (DT) events. In the final states of $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})$, the $\Sigma^{-}$ hyperon is detected through $\Sigma^{-}\rightarrow n\pi^{-}$. As the neutron is not reconstructed in this analysis, we deduce its kinematic properties by four-momentum conservation. The absolute branching fraction (BF) of $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})$ is derived from the probability of detecting the DT signals in the ST sample. Hence, this method provides a clean and straightforward BF measurement that is independent of the number of $\Lambda^{+}_{c}\bar{\Lambda}^{-}_{c}$ events produced.

2 BESIII Detector and Monte Carlo Simulation

BESIII [6] is a cylindrical detector with a coverage of 93% of the full $4\pi$ solid angle. It consists of a Helium-gas based main drift chamber (MDC), a plastic scintillator time-of-flight (TOF) system, a CsI (Tl) electromagnetic calorimeter (EMC), a superconducting solenoid providing a 1.0 T magnetic field, and a muon detection system in the iron flux return of the magnet. The charged particle momentum resolution is 0.5% at a transverse momentum of 1 GeV/$c$. The photon energy resolution at 1 GeV is 2.5% in the central barrel region and 5.0% in the two end caps. More details about the design and performance of the detector are given in Ref. [6].

A GEANT4-based [8] Monte Carlo (MC) simulation package, which includes the geometric description of the detector and the detector response, is used to determine the detection efficiency and to estimate the potential backgrounds. MC samples of the signal mode $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})$, together with a $\bar{\Lambda}^{-}_{c}$ decaying to specified ST modes, are generated with KKMC [9] and EVTGEN [10], taking into account initial-state radiation (ISR) [11] and final-state radiation [12] effects. The $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})$ decay is simulated by reweighting the phase-space-generated MC events to approximate observed kinematic distributions in data. To understand potential background contributions, an inclusive MC sample is used. It includes generic $\Lambda_{c}^{+}\bar{\Lambda}_{c}^{-}$ events, $D_{(s)}^{(*)}\bar{D}_{(s)}^{(*)}+X$ production, ISR return to the charmonium states at lower masses and continuum $q\bar{q}$ processes. Previously measured decay modes of the $\Lambda_{c}$, $\psi$ and $D_{(s)}$ are simulated with EVTGEN, using BFs from the Particle Data Group (PDG) [3]. The unknown decays of the $\psi$ states are generated with LUNDCHARM [13].

3 Analysis

The ST and DT selection technique that is used in our analysis follows closely the one used and described in Ref. [14]. We reconstruct the $\bar{\Lambda}^{-}_{c}$ baryons in the eleven hadronic decay modes listed in Table 1. Intermediate particles are reconstructed through their decays $K^{0}_{S}\rightarrow\pi^{+}\pi^{-}$, $\bar{\Lambda}\rightarrow\bar{p}\pi^{+}$, $\bar{\Sigma}^{0}\rightarrow\gamma\bar{\Lambda}$ with $\bar{\Lambda}\rightarrow\bar{p}\pi^{+}$, $\bar{\Sigma}^{-}\rightarrow\bar{p}\pi^{0}$, and $\pi^{0}\rightarrow\gamma\gamma$. The selection criteria for the proton, kaon, pion, $\pi^{0}$, $K^{0}_{S}$ and $\bar{\Lambda}$ candidates used in the reconstruction of the ST signals are described in Ref. [14].

The ST $\bar{\Lambda}^{-}_{c}$ signals are identified using the beam-energy-constrained mass, $M_{\rm BC}=\sqrt{E^{2}_{\rm beam}-|\vec{p}_{\bar{\Lambda}^{-}_{c}}|^{2}}$, where $E_{\rm beam}$ is the beam energy and $\vec{p}_{\bar{\Lambda}^{-}_{c}}$ is the momentum of the $\bar{\Lambda}^{-}_{c}$ candidate in the rest frame of the initial $e^{+}e^{-}$ system [15]. To improve the signal purity, the energy difference $\Delta E=E_{\rm beam}-E_{\bar{\Lambda}^{-}_{c}}$ for each candidate is required to be within approximately $\pm 3\sigma$ of the $\Delta E$ signal peak position, where $\sigma$ is the $\Delta E$ resolution and $E_{\bar{\Lambda}^{-}_{c}}$ is the reconstructed $\bar{\Lambda}^{-}_{c}$ energy. Table 1 shows the mode-dependent $\Delta E$ requirements and the ST yields in the $M_{\rm BC}$ signal region $(2.280,2.296)$ GeV/$c^{2}$, which are obtained by fits to the $M_{\rm BC}$ distributions. See Ref. [14] for more details. The total ST yield is $N^{\rm tot}_{\bar{\Lambda}^{-}_{c}}=14415\pm 159$, where the uncertainty is statistical only.

Candidates for the decay $\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})$ with $\Sigma^{-}\rightarrow n\pi^{-}$ are reconstructed from the tracks not used in the ST $\bar{\Lambda}^{-}_{c}$ reconstruction. It is required that there are only three charged tracks in the system recoiling against the $\bar{\Lambda}_{c}^{-}$ satisfying $|\cos\theta|<0.93$, where $\theta$ is the polar angle with respect to the beam direction. For the two $\pi^{+}$ candidates from the $\Lambda_{c}^{+}$, the distances of closest approach to the interaction point must be within $\pm 10$ cm along the beam direction and within 1 cm in the perpendicular plane, while the $\pi^{-}$ candidate from $\Sigma^{-}$ decay is not subjected to this requirement. Identification of charged tracks is performed by combining the $dE/dx$ information from the MDC and the time of flight measured in the TOF to obtain the probability $\mathcal{L}_{h}$ for each hadron type $h$. The three charged pions must satisfy $\mathcal{L}_{\pi}>\mathcal{L}_{K}$. Photon candidates are reconstructed from isolated clusters in the EMC in the regions $|\cos\theta|\leq 0.80$ (barrel) and $0.86\leq|\cos\theta|\leq 0.92$ (end cap). The deposited energy of a neutral cluster is required to be larger than 25 (50) MeV in the barrel (end cap) region, and the angle between the photon candidate and the nearest charged track must be larger than 10∘. To suppress electronic noise and energy deposits unrelated to the event, the difference between the EMC time and the event start time is required to be within $(0,700)$ ns. To reconstruct $\pi^{0}$ candidates, the invariant mass of photon pairs is required to be within (0.110, 0.155) GeV/$c^{2}$ and, as a second step, a kinematic fit is implemented to constrain the $\gamma\gamma$ invariant mass to the nominal $\pi^{0}$ mass [3].

The kinematic variable

[TABLE]

is computed to characterize the reconstructed mass of the undetected neutron, where $E_{\pi^{+}\pi^{+}\pi^{-}(\pi^{0})}$ is the energy of the $\pi^{+}\pi^{+}\pi^{-}(\pi^{0})$ combination and $\overrightarrow{p}_{\pi^{+}\pi^{+}\pi^{-}(\pi^{0})}$ is the three-momentum of the $\pi^{+}\pi^{+}\pi^{-}(\pi^{0})$ combination. The expected momentum $\vec{p}_{\Lambda_{c}^{+}}$ of the $\Lambda_{c}^{+}$ is calculated by $\vec{p}_{\Lambda_{c}^{+}}=-\hat{p}_{\rm tag}\sqrt{E_{\rm beam}^{2}-m^{2}_{\Lambda^{+}_{c}}},$ where $\hat{p}_{\rm tag}$ is the direction of the momentum of the ST $\bar{\Lambda}_{c}^{-}$ candidate and $m_{\Lambda_{c}^{+}}$ is the mass of the $\Lambda^{+}_{c}$ taken from the PDG [3]. Similarly, we can construct the variable

[TABLE]

to represent the reconstructed mass of the $\Sigma^{-}$.

The distributions of $M_{n}$ versus $M_{n\pi^{-}}$ for the $\Lambda^{+}_{c}\to\Sigma^{-}\pi^{+}\pi^{+}$ and $\Lambda^{+}_{c}\to\Sigma^{-}\pi^{+}\pi^{+}\pi^{0}$ candidates in data are shown in Figs. 1 (a) and (b), respectively, where clusters corresponding to signal decays are evident. To improve the resolution of the signal mass, as well as to better handle the backgrounds around the $\Sigma^{-}$ and neutron mass regions, we determine the signal yields from the distribution of the mass difference $M_{n\pi^{-}}-M_{n}$, since $M_{n\pi^{-}}$ and $M_{n}$ are highly correlated. Based on a study of the inclusive MC samples, no peaking backgrounds are expected for these two channels. We perform an unbinned maximum likelihood fit to the $M_{n\pi^{-}}-M_{n}$ spectra, as shown in Figs. 1 (c) and (d). In the fits, the signals are described by non-parametric functions extracted from the signal MC convoluted with a Gaussian function accounting for the resolution difference between data and MC, while the background shapes are described with a second-order polynomial function. The width of the Gaussian is left free in the fit, while its mean is fixed to zero. From the fits, we find the DT signal yields $N^{\rm obs}_{\Sigma^{-}\pi^{+}\pi^{+}}=161\pm 15$ and $N^{\rm obs}_{\Sigma^{-}\pi^{+}\pi^{+}\pi^{0}}=88\pm 14$, where the uncertainties are statistical only. Backgrounds from non-$\Lambda_{c}^{+}$ decays are estimated by examining the ST candidates in the $M_{\rm BC}$ sideband $(2.252,2.272)$ GeV/$c^{2}$ in data. The backgrounds from non-$\Lambda_{c}^{+}$ decays are found to be negligible.

The absolute BFs for $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}$ and $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}\pi^{0}$ are determined by

[TABLE]

where $\varepsilon_{\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})}$ is the detection efficiency for the $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})$ decay with $\Sigma^{-}\to n\pi^{-}$. The intermediate decay branching fraction of $\Sigma^{-}\to n\pi^{-}$ is included in the denominator of Eq. (1). For each ST mode $i$, the efficiency $\varepsilon^{i}_{\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})}$ is obtained by dividing the DT efficiency $\varepsilon^{i}_{{\rm tag},\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})}$ by the ST efficiency $\varepsilon^{i}_{\rm tag}$. After weighting $\varepsilon^{i}_{\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})}$ by the mode-by-mode ST yields in data, we find the overall average efficiencies $\varepsilon_{\Sigma^{-}\pi^{+}\pi^{+}}=(61.8\pm 0.4)\%$ and $\varepsilon_{\Sigma^{-}\pi^{+}\pi^{+}\pi^{0}}=(29.0\pm 0.2)\%$, where the branching fraction for $\pi^{0}\rightarrow\gamma\gamma$ is included. Substituting the values of $N^{\rm obs}_{\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})}$, $N^{\rm tot}_{\bar{\Lambda}^{-}_{c}}$, $\varepsilon_{\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})}$ and $\mathcal{B}(\Sigma^{-}\rightarrow n\pi^{-})$ in Eq. (1), we obtain $\mathcal{B}(\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+})=(1.81\pm 0.17\pm 0.09)\%$ and $\mathcal{B}(\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}\pi^{0})=(2.11\pm 0.33\pm 0.14)\%$, where the first uncertainties are statistical, and the second are systematic, as described below.

With the DT technique, the BF measurement is insensitive to uncertainty in the ST efficiencies. The systematic uncertainties in measuring $\mathcal{B}(\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+})$ and $\mathcal{B}(\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}\pi^{0})$ mainly arise from the efficiencies of $\pi$ detection and identification, fits to the $M_{n\pi^{-}}-M_{n}$ distributions and the signal modelling in the MC simulation. The systematic uncertainties in the $\pi^{\pm}$ tracking and identification are both determined to be 1.0% by studying a set of samples of $e^{+}e^{-}\rightarrow\pi^{+}\pi^{-}\pi^{+}\pi^{-}$, $e^{+}e^{-}\rightarrow K^{+}K^{-}\pi^{+}\pi^{-}$ and $e^{+}e^{-}\rightarrow p\bar{p}\pi^{+}\pi^{-}$ obtained from data with c.m. energy above 4.0 GeV. The $\pi^{0}$ reconstruction efficiency is validated by analyzing DT events with $\bar{D}^{0}\rightarrow K^{+}\pi^{-}$ or $K^{+}\pi^{-}\pi^{0}$ versus $D^{0}\rightarrow K^{-}\pi^{+}\pi^{0}$ [16]. The difference of the $\pi^{0}$ reconstruction efficiencies between data and MC simulations is estimated to be 2.0%. The uncertainty from the fit to the $M_{n\pi^{-}}-M_{n}$ distribution is evaluated by checking the relative changes of $N^{\rm obs}_{\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})}$ with different choices for signal shapes (double Gaussian function), background shapes (first-order polynomial function, third-order polynomial function and a MC-derived background shape) and fit ranges ((0.19, 0.34) GeV/$c^{2}$). The uncertainty in modelling the signal process is obtained by varying the reweighting factors of the observed kinematic variables within their statistical uncertainties and extracting the difference of the resultant efficiencies. The difference is estimated to be 2.0% for the studied channels and is taken as the systematic uncertainty due to the signal modelling. In addition, there are systematic uncertainties in obtaining $N^{\rm tot}_{\bar{\Lambda}_{c}^{-}}$ evaluated by using alternative signal shapes in the fits to the $M_{\rm BC}$ spectra [14], resulting in an uncertainty of 1.0%, and in the statistical limitation of the MC samples, which is estimated to be 0.6 (0.7)% for $\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}(\pi^{0})$. The uncertainties from the BFs of $\Sigma^{-}\rightarrow n\pi^{-}$ and $\pi^{0}\rightarrow\gamma\gamma$ are negligible. All of the above systematic uncertainties are summarized in Table 2, and the total uncertainties are evaluated to be 5.2% and 6.4% for $\mathcal{B}(\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+})$ and $\mathcal{B}(\Lambda_{c}^{+}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}\pi^{0})$, respectively, by combining all items in quadrature.

4 Summary

Based on an $e^{+}e^{-}$ collision data sample with an integrated luminosity of 567 pb*-1* taken at $\sqrt{s}=4.6$ GeV with the BESIII detector, we report the first observation of the decay $\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}\pi^{0}$ and the first absolute BF measurement for $\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}$. The results are $\mathcal{B}(\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+})=(1.81\pm 0.17\pm 0.09)\%$ and $\mathcal{B}(\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+}\pi^{0})=(2.11\pm 0.33\pm 0.14)\%$, where the first uncertainties are statistical and the second are systematic.

Our result for $\mathcal{B}(\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+})$ is consistent with and more precise than the previous result [3]. BESIII measured the BF of the isospin symmetric channel $\mathcal{B}(\Lambda^{+}_{c}\rightarrow\Sigma^{+}\pi^{+}\pi^{-})=(4.25\pm 0.24\pm 0.20)\%$ [17]. This allows us to determine the ratio $\mathcal{B}(\Lambda^{+}_{c}\rightarrow\Sigma^{-}\pi^{+}\pi^{+})/\mathcal{B}(\Lambda^{+}_{c}\rightarrow\Sigma^{+}\pi^{+}\pi^{-})=0.42\pm 0.05\pm 0.02$, where the first uncertainty is statistical and the second systematic. The statistical uncertainty of the ratio dominates, as many common systematic uncertainties cancel. This is consistent with and more precise than the value previously measured by the E687 Collaboration $(0.53\pm 0.15\pm 0.07)$ [18].

5 Acknowledgments

The BESIII collaboration thanks the staff of BEPCII and the IHEP computing center for their strong support. This work is supported in part by National Key Basic Research Program of China under Contract No. 2015CB856700; National Natural Science Foundation of China (NSFC) under Contracts Nos. 11125525, 11235011, 11275266, 11305180, 11322544, 11322544, 11335008, 11425524, 11505010; the Chinese Academy of Sciences (CAS) Large-Scale Scientific Facility Program; the CAS Center for Excellence in Particle Physics (CCEPP); Joint Large-Scale Scientific Facility Funds of the NSFC and CAS under Contracts Nos. U1332201, U1532257, U1532258; CAS under Contracts Nos. KJCX2-YW-N29, KJCX2-YW-N45, QYZDJ-SSW-SLH003; 100 Talents Program of CAS; National 1000 Talents Program of China; INPAC and Shanghai Key Laboratory for Particle Physics and Cosmology; German Research Foundation DFG under Contracts Nos. Collaborative Research Center CRC 1044, FOR 2359; Istituto Nazionale di Fisica Nucleare, Italy; Joint Large-Scale Scientific Facility Funds of the NSFC and CAS; Koninklijke Nederlandse Akademie van Wetenschappen (KNAW) under Contract No. 530-4CDP03; Ministry of Development of Turkey under Contract No. DPT2006K-120470; National Science and Technology fund; The Swedish Resarch Council; U.S. Department of Energy under Contracts Nos. DE-FG02-05ER41374, DE-SC-0010118, DE-SC-0010504, DE-SC-0012069; University of Groningen (RuG) and the Helmholtzzentrum fuer Schwerionenforschung GmbH (GSI), Darmstadt; WCU Program of National Research Foundation of Korea under Contract No. R32-2008-000-10155-0. This paper is also supported by Beijing municipal government under Contract Nos. KM201610017009, 2015000020124G064.