include/Tauola/f_Decay.h

#ifndef _f_Decay_h_included_
#define _f_Decay_h_included_
 
/**
 * This file contains an interface between the C++ code and TAUOLA
 * FORTRAN routines for decaying taus. TauolaDecay() should be used
 * by C++ code. This call the dexay_ or dekay_ routines defined in
 * tauola.f
 *
 * @author Nadia Davidson
 * @date 17 June 2008
 */
 
#include <vector>
#include <iostream>
 
namespace Tauolapp
{
 
extern "C" {
 
  extern struct{//positions of taus in the LUND common block
    int npa;
    int npb;
  } taupos_;
 
  //extern void dexay_(int *state, double pol[4]);
  extern void dekay_(int *state, double pol[4]);
 
  extern void taupi0_(double pp[4],int *k);
  extern void tauk0s_(double pp[4],int *k);
  extern void taueta_(double pp[4],int *k);
}
 
/** Invokes DEKAY with "1" or "2" to get the polarization information. */
void TauolaDecay(int sign_type, double *polx, double *poly,
                 double *polz, double *poln);
 
/** Invokes DEKAY with "11" or "12" to produce the decay. */
void TauolaWriteDecayToEventRecord(int sign_type);
 
} // namespace Tauolapp
#endif