Abstract
A high statistics data sample of the decays of \(K^+\) mesons to three charged particles was accumulated by the OKA experiment in 2012 and 2013. This allowed to select a clean sample of about 450 events with \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) decays with the energy of the photon in the kaon rest frame greater than 30 MeV. The measured branching fraction of the \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \), with \(E_{\gamma }^{*} > 30\ \hbox {MeV}\) is equal to \((0.71 \pm 0.05) \times 10^{-5}\). The measured differential branching fraction over photon energy is compared with the prediction of the chiral perturbation theory to \({\mathcal {O}}(p^{4})\). A search for an up-down asymmetry of the photon with respect to the hadronic system decay plane is also performed.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
1 Introduction
The present experimental status of \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) is rather meagre. It was observed in one experiment on statistics of 7 events [1]. The photon energies in these events were low, that did not allow to search for deviations from a simple QED process of photon emission. The measured value of the \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) branching ratio [2] is \((1.04\pm 0.31) \times 10^{-4}\) for \(E_{\gamma }^{*} > 5\ \hbox {MeV}\). In the present analysis we have a possibility for more detailed study of this decay using larger data sample collected by the OKA experiment. This decay has certain interest for the theory, in particular, for the chiral perturbation theory [3,4,5,6] (CHPT), which gives appropriate framework for analysis of such decays. To the lowest order in an expansion in momenta and meson masses, this radiative decay is completely determined [7] by the non-radiative amplitude for \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\). However, at next-to-leading order, a full-fledged CHPT calculation of nonleptonic weak amplitudes of \(\hbox {O}(p^{4})\) is required. Such a calculation is done in [8]. It is interesting to confront these results to the experimental data. In a decay in the beauty sector \(B^{+}\rightarrow K^{+}\pi ^{+}\pi ^{-}\gamma \) the LHCb experiment has found significant up-down asymmetry of the photon with respect to the hadronic system decay plane [9]. This observable was proposed in [10], it is both P and T-odd. We perform an analogous study for the radiative kaon decay \(K\rightarrow 3\pi \gamma \) to search for new physics effects.
2 Separated kaon beam and OKA experiment
The OKA experiment makes use of a secondary hadron beam at the U-70 Proton Synchrotron of NRC “Kurchatov Institute” – IHEP, Protvino, with enhanced fraction of kaons obtained by RF-separation with Panofsky scheme [11]. The OKA setup, Fig. 1, is a double magnetic spectrometer complemented by electromagnetic and hadron calorimeters and a Decay Volume. The first magnetic spectrometer, consisting of the magnet \(\hbox {M}_{1}\) and surrounding 1 mm pitch PC’s (\(\hbox {BPC}_{(1Y)}\), \(\hbox {BPC}_{(2Y,2X)}\), \(\hbox {BPC}_{(3X,3Y)}\), \(\hbox {BPC}_{(4X,4Y)}\)) serves for the beam momentum measurement. It is supplemented by two gas threshold Cherenkov counters Č\(_{1}\), Č\(_{2}\) for kaon identification and by beam scintillation counters \(\hbox {S}_{1}\), \(\hbox {S}_{2}\), \(\hbox {S}_{3}\), \(\hbox {S}_{4}\). The 11 m long Decay Volume (DV) filled with helium contains 11 rings of guard system (GS). To reinforce GS, a gamma detector (BGD), made of lead glass blocks located behind the DV is used as a veto at large angles, while low angle particles pass through a central opening. The wide aperture spectrometric magnet, \(\hbox {SM}_{(SP40A)}\), with a field integral of \(\sim \ 1\ \hbox {Tm}\) serves as a spectrometer for the charged decay products together with corresponding tracking chambers: 2 mm PC’s (\(\hbox {PC}_{1,\ldots ,8}\)), 9 mm diameter straw tubes \(\hbox {ST}_{(1,2,3)}\) and 30 mm diameter drift tubes \(\hbox {DT}_{1,2}\). The matrix hodoscope \(\hbox {HODO}_{(matrix)}\) is composed of 252 scintillator tiles with WLS+SiPM readout. It is used in the trigger, improves time resolution and links x–y projections of a track. Two scintillation counters \(\hbox {S}_{bk1}\), \(\hbox {S}_{bk2}\) serve to suppress undecayed beam particles. At the end of the OKA setup there are two calorimeters: electromagnetic (GAMS-2000) made of lead glass blocks and a hadron one (\(\hbox {HCAL}_{(GDA)}\)) – 100 iron-scintillator sandwiches. Finally, four partially overlapping muon counters are located downstream the HCAL. The data acquisition system of the OKA setup [12] operates at \(\sim \ 25\ \hbox {kHz}\) event rate with the mean event size of \(\sim \ 4\ \hbox {kByte}\). The detailes of the setup and the beam can be found elsewhere [13].
3 The data and the analysis procedure
The study of the decay \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) is performed with the data set accumulated in 2012 and 2013 runs with 17.7 GeV/c beam momentum. Two main triggers were used,both include the same component which selects beam kaons decaying inside the OKA setup. It requires a coincidence of 4 beam scintillation counters, a combination of two Cerenkov’s (Č\(_{1}\) sees pions, Č\(_{2}\) pions and kaons) and, finally, anticoincidence of two scintillation counters \(\hbox {S}_{bk1}\), \(\hbox {S}_{bk2}\), located on the beam axis after the magnet to suppress undecayed beam particles \(\texttt {Tr}_{\texttt {Kdecay}}=\texttt {S}_{\texttt {1}} {\cdot }{} \texttt {S}_{\texttt {2}}{\cdot }{} \texttt {S}_{\texttt {3}} {\cdot }{} \texttt {S}_{\texttt {4}}{\cdot }\overline{\check{\texttt {C}}}_{\texttt {1}} {\cdot }\check{\texttt {C}}_{\texttt {2}}{\cdot }\overline{\texttt {S}}_{\texttt {bk}}\). Then the first trigger requires an energy deposition in GAMS-2000 e.m. calorimeter higher than 2.5 GeV to suppress the dominating \(K^{+}\rightarrow \mu ^{+}\nu \) decay: \(\texttt {Tr}_{\texttt {GAMS}} =\texttt {Tr}_{\texttt {Kdecay}} \cdot (\texttt {E}_{\texttt {GAMS}}>\texttt {2.5})\). The second one, \(\texttt {Tr}_{\texttt {HODO}}=\texttt {Tr}_{\texttt {Kdecay}}{\cdot } {(\texttt {2} \le \texttt {Mult} \le \texttt {4})}\), includes additionally a requirement on multiplicity in the Matrix hodoscope, which allows to efficiently select events with 3 charged particles. The beam intensity, counted by the coincidence of 4 counters (\(\texttt {S}_{\texttt {1}}\cdot \texttt {S}_{\texttt {2}} \cdot \texttt {S}_{\texttt {3}}\cdot \texttt {S}_{\texttt {4}}\)) was \(\sim 2\times 10^{6}\) per spill, the fraction of kaons in the beam was \(\sim \ 12.5\%\), i.e. the kaon intensity was \(\sim 250\,\hbox {k}/\hbox {spill}\). The total number of kaons entering the DV corresponds to \(\sim 3.4\times 10^{10}\).
The Monte Carlo statistics was generated with Geant3 [14] program which includes a realistic description of the setup. The Monte Carlo events are passed through full OKA reconstruction procedures. Monte Carlo sample for the signal events is using as the input the weight, proportional to the square of the absolute value of the leading order matrix element for the decay \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \), given in [7]. For the estimation of the background contribution to the selected data set, a sample of the Monte Carlo events with six main decay channels of charged kaon ( \(\pi ^+\pi ^0\), \(\pi ^+\pi ^0\pi ^0\), \(\pi ^+\pi ^-\pi ^+\), \(\mu ^+\nu \), \(\pi ^0\mu ^+\nu \), \(\pi ^0e^+\nu \)) mixed accordingly to the branching fractions, with the total statistics about equal to that of the recorded data sample is used. The weights for the 3-body decays are calculated from the data, presented in PDG [2]. The same sample serves for the \(K\rightarrow 3\pi \) normalization channel efficiency calculations.
3.1 Event selection
To select \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) decay channel in off-line analysis a set of requirements is applied:
-
the momentum of the beam track is required to be measured;
-
the number of the secondary charged tracks is required to be equal to three and the net charge of them is equal to \(+1\);
-
the decay vertex should have good \(\chi ^2\) and should be with a margin inside the DV;
-
the charged tracks are not identified as electrons in the GAMS-2000 electromagnetic calorimeter;
-
the missing mass squared to each positive pion \(M^{2}_{miss}(\pi ^+)=(P_{K^{+}}-P_{\pi ^{+}})^2\) should be greater than \(0.07\ \hbox {GeV}^2\);
-
the event should contain one and only one photon with the energy \(E_\gamma >0.5\ \hbox {GeV}\);
-
the invariant mass of the photon with each pion (\(M(\pi \gamma )\)) should be greater than 0.17 GeV;
-
the square of the transverse momentum of the \(3\pi \gamma \) system is less than \(0.001\ \hbox {GeV}^2\);
-
the ratio of the momentum of the \(3\pi \gamma \) system to the momentum of the beam track should be within the range 0.95–1.05.
The cut on the square of the missing mass to each positive pion \(M^{2}_{miss}(\pi ^+)\) is used for the suppression of the background from the decay \(K^{+}\rightarrow \pi ^{+}\pi ^{0}\) with \(\pi ^{0}\rightarrow e^{+}e^{-}\gamma \). The main source of the background for the \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) decay is the decay of kaon to three charged pions when the pions produce hadron showers in the electromagnetic calorimeter and, because of fluctuations, a part of a shower is not associated with a charged track by the reconstruction program. To suppress such kind of background a cut on \(M(\pi \gamma )\) is done. The reconstructed momentum of the \(3\pi \gamma \) system should be equal to the momentum of the beam particle. This motivates the last two cuts listed above. For the illustration, the distributions of the square of the transverse momentum of the \(3\pi \gamma \) system and the ratio of absolute value of its momentum to the momentum of the beam particle after all the previous cuts is shown in Figs. 2 and 3. The invariant mass distribution of the \(3\pi \gamma \) system after application of all the cuts listed above is shown in Fig. 4. We see clear separation of the signal, peaking around the nominal value of the kaon PDG [2] mass and the background concentrating at higher masses. The number of events in the signal region is about 450. As almost all the photons in the selected \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) events have the energy in the kaon rest frame greater than 0.03 GeV, in the following analysis we apply an additional explicit cut \(E_{\gamma }^{*} > 0.03\ \hbox {GeV}\).
3.2 Measurement of the branching ratio of the \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) decay
The invariant mass distribution of the \(3\pi \gamma \) system in the signal region is shown in Fig. 5. The Monte Carlo distributions for the main sources of the background and the signal events are shown together with the data. We define the mass range of 0.486–0.504 GeV as a signal region. The number of the decays of \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) is determined as the difference of the number of the data events in the signal region and the expected number of events from the background. The expected background contribution to the signal region is determined using Monte Carlo events for six main channels of the kaon decays described in the above section. To determine the branching fraction of the decay \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) we use the decay \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\) as the normalization channel. It is triggered with the same multiplicity trigger (\(\mathtt{Tr_{HODO}}\)) as the signal and selected with the same criteria as for the pions of the \(3\pi \gamma \) decay, thus we expect the cancellation of many systematic effects related to pions in the ratio of the branching fractions \(\hbox {Br}(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma )/ \hbox {Br}(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+})\). The invariant mass spectrum of the system of three charged pions for the same selection criteria as for the charged part of the \(3\pi \gamma \) decay is shown in Fig. 6.
The branching fraction of the decay \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) is determined in the following way:
where \(Br(3\pi )_{PDG}\) is the branching fraction of the decay \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\) from PDG [2]; \(\epsilon (3\pi )\) and \(\epsilon (3\pi \gamma )\) are the reconstruction efficiencies for the decays \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\) and \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \), determined from the Monte Carlo; \(N(3\pi \gamma )_{D}\) and \(N(3\pi )_{D}\) are the numbers of the decays of \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) and \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\) in the data. Efficiencies \(\epsilon (3\pi )\) and \(\epsilon (3\pi \gamma )\), determined from the Monte Carlo, are \(0.12\pm 0.002\) and \(0.024\pm 0.001\), respectively. The obtained result is
The main source of the systematic error is the uncertainty in the estimate of the background contribution to the signal region (see Figs. 4 and 5). For the estimation of this uncertainty we varied the normalization of the background distribution. This leads to the systematic uncertainty of about \(0.027 \times 10^{-5}\) in the branching fraction. Another source of systematic error is cut’s related. To estimate that, the cuts were varied within experimental resolution around nominal value and the branching was recalculated. The maximal change in the branching was considered as the systematic error. Adding all the systematic errors in quadrature gives a total systematic error of \(0.03 \times 10^{-5}\). The final result is:
The theory prediction in the framework of the CHPT is \(Br_{th}=(0.665\pm 0.005)\times 10^{-5}\) [8].
3.3 Measurement of the differential branching fraction over photon energy
To perform a measurement of the differential branching fraction over photon energy, we split the data sample into four parts with photon energies in the rest frame of the kaon lying in the ranges 0.03–0.04, 0.04–0.05, 0.05–0.06 and 0.06–0.07 GeV. After that we apply the procedure of the previous section to each of the four data samples. The \(3\pi \gamma \) invariant mass spectra for the listed above ranges of the photon energy are shown in Fig. 7 together with expected signal and background contributions.
The obtained values of the branching fractions for the given energy ranges are listed in the Table 1 and shown in the Fig. 8. For comparison, we also give the \(\hbox {O}(p^{4})\) CHPT predictions from [8].
4 Search for a photon up-down asymmetry
In Fig. 9 we show the distribution of the cosine of the angle of the photon direction with respect to the pion system decay plane in the 3-pion rest frame:
where \(n_\gamma \) is the unit vector of photon direction in the 3-pion rest frame, \([p_{f}(\pi ) \times p_{s}(\pi )]\) is the vector product of the momenta of the fastest and slowest pions in the same frame. For comparison, we show the same distribution for the Monte Carlo signal events with the matrix element from [8]. The observed asymmetry in the data is:
The main source of the systematic error comes from the experimental resolution in \(cos(\theta )\). The alternative choice to define the decay plane is to use the fastest and slowest \(\pi ^{+}\). This results in \( A^{+}= -0.04 \pm 0.05(stat.) \pm 0.03(syst.)\). In both cases the observed asymmetry is consistent with zero within errors.
5 Conclusions
The decay \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) is studied on statistics of 450 events. The measured branching fraction is \(Br(3\pi \gamma ) = (0.71\pm 0.05)\times 10^{-5}\) for \(E_\gamma ^{*} > 0.03\ \hbox {GeV}\). The differential branching fraction over photon energy is also determined. The measured branching fractions agree well within the errors with the calculations in the framework of the chiral perturbation theory. The measured up-down asymmetry of the photon with respect to the decay plane of the hadronic system is \(0.03\pm 0.06\). No sign of P and T-odd effects is observed within present statistics.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: OKA data will be released after the intended period of operation, and corresponding publication of results.]
References
V.V. Barmin et al., Sov. J. Nucl. Phys. 50, 421–423 (1989)
C. Patrignani et al., Chin. Phys. C 40, 100001 (2016)
S. Weinberg, Phys. A 96, 327–340 (1979)
J. Gasser, H. Leutwyler, Ann. Phys. (N.Y.) 158, 142 (1984)
J. Gasser, H. Leutwyler, Nucl. Phys. B 250, 465 (1985)
H. Leutwyler, Ann. Phys. (N.Y.) 235, 165 (1994)
G. D’Ambrosio, G. Ecker, G. Isidori, H. Neufeld, Phys. Lett. B 380, 165–170 (1996)
G. D’Ambrosio, G. Ecker, G. Isidori, H. Neufeld, Z. Phys. C 76(2), 301–310 (1997)
R. Aaij et al. (LHCb Collaboration), Phys. Rev. Lett., 112, 1–8 (2014), (Article 161801)
M. Gronau, Y. Grossman, D. Pirjol, A. Ryd, Phys. Rev. Lett. 88, 051802 (2002)
V.I. Garkusha et al., IHEP preprint, IHEP (2003–4)
S.V. Donskov et al., Instrum. Exp. Tech. 59(4), 519–526 (2016)
A.S. Sadovsky et al., Eur. Phys. J. C 78, 92 (2018)
R. Brun et al., GEANT 3 Manual, CERN Program Library Long Writeup W5013, (1994)
Acknowledgements
We express our gratitude to our colleagues in the accelerator department for the good performance of the U-70 during data taking; to colleagues from the beam department for the stable operation of the 21K beam line, including RF-deflectors, and to colleagues from the engineering physics department for the operation of the cryogenic system of the RF-deflectors. The work is supported by the Russian Fund for Basic Research, Grant N18-02-00179A.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funded by SCOAP3
About this article
Cite this article
Shapkin, M.M., Akimenko, S.A., Artamonov, A.V. et al. Study of the decay \(K^{+}\rightarrow \pi ^{+}\pi ^{-}\pi ^{+}\gamma \) in the OKA experiment. Eur. Phys. J. C 79, 296 (2019). https://doi.org/10.1140/epjc/s10052-019-6797-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-019-6797-1