[u/mrichter/AliRoot.git] / GEANT321 / geant321 / nuclev.inc

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:20:41  cernlib
* Geant
*
*
#ifndef CERNLIB_GEANT321_NUCLEV_INC
#define CERNLIB_GEANT321_NUCLEV_INC
*
*
* nuclev.inc
*
*
*=== nuclev ===========================================================*
*
*----------------------------------------------------------------------*
*                                                                      *
*     NUClear LEVel common:                                            *
*                                                                      *
*     Created on  10 january 1992  by    Alfredo Ferrari & Paola Sala  *
*                                                   Infn - Milan       *
*                                                                      *
*     Last change on 26-mar-93     by    Alfredo Ferrari               *
*                                                                      *
*     Included in the following routines:                              *
*                                                                      *
*                        BERTTP                                        *
*                        BIMSEL                                        *
*                        NCLVIN                                        *
*                        NCLVGT                                        *
*                        NCLVST                                        *
*                        NUCNUC                                        *
*                        NWISEL                                        *
*                        PFNCLV                                        *
*                        PIOABS                                        *
*                                                                      *
*     Description of the variables (nuclev):                           *
*                                                                      *
*                   paenuc(k,i) = pairing energies (GeV), for i=1 pro- *
*                                 ton pairing energy for Z=k, for i=2  *
*                                 neutron pairing energy for N=k       *
*                   shenuc(k,i) = shell energies (GeV), for i=1 pro-   *
*                                 ton shell energy for Z=k, for i=2    *
*                                 neutron shell energy for N=k         *
*                     defrmi(i) = level spacing at the Fermi level ac- *
*                                 cording to a Fermi gas               *
*                     defmag(i) = extra energy needed to overcome the  *
*                                 gap between the closed shell and the *
*                                 next level for nuclei with magic     *
*                                 numbers of protons or neutrons       *
*                   cumrad(k,i) = sum for j=1,k of 1 /Rj where Ri is   *
*                                 the maximum radius at which can be   *
*                                 found the jth nucleon of type i      *
*                                 (i=1 proton, i=2 neutron)            *
*                        ruspro = maximum radius at which protons can  *
*                                 still be found                       *
*                        rusneu = maximum radius at which neutrons can *
*                                 still be found                       *
*                        juspro = index for flagging that the jth      *
*                                 proton has been already used if      *
*                                 = inuclv                             *
*                        jusneu = index for flagging that the jth      *
*                                 neutron has been already used if     *
*                                 = inuclv                             *
*                        ntapro = number of protons in the target      *
*                        ntaneu = number of neutrons in the target     *
*                        navpro = number of available proton levels in *
*                                 the target                           *
*                        navneu = number of available neutron levels   *
*                                 in the target                        *
*                        nlspro = number of protons in the last level  *
*                        nlsneu = number of neutrons in the last level *
*                        ncopro = number of protons in the core        *
*                        nconeu = number of neutrons in the core       *
*                        nskpro = number of protons in the skin        *
*                        nskneu = number of neutrons in the skin       *
*                        nhapro = number of protons in the halo        *
*                        nhaneu = number of neutrons in the halo       *
*                        nuspro = number of "used" protons             *
*                        nusneu = number of "used" neutrons            *
*                        jmxpro = index  of the highest energy avai-   *
*                                 lable proton                         *
*                        jmxneu = index  of the highest energy avai-   *
*                                 lable neutron                        *
*                        nprnuc = number of presently used nucleons    *
*                                 (max 2)                              *
*                        iprnuc = type of the presently considered     *
*                                 nucleon                              *
*                        jprnuc = index of the presently considered    *
*                                 nucleon                              *
*                        magnum = magic numbers                        *
*                        magpro = last (possibly not filled) proton    *
*                                 shell                                *
*                        magneu = last (possibly not filled) neutron   *
*                                 shell                                *
*                     mgspro(k) = status of the kth proton shell:      *
*                                 -2 not filled but paired             *
*                                 -1 not filled and unpaired           *
*                                  0 filled (paired of course)         *
*                                 +1 originally  filled but now unpai- *
*                                    red                               *
*                     mgsneu(k) = status of the kth neutron shell      *
*                        inuclv = index of the present event           *
*                                                                      *
*----------------------------------------------------------------------*
*
*  NUClear LEVel
      COMMON / FKNULE / PAENUC (200,2), SHENUC (200,2), DEFRMI (2),
     &                  DEFMAG (2), ENNCLV (160,2), RANCLV (160,2),
     &                  CUMRAD (0:160,2), RUSNUC (2),
     &                  ENPLVL (114), ENNLVL(164), JUSNUC (160,2),
     &                  NTANUC (2), NAVNUC (2), NLSNUC (2), NCONUC (2),
     &                  NSKNUC (2), NHANUC (2), NUSNUC (2), JMXNUC (2),
     &                  IPRNUC (2), JPRNUC (2), MAGNUM (8), MAGNUC (2),
     &                  MGSNUC (8,2), MGSSNC (25,2), NSBSHL (2),
     &                  NPRNUC, INUCLV
      DIMENSION JUSPRO (160), JUSNEU (160), MGSPRO (8), MGSNEU (8),
     &          MGSSPR (19) , MGSSNE (25)
      EQUIVALENCE ( RUSNUC (1), RUSPRO )
      EQUIVALENCE ( RUSNUC (2), RUSNEU )
      EQUIVALENCE ( JUSNUC (1,1), JUSPRO (1) )
      EQUIVALENCE ( JUSNUC (1,2), JUSNEU (1) )
      EQUIVALENCE ( MGSNUC (1,1), MGSPRO (1) )
      EQUIVALENCE ( MGSNUC (1,2), MGSNEU (1) )
      EQUIVALENCE ( MGSSNC (1,1), MGSSPR (1) )
      EQUIVALENCE ( MGSSNC (1,2), MGSSNE (1) )
      EQUIVALENCE ( NTANUC (1), NTAPRO )
      EQUIVALENCE ( NTANUC (2), NTANEU )
      EQUIVALENCE ( NAVNUC (1), NAVPRO )
      EQUIVALENCE ( NAVNUC (2), NAVNEU )
      EQUIVALENCE ( NLSNUC (1), NLSPRO )
      EQUIVALENCE ( NLSNUC (2), NLSNEU )
      EQUIVALENCE ( NCONUC (1), NCOPRO )
      EQUIVALENCE ( NCONUC (2), NCONEU )
      EQUIVALENCE ( NSKNUC (1), NSKPRO )
      EQUIVALENCE ( NSKNUC (2), NSKNEU )
      EQUIVALENCE ( NHANUC (1), NHAPRO )
      EQUIVALENCE ( NHANUC (2), NHANEU )
      EQUIVALENCE ( NUSNUC (1), NUSPRO )
      EQUIVALENCE ( NUSNUC (2), NUSNEU )
      EQUIVALENCE ( JMXNUC (1), JMXPRO )
      EQUIVALENCE ( JMXNUC (2), JMXNEU )
      EQUIVALENCE ( MAGNUC (1), MAGPRO )
      EQUIVALENCE ( MAGNUC (2), MAGNEU )
 

#endif
Commit	Line	Data
fe4da5cc	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