ISAJET/code/sigtc2.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE SIGTC2
   3 C
   4 C          Compute the techni-rho decay distribution cross section
   5 C          D(SIGMA)/D(QMW**2)D(YW)D(OMEGA)
   6 C          for the specified jet types. This is trivial but done for
   7 C          compatibility with Drell-Yan and Higgs.
   8 C
   9 #if defined(CERNLIB_IMPNONE)
  10       IMPLICIT NONE
  11 #endif
  12 #include "isajet/itapes.inc"
  13 #include "isajet/const.inc"
  14 #include "isajet/jetpar.inc"
  15 #include "isajet/jetsig.inc"
  16 #include "isajet/pjets.inc"
  17 #include "isajet/wsig.inc"
  18 #include "isajet/tcpar.inc"
  19 C
  20       REAL AM12,AM22,ANGFAC,S,T,U
  21       EQUIVALENCE (S,SHAT),(T,THAT),(U,UHAT)
  22 C
  23 C          Angfac is (1-z**2), and is determined in terms of S,T,U.
  24 C          Note that both rho+- and rho0 are always elastic.
  25       AM12=PJETS(5,1)**2
  26       AM22=PJETS(5,2)**2
  27       ANGFAC=4.*(T*U-AM12*AM22)/((S-AM12-AM22)**2-4.*AM12*AM22)
  28 C          Differential cross section
  29       SIGLLQ=SIGEVT*ANGFAC*3./(8.*PI)
  30       RETURN
  31       END