]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | SUBROUTINE MGINIT | |
3 | C | |
4 | C Initialize common blocks for MadGraph code in ISAJET | |
5 | C Note the QCD coupling constant is g=1. | |
6 | C | |
7 | #if defined(CERNLIB_IMPNONE) | |
8 | IMPLICIT NONE | |
9 | #endif | |
10 | C | |
11 | #include "isajet/itapes.inc" | |
12 | #include "isajet/sstype.inc" | |
13 | #include "isajet/mgcoms.inc" | |
14 | C | |
15 | INTEGER I | |
16 | REAL AMGMW | |
17 | REAL*8 SW2 | |
18 | C | |
19 | C Fermion masses and widths | |
20 | C | |
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 | |
35 | 100 CONTINUE | |
36 | C | |
37 | C Boson masses and widths | |
38 | C | |
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) | |
48 | C | |
49 | C Calls to Helas routines to set couplings | |
50 | C | |
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)) | |
56 | 110 CONTINUE | |
57 | C | |
58 | C QCD couplings | |
59 | C | |
60 | G = 1D0 | |
61 | GG(1)=-G | |
62 | GG(2)=-G | |
63 | RETURN | |
64 | END |