HERWIG/jimmy/jimmy/qcdres.F

   1       DOUBLE PRECISION FUNCTION QCDRES( Z )
   2 C     ---------------------------------------------------------------------
   3 C     Purpose: To calculate the inclusive hadron(or photon)+hadron(or
   4 C     photon) => JETS
   5 C     ---------------------------------------------------------------------
   6 #include "herwig65.inc"
   7 #include "jimmy.inc"
   8
   9       INTEGER I, MINPTS, NDIM, MAXPTS, LENWRK
  10
  11 C     Set parameters for integration
  12       PARAMETER(NDIM=3,MAXPTS=10000*NDIM,
  13      *   LENWRK=(NDIM+2)*(1+MAXPTS/(2**NDIM+2*NDIM*NDIM+2*NDIM+1)))
  14       DOUBLE PRECISION Z, EPS
  15       EXTERNAL JMKERN
  16       REAL A(NDIM), B(NDIM)
  17       REAL RESULT, ERROR
  18
  19 C     INITIALISE INTEGRATION LIMITS (ALL 0 TO 1)
  20       DO 10 I=1, NDIM
  21         A(I)=0.0
  22         B(I)=1.0
  23   10  CONTINUE
  24
  25       MINPTS=0
  26       EPS=1.0D-3
  27 C     Store the current z value
  28       YGAMMA=Z
  29
  30       FN_TYPE=101
  31       CALL PARTN( NDIM, A, B, 1.0, 10000 )
  32       CALL INTGRL( NDIM, 0, 200, RESULT, ERROR )
  33       QCDRES=DBLE(RESULT)
  34
  35 C     Re-include the PHAD factor (JMKERN divides out the probability of the
  36 C     photon becoming a hadronic = 1/phad)
  37       QCDRES = QCDRES/PHAD
  38
  39       RETURN
  40
  41       END
  42
  43
  44
  45
  46
  47
  48