[u/mrichter/AliRoot.git] / TEvtGen / 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
Commit	Line	Data
0ca57c2f	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