TauolaWrapper.h

Go to the documentation of this file.

#ifndef gen_ExternalDecays_TauolaWrapper_h
#define gen_ExternalDecays_TauolaWrapper_h

/********** TauolaWrapper
 *
 * Wrapper to Fortran functions in TAUOLA tau decay library
 *
 * Christian Veelken
 * 04/17/07
 *
 * Modified to contain access to individual decays (COMMON TAUBRA), and
 * force polarization like in case of particle gun taus (SUBROUTINE DEXAY
 * and COMMON MOMDEC)
 * 23.2.2009/S.Lehti
 *
 * Modified to remove everything related to pythia6 or pretauola;
 * all pythia6-related functionalities and/or ties are moving back
 * to Pythia6Interface;
 * the concept will ensure transparent use of Tauola with ANY
 * multi-purpose generator
 * J.V.Yarba, Feb.26, 2009
 *
 ***************************************/

//
//-------------------------------------------------------------------------------
//

// main function(s) of TAUOLA/pretauola tau decay library

extern "C" {
  void tauola_(int*, int*);
  void tauola_srs_(int*,int*);
  void taurep_(int*);
  void ranmar_(float*,int*);
  void rmarin_(int*, int*, int*);
}
#define tauola tauola_

void inline call_tauola (int mode, int polarization) { tauola(&mode, &polarization); }

extern "C" {
  extern void dexay_(int*, float[4]);
}
#define dexay dexay_

void inline call_dexay (int mode, float polarization[4]) { dexay(&mode, polarization); }

//
//-------------------------------------------------------------------------------
//

// common block with steering parameters for CMS specific Fortran interface to TAUOLA

extern "C" {
  extern struct {
    int pjak1;
    int pjak2;
    int mdtau;
  } ki_taumod_;
}
#define ki_taumod ki_taumod_

extern "C" {
  extern struct {
    int jak1;
    int jak2;
    int itdkrc;
    int ifphot;
    int ifhadm;
    int ifhadp;
  } libra_ ;
}
#define libra libra_

extern "C" {
  extern struct {
    float gamprt[30];
    int jlist[30];
    int nchan;
  } taubra_;
}
#define taubra taubra_

extern "C" {
  extern struct {
    double q1[4];
    double q2[4];
    double p1[4];
    double p2[4];
    double p3[4];
    double p4[4];
  } momdec_;
}
#define momdec momdec_

extern "C" {
  extern struct {
    int np1;
    int np2;
  } taupos_;
}
#define taupos taupos_

#endif