f_Decay.h

00001 #ifndef _f_Decay_h_included_
00002 #define _f_Decay_h_included_
00003 
00004 /**
00005  * This file contains an interface between the C++ code and TAUOLA
00006  * FORTRAN routines for decaying taus. TauolaDecay() should be used
00007  * by C++ code. This call the dexay_ or dekay_ routines defined in
00008  * tauola.f
00009  *
00010  * @author Nadia Davidson
00011  * @date 17 June 2008
00012  */
00013 
00014 #include <vector>
00015 #include <iostream>
00016 
00017 namespace Tauolapp
00018 {
00019 
00020 extern "C" {
00021 
00022   extern struct{//positions of taus in the LUND common block
00023     int npa;
00024     int npb;
00025   } taupos_;
00026 
00027   //extern void dexay_(int *state, double pol[4]);
00028   extern void dekay_(int *state, double pol[4]);
00029 
00030   extern void taupi0_(double pp[4],int *k);
00031   extern void tauk0s_(double pp[4],int *k);
00032   extern void taueta_(double pp[4],int *k);
00033 }
00034 
00035 /** Invokes DEKAY with "1" or "2" to get the polarization information. */
00036 void TauolaDecay(int sign_type, double *polx, double *poly,
00037                  double *polz, double *poln);
00038 
00039 /** Invokes DEKAY with "11" or "12" to produce the decay. */
00040 void TauolaWriteDecayToEventRecord(int sign_type);
00041 
00042 } // namespace Tauolapp
00043 #endif