[u/mrichter/AliRoot.git] / GEANT321 / neutron / fisrnf.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:21:56  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/04 23/02/95  14.46.01  by  S.Giani
*-- Author :
C*********************************************************************
      FUNCTION FISRNF(A,B)
C*********************************************************************
C Sample secondary fission neutron energy from Watt spectrum
C taken from ORNL/TM-7631
C CZ 3/11/92
      DIMENSION RNDM(3)
C
      CALL GRNDM(RNDM,3)
      Z=SQRT(-ALOG(RNDM(1)))
      S=6.28319*RNDM(2)
      ALOGR3=ALOG(RNDM(3))
      X=SQRT(A*B)/2.
      E1=A*((Z*COS(S)+X)**2-ALOGR3)
C--  E2=A*((Z*SIN(S)+X)**2-ALOGR3)
C distribution of E1 and E2 are identical
      FISRNF = E1
      RETURN
      END
Commit	Line	Data
fe4da5cc	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((ZCOS(S)+X)**2-ALOGR3)
	26	C-- E2=A((ZSIN(S)+X)**2-ALOGR3)
	27	C distribution of E1 and E2 are identical
	28	FISRNF = E1
	29	RETURN
	30	END