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