f_Decay.c

00001 #include "f_Decay.h"
00002 
00003 namespace Tauolapp
00004 {
00005 
00006 void TauolaDecay(int sign_type, double *polx, double* poly, 
00007                  double * polz, double * poln){
00008 
00009   taupos_.npa=1; //tau position in particle list
00010   taupos_.npb=1;
00011 
00012   double pol[4]={0};
00013   dekay_(&sign_type, pol);
00014   *polx=pol[0];
00015   *poly=pol[1];
00016   *polz=pol[2];
00017   *poln=pol[3];
00018 
00019   //  std::cout << "Polarimetric: "<<pol[0]<<","<<pol[1]
00020   //            <<","<<pol[2]<<","<<pol[3]<<std::endl;
00021   
00022 
00023 }
00024 
00025 void TauolaWriteDecayToEventRecord(int sign_type){
00026   taupos_.npa=1; //tau position in particle list
00027   taupos_.npb=1;
00028 
00029   double pol[4]={0};
00030 
00031   sign_type+=10;
00032   dekay_(&sign_type, pol); //write to event record
00033 
00034 }
00035 
00036 } // namespace Tauolapp