GEANT321/gheisha/cinema.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:03  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.39  by  S.Giani
  11 *-- Author :
  12       FUNCTION CINEMA(EK1)
  13 C
  14 C *** INELASTICITY IN NUCLEAR INTERACTIONS AS A FUNCTION ***
  15 C *** OF ATOMIC NUMBER ATNO2 AND KINETIC ENERGY EK1 ***
  16 C *** NVE 12-JUL-1988 CERN GENEVA ***
  17 C
  18 C ORIGIN : H.FESEFELDT (14-OCT-1987)
  19 C
  20 C THE FUNCTIONAL DEPENDENCE AND THE PARAMETERS HAVE BEEN OBTAINED
  21 C BY STUDY OF VARIOUS NUCLEAR STRUCTURE MODELS.
  22 C BUT: IT IS OF COURSE AN INTERPOLATION AS FUNCTION OF ATOMIC
  23 C      NUMBER, FOR CERTAIN NUCLEI A DIFFERENT DESCRIPTION MAY BE
  24 C      MORE ADEQUATE. DETAILED TESTS HAVE BEEN PERFORMED FOR
  25 C      FE, CU, PB ,U AND SOME MIXTURES LIKE NAI, BGO, CONCRETE.
  26 C
  27 #include "geant321/mxgkgh.inc"
  28 #include "geant321/s_curpar.inc"
  29 #include "geant321/s_result.inc"
  30 #include "geant321/s_mat.inc"
  31 #include "geant321/limits.inc"
  32 C
  33       CINEMA=0.0
  34       ND=IND+1
  35       ALA=LOG(ATNO2)
  36       ALEK1=LOG(EK1)
  37       SIG1=0.50
  38       SIG2=0.50
  39       EM=0.2390+0.0408*ALA**2
  40       IF (EM. GT. 1.0) EM=1.0
  41       CINEM=0.0019*ALA**3
  42       IF(CINEM.GT.0.15) CINEM=0.15
  43       IF (PARMAT(ND,10) .GE. 0.01) CINEM=CINEM*PARMAT(ND,10)
  44 C
  45       IF (ALEK1 .GT. EM) GO TO 1
  46 C
  47       CORR=-(ALEK1-EM)**2/(2.0*SIG1**2)
  48       IF (CORR .LT. EXPXL) CORR=EXPXL
  49       IF (CORR .GT. EXPXU) CORR=EXPXU
  50       DUM1=-EK1*CINEM
  51       DUM2=ABS(DUM1)
  52       DUM3=EXP(CORR)
  53       CINEMA=0.0
  54       IF (DUM2 .GE. 1.0) CINEMA=DUM1*DUM3
  55       IF ((DUM2 .LT. 1.0) .AND. (DUM3 .GT. 1.0E-10)) CINEMA=DUM1*DUM3
  56       GO TO 2
  57 C
  58  1    CONTINUE
  59       CORR=-(ALEK1-EM)**2/(2.0*SIG2**2)
  60       IF (CORR .LT. EXPXL) CORR=EXPXL
  61       IF (CORR .GT. EXPXU) CORR=EXPXU
  62       DUM1=-EK1*CINEM
  63       DUM2=ABS(DUM1)
  64       DUM3=EXP(CORR)
  65       CINEMA=0.0
  66       IF (DUM2 .GE. 1.0) CINEMA=DUM1*DUM3
  67       IF ((DUM2 .LT. 1.0) .AND. (DUM3 .GT. 1.0E-10)) CINEMA=DUM1*DUM3
  68 C
  69  2    CONTINUE
  70       IF (CINEMA .LT. -EK1) CINEMA=-EK1
  71 C
  72       END