[u/mrichter/AliRoot.git] / DPMJET / flukapro / (NUCSFT)

*$ CREATE NUCSFT.ADD
*COPY NUCSFT
*
*=== nucsft ===========================================================*
*
*----------------------------------------------------------------------*
*                                                                      *
*     NUClear Symmetrized Fermi/woods-saxon density and potential      *
*     tabulations                                                      *
*                                                                      *
*     Last change   on   20-jun-94   by   Alfredo Ferrari, INFN-Milan  *
*                                                                      *
*          Vparvp, cparvp, aparvp (bparvp auxiliary parameter          *
*          for Cparvp calculations):                                   *
*                                                                      *
*          V0 (r)    = Vparvp * sinh (cparvp/aparvp)                   *
*                    / ( cosh (cparvp/aparvp) + cosh (r/aparvp) )      *
*          expvca    = exp ( -cparvp / aparvp )                        *
*          expvra    = exp ( rcursf / aparvp )                         *
*          bnvpsf    = "average" binding energy used when generating   *
*                      the potential tabulations                       *
*          v0cesf    = well depth at r=0 for the "average" potential   *
*          evsfmx    = Max. Fermi energy for the ip nucleon potential  *
*                      evsfmx = |v0cesf| - bnvpsf                      *
*          pvsfmx    = Max. Fermi momentum for proton/neutron          *
*          pfsksf    = Fermi momentum for proton/neutron computed at   *
*                      radiu1 defined ad the (mrsfbn-5)th bin limit    *
*          vcoumx    = maximum height of the Coulomb barrier           *
*                      (for Z=Zvpcou), modified during tracking to     *
*                      get the real barrier                            *
*          vcousv    = maximum height of the Coulomb barrier           *
*                      (for Z=Zvpcou)                                  *
*          v0ccou    = Coulomb potential at R=Radtot                   *
*          zvpcou    = Z for which the barrier has been computed       *
*          zbourd    = Z reduction factor to be used for bound nucleons*
*                      when computing Coulomb effects                  *
*          ravpsf    = radius at which the "average" nuclear potential *
*                      stops                                           *
*          ravcou    = radius at which the Coulomb potential stops     *
*          ravpot    = radius at which the nuclear potential stops     *
*          dradsf    = radial increment used from radiu0 up to radtot  *
*          dravsf    = radial increment used from radtot up to ravpot  *
*          ivsfmn    = index of the innermost potential shell which    *
*                      can be reached by the current particle accord-  *
*                      ing to its angular momentum                     *
*          ivsfmx    = index of the outermost potential shell which    *
*                      can be reached by the current particle accord-  *
*                      ing to its energy and angular momentum          *
*          ivcomx    = bin index of the maximum height of the Coulomb  *
*                      barrier                                         *
*          isfint    = current tabulation index                        *
*          ibsfmn    = index of the radial bin of the the minimum      *
*                      approach point                                  *
*          ibsfmx    = index of the radial bin of the the maximum      *
*                      approach point                                  *
*                                                                      *
*     The actual potential is given by (Preeqmvax version with no      *
*     energy dependence of the potential):                             *
*     for r < radtot (note V_nuc < 0):                                 *
*       neutral particles :                                            *
*          V_true(r) = V_nuc(r,ip) + bnvpsf (ip) - bn_true (ip)        *
*       charged particles :                                            *
*          V_true(r) = V_nuc(r,ip) + bnvpsf (ip) - bn_true (ip)        *
*                    + Zproj * V_coul(r)                               *
*     for ravpsf > r > radtot (note V_nuc < 0):                        *
*          V_true(r) = V_nuc(r,ip) / bnvpsf (ip) x bn_true (ip)        *
*                    + Znow / Zvpcou x Zproj x V_coul(r)               *
*     for r > ravpsf :                                                 *
*          V_true(r) = Znow / Zvpcou x Zproj x V_coul(r)               *
*     Please note that Ekfcen(ip)=Evsfmx(ip), for ip=1,2               *
*     For the actual potential with the Preeqmvax version with the     *
*     energy dependence of the potential, see the relevant routines    *
*                                                                      *
*----------------------------------------------------------------------*
*
      COMMON / NUCSFT / VPARVP (2), CPARVP (2), APARVP (2), BPARVP (2),
     &                  EXPVCA (2), EXPVRA (2), RAVPSF (2), BNVPSF (2),
     &                  V0CESF (2), EVSFMX (2), VSFHLP (2), PVSFMX (2),
     &                  PFSKSF (2), RAVCOU, RAVPOT, VCOUMX, V0CCOU,
     &                  ZVPCOU, ZBOURD, VPSFA0, VPSFA1, VPSFAE, DRADSF,
     &                  DRAVSF, VCOUSV,
     &                  IVSFMN, IVSFMX, IVCOMX, ISFINT, IBSFMN, IBSFMX
Commit	Line	Data
d30b8254	1	*$ CREATE NUCSFT.ADD
	2	*COPY NUCSFT
	3	*
	4	=== nucsft ===========================================================
	5	*
	6	----------------------------------------------------------------------
	7	* *
	8	* NUClear Symmetrized Fermi/woods-saxon density and potential *
	9	* tabulations *
	10	* *
	11	* Last change on 20-jun-94 by Alfredo Ferrari, INFN-Milan *
	12	* *
	13	* Vparvp, cparvp, aparvp (bparvp auxiliary parameter *
	14	* for Cparvp calculations): *
	15	* *
	16	* V0 (r) = Vparvp * sinh (cparvp/aparvp) *
	17	* / ( cosh (cparvp/aparvp) + cosh (r/aparvp) ) *
	18	* expvca = exp ( -cparvp / aparvp ) *
	19	* expvra = exp ( rcursf / aparvp ) *
	20	* bnvpsf = "average" binding energy used when generating *
	21	* the potential tabulations *
	22	* v0cesf = well depth at r=0 for the "average" potential *
	23	* evsfmx = Max. Fermi energy for the ip nucleon potential *
	24	* evsfmx = \|v0cesf\| - bnvpsf *
	25	* pvsfmx = Max. Fermi momentum for proton/neutron *
	26	* pfsksf = Fermi momentum for proton/neutron computed at *
	27	* radiu1 defined ad the (mrsfbn-5)th bin limit *
	28	* vcoumx = maximum height of the Coulomb barrier *
	29	* (for Z=Zvpcou), modified during tracking to *
	30	* get the real barrier *
	31	* vcousv = maximum height of the Coulomb barrier *
	32	* (for Z=Zvpcou) *
	33	* v0ccou = Coulomb potential at R=Radtot *
	34	* zvpcou = Z for which the barrier has been computed *
	35	* zbourd = Z reduction factor to be used for bound nucleons*
	36	* when computing Coulomb effects *
	37	* ravpsf = radius at which the "average" nuclear potential *
	38	* stops *
	39	* ravcou = radius at which the Coulomb potential stops *
	40	* ravpot = radius at which the nuclear potential stops *
	41	* dradsf = radial increment used from radiu0 up to radtot *
	42	* dravsf = radial increment used from radtot up to ravpot *
	43	* ivsfmn = index of the innermost potential shell which *
	44	* can be reached by the current particle accord- *
	45	* ing to its angular momentum *
	46	* ivsfmx = index of the outermost potential shell which *
	47	* can be reached by the current particle accord- *
	48	* ing to its energy and angular momentum *
	49	* ivcomx = bin index of the maximum height of the Coulomb *
	50	* barrier *
	51	* isfint = current tabulation index *
	52	* ibsfmn = index of the radial bin of the the minimum *
	53	* approach point *
	54	* ibsfmx = index of the radial bin of the the maximum *
	55	* approach point *
	56	* *
	57	* The actual potential is given by (Preeqmvax version with no *
	58	* energy dependence of the potential): *
	59	* for r < radtot (note V_nuc < 0): *
	60	* neutral particles : *
	61	* V_true(r) = V_nuc(r,ip) + bnvpsf (ip) - bn_true (ip) *
	62	* charged particles : *
	63	* V_true(r) = V_nuc(r,ip) + bnvpsf (ip) - bn_true (ip) *
	64	* + Zproj * V_coul(r) *
65	* for ravpsf > r > radtot (note V_nuc < 0): *
66	* V_true(r) = V_nuc(r,ip) / bnvpsf (ip) x bn_true (ip) *
67	* + Znow / Zvpcou x Zproj x V_coul(r) *
68	* for r > ravpsf : *
69	* V_true(r) = Znow / Zvpcou x Zproj x V_coul(r) *
70	* Please note that Ekfcen(ip)=Evsfmx(ip), for ip=1,2 *
71	* For the actual potential with the Preeqmvax version with the *
72	* energy dependence of the potential, see the relevant routines *
73	* *
74	----------------------------------------------------------------------
75	*
76	COMMON / NUCSFT / VPARVP (2), CPARVP (2), APARVP (2), BPARVP (2),
77	& EXPVCA (2), EXPVRA (2), RAVPSF (2), BNVPSF (2),
78	& V0CESF (2), EVSFMX (2), VSFHLP (2), PVSFMX (2),
79	& PFSKSF (2), RAVCOU, RAVPOT, VCOUMX, V0CCOU,
80	& ZVPCOU, ZBOURD, VPSFA0, VPSFA1, VPSFAE, DRADSF,
81	& DRAVSF, VCOUSV,
82	& IVSFMN, IVSFMX, IVCOMX, ISFINT, IBSFMN, IBSFMX
83