]>
Commit | Line | Data |
---|---|---|
0ca57c2f | 1 | #ifndef _f_Decay_h_included_ |
2 | #define _f_Decay_h_included_ | |
3 | ||
4 | /** | |
5 | * This file contains an interface between the C++ code and TAUOLA | |
6 | * FORTRAN routines for decaying taus. TauolaDecay() should be used | |
7 | * by C++ code. This call the dexay_ or dekay_ routines defined in | |
8 | * tauola.f | |
9 | * | |
10 | * @author Nadia Davidson | |
11 | * @date 17 June 2008 | |
12 | */ | |
13 | ||
14 | #include <vector> | |
15 | #include <iostream> | |
16 | ||
17 | namespace Tauolapp | |
18 | { | |
19 | ||
20 | extern "C" { | |
21 | ||
22 | extern struct{//positions of taus in the LUND common block | |
23 | int npa; | |
24 | int npb; | |
25 | } taupos_; | |
26 | ||
27 | //extern void dexay_(int *state, double pol[4]); | |
28 | extern void dekay_(int *state, double pol[4]); | |
29 | ||
30 | extern void taupi0_(double pp[4],int *k); | |
31 | extern void tauk0s_(double pp[4],int *k); | |
32 | extern void taueta_(double pp[4],int *k); | |
33 | } | |
34 | ||
35 | /** Invokes DEKAY with "1" or "2" to get the polarization information. */ | |
36 | void TauolaDecay(int sign_type, double *polx, double *poly, | |
37 | double *polz, double *poln); | |
38 | ||
39 | /** Invokes DEKAY with "11" or "12" to produce the decay. */ | |
40 | void TauolaWriteDecayToEventRecord(int sign_type); | |
41 | ||
42 | } // namespace Tauolapp | |
43 | #endif |