ISAJET/code/xzzww.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE XZZWW
   3 C
   4 C          SET UP Z0 Z0 -> W+ W- AMPLITUDES AS RATIONAL FUNCTIONS OF Z
   5 C
   6 C          RE(F(Z,L)) = SUM(I,J)(ANWWWW(I+1,J,L)*Z**I
   7 C                                  /(ADWWWW(1,J)+ADWWWW(2,J)*Z))
   8 C          IM(F(Z,L)) = AIWWWW(L)   (INDEPENDENT OF Z)
   9 C          J LABELS PIECES WITH SAME DENOMINATOR.
  10 C          L=1 FOR 0,0; L=2 FOR 1,-1; L=3 FOR 1,1; L=4 FOR 0,1
  11 C
  12 C          *NOTE* A FACTOR OF SIN(THETA)/SQRT(2) IS REMOVED FROM F01
  13 C
  14 #include "isajet/const.inc"
  15 #include "isajet/jetpar.inc"
  16 #include "isajet/wcon.inc"
  17 #include "isajet/hcon.inc"
  18 #if defined(CERNLIB_DOUBLE)
  19       DOUBLE PRECISION WM,ZM,ZM2,ZM3,ZM4,ZM5,ZM6,HM,HM2,HM3,HM4,HG,HG2
  20      $,PROPH,RTS,S,S2,S3,SW,QQ0,QQI,QQF
  21 #endif
  22 C
  23 C          USE UNITS OF WM TO AVOID LARGE NUMBERS - NOTE ANWWWW/ADWWWW
  24 C          AND AIWWWW ARE DIMENSIONLESS
  25       WM=WMASS(2)
  26       ZM=WMASS(4)/WM
  27       ZM2=ZM**2
  28       ZM3=ZM**3
  29       ZM4=ZM**4
  30       ZM5=ZM**5
  31       ZM6=ZM**6
  32       HM=HMASS/WM
  33       HM2=HM**2
  34       HM3=HM**3
  35       HM4=HM**4
  36       HG=HGAM/WM
  37       HG2=HG**2
  38       RTS=QMW/WM
  39       S=RTS**2
  40       S2=S**2
  41       S3=S**3
  42       PROPH=(S-HM2)**2+(HM*HG)**2
  43 C          CORRECT SIGN OF HIGGS AMPLITUDE.
  44       PROPH=-PROPH
  45 C
  46       CW=1./ZM
  47       CW2=CW**2
  48       SW2=1.-CW2
  49       SW=SQRT(SW2)
  50       QQ0=.5*RTS
  51       QQI=.5*SQRT(S-4.*ZM2)
  52       QQF=.5*SQRT(S-4.)
  53       GSQ=4.*PI*ALFA/SW2
  54 C
  55 C          FROM ZZWW3.EX
  56       ANWWWW(1,1,1) = -1.60E+01 * ((HM2 * ZM3) / (CW * PROPH))
  57      $ + 1.60E+01 * ((S * ZM3) / (CW * PROPH) ) - 8.00E+00 * ((S2
  58      $ * ZM) / (CW * PROPH)) - 8.00E+00 * ((S2 * ZM3) / (CW * PROPH))
  59      $ + 4.00E+00 * ((S3 * ZM) / (CW * PROPH)) + 8.00E+00 * ((HM2 * S
  60      $ * ZM) / (CW * PROPH)) + 8.00E+00 * ((HM2 * S * ZM3) / (CW
  61      $ * PROPH)) - 4.00E+00 * ((HM2 * S2 * ZM) / (CW * PROPH))
  62      $ + 8.00E+00 * CW2 * S - 6.00E+00 * CW2 * S2 + 8.00E+00 * CW2
  63      $ * S * ZM2
  64       ANWWWW(1,1,2) = 4.00E+00 * CW2 * S
  65       ANWWWW(1,1,3) = -1.60E+01 * ((HM2 * ZM3) / (CW * PROPH))
  66      $ + 1.60E+01 * ((S * ZM3) / (CW * PROPH) ) - 8.00E+00 * ((S2
  67      $ * ZM) / (CW * PROPH)) + 8.00E+00 * ((HM2 * S * ZM) / (CW
  68      $ * PROPH)) + 1.20E+01 * CW2 * S - 3.20E+01 * CW2 * ZM2
  69       ANWWWW(1,1,4) = 0.00E+00
  70       ANWWWW(1,2,1) = -4.00E+00 * CW2 * S2 + 3.00E+00 * CW2 * S3
  71      $ - 9.60E+01 * CW2 * ZM4 + 3.20E+01 * CW2 * ZM6 + 8.00E+00 * CW2
  72      $ * S * ZM2 + 1.60E+01 * CW2 * S * ZM4 - 8.00E+00 * CW2 * S
  73      $ * ZM6 - 1.00E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4
  74       ANWWWW(1,2,2) = -4.00E+00 * CW2 * S - 2.00E+00 * CW2 * S2
  75      $ + 6.40E+01 * CW2 * ZM4 + 1.20E+01 * CW2 * S * ZM2
  76       ANWWWW(1,2,3) = 4.00E+00 * CW2 * S - 6.00E+00 * CW2 * S2
  77      $ - 6.40E+01 * CW2 * ZM4 + 2.00E+01 * CW2 * S * ZM2
  78       ANWWWW(1,2,4) = 6.40E+01 * CW2 * QQ0 * QQF * QQI + 1.60E+01
  79      $ * CW2 * QQ0 * QQF * QQI * S + 9.60E+01 * CW2 * QQ0 * QQF * QQI
  80      $ * ZM2 + 3.20E+01 * CW2 * QQ0 * QQF * QQI * ZM4
  81       ANWWWW(1,3,1) = -4.00E+00 * CW2 * S2 + 3.00E+00 * CW2 * S3
  82      $ - 9.60E+01 * CW2 * ZM4 + 3.20E+01 * CW2 * ZM6 + 8.00E+00 * CW2
  83      $ * S * ZM2 + 1.60E+01 * CW2 * S * ZM4 - 8.00E+00 * CW2 * S
  84      $ * ZM6 - 1.00E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4
  85       ANWWWW(1,3,2) = -4.00E+00 * CW2 * S - 2.00E+00 * CW2 * S2
  86      $ + 6.40E+01 * CW2 * ZM4 + 1.20E+01 * CW2 * S * ZM2
  87       ANWWWW(1,3,3) = 4.00E+00 * CW2 * S - 6.00E+00 * CW2 * S2
  88      $ - 6.40E+01 * CW2 * ZM4 + 2.00E+01 * CW2 * S * ZM2
  89       ANWWWW(1,3,4) = -6.40E+01 * CW2 * QQ0 * QQF * QQI - 1.60E+01
  90      $ * CW2 * QQ0 * QQF * QQI * S - 9.60E+01 * CW2 * QQ0 * QQF * QQI
  91      $ * ZM2 - 3.20E+01 * CW2 * QQ0 * QQF * QQI * ZM4
  92       ANWWWW(1,4,1) = 0.00E+00
  93       ANWWWW(1,4,2) = 0.00E+00
  94       ANWWWW(1,4,3) = 0.00E+00
  95       ANWWWW(1,4,4) = 0.00E+00
  96       ANWWWW(2,1,1) = 0.00E+00
  97       ANWWWW(2,1,2) = 0.00E+00
  98       ANWWWW(2,1,3) = 0.00E+00
  99       ANWWWW(2,1,4) = 8.00E+00 * CW2 * QQ0 * S
 100       ANWWWW(2,2,1) = 4.80E+01 * CW2 * QQF * QQI * S - 2.00E+01 * CW2
 101      $ * QQF * QQI * S2 + 6.40E+01 * CW2 * QQF * QQI * ZM2 + 9.60E+01
 102      $ * CW2 * QQF * QQI * S * ZM2 - 1.60E+01 * CW2 * QQF * QQI * S
 103      $ * ZM4
 104       ANWWWW(2,2,2) = 8.00E+00 * CW2 * QQF * QQI * S
 105       ANWWWW(2,2,3) = 5.60E+01 * CW2 * QQF * QQI * S
 106       ANWWWW(2,2,4) = -8.00E+00 * CW2 * QQ0 * S2 + 1.92E+02 * CW2
 107      $ * QQ0 * ZM2 + 1.28E+02 * CW2 * QQ0 * ZM4 - 8.00E+00 * CW2
 108      $ * QQ0 * S * ZM2 - 8.00E+00 * CW2 * QQ0 * S * ZM4
 109       ANWWWW(2,3,1) = -4.80E+01 * CW2 * QQF * QQI * S + 2.00E+01 * CW2
 110      $ * QQF * QQI * S2 - 6.40E+01 * CW2 * QQF * QQI * ZM2
 111      $ - 9.60E+01 * CW2 * QQF * QQI * S * ZM2 + 1.60E+01 * CW2 * QQF
 112      $ * QQI * S * ZM4
 113       ANWWWW(2,3,2) = -8.00E+00 * CW2 * QQF * QQI * S
 114       ANWWWW(2,3,3) = -5.60E+01 * CW2 * QQF * QQI * S
 115       ANWWWW(2,3,4) = -8.00E+00 * CW2 * QQ0 * S2 + 1.92E+02 * CW2
 116      $ * QQ0 * ZM2 + 1.28E+02 * CW2 * QQ0 * ZM4 - 8.00E+00 * CW2
 117      $ * QQ0 * S * ZM2 - 8.00E+00 * CW2 * QQ0 * S * ZM4
 118       ANWWWW(2,4,1) = 0.00E+00
 119       ANWWWW(2,4,2) = 0.00E+00
 120       ANWWWW(2,4,3) = 0.00E+00
 121       ANWWWW(2,4,4) = 0.00E+00
 122       ANWWWW(3,1,1) = 2.00E+00 * CW2 * S2
 123       ANWWWW(3,1,2) = -4.00E+00 * CW2 * S
 124       ANWWWW(3,1,3) = 4.00E+00 * CW2 * S
 125       ANWWWW(3,1,4) = 0.00E+00
 126       ANWWWW(3,2,1) = 3.20E+01 * CW2 * S - 1.60E+01 * CW2 * S2 + CW2
 127      $ * S3 + 9.60E+01 * CW2 * S * ZM2 + 3.20E+01 * CW2 * S * ZM4
 128      $ - 2.20E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4
 129       ANWWWW(3,2,2) = 4.00E+00 * CW2 * S + 2.00E+00 * CW2 * S2
 130      $ - 6.40E+01 * CW2 * ZM4 - 1.20E+01 * CW2 * S * ZM2
 131       ANWWWW(3,2,3) = 2.80E+01 * CW2 * S - 1.00E+01 * CW2 * S2
 132      $ + 6.40E+01 * CW2 * ZM4 + 1.20E+01 * CW2 * S * ZM2
 133       ANWWWW(3,2,4) = 1.60E+01 * CW2 * QQ0 * QQF * QQI * S
 134       ANWWWW(3,3,1) = 3.20E+01 * CW2 * S - 1.60E+01 * CW2 * S2 + CW2
 135      $ * S3 + 9.60E+01 * CW2 * S * ZM2 + 3.20E+01 * CW2 * S * ZM4
 136      $ - 2.20E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4
 137       ANWWWW(3,3,2) = 4.00E+00 * CW2 * S + 2.00E+00 * CW2 * S2
 138      $ - 6.40E+01 * CW2 * ZM4 - 1.20E+01 * CW2 * S * ZM2
 139       ANWWWW(3,3,3) = 2.80E+01 * CW2 * S - 1.00E+01 * CW2 * S2
 140      $ + 6.40E+01 * CW2 * ZM4 + 1.20E+01 * CW2 * S * ZM2
 141       ANWWWW(3,3,4) = -1.60E+01 * CW2 * QQ0 * QQF * QQI * S
 142       ANWWWW(3,4,1) = 0.00E+00
 143       ANWWWW(3,4,2) = 0.00E+00
 144       ANWWWW(3,4,3) = 0.00E+00
 145       ANWWWW(3,4,4) = 0.00E+00
 146       ANWWWW(4,1,1) = 0.00E+00
 147       ANWWWW(4,1,2) = 0.00E+00
 148       ANWWWW(4,1,3) = 0.00E+00
 149       ANWWWW(4,1,4) = 0.00E+00
 150       ANWWWW(4,2,1) = 4.00E+00 * CW2 * QQF * QQI * S2
 151       ANWWWW(4,2,2) = -8.00E+00 * CW2 * QQF * QQI * S
 152       ANWWWW(4,2,3) = 8.00E+00 * CW2 * QQF * QQI * S
 153       ANWWWW(4,2,4) = 0.00E+00
 154       ANWWWW(4,3,1) = -4.00E+00 * CW2 * QQF * QQI * S2
 155       ANWWWW(4,3,2) = 8.00E+00 * CW2 * QQF * QQI * S
 156       ANWWWW(4,3,3) = -8.00E+00 * CW2 * QQF * QQI * S
 157       ANWWWW(4,3,4) = 0.00E+00
 158       ANWWWW(4,4,1) = 0.00E+00
 159       ANWWWW(4,4,2) = 0.00E+00
 160       ANWWWW(4,4,3) = 0.00E+00
 161       ANWWWW(4,4,4) = 0.00E+00
 162 C
 163       ADWWWW(1,1) = 1.00E+00
 164       ADWWWW(1,2) = S - 2.00E+00 * ZM2
 165       ADWWWW(1,3) = S - 2.00E+00 * ZM2
 166       ADWWWW(1,4) = 1.00E+00
 167       ADWWWW(2,1) = 0.00E+00
 168       ADWWWW(2,2) = -1.00E+00 * S + 4.00E+00 * (S / (S + 4.00E+00
 169      $ * QQF * QQI)) - 1.60E+01 * (ZM2 / (S + 4.00E+00 * QQF * QQI))
 170      $ + 4.00E+00 * ((S * ZM2) / (S + 4.00E+00 * QQF * QQI))
 171       ADWWWW(2,3) = S - 4.00E+00 * (S / (S + 4.00E+00 * QQF * QQI))
 172      $ + 1.60E+01 * (ZM2 / (S + 4.00E+00 * QQF * QQI)) - 4.00E+00
 173      $ * ((S * ZM2) / (S + 4.00E+00 * QQF * QQI))
 174       ADWWWW(2,4) = 0.00E+00
 175 C
 176       AIWWWW(1) = 1.60E+01 * ((HG * HM * ZM3) / (CW * PROPH))
 177      $ - 8.00E+00 * ((HG * HM * S * ZM) / (CW * PROPH))
 178      $ - 8.00E+00 * ((HG * HM * S * ZM3) / (CW * PROPH)) + 4.00E+00
 179      $ * ((HG * HM * S2 * ZM) / (CW * PROPH))
 180       AIWWWW(2) = 0.00E+00
 181       AIWWWW(3) = 1.60E+01 * ((HG * HM * ZM3) / (CW * PROPH))
 182      $ - 8.00E+00 * ((HG * HM * S * ZM) / (CW * PROPH))
 183       AIWWWW(4) = 0.00E+00
 184 C
 185 C          RESTORE MISSING FACTORS
 186       DO 100 J=1,4
 187       AIWWWW(J)=AIWWWW(J)*GSQ/(16.*ZM2)
 188       DO 100 I=1,4
 189       DO 110 K=1,4
 190 110   ANWWWW(K,I,J)=ANWWWW(K,I,J)*GSQ/(16.*ZM2)
 191 100   CONTINUE
 192 C
 193       RETURN
 194       END