GEANT321/neutron/fisrnf.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:56  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/04 23/02/95  14.46.01  by  S.Giani
  11 *-- Author :
  12 C*********************************************************************
  13       FUNCTION FISRNF(A,B)
  14 C*********************************************************************
  15 C Sample secondary fission neutron energy from Watt spectrum
  16 C taken from ORNL/TM-7631
  17 C CZ 3/11/92
  18       DIMENSION RNDM(3)
  19 C
  20       CALL GRNDM(RNDM,3)
  21       Z=SQRT(-ALOG(RNDM(1)))
  22       S=6.28319*RNDM(2)
  23       ALOGR3=ALOG(RNDM(3))
  24       X=SQRT(A*B)/2.
  25       E1=A*((Z*COS(S)+X)**2-ALOGR3)
  26 C--  E2=A*((Z*SIN(S)+X)**2-ALOGR3)
  27 C distribution of E1 and E2 are identical
  28       FISRNF = E1
  29       RETURN
  30       END