TEvtGen/Tauola/f_Decay.cxx

   1 #include "f_Decay.h"
   2
   3 namespace Tauolapp
   4 {
   5
   6 void TauolaDecay(int sign_type, double *polx, double* poly,
   7                  double * polz, double * poln){
   8
   9   taupos_.npa=1; //tau position in particle list
  10   taupos_.npb=1;
  11
  12   double pol[4]={0};
  13   dekay_(&sign_type, pol);
  14   *polx=pol[0];
  15   *poly=pol[1];
  16   *polz=pol[2];
  17   *poln=pol[3];
  18
  19   //  std::cout << "Polarimetric: "<<pol[0]<<","<<pol[1]
  20   //            <<","<<pol[2]<<","<<pol[3]<<std::endl;
  21
  22
  23 }
  24
  25 void TauolaWriteDecayToEventRecord(int sign_type){
  26   taupos_.npa=1; //tau position in particle list
  27   taupos_.npb=1;
  28
  29   double pol[4]={0};
  30
  31   sign_type+=10;
  32   dekay_(&sign_type, pol); //write to event record
  33
  34 }
  35
  36 } // namespace Tauolapp