ISAJET/code/xwwzz.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE XWWZZ
   3 C
   4 C          SET UP W+ W- -> Z0 Z0 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.)
  52       QQF=.5*SQRT(S-4.*ZM2)
  53       GSQ=4.*PI*ALFA/SW2
  54 C
  55 C          FROM WWZZ3.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 * ZM2
  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      $ * ZM3) / (CW * PROPH)) + 8.00E+00 * ((HM2 * S * ZM3) / (CW
  68      $ * PROPH)) - 3.20E+01 * CW2 * ZM2 + 1.20E+01 * CW2 * S * 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) = 6.40E+01 * CW2 * ZM2 + 1.20E+01 * CW2 * S * ZM2
  75      $ - 4.00E+00 * CW2 * S * ZM4 - 2.00E+00 * CW2 * S2 * ZM2
  76       ANWWWW(1,2,3) = -6.40E+01 * CW2 * ZM2 + 2.00E+01 * CW2 * S * ZM2
  77      $ + 4.00E+00 * CW2 * S * ZM4 - 6.00E+00 * CW2 * S2 * ZM2
  78       ANWWWW(1,2,4) = 1.92E+02 * CW2 * QQ0 * QQF * QQI * ZM
  79      $ - 3.20E+01 * CW2 * QQ0 * QQF * QQI * ZM3 + 3.20E+01 * CW2
  80      $ * QQ0 * QQF * QQI * ZM5 + 1.60E+01 * CW2 * QQ0 * QQF * QQI * S
  81      $ * ZM
  82       ANWWWW(1,3,1) = -4.00E+00 * CW2 * S2 + 3.00E+00 * CW2 * S3
  83      $ - 9.60E+01 * CW2 * ZM4 + 3.20E+01 * CW2 * ZM6 + 8.00E+00 * CW2
  84      $ * S * ZM2 + 1.60E+01 * CW2 * S * ZM4 - 8.00E+00 * CW2 * S
  85      $ * ZM6 - 1.00E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4
  86       ANWWWW(1,3,2) = 6.40E+01 * CW2 * ZM2 + 1.20E+01 * CW2 * S * ZM2
  87      $ - 4.00E+00 * CW2 * S * ZM4 - 2.00E+00 * CW2 * S2 * ZM2
  88       ANWWWW(1,3,3) = -6.40E+01 * CW2 * ZM2 + 2.00E+01 * CW2 * S * ZM2
  89      $ + 4.00E+00 * CW2 * S * ZM4 - 6.00E+00 * CW2 * S2 * ZM2
  90       ANWWWW(1,3,4) = -1.92E+02 * CW2 * QQ0 * QQF * QQI * ZM
  91      $ + 3.20E+01 * CW2 * QQ0 * QQF * QQI * ZM3 - 3.20E+01 * CW2
  92      $ * QQ0 * QQF * QQI * ZM5 - 1.60E+01 * CW2 * QQ0 * QQF * QQI * S
  93      $ * ZM
  94       ANWWWW(1,4,1) = 0.00E+00
  95       ANWWWW(1,4,2) = 0.00E+00
  96       ANWWWW(1,4,3) = 0.00E+00
  97       ANWWWW(1,4,4) = 0.00E+00
  98       ANWWWW(2,1,1) = 0.00E+00
  99       ANWWWW(2,1,2) = 0.00E+00
 100       ANWWWW(2,1,3) = 0.00E+00
 101       ANWWWW(2,1,4) = 8.00E+00 * CW2 * QQ0 * S * ZM
 102       ANWWWW(2,2,1) = 4.80E+01 * CW2 * QQF * QQI * S - 2.00E+01 * CW2
 103      $ * QQF * QQI * S2 + 6.40E+01 * CW2 * QQF * QQI * ZM2 + 9.60E+01
 104      $ * CW2 * QQF * QQI * S * ZM2 - 1.60E+01 * CW2 * QQF * QQI * S
 105      $ * ZM4
 106       ANWWWW(2,2,2) = 8.00E+00 * CW2 * QQF * QQI * S * ZM2
 107       ANWWWW(2,2,3) = 5.60E+01 * CW2 * QQF * QQI * S * ZM2
 108       ANWWWW(2,2,4) = 1.28E+02 * CW2 * QQ0 * ZM + 1.92E+02 * CW2 * QQ0
 109      $ * ZM3 - 3.20E+01 * CW2 * QQ0 * S * ZM + 2.40E+01 * CW2 * QQ0
 110      $ * S * ZM3 - 8.00E+00 * CW2 * QQ0 * S * ZM5 - 8.00E+00 * CW2
 111      $ * QQ0 * S2 * ZM
 112       ANWWWW(2,3,1) = -4.80E+01 * CW2 * QQF * QQI * S + 2.00E+01 * CW2
 113      $ * QQF * QQI * S2 - 6.40E+01 * CW2 * QQF * QQI * ZM2
 114      $ - 9.60E+01 * CW2 * QQF * QQI * S * ZM2 + 1.60E+01 * CW2 * QQF
 115      $ * QQI * S * ZM4
 116       ANWWWW(2,3,2) = -8.00E+00 * CW2 * QQF * QQI * S * ZM2
 117       ANWWWW(2,3,3) = -5.60E+01 * CW2 * QQF * QQI * S * ZM2
 118       ANWWWW(2,3,4) = 1.28E+02 * CW2 * QQ0 * ZM + 1.92E+02 * CW2 * QQ0
 119      $ * ZM3 - 3.20E+01 * CW2 * QQ0 * S * ZM + 2.40E+01 * CW2 * QQ0
 120      $ * S * ZM3 - 8.00E+00 * CW2 * QQ0 * S * ZM5 - 8.00E+00 * CW2
 121      $ * QQ0 * S2 * ZM
 122       ANWWWW(2,4,1) = 0.00E+00
 123       ANWWWW(2,4,2) = 0.00E+00
 124       ANWWWW(2,4,3) = 0.00E+00
 125       ANWWWW(2,4,4) = 0.00E+00
 126       ANWWWW(3,1,1) = 2.00E+00 * CW2 * S2
 127       ANWWWW(3,1,2) = -4.00E+00 * CW2 * S * ZM2
 128       ANWWWW(3,1,3) = 4.00E+00 * CW2 * S * ZM2
 129       ANWWWW(3,1,4) = 0.00E+00
 130       ANWWWW(3,2,1) = 3.20E+01 * CW2 * S - 1.60E+01 * CW2 * S2 + CW2
 131      $ * S3 + 9.60E+01 * CW2 * S * ZM2 + 3.20E+01 * CW2 * S * ZM4
 132      $ - 2.20E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4
 133       ANWWWW(3,2,2) = -6.40E+01 * CW2 * ZM2 - 1.20E+01 * CW2 * S
 134      $ * ZM2 + 4.00E+00 * CW2 * S * ZM4 + 2.00E+00 * CW2 * S2 * ZM2
 135       ANWWWW(3,2,3) = 6.40E+01 * CW2 * ZM2 + 1.20E+01 * CW2 * S * ZM2
 136      $ + 2.80E+01 * CW2 * S * ZM4 - 1.00E+01 * CW2 * S2 * ZM2
 137       ANWWWW(3,2,4) = 1.60E+01 * CW2 * QQ0 * QQF * QQI * S * ZM
 138       ANWWWW(3,3,1) = 3.20E+01 * CW2 * S - 1.60E+01 * CW2 * S2 + CW2
 139      $ * S3 + 9.60E+01 * CW2 * S * ZM2 + 3.20E+01 * CW2 * S * ZM4
 140      $ - 2.20E+01 * CW2 * S2 * ZM2 + 2.00E+00 * CW2 * S2 * ZM4
 141       ANWWWW(3,3,2) = -6.40E+01 * CW2 * ZM2 - 1.20E+01 * CW2 * S
 142      $ * ZM2 + 4.00E+00 * CW2 * S * ZM4 + 2.00E+00 * CW2 * S2 * ZM2
 143       ANWWWW(3,3,3) = 6.40E+01 * CW2 * ZM2 + 1.20E+01 * CW2 * S * ZM2
 144      $ + 2.80E+01 * CW2 * S * ZM4 - 1.00E+01 * CW2 * S2 * ZM2
 145       ANWWWW(3,3,4) = -1.60E+01 * CW2 * QQ0 * QQF * QQI * S * ZM
 146       ANWWWW(3,4,1) = 0.00E+00
 147       ANWWWW(3,4,2) = 0.00E+00
 148       ANWWWW(3,4,3) = 0.00E+00
 149       ANWWWW(3,4,4) = 0.00E+00
 150       ANWWWW(4,1,1) = 0.00E+00
 151       ANWWWW(4,1,2) = 0.00E+00
 152       ANWWWW(4,1,3) = 0.00E+00
 153       ANWWWW(4,1,4) = 0.00E+00
 154       ANWWWW(4,2,1) = 4.00E+00 * CW2 * QQF * QQI * S2
 155       ANWWWW(4,2,2) = -8.00E+00 * CW2 * QQF * QQI * S * ZM2
 156       ANWWWW(4,2,3) = 8.00E+00 * CW2 * QQF * QQI * S * ZM2
 157       ANWWWW(4,2,4) = 0.00E+00
 158       ANWWWW(4,3,1) = -4.00E+00 * CW2 * QQF * QQI * S2
 159       ANWWWW(4,3,2) = 8.00E+00 * CW2 * QQF * QQI * S * ZM2
 160       ANWWWW(4,3,3) = -8.00E+00 * CW2 * QQF * QQI * S * ZM2
 161       ANWWWW(4,3,4) = 0.00E+00
 162       ANWWWW(4,4,1) = 0.00E+00
 163       ANWWWW(4,4,2) = 0.00E+00
 164       ANWWWW(4,4,3) = 0.00E+00
 165       ANWWWW(4,4,4) = 0.00E+00
 166 C
 167       ADWWWW(1,1) = 1.00E+00
 168       ADWWWW(1,2) = S - 2.00E+00 * ZM2
 169       ADWWWW(1,3) = S - 2.00E+00 * ZM2
 170       ADWWWW(1,4) = 1.00E+00
 171       ADWWWW(2,1) = 0.00E+00
 172       ADWWWW(2,2) = -1.00E+00 * S + 4.00E+00 * (S / (S + 4.00E+00
 173      $ * QQF * QQI)) - 1.60E+01 * (ZM2 / (S + 4.00E+00 * QQF * QQI))
 174      $ + 4.00E+00 * ((S * ZM2) / (S + 4.00E+00 * QQF * QQI))
 175       ADWWWW(2,3) = S - 4.00E+00 * (S / (S + 4.00E+00 * QQF * QQI))
 176      $ + 1.60E+01 * (ZM2 / (S + 4.00E+00 * QQF * QQI)) - 4.00E+00
 177      $ * ((S * ZM2) / (S + 4.00E+00 * QQF * QQI))
 178       ADWWWW(2,4) = 0.00E+00
 179 C
 180       AIWWWW(1) = 1.60E+01 * ((HG * HM * ZM3) / (CW * PROPH))
 181      $ - 8.00E+00 * ((HG * HM * S * ZM) / (CW * PROPH)) - 8.00E+00
 182      $ * ((HG * HM * S * ZM3) / (CW * PROPH)) + 4.00E+00 * ((HG * HM
 183      $ * S2 * ZM) / (CW * PROPH))
 184       AIWWWW(2) = 0.00E+00
 185       AIWWWW(3) = 1.60E+01 * ((HG * HM * ZM3) / (CW * PROPH))
 186      $ - 8.00E+00 * ((HG * HM * S * ZM3) / (CW * PROPH))
 187       AIWWWW(4) = 0.00E+00
 188 C
 189 C          RESTORE MISSING FACTORS
 190       DO 100 J=1,4
 191       AIWWWW(J)=AIWWWW(J)*GSQ/(16.*ZM2)
 192       DO 100 I=1,4
 193       DO 110 K=1,4
 194 110   ANWWWW(K,I,J)=ANWWWW(K,I,J)*GSQ/(16.*ZM2)
 195 100   CONTINUE
 196 C
 197       RETURN
 198       END