This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / fluka / geta.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:19:56  cernlib
6 * Geant
7 *
8 *
9 #include "geant321/pilot.h"
10 *CMZ :  3.21/02 29/03/94  15.41.42  by  S.Giani
11 *-- Author :
12 *$ CREATE GETA.FOR
13 *COPY GETA
14 *
15 *=== geta =============================================================*
16 *
17       FUNCTION GETA ( U, JZ, JN, MODE, IS, AOGMAX,
18      &                                 AOGMIN )
19  
20 #include "geant321/dblprc.inc"
21 #include "geant321/dimpar.inc"
22 #include "geant321/iounit.inc"
23 *
24 *----------------------------------------------------------------------*
25 *                                                                      *
26 *     Created on 18 january 1993   by    Alfredo Ferrari & Paola Sala  *
27 *                                                   Infn - Milan       *
28 *                                                                      *
29 *     Last change on 28-jan-93     by    Alfredo Ferrari               *
30 *                                                                      *
31 *----------------------------------------------------------------------*
32 *
33 #include "geant321/eva0.inc"
34 #if !defined(CERNLIB_SINGLE)
35       PARAMETER ( ASMTOG = 6.0D0 / PIPIPI**2 )
36 #endif
37 #if defined(CERNLIB_SINGLE)
38       PARAMETER ( ASMTOG = 6.0e0 / PIPIPI**2 )
39 #endif
40       LOGICAL LASMLL
41 *
42       LASMLL = .TRUE.
43       IF ( JZ .LE. 0 .OR. JN .LE. 0 ) THEN
44          GETA = ( JZ + JN) / B0
45          AOGMAX = GETA
46          AOGMIN = GETA
47          RETURN
48       END IF
49       GO TO 5
50       ENTRY GETG ( U, JZ, JN, MODE, IS, AOGMAX, AOGMIN )
51       LASMLL = .FALSE.
52     5 CONTINUE
53       JA = JN + JZ
54       AA = JA
55       ZZ = JZ
56       GO TO (10,20,30), MODE
57 *  +-------------------------------------------------------------------*
58 *  |  Standard EVAP parametrization for the level density
59    10 CONTINUE
60          TEMP = AA * ( ONEONE + Y0 * ( ( AA - TWOTWO*ZZ ) / AA )**2) /B0
61          AOGMAX = TEMP
62          AOGMIN = TEMP
63       GO TO 40
64 *  |
65 *  +-------------------------------------------------------------------*
66 *  |  Gilbert & Cameron level density (with Z and A dependent correction
67 *  |  smoothly going into the Ignatyuk high energy limit, according
68 *  |  to the Ignatyuk interpolation formula
69    20 CONTINUE
70          STOP 'GETA:20'
71 ***      GO TO 40
72 *  |
73 *  +-------------------------------------------------------------------*
74 *  |  Julich zero energy A-dependent level density:
75    30 CONTINUE
76          STOP 'GETA:30'
77    40 CONTINUE
78 *  |
79 *  +-------------------------------------------------------------------*
80 *  +-------------------------------------------------------------------*
81 *  |  geta = a (U) (aogmax = a_max)
82       IF ( LASMLL ) THEN
83          GETA   = TEMP
84 *  |
85 *  +-------------------------------------------------------------------*
86 *  |  getg = g (U) (aogmax = g_max)
87       ELSE
88          GETG   = ASMTOG * TEMP
89          AOGMAX = ASMTOG * AOGMAX
90          AOGMIN = ASMTOG * AOGMIN
91       END IF
92 *  |
93 *  +-------------------------------------------------------------------*
94 *=== End of Function geta =============================================*
95       RETURN
96       END