PhokharaDef.h

Go to the documentation of this file.

#ifndef PHOKHARA__H
#define PHOKHARA__H
 
#include "cfortran/cfortran.h"
#include <string.h>
 
// common blocks
typedef struct {
  double r, i;
} COMPLEX16;
 
// --- couplings, masses, momenta ---
typedef struct {
  double    pi, gev2nbarn, gev2pbarn, alpha, me, mmu, mpi, mpi0, Sp, mp;
  double    ebeam, ver_s, rmass[10], momenta[4][15];
  COMPLEX16 cvac_s;
  double    mnt, mKp, mKn;
  COMPLEX16 cvac_s_pi, cvac_s_phi;
  double    mlamb;
} CTES_DEF;
 
#define CTES COMMON_BLOCK( CTES_DEF, ctes )
COMMON_BLOCK_DEF( CTES_DEF, CTES );
 
// --- cuts ---
typedef struct {
  double Emin, gmin, phot1cut, phot2cut, pi1cut, pi2cut;
  double piphcut, accecut, q2min, w, q2_max_c, q2_min_c;
} CUTS_DEF;
 
#define CUTS COMMON_BLOCK( CUTS_DEF, cuts )
COMMON_BLOCK_DEF( CUTS_DEF, CUTS );
 
// --- Kaons masses and gamms
typedef struct {
  double m_phi0_Kp, g_phi0_Kp, m_phi1_Kp, g_phi1_Kp;
  double m_om0_Kp, g_om0_Kp, m_om1_Kp, g_om1_Kp, m_om2_Kp, g_om2_Kp;
  double m_rho0_Kp, g_rho0_Kp, m_rho1_Kp, g_rho1_Kp, m_rho2_Kp, g_rho2_Kp;
} MASS_GAMM_KP_DEF;
 
#define MASS_GAMM_KP COMMON_BLOCK( MASS_GAMM_KP_DEF, mass_gamm_Kp )
COMMON_BLOCK_DEF( MASS_GAMM_KP_DEF, MASS_GAMM_KP );
 
// --- Pions pi+pi- masses and gamms
typedef struct {
  double m_rho0_pion, g_rho0_pion, m_rho1_pion;
  double g_rho1_pion, m_rho2_pion, g_rho2_pion, m_rho3_pion, g_rho3_pion;
  double m_om0_pion, g_om0_pion;
} MASS_GAMM_PION_DEF;
 
#define MASS_GAMM_PION COMMON_BLOCK( MASS_GAMM_PION_DEF, mass_gamm_pion )
COMMON_BLOCK_DEF( MASS_GAMM_PION_DEF, MASS_GAMM_PION );
 
// --- Lambdas
typedef struct {
  double alpha_lamb, ratio_lamb;
} LAMBDA_PAR_DEF;
 
#define LAMBDA_PAR COMMON_BLOCK( LAMBDA_PAR_DEF, lambda_par )
COMMON_BLOCK_DEF( LAMBDA_PAR_DEF, LAMBDA_PAR );
 
// --- histograms ---
typedef struct {
  char val[20];
} titledef;
 
typedef struct {
  DOUBLE_PRECISION xlow[20], xup[20];
  int              bins[20];
  titledef         title[20];
} HISTO_DEF;
 
#define HISTO COMMON_BLOCK( HISTO_DEF, histo )
COMMON_BLOCK_DEF( HISTO_DEF, HISTO );
 
// --- Maxima ---
typedef struct {
  DOUBLE_PRECISION Mmax[2], gross[2], klein[2], tr[2], count[2];
  int              iprint;
} MAXIMA_DEF;
 
#define MAXIMA COMMON_BLOCK( MAXIMA_DEF, maxima )
COMMON_BLOCK_DEF( MAXIMA_DEF, MAXIMA );
 
// --- flags --
typedef struct {
  int nlo, pion, fsr, fsrnlo, ivac, FF_pion, f0_model, FF_kaon, narr_res;
} FLAGS_DEF;
 
#define FLAGS COMMON_BLOCK( FLAGS_DEF, flags )
COMMON_BLOCK_DEF( FLAGS_DEF, FLAGS );
 
PROTOCCALLSFSUB2( RLXDINIT, rlxdinit, INT, INT )
#define RLXDINIT( LUXURY, SEED ) CCALLSFSUB2( RLXDINIT, rlxdinit, INT, INT, LUXURY, SEED )
 
PROTOCCALLSFSUB1( RLXDRESETF, rlxdresetf, INTV )
#define RLXDRESETF( SEED ) CCALLSFSUB1( RLXDRESETF, rlxdresetf, INTV, SEED )
 
PROTOCCALLSFSUB0( INPUT, input )
#define INPUT() CCALLSFSUB0( INPUT, input )
 
// PROTOCCALLSFSUB3(INPUT,input,PLONG,PINT,PSTRING);
// #define INPUT(NGES,NM,OUTFILE) CCALLSFSUB3(INPUT,input,PLONG,PINT,PSTRING,NGES,NM,OUTFILE)
 
// PROTOCCALLSFSUB0(INITHISTO,inithisto);
// #define INITHISTO() CCALLSFSUB0(INITHISTO,inithisto)
 
// PROTOCCALLSFSUB0(ENDHISTO,endhisto);
// #define ENDHISTO() CCALLSFSUB0(ENDHISTO,endhisto)
 
// PROTOCCALLSFSUB0(WRITEEVENT,writeevent);
// #define WRITEEVENT() CCALLSFSUB0(WRITEEVENT,writeevent)
 
PROTOCCALLSFSUB2( RANLXDF, ranlxdf, DOUBLEV, INT )
#define RANLXDF( AR, VAL ) CCALLSFSUB2( RANLXDF, ranlxdf, DOUBLEV, INT, AR, VAL )
 
PROTOCCALLSFSUB7( GEN_1PH, gen_1ph, INT, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE )
#define GEN_1PH( I, QQMIN, QQMAX, COS1MIN, COS1MAX, COS3MIN, COS3MAX )                        \
  CCALLSFSUB7( GEN_1PH, gen_1ph, INT, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE,   \
               I, QQMIN, QQMAX, COS1MIN, COS1MAX, COS3MIN, COS3MAX )
 
PROTOCCALLSFSUB8( GEN_2PH, gen_2ph, INT, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE,
                  PDOUBLE )
#define GEN_2PH( I, QQMIN, COS1MIN, COS1MAX, COS2MIN, COS2MAX, COS3MIN, COS3MAX )             \
  CCALLSFSUB8( GEN_2PH, gen_2ph, INT, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE,   \
               PDOUBLE, I, QQMIN, COS1MIN, COS1MAX, COS2MIN, COS2MAX, COS3MIN, COS3MAX )
 
#endif