Pluto: a Monte Carlo simulation tool for hadronic physics

Pluto is a Monte-Carlo event generator designed for hadronic interactions from Pion production threshold to intermediate energies of a few GeV per nucleon, as well as for studies of heavy ion reactions. The package is entirely based on ROOT, without the need of additional packages, and uses the embedded C++ interpreter of ROOT to control the event production. The generation of events based on a single reaction chain and the storage of the resulting particle objects can be done with a few lines of a ROOT-macro. However, the complete control of the package can be taken over by the steering macro and user-defined models may be added without a recompilation of the framework. Multi-reaction cocktails can be facilitated as well using either mass-dependent or user-defined static branching ratios. The included physics uses resonance production with mass-dependent Breit-Wigner sampling. The calculation of partial and total widths for resonances producing unstable particles is performed recursively in a coupled-channel approach. Here, particular attention is paid to the electromagnetic decays, motivated by the physics program of HADES. The thermal model supports 2-component thermal distributions, longitudinal broadening, radial blast, direct and elliptic flow, and impact-parameter sampled multiplicities. The interface allows angular distribution models (e.g. for the primary meson emission) to be attached by the user as well as descriptions of multi-particle correlations using decay chain templates. The exchange of mass sampling or momentum generation models is also possible. The first feature allows for consistent coupled-channel calculations, needed for a correct description of hadronic interactions. For elementary reactions, angular distribution models for selected channels are already part of the framework, based on parameterizations of existing data. This report gives an overview of the design of the package, the included models and the user interface

Amplitude analysis and branching fraction measurement of Ds+→K+K−π+D_s^+ → K^+ K^− π^+

Utilizing a data set corresponding to an integrated luminosity of 6.32~$\rm fb^{-1}$, recorded by the BESIII detector at center-of-mass energies between 4.178 and 4.226~GeV, we perform an amplitude analysis of the decay $D_{s}^{+} \to K_{S}^{0}\pi^{+}\pi^{0}$ and determine the relative fractions and phase differences of different intermediate processes, which include $K_{S}^{0}\rho(770)^{+}$, $K_{S}^{0}\rho(1450)^{+}$, $K^{*}(892)^{0}\pi^{+}$, $K^{*}(892)^{+}\pi^{0}$, and $K^{*}(1410)^{0}\pi^{+}$. Using a double-tag technique, and making an efficiency correction that relies on our knowledge of the phase-space distribution of the decays coming from the amplitude analysis, the absolute branching fraction is measured to be $\mathcal{B}(D_{s}^{+} \to K_{S}^{0}\pi^{+}\pi^{0})=(5.43\pm0.30_{\text{stat}}\pm 0.15_{\text{syst}})\times 10^{-3}$

Partial wave analysis of J/ψ→γη′η′

Using a sample of (10.09±0.04)×109 J/ψ events collected with the BESIII detector, a partial wave analysis of J/ψ→γη′η′ is performed.The masses and widths of the observed resonances and their branching fractions are reported. The main contribution is from J/ψ→γf0(2020) with f0(2020)→η′η′, which is found with a significance of greater than 25σ. The product branching fraction B(J/ψ → γf0(2020))⋅B(f0(2020) → η′η′ is measured to be (2.63±0.06(stat.) + 0.31−0.46(syst.))×10−4

Search for the decay D0→π0ν¯ν

We present the first experimental search for the rare charm decay D0→π0ν¯ν. It is based on an e+e− collision sample consisting of 10.6×10^6 pairs of D0¯D0 mesons collected by the BESIII detector at √s=3.773 GeV, corresponding to an integrated luminosity of 2.93 fb^−1. A data-driven method is used to ensure the reliability of the background modeling. No significant D0→π0ν¯ν signal is observed in data and an upper limit of the branching fraction is set to be 2.1×10^-4 at the 90% confidence level. This is the first experimental constraint on charmed-hadron decays into dineutrino final states

Amplitude analysis and branching fraction measurement of D+s→K−K+π+π0

Relative fractions and phases of the intermediate decays are determined. With the detection efficiency estimated by the results of the amplitude analysis, the branching fraction of Dþ s → K−Kþπþπ0 decay is measured to be ð5.42 0.10stat 0.17systÞ%

Measurements of the center-of-mass energies of e+e− collisions at BESIII

During the 2016-17 and 2018-19 running periods, the BESIII experiment collected 7.5~fb−1 of e+e− collision data at center-of-mass energies ranging from 4.13 to 4.44~GeV. These data samples are primarily used for the study of excited charmonium and charmoniumlike states. By analyzing the di-muon process e+e−→(γISR/FSR)μ+μ−, we measure the center-of-mass energies of the data samples with a precision of 0.6 MeV. Through a run-by-run study, we find that the center-of-mass energies were stable throughout most of the data-taking period

Measurements of the center-of-mass energies of e+e- collisions at BESIII

During the 2016-17 and 2018-19 running periods, the BESIII experiment collected 7.5 fb -1 of e+e− collision data at center-of-mass energies ranging from 4.13 to 4.44 GeV. These data samples are primarily used for the study of excited charmonium and charmoniumlike states. By analyzing the di-muon process e+e− (γISR/FSR)µ -> +µ-, we measure the center-of-mass energies of the data samples with a precision of 0.6 MeV. Through a run-by-run study, we find that the center-of-mass energies were stable throughout most of the data-collection period

Observation of the decay D_(s)⁺ → ωπ⁺η

Using 7.33 fb−1 of e+e− collision data collected by the BESIII detector at center-of-mass energies between 4.128 and 4.226~GeV, we observe for the first time the decay D±s→ωπ±η with a statistical significance of 7.6σ. The measured branching fraction of this decay is (0.54±0.12±0.04)%, where the first uncertainty is statistical and the second is systematic

Measurement of cross section for e+⁢e−→Ξ−⁢¯Ξ+ near threshold at BESIII

The Born cross sections and effective form factors for process +⁢−→Ξ−⁢¯Ξ+ are measured at eight center-of-mass energies between 2.644 and 3.080 GeV, using a total integrated luminosity of 363.9  pb−1 +⁢− collision data collected with the BESIII detector at BEPCII. After performing a fit to the Born cross section of +⁢−→Ξ−⁢¯Ξ+, no significant threshold effect is observed