ISAJET/code/xwwww.F

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