ISAJET/code/mginit.F

   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