TEvtGen/Tauola/f_Decay.h

   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