This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / fluka / eexlvl.F
1 *
2 * $Id$
3 *
4 * $Log$
5 * Revision 1.1.1.1  1995/10/24 10:19:55  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 EEXLVL.FOR
13 *COPY EEXLVL
14 *
15 *=== eexlvl ===========================================================*
16 *
17       SUBROUTINE EEXLVL ( JA, JZ, EEX1ST, EEX2ND, EEXCON )
18  
19 #include "geant321/dblprc.inc"
20 #include "geant321/dimpar.inc"
21 #include "geant321/iounit.inc"
22 *
23 *----------------------------------------------------------------------*
24 *                                                                      *
25 *     Created on 06 december 1991  by    Alfredo Ferrari & Paola Sala  *
26 *                                                   Infn - Milan       *
27 *                                                                      *
28 *     Last change on 28-apr-92     by    Alfredo Ferrari               *
29 *                                                                      *
30 *                                                                      *
31 *----------------------------------------------------------------------*
32 *
33 #include "geant321/eva0.inc"
34 *
35       SAVE SQATAR, JAOLD
36       DATA JAOLD / 0 /
37 *
38       IF ( JA .EQ. JZ ) THEN
39          EEXCON = 0.D+00
40          EEX1ST = 0.D+00
41          EEX2ND = 0.D+00
42          RETURN
43       END IF
44       IF ( JZ .EQ. 0 .OR. JA .EQ. JZ ) THEN
45          EEXCON = 0.D+00
46       ELSE
47          EEXCON = 1.D-03 * ( CAM4 (JZ) + CAM5 (JA-JZ) )
48       END IF
49 *  **** Very tentative selection of Eex1st, Eex2nd, based on pairing
50 *       energies according to delta = 12 MeV / A^1/2 ****
51       IZODD  = 1 - MOD ( JZ, 2 )
52       INODD  = 1 - MOD ( JA-JZ, 2 )
53       IODD   = IZODD + INODD
54 *  +-------------------------------------------------------------------*
55 *  |  Even-even nucleus
56       IF ( IODD .GE. 2 ) THEN
57 *  |  +----------------------------------------------------------------*
58 *  |  |
59          IF ( JA .NE. JAOLD ) THEN
60             JAOLD  = JA
61             DJA = JA
62             SQATAR = SQRT ( DJA )
63          END IF
64 *  |  |
65 *  |  +----------------------------------------------------------------*
66          EEX1ST = 12.0D-03 / SQATAR
67          EEX2ND = 2.D+00 * EEX1ST
68 *  |
69 *  +-------------------------------------------------------------------*
70 *  |  even-odd nucleus
71       ELSE IF ( IODD .GT. 0 ) THEN
72 *  |  +----------------------------------------------------------------*
73 *  |  |
74          IF ( JA .NE. JAOLD ) THEN
75             JAOLD  = JA
76             DJA    = JA
77             SQATAR = SQRT ( DJA )
78          END IF
79 *  |  |
80 *  |  +----------------------------------------------------------------*
81          EEX1ST = 12.0D-03 / SQATAR
82          EEX2ND = EEX1ST
83 *  |
84 *  +-------------------------------------------------------------------*
85 *  |  odd-odd nucleus
86       ELSE
87          EEX2ND = 0.D+00
88          EEX1ST = 0.D+00
89       END IF
90 *  |
91 *  +-------------------------------------------------------------------*
92       RETURN
93 *=== End of subroutine eexlvl =========================================*
94       END