This commit was generated by cvs2svn to compensate for changes in r7641,
[u/mrichter/AliRoot.git] / DPMJET / flukapro / (NUCSFT)
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