C++ Interface to Tauola
Description of Tauola Interface in C++
- Nadia Davidson, Gizo Nanava, Tomasz Przedzinski, Elzbieta Richter-Was, Zbigniew Was
The source code and documentation for release 1.1.1a
- NOTE: for downloads Nov 20-27 2012 please check note for version 1.1.1.
- Note that LCG/Genser Generator Services Subproject distributes compiled, platform adapted tar balls of our programs, in particular TAUOLA.1.0.2, TAUOLA.1.0.4 and TAUOLA.1.0.7.
The source code and documentation are updated daily from the repository as well. The following files are provided for download of the developement version:
The source code and documentation for release 1.0. The following files are provided for download:
The source code and documentation for release 1.0.2
The source code and documentation for release 1.0.4
The source code and documentation for release 1.0.5
The source code and documentation for release 1.0.6
The source code and documentation for release 1.0.7
The source code and documentation for release 1.1.0
The source code and documentation for release 1.1.1
- NOTE: this version was erroneously made public for three days (Nov 20-22) and lacks numerical protecting code against rounding errors for TauSpinner operation on events with bremsstrahlung in leptonic decays. Version 1.1.1a should be used instead. For the period of time Nov 23 - Nov 27 this (1.1.1) version number was attributed to the code of version 1.1.1a but svn label was not updated at that time.
At present (since Feb 2 2010) the C++ interface functionality for TAUOLA is complete. Longitudinal spin correlation are checked using MC-TESTER now. Transverse spin correlations and full functionality of TAUOLA Universal Interface are coded. Genuine electroweak corrections for processes mediated by Z/gamma are implemented as well.
The tar file contains the C++ interface along with the source code for tauola itself (as available from old web page; version Oct 11 2005). The developement version contains the latest source code for the interface from our subversion repository. Note that revision numbers, dates and other info for this development version can be found in the revision information file listed above.
At present, the tar ball includes everything that is needed for installation. The user is advised to use HepMC 2.04+ libraries as described below.
For compilation, and to run the simple example, the interface requires:
For further examples, one need to also install:
- ROOT v5.18 or later
- PYTHIA 8.1 or later. PYTHIA must be compiled with HepMC 2 so that the PYTHIA library hepmcinterface exists.
- MC-TESTER v1.24 or later. Do not forget to compile the additional HepMC library libHepMCEvent as well.
Examples of TAUOLA C++ combined with PHOTOS and its C++ interface are available from:
In order to compile the TAUOLA C++ interface:
- Execute './configure' with additional command line options:
- '--with-hepmc=<path> ' provides the path to HepMC installation directory. One can set HEPMCLOCATION variable instead of using this directive. This options is required for interface to compile. To compile without HepMC use '--without-hepmc'
- '--prefix=<path>' provides the installation path. The 'include' and 'lib' directories will be copied there if 'make install' is executed later. If none has been provided, the default directory for installation is '/usr/local'.
- Execute 'make'
- Optionally, execute 'make install' to copy files to the directory provided during configuration.
After compiling the 'tauola-fortran' part, the TAUOLA C++ interface will be compiled and the '/lib' and '/include' directories will contain the appropriate library and include files.
In order to compile the examples, enter the 'examples' directory, and:
- execute './configure' to determine which examples can be compiled. Additional paths can be provided as command line options:
- '--with-pythia8=<path>' provides the path to the Pythia8 installation directory. One can set the PYTHIALOCATION variable instead of using this directive. This path is required for all additional examples and tests.
- '--with-mc-tester=<path>' provides the path to the MC-Tester installation directory (libHepMCEvent must be compiled as well, check MC-Tester documentation for more details). One can set the MCTESTERLOCATION variable instead of using this directive. This path is required for all additional examples and tests. It is assumed that using this option also implies that ROOT has already been installed (since it's required by MC-TESTER). The location of its binaries should be listed in the PATH variable.
- execute 'make'
Note that for examples working with PYTHIA 8.1, the PYTHIA8DATA global variable must be set (refer to instructions provided during configuration). Similarly, for examples in the examples/testing directory to work, the MCTESTERLOCATION global variable must be set. If neither PYTHIA nor MC-TESTER are present, only the simple example will be provided. The '/examples' directory will contain the compiled example files.
If you prefer not to use configuration scripts, you may find README-NO-CONFIG.txt useful.
In order to run some more specific tests both PYTHIA and MC-TESTER must be installed.
- Compile TAUOLA C++ interface as well as examples.
- Check that the appropriate system variables are set: normally set by the script configure.paths.sh [.csh] (the configuration step mentions this script).
- Enter the /examples/testing directory. Modify test.inc if needed.
- Enter the selected directory and execute 'make'.
The appropriate .root files as well as .pdf files generated by MC-TESTER will be created inside the chosen directory. You can execute 'make clobber' to clean the directory. You can also execute 'make' inside the 'TAUOLA/examples/testing' directory to run all available tests one after another.
Electroweak corrections may affect spin correlations between tau+ and tau- in a significant way. This is the case at high energies, far above WW threshold. At present we use the alpha scheme for electroweak calculations, this choice may not be optimal for the cross section weight aiming at implementation of electroweak corrections. For this, the end scheme should be consistent with the choice of the host generator. On the other hand this is not the problem for spin correlations, our prime interest. In addition, if tables are used, transverse spin correlations for processes mediated by Z and gamma are taken into account.
Changes to the tables generated by the SANC module can be implemented by modifying the interface located in the '/SANC' directory. Details regarding the structure and initialization variables of the interface are described in the documentation. In order to generate new tables:
- execute 'make' in the '/SANC' directory to compile the library and tools.
- Adjust the initialization in the SANC tables calculation if needed.
- execute 'make tables' in the '/SANC' directory to generate the tables.
- move tables to the directory from which your main program is executed.
Design of Interface. Components are described in more detail below.
The simplified version of the main program structure including the tauola decay routines would consist of:
- initialize the Monte Carlo generator and all other generators / analysis tools.
- set the appropriate parameters for the Tauola interface.
- invoke the Tauola::initialize() routine.
- for each event to be generated:
- generate an event using any Monte Carlo generator
- convert a generated event to the HepMC event record
- create a TauolaHepMCEvent from a HepMC event object
- invoke decayTaus() for the TauolaHepMCEvent object. If a tau has already been decayed by a MC generator, use undecayTaus() first.
- the HepMC event is now prepared and can be used in the analysis
- finalize the analysis
- print any summary output if needed and exit
The following is a more detailed example point by point of how a tau, or set of taus, are decayed using the interface with the HepMC event record implementation.
Starting from the main program:
- At the initialization step, appropriate configuration routines via the 'Tauola' class methods can be executed.
- The Tauola::initialize() method must be invoked, which calls the Initialize() methods in f_Init, which in turn calls the initilization routines of TAUOLA.
In the main event loop, after the HepMC event is filled by a Monte Carlo generator:
- An appropriate implementation of the abstract class TauolaEvent - TauolaHepMCEvent object must be created, taking a HepMC::GenEvent object as its parameter.
- The method decayTaus() for the created TauolaHepMCEvent object must be invoked. If the Monte Carlo generator has already performed the decays of a tau, the undecayTaus() method can be executed just before the decayTaus() method.
After executing decayTaus(), the following procedure takes place:
- The HepMC event record is traversed and a list of stable taus in the event is created.
- the taus from the list are paired with other decay products (e.g. tau neutrino or second tau). If the pair consists of two taus, the second tau is removed from the previous list.
- TauolaParticlePairs are then created from these pairings.
- For each pair, the spin density matrix is calculated using information about the production process.
- The decayTauPairs() method is called for each pair. Within this method:
- For each tau in the pair, decay() is called.
- the decay method in TauolaParticle is responsible for adding itself (the tau) to the DecayList in the first position. It then calls decay() in f_Decay which in turn calls the dekay_() routine of TAUOLA.
- dekay_() in TAUOLA generates a set of daughter particles with momentum as though the tau was decayed in its rest frame.
- The dekay_() also provides a polarimetric vector, which is used in combination with spin density matrix to calculate the spin weight. It is used to determine whether the decay is accepted or not.
- If rejected, the list of decay products is cleared and the pair is decayed anew. This way unweighting of spin effects is performed.
- Once accepted, for each tau in the pair the routine addDecayToEventRecord() is executed, which in turn calls TAUOLA dekay_() with an option set to write to the event record.
- For each daughter TAUOLA calls filhep_() in the file f_FilHep.c.
- filhep_() adds a new TauolaHepMCParticle to DecayList with all necessary information taken from TAUOLA.
- The DecayList data structure translates the parent index given by TAUOLA to a pointer to the TauolaParticle(TauolaHepMCParticle) object.
- The parents are set via the TauolaParticle methods (setMother, etc.) By setting the mothers it is automatically assigned a barcode and added to the HepMC event.
- The particle is boosted by the tau's momentum.
- The position of the vertex of the particle is modified accordingly to the tau lifetime.
- Control is then returned to TauolaHepMCEvent which proceeds onto the next pair.
- When all pairs are decayed, the TauolaHepMCEvent converts the HepMC::GenEvent to the units selected by the user (GeV/MeV mm/cm) and finalizes execution. The loop in the main program proceeds to the next event.
This section in future will describe examples of events, particularly unusual ones, which the TAUOLA Interface will nonetheless need to be able to handle, but which were not discussed at the time of writing the documentation. At present this point is basically empty, except one example. Such type of events is already discussed in the documentation.
See software documentation in pdf form.
Last update 19 November 2012.