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