]>
Commit | Line | Data |
---|---|---|
ef94df36 | 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 |