ISAJET/code/sigwh.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE SIGWH
   3 C
   4 C          Calculate d(sigma)/d(pt**2)d(y1)d(y2) for WH and ZH
   5 C          associated production.
   6 C
   7 C          SIGMA    = cross section summed over types allowed by
   8 C                     JETTYPE cards.
   9 C          SIGS(I)  = partial cross section for I1 + I2 --> I3 + I4
  10 C          INOUT(I) = IOPAK**3*I4 + IOPAK**2*I3 + IOPAK*I2 +I1
  11 C
  12 C          Extra factor of 1/2 needed for nonidentical final jets.
  13 C          Y=-log(tan(theta/2)) gives jacobean P1*P2/E1*E2
  14 C
  15 C
  16 #if defined(CERNLIB_IMPNONE)
  17       IMPLICIT NONE
  18 #endif
  19 #include "isajet/itapes.inc"
  20 #include "isajet/const.inc"
  21 #include "isajet/jetpar.inc"
  22 #include "isajet/jetsig.inc"
  23 #include "isajet/primar.inc"
  24 #include "isajet/q1q2.inc"
  25 #include "isajet/qcdpar.inc"
  26 #include "isajet/wcon.inc"
  27 #include "isajet/hcon.inc"
  28 #include "isajet/xmssm.inc"
  29 C
  30       REAL X(2)
  31       EQUIVALENCE (X(1),X1)
  32       EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT)
  33       REAL SIG,S,T,U,FAC,AMW,AMZ,AMW2,AMZ2,E1,E2
  34       REAL QFCN,STRUC,SIGHW
  35       REAL PROPZ,PROPW,GV(2),GA(2),AMH,AMH2,GAMW,GAMZ
  36       INTEGER IS2UD(25),IQ,IH,I,IQ1,IQ2,IFLQ
  37       SAVE IS2UD
  38 C
  39 C          IS2UD: Susy jettype -> u/d code
  40       DATA IS2UD/0,1,1,2,2,2,2,1,1,2,2,1,1,1,1,2,2,2,2,1,1,2,2,1,1/
  41
  42 C          Functions
  43       QFCN(IQ,IH)=STRUC(X(IH),QSQ,IQ,IDIN(IH))/X(IH)
  44 C
  45       IF (GOMSSM) THEN
  46         CALL SIGWHS
  47         RETURN
  48       END IF
  49 C          Initialize
  50       DO 10 I=1,MXSIGS
  51 10    SIGS(I)=0.
  52       SIGMA=0.
  53       NSIGS=0
  54 C
  55       AMW=WMASS(2)
  56       AMW2=AMW**2
  57       AMZ=WMASS(4)
  58       AMZ2=AMZ**2
  59       AMH=HMASS
  60       AMH2=AMH**2
  61       GAMW=WGAM(2)
  62       GAMZ=WGAM(4)
  63       GV(1)=.25-2*SIN2W/3.
  64       GV(2)=-.25+SIN2W/3.
  65       GA(1)=-.25
  66       GA(2)=.25
  67 C
  68 C          WH production via W-*
  69 C
  70       IF (GOQ(28,1).AND.GOQ(30,2)) THEN
  71           CALL TWOKIN(0.,0.,AMW,AMH)
  72           IF(X1.GE.1..OR.X2.GE.1.) GO TO 100
  73           E1=SQRT(P(1)**2+AMW**2)
  74           E2=SQRT(P(2)**2+AMH**2)
  75           FAC=1./(12.*PI*S**2)
  76           FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
  77           PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
  78           SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
  79      $          PROPW*TBRWW(3,1)
  80           SIG=.5*SIGHW*FAC*QFCN(3,1)*QFCN(4,2)
  81           CALL SIGFIL(SIG,3,4,28,30)
  82           SIG=.5*SIGHW*FAC*QFCN(4,1)*QFCN(3,2)
  83           CALL SIGFIL(SIG,4,3,28,30)
  84           SIG=.5*SIGHW*FAC*QFCN(9,1)*QFCN(6,2)
  85           CALL SIGFIL(SIG,9,6,28,30)
  86           SIG=.5*SIGHW*FAC*QFCN(6,1)*QFCN(9,2)
  87           CALL SIGFIL(SIG,6,9,28,30)
  88 100       CONTINUE
  89       END IF
  90 C
  91       IF (GOQ(30,1).AND.GOQ(28,2)) THEN
  92           CALL TWOKIN(0.,0.,AMH,AMW)
  93           IF(X1.GE.1..OR.X2.GE.1.) GO TO 110
  94           E1=SQRT(P(1)**2+AMH**2)
  95           E2=SQRT(P(2)**2+AMW**2)
  96           FAC=1./(12.*PI*S**2)
  97           FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
  98           PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
  99           SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
 100      $          PROPW*TBRWW(3,2)
 101           SIG=.5*SIGHW*FAC*QFCN(3,1)*QFCN(4,2)
 102           CALL SIGFIL(SIG,3,4,30,28)
 103           SIG=.5*SIGHW*FAC*QFCN(4,1)*QFCN(3,2)
 104           CALL SIGFIL(SIG,4,3,30,28)
 105           SIG=.5*SIGHW*FAC*QFCN(9,1)*QFCN(6,2)
 106           CALL SIGFIL(SIG,9,6,30,28)
 107           SIG=.5*SIGHW*FAC*QFCN(6,1)*QFCN(9,2)
 108           CALL SIGFIL(SIG,6,9,30,28)
 109 110       CONTINUE
 110       END IF
 111 C
 112 C
 113 C          WH production via W+*
 114 C
 115       IF (GOQ(27,1).AND.GOQ(30,2)) THEN
 116           CALL TWOKIN(0.,0.,AMW,AMH)
 117           IF(X1.GE.1..OR.X2.GE.1.) GO TO 120
 118           E1=SQRT(P(1)**2+AMW**2)
 119           E2=SQRT(P(2)**2+AMH**2)
 120           FAC=1./(12.*PI*S**2)
 121           FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
 122           PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
 123           SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
 124      $          PROPW*TBRWW(2,1)
 125           SIG=.5*SIGHW*FAC*QFCN(2,1)*QFCN(5,2)
 126           CALL SIGFIL(SIG,2,5,27,30)
 127           SIG=.5*SIGHW*FAC*QFCN(5,1)*QFCN(2,2)
 128           CALL SIGFIL(SIG,5,2,27,30)
 129           SIG=.5*SIGHW*FAC*QFCN(8,1)*QFCN(7,2)
 130           CALL SIGFIL(SIG,8,7,27,30)
 131           SIG=.5*SIGHW*FAC*QFCN(7,1)*QFCN(8,2)
 132           CALL SIGFIL(SIG,7,8,27,30)
 133 120       CONTINUE
 134       END IF
 135 C
 136       IF (GOQ(30,1).AND.GOQ(27,2)) THEN
 137           CALL TWOKIN(0.,0.,AMH,AMW)
 138           IF(X1.GE.1..OR.X2.GE.1.) GO TO 130
 139           E1=SQRT(P(1)**2+AMH**2)
 140           E2=SQRT(P(2)**2+AMW**2)
 141           FAC=1./(12.*PI*S**2)
 142           FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
 143           PROPW=(S-AMW**2)**2+AMW**2*GAMW**2
 144           SIGHW=GF**2*AMW**8*(S/AMW2+(1.-T/AMW2)*(1.-U/AMW2))/
 145      $          PROPW*TBRWW(2,2)
 146           SIG=.5*SIGHW*FAC*QFCN(2,1)*QFCN(5,2)
 147           CALL SIGFIL(SIG,2,5,30,27)
 148           SIG=.5*SIGHW*FAC*QFCN(5,1)*QFCN(2,2)
 149           CALL SIGFIL(SIG,5,2,30,27)
 150           SIG=.5*SIGHW*FAC*QFCN(8,1)*QFCN(7,2)
 151           CALL SIGFIL(SIG,8,7,30,27)
 152           SIG=.5*SIGHW*FAC*QFCN(7,1)*QFCN(8,2)
 153           CALL SIGFIL(SIG,7,8,30,27)
 154 130       CONTINUE
 155       END IF
 156 C
 157 C          ZH production via Z*
 158 C
 159       IF (GOQ(29,1).AND.GOQ(30,2)) THEN
 160           CALL TWOKIN(0.,0.,AMZ,AMH)
 161           IF(X1.GE.1..OR.X2.GE.1.) GO TO 200
 162           E1=SQRT(P(1)**2+AMZ2)
 163           E2=SQRT(P(2)**2+AMH2)
 164           FAC=1./(3.*PI*S**2)
 165           FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
 166           PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
 167           DO 210 IQ1=2,11
 168             IFLQ=IS2UD(IQ1)
 169             IQ2=MATCH(IQ1,4)
 170             IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 210
 171             SIG=GF**2*AMZ**8*(GV(IFLQ)**2+GA(IFLQ)**2)*
 172      $     (S/AMZ2+(1.-T/AMZ2)*(1.-U/AMZ2))/PROPZ*TBRWW(4,1)
 173             SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
 174             CALL SIGFIL(SIG,IQ1,IQ2,29,30)
 175 210         CONTINUE
 176 200       CONTINUE
 177       END IF
 178 C          HZ production via Z*
 179 C
 180       IF (GOQ(30,1).AND.GOQ(29,2)) THEN
 181           CALL TWOKIN(0.,0.,AMH,AMZ)
 182           IF(X1.GE.1..OR.X2.GE.1.) GO TO 220
 183           E1=SQRT(P(1)**2+AMH2)
 184           E2=SQRT(P(2)**2+AMZ2)
 185           FAC=1./(3.*PI*S**2)
 186           FAC=FAC*S/SCM*(P(1)*P(2)/(E1*E2))*UNITS
 187           PROPZ=(S-AMZ**2)**2+AMZ**2*GAMZ**2
 188           DO 230 IQ1=2,11
 189             IFLQ=IS2UD(IQ1)
 190             IQ2=MATCH(IQ1,4)
 191             IF (IQ2.EQ.0.OR.IQ2.GE.12) GO TO 230
 192             SIG=GF**2*AMZ**8*(GV(IFLQ)**2+GA(IFLQ)**2)*
 193      $     (S/AMZ2+(1.-T/AMZ2)*(1.-U/AMZ2))/PROPZ*TBRWW(4,2)
 194             SIG=.5*SIG*FAC*QFCN(IQ1,1)*QFCN(IQ2,2)
 195             CALL SIGFIL(SIG,IQ1,IQ2,30,29)
 196 230         CONTINUE
 197 220       CONTINUE
 198       END IF
 199       RETURN
 200       END