EvtPhokharaDef.hh

Go to the documentation of this file.

#ifndef PHOKHARA__H
#define PHOKHARA__H
 
#include "../EvtGenBase/EvtRandom.hh"
#include "cfortran/cfortran.h"
#include <string.h>
//
// common block is easy to hide bugs, please take care of it.
//
// common blocks
typedef struct {
  double r, i;
} COMPLEX16;
 
// --- couplings, masses, momenta ---
extern "C" struct {
  double    pi, gev2nbarn, gev2pbarn, alpha, me, mmu, mpi, mpi0, Sp, mp;
  double    ebeam, ver_s, rmass[10], momenta[4][16];
  COMPLEX16 cvac_s;
  double    mnt, mKp, mKn;
  COMPLEX16 cvac_s_pi, cvac_s_pi1;
  double    mlamb, meta, GFermi, sinthW2; // new Added
} ctes_;
 
// --- cuts ---
extern "C" struct {
  double Emin, gmin, phot1cut, phot2cut, pi1cut, pi2cut;
  double piphcut, accecut, q2min, w, q2_max_c, q2_min_c;
} cuts_;
 
// --- Kaons masses and gamms
extern "C" 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_;
 
// --- Pions pi+pi- masses and gamms
extern "C" 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_;
 
// --- Lambdas
extern "C" struct {
  double alpha_lamb, ratio_lamb;
} lambda_par_;
 
// --- histograms ---
typedef struct {
  char val[20];
} titledef;
 
extern "C" struct {
  DOUBLE_PRECISION xlow[20], xup[20];
  int              bins[20];
  titledef         title[20];
} histo_;
 
// --- Maxima ---
extern "C" struct {
  DOUBLE_PRECISION Mmax[3], gross[3], klein[3], tr[3], count[3];
  int              iprint;
} maxima_;
 
// --- flags --
extern "C" struct {
  int nlo, ph0, pion, fsr, fsrnlo, ivac, FF_pion, f0_model, FF_kaon, narr_res, FF_pp, chi_sw,
      chi_pion, be_r, FF_Pgg,
      nlo2; // new Added
} flags_;
 
// phokhara10.0 newly implemented
extern "C" struct {
  double hv1, A1pi0, A2pi0, Fsmv2Hv2, A1eta, A1etaP, fsmv1, H_om;
  double F_om, F_phi, fchir, A_pi_phiom, A_eta_phiom;
  double GammaTotalrho_c, GammaTotalrho_c2, GammaTotalrho_c3;
  double GammaTotalomega_c, GammaTotalomega_c2, GammaTotalomega_c3;
  double GammaTotalphi_c, GammaTotalphi_c2, GammaTotalphi_c3;
  double Mrho_c, Momega_c, Mphi_c;
  double Mrho_pr, Momega_pr, Mphi_pr;
  double Mrho_dpr, Momega_dpr, Mphi_dpr, metaf, metaP, mpi0f;
  double Cq_c, Cs_c, Cq_P, Cs_P, A3pi0, A3eta, A3etaP;
} pscalars_;
 
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 )
 
PROTOCCALLSFSUB5( GEN_0PH, gen_0ph, INT, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE )
#define GEN_0PH( I, QQMIN, SP, COS3MIN, COS3MAX )                                             \
  CCALLSFSUB5( GEN_0PH, gen_0ph, INT, PDOUBLE, PDOUBLE, PDOUBLE, PDOUBLE, I, QQMIN, SP,       \
               COS3MIN, COS3MAX )
 
#endif