This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / 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