GEANT321/fluka/fpower.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:19:56  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.42  by  S.Giani
  11 *-- Author :
  12 *$ CREATE FPOWER.FOR
  13 *COPY FPOWER
  14 *
  15 *=== fpower ===========================================================*
  16 *
  17       FUNCTION FPOWER ( IPOWER, NU, ANUAV )
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 *
  23 *----------------------------------------------------------------------*
  24 *                                                                      *
  25 *     Created on   05  may  1992   by    Alfredo Ferrari & Paola Sala  *
  26 *                                                   Infn - Milan       *
  27 *                                                                      *
  28 *     Last change on 10-jun-92     by    Alfredo Ferrari               *
  29 *                                                                      *
  30 *                                                                      *
  31 *----------------------------------------------------------------------*
  32 *
  33       PARAMETER ( TWOLOG = 0.6931471805599453D+00 )
  34 *
  35       IF ( IPOWER .NE. 13 ) STOP 'FPOWER'
  36          IF ( ANUAV .LT. 2.D+00 ) THEN
  37             AKAPP0 = -0.93699278D+00*ANUAV**2 + 3.9082338D+00
  38          ELSE
  39             AKAPP0 = 4.2019D+00 / ANUAV**2.2835973D+00
  40          END IF
  41          IF ( NU .LE. 1 ) THEN
  42             FPOWER = AKAPP0
  43          ELSE
  44             AFACT  = ( NU - 1.D+00 ) / ( ANUAV + NU - 1.D+00 )
  45 *  |  This correction (1 card) has been put by A.F. on 29-12-92, trying
  46 *  |  to widen a bit the residual nuclei distribution without
  47 *  |  spoiling the good agreement with the grey particle distribution
  48             AFACT  = AFACT / SQRT (ANUAV)
  49             AKAPPA = ( 1.D+00 - AFACT ) * AKAPP0 + AFACT
  50             APOWER = ANUAV / NU + 1.D+00 * ANUAV / ( ANUAV + 1.D+00 )
  51             DDNU   = NU
  52             FPOWER = APOWER + LOG ( AKAPPA ) / LOG ( DDNU )
  53          END IF
  54 *=== End of function fpower ===========================================*
  55       RETURN
  56       END