Mini-Proceedings, 16th meeting of the Working Group on Radiative Corrections and MC Generators for Low Energies
18 - 19 November, Laboratori Nazionali di Frascati, Italy
Henryk Czyz (Katowice), Pere Masjuan (Mainz), Graziano Venanzoni (Frascati)
The mini-proceedings of the 16 Meeting of the ”Working Group on Radiative Corrections and MonteCarlo Generators for Low Energies” held in Frascati, 18 - 19 November, are presented. These meetings, started in 2006, have as aim to bring together experimentalists and theoreticians working in the fields of meson transition form factors, hadronic contributions to the anomalous magnetic moment of the leptons, and the effective fine structure constant. The development of MonteCarlo generators and Radiative Corrections for precision and -lepton physics are also covered.
The web page of the conference:
contains the presentations.
We acknowledge the support and hospitality of the Laboratori Nazionali di Frascati.
- 1 Introduction to the Workshop
- A Summaries of the talks
- A.2 Current status of luminosity measurement with the CMD-3 detector at the VEPP-2000 ee collider
- A.3 The role of experimental data as input information for precise hadronic calculations: muon , rare decays, and mixing parameters
- A.4 On the positronium contribution to the electron
- A.5 Hadronic light-by-light scattering in the muon : a dispersive approach
- A.6 Primary Monte-Carlo generator of the process for the CMD-3 experiment
- A.7 Automation of the leading order calculations for hadrons
- A.8 MCGPJ for the processes for experiments with CMD-3 detector at the VEPP-2000 collider
- A.9 and production at colliders - preliminary results
- A.10 Nucleon form factors in PHOKHARA
- A.11 Current status of Monte Carlo generator Tauola
- B List of participants
1 Introduction to the Workshop
H. Czyż and G. Venanzoni
Institute of Physics, University of Silesia, 40007 Katowice, Poland
Laboratori Nazionali di Frascati dellÕINFN, 00044 Frascati, Italy
The importance of continuous and close collaboration between the experimental and theoretical groups is crucial in the quest for precision in hadronic physics. This is the reason why the Working Group on “Radiative Corrections and Monte Carlo Generators for Low Energies” (Radio MonteCarLow) was formed a few years ago bringing together experts (theorists and experimentalists) working in the field of low-energy physics and partly also the community. Its main motivation was to understand the status and the precision of the Monte Carlo generators used to analyze the hadronic cross section measurements obtained as well with energy scans as with radiative return, to determine luminosities, and whatever possible to perform tuned comparisons, i.e. comparisons of MC generators with a common set of input parameters and experimental cuts. This main effort was summarized in a report published in 2010 . During the years the WG structure has been enriched of more physics items and now it includes seven subgroups: Luminosity, R-measurement, ISR, Hadronic VP and Delta alpha, gamma-gamma physics, FSR models, tau.
During the workshop the last achievements of each subgroups have been presented. The present accuracy and the future prospects of MC generators for into leptonic, , and hadronic final states have been reviewed. The recent evaluation of the positronium contribution to the electron and the role of experimental data to the hadronic LO and Light-by-Light NLO contributions to the of the muon have been discussed. New results from CMD3 and BESIII experiments have been presented. Finally the status of HPrecisionNet work package of the networking program HPH to Horizon 2020, was presented.
The workshop was held from the 18 to the 19 November, at the Laboratori Nazionali di Frascati dellÕINFN, Italy.
Webpage of the conference is
where detailed program and talks can be found.
All the information on the WG can be found at the web page:
Appendix A Summaries of the talks
a.1 Present accuracy and future prospects of Monte Carlo generators for Bhabha and
C.M. Carloni Calame
Dipartimento di Fisica, Università di Pavia, Via A. Bassi 6, 27100 Pavia, Italy
The knowledge of the luminosity is a key ingredient for any measurement at machines. The usual strategy to calculate it is through the relation , where is the theoretical cross section of a QED process and the number of events. QED processes are the best choice because of their clean signal, low background and the possibility to push the theoretical accuracy up to the level or better. The latter requires the inclusion of the relevant radiative corrections (RCs) and their implementation into Monte Carlo (MC) event generators (EGs) to easily account for experimental event selection criteria.
Modern EGs used for luminometry simulate Bhabha, and (or a sub-set of them), including the exact NLO QED corrections and/or a leading-log (LL) approximation of higher-order (h.o.) effects [5, 6, 7, 8, 9]. The consistent inclusion of NLO and h.o. LL contributions is mandatory in view of the required theoretical accuracy.
Focusing on Bhabha scattering, in order to estimate the theoretical error of the EGs, it is extremely important to perform tuned comparisons among them, to assess the technical precision and have an idea of the accuracy of the included corrections, usually implemented according to different approaches. This has been done and reported in , from which tab. 1 has been extracted.
In general, it is found that the different MC EGs predict cross sections which differ by at most when including NLO and h.o. LL corrections. This is in fair agreement with the accuracy of the different approaches estimated by the authors.
A further step to put on firmer ground the theoretical error is to compare with exact NNLO results, which have been calculated for Bhabha scattering by various groups in the last years (see references in ). NNLO RCs are partly included in EGs and, once extracted, their NNLO contributions can be consistently compared with exact calculations. This has been done for example in the paper of  and in .
In Tab. 2 (adapted from ) the total error “budget” is summarized for typical luminometry conditions at flavor factories. The main conclusion to draw from it is that RCs currently implemented into MC EGs allow to reach a theoretical precision up to the level.
|Source of error (%)||factories||= 3.5 GeV||factories|
It has to be mentioned that the error induced by vacuum polarization (VP) corrections is driven and dominated by experimental errors. At energies around the narrow resonances (such as ), VP errors might be larger than in tab. 2 and a dedicated study is needed.
A similar picture ought to be valid also for final state, with the added advantage that, at least up to NLO, VP corrections do not contribute to the cross section. Nevertheless, a careful estimate of the theoretical error in this case has not been performed yet and would be of high interest.
A possible improvement of the theoretical accuracy, if going beyond the level is needed at all, would be the inclusion of the full NNLO results into the MC EGs, which is a non trivial but feasible task.
I’d like to thank H. Czyz and G. Venanzoni for the kind invitation and the organization of a really stimulating workshop.
a.2 Current status of luminosity measurement with the CMD-3 detector at the VEPP-2000 ee collider
Budker Institute of Nuclear Physics, SB RAS, Novosibirsk, 630090, Russia
Novosibirsk State University, Novosibirsk, 630090, Russia
Since December 2010 the CMD-3  detector has taken data at the electron-positron collider VEPP-2000 . The collected data sample corresponds to an integrated luminosity of 60 pb in the c.m. energy from 0.32 up to 2 GeV.
The luminosity is a key part in many experiments which study the hadronic cross sections at colliders. As a rule, the systematic error of the luminosity determination represents one of the largest sources of uncertainty which can cause significant reduction of the hadronic cross sections accuracy. Therefore it is very important to have several well known QED processes such as , , to determine the luminosity. The combined application of them will help to better understand and estimate a real systematic accuracy of the luminosity. The combined application of them will help to better understand and estimate a real systematic accuracy of the luminosity.
The energy range from 1 to 2 GeV was scanned up and down with a step of 50 MeV. At each energy point the integrated luminosity 500 nb was collected. During the scan down the energy points, at which the data were collected, have been shifted to the previous one by 25 MeV. The data were collected at an average luminosity scm. At the highest energies the peak luminosity reached the values about scm and was restricted by the positron storage rate in the booster. The project luminosity scm will be provided only with start of operating of new positron injection facility in 2015. The beam energy has been monitored ( 0.5 MeV) by measuring the current in dipole magnets of the main ring. The time period of this run was extend from January to June 2011. In 2012 the luminosity was measured at 16 energy points from 1.32 GeV to 1.98 GeV and collected luminosity was about 14 pb.
In 2013 the energy range from 0.32 GeV to 1 GeV was scan with the 10 MeV step. The integrated luminosities about 8.3 and 8.4 pb were collected around and mesons. Over the 2013 year the integrated luminosity 25 pb has been collected.
The sample of collinear events , , , and cosmic background were selected for luminosity determination. The two-dimensional plot of energy deposition in calorimeters for these events is presented in Fig.2 for the beam energy 950 MeV. It is clear seen that Bhabha events are distributed predominantly at the upper right corner whereas other particles are concentrated in the bottom left one. Thus, the integrated luminosity can be determined by the well selected Bhabha events.
To select events the information about they energy deposition in calorimeters is used. It is seen that the signal events are concentrated as a cluster of dots in upper-right corner of this plot. At the same time two train seen - concentration of dots in two mutually perpendicular directions due to ISR. The events of this sample should have energy deposition inside interval: . Unfortunately the small part of the Bhabha events can seep under the central peak and imitate events. To exclude such events the additional condition was applied - the Z-chamber sectors associated with clusters must be triggered. Visual scan of the remaining events proofed - there are not Bhabha events under central. Unfortunately this condition delete some events due to albedo coming from showers. The fraction of such events amounts to and as a result we should include correction about to restore the number of events.
The luminosities ratio determined with use of two processes vs energy is presented in Fig. 2 and in Fig. 4, where only statistical errors are shown. The blue circles correspond to the scan up, whereas red circles - scan down. The horizontal line is a fit for this ratio for scan down. In this case the relative difference between luminosities is in average %. However, at the beginning of the run the difference was 3% and explained by hardware problems and the quality of inter-calibration of the detector subsystems. Collecting all the main sources which contribute to systematic error of the luminosity, we estimate the current accuracy as 1% while. The first energy scan below 1 GeV was performed at VEPP-2000 during the season of 2013. The preliminary results of the luminosity measurement are shown in Fig. 4. The already collected statistics is higher than that in the previous CMD-2 experiment and at the level or better than in BaBar and KLOE experiments. One of the tests in this analysis is to measure the cross section of the process at low energy, where particles separation is possible using only momentum information from DC. Preliminary results of this test are consistent with the QED prediction as it is seen in Fig. 4. The radiative corrections (RC) to this cross section with photon jets radiation in collinear regions were taken into account according to  and their accuracy is better than .
a.3 The role of experimental data as input information for precise hadronic calculations: muon , rare decays, and mixing parameters
PRISMA Cluster of Excellence, Institut für Kernphysik, Johannes Gutenberg-Universtät, Mainz D-55099, Germany
One of the open questions concerning the Hadronic Light-by-Light scattering contribution to the muon (HLBL) is the role of experimental data.
Part of the difficulty of including experimental data in the HLBL is due to the particular framework where the main calculations are done , the large- of QCD . In such limit, one uses the resonance saturation scheme to reproduce the pseudoscalar transition form factor (TFF) that appears in the dominant piece of the HLBL, the pseudscalar-exchange contribution . The main inputs are, then, the pion decay constant and the values of the resonance masses. On top, even though data on the TFF are willing to be included, one still faces the problem on how to link the different kinematic regimes between the experiment for the TFF (low-energy region at time- and space-like, together with intermediate energies at space like) and the kinematics for the pseudoscalar-exchange diagram (whole space-like energy region, from origin of energies up to infinity) .
In this talk I summarized our attempt to provide an answer to that question in a model-independent fashion 
The method proposed can be summarized as follows:
our attempt is, indeed, a method, not a model.
it shall be simple, easy to understand, to apply and reproduce, in contrast to more involved procedures such as dispersion relations.
it may contain approaches (to say, improvable), but not assumptions (not improvable).
it should be systematic, easy to update with new experimental data but also it should provide a systematic error, a pure error from the method itself.
finally, it should be predictive meaning checkable
The method proposed is based on the mathematical theory of Padé approximants (PA). It was pointed out in Ref.  that, in the large- framework, the resonance saturation scheme employed in the HLBL can be understood from the theory of PA to meromorphic functions , where one can compute the desired quantities in a model-independent way. Also, the analytical properties of the TFF indicate that the convergence of the PA is guaranteed at the energies we are working with. Altogether defines the systematics of our method (more iterations should give better approximation) and ascribe an error to that .
To exemplify the advantages of our method, we considered a model from Ref. . As stated before, the inputs for the model can have two different sources: first, a pure theoretical origin based on large- and chiral limits (inputs are resonance masses within the half-width rule  and the meson decay constant in the chiral limit ); second, a reconstruction of the models based on a matching with the TFF low-energy constants , i.e, á la PA [24, 25, 26] minimizing in such a way the model dependence (see  for details). The former yields a final error for the contribution to HLBL to ( from and from the masses). The later, the Padé method, yield a similar provided that the error on the slope ( on curvature) implies an error of () in the pion contribution; the impact of is dramatic since its error implies a error on HLBL. Interestingly enough, the central value is higher, driving non compatible results! The PA method, predictable , can accommodate space- and time-like data. Provides also a rule-of-thumb for estimating the impact of experimental uncertainties, a point never discussed before.
Notice, nevertheless, that the standard procedure [17, 27] to treat the TFF is through a factorization approach, e.g., where is the measured quantity. The impact of such approach is not negligible (see my contribution in ).
In conclusion, we remark the important role of experimental data to determine the dominant pieces of the HLBL (i.e., ). We argue that the way of including such information should be based on PA which provides a systematic error and a simple rule for estimating the impact of experimental uncertainties, both from the space- and the time-like . We notice, finally, that the errors discussed above have been unfortunately ignored in the main reviews (no error for or resonance masses have been properly estimated, neither the possibility to match with experimental low-energy description of the TFF) and that posses a warrant on the reliability of the current error estimates for the HLBL. Similar discussions concerning as well as mixing angles were also addressed.
As discussed during the meeting, to extract resonance poles using PA, see Ref..
a.4 On the positronium contribution to the electron
M. Fael and M. Passera
Albert Einstein Center for Fundamental Physics, Institute for Theoretical Physics, University of Bern, CH-3012 Bern, Switzerland
Istituto Nazionale Fisica Nucleare, Sezione di Padova, I-35131 Padova, Italy
The leading contribution of positronium, the bound state, to the anomalous magnetic moment of the electron has been computed in Ref. . The result of this calculation,
where and is the fine-structure constant, is of the same order of as the perturbative QED five-loop contribution  and comparable with the present experimental uncertainty . As it seems reasonable to expect a reduction of to a part in (or better) in ongoing efforts to improve this measurement, and work is in progress to reduce the error induced in the theoretical prediction for by the uncertainty of [36, 37], a test of the electron -2 at the level of (or below) is a goal that may be achieved not too far in the future. This will bring to play a pivotal role in probing new physics and allow to test whether the long-standing 3–4 discrepancy in the muon -2 also manifests itself in the electron one .
Recently the authors of Ref.  pointed out the presence of the continuum nonperturbative contribution
arising from the region right above the threshold, which corresponds to scattering states with the exchange of Coulomb photons. Comparing Eqs. (1) and (2) they showed that this additional nonperturbative contribution cancels one-half of that of the positronium poles. The question is therefore how to deal with the remaining half: should one add it to the perturbative five-loop QED result of Ref. ? Reference  argued that this remaining term is already contained in the perturbative contribution to computed in Ref.  and, therefore, it should not be added to it. On the other hand, in Ref.  it was claimed that positronium contributes to only through diagrams of or higher. Also, on more general grounds , Ref.  argued that does not exist.
In order to clarify this issue, in Ref.  we used the closed form for the QED vacuum polarization function near the threshold of Refs. [41, 44] to verify that the total (positronium poles plus continuum) nonperturbative contribution to arising from the threshold region is equal to . Then, using the analytic QED vacuum polarization at four-loop of Ref. , we showed explicitly that the perturbative five-loop calculation of of Ref.  does indeed contain the remaining term , in agreement with the arguments of Ref. . We also showed that this term arises from the class I(i) of five-loop diagrams of Ref.  containing only one closed electron loop.
In conclusion, we showed explicitly that there is no additional contribution of QED bound states to beyond perturbation theory.
Acknowledgments The work of M.F. is supported by the Swiss National Science Foundation. M.P. thanks the Department of Physics and Astronomy of the University of Padova for its support. His work was supported in part by the PRIN 2010-11 of the Italian MIUR and by the European Program INVISIBLES (PITN-GA-2011-289442).
a.5 Hadronic light-by-light scattering in the muon : a dispersive approach
Institut für Kernphysik, Technische Universität
ExtreMe Matter Institute EMMI, GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany
The uncertainty in the Standard-Model prediction for the anomalous magnetic moment of the muon is dominated by strong interactions. While hadronic vacuum polarization is intimately related to by means of a dispersion integral, a similarly data-driven approach has only recently been suggested for hadronic light-by-light scattering (HLbL) (see [47, 48] for even higher-order hadronic contributions). Our framework [49, 50] exploits the analytic structure of the HLbL tensor, concentrating on pseudoscalar poles and two-meson intermediate states, which dominate at low energies.
The calculation of the pion transition form factor  can be understood as a generalization of existing analyses for  and , which provide access to the form factor at fixed isoscalar virtualities and , respectively. However, the normalization of the amplitude for a given cannot be predicted, but needs to be inferred from experiment, in case of by means of the Wess–Zumino–Witten anomaly, in case of via the decay width. For general , the normalization is extracted from cross-section data, providing a prediction for without adjusting further parameters. So far, the phenomenological analysis of the singly-virtual form factor has been carried out, including accurate predictions for the slope of the form factor and its analytic continuation into the space-like region.
A crucial step in the derivation of our dispersive formalism  concerns the construction of a suitable basis for the HLbL tensor, in such a way that the coefficient functions are free of kinematic singularities [63, 64, 65]. Moreover, contributions involving double-spectral regions need to be considered separately, so that the sQED pion loop augmented with pion vector form factors (FsQED), as identified on the level of the Mandelstam representation, is evaluated based on Feynman loop integrals. In the talk, also a first numerical evaluation of -wave intermediate states was presented. Despite the double-spectral regions and solely based on -waves the FsQED contribution can be reproduced at the level. Including rescattering in the partial waves in a simplified formalism that involves a Born-term left-hand cut and a finite matching point below the threshold, we find that the sum of and rescattering contributes and, taken together with FsQED, to HLbL scattering in the muon .
a.6 Primary Monte-Carlo generator of the process for the CMD-3 experiment
Budker Institute of Nuclear Physics and Novosibirsk State University,
Russia, 630090, Novosibirsk
Electron-positron collider VEPP-2000  has been operating in Budker Institute of Nuclear Physics since 2010. Center-of-mass energy () range covered by the collider is from threshold of hadron production and up to 2 GeV. Special optics, so called “round beams”, used in the collider construction, allowed to obtain luminosity cms at 1.8 GeV.
The general purpose detector CMD-3 has been described in detail elsewhere . Its tracking system consists of a cylindrical drift chamber (DC)  and double-layer multiwire proportional Z-chamber, both also used for a trigger, and both inside a thin (0.2 X) superconducting solenoid with a field of 1.3 T. The liquid xenon (LXe) barrel calorimeter with 5.4 X thickness has fine electrode structure, providing 1-2 mm spatial resolution , and shares the cryostat vacuum volume with the superconducting solenoid. The barrel CsI crystal calorimeter with thickness of 8.1 X is placed outside the LXe calorimeter, and the end-cap BGO calorimeter with a thickness of 13.4 X is placed inside the solenoid . The luminosity is measured using events of Bhabha scattering at large angles .
Physics program of the CMD-3 experiment includes the study of the multi-hadron production. The cross section measurement of the process in GeV has been already published . §± Preliminary results the study 2 final state has been reported .
The study of intermediate states which lead to 2 final state is essential to correctly describe the angular correlations between the particles and determine the registration efficiency of the process under study. As it was reported at  the intermediate states , and allow satisfactorily describe mass and angular distributions of the 2 production in GeV.
But for higher it is not possible to describe signal, seen in three-pion mass distribution of the 2, by contribution either or , because corresponding cross sections are small enough. We supposed that can be explained by the process with dominant decay of into . The primary Monte-Carlo generator for the process has been created out and implemented into the CMD-3 experiment Monte-Carlo simulation package. Using the generator the signal of the process has been observed in the experimental data for 2 final state.
However, it was not possible to describe 2- 3 and 4 mass distributions as well as angular correlations for 2 final state at GeV by contributions of , and intermediate states. Experimental 4 mass spectra demonstrate presence of narrow state. One of the candidate for this state is with dominant decay into and decays into 2 and . So, the next step in the study of the 2 dynamics will be creation of the primary Monte-Carlo generator of the process .
a.7 Automation of the leading order calculations for hadrons
Institute of Physics, University of Silesia, ul. Uniwersytecka 4, PL-40 007
After some modifications, carlomat [75, 76], a program for automatic computation of the leading order (LO) cross sections of multiparticle reactions, that was originally dedicated mainly to description of the processes of production and decay of heavy particles such as top quarks, the Higgs boson, or electroweak gauge bosons, can be used to obtain predictions for in the framework of effective models. At low energies, the hadronic final states consist mostly of pions, kaons, or nucleons which can be accompanied by one or more photons, or light fermion pairs such as , or . Some effective models which can be useful in this context, including the scalar electrodynamics (sQED) and the interaction with operators of dimension up to 5, were already implemented in version 2 of the program .
The effective Lagrangian of the interaction has the following form :
where the couplings , , , can be complex in general. The electromagnetic (EM) interaction of spin 1/2 nucleons has a similar form:
The form factors and , where , were adopted from PHOKARA , thus making possible Monte Carlo (MC) simulations of processes involving the EM interaction of nucleons.
At low energies, can be treated as point like particles and their EM interaction can be effectively described in the framework of sQED  the interaction vertices examples of which are shown in Fig. 5.
Another step toward better description of at low energies is the inclusion of the Feynman rules of the Resonance Chiral Perturbation Theory (RChPT). The interaction vertices and particle mixing terms of RChPT that can be relevant in this context were provided by Fred Jegerlehner . Some examples of them are shown in Figs. 6 and 7. The implementation of the triple and quartic interaction vertices was more or less straightforward, as it just required writing a few new subroutines for computation of the helicity amplitudes involving the Lorentz tensors that are different from those of the sQED vertices. The couplings , , , , and are currently set either to 1 or . However, implementation of the particle mixing is more challenging, because it must be added at the stage, where the topologies of diagrams which, in carlomat, contain only triple and quartic vertices, are confronted with the Feynman rules. This required substantial changes in the code generating part of the program.
To illustrate how the program works, consider the process . Taking into account the Feynman rules of the standard model and the rules of Figs. 5, 6 and 7, carlomat generates the gauge invariant matrix element, which receives contributions from 903 LO Feynman diagrams, together with a dedicated multichannel phase space integration routine in just a few seconds. A computation of the total cross section, including any number of differential distributions, which is performed as the next step, takes several dozen seconds or several minutes time, dependent on the desired precision of the MC integration.
This project was supported in part with financial resources of the Polish National Science Centre (NCN) under grant decision No. DEC-2011/03/B/ST6/01615.
a.8 MCGPJ for the processes for experiments with CMD-3 detector at the VEPP-2000 collider
G.V. Fedotovich, V.L. Ivanov, D.N. Shemyakin
Budker Institute of Nuclear Physics, SB RAS, Novosibirsk, 630090, Russia
Novosibirsk State University, Novosibirsk, 630090, Russia
The hadronic contribution to of muon, coming from VEPP-2000 energy range, is about . One of the aims of the experiments with CMD-3 detector is to measure the main multihadrons cross sections (MHCS) with systematic uncertainty smaller than ( in ). The systematic accuracy of the cross section for the channel at least better than is required, as it is seen from the “bench” estimation. It is well known that the hadronic contribution to anomalous magnetic moment of muon (AMM) is about : . The aim of the new FNAL and JPARC experiments to measure the of muon is to improve the previous BNL result by a factor of 4 and to achieve the accuracy of - good test of the SM. It is obvious that the systematic uncertainly of radiative corrections (RC) for the MHCS should be smaller than .
Previous experience of the studying of channel at CMD-2 confirms, that final state radiation (FSR) contributes to the cross section (CS) at the level of . Obviously, that for the multihadron channels the contribution of FSR to CS will be smaller than . In the scale of for expected systematic accuracy we can neglect by FSR and consider photon jets radiation in collinear regions only. After many discussions with our experts in Dubna (JINR) we chose the following strategy, which will be described using the channels and as the examples.
To select the clean signal events at first step four charged particles with zero net charge are selected using information from the drift chamber (DC). At the second step after the procedure of separation of kaons and pions (using information ) we calculate the values: and . Two dimensional plot vs for selected events is shown in Fig.1. The events of the signal process are inside the rectangle. The analysis of these events revealed that at least four intermediate states exist: , , , . The simplest model was chosen to describe (more or less) correctly the experimental angular and momentum distributions. In Figures 2-5 the results of simulation vs experiment are presented. To increase the statistics we combine the data collected at several energy points. Red points correspond to experiment, black points - to simulation according to ”realistic” model mentioned above, blue points - to simulation according to the phase space. When the events are selected and the dynamics of their production is defined we are able to calculate the detection efficiency and to determine the visible CS. In order to calculate the RC we calculate the following integral with Structure Functions , which describe photon jets radiation in the collinear regions:
At first iteration in the integral we use born CS, measured by BABAR (if there were not previous measurements of the born CS, then instead of the latter the visible CS should be used). This procedure is repeated for several iterations while RC does not become stable inside a corridor of . The results of such calculation are plotted in Fig.6. The uncertainties of the RC are caused by the uncertainty of the form of the CS.
This final state is almost completely produced by the mechanism. In our analysis we do not reconstruct the from its decay products and do not use the information from calorimeters. This approach allows us to use all the modes of decay and thus enlarges the statistics, but complicates background subtraction. The main sources of background here are and (especially their intermediate mechanism).
Using the information from DC we calculate the parameter :
which represents the total energy of the final particles minus twice beam energy in assumption that the missing particle is the -meson. The distribution of this parameter for the has a peak around zero for any energy point and is used for defining the number of selected signal events . The shape of the distribution of this parameter for background is determined from MC. The results of the fitting of the experimental distributions are summarized in the Fig.7.
Having gotten the and the detection efficiency we calculate the visible CS. The calculation of the RC is performed according the procedure, described above and shown in Fig.8.
In summary, we pointed out that until the precision of the multihadron cross sections can be studied without taking into account the FSR, i.e. considering photon jets radiation in collinear regions only. As the examples of such approach we have considered the calculation of RC for the processes and , which are under study at CMD-3.
a.9 and production at colliders - preliminary results
S. Tracz, H. Czyż, P. Kisza
Institute of Physics, University of Silesia, Katowice, Poland
With the improving luminosity of colliders, the search for a production of , and states become possible. The production of these states goes through two virtual photons. The amplitude describing creation of the state through reaction , going through loop diagram , is proportional to electron mass and thus highly suppressed. All states can be however produced through , , reaction.
Measurements of the cross sections of the reaction will allow to measure the electronic widths () of the and resonances. Combined with measurements of the differential cross section of the reactions they will allow for detailed tests of models describing these charmonium bound states.
Expected range of have been calculated inside two models, quarkonium model and vector dominance model already in . Within the quarkonium model the amplitudes describing coupling of two virtual photons to states depend on binding energy and the derivative of the wave function at the origin. This model predicts also that only some of the allowed amplitudes contribute. For from two allowed amplitudes only one contributes. For from three independent amplitudes one gives contribution. While for from five possible amplitudes only one gives contribution.
For production of states one can concentrate on a selected final state which is easy to be observed experimentally, mainly the decay of into , where subsequently decays into pair of muons. For the the process has to be considered together with radiative return process, which is a non-reducible background (see Figure 16 for diagrams). Furthermore one should take into account the unknown relative phase of the signal and the background amplitudes, which could give an interesting interference pattern.
The binding energies and the derivative of the wave functions at the origin for can be extracted from known  values of and Assuming that the binding energies are different for each state and the derivative of the wave functions are equal, we have performed four fits. First fit was done to the data for and , second to the data for and , third to the data for and and fourth one to all data. The obtained results show that it is impossible to fit simultaneously the data for the states and . Thus the quarkonium model is not able to accommodate these data, even if the discrepancy is not dramatic. In the case of the global fit only width does not fit well. One has to remember however that the model is non-relativistic, while obtained binding energies and hence velocities of quarks are large. Using the fitted parameters we have made predictions of electronic widths. They are not greater than , which is at the limit of BES-III sensitivity.
To examine a possibility of studies of amplitudes at meson factories, we have calculated the cross section of the reaction within the same model. We have assumed integrated luminosities at energy 4.23 GeV(BES-III), (BaBar), (BELLE) and (BELLE-2). In the last three we have used energy 10.56 GeV. For BES-III energy and luminosity the event rates are too small to be observed. The biggest expected rate is for with 20 expected events if the final electron and positron are not tagged. For BaBar, BELLE and BELLE-2 the expected event rates with no electron-positron tagging are 1700,3300,160000 for , 19000,36000,1800000 for and 26000,50000,2500000 for . For BELLE-2, even if one tags both electron and positron within the angular range , the expected event rates are big enough to be observed and are equal to 7500, 400000,100000 for respectively. The stricking difference between and seen here comes from the fact that for the contribution from real photons is equal to zero and thus for small photons virtualities the amplitude is small.
The modes of production of with the subsequent decay described above have been implemented in PHOKHARA Monte Carlo generator and the production and decay mode was added in EKHARA Monte Carlo generator. Preliminary results shown above indicate that amplitudes can be studied in existing or near future experiments. More detailed analysis will be presented in .
Work supported in part by the Polish National Science Centre, grant number DEC-2012/07/B/ST2/03867 and German Research Foundation DFG under Contract No. Collaborative Research Center CRC-1044. Szymon Tracz and Patrycja Kisza are supported by the Forszt project co-financed by EU from the European Social Fund.
a.10 Nucleon form factors in PHOKHARA
H. Czyż, S. Tracz
Institute of Physics, University of Silesia, PL-40007 Katowice, Poland
The electromagnetic nucleon form factors were studied from the conception of particle physics [87, 88, 89, 90]. and yet a lot has to be done to built a model which meets requirements in the era of precision hadronic physics. The not expected developments in this field show that we still have much to learn. The measurements of the ratio of the electric and magnetic form factors of the proton with two different methods were giving different results . The two-photon exchange radiative corrections explaining to large extent this difference (see  for review) turned out to be unexpectedly large. Moreover their modelling goes beyond the nucleon form factor modelling adding to the complexity of the problem. In addition, usually the models are built separately in the space-like (see  for review) or time-like regions (see  for review). As one expects that each of the form factors is a unique analytic function valid in both space-like and time-like regions this attitude has to be changed and a model describing both regions has to be constructed. A step towards such a model was done in . The model describes well data from both space-like and time-like regions. The form factors are normalised properly at zero invariant mass and, by construction, have correct behaviour  at large invariant masses. Yet the model is far from being satisfactory and further studies of photon-nucleon interactions are necessary. It is also clear that the progress can be achieved only through close collaboration of experimental and theory groups. Careful studies of charge and/or forward-backward asymmetries in , processes together with angular distributions in and processes should allow for disentangling of the two-photon exchange contributions from the form factors.
Model testing is simplified if it is implemented into a Monte Carlo event generator. Such a generator serves also for other purposes like calculations of acceptance and/or efficiency corrections. For the radiative return (called also ISR) method such a tool was developed some time ago . Nucleon final states were implemented in it already in  and the nucleon form factors were updated in , where also the modelling of the final state radiative corrections was addressed. Recently  also a possibility of generation of the process hadrons, useful for scan experiments, was added.
Work supported in part by the Polish National Science Centre, grant number DEC-2012/07/B/ST2/03867 and German Research Foundation DFG under Contract No. Collaborative Research Center CRC-1044. Szymon Tracz is supported by the Forszt project co-financed by EU from the European Social Fund.
a.11 Current status of Monte Carlo generator Tauola
Institute of Nuclear Physics PAN, ul. Radzikowskiego 152, 31-342 Krakow, Poland
Kharkov Institute of Physics and Technology, Akademicheskaya,1, 61108 Kharkov, Ukraine
The Monte Carlo generator TAUOLA , which is used to simulate tau lepton decays, is a computing project with a rather long history. It started in the years 80’s and is under development still now. The main problem in the theoretical description of the hadronic decay modes of the tau lepton is a lack of a theory coming from the first principle in the energy region from the threshold till the tau mass. The hadronic currents in the first version of TAUOLA as well as the subsequent internal versions of the code used by both Aleph and Cleo collaborations were based on the Vector Meson Dominance approach. However, that approach is able to reproduce only the leading-order results of Chiral Perturbation Theory (ChPT). As was shown in the case of the mode this approach spoils the Wess-Zumino anomaly normalization, that appears at in ChPT . An alternative approach is to include the lightest resonances as active degrees of freedom in the theory. This can be done by adding resonance fields to the ChPT Lagrangian, without any dynamical assumption, leading to Resonance Chiral Lagrangian approach [102, 103]. The RChL approach succeeds in reproducing low energy results, predicted by ChPT, at least at the next-to-leading order, and also complies with QCD high energy constraints.
Up to now the RChL currents for the main two-meson (final states with two pion, pion-kaon, two kaons) and three-pseudoscalar (three pion, two kaon-one pion) decay modes have been installed into TAUOLA. This set covers more than of the tau lepton hadronic decay width. The implementation of the currents, the related technical tests as well as the necessary theoretical concepts are documented in .
To get numerical values of the model parameters one has to fit the theoretically predicted spectra to those measured in experiments. We started with the decay. The first fit to the preliminary BaBar data  allowed to extract numerical values of the model parameters and demonstrated satisfactory agreement with the three pion invariant mass spectrum and a mismatch in the low energy tail of the two pion invariant mass distributions . Improvement has been achieved by adding the resonance contribution in the hadronic current . To fit the data we used the MINUIT package through the ROOT framework and the fit result is presented in Fig. 1 of . For the numerical values of the model parameters see Table 1 in . The goodness of the fit is quantified by . To compare with the previous result  we have estimated the value using the combined statistical and systematic uncertainties and we find that is eight times better than in .
The following tests have been done to validate our fitting procedure: asserting statistical errors and correlation coefficients between the model parameters, verifying that the obtained result is a global minimum and performing studies of the systematical errors. A detailed description as well as the results of the tests are presented in .
The fitting procedure for the three pion mode has been generalized to an arbitrary three meson final state. Thus the new fitting framework allows to perform fits for arbitrary tau decay mode, using either Fortran or C++ code. To test it, we have first reproduced the results for the three pion decay mode described above and then fitted the currents to the BaBar preliminary data . For the moment we have used some simplifications for the for current which will be removed in the final version of the fitting strategy. In addition, fits for decay mode has been added allowing to fit the RChL two pion form factor to the Belle parameterization. Next step will be to study the stability of the generalized fitting procedure and to include the experimental errors.
Another important task concerned the C++ interface for the decay channels of TAUOLA generator. The interface allows to add, substitute or modify TAUOLA decay modes using either C++ or FORTRAN code. The interface is ready to be used in Tauola++ project (http://tauolapp.web.cern.ch/tauolapp/) as well as in the FORTRAN environments, where TAUOLA is still being used. Validation of all BaBar currents introduced along with the C++ environment has been finished and extensive testing of the environment has been performed.
This research was supported in part by Foundation of Polish Science grant POMOST/2013-7/12, that is co-financed from European Union, Regional Development Fund and from funds of Polish National Science Centre under decisions DEC-2011/03/B/ST2/00107.
Appendix B List of participants
Antonio Anastasi, University of Messina and LNF, firstname.lastname@example.org
Marcin Berlowski, National Center for Nuclear Research, Marcin.Berlowski@fuw.edu.pl
Bo Cao, University of Uppsala, email@example.com
Carlo M. Carloni Calame, University of Pavia, firstname.lastname@example.org
Henryk Czyż, University of Silesia, email@example.com
Achim Denig, Universtät Mainz, firstname.lastname@example.org
Gennadiy V. Fedotovich, BINP and NSU, G.V.Fedotovich@inp.nsk.su
Martin Hoferichter, University of Bern, email@example.com
Karol Kołodziej, University of Silesia, firstname.lastname@example.org
Andrzej Kupsc, University of Uppsala, Andrzej.Kupsc@physics.uu.se
Peter Lukin, BINP and NSU, P.A.Lukin@inp.nsk.su
Pere Masjuan, Universtät Mainz, email@example.com
Dario Moricciani, INFN - Roma Tor Vergata, firstname.lastname@example.org
Massimo Passera, INFN Padova, email@example.com
Olga Shekhovtsova, IFJ PAN, Olga.Shekhovtsova@lnf.infn.it
Szymon Tracz, University of Silesia, firstname.lastname@example.org
Graziano Venanzoni, LNF, Graziano.Venanzoni@lnf.infn.it
Ping Wang, Institute of High-energy Physics, Beijing, email@example.com
- S. Actis et al. [Working Group on Radiative Corrections and Monte Carlo Generators for Low Energies Collaboration], Eur. Phys. J. C 66 (2010) 585 [arXiv:0912.0749 [hep-ph]].
- P. Masjuan, G. Venanzoni, H. Czyż, A. Denig, M. Vanderhaeghen, G. Venanzoni, A. Denig and S. Eidelman et al., arXiv:1306.2045 [hep-ph].
- H. Czyż, S. Eidelman, G. V. Fedotovich, A. Korobov, S. E. Mller, A. Nyffeler, P. Roig and O. Shekhovtsova et al., arXiv:1312.0454 [hep-ph].
- J. J. van der Bij, H. Czyż, S. Eidelman, G. Fedotovich, T. Ferber, V. Ivanov, A. Korobov and Z. Liu et al., arXiv:1406.4639 [hep-ph].
C. M. Carloni Calame et al., Nucl. Phys. B 584 (2000) 459
C. M. Carloni Calame et al., Nucl. Phys. Proc. Suppl. 131 (2004) 48 [hep-ph/0312014]
G. Balossini et al., Nucl. Phys. B 758 (2006) 227 [hep-ph/0607181]
G. Balossini et al., Phys. Lett. B 663 (2008) 209 [arXiv:0801.3360 [hep-ph]]
- S. Jadach, W. Placzek, B. F. L. Ward, Phys. Lett. B 390 (1997) 298 [hep-ph/9608412]
- A. B. Arbuzov et al., Eur. Phys. J. C 46 (2006) 689 [hep-ph/0504233]
- E. Drago, G. Venanzoni, INFN-AE-97-48 (1997)
F. A. Berends, R. Kleiss,
Nucl. Phys. B 228 (1983) 537
F. A. Berends, R. Kleiss, Nucl. Phys. B 186, 22 (1981)
- S. Actis et al., Eur. Phys. J. C 66 (2010) 585 [arXiv:0912.0749 [hep-ph]].
- C. M. Carloni Calame et al., JHEP 1107 (2011) 126 [arXiv:1106.3178 [hep-ph]].
- B.I.Khazin et al., Nucl.Phys.B, Proc. Suppl. 181-182, 376 (2008).
V.V. Danilov et al., Proceedings of EPAC96, Barcelona,
I.A. Koop, Nucl. Phys. B (Proc. Suppl.) 181-182, 371 (2008).
- G.Grawford et al., NIM A345, 429 (1994).
- A.Arbuzov et al., EPJ C46, 689 (2006).
- S.Actis et al., Eur. Phys. J. C66, 585 (2010).
- F. Jegerlehner and A. Nyffeler, Phys. Rept. 477, 1 (2009) [arXiv:0902.3360 [hep-ph]].
- G. ’t Hooft, Nucl. Phys. B 72, 461 (1974).
- P. Masjuan and P. Sanchez-Puertas, In preparation.
- P. Masjuan, arXiv:1411.6397 [hep-ph].
- See M. Hoferichter in this proceedings.
- P. Masjuan, Phys. Rev. D 86, 094021 (2012) [arXiv:1206.2549 [hep-ph]]; P. Masjuan and M. Vanderhaeghen, arXiv:1212.0357 [hep-ph].
- R. Escribano, P. Masjuan and P. Sanchez-Puertas, Phys. Rev. D 89 (2014) 3, 034014 [arXiv:1307.2061 [hep-ph]].
- P. Masjuan and S. Peris, JHEP 0705, 040 (2007) [arXiv:0704.1247 [hep-ph]]; Phys. Lett. B 663, 61 (2008) [arXiv:0801.3558 [hep-ph]].
- G. A. Baker and P. Graves-Morris, Encyclopedia of Mathematics and its Applications, Cambridge Univ. Press, 1996; P. Masjuan Queralt, arXiv:1005.5683 [hep-ph].
- P. Masjuan, et al., Phys. Lett. B 668 (2008) 14 [arXiv:0805.3291 [hep-ph]]; P. Masjuan and S. Peris, Phys. Lett. B 686, 307 (2010) [arXiv:0903.0294 [hep-ph]].
- M. Knecht and A. Nyffeler, Phys. Rev. D 65, 073034 (2002) [hep-ph/0111058].
- P. Masjuan, E. Ruiz Arriola and W. Broniowski, Phys. Rev. D 85, 094006 (2012) [arXiv:1203.4782 [hep-ph]]; Phys. Rev. D 87, 014005 (2013) [arXiv:1210.0760 [hep-ph]]; Phys. Rev. D 87 (2013) 11, 118502 [arXiv:1305.3493 [hep-ph]].
- G. Ecker, et al., Eur. Phys. J. C 74 (2014) 2748 [arXiv:1310.8452 [hep-ph]].
- P. Aguar-Bartolome et al. [A2 Collaboration], Phys. Rev. C 89 (2014) 4, 044608 [arXiv:1309.5648 [hep-ex]].
- P. Adlarson et al., arXiv:1412.5451 [nucl-ex].
- P. Masjuan et al., Eur. Phys. J. C 73 (2013) 2594 [arXiv:1306.6308 [hep-ph]]; Phys. Rev. D 90 (2014) 9, 097901 [arXiv:1410.2397 [hep-ph]].
- G. Mishima, arXiv:1311.7109 [hep-ph].
- T. Aoyama, M. Hayakawa, T. Kinoshita and M. Nio, Phys. Rev. Lett. 109 (2012) 111807 [arXiv:1205.5368 [hep-ph]].
- D. Hanneke, S. Fogwell and G. Gabrielse, Phys. Rev. Lett. 100 (2008) 120801.
- R. Bouchendira et al., Phys. Rev. Lett. 106 (2011) 080801.
- F. Terranova, G. M. Tino, Phys. Rev. A 89 (2014) 052118 [arXiv:1312.2346 [hep-ex]].
- G. F. Giudice, P. Paradisi and M. Passera, JHEP 1211 (2012) 113 [arXiv:1208.6583 [hep-ph]].
- K. Melnikov, A. Vainshtein and M. Voloshin, Phys. Rev. D 90 (2014) 017301 [arXiv:1402.5690 [hep-ph]].
- M. Hayakawa, arXiv:1403.0416 [hep-ph].
- M. A. Braun, Zh. Eksp. Teor. Fiz. 54 (1968) 1220 [Sov. Phys. JETP. 27 (1968) 652].
- M. I. Eides, Phys. Rev. D 90 (2014) 057301 [arXiv:1402.5860 [hep-ph]].
- M. Fael and M. Passera, Phys. Rev. D 90 (2014) 056004 [arXiv:1402.1575 [hep-ph]].
- R. Barbieri, P. Christillin and E. Remiddi, Phys. Rev. A 8 (1973) 5, 2266.
- P. A. Baikov, A. Maier and P. Marquard, Nucl. Phys. B 877 (2013) 647 [arXiv:1307.6105 [hep-ph]].
- T. Aoyama, M. Hayakawa, T. Kinoshita and M. Nio, Phys. Rev. D 83 (2011) 053003 [arXiv:1012.5569 [hep-ph]].
- A. Kurz, T. Liu, P. Marquard and M. Steinhauser, Phys. Lett. B 734 (2014) 144 [arXiv:1403.6400 [hep-ph]].
- G. Colangelo, M. Hoferichter, A. Nyffeler, M. Passera and P. Stoffer, Phys. Lett. B 735 (2014) 90 [arXiv:1403.7512 [hep-ph]].
- G. Colangelo, M. Hoferichter, M. Procura and P. Stoffer, JHEP 1409 (2014) 091 [arXiv:1402.7081 [hep-ph]].
- G. Colangelo, M. Hoferichter, B. Kubis, M. Procura and P. Stoffer, Phys. Lett. B 738 (2014) 6 [arXiv:1408.2517 [hep-ph]].
- V. Pauk and M. Vanderhaeghen, arXiv:1409.0819 [hep-ph].
- T. Blum, S. Chowdhury, M. Hayakawa and T. Izubuchi, arXiv:1407.2923 [hep-lat].
- M. Hoferichter, B. Kubis, S. Leupold, F. Niecknig and S. P. Schneider, Eur. Phys. J. C 74 (2014) 3180 [arXiv:1410.4691 [hep-ph]].
- R. García-Martín and B. Moussallam, Eur. Phys. J. C 70 (2010) 155 [arXiv:1006.5373 [hep-ph]].
- M. Hoferichter, D. R. Phillips and C. Schat, Eur. Phys. J. C 71 (2011) 1743 [arXiv:1106.4147 [hep-ph]].
- B. Moussallam, Eur. Phys. J. C 73 (2013) 2539 [arXiv:1305.3143 [hep-ph]].
- M. Hoferichter, G. Colangelo, M. Procura and P. Stoffer, arXiv:1309.6877 [hep-ph].
- L. Y. Dai and M. R. Pennington, Phys. Rev. D 90 (2014) 036004 [arXiv:1404.7524 [hep-ph]].
- C. Hanhart, A. Kupść, U.-G. Meißner, F. Stollenwerk and A. Wirzba, Eur. Phys. J. C 73 (2013) 2668 [arXiv:1307.5654 [hep-ph]].
- M. Hoferichter, B. Kubis and D. Sakkas, Phys. Rev. D 86 (2012) 116009 [arXiv:1210.6793 [hep-ph]].
- F. Niecknig, B. Kubis and S. P. Schneider, Eur. Phys. J. C 72 (2012) 2014 [arXiv:1203.2501 [hep-ph]].
- S. P. Schneider, B. Kubis and F. Niecknig, Phys. Rev. D 86 (2012) 054013 [arXiv:1206.3098 [hep-ph]].
- W. A. Bardeen and W. K. Tung, Phys. Rev. 173 (1968) 1423 [Erratum-ibid. D 4 (1971) 3229].
- R. Tarrach, Nuovo Cim. A 28 (1975) 409.
- R. A. Leo, A. Minguzzi and G. Soliani, Nuovo Cim. A 30 (1975) 270.
V.V. Danilov et al., Proceedings of EPAC96, Barcelona,
I.A. Koop, Nucl. Phys. B (Proc. Suppl.) 181-182, 371 (2008).
- B.I. Khazin, Nucl. Phys. B (Proc. Suppl.) 181-182, 376 (2008).
- F. Grancagnolo et al., Nucl. Instr. and Meth. A 623, 114 (2010).
- A.V. Anisyonkov et al., Nucl. Instr. and Meth. A 598, 266 (2009)
- D. Epifanov (CMD-3 Collaboration), J. Phys. Conf. Ser. 293, 012009 (2011).
- R.R. Akhmetshin et al., Nucl. Phys. B (Proc. Suppl.) 225-227, 69 (2012).
- R.R. Akhmetshin et al., Phys.Lett. B723 82 (2013).
- P.A. Lukin et al., EPJ Web Conf. 81 02010 (2014).
- J.J. van der Bij et al., arXiv:1406.4639 [hep-ph].
- K. Kołodziej, Comput. Phys. Commun. 180 (2009) 1671, [arXiv:0903.3334].
- K. Kołodziej, Comput. Phys. Commun. 185 (2014) 323, [arXiv:1305.5096].
- G.L. Kane, G.A. Ladinsky, C.-P. Yuan, Phys. Rev. D45 (1992) 124.
- S. Tracz, H. Czyż, Acta Phys. Pol. B44 (2013) 2281.
- F. Jegerlehner, The Anomalous Magnetic Moment of the Muon, Springer, Nov. 2007, ISBN 9783540726333.
- M. Benayoun, P. David, L. DelBuono, F. Jegerlehner, Eur. Phys. J. C73 (2013) 2453.
- A.B. Arbuzov (Dubna, JINR), G.V. Fedotovich, F.V. Ignatov (Novosibirsk, IYF) et al. Eur.Phys.J. C46 (2006) 689-703
- D.N. Shemyakin (Novosibirsk, IYF & Novosibirsk State U.), R.R Akhmetshin, A. Anisenkov (Novosibirsk, IYF) et al. PoS Hadron2013 (2013) 138.
- E.A. Kuraev, V.S. Fadin. Yad.Fiz. 41 (1985) 733-742.
- J. H. Kuhn, J. Kaplan and E. G. O. Safiani, Nucl. Phys. B 157 (1979) 125.
- K. A. Olive et al. [Particle Data Group Collaboration], Chin. Phys. C 38 (2014) 090001.
- H. Czyz, J. Kuhn, Sz. Tracz, P. Kisza (in preparation)
- I. Estermann, R. Frisch, and O. Stern Nature 132 (1933) 169.
- M. N. Rosenbluth, Phys. Rev. 79 (1950) 615.
- R. W. Mcallister and R. Hofstadter, Phys. Rev. 102 (1956) 851.
- R. Hofstadter, Rev. Mod. Phys. 28 (1956) 214.
- J. Arrington, Phys. Rev. C 68 (2003) 034325 [nucl-ex/0305009].
- C. E. Carlson and M. Vanderhaeghen, Ann. Rev. Nucl. Part. Sci. 57 (2007) 171 [hep-ph/0701272 [HEP-PH]].
- C. F. Perdrisat, V. Punjabi and M. Vanderhaeghen, Prog. Part. Nucl. Phys. 59 (2007) 694 [hep-ph/0612014].
- A. Denig and G. Salme, Prog. Part. Nucl. Phys. 68 (2013) 113 [arXiv:1210.4689 [hep-ex]].
- H. Czyz, J. H. Kuhn and S. Tracz, arXiv:1407.7995 [hep-ph].
- G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22 (1980) 2157.
- G. Rodrigo, H. Czyz, J. H. Kuhn and M. Szopa, Eur. Phys. J. C 24 (2002) 71 [hep-ph/0112184].
- H. Czyz, J. H. Kuhn, E. Nowak and G. Rodrigo, Eur. Phys. J. C 35 (2004) 527 [hep-ph/0403062].
- H. Czyz, M. Gunia and J. H. Kuhn, JHEP 1308 (2013) 110 [arXiv:1306.1985, arXiv:1306.1985 [hep-ph]].
- S. Jadach, Z. Was, R. Decker and J. H. Kuhn, Comput. Phys. Commun. 76 (1993) 361.
- T. E. Coan et al. [CLEO Collaboration], Phys. Rev. Lett. 92 (2004) 232001 [hep-ex/0401005].
- G. Ecker, J. Gasser, A. Pich and E. de Rafael, Nucl. Phys. B 321 (1989) 311.
- G. Ecker, J. Gasser, H. Leutwyler, A. Pich and E. de Rafael, Phys. Lett. B 223 (1989) 425.
- O. Shekhovtsova, T. Przedzinski, P. Roig and Z. Was, Phys. Rev. D 86 (2012) 113008 [arXiv:1203.3955 [hep-ph]].
- I. M. Nugent, T. Przedzinski, P. Roig, O. Shekhovtsova and Z. Was, Phys. Rev. D 88 (2013) 9, 093012 [arXiv:1310.1053 [hep-ph]].
- O. Shekhovtsova, I. M. Nugent, T. Przedzinski, P. Roig and Z. Was, arXiv:1301.1964 [hep-ph].
- I. M. Nugent, arXiv:1301.7105 [hep-ex].