[u/mrichter/AliRoot.git] / THerwig / Herwig6Calls.h

#ifndef ROOT_Herwig6Calls
#define ROOT_Herwig6Calls

// Translation of Fortran commons from the Herwig6
// f77 program into c++ structures to be used in ROOT
// and declaration of Fortran functions as extern
// C functions to be called from the class Herwig6
// author: j. g. contreras jgcn@moni.mda.cinvestav.mx
// date: december 22, 2000

struct dcpx {double dr,di;};

int const NMXHEP = 4000;

struct Hepevt_t {
  int NEVHEP;
  int NHEP;
  int ISTHEP[NMXHEP];
  int IDHEP[NMXHEP];
  int JMOHEP[NMXHEP][2];
  int JDAHEP[NMXHEP][2];
  double PHEP[NMXHEP][5];
  double VHEP[NMXHEP][4];
};

struct Hwbeam_t {
  int IPART1;
  int IPART2;
};

struct Hwbmch_t {
  char PART1[8];
  char PART2[8];
};

struct Hwproc_t {
  double EBEAM1;
  double EBEAM2;
  double PBEAM1;
  double PBEAM2;
  int    IPROC;
  int    MAXEV;
};

struct Hwpram_t {
  double AFCH[2][16];
  double ALPHEM;
  double B1LIM;
  double BETAF;
  double BTCLM;
  double CAFAC;
  double CFFAC;
  double CLMAX;
  double CLPOW;
  double CLSMR[2];
  double CSPEED;
  double ENSOF;
  double ETAMIX;
  double F0MIX;
  double F1MIX;
  double F2MIX;
  double GAMH;
  double GAMW;
  double GAMZ;
  double GAMZP;
  double GEV2NB;
  double H1MIX;
  double PDIQK;
  double PGSMX;
  double PGSPL[4];
  double PHIMIX;
  double PIFAC;
  double PRSOF;
  double PSPLT[2];
  double PTRMS;
  double PXRMS;
  double QCDL3;
  double QCDL5;
  double QCDLAM;
  double QDIQK;
  double QFCH[16];
  double QG;
  double QSPAC;
  double QV;
  double SCABI;
  double SWEIN;
  double TMTOP;
  double VFCH[2][16];
  double VCKM[3][3];
  double VGCUT;
  double VQCUT;   
  double VPCUT;
  double ZBINM;
  double EFFMIN;
  double OMHMIX;
  double ET2MIX;
  double PH3MIX;
  double GCUTME;
  int    IOPREM;
  int    IPRINT;
  int    ISPAC;
  int    LRSUD;
  int    LWSUD;
  int    MODPDF[2];
  int    NBTRY;
  int    NCOLO;
  int    NCTRY;
  int    NDTRY;
  int    NETRY;
  int    NFLAV;
  int    NGSPL;
  int    NSTRU;
  int    NSTRY;
  int    NZBIN;
  int    IOP4JT[2];
  int    NPRFMT;
  int AZSOFT;
  int AZSPIN;
  int CLDIR[2];
  int HARDME;
  int NOSPAC;
  int PRNDEC;
  int PRVTX;
  int SOFTME;
  int ZPRIME;
  int PRNDEF;
  int PRNTEX;
  int PRNWEB;
};

struct Hwprch_t {
  char AUTPDF[2][20];
  char BDECAY[4];
};

int const NMXPAR = 500;

struct Hwpart_t {
  int  NEVPAR;
  int  NPAR;
  int  ISTPAR[NMXPAR];
  int  IDPAR[NMXPAR];
  int  JMOPAR[NMXPAR][2];
  int  JDAPAR[NMXPAR][2];
  double  PPAR[NMXPAR][5];
  double  VPAR[NMXPAR][4];
};

struct Hwparp_t {
  double DECPAR[NMXPAR][2];
  double PHIPAR[NMXPAR][2];
  double RHOPAR[NMXPAR][2];
  int TMPAR[NMXPAR];
};

int const MODMAX = 50;

struct Hwbosc_t {
  double  ALPFAC;
  double  BRHIG[12];
  double  ENHANC[12];
  double  GAMMAX;
  double  RHOHEP[NMXHEP][3];
  int     IOPHIG;
  int     MODBOS[MODMAX];
};

struct Hwparc_t {
  int     JCOPAR[NMXPAR][4];
};

struct Hwbrch_t {
  double ANOMSC[2][2];
  double HARDST;
  double PTINT[2][3];
  double XFACT;
  int    INHAD;
  int    JNHAD;
  int    NSPAC[7];
  int    ISLENT;
  int    BREIT;
  int    FROST;
  int    USECMF;
};

struct Hwevnt_t {
  double AVWGT;
  double EVWGT;
  double GAMWT;
  double TLOUT;
  double WBIGST;
  double WGTMAX;
  double WGTSUM;
  double WSQSUM;
  int    IDHW[NMXHEP];
  int    IERROR;
  int    ISTAT;
  int    LWEVT;
  int    MAXER;
  int    MAXPR;
  int    NOWGT;
  int    NRN[2];
  int    NUMER;
  int    NUMERU;
  int    NWGTS;
  int    GENSOF;
};

struct Hwhard_t {
  double ASFIXD;
  double CLQ[6][7];
  double COSS;
  double COSTH;
  double CTMAX;
  double DISF[2][13];
  double EMLST;
  double EMMAX;
  double EMMIN;
  double EMPOW;
  double EMSCA;
  double EPOLN[3];
  double GCOEF[7];
  double GPOLN;
  double OMEGA0;
  double PHOMAS;
  double PPOLN[3];
  double PTMAX;
  double PTMIN;
  double PTPOW;
  double Q2MAX;
  double Q2MIN;
  double Q2POW;
  double Q2WWMN;
  double Q2WWMX;
  double QLIM;
  double SINS;
  double THMAX;
  double Y4JT;
  double TMNISR;
  double TQWT;
  double XX[2];
  double XLMIN;
  double XXMIN;
  double YBMAX;
  double YBMIN;
  double YJMAX;
  double YJMIN;
  double YWWMAX;
  double YWWMIN;
  double WHMIN;
  double ZJMAX;
  double ZMXISR;
  int    IAPHIG;
  int    IBRN[2];
  int    IBSH;
  int    ICO[10];
  int    IDCMF;
  int    IDN[10];
  int    IFLMAX;
  int    IFLMIN;
  int    IHPRO;
  int    IPRO;
  int    MAPQ[10];
  int    MAXFL;
  int    BGSHAT;
  int    COLISR;
  int    FSTEVT;
  int    FSTWGT;
  int    GENEV;
  int    HVFCEN;
  int    TPOL;
  int     DURHAM;
};

int const NMXRES = 500;

struct Hwprop_t {
  double RLTIM[NMXRES+1];
  double RMASS[NMXRES+1];
  double RSPIN[NMXRES+1];
  int    ICHRG[NMXRES+1];
  int    IDPDG[NMXRES+1];
  int    IFLAV[NMXRES+1];
  int    NRES;
  int    VTOCDK[NMXRES+1];
  int    VTORDK[NMXRES+1];
  int    QORQQB[NMXRES+1];
  int    QBORQQ[NMXRES+1];
};

struct Hwunam_t {
  char RNAME[NMXRES+1][8];
  char TXNAME[NMXRES+1][2][37];
};

int const NMXDKS = 4000;
int const NMXMOD = 200;

struct Hwupdt_t {
  double BRFRAC[NMXDKS];
  double CMMOM[NMXDKS];
  double DKLTM[NMXRES];
  int    IDK[NMXDKS];
  int    IDKPRD[NMXDKS][5];
  int    LNEXT[NMXDKS];
  int    LSTRT[NMXRES];
  int    NDKYS;
  int    NME[NMXDKS];
  int    NMODES[NMXRES];
  int    NPRODS[NMXDKS];
  int    DKPSET;
  int    RSTAB[NMXRES+1];
};

struct Hwuwts_t {
  double REPWT[5][4][4];
  double SNGWT;
  double DECWT;
  double QWT[3];
  double PWT[12];
  double SWTEF[NMXRES];
};

int const NMXCDK = 4000;

struct Hwuclu_t {
  double CLDKWT[NMXCDK];
  double CTHRPW[12][12];
  double PRECO;
  double RESN[12][12];
  double RMIN[12][12];
  int    LOCN[12][12];
  int    NCLDK[NMXCDK];
  int    NRECO;
  int    CLRECO;
};

struct Hwdist_t {
  double EXAG;
  double GEV2MM;
  double HBAR;
  double PLTCUT;
  double VMIN2;
  double VTXPIP[5];
  double XMIX[2];
  double XMRCT[2];
  double YMIX[2];
  double YMRCT[2];
  int    IOPDKL;
  int    MAXDKL;
  int    MIXING;
  int    PIPSMR;
};

int const NMXQDK=20;

struct Hwqdks_t {
  double VTXQDK[NMXQDK][4];
  int    IMQDK[NMXQDK];
  int    LOCQ[NMXQDK];
  int    NQDK;
};

int const NMXSUD = 1024;

struct Hwusud_t {
  double ACCUR;
  double QEV[6][NMXSUD];
  double SUD[6][NMXSUD];
  int    INTER;
  int    NQEV;
  int    NSUD;
  int    SUDORD;
};

struct Hwsusy_t {
  double TANB;
  double ALPHAH;
  double COSBPA;
  double SINBPA;
  double COSBMA;
  double SINBMA;
  double COSA;
  double SINA;
  double COSB;
  double SINB;
  double COTB;
  double ZMIXSS[4][4];
  double ZMXNSS[4][4];
  double ZSGNSS[4]; 
  double LFCH[16];
  double RFCH[16];
  double SLFCH[4][16];
  double SRFCH[4][16]; 
  double WMXUSS[2][2];
  double WMXVSS[2][2]; 
  double WSGNSS[2];
  double QMIXSS[2][2][6];
  double LMIXSS[2][2][6];
  double THETAT;
  double THETAB;
  double THETAL;
  double ATSS;
  double ABSS;
  double ALSS;
  double MUSS;
  double FACTSS;
  double GHWWSS[3];
  double GHZZSS[3];
  double GHDDSS[4];
  double GHUUSS[4];
  double GHWHSS[3];
  double GHSQSS[2][2][6][4];
  double XLMNSS;
  double RMMNSS;
  double DMSSM;
  double SENHNC[24];
  double SSPARITY;
  int    SUSYIN;
};

struct Hwrpar_t {
  double LAMDA1[3][3][3];
  double LAMDA2[3][3][3];
  double LAMDA3[3][3][3];
  int    HRDCOL[5][2];
  int    RPARTY;
  int    COLUPD;
};

struct Hwminb_t {
  double PMBN1;
  double PMBN2;
  double PMBN3;
  double PMBK1;
  double PMBK2;
  double PMBM1;
  double PMBM2;
  double PMBP1;
  double PMBP2;
  double PMBP3;
};

int const NMXCL = 500;

struct Hwclus_t {
  double PPCL[NMXCL][5];
  int    IDCL[NMXCL];
  int    NCL;
};

// herwig 6.507

struct Hwgrav_t {
  double GRVLAM;
  double EMGRV;
  double GAMGRV;
};

struct Hw6202_t {
  double VIPWID[3];
  double DXRCYL;
  double DXZMAX;
  double DXRSPH;
  int    WZRFR;
  int    FIX4JT;
  int    IMSSM;
  int    IHIGGS;
  int    PARITY;
  int    LRSUSY;
};

struct Hw6203_t {
  double ABWGT;
  double ABWSUM;
  double AVABW;
  int    NNEGWT;
  int    NNEGEV;
  int    NEGWTS;
};

int const IMAXCH = 20;

struct Hw6300_t {
  double MJJMIN;
  double CHNPRB[IMAXCH];
  int    IOPSTP;
  int    IOPSH;
  int    OPTM;
  int    CHON[IMAXCH];
};


int const NXMRS = 49;
int const NQMRS = 37;
int const NPMRS = 8;

struct Hwpmrs_t {
  double FMRS[NQMRS+1][NXMRS][NPMRS][3];
};

struct Hwcirc_t {
  int	   CIRCOP;
  int	   CIRCAC;
  int	   CIRCVR;
  int	   CIRCRV;
  int	   CIRCCH;
};

int const NCFMAX = 3;
int const NMODE2 = 500;
int const NMODE3 = 500;
int const NDIAGR = 8;
int const NMODEB = 50;
int const NMODE4 = 4;

struct Hwdspb_t {
  double ABMODE[NMODEB][2];
  double BBMODE[NMODEB][12][2];
  double PBMODE[NMODEB][12];
  double WTBMAX[NMODEB][12];
  int    IDBPRT[NMODEB];
  int    IBDRTP[NMODEB];
  int    IBMODE[NMODEB];
  int    NBMODE;
};

struct Hwdsp2_t {
  double A2MODE[NMODE2][2];
  double P2MODE[NMODE2];
  double WT2MAX[NMODE2];
  int    ID2PRT[NMODE2];
  int    I2DRTP[NMODE2];
  int    N2MODE;
};

struct Hwdsp3_t {
  double A3MODE[NMODE3][NDIAGR][2];
  double B3MODE[NMODE3][NDIAGR][2];
  double P3MODE[NMODE3];
  double WT3MAX[NMODE3];
  double SPN3CF[NMODE3][NCFMAX][NCFMAX];
  int    ID3PRT[NMODE3];
  int    I3MODE[NDIAGR][NMODE3];
  int    I3DRTP[NDIAGR][NMODE3];
  int    N3MODE;
  int    NDI3BY[NMODE3];
  int    N3NCFL[NMODE3];
  int    I3DRCF[NDIAGR][NMODE3];
};

struct Hwdsp4_t {
  double A4MODE[NMODE4][12][2];
  double B4MODE[NMODE4][12][2];
  double P4MODE[NMODE4][12][12];
  double WT4MAX[NMODE4][12][12];
  int    ID4PRT[NMODE4];
  int    I4MODE[NMODE4][2];
  int    N4MODE;
};

struct Hwdspn_t {
  int    NDECSY;
  int    NSEARCH;
  int    LRDEC;
  int    LWDEC;
  int    SYSPIN;
  int    THREEB;
  int    FOURB;
  char   TAUDEC[6];
};

int const NMXSPN = 50;

struct Hwspin_t {
  dcpx   MESPN[NMXSPN][NCFMAX][2][2][2][2];
  dcpx   RHOSPN[NMXSPN][2][2];
  double SPNCFC[NMXSPN][NCFMAX][NCFMAX];
  int    IDSPN[NMXSPN];
  int    JMOSPN[NMXSPN];
  int    JDASPN[NMXSPN][2];
  int    NSPN;
  int    ISNHEP[NMXHEP];
  int    NSNTRY;
  int    DECSPN[NMXSPN];
  int    NCFL[NMXSPN];
  int    SPCOPT;
};

struct Hwstau_t {
  int    JAK1;
  int    JAK2;
  int    ITDKRC;
  int    IFPHOT;
};

int const MAXHRP = 100;

struct Hwgupr_t {
  double LHWGT[MAXHRP];
  double LHWGTS[MAXHRP];
  double LHXSCT[MAXHRP];
  double LHXERR[MAXHRP];
  double LHXMAX[MAXHRP];
  double LHMXSM;
  int    LHIWGT[MAXHRP];
  int    LHNEVT[MAXHRP];
  int    ITYPLH;
  int    LHSOFT;
  int    LHGLSF;    
};

struct Hw6500_t {
  int    PRESPL;
};

struct Hw6504_t {
  int    ITOPRD;
};

struct Hw6506_t {
  double PDFX0;
  double PDFPOW;
};


extern "C" {
  void  hwigin_();
  void  hwuinc_();
  void  hwusta_(char * name, int);
  void  hweini_();
  void  hwuine_();
  void  hwepro_();
  void  hwbgen_();
  void  hwdhob_();
  void  hwcfor_();
  void  hwcdec_();
  void  hwdhad_();
  void  hwdhvy_();
  void  hwmevt_();
  void  hwufne_();
  void  hwefin_();
}

// subroutines to be call by JIMMY

extern "C" {
  void  jminit_();
  void  jimmin_();
  void  jmefin_();
}


#endif
Commit	Line	Data
	1	#ifndef ROOT_Herwig6Calls
	2	#define ROOT_Herwig6Calls
	3
	4	// Translation of Fortran commons from the Herwig6
	5	// f77 program into c++ structures to be used in ROOT
	6	// and declaration of Fortran functions as extern
	7	// C functions to be called from the class Herwig6
	8	// author: j. g. contreras jgcn@moni.mda.cinvestav.mx
	9	// date: december 22, 2000
	10
	11	struct dcpx {double dr,di;};
	12
	13	int const NMXHEP = 4000;
	14
	15	struct Hepevt_t {
	16	int NEVHEP;
	17	int NHEP;
	18	int ISTHEP[NMXHEP];
	19	int IDHEP[NMXHEP];
	20	int JMOHEP[NMXHEP][2];
	21	int JDAHEP[NMXHEP][2];
	22	double PHEP[NMXHEP][5];
	23	double VHEP[NMXHEP][4];
	24	};
	25
	26	struct Hwbeam_t {
	27	int IPART1;
	28	int IPART2;
	29	};
	30
	31	struct Hwbmch_t {
	32	char PART1[8];
	33	char PART2[8];
	34	};
	35
	36	struct Hwproc_t {
	37	double EBEAM1;
	38	double EBEAM2;
	39	double PBEAM1;
	40	double PBEAM2;
	41	int IPROC;
	42	int MAXEV;
	43	};
	44
	45	struct Hwpram_t {
	46	double AFCH[2][16];
	47	double ALPHEM;
	48	double B1LIM;
	49	double BETAF;
	50	double BTCLM;
	51	double CAFAC;
	52	double CFFAC;
	53	double CLMAX;
	54	double CLPOW;
	55	double CLSMR[2];
	56	double CSPEED;
	57	double ENSOF;
	58	double ETAMIX;
	59	double F0MIX;
	60	double F1MIX;
	61	double F2MIX;
	62	double GAMH;
	63	double GAMW;
	64	double GAMZ;
	65	double GAMZP;
	66	double GEV2NB;
	67	double H1MIX;
	68	double PDIQK;
	69	double PGSMX;
	70	double PGSPL[4];
	71	double PHIMIX;
	72	double PIFAC;
	73	double PRSOF;
	74	double PSPLT[2];
	75	double PTRMS;
	76	double PXRMS;
	77	double QCDL3;
	78	double QCDL5;
	79	double QCDLAM;
	80	double QDIQK;
	81	double QFCH[16];
	82	double QG;
	83	double QSPAC;
	84	double QV;
	85	double SCABI;
	86	double SWEIN;
	87	double TMTOP;
	88	double VFCH[2][16];
	89	double VCKM[3][3];
	90	double VGCUT;
	91	double VQCUT;
	92	double VPCUT;
	93	double ZBINM;
	94	double EFFMIN;
	95	double OMHMIX;
	96	double ET2MIX;
	97	double PH3MIX;
	98	double GCUTME;
	99	int IOPREM;
	100	int IPRINT;
	101	int ISPAC;
	102	int LRSUD;
	103	int LWSUD;
	104	int MODPDF[2];
	105	int NBTRY;
	106	int NCOLO;
	107	int NCTRY;
	108	int NDTRY;
	109	int NETRY;
	110	int NFLAV;
	111	int NGSPL;
	112	int NSTRU;
	113	int NSTRY;
	114	int NZBIN;
	115	int IOP4JT[2];
	116	int NPRFMT;
	117	int AZSOFT;
	118	int AZSPIN;
	119	int CLDIR[2];
	120	int HARDME;
	121	int NOSPAC;
	122	int PRNDEC;
	123	int PRVTX;
	124	int SOFTME;
	125	int ZPRIME;
	126	int PRNDEF;
	127	int PRNTEX;
	128	int PRNWEB;
	129	};
	130
	131	struct Hwprch_t {
	132	char AUTPDF[2][20];
	133	char BDECAY[4];
	134	};
	135
	136	int const NMXPAR = 500;
	137
	138	struct Hwpart_t {
	139	int NEVPAR;
	140	int NPAR;
	141	int ISTPAR[NMXPAR];
	142	int IDPAR[NMXPAR];
	143	int JMOPAR[NMXPAR][2];
	144	int JDAPAR[NMXPAR][2];
	145	double PPAR[NMXPAR][5];
	146	double VPAR[NMXPAR][4];
	147	};
	148
	149	struct Hwparp_t {
	150	double DECPAR[NMXPAR][2];
	151	double PHIPAR[NMXPAR][2];
	152	double RHOPAR[NMXPAR][2];
	153	int TMPAR[NMXPAR];
	154	};
	155
	156	int const MODMAX = 50;
	157
	158	struct Hwbosc_t {
	159	double ALPFAC;
	160	double BRHIG[12];
	161	double ENHANC[12];
	162	double GAMMAX;
	163	double RHOHEP[NMXHEP][3];
	164	int IOPHIG;
	165	int MODBOS[MODMAX];
	166	};
	167
	168	struct Hwparc_t {
	169	int JCOPAR[NMXPAR][4];
	170	};
	171
	172	struct Hwbrch_t {
	173	double ANOMSC[2][2];
	174	double HARDST;
	175	double PTINT[2][3];
	176	double XFACT;
	177	int INHAD;
	178	int JNHAD;
	179	int NSPAC[7];
	180	int ISLENT;
	181	int BREIT;
	182	int FROST;
	183	int USECMF;
	184	};
	185
	186	struct Hwevnt_t {
	187	double AVWGT;
	188	double EVWGT;
	189	double GAMWT;
	190	double TLOUT;
	191	double WBIGST;
	192	double WGTMAX;
	193	double WGTSUM;
	194	double WSQSUM;
	195	int IDHW[NMXHEP];
	196	int IERROR;
	197	int ISTAT;
	198	int LWEVT;
	199	int MAXER;
	200	int MAXPR;
	201	int NOWGT;
	202	int NRN[2];
	203	int NUMER;
	204	int NUMERU;
	205	int NWGTS;
	206	int GENSOF;
	207	};
	208
	209	struct Hwhard_t {
	210	double ASFIXD;
	211	double CLQ[6][7];
	212	double COSS;
	213	double COSTH;
	214	double CTMAX;
	215	double DISF[2][13];
	216	double EMLST;
	217	double EMMAX;
	218	double EMMIN;
	219	double EMPOW;
	220	double EMSCA;
	221	double EPOLN[3];
	222	double GCOEF[7];
	223	double GPOLN;
	224	double OMEGA0;
	225	double PHOMAS;
	226	double PPOLN[3];
	227	double PTMAX;
	228	double PTMIN;
	229	double PTPOW;
	230	double Q2MAX;
	231	double Q2MIN;
	232	double Q2POW;
	233	double Q2WWMN;
	234	double Q2WWMX;
	235	double QLIM;
	236	double SINS;
	237	double THMAX;
	238	double Y4JT;
	239	double TMNISR;
	240	double TQWT;
	241	double XX[2];
	242	double XLMIN;
	243	double XXMIN;
	244	double YBMAX;
	245	double YBMIN;
	246	double YJMAX;
	247	double YJMIN;
	248	double YWWMAX;
	249	double YWWMIN;
	250	double WHMIN;
	251	double ZJMAX;
	252	double ZMXISR;
	253	int IAPHIG;
	254	int IBRN[2];
	255	int IBSH;
	256	int ICO[10];
	257	int IDCMF;
	258	int IDN[10];
	259	int IFLMAX;
	260	int IFLMIN;
	261	int IHPRO;
	262	int IPRO;
	263	int MAPQ[10];
	264	int MAXFL;
	265	int BGSHAT;
	266	int COLISR;
	267	int FSTEVT;
	268	int FSTWGT;
	269	int GENEV;
	270	int HVFCEN;
	271	int TPOL;
	272	int DURHAM;
	273	};
	274
	275	int const NMXRES = 500;
	276
	277	struct Hwprop_t {
	278	double RLTIM[NMXRES+1];
	279	double RMASS[NMXRES+1];
	280	double RSPIN[NMXRES+1];
	281	int ICHRG[NMXRES+1];
	282	int IDPDG[NMXRES+1];
	283	int IFLAV[NMXRES+1];
	284	int NRES;
	285	int VTOCDK[NMXRES+1];
	286	int VTORDK[NMXRES+1];
	287	int QORQQB[NMXRES+1];
	288	int QBORQQ[NMXRES+1];
	289	};
	290
	291	struct Hwunam_t {
	292	char RNAME[NMXRES+1][8];
	293	char TXNAME[NMXRES+1][2][37];
	294	};
	295
	296	int const NMXDKS = 4000;
	297	int const NMXMOD = 200;
	298
	299	struct Hwupdt_t {
	300	double BRFRAC[NMXDKS];
	301	double CMMOM[NMXDKS];
	302	double DKLTM[NMXRES];
	303	int IDK[NMXDKS];
	304	int IDKPRD[NMXDKS][5];
	305	int LNEXT[NMXDKS];
	306	int LSTRT[NMXRES];
	307	int NDKYS;
	308	int NME[NMXDKS];
	309	int NMODES[NMXRES];
	310	int NPRODS[NMXDKS];
	311	int DKPSET;
	312	int RSTAB[NMXRES+1];
	313	};
	314
	315	struct Hwuwts_t {
	316	double REPWT[5][4][4];
	317	double SNGWT;
	318	double DECWT;
	319	double QWT[3];
	320	double PWT[12];
	321	double SWTEF[NMXRES];
	322	};
	323
	324	int const NMXCDK = 4000;
	325
	326	struct Hwuclu_t {
	327	double CLDKWT[NMXCDK];
	328	double CTHRPW[12][12];
	329	double PRECO;
	330	double RESN[12][12];
	331	double RMIN[12][12];
	332	int LOCN[12][12];
	333	int NCLDK[NMXCDK];
	334	int NRECO;
	335	int CLRECO;
	336	};
	337
	338	struct Hwdist_t {
	339	double EXAG;
	340	double GEV2MM;
	341	double HBAR;
	342	double PLTCUT;
	343	double VMIN2;
	344	double VTXPIP[5];
	345	double XMIX[2];
	346	double XMRCT[2];
	347	double YMIX[2];
	348	double YMRCT[2];
	349	int IOPDKL;
	350	int MAXDKL;
	351	int MIXING;
	352	int PIPSMR;
	353	};
	354
	355	int const NMXQDK=20;
	356
	357	struct Hwqdks_t {
	358	double VTXQDK[NMXQDK][4];
	359	int IMQDK[NMXQDK];
	360	int LOCQ[NMXQDK];
	361	int NQDK;
	362	};
	363
	364	int const NMXSUD = 1024;
	365
	366	struct Hwusud_t {
	367	double ACCUR;
	368	double QEV[6][NMXSUD];
	369	double SUD[6][NMXSUD];
	370	int INTER;
	371	int NQEV;
	372	int NSUD;
	373	int SUDORD;
	374	};
	375
	376	struct Hwsusy_t {
	377	double TANB;
	378	double ALPHAH;
	379	double COSBPA;
	380	double SINBPA;
	381	double COSBMA;
	382	double SINBMA;
	383	double COSA;
	384	double SINA;
	385	double COSB;
	386	double SINB;
	387	double COTB;
	388	double ZMIXSS[4][4];
	389	double ZMXNSS[4][4];
	390	double ZSGNSS[4];
	391	double LFCH[16];
	392	double RFCH[16];
	393	double SLFCH[4][16];
	394	double SRFCH[4][16];
	395	double WMXUSS[2][2];
	396	double WMXVSS[2][2];
	397	double WSGNSS[2];
	398	double QMIXSS[2][2][6];
	399	double LMIXSS[2][2][6];
	400	double THETAT;
	401	double THETAB;
	402	double THETAL;
	403	double ATSS;
	404	double ABSS;
	405	double ALSS;
	406	double MUSS;
	407	double FACTSS;
	408	double GHWWSS[3];
	409	double GHZZSS[3];
	410	double GHDDSS[4];
	411	double GHUUSS[4];
	412	double GHWHSS[3];
	413	double GHSQSS[2][2][6][4];
	414	double XLMNSS;
	415	double RMMNSS;
	416	double DMSSM;
	417	double SENHNC[24];
	418	double SSPARITY;
	419	int SUSYIN;
	420	};
	421
	422	struct Hwrpar_t {
	423	double LAMDA1[3][3][3];
	424	double LAMDA2[3][3][3];
	425	double LAMDA3[3][3][3];
	426	int HRDCOL[5][2];
	427	int RPARTY;
	428	int COLUPD;
	429	};
	430
	431	struct Hwminb_t {
	432	double PMBN1;
	433	double PMBN2;
	434	double PMBN3;
	435	double PMBK1;
	436	double PMBK2;
	437	double PMBM1;
	438	double PMBM2;
	439	double PMBP1;
	440	double PMBP2;
	441	double PMBP3;
	442	};
	443
	444	int const NMXCL = 500;
	445
	446	struct Hwclus_t {
	447	double PPCL[NMXCL][5];
	448	int IDCL[NMXCL];
	449	int NCL;
	450	};
	451
	452	// herwig 6.507
	453
	454	struct Hwgrav_t {
	455	double GRVLAM;
	456	double EMGRV;
	457	double GAMGRV;
	458	};
	459
	460	struct Hw6202_t {
	461	double VIPWID[3];
	462	double DXRCYL;
	463	double DXZMAX;
	464	double DXRSPH;
	465	int WZRFR;
	466	int FIX4JT;
	467	int IMSSM;
	468	int IHIGGS;
	469	int PARITY;
	470	int LRSUSY;
	471	};
	472
	473	struct Hw6203_t {
	474	double ABWGT;
	475	double ABWSUM;
	476	double AVABW;
	477	int NNEGWT;
	478	int NNEGEV;
	479	int NEGWTS;
	480	};
	481
	482	int const IMAXCH = 20;
	483
	484	struct Hw6300_t {
	485	double MJJMIN;
	486	double CHNPRB[IMAXCH];
	487	int IOPSTP;
	488	int IOPSH;
	489	int OPTM;
	490	int CHON[IMAXCH];
	491	};
	492
	493
	494	int const NXMRS = 49;
	495	int const NQMRS = 37;
	496	int const NPMRS = 8;
	497
	498	struct Hwpmrs_t {
	499	double FMRS[NQMRS+1][NXMRS][NPMRS][3];
	500	};
	501
	502	struct Hwcirc_t {
	503	int CIRCOP;
	504	int CIRCAC;
	505	int CIRCVR;
	506	int CIRCRV;
	507	int CIRCCH;
	508	};
	509
	510	int const NCFMAX = 3;
	511	int const NMODE2 = 500;
	512	int const NMODE3 = 500;
	513	int const NDIAGR = 8;
	514	int const NMODEB = 50;
	515	int const NMODE4 = 4;
	516
	517	struct Hwdspb_t {
	518	double ABMODE[NMODEB][2];
	519	double BBMODE[NMODEB][12][2];
	520	double PBMODE[NMODEB][12];
	521	double WTBMAX[NMODEB][12];
	522	int IDBPRT[NMODEB];
	523	int IBDRTP[NMODEB];
	524	int IBMODE[NMODEB];
	525	int NBMODE;
	526	};
	527
	528	struct Hwdsp2_t {
	529	double A2MODE[NMODE2][2];
	530	double P2MODE[NMODE2];
	531	double WT2MAX[NMODE2];
	532	int ID2PRT[NMODE2];
	533	int I2DRTP[NMODE2];
	534	int N2MODE;
	535	};
	536
	537	struct Hwdsp3_t {
	538	double A3MODE[NMODE3][NDIAGR][2];
	539	double B3MODE[NMODE3][NDIAGR][2];
	540	double P3MODE[NMODE3];
	541	double WT3MAX[NMODE3];
	542	double SPN3CF[NMODE3][NCFMAX][NCFMAX];
	543	int ID3PRT[NMODE3];
	544	int I3MODE[NDIAGR][NMODE3];
	545	int I3DRTP[NDIAGR][NMODE3];
	546	int N3MODE;
	547	int NDI3BY[NMODE3];
	548	int N3NCFL[NMODE3];
	549	int I3DRCF[NDIAGR][NMODE3];
	550	};
	551
	552	struct Hwdsp4_t {
	553	double A4MODE[NMODE4][12][2];
	554	double B4MODE[NMODE4][12][2];
	555	double P4MODE[NMODE4][12][12];
	556	double WT4MAX[NMODE4][12][12];
	557	int ID4PRT[NMODE4];
	558	int I4MODE[NMODE4][2];
	559	int N4MODE;
	560	};
	561
	562	struct Hwdspn_t {
	563	int NDECSY;
	564	int NSEARCH;
	565	int LRDEC;
	566	int LWDEC;
	567	int SYSPIN;
	568	int THREEB;
	569	int FOURB;
	570	char TAUDEC[6];
	571	};
	572
	573	int const NMXSPN = 50;
	574
	575	struct Hwspin_t {
	576	dcpx MESPN[NMXSPN][NCFMAX][2][2][2][2];
	577	dcpx RHOSPN[NMXSPN][2][2];
	578	double SPNCFC[NMXSPN][NCFMAX][NCFMAX];
	579	int IDSPN[NMXSPN];
	580	int JMOSPN[NMXSPN];
	581	int JDASPN[NMXSPN][2];
	582	int NSPN;
	583	int ISNHEP[NMXHEP];
	584	int NSNTRY;
	585	int DECSPN[NMXSPN];
	586	int NCFL[NMXSPN];
	587	int SPCOPT;
	588	};
	589
	590	struct Hwstau_t {
	591	int JAK1;
	592	int JAK2;
	593	int ITDKRC;
	594	int IFPHOT;
	595	};
	596
	597	int const MAXHRP = 100;
	598
	599	struct Hwgupr_t {
	600	double LHWGT[MAXHRP];
	601	double LHWGTS[MAXHRP];
	602	double LHXSCT[MAXHRP];
	603	double LHXERR[MAXHRP];
	604	double LHXMAX[MAXHRP];
	605	double LHMXSM;
	606	int LHIWGT[MAXHRP];
	607	int LHNEVT[MAXHRP];
	608	int ITYPLH;
	609	int LHSOFT;
	610	int LHGLSF;
	611	};
	612
	613	struct Hw6500_t {
	614	int PRESPL;
	615	};
	616
	617	struct Hw6504_t {
	618	int ITOPRD;
	619	};
	620
	621	struct Hw6506_t {
	622	double PDFX0;
	623	double PDFPOW;
	624	};
	625
	626
	627
	628	extern "C" {
	629	void hwigin_();
	630	void hwuinc_();
	631	void hwusta_(char * name, int);
	632	void hweini_();
	633	void hwuine_();
	634	void hwepro_();
	635	void hwbgen_();
	636	void hwdhob_();
	637	void hwcfor_();
	638	void hwcdec_();
	639	void hwdhad_();
	640	void hwdhvy_();
	641	void hwmevt_();
	642	void hwufne_();
	643	void hwefin_();
	644	}
	645
	646	// subroutines to be call by JIMMY
	647
	648	extern "C" {
	649	void jminit_();
	650	void jimmin_();
	651	void jmefin_();
	652	}
	653
	654
	655	#endif
	656
	657
	658