TEvtGen/Tauola/f_InitTauola.cxx

   1 #include "f_InitTauola.h"
   2 #include "TauolaLog.h"
   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