This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / giface / gpfisi.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:21:15  cernlib
6 * Geant
7 *
8 *
9 #include "geant321/pilot.h"
10 *CMZ :  3.21/02 29/03/94  15.41.38  by  S.Giani
11 *-- Author :
12       SUBROUTINE GPFISI
13 C
14 C *** CALCULATION OF PHOTO-FISSION AND PHOTO-ABSORBTION         ***
15 C *** CROSS SECTIONS IN CURRENT MATERIAL                        ***
16 C *** HMF 25-AUG-1989 RWTH AACHEN / NVE 11-MAY-1990 CERN GENEVA ***
17 C
18 C CALLED BY : GPHYSI
19 C ORIGIN    : F.CARMINATI, H.FESEFELDT
20 C
21 #include "geant321/gcbank.inc"
22 #include "geant321/gcjloc.inc"
23 #include "geant321/gcmulo.inc"
24 #include "geant321/gcmate.inc"
25 #include "geant321/gconsp.inc"
26 #include "geant321/gctrak.inc"
27 C
28       DIMENSION CSGN(19),CSGF(19),EKFISS(19)
29       SAVE CSGN,CSGF,EKFISS
30 C
31 C     MEASURED CROSS SECTIONS ON U(238) IN MILLBARN, ENERGIES IN GEV.
32 C
33       DATA CSGN  /   0.,  12.,  20.,  46., 115., 260., 450., 520.,
34      +             400., 270., 200., 150., 120.,  80.,  55.,  27.,
35      +              15.,   1.,   0./
36       DATA CSGF  /   0.,   0.,   0.,  14.,  30.,  65., 100., 130.,
37      +             150., 180., 180., 150., 120.,  80.,  55.,  27.,
38      +              15.,   1.,   0./
39       DATA EKFISS/0.004,0.005,0.006,0.007,0.008,0.009,0.010,0.011,
40      +            0.012,0.013,0.014,0.015,0.016,0.018,0.020,0.025,
41      +            0.030,0.035,0.040/
42 C
43       SIG=0.
44       SIGFIS=0.
45       IF(JPFIS.LE.0)                  GOTO 99
46       IF(ELOW(IEKBIN).LE.EKFISS(1))   GOTO 90
47       IF(ELOW(IEKBIN).GE.EKFISS(19))  GOTO 90
48       IF(A.LT.230..OR.A.GT.240.)      GOTO 90
49 C
50 C     A-DEPENDENCE, MEASURED AT THE MAXIMUM OF THE RESONANCE
51 C
52       Z43BA=-67.0+38.7*Z**(4./3.)/A
53       IF(Z43BA.LE.0.)                 GOTO 90
54 C
55 C     ENERGY DEPENDENCE, AS TABULATED IN THE DATA STATEMENTS
56 C
57       DO 10 I=2,19
58         IF(ELOW(IEKBIN).LT.EKFISS(I)) GOTO 20
59    10 CONTINUE
60       I=19
61    20 RAT = (EKFISS(I)-ELOW(IEKBIN))/(EKFISS(I)-EKFISS(I-1))
62 C
63       SIGRES = CSGN(I-1)*RAT + CSGN(I)*(1.-RAT)
64       SIGRES = SIGRES*Z43BA/0.53
65       SIGRES = SIGRES*0.0006022045*DENS/A
66 C
67       SIGFIS = CSGF(I-1)*RAT + CSGF(I)*(1.-RAT)
68       SIGFIS = SIGFIS*Z43BA/0.53
69       SIGFIS = SIGFIS*0.0006022045*DENS/A
70 C
71       SIG    = SIGRES+SIGFIS
72 C
73    90 IF(SIG.GT.0.) THEN
74         Q(JPFIS+IEKBIN) = 1./SIG
75       ELSE
76         Q(JPFIS+IEKBIN) = BIG
77       ENDIF
78 C
79       IF(SIGFIS.GT.0.) THEN
80         Q(JPFIS+NEK1+IEKBIN) = 1./SIGFIS
81       ELSE
82         Q(JPFIS+NEK1+IEKBIN) = BIG
83       ENDIF
84 C
85   99  END