## RECENT NA48/2 AND NA62 RESULTS

## D. MADIGOZHIN a

Joint Institute for Nuclear Research, Dubna, Russia E-mail: madigo@mail.cern.ch

The NA48/2 Collaboration at CERN has accumulated and analysed unprecedented statistics of rare kaon decays in the K e 4 modes: K e 4 (+ -) ( K ± → π + π - ± e ν ) and K e 4 (00) ( K ± → π π e 0 0 ± ν ) with nearly one percent background contamination. It leads to the improved measurement of branching fractions and detailed form factor studies. New final results from the analysis of 381 K ± → π ± γγ rare decay candidates collected by the NA48/2 and NA62 experiments at CERN are presented. The results include a decay rate measurement and fits to Chiral Perturbation Theory (ChPT) description.

## 1 Introduction

The main goal of NA48/2 experiment was the search for CP-violating asymmetry in K ± → 3 π ± decays 1 . Also it has provided in 2003-2004 a large data sample for charged kaon rare decay studies. In 2007-2008, the NA62 experiment 2 ( R K phase) has collected a large data sample with the same detector 3 but modified beam line. Two simultaneous K + and K -beams were produced by 400 GeV/ c protons on a beryllium target. Particles of opposite charge with a central momentum of 60 GeV/ c and a momentum band of ± 3 8% ( . rms ) (74 GeV/ c ± 1 9% for . NA62 R K phase) were selected by the system of magnets and collimators. Charged products of K ± decays were measured by the magnetic spectrometer consisting of four drift chambers (DCH1-DCH4) and a dipole magnet located between DCH2 and DCH3. The spectrometer was

a for the NA48/2 and NA62 Collaborations: F. Ambrosino, A. Antonelli, G. Anzivino, R. Arcidiacono, W. Baldini, S. Balev, J.R. Batley, M. Behler, S. Bifani, C. Biino, A. Bizzeti, B. Bloch-Devaux, G. Bocquet, V. Bolotov, F. Bucci, N. Cabibbo, M. Calvetti, N. Cartiglia, A. Ceccucci, P. Cenci, C. Cerri, C. Cheshkov, J.B. Ch` eze, M. Clemencic, G. Collazuol, F. Costantini, A. Cotta Ramusino, D. Coward, D. Cundy, A. Dabrowski, G. D'Agostini, P. Dalpiaz, C. Damiani, H. Danielsson, M. De Beer, G. Dellacasa, J. Derr´ e, H. Dibon, D. Di Filippo, L. DiLella, N. Doble, V. Duk, J. Engelfried, K. Eppard, V. Falaleev, R. Fantechi, M. Fidecaro, L. Fiorini, M. Fiorini, T. Fonseca Martin, P.L. Frabetti, A. Fucci, S. Gallorini, L. Gatignon, E. Gersabeck, A. Gianoli, S. Giudici, A. Gonidec, E. Goudzovski, S. Goy Lopez, E. Gushchin, B. Hallgren, M. Hita-Hochgesand, M. Holder, P. Hristov, E. Iacopini, E. Imbergamo, M. Jeitler, G. Kalmus, V. Kekelidze, K. Kleinknecht, V. Kozhuharov, W. Kubischta, V. Kurshetsov, G. Lamanna, C. Lazzeroni, M. Lenti, E. Leonardi, L. Litov, D. Madigozhin, A. Maier, I. Mannelli, F. Marchetto, G. Marel, M. Markytan, P. Marouelli, M. Martini, L. Masetti, P. Massarotti, E. Mazzucato, A. Michetti, I. Mikulec, M. Misheva, N. Molokanova, E. Monnier, U. Moosbrugger, C. Morales Morales, M. Moulson, S. Movchan, D.J. Munday, M. Napolitano, A. Nappi, G. Neuhofer, A. Norton, T. Numao, V. Obraztsov, V. Palladino, M. Patel, M. Pepe, A. Peters, F. Petrucci, M.C. Petrucci, B. Peyaud, R. Piandani, M. Piccini, G. Pierazzini, I. Polenkevich, I. Popov, Yu. Potrebenikov, M. Raggi, B. Renk, F. Reti`re, P. Riedler, e A. Romano, P. Rubin, G. Ruggiero, A. Salamon, G. Saracino, M. Savri´ e, M. Scarpa, V. Semenov, A. Sergi, M. Serra, M. Shieh, S. Shkarovskiy, M.W. Slater, M. Sozzi, T. Spadaro, S. Stoynev, E. Swallow, M. Szleper, M. Valdata-Nappi, P. Valente, B. Vallage, M. Velasco, M. Veltri, S. Venditti, M. Wache, H. Wahl, A. Walker, R. Wanke, L. Widhalm, A. Winhart, R. Winston, M.D. Wood, S.A. Wotton, O. Yushchenko, A. Zinchenko, M. Ziolkowski.

followed by a scintillator hodoscope. A Liquid Krypton calorimeter (LKr) was used to measure the energy of electrons and photons.

## 2 K + -e 4 decay

Kinematics of the K ± → π + π - ± e ν ( K + -e 4 ) decay is defined by five variables 4 : the squared invariant masses of dipion ( S π ) and dilepton ( S e ), the angle θ π of π ± in dipion rest frame with respect to the flight direction of dipion in the kaon center of mass system, the similar angle θ e of e ± in dilepton rest frame, and the angle φ between dipion and dilepton planes.

K e 4 decay amplitude is a product of the leptonic weak current and (V-A) hadronic current, described in terms of three (F,G,R) axial-vector and one (H) vector complex form factors. Sensitivity of K e 4 decay matrix element to R form factor is negligible due to the small mass of electron. Form-factors may be developed in a partial wave expansion:

$$F = F _ { s } e ^ { i \delta _ { f _ { s } } } + F _ { p } e ^ { i \delta _ { f _ { p } } } & c o s \theta _ { \pi } + F _ { d } e ^ { i \delta _ { f _ { d } } } c o s ^ { 2 } \theta _ { \pi } + \dots \\ G = G _ { p } e ^ { i \delta _ { g _ { p } } } + G _ { d } e ^ { i \delta _ { g _ { d } } } c o s \theta _ { \pi } + \dots \\ H = H _ { p } e ^ { i \delta _ { h _ { p } } } + H _ { d } e ^ { i \delta _ { h _ { d } } } c o s \theta _ { \pi } + \dots$$

Limiting the expansion to S- and P-waves and considering a unique phase δ p for all P-wave form factors in absence of CP violating weak phases, one will obtain the decay probability, that depends only on the real form factor magnitudes F , F s p , G p , H p , a single phase shift δ = δ s -δ p and kinematic variables. The form factors can be developed in a series expansion of the dimensionless invariants q 2 = ( S / m π 4 2 π ) -1 and S / m e 4 2 5 π . Two slope and one curvature terms are sufficient to describe the F s form factor variation within the available statistics ( F /f s s = 1+( f /f ′ s s ) q 2 +( f ′′ s /f s ) q 4 +( f /f ′ e s ) S / m e 4 2 π ), while two terms are enough to describe the G p form factor ( G /f p s = g /f p s +( g /f ′ p s ) q 2 ), and two constants - to describe the F p and H p form factors.

Hadronic form factors for the S- and P-waves have been obtained by NA48/2 concurrently with the phase difference between the S- and P-wave states of ππ system, leading to the precise determination of a 0 0 and a 0 2 , the I=0 and I=2 S-wave ππ scattering lengths 7 .

A high precision measurement of K + -e 4 form factors and branching fraction has been published by NA48/2 few years later 8 9 , . K + -e 4 decay rate was measured relative to K ± → π + π -π ± ( K + -3 π ) normalization channel. A track of charged particle with a momentum p &gt; 2 75 GeV . and 0 9 . &lt; E/p &lt; 1 1 was identified as . e ± , while the track with p &gt; 5 GeV and E/p &lt; 0 8 was . regarded as π ± . A dedicated linear discriminant variable based on shower properties has been applied to reject events with one misidentified pion. To suppress K + -3 π background, the vertex invariant mass M 3 π in the π + π -π ± hypothesis and its transverse momentum p t were required to be outside an ellipse centered at PDG kaon mass 10 and zero transversal momentum, with semi-axes of 20 MeV/c 2 and 35 MeV/c , respectively.

The squared missing mass was required to be &gt; 0 04 ( . Gev/c 2 ) 2 to reject π ± π 0 decays with a subsequent π 0 → e + e -γ process. The invariant mass of e + e -system was required to be &gt; 0 03 . GeV/c 2 in order to reject photon conversions.

For the normalization channel K + -3 π the M 3 π and p t were inside a smaller ellipse with semiaxes 12 MeV/c 2 and 25 MeV/c , respectively. A sample of about 1.11 million K + -e 4 candidates and about 19 millions of prescaled K 3 π candidates were selected from data recorded in 2003-2004.

Two main background sources are known: K ± → π + π -π ± decays with subsequent π → eν decay or a pion mis-identified as an electron; and K ± → π 0 ( π 0 ) π ± with subsequent π 0 → e + e -γ decay with undetected photons and an electron mis-identified as a pion. Their admixture in the signal events is estimated to be below 1%.

A detailed GEANT3-based 11 Monte Carlo simulation was used to take into account full detector geometry, DCH alignment, local inefficiencies and beam properties.

The resulting K + -e 4 branching fraction 9 BR K ( + -e 4 ) = (4 257 . ± 0 004 . stat ± 0 016 . syst ± 0 031 . ext )10 -5 is 3 times more precise than available PDG value 10 . It has been used to extract the common normalization form factor f s 9 .

## 3 K 00 e 4 decay

The K 00 e 4 rate is measured relative to the K ± → π π π 0 0 ± ( K 00 3 π ) normalization channel. These two modes are collected using the same trigger and with a similar event selections. The separation between them occurs only at a later stage.

Events with at least four γ , detected by LKr, and at least one track, reconstructed from spectrometer data, were regarded as K 00 e 4 or K 00 3 π candidates. Every combination of 4 reconstructed γ with energies E &gt; 3 GeV was considered as a possible pair of π 0 decays. Reconstructed longitudinal positions Z 1 and Z 2 of both π 0 → 2 γ decay candidates were required to coincide within 500 cm , with their average position Z n = ( Z 1 + Z 2 ) / 2 inside the fiducial volume 106 m long.

Decay longitudinal position Z ch , assigned to the track, was defined by the closest distance approach between the track and the beam axis. Combined vertex, composed of four LKr clusters and one charged track with momentum p &gt; 5 GeV , was required to have the difference | Z n -Z ch | less than 800 cm . If several combinations satisfy the vertex criteria, the case of minimum ( Z 1 -Z 2 σ n ) 2 +( Z n -Z ch σ c ) 2 has been chosen, where σ n and σ c are the Z n -dependent widths of corresponding distributions.

A track was preliminarily identified as e ± , if it has an associated LKr cluster with E/p between 0.9 and 1.1, otherwise π ± was assumed at the first stage. Further suppression of pions mis-identified as electrons is obtained by means of discriminant variable which is a linear combination of E/p , shower width and energy weighted track-to-cluster distance at LKr front face.

K 00 e 4 and K 00 3 π decays were discriminated by means of elliptic cuts in the ( M π π π 0 0 ± , p t ) plane, where M π π π 0 0 ± is the invariant mass of combined vertex in the K 0 3 π hypothesis, and p t is the transversal momentum. Elliptic cut separates about 94 million K 00 3 π normalization events from about 65000 K 00 e 4 candidates. Residual fake-electron background is about 0.65% of K 00 e 4 amount. Background from K 00 3 π with the subsequent π ± → e ± ν is 0.12% of the signal, and the accidentalrelated background is about 0.23%. It gives in total 1% of background admixture.

For the case of K ± → π π e 0 0 ± ν ( K 00 e 4 ) decay, due to restrictions of symmetry, matrix element doesn't depend on θ π and φ angles. It is parametrized in terms of the only formfactor F s , that may depend on S π and S e . Form factor F s was extracted from the fit of events distribution on ( S , S e π ) plane, taking into account the acceptance, calculated from MC simulation. The following empirical parameterization has been used: F /f s s = 1 + ( f /f ′ s s ) q 2 +( f ′′ s /f s ) q 4 +( f /f ′ e s ) S / m e 4 2 π for q 2 &gt; 0 and F /f s s = 1 + d √ | q 2 / (1 + q 2 ) | +( f /f ′ e s ) S / m e 4 2 π for q 2 &lt; 0.

The results are in agreement with NA48/2 K + -e 4 analysis described above: f /f ′ s s = 0 149 . ± 0 033 . stat ± 0 014 . syst , f ′′ s /f s = -0 070 . ± 0 039 . stat ± 0 013 . syst , f /f ′ e s = 0 113 . ± 0 022 . stat ± 0 007 . syst , d = -0 256 . ± 0 049 . ± 0 016 . syst .

The obtained form factor was used to obtain the final result of branching fraction measurement: Br K ( 00 e 4 ) = (2 552 . ± 0 010 . stat ± 0 010 . syst ± 0 032 . ext )10 -5 . It is 10 times more precise, than PDG corresponding value 10 . Systematic error includes the contributions from background, simulation statistical error, sensitivity to form factor, radiation correction, trigger efficiency and beam geometry. External error comes from uncertainty of normalization channel K 00 3 π branching fraction.

Below the threshold of S π = (2 m π ± ) 2 the measured K 00 e 4 decay form factor shows a deficit of events, that is well described by the present empirical parameterization (Fig. 1). It is similar to the effect of π + π -→ π π 0 0 rescattering in K ± → π π π 0 0 ± decay (cusp effect 6 ), investigated by NA48/2 collaboration earlier 7 on the basis of ChPT formulations.

Figure 1 K 00 e 4 normalized form factor squared as a function of q 2 . The line corresponds to the adopted empirical fit. The arrow points to the 2 m π threshold.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0585v1-figure-44cf31c2-13e5-4c27-83a2-2a5c599fdd33.jpg)

## 4 K ± → π ± γγ decay

In the ChPT framework, the K ± → π ± γγ decay receives two non-interfering contributions at lowest non-trivial order O ( p 4 ): the pion and kaon loop amplitude depending on an unknown O (1) constant ˆ representing the total contribution of the counterterms, and the c pole amplitude 12 .

New measurements of this decay have been performed using data collected during a 3-day special NA48/2 run in 2004 and a 3-month NA62 run in 2007. Signal events are selected in the region of z = ( m γγ /m K ) 2 &gt; 0 2 . to reject the K ± → π ± π 0 background peaking at z = 0 075. . 149 (232) decays candidates are observed in the 2004 (2007) data set, with backgrounds contaminations of 10.4% (7.5%) from K ± → π ± π 0 ( π 0 )( γ ) decays with merged photon clusters in the electromagnetic calorimeter.

The values of ˆ in the frameworks of the ChPT c O ( p 4 ) and O ( p 6 ) parameterizations 13 as well as branching ratio have been measured using likelihood fits to the data. The main systematic effect is due to the background uncertainty. Uncertainties related to trigger, particle identification, acceptance and accidental effects found to be negligible. The final combined results based on 2004 and 2007 runs data 14 15 , are: ˆ for c O ( p 4 ) fit = 1 72 . ± 0 20 . stat ± 0 06 . syst ; ˆ for c O ( p 6 ) fit = 1 86 . ± 0 23 . stat ± 0 11 . syst ; branching fraction Br K ( πγγ ) for O ( p 6 ) fit = (1 003 . ± 0 056) . × 10 -6 . The model-independent branching ratio for z &gt; 0 2 is . equal to (0 965 . ± 0 063) . × 10 -6 . New results are in agreement with the earlier (based on 31 events) BNL E787 16 ones.

## References

- 1. J.R.Batley, et al., Eur. Phys. J. C 52 (2007) 875-891.
- 2. G. Anelli, et al., CERN-SPSC-2005-013.
- 3. V. Fanti, et al., Nucl.Instrum.Meth. A574 (2007) 433-471.
- 4. N. Cabibbo, A. Maksymowicz, Phys.Rev. 137 (1965) B438-B443.
- 5. G. Amoros, J. Bijnens, J.Phys.G G25 (1999) 1607-1622.
- 6. N. Cabibbo, Phys.Rev.Lett. 93 (2004) 121801.
- 7. J. Batley, et al., Eur.Phys.J. C64 (2009) 589-608.
- 8. J. Batley, et al., Eur.Phys.J. C70 (2010) 635-657.
- 9. J. Batley, et al., Phys.Lett. B715 (2012) 105-115.
- 10. J. Beringer, et al., Phys.Rev. D86 (2012) 010001.
- 11. R. Brun, F. Carminati, S. Giani, CERN-W-5013.
- 12. G. Ecker, A. Pich and E. de Rafael, Nucl. Phys. B303 (1988) 665.
- 13. G. D'Ambrosio and J. Portol´s, Phys. Lett. e B386 (1996) 403.
- 14. J. Batley, et al., Phys. Lett. B730 (2014) 141.
- 15. J. Batley, et al., Phys. Lett. B732 (2014) 65.
- 16. P. Kitching et al., Phys. Rev. Lett. 79 (1997) 4079.
F. Ambrosino A. Antonelli G. Anzivino R. Arcidiacono W. Baldini S. Balev J. R. Batley M. Behler S. Bifani C. Biino A. Bizzeti B. Bloch-Devaux G. Bocquet V. Bolotov F. Bucci N. Cabibbo M. Calvetti N. Cartiglia A. Ceccucci P. Cenci C. Cerri C. Cheshkov J. B. Cheze M. Clemencic G. Collazuol F. Costantini A. Cotta Ramusino D. Coward D. Cundy A. Dabrowski G. D'Agostini P. Dalpiaz C. Damiani H. Danielsson M. De Beer G. Dellacasa J. Derre H. Dibon D. Di Filippo L. DiLella N. Doble V. Duk J. Engelfried K. Eppard V. Falaleev R. Fantechi M. Fidecaro L. Fiorini M. Fiorini T. Fonseca Martin P. L. Frabetti A. Fucci S. Gallorini L. Gatignon E. Gersabeck A. Gianoli S. Giudici A. Gonidec E. Goudzovski S. Goy Lopez E. Gushchin B. Hallgren M. Hita-Hochgesand M. Holder P. Hristov E. Iacopini E. Imbergamo M. Jeitler G. Kalmus V. Kekelidze K. Kleinknecht V. Kozhuharov W. Kubischta V. Kurshetsov G. Lamanna C. Lazzeroni M. Lenti E. Leonardi L. Litov D. Madigozhin A. Maier I. Mannelli F. Marchetto G. Marel M. Markytan P. Marouelli M. Martini L. Masetti P. Massarotti E. Mazzucato A. Michetti I. Mikulec M. Misheva N. Molokanova E. Monnier U. Moosbrugger C. Morales Morales M. Moulson S. Movchan D. J. Munday M. Napolitano A. Nappi G. Neuhofer A. Norton T. Numao V. Obraztsov V. Palladino M. Patel M. Pepe A. Peters F. Petrucci M. C. Petrucci B. Peyaud R. Piandani M. Piccini G. Pierazzini I. Polenkevich I. Popov Yu. Potrebenikov M. Raggi B. Renk F. Retiere P. Riedler A. Romano P. Rubin G. Ruggiero A. Salamon G. Saracino M. Savrie M. Scarpa V. Semenov A. Sergi M. Serra M. Shieh S. Shkarovskiy M. W. Slater M. Sozzi T. Spadaro S. Stoynev E. Swallow M. Szleper M. Valdata-Nappi P. Valente B. Vallage M. Velasco M. Veltri S. Venditti M. Wache H. Wahl A. Walker R. Wanke L. Widhalm A. Winhart R. Winston M. D. Wood S. A. Wotton O. Yushchenko A. Zinchenko M. Ziolkowski
Recent NA48/2 and NA62 results
The NA48/2 Collaboration at CERN has accumulated and analysed unprecedented
statistics of rare kaon decays in the $K_{e4}$ modes: $K_{e4}(+-)$ ($K^\pm \to
\pi^+ \pi^- e^\pm \nu$) and $K_{e4}(00)$ ($K^\pm \to \pi^0 \pi^0 e^\pm \nu$)
with nearly one percent background contamination. It leads to the improved
measurement of branching fractions and detailed form factor studies. New final
results from the analysis of 381 $K^\pm \to \pi^\pm \gamma \gamma$ rare decay
candidates collected by the NA48/2 and NA62 experiments at CERN are presented.
The results include a decay rate measurement and fits to Chiral Perturbation
Theory (ChPT) description.