JIMMY first commit.
[u/mrichter/AliRoot.git] / HERWIG / jimmy / jimmy / eikon.F
1       DOUBLE PRECISION FUNCTION EIKON( DUMMY )
2 C     --------------------------------------------------------------------
3 C     Purpose: To eikonalise the total hadronic cross section pass the
4 C     cross section to be eikonalised (XSECN) and the routine returns
5 C     the eikonalised version (EIKON)
6 C     --------------------------------------------------------------------
7       IMPLICIT NONE
8
9 #include "jimmy.inc"
10
11       DOUBLE PRECISION EPS, XSECN, DUMMY
12       INTEGER NDIM, MAXPTS, MINPTS, LENWRK, I
13       PARAMETER(NDIM=2,MAXPTS=10000*NDIM,
14      & LENWRK=(NDIM+2)*(1+MAXPTS/(2**NDIM+2*NDIM*NDIM+2*NDIM+1)))
15       REAL A(NDIM), B(NDIM), RESULT, ERROR
16
17       COMMON/BLOCK1/XSECN               !COMMON local to eikfun
18
19       XSECN = DUMMY
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.D-3
27       FN_TYPE = 3
28       CALL PARTN( NDIM, A, B, 1.0, 10000 )
29       CALL INTGRL( NDIM, 0, 200, RESULT, ERROR )
30       EIKON = DBLE(RESULT)
31
32       RETURN
33       END
34
35