]>
Commit | Line | Data |
---|---|---|
e1938fed | 1 | #include "f_InitTauola.h" |
2 | #include "TauolaLog.h" | |
0ca57c2f | 3 | #include <iostream> |
4 | using std::endl; | |
5 | ||
6 | namespace Tauolapp | |
7 | { | |
8 | ||
9 | void f_interface_tauolaInitialize(int pdg_id, int firstDecayMode, | |
10 | int secondDecayMode, bool rad, | |
11 | double rad_cut_off, double iniphy){ | |
12 | ||
13 | double pol[4]={0,0,0,0}; //dummy | |
14 | ||
15 | jaki_.jak1=firstDecayMode; | |
16 | jaki_.jak2=secondDecayMode; | |
17 | ||
18 | taurad_.itdkrc=rad; | |
19 | taurad_.xk0dec=rad_cut_off; //<-this one was causing nan's | |
20 | idfc_.idff=pdg_id; | |
21 | ||
22 | inimas_(); | |
23 | initdk_(); | |
24 | float iniphy_param = iniphy; | |
25 | iniphy_(&iniphy_param); | |
26 | ||
27 | int init_state=-1; | |
28 | // dexay_(&init_state,pol); | |
29 | dekay_(&init_state,pol); | |
30 | } | |
31 | ||
32 | double f_getTauMass(){ | |
33 | return (double) parmas_.amtau; | |
34 | } | |
35 | ||
36 | void f_interface_tauolaInitialise(int pdg_id, int firstDecayMode, | |
37 | int secondDecayMode, bool rad, | |
38 | double rad_cut_off, double iniphy) | |
39 | { | |
40 | Log::Warning() <<"Deprecated routine 'f_interface_tauolaInitialise'"<<endl; | |
41 | Log::Warning(0)<<"Use 'f_interface_tauolaInitialize' instead."<<endl; | |
42 | ||
43 | f_interface_tauolaInitialize(pdg_id, firstDecayMode, | |
44 | secondDecayMode, rad, | |
45 | rad_cut_off, iniphy); | |
46 | ||
47 | // Deprecated routines: initialise, setInitialisePhy, | |
48 | // f_interface_tauolaInitialise | |
49 | } | |
50 | ||
51 | } // namespace Tauolapp |