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