# PHOKHARA 7.0 Monte Carlo generator:
the narrow resonances implementation and new pion and kaon form factors.^{†}^{†}thanks:
Supported in part by Polish Ministry of Science and High Education
from budget for science for years 2010-2013: grant number N N202 102638 and
the European Community-Research
Infrastructure Integrating Activity
“Study of Strongly Interacting Matter”
(acronym HadronPhysics2, Grant Agreement n. 227431)
under the Seventh Framework Program of EU.

###### Abstract

Experiments at high luminosity electron-positron colliders allow to study the kaon and the pion form factors in the time-like region up to high energies. Also the kaon and the pion pair production at and around the narrow resonances and can be investigated, with the interference between electromagnetic and hadronic amplitudes as one of the most interesting phenomenas. Parameterisations of charged and neutral kaon as well as pion form factors, which lead to an improved description of the data in the region of large invariant masses of the meson pair, are presented. They are implemented into the Monte Carlo generator PHOKHARA, together with the hadronic couplings of charged and neutral kaons to and . The physics case as well as details of the implementation are discussed.

## 1 Introduction

During last years meson factories, using the radiative return method [3, 4], contributed a lot to the error reduction of the hadronic contributions to the muon anomalous magnetic moment and the running electromagnetic fine structure constant [1, 2].

An important tool in these analyses is a Monte Carlo generator which simulates all the measured reactions. To meet these requirements, the generator EVA was developed [4, 5], which is based on leading order matrix elements combined with structure function methods for an improved treatment of initial state radiation. To improve the accuracy of these simulations the complete next-to-leading order (NLO) QED corrections were evaluated [6, 7] and implemented into the generator PHOKHARA [8, 9, 10, 11, 12, 13, 14, 15, 16, 17]. For a recent review of theoretical and experimental results see e.g. [18].

In the last paper, resulting in the implementation in the PHOKHARA generator [17], the possibility of studies of the narrow resonances and at B-meson factories was investigated in details profiting from the previous work on the narrow resonances [19]. In this contribution the results of [19, 17] are summarised and complemented with additional material not presented there.

The new models of pion and kaons form factors are presented in Section 2. In Section 3 the modelling of the narrow resonances is shortly sketched. In Section 4 the implementation of the narrow resonances in PHOKHARA is described and a discussion of the importance of the FSR NLO radiative corrections is presented. A short summary follows in Section 5.

## 2 Modeling of the pion and kaon form factors

For a realistic generation of the processes , and models for the electromagnetic pion and kaons form factors are required. In PHOKHARA 5.0 [14, 13] and 6.0 [15] models presented in [20] were implemented. They were published before the CLEO- c measurement of the form factor in the vicinity of the resonance [21] and underestimate the experimental results significantly. Similarly the model predictions at are lower than the pion form factor calculated in [22] from and decay rates. In [17] an effort was made to get better description of the form form factors, especially at high meson pairs invariant masses. In combination with the and hadronic amplitudes parameterisations, obtained in [19] on the basis of new experimental data, it allows for a description of the pion and the kaon pairs production reflecting all existing experimental data.

The models proposed in [17] are generalisations of the models used in [20] and for both, pions and kaons, they are inspired by the dual QCD model of the pion form factor [23]. The infinite towers of , and radial excitations, present in this models, are essential for getting the right behaviour of the form factors at high invariant masses of the meson pairs.

The ansatz for the pion form factor reads

(1) | |||||

with the parameters of the first six radial excitations fitted to the data and others taken as assumed in the original model [23]. For details we refer the reader to [17]. In the original model the coupling constants are real, but it is impossible to fit the data using real couplings, even if one allows that they are different from their model values. Thus they were allowed to be complex, keeping the normalisation (sum of all couplings ) unchanged. With these assumptions the fit is very good . This might be an indication that the radial excitations decay to the same final states, as in this case the mixing between them can generate complex couplings after diagonalisation of the mass matrix. This subject was not studied in details as the experimental knowledge about the decay modes of the mesons is very poor with the exception of [24].

The ansatz for the kaon form factors read

(2) | |||||

with

(3) |

Two versions of the model were used, one with an additional assumption (the ’constrained’ model) and one without this assumption (the ’unconstrained’ model). , , indicate the highest radial resonance used in the fit. For the resonances above these values their model parameters were used. , for the constrained model and , , for the unconstrained model. The unconstrained model fits the data better [17], but the quality of the second fit is also reasonable. It is interesting to observe that the models give distinct predictions in the region were there are no data available (Fig. 1), even if in the region were data are available their behaviour is very similar.

It is also of interest to see how much important is the infinite tower of the resonances in the form factors. In Fig. 2 the relative difference of the form factors without (’no tail’) and with the infinite tower of resonances is shown. In the latter case only resonance contributions are included. As one could guess, the contributions of the higher resonances to the form factors are more pronounced for higher invariant masses, but they are not negligible even in the lower invariant mass range, where the fitted resonances are situated. A similar behaviour of the pion form factor was observed in [17].

## 3 Description of narrow resonances

The narrow resonance (only , and were considered) contributions were included [17] into the Monte Carlo generator PHOKHARA for the muon, the charged pion and the kaon (both charged and neutral) pair production through the following substitution , where

(4) | |||||

(5) |

For and , (it was assumed that the narrow resonances do not decay directly into and ). The contributions to the kaon pair production are included in the kaon form factor, hence , while for and . The notation and the detailed description of the narrow resonance contribution to the amplitude can be found in [19] together with numerical values of the couplings.

## 4 Phokhara 7.0

The new form factors described in Section 2 and the couplings of , , and pairs to narrow resonances described shortly in Section 3 were implemented into a Monte Carlo event generator PHOKHARA. It’s new release (PHOKHARA 7.0) contains also the implementation of the 4-pion hadronic currents as described in [16], not included in the previous releases of the code. The code can be found on the web page http://ific.uv.es/~rodrigo/phokhara/ , where also a user guide is provided.

The variance reduction through standard change of variables using the Breit-Wigner shape of resonances was enough to allow for a reasonable acceptance rate. Other technical problems with the treatment of the narrow resonances and their solution were discussed at length in [17, 30] and will not be repeated here.

A typical invariant mass spectrum of the final meson or muon pair in the vicinity of a narrow resonance, as seen in a radiative return experiment, resembles more a resolution curve of a detector then the shape of the resonance. A typical picture is shown in Fig. 3, where the generated mass spectrum of the pair was convoluted with a Gaussian with a detector resolution of the BaBar detector [31].

For the resonance the smearing effects to large extend dilute the influence of the resonance on the invariant mass spectrum, as shown in Fig. 4 for the charged kaon pair invariant mass distribution. Thus to study it with the radiative return method fairly large statistics is required and it is difficult to expect that the method can compete with the accuracy of scan experiments.

The contributions from the to all studied invariant mass spectra is sizable and one can study the decays with a good accuracy. Thus a legitimate question is how big are the contributions coming from the final state photon radiation (FSR). If the detector smearing effects are not taken into account, the distortions of the meson (muon) pair invariant mass distribution by the FSR effects, as compared to ISR, are huge as shown in Figs. 5 and 6. They are still sizable when one takes into account the smearing effects, as evident from Figs. 7 and 8. Moreover they do depend on the event selection [30]. However in the integrated cross section to large extend the corrections cancel. For the muon pair production pb and pb, while for the charged kaon pair production nb and nb. Thus when one uses only the integrated spectra the FSR corrections can be safely neglected.

## 5 Summary

The implementation into the PHOKHARA Monte Carlo event generator of the narrow resonance contributions is presented and the possibility of studies of the narrow resonances at B factories is discussed.

## References

- [1] Andreas Hoecker, these proceedings
- [2] Thomas Teubner, these proceedings
- [3] Min-Shih Chen and P. M. Zerwas, Phys. Rev. D 11 (1975) 58.
- [4] S. Binner, J. H. Kühn and K. Melnikov, Phys. Lett. B 459 (1999) 279 [hep-ph/9902399].
- [5] H. Czyż and J. H. Kühn, Eur. Phys. J. C 18 (2001) 497 [hep-ph/0008262].
- [6] G. Rodrigo, A. Gehrmann-De Ridder, M. Guilleaume and J. H. Kühn, Eur. Phys. J. C 22 (2001) 81 [hep-ph/0106132].
- [7] J. H. Kühn and G. Rodrigo, Eur. Phys. J. C 25 (2002) 215 [hep-ph/0204283].
- [8] G. Rodrigo, H. Czyż, J.H. Kühn and M. Szopa, Eur. Phys. J. C 24 (2002) 71 [hep-ph/0112184].
- [9] H. Czyż, A. Grzelińska, J. H. Kühn and G. Rodrigo, Eur. Phys. J. C 27 (2003) 563 [hep-ph/0212225].
- [10] H. Czyż, A. Grzelińska, J. H. Kühn and G. Rodrigo, Eur. Phys. J. C 33 (2004) 333 [hep-ph/0308312].
- [11] H. Czyż, J. H. Kühn, E. Nowak and G. Rodrigo, Eur. Phys. J. C 35 (2004) 527 [hep-ph/0403062].
- [12] H. Czyż, A. Grzelińska, J. H. Kühn and G. Rodrigo, Eur. Phys. J. C 39 (2005) 411 [hep-ph/0404078].
- [13] H. Czyż, A. Grzelińska and J. H. Kühn, Phys. Lett. B 611 (2005) 116 [hep-ph/0412239].
- [14] H. Czyż, A. Grzelińska, J. H. Kühn and G. Rodrigo, Eur. Phys. J. C 47 (2006) 617 [arXiv:hep-ph/0512180].
- [15] H. Czyż, A. Grzelińska and J. H. Kühn, Phys. Rev. D 75, 074026 (2007) [arXiv:hep-ph/0702122].
- [16] H. Czyż, J. H. Kühn and A. Wapienik, Phys. Rev. D 77, 114005 (2008) [arXiv:0804.0359 [hep-ph]].
- [17] H. Czyż, A. Grzelińska and J. H. Kühn, Phys. Rev. D 81 (2010) 094014 [arXiv:1002.0279 [hep-ph]].
- [18] S. Actis et al., arXiv:0912.0749 [hep-ph].
- [19] H. Czyż and J. H. Kühn, Phys. Rev. D 80 (2009) 034035 [arXiv:0904.0515 [hep-ph]].
- [20] C. Bruch, A. Khodjamirian and J.H. Kühn, Eur. Phys. J. C 39 (2005) 41, [hep-ph/0409080].
- [21] T. K. Pedlar et al. [CLEO Collaboration], Phys. Rev. Lett. 95, 261803 (2005) [arXiv:hep-ex/0510005].
- [22] J. Milana, S. Nussinov and M. G. Olsson, Phys. Rev. Lett. 71, 2533 (1993) [arXiv:hep-ph/9307233].
- [23] C. A. Dominguez, Phys. Lett. B 512, 331 (2001) [arXiv:hep-ph/0102190].
- [24] K. Nakamura et al. [ Particle Data Group Collaboration ], J. Phys. G G37 (2010) 075021.
- [25] F. Mane, D. Bisello, J. C. Bizot, J. Buon, A. Cordier and B. Delcourt, Phys. Lett. B 99, 261 (1981).
- [26] M. N. Achasov et al., Phys. Rev. D 63, 072002 (2001).
- [27] R. R. Akhmetshin et al., Phys. Lett. B 551, 27 (2003) [arXiv:hep-ex/0211004].
- [28] R. R. Akhmetshin et al. [CMD-2 Collaboration], Phys. Lett. B 578, 285 (2004) [arXiv:hep-ex/0308008].
- [29] M. N. Achasov et al., J. Exp. Theor. Phys. 103, 720 (2006) [Zh. Eksp. Teor. Fiz. 103, 831 (2006)] [arXiv:hep-ex/0606057].
- [30] H. Czyż and A. Grzelińska, Chinese Physics C 34 (2010) 855.
- [31] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 69 (2004) 011103 [arXiv:hep-ex/0310027].