GEANT321/peanut/fekfnc.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:22:00  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.46  by  S.Giani
  11 *-- Author :
  12 *$ CREATE FEKFNC.FOR
  13 *COPY FEKFNC
  14 *
  15 *=== fekfnc ===========================================================*
  16 *
  17       FUNCTION FEKFNC (PF,I)
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 *
  23 #include "geant321/nucgeo.inc"
  24 #include "geant321/nucdat.inc"
  25 *
  26       IF ( PFRCEN (I) - PF .LE. ANGLGB ) THEN
  27          FEKFNC = EKFCEN (I)
  28       ELSE
  29          FEKFNC = SQRT ( AMNUSQ (I) + PF**2 ) - AMNUCL (I)
  30       END IF
  31       RETURN
  32       END