M2: Accomplishments of the second funding period

  • η-η′ mixing in large-Nc chiral perturbation theory
    We performed a calculation of the η-η′ mixing in the framework of large-Nc chiral perturbation theory. A general expression for the η-η′ mixing at next-to-next-to-leading order (NNLO) was derived, including higher-derivative terms up to fourth order in the four momentum, kinetic and mass terms. In addition, the axial-vector decay constants of the η-η′ system were determined at NNLO. The numerical analysis of the results was performed successively at LO, NLO, and NNLO.
  • Quantum corrections to the chiral anomaly in large-Nc chiral perturbation theory
    We calculated the anomalous decays η(′) → γ(∗)γ(∗) and η(′) → π+πγ(∗) at the one-loop level up to next-to-next-to-leading order (NNLO) in large-Nc chiral perturbation theory. Both the decays to real photons and the decays involving virtual photons, providing access to the substructure of the mesons, have been discussed. The results were numerically evaluated successively at LO, NLO, and NNLO. The appearing low-energy constants were determined through fits to the available experimental data. In the case of η(′) → γ(∗)γ(∗) we investigated the decay widths to real photons, the widths of η(′) → γl+l, where l = e,μ, and the single-virtual transition form factors. The considered observables of the decays η(′) → π+π−γ(∗) are the spectra of the decays involving a real photon, η(′) → π+πγ, as well as the spectra of η(′) → π+πl+l, where l = e,μ, with respect to the invariant masses of the π+π and l+l systems.
  • Dalitz plot analysis of the decay η′ → ηπ0π0
    We performed a Dalitz plot analysis of the decay η′ → ηπ0π0 using about 124,000 signal decays collected with the Crystal Ball/TAPS setup and the Endpoint Tagger (EPT) in 2014. The A2 dataset analysed in this paper is the largest sample of η′ → ηπ0π0 decays so far and results in the most precise to date determination of the parameters describing the matrix element for this decay.
    The large dataset also allows, for the first time, the observation of a structure below the π+π threshold, which is in good agreement with the cusp that was predicted based on the ππ scattering length combination, a0 − a2, extracted from K → 3π decays.
  • Measurement of the Dalitz plot distribution of η′ → π0π0π0
    [PhD project of Martin Wolfes.] Using the EPT data set of A2, we attempt to reconstruct the decay η′ → π0π0π0, with the goal to perform a Dalitz plot analysis. Initial expectations based on the number of produced η′, the branching fraction, and a signal Monte Carlo study indicated that a signal yield comparable to the BESIII sample of 1900 reconstructed η′ → π0π0π0 would be within reach.
    However, performing the analysis we had to realise that a clean extraction of the signal decay is not possible due to the presence of a large, irreducible, background from 3π0 production which does not proceed via a primary η′. This analysis therefore turns out to be not feasible within the scope of a PhD thesis.
    Instead, a measurement of the cross section γp → Σ+K0 with Σ+ → pπ0 and K0S → π0π0 is being performed. The analysis sees about 52000 reconstructed signal decays; a measurement of the differential cross section in several bins of Eγ is therefore possible. Results are expected for 2018.
  • Measurement of the branching fraction η′ → ωγ
    [This measurement is the PhD project of Andreas Neiser.] The plan is to measure the branching fraction of η′ → ωγ → π0γ relative to the abundant reference channel η′ → γγ. Signal yields of about 50,000 reference and 1300 signal decays are found in the full endpoint tagger data set. However, there is a systematic difference between data and simulation in the signal efficiency, so that the preliminary result for the branching fraction comes out significantly smaller than the PDG average. A publication using an alternative analysis chain is under preparation.
  • Measurement of the branching fraction ω → ηγ
    [PhD project of Oliver Steffen.] The decay ω → ηγ was planned to be measured using the decay ω → π0γ as a reference, using the approximately 60 × 106 ω mesons produced in the EPT beam times. The reference channel is cleanly reconstructed with the expected rates, but identification of the signal decay is seriously hampered by the presence of a large background from the processes γp → ηπ0p, in which one photon from the π0 is not detected. Additional backgrounds are also present, and the distribution of m(ηγ) is not well described by simulation, making a determination of the signal yield very unreliable.
    As an alternative result, the differential cross section for ω photoproduction is measured in bins of the incident photon energy; this agrees within uncertainties with previous results.
  • Study of charmonium decays to three pions
    We developed and applied a data selection for the channels J/ψ → π+π−π0 and ψ′ → π+π−π0 on the full BESIII data sets. Fits to these data using Veneziano amplitudes in close collaboration with A. Szczepaniak (Indiana University) were started and show a very good description of the J/ψ data; more work is needed for a thorough understanding of the ψ′ data.
  • Dalitz plot analysis of D0 → KS0K+K
    We finalised the Dalitz plot analysis and measurement of the branching fraction of the decay D0 →K0SK+Kusing the full 2.92fb−1 data sample taken at √s=3.773 GeV.The branching fraction result is the most precise measurement to date. The uncertainty is dominated by systematic uncertainties and improves on the PDG average by a factor of 1.7. A journal publication is under internal review within BESIII; submission is expected in the first quarter of 2018.
  • Exotic Charmonium states
    Within the last two years we performed the search for the two isospin singlet states ηcη in the reaction channel e+e → ηcηππ and Jψη in the reaction channel e+e → Jψ ηη with data of the BESIII experiment. The first step was the optimization of the extraction of events with the correct signature in the detector.
    To extract the ηcη events we had to sum over 16 different ηc decay channels in order to get a reasonable amount of ηc candidates. Concerning the background contributions we were not able to describe the remaining background distribution with the recent models. Therefore we decided to describe the remaining contribution in the data directly by fitting a combined model of a signal and a background function simultaneously to each of the 16 ηc decay channels. After the study of the systematic uncertainties the upper limit for the production of this reaction channel was calculated for the three energies at 4.23, 4.26 and 4.60 GeV. For this analysis the review process at BESIII already started. A review committee already gave the second round of questions.
    For the Jψη case extensive studies were performed in order to subtract the background contributions. Not for all contributions measured cross sections exist, therefore a rough analysis was done for the missing channels to determine at least the upper limit of those channels. Only very few events survive for the 5 energy points 4.23, 4.26, 4.36, 4.42 and 4.60 GeV. Therefore we calculated the upper limits for each energy and took into account all systematic uncertainties. Next step here is writing the memo to start the review of this analysis.
    In parallel we started the analyses e+e → (DD ̄*)+π and e+e → (DD ̄*)0π0 at 4.42 GeV. The selection of these channels is almost finished. The two independent analyses were cross-checked in detail and differences understood.