HERWIG/jimmy/jimmy/jmsn.F

   1       DOUBLE PRECISION FUNCTION JMSN( DUMMY1, DUMMY2 )
   2 C --------------------------------------------------------------------
   3 C  Purpose: To calculate the cross section for N and
   4 C           only N scatters.
   5 C --------------------------------------------------------------------
   6
   7       IMPLICIT NONE
   8
   9 #include "jimmy.inc"
  10
  11
  12       INTEGER N, DUMMY2
  13       DOUBLE PRECISION JMSNFN, EPS, XSECN, DUMMY1
  14       INTEGER NDIM, MAXPTS, MINPTS, LENWRK, I
  15       PARAMETER(NDIM=2,MAXPTS=10000*NDIM,
  16      *   LENWRK=(NDIM+2)*(1+MAXPTS/(2**NDIM+2*NDIM*NDIM+2*NDIM+1)))
  17
  18       REAL A(NDIM), B(NDIM), RESULT, ERROR
  19
  20       EXTERNAL JMSNFN
  21       COMMON/BLOCK2/ XSECN, N       ! COMMON local to JMSNFN
  22
  23       XSECN = DUMMY1
  24       N = DUMMY2
  25       DO 10 I=1,NDIM
  26         A(I)=0.0
  27         B(I)=1.0
  28   10  CONTINUE
  29
  30       MINPTS = 0
  31       EPS    = 1.D-3
  32       FN_TYPE = 4
  33       CALL PARTN( NDIM, A, B, 1.0, 10000 )
  34       CALL INTGRL( NDIM, 0, 200, RESULT, ERROR )
  35       IF (JMBUG.GT.0) write(*,*) 'jmsn: result, error',result,error
  36       JMSN = DBLE(RESULT)
  37
  38       RETURN
  39       END