GEANT321/geant321/nuclev.inc

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:41  cernlib
   6 * Geant
   7 *
   8 *
   9 #ifndef CERNLIB_GEANT321_NUCLEV_INC
  10 #define CERNLIB_GEANT321_NUCLEV_INC
  11 *
  12 *
  13 * nuclev.inc
  14 *
  15 *
  16 *=== nuclev ===========================================================*
  17 *
  18 *----------------------------------------------------------------------*
  19 *                                                                      *
  20 *     NUClear LEVel common:                                            *
  21 *                                                                      *
  22 *     Created on  10 january 1992  by    Alfredo Ferrari & Paola Sala  *
  23 *                                                   Infn - Milan       *
  24 *                                                                      *
  25 *     Last change on 26-mar-93     by    Alfredo Ferrari               *
  26 *                                                                      *
  27 *     Included in the following routines:                              *
  28 *                                                                      *
  29 *                        BERTTP                                        *
  30 *                        BIMSEL                                        *
  31 *                        NCLVIN                                        *
  32 *                        NCLVGT                                        *
  33 *                        NCLVST                                        *
  34 *                        NUCNUC                                        *
  35 *                        NWISEL                                        *
  36 *                        PFNCLV                                        *
  37 *                        PIOABS                                        *
  38 *                                                                      *
  39 *     Description of the variables (nuclev):                           *
  40 *                                                                      *
  41 *                   paenuc(k,i) = pairing energies (GeV), for i=1 pro- *
  42 *                                 ton pairing energy for Z=k, for i=2  *
  43 *                                 neutron pairing energy for N=k       *
  44 *                   shenuc(k,i) = shell energies (GeV), for i=1 pro-   *
  45 *                                 ton shell energy for Z=k, for i=2    *
  46 *                                 neutron shell energy for N=k         *
  47 *                     defrmi(i) = level spacing at the Fermi level ac- *
  48 *                                 cording to a Fermi gas               *
  49 *                     defmag(i) = extra energy needed to overcome the  *
  50 *                                 gap between the closed shell and the *
  51 *                                 next level for nuclei with magic     *
  52 *                                 numbers of protons or neutrons       *
  53 *                   cumrad(k,i) = sum for j=1,k of 1 /Rj where Ri is   *
  54 *                                 the maximum radius at which can be   *
  55 *                                 found the jth nucleon of type i      *
  56 *                                 (i=1 proton, i=2 neutron)            *
  57 *                        ruspro = maximum radius at which protons can  *
  58 *                                 still be found                       *
  59 *                        rusneu = maximum radius at which neutrons can *
  60 *                                 still be found                       *
  61 *                        juspro = index for flagging that the jth      *
  62 *                                 proton has been already used if      *
  63 *                                 = inuclv                             *
  64 *                        jusneu = index for flagging that the jth      *
  65 *                                 neutron has been already used if     *
  66 *                                 = inuclv                             *
  67 *                        ntapro = number of protons in the target      *
  68 *                        ntaneu = number of neutrons in the target     *
  69 *                        navpro = number of available proton levels in *
  70 *                                 the target                           *
  71 *                        navneu = number of available neutron levels   *
  72 *                                 in the target                        *
  73 *                        nlspro = number of protons in the last level  *
  74 *                        nlsneu = number of neutrons in the last level *
  75 *                        ncopro = number of protons in the core        *
  76 *                        nconeu = number of neutrons in the core       *
  77 *                        nskpro = number of protons in the skin        *
  78 *                        nskneu = number of neutrons in the skin       *
  79 *                        nhapro = number of protons in the halo        *
  80 *                        nhaneu = number of neutrons in the halo       *
  81 *                        nuspro = number of "used" protons             *
  82 *                        nusneu = number of "used" neutrons            *
  83 *                        jmxpro = index  of the highest energy avai-   *
  84 *                                 lable proton                         *
  85 *                        jmxneu = index  of the highest energy avai-   *
  86 *                                 lable neutron                        *
  87 *                        nprnuc = number of presently used nucleons    *
  88 *                                 (max 2)                              *
  89 *                        iprnuc = type of the presently considered     *
  90 *                                 nucleon                              *
  91 *                        jprnuc = index of the presently considered    *
  92 *                                 nucleon                              *
  93 *                        magnum = magic numbers                        *
  94 *                        magpro = last (possibly not filled) proton    *
  95 *                                 shell                                *
  96 *                        magneu = last (possibly not filled) neutron   *
  97 *                                 shell                                *
  98 *                     mgspro(k) = status of the kth proton shell:      *
  99 *                                 -2 not filled but paired             *
 100 *                                 -1 not filled and unpaired           *
 101 *                                  0 filled (paired of course)         *
 102 *                                 +1 originally  filled but now unpai- *
 103 *                                    red                               *
 104 *                     mgsneu(k) = status of the kth neutron shell      *
 105 *                        inuclv = index of the present event           *
 106 *                                                                      *
 107 *----------------------------------------------------------------------*
 108 *
 109 *  NUClear LEVel
 110       COMMON / FKNULE / PAENUC (200,2), SHENUC (200,2), DEFRMI (2),
 111      &                  DEFMAG (2), ENNCLV (160,2), RANCLV (160,2),
 112      &                  CUMRAD (0:160,2), RUSNUC (2),
 113      &                  ENPLVL (114), ENNLVL(164), JUSNUC (160,2),
 114      &                  NTANUC (2), NAVNUC (2), NLSNUC (2), NCONUC (2),
 115      &                  NSKNUC (2), NHANUC (2), NUSNUC (2), JMXNUC (2),
 116      &                  IPRNUC (2), JPRNUC (2), MAGNUM (8), MAGNUC (2),
 117      &                  MGSNUC (8,2), MGSSNC (25,2), NSBSHL (2),
 118      &                  NPRNUC, INUCLV
 119       DIMENSION JUSPRO (160), JUSNEU (160), MGSPRO (8), MGSNEU (8),
 120      &          MGSSPR (19) , MGSSNE (25)
 121       EQUIVALENCE ( RUSNUC (1), RUSPRO )
 122       EQUIVALENCE ( RUSNUC (2), RUSNEU )
 123       EQUIVALENCE ( JUSNUC (1,1), JUSPRO (1) )
 124       EQUIVALENCE ( JUSNUC (1,2), JUSNEU (1) )
 125       EQUIVALENCE ( MGSNUC (1,1), MGSPRO (1) )
 126       EQUIVALENCE ( MGSNUC (1,2), MGSNEU (1) )
 127       EQUIVALENCE ( MGSSNC (1,1), MGSSPR (1) )
 128       EQUIVALENCE ( MGSSNC (1,2), MGSSNE (1) )
 129       EQUIVALENCE ( NTANUC (1), NTAPRO )
 130       EQUIVALENCE ( NTANUC (2), NTANEU )
 131       EQUIVALENCE ( NAVNUC (1), NAVPRO )
 132       EQUIVALENCE ( NAVNUC (2), NAVNEU )
 133       EQUIVALENCE ( NLSNUC (1), NLSPRO )
 134       EQUIVALENCE ( NLSNUC (2), NLSNEU )
 135       EQUIVALENCE ( NCONUC (1), NCOPRO )
 136       EQUIVALENCE ( NCONUC (2), NCONEU )
 137       EQUIVALENCE ( NSKNUC (1), NSKPRO )
 138       EQUIVALENCE ( NSKNUC (2), NSKNEU )
 139       EQUIVALENCE ( NHANUC (1), NHAPRO )
 140       EQUIVALENCE ( NHANUC (2), NHANEU )
 141       EQUIVALENCE ( NUSNUC (1), NUSPRO )
 142       EQUIVALENCE ( NUSNUC (2), NUSNEU )
 143       EQUIVALENCE ( JMXNUC (1), JMXPRO )
 144       EQUIVALENCE ( JMXNUC (2), JMXNEU )
 145       EQUIVALENCE ( MAGNUC (1), MAGPRO )
 146       EQUIVALENCE ( MAGNUC (2), MAGNEU )
 147
 148
 149 #endif