More volume overlaps corrected
[u/mrichter/AliRoot.git] / ISAJET / code / mginit.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 SUBROUTINE MGINIT
3C
4C Initialize common blocks for MadGraph code in ISAJET
5C Note the QCD coupling constant is g=1.
6C
7#if defined(CERNLIB_IMPNONE)
8 IMPLICIT NONE
9#endif
10C
11#include "isajet/itapes.inc"
12#include "isajet/sstype.inc"
13#include "isajet/mgcoms.inc"
14C
15 INTEGER I
16 REAL AMGMW
17 REAL*8 SW2
18C
19C Fermion masses and widths
20C
21 FMASS(1) = AMGMW(IDE,1)
22 FMASS(2) = 0D0
23 FMASS(3) = AMGMW(IDUP,1)
24 FMASS(4) = AMGMW(IDDN,1)
25 FMASS(5) = AMGMW(IDMU,1)
26 FMASS(6) = 0D0
27 FMASS(7) = AMGMW(IDCH,1)
28 FMASS(8) = AMGMW(IDST,1)
29 FMASS(9) = AMGMW(IDTAU,1)
30 FMASS(10)= 0D0
31 FMASS(11)= AMGMW(IDTP,1)
32 FMASS(12)= AMGMW(IDBT,1)
33 DO 100 I=1,12
34 FWIDTH(I)=0D0
35100 CONTINUE
36C
37C Boson masses and widths
38C
39 AMASS=0D0
40 AWIDTH=0D0
41 WMASS=AMGMW(IDW,1)
42 WWIDTH=AMGMW(IDW,2)
43 ZMASS=AMGMW(IDZ,1)
44 ZWIDTH=AMGMW(IDZ,2)
45 HMASS=AMGMW(IDH,1)
46 HWIDTH=AMGMW(IDH,2)
47 SW2=AMGMW(1,3)
48C
49C Calls to Helas routines to set couplings
50C
51 CALL COUP1X(SW2,GW,GWWA,GWWZ)
52 CALL COUP2X(SW2,GAL,GAU,GAD,GWF,GZN,GZL,GZU,GZD,G1)
53 CALL COUP3X(SW2,ZMASS,HMASS,GWWH,GZZH,GHHH,GWWHH,GZZHH,GHHHH)
54 DO 110 I=1,12
55 CALL COUP4X(SW2,ZMASS,FMASS(I),GCHF(1,I))
56110 CONTINUE
57C
58C QCD couplings
59C
60 G = 1D0
61 GG(1)=-G
62 GG(2)=-G
63 RETURN
64 END