This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / GEANT321 / geant321 / nucstf.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_NUCSTF_INC
10 #define CERNLIB_GEANT321_NUCSTF_INC
11 *
12 *
13 * nucstf.inc
14 *
15 *
16 *=== nucstf ===========================================================*
17 *
18 *----------------------------------------------------------------------*
19 *                                                                      *
20 *     NUClear STatement Functions:                                     *
21 *                                                                      *
22 *     Created on  14 january 1992  by    Alfredo Ferrari & Paola Sala  *
23 *                                                   Infn - Milan       *
24 *                                                                      *
25 *     Last change on 04-feb-92     by    Alfredo Ferrari               *
26 *                                                                      *
27 *     Included in the following routines:                              *
28 *                                                                      *
29 *                        FPFRNC                                        *
30 *                        FRADNC                                        *
31 *                        FRHINC                                        *
32 *                        FRHONC                                        *
33 *                        SBCOMP                                        *
34 *                                                                      *
35 *----------------------------------------------------------------------*
36 *
37 *  NUClear STatement Functions: nucleus with a costant core and
38 *  linearly decreasing skin and halo
39 *  rho (r):
40 *     Statement functions:
41       RHCORE (R) = RHOCEN
42       RHSKIN (R) = RHOCEN * ( RADSKN - R ) / SKNEFF
43       RHHALO (R) = RHOSKN * ( RADTOT - R ) / HALODP
44 *  other
45 *     RHCORE (R) = RHOCEN
46 *  Skin with a linear and a cubic term:
47 *     RHSKIN (R) = RHOCEN * ( ( RADSKN - R ) / SKNEFF + SK3PAR * ( R
48 *    &           - CPARWS )**3 )
49 *  Skin with a linear, a cubic and a quartic term:
50 *     RHSKIN (R) = RHOCEN * ( ( RADSKN - R ) / SKNEFF
51 *    &           + ( R - CPARWS )**3 * ( SK3PAR + SK4PAR
52 *    &           * ( R - CPARWS ) ) )
53 *     RHHALO (R) = RHOCEN * HABPAR * ( RADTOT - R )**1.5D+00
54 * / r      2
55 * | 4 pi r' rho (r') dr':
56 * / 0
57 *     Statement functions:
58       RICORE (R) = 4.D+00 * PI / 3.D+00 * RHOCEN * R**3
59       RISKIN (R) = RICORE (RADIU0) + 4.D+00 * PI / 3.D+00 * RHOCEN
60      &           * ( RADSKN / SKNEFF * ( R**3 - RADIU0**3 )
61      &           - 0.75D+00 * ( R**4 - RADIU0**4 ) / SKNEFF )
62       RIHALO (R) = RISKIN (RADIU1) + 4.D+00 * PI / 3.D+00 * RHOSKN
63      &           * ( RADTOT / HALODP * ( R**3 - RADIU1**3 )
64      &           - 0.75D+00 * ( R**4 - RADIU1**4 ) / HALODP )
65 *  Other:
66 *     RICORE (R) = 4.D+00 * PI / 3.D+00 * RHOCEN * R**3
67 *  Skin with a linear and a cubic term:
68 *     RISKIN (R) = RICORE (RADIU0) + 4.D+00 * PI / 3.D+00 * RHOCEN
69 *    &           * ( ( RADSKN / SKNEFF - CPARWS**3 * SK3PAR )
70 *    &           * ( R**3 - RADIU0**3 ) - 0.75D+00 * ( 1.D+00 / SKNEFF
71 *    &           - 3.D+00 * SK3PAR * CPARWS**2 ) * ( R**4 - RADIU0**4 )
72 *    &           - 1.8D+00 * SK3PAR * CPARWS * ( R**5 - RADIU0**5 )
73 *    &           + 0.5D+00 * SK3PAR * ( R**6 - RADIU0**6 ) )
74 *  Skin with a linear, a cubic and a quartic term:
75 *     RISKIN (R) = RICORE (RADIU0) + 4.D+00 * PI / 3.D+00 * RHOCEN
76 *    &           * ( ( RADSKN / SKNEFF - CPARWS**3 * ( SK3PAR - CPARWS
77 *    &           * SK4PAR ) ) * ( R**3 - RADIU0**3 ) - 0.75D+00
78 *    &           * ( 1.D+00 / SKNEFF - CPARWS**2 * ( 3.D+00 * SK3PAR
79 *    &           - 4.D+00 * SK4PAR * CPARWS ) ) * ( R**4 - RADIU0**4 )
80 *    &           - 1.8D+00 * CPARWS * ( SK3PAR - 2.D+00 * SK4PAR
81 *    &           * CPARWS ) * ( R**5 - RADIU0**5 )
82 *    &           - ( 2.D+00 * SK4PAR * CPARWS - 0.5D+00 * SK3PAR )
83 *    &           * ( R**6 - RADIU0**6 ) + 3.D+00 / 7.D+00 * SK4PAR
84 *    &           * ( R**7 - RADIU0**7 ) )
85 *     RIHALO (R) = RISKIN (RADIU1) + 8.D+00 * PI * RHOCEN
86 *    &           * HABPAR * ( ( 0.2D+00 * RADTOT**2 - 2.D+00 / 7.D+00
87 *    &           * RADTOT * HALODP + HALODP**2 / 9.D+00 )
88 *    &           * ( SQRT (HALODP) )**5 - ( 0.2D+00 * RADTOT**2
89 *    &           - 2.D+00 / 7.D+00 * RADTOT * ( RADTOT - R )
90 *    &           + ( RADTOT - R )**2 / 9.D+00 )
91 *    &           * ( RADTOT - R )**2.5D+00 )
92 *  r (rho):
93 *     Statement functions:
94       RACORE (RHO) = RADIU0
95       RASKIN (RHO) = RADSKN - SKNEFF * RHO / RHOCEN
96       RAHALO (RHO) = RADTOT - HALODP * RHO / RHOSKN
97 *  p_f (rho,i), i=1 proton, 2 neutron:
98 *     Statement functions:
99       PFCORE (RHO,I) = PFRCEN (I)
100       PFSKIN (RHO,I) = PFRCEN (I) * (RHO/RHOCEN)**0.3333333333333333D+00
101       PFHALO (RHO,I) = PFRCEN (I) * (RHO/RHOCEN)**0.3333333333333333D+00
102 *  end NUClear STatement Functions
103  
104
105 #endif