f_Decay.h

#ifndef _f_Decay_h_included_
#define _f_Decay_h_included_

#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);
}

void TauolaDecay(int sign_type, double *polx, double *poly,
                 double *polz, double *poln);

void TauolaWriteDecayToEventRecord(int sign_type);

} // namespace Tauolapp
#endif