This commit was generated by cvs2svn to compensate for changes in r7641,
[u/mrichter/AliRoot.git] / DPMJET / flukapro / (NUCSFX)
1 *$ CREATE NUCSFX.ADD
2 *COPY NUCSFX
3 *
4 *=== nucsfx ===========================================================*
5 *
6 *----------------------------------------------------------------------*
7 *                                                                      *
8 *     NUClear Symmetrized Fermi/woods-saXon                            *
9 *                                                                      *
10 *          rh0sfx, cparsf, aparsf:                                     *
11 *                                                                      *
12 *          rho (r) = rh0sfx * sinh (cparfx/aparfx)                     *
13 *                  / ( cosh (cparfx/aparfx) + cosh (r/aparfx) )        *
14 *                                                                      *
15 *          radsfx(i) = max ( 0, asfhlf x (nsfxbf+i) + cparsf )         *
16 *          rhosfx(i) = rho (radsfx(i))                                 *
17 *          rhpsfx(i) = d rho / dr |radsfx(i)                           *
18 *                      / radsfx(i)                                     *
19 *          rhisfx(i) = | 4 pi r^2 rho dr                               *
20 *                      / radsfx (i-1)                                  *
21 *          exisfx(i) = exp ((radsfx(i)-cparsf)/asfhlf) =               *
22 *                      exp((nsfxbf+i)/2)                               *
23 *          asfhlf    = aparsf / 2                                      *
24 *          expcoa    = exp ( -cparsf / aparsf )                        *
25 *          exproa    = exp ( rcursf / aparsf )                         *
26 *          x,y,zcursf= current position                                *
27 *          u,v,wcursf= current direction                               *
28 *          ecursf    = current (total) energy                          *
29 *          pcursf    = current momentum                                *
30 *          vpcrsf    = current potential                               *
31 *          rcursf    = current radius                                  *
32 *          lingng    = logical flag to signal if the current particle  *
33 *                      is approaching the nucleus centre               *
34 *                                                                      *
35 *----------------------------------------------------------------------*
36 *
37       PARAMETER ( NINSFX = +21 )
38       PARAMETER ( NSFXBF = -11 )
39       LOGICAL LINGNG
40       COMMON / NUCSFX / RHOSFX (0:NINSFX), RHPSFX (0:NINSFX),
41      &                  RADSFX (0:NINSFX), RHISFX (NINSFX),
42      &                  EXISFX (0:NINSFX), RH0SFX, CPARSF, APARSF,
43      &                  ASFHLF, EXPROA, XCURSF, YCURSF, ZCURSF, UCURSF,
44      &                  VCURSF, WCURSF, ECURSF, PCURSF, VPCRSF, RCURSF,
45      &                  RHOCUR, RPCRSF, R13CUR, ISFX  , LINGNG
46       EQUIVALENCE ( EXPCOA, EXISFX (0) )
47