ISAJET/code/sigh2.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE SIGH2
   3 C
   4 C          COMPUTE THE WEINBERG-SALAM HIGGS CROSS SECTION
   5 C          D(SIGMA)/D(QMW**2)D(YW)D(OMEGA)
   6 C          FOR THE SPECIFIED JET TYPES. TRIVIAL EXCEPT FOR W W FUSION,
   7 C          WHICH HAS INTERFERENCE WITH W W SCATTERING.
   8 C
   9 #include "isajet/itapes.inc"
  10 #include "isajet/qcdpar.inc"
  11 #include "isajet/jetpar.inc"
  12 #include "isajet/pjets.inc"
  13 #include "isajet/primar.inc"
  14 #include "isajet/q1q2.inc"
  15 #include "isajet/jetsig.inc"
  16 #include "isajet/wsig.inc"
  17 #include "isajet/qsave.inc"
  18 #include "isajet/wcon.inc"
  19 #include "isajet/const.inc"
  20 #include "isajet/hcon.inc"
  21 C
  22       DIMENSION X(2),LISTJ(29),WTHELI(4)
  23       EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT),(X(1),X1)
  24 #if defined(CERNLIB_DOUBLE)
  25       DOUBLE PRECISION C,TERM,SUM,DENOM,ZCM
  26 #endif
  27 C
  28 C          WTHELI ARE WEIGHTS OF HELICITY AMPLITUDES IN SIGMA.
  29       DATA WTHELI/1.,2.,2.,4./
  30       DATA LISTJ/
  31      $9,1,-1,2,-2,3,-3,4,-4,5,-5,6,-6,
  32      $11,-11,12,-12,13,-13,14,-14,15,-15,16,-16,
  33      $10,80,-80,90/
  34 C
  35 C          QUARK OR GLUON FUSION TO HIGGS
  36 C
  37       IF(INITYP(1).LE.25) THEN
  38         SIGLLQ=SIGEVT/(4.*PI)
  39         RETURN
  40       ENDIF
  41 C
  42 C          W+W FUSION AND W+W->W+W IN EFFECTIVE W APPROXIMATION.
  43 C
  44 C          KINEMATICS
  45       IFL1=LISTJ(JETTYP(1))
  46       IFL2=LISTJ(JETTYP(2))
  47       IFIN1=LISTJ(INITYP(1))
  48       IFIN2=LISTJ(INITYP(2))
  49       WMF=AMASS(IFL1)
  50       WMI=AMASS(IFIN1)
  51       PINPF=SQRT((S-4.*WMI**2)*(S-4.*WMF**2))
  52       ZCM=(.5*S+T-WMI**2-WMF**2)/(.5*PINPF)
  53 C          RESET COEFFICIENTS FOR SELECTED PROCESS
  54       IABSI=IABS(IFIN1)
  55       IABSF=IABS(IFL1)
  56       IF(IABSI.EQ.80) THEN
  57         IF(IABSF.EQ.80) THEN
  58           CALL XWWWW
  59         ELSE
  60           CALL XWWZZ
  61         ENDIF
  62       ELSE
  63         IF(IABSF.EQ.80) THEN
  64           CALL XZZWW
  65         ELSE
  66           CALL XZZZZ
  67         ENDIF
  68       ENDIF
  69 C          SUM CROSS SECTION TERMS. I,J RUN OVER AMPLITUDE TERMS.
  70 C          L RUNS OVER HELICITY STATES. N RUNS OVER POWERS.
  71 C          REMEMBER THAT L=4 IS MISSING SIN(THETA)/SQRT(2)
  72       SUM=0.
  73       DO 111 I=1,4
  74       DO 111 J=I,4
  75       DENOM=1./((ADWWWW(1,I)+ADWWWW(2,I)*ZCM)
  76      $*(ADWWWW(1,J)+ADWWWW(2,J)*ZCM))
  77         DO 112 L=1,4
  78         TERM=0.
  79           DO 113 N=0,6
  80           C=0.
  81           N1=MAX(N-3,0)
  82           N2=MIN(3,N)
  83           DO 114 K=N1,N2
  84 114       C=C+ANWWWW(K+1,I,L)*ANWWWW(N-K+1,J,L)
  85           C=C*WTHELI(L)
  86           IF(J.NE.I) C=2.*C
  87           TERM=TERM+C*ZCM**N
  88 113       CONTINUE
  89         IF(L.EQ.4) TERM=TERM*(1.-ZCM**2)/2.
  90         TERM=TERM*DENOM
  91         SUM=SUM+TERM
  92 112     CONTINUE
  93 111   CONTINUE
  94 C          ADD IMAGINARY PART SQUARED.
  95       SUM=SUM+WTHELI(1)*AIWWWW(1)**2+WTHELI(2)*AIWWWW(2)**2
  96      $+WTHELI(3)*AIWWWW(3)**2+WTHELI(4)*AIWWWW(4)**2
  97 C          CROSS SECTION. NOTE D(OMEGA)=2.*PI*D(Z)
  98       SIG0=SUM/(64.*PI**2*S*SCM)*UNITS
  99       SIG0=SIG0*TBRWW(JETTYP(1)-25,1)*TBRWW(JETTYP(2)-25,2)
 100 C          SYMMETRY FACTOR
 101       IF(IABSF.EQ.90) SIG0=.5*SIG0
 102       SIGLLQ=SIG0*QSAVE(INITYP(1),1)*QSAVE(INITYP(2),2)
 103       RETURN
 104       END