TAUOLA, TAUOLA universal interface PHOTOS and MC-TESTER: Status Report
The status of Monte Carlo programs for the simulation of -lepton production and decay in high-energy accelerator experiments is reviewed. No basic changes in the organization of the programs were necessary since the previous -lepton conference in 2006. Minor in size, but practical, extensions for the use of the programs in high precision mixed languages enviromens are being introduced: (i) For the TAUOLA Monte Carlo generator of -lepton decays, automated and simultaneous use of many versions of form-factors for the calculation of optional weights for fits was discussed. A pilot example of its use was presented. (ii) New tests for PHOTOS Monte Carlo for QED bremsstrahlung in decays were shown. (iii) Prototype version of the TAUOLA universal interface based on HepMC (the C++ event record) was mentioned. Its tests with the help of MC-TESTER were discussed.
Presented results illustrate the status of the projects performed in collaboration with Vladimir Cherepanov, Nadia Davidson, Piotr Golonka, Gizo Nanava, Tomasz Przedziński Elżbieta Richter-W¸as and others.
Presented at International workshop on Tau Lepton Physics, TAU08 Novosibirsk, Russia September, 2008
The TAUOLA package [1, 2, 3, 4] for the simulation of -lepton decays and PHOTOS [5, 6] for the simulation of radiative corrections in decays, are computing projects with a rather long history. Written and maintained by well-defined authors, they nonetheless migrated into a wide range of applications where they became ingredients of complicated simulation chains. As a consequence, a large number of different versions are presently in use. From the algorithmic point of view, they often differ only in a few small details, but incorporate many specific results from distinct -lepton measurements or phenomenological projects. Such versions were mainly maintained (and will remain so) by the experiments taking precision data on leptons. On the other hand, many new applications were developed recently, often requiring a program interface to other packages (e.g. generating events for LHC, LC, Belle or BaBar physics processes). The programs structure, prepared for the convenience of users, was presented during previous conferences, and we will not repeat it here.
This time, let us concentrate on physics oriented results and special techniques being developped and used for fits, over the years. They may be worth some attention now, when high precision, high statistic data are being analyzed. FORTRAN to C++ migration is not as simple as one may think. On one hand it is generally believed that large projects are easier to maintain in more modern programming frameworks. On the other hand, FORTRAN software can not be abandoned. It usually carries a lot of physics expertise. It is also impractical to replace the software environment of mature experiments with something new. Fortunately, co-existence of the two languages is not a problem, at least not from the software point of view.
Our presentation is organized as follows. Section 2 is devoted to the discussion of optional weights in TAUOLA and their use for fits to the data already at the level of comparison with raw data. In section 3 we present some new results for the simulation with PHOTOS of radiative corrections in decays. Section 4 is devoted to MC-TESTER; the program which can be used for semi-automatic comparisons of simulation samples originating from different programs. The prototype, C++ based, TAUOLA universal interface is mentioned in that Section too.
Because of the limited space of the contribution, and sizable amount of other physically interesting results, some of them will be excluded from conference proceedings. Let us hope that they will find place in future works, possibly with collaborators mentioned in the Abstract. For these works, the present paper may serve as an announcement.
2 Optional weights in Tauola Monte Carlo
Physics of lepton decays requires sophisticated strategies for the confrontation of phenomenological models with experimental data. On one side high statistics experimental samples are collected, and the obtained precision is high, on the other hand, there is a significant cross-contamination between distinct decay channels. Starting from a certain precision level all channels need to be analyzed simultaneously. Change of parametrization for one channel contributiong to the background may be important for the fit of another one. This situation leads to a complex configuatin where a multitude of parameters needs to be simultaneously confronted with a multitude of observables. One has to keep in mind that the models used to obtain distributionbs in the fits may require refinements or even substantial rebuilds. The isospin relation to distribution measured in collistions at low energies offer additional control, but may be the source of even more complex systematic error considerations. Effects of QED bremsstrahlung need to be taken ito account as well. It may affect selection criteria and background contaminations in quite complex and unexpected ways.
From the statistical point of view it is best to resolve such system in one
automated step using the method
such as [7, 8] for example.
This can be of course very dangerous from
the point of view of systematic error control. We will return to this point
later. At this moment let us discuss necessary adaptations for TAUOLA
for such methods to be used. In fact the principle of the changes is quite
simple. The method is in use since long time and in different situations,
see eg. .
It is enough, to calculate, for each generated event (for each
present in it decay of and separately) alternative weights; the ratios
of matrix element squared obtained with new currents provided by the user,
and the one actually used in generation are then the vector of weights.
If the currents
can be represented as a linear combination of elementary ones, then
the weight for arbitary values
of coefficients is a bilinear form in the relative coefficients
The first example of the use of the method, is demonstrated in fig. 1. At this step of the work no detector effects were included. Experimental data were provided as a one-dimensional histogram. Simulation of production process was not performed. One can see, that the theoretical model provides double peak which is absent in the data. This is because of an interplay between intermediate state resonance parametrization and phase space jacobian. The same parametrization was used for 3-scalar and 2-scalar systems. This points to the necessity of further investigations of unitarity constraints on the running width or better control of the data.
Options how to handle, in experimental software, information for weight calculation is under study now. In principle, information stored in the event record should be enough. This may not be a perfect match in all cases. For example, effects of the boost from laboratory frame to rest-frame of lepton and back, may result in rounding errors. In some cases the kinematical configuration may be changed because of QED bremsstrahlung generated by PHOTOS. The order of final state particles stored in the event record may not be the same as used internally in TAUOLA at the event construction step. This of course can be cross checked on the basis of inspection of a smaller sample and without detector simulation part being used. Such tests are nonetheless necessary, but may need less time than would be the case if collaboration software would be rebuild. From that point of view, the easiest would be to calculate optional weights from information already available in experiments production tapes, and only for those events where bremsstrahlung effects are not generated at all. This can be correted at the next step of form-factor iteration. Related slow down of the convergence should not be large. Possibly two steps will be enough, the second one for the verification. It will be necessary for control of spin correlations between the two decaying -s as well. It is technically much simpler not to include changes into spin correlations. In such a case the vector of weights would depend on kinematical configurations of , and of the production process simultaneously. The optimal choice require further interaction with collaborations. Technical side of the activity will depend on the way the user want to calculate the hadronic current. It can be easily calculated as external generalized function communicating with TAUOLA at the level of linker or through the fifo pipes. They can be provided as FORTRAN or eg. C++ code and have well defined and easy to control input/output for modifications introduced into TAUOLA.
Int the past, and for a good reason, significant effort was devoted to strategies where the methotds as indicated here can be avoided. For example in  a system of projections was proposed for the 3-scalar decays. If such a method can be used it is possible to obtain from the data just one of the contribution to the current, which is then parametrized by a single scalar function of two invariants only. This is of course on the expense of statistical strength on the measurement and abstraction of experimental sytematic error. Such methods, easy to apply for final states of two or three scalars, are less practical for higher multiplicities. Nonetheless they can be of great help to control the analysis and to verify which aspects of the models require improvements.
3 Photos and NLO effects in and decays
The changes introduced over the last two years into PHOTOS Monte Carlo program itself were rather small. On the other hand, the complexity of the work on its theoretical foundation was large and matches neither size nor the main purpose of the present talk (some of the related topics are revieved in ). Previous tests of two body decays of the into a pair of charged leptons  and scalar into a pair of scalars  was recently supplemented  with the study of . The study of the process for is on-going . In all of these cases universal kernel of PHOTOS was replaced with the one matching exact first order matrix element. In this way terms for NLO/NLL level are implemented. As usual, algorithm covers the full multiphoton phase-space and is exact in the infrared region of the phase space. This is rather unusual for the NLL compatible algorithms. One should not forget that PHOTOS generates weight one events and does not require any factorization scales. There is full overlap of the phase space where hard matrix element is used with the one for iterative photon emissions. All interference effects are implemented with the help of internal weights.
The results of all tests of PHOTOS with a NLO kernel are at sub-permille level. No differences with benchmarks can be observed, even for samples of events. For the comparisons when simpler physics assumptions are used, differences were sub-permille too, if compared to total rates (or were matching precision of programs used for tests).
This is very encouraging, and points to the possible extension of the approach outside of the QED (scalar QED) only. In particular to the domain of QCD or if phenomenological lagrangians for interactions of photons need to be applied. For that work to be completed spin amplitudes need to be further studied. Let us point here to the paper .
The discussed here refinements only indirectly affect practical side of simulations for physics. Changes in the kernels necessary for NLO may remain as options for tests only. They are available from the PHOTOS web page  but are not recomended for wider use. An example of such new tests is given in fig. 2b.
4 Mc-Tester and C++ TAUOLA universal interface
In the development of packages such as TAUOLA or PHOTOS, question of tests and appriopriate relations to users applications are essential for their usefulness. In fact, users applications may be much larger in size and human effort than the programs discussed here. Good example of such ‘user applications’ are complete environments to fit experimental data and control detector uncertainties at the same time. As a consequence our programs need to work well with FORTRAN where HEPEVT event record is used and C++ where HepMC  is used instead. For the sake of automatization of tests MC-TESTER was prepared . Recently, new functionalities were introduced into the package . In particular, it works now with the HepMC event record, the standard of C++ programs. The complete set-up for benchmarking the interfaces, such as interface between -lepton production and decay, including QED bremsstrahlung effects is prepared. The example is chosen to illustrate the new options introduced into the program and novel ways of its use. From the technical perspective,  documents program updates and supplements previous documentation .
As in the past, test consists of two steps. Distinct Monte Carlo programs are run separately; events are searched for the decays of the chosen object, and information is stored by MC-TESTER. Then, at the analysis step, information from a pair of such runs may be compared and represented in the form of tables and plots.
One may have the impression that the modifications introduced into the release are minor and consist of simple improvements in the graphical representation of the output and purely technical reorganization thanks to the use of C++ lists.
To some degree this is true, but on the other hand the changes were introduced because of pressure from applications. Let us show how program modifications can be used for non-trivial practical applications.
It is quite common that information stored in the event record is too large. For example individual soft photons which remain undetectable are present. Not only they do not influence the detector response at all, but they exhibit technical aspects of eg. infrared regulators of QED bremsstrahlung. In response, MC-TESTER should ignore those photons while analyzing decays (or group together with other particles). Otherwise, comparisons of different Monte Carlo programs would be dominated by the technical aspect of the implementation of the infrared regulator; MC-TESTER operation need to be adopted for this, see eg. [20, 12].
Another example where the event tree may need to be simplified for validation is
if spin correlations are appropriately introduced into various production
processes. Let us use as an example
Final state activities will lead to final states of decays where the -pair is accompanied by bremsstrahlung photons or soft hadrons. Such distribution are usually strongly peaked. In fig. 2b we provide an example of such a distribution but for the comparison of PHOTOS with exact marix element for the decay .
With the above brief examples we have demonstrated that MC-TESTER may be useful for tests of libraries of particles decays, as well as for tests of their interfaces. The work on HepMC based TAUOLA universal interface is on-going. A preliminary version is already available , but from physics point of view it is still less complete than the FORTRAN version .
5 Summary and future possibilities
The status of the computer programs for the decay of leptons TAUOLA and associated projects TAUOLA universal interface and MC-TESTER was reviewed. The high-precision version of PHOTOS for radiative corrections was presented too.
New results for PHOTOS and decay of were mentioned. For this channel complete next-to-leading order effects can now be simulated. However, for most of the applications these effects are not necessary, leaving the standard modular version of PHOTOS sufficient. This exercise is important not only for photonic bremsstrahlung. It helps to understand better questions related with matching hard emission matrix elements with parton showers without necessity to introduce any boundaries within phase space. The important result of the above work is that the path to include electromagnetic form-factors of the particles participating in decay is now open for future fits to the data. These form-factor effects may be significantly larger and physically more justified than complete next-to-leading order effects of scalar QED as in meson decays.
The presentation of the TAUOLA general-purpose interface was brief. The C++ version exist, but is still physically less refined than the FORTRAN version. Examples for its use can be found on the continuously modified web page . The new version of MC-TESTER  is already public now. It not only works with HepMC  of C++ but enables user defined tests.
Distinct versions of the TAUOLA library for lepton decay, and of PHOTOS for radiative corrections in decays, are now in use. The principles how to use the distribution package are presented in refs. [4, 23], on the other hand, development of new and recovery of codes for old strategies of fits is becoming available. This includes strategy to use weighted events, and the projection operators  for the three scalars final states. Let us hope that the present talk will provide some stimulation for the larger work in the near future.
Discussions with members of the Belle and BaBar collaborations are also acknowledged. Exchange of e-mails and direct discussions with S. Banerjee, S. Eidelman, H. Hayashii, K. Inami, A. Korchin, J. H. Kühn and O. Shekhovtsova was a valuable input to present and future steps in project development.
- We have applied such method in one of our other programs  for the studies of anomalous couplings of and LEP data.
- Spin correlations in decays of into lepton(s) are nearly identical.
- S. Jadach, J. H. Kuhn, and Z. Was, Comput. Phys. Commun. 64 (1990) 275.
- M. Jezabek, Z. Was, S. Jadach, and J. H. Kuhn, Comput. Phys. Commun. 70 (1992) 69.
- S. Jadach, Z. Was, R. Decker, and J. H. Kuhn, Comput. Phys. Commun. 76 (1993) 361–380.
- P. Golonka et al., Comput. Phys. Commun. 174 (2006) 818–835, hep-ph/0312240.
- E. Barberio, B. van Eijk, and Z. Was, Comput. Phys. Commun. 66 (1991) 115.
- E. Barberio and Z. Was, Comput. Phys. Commun. 79 (1994) 291–308.
- A. Hoecker, P. Speckmayer, J. Stelzer, F. Tegenfeldt, H. Voss, K. Voss, et al., http://tmva.sourceforge.net/.
- A. Hocker et al., PoS ACAT (2007) 040, physics/0703039.
- A. Jacholkowska, J. Kalinowski, and Z. Was, Comput. Phys. Commun. 124 (2000) 238–242, hep-ph/9905225.
- J. H. Kuhn and E. Mirkes, Z. Phys. C56 (1992) 661–672.
- Z. Was, Acta Phys. Polon. B39 (2008) 1761, 0807.2775.
- P. Golonka and Z. Was, Eur. Phys. J. C50 (2007) 53–62, hep-ph/0604232.
- G. Nanava and Z. Was, Eur. Phys. J. C51 (2007) 569–583, hep-ph/0607019.
- P. Golonka, G. Nanava, and Z. Was, Tests of PHOTOS Hard Bremsstrahlung, http://mc-tester.web.cern.ch/MC-TESTER/PHOTOS-MCTESTER/.
- Z. Was and Q. Xu, In preparation.
- A. van Hameren and Z. Was, 0802.2182.
- M. Dobbs and J. B. Hansen, Comput. Phys. Commun. 134 (2001) 41–46, https://savannah.cern.ch/projects/hepmc/.
- P. Golonka, T. Pierzchala, and Z. Was, Comput. Phys. Commun. 157 (2004) 39–62, hep-ph/0210252.
- N. Davidson, P. Golonka, T. Przedzinski, and Z. Was, 0812.3215.
- P. Golonka and Z. Was, Eur. Phys. J. C45 (2006) 97–107, hep-ph/0506026.
- N. Davidson, T. Przedzinski, E. Richter-Was, and Z. Was, Presently the code is evolving quickly. Interested pilot users can download code from http://www.ph.unimelb.edu.au/ ~ndavidson/tauola/doxygen/index.html It is automatically updated daily from our source code repository. Version information can be found in a text file contained in the distribution.
- T. Pierzchała, E. Richter-Wa̧s, Z. Wa̧s, and M. Worek, Acta Phys. Polon. B32 (2001) 1277–1296, hep-ph/0101311.
- Z. Was and P. Golonka, Nucl. Phys. Proc. Suppl. 144 (2005) 88–94, hep-ph/0411377.