[u/mrichter/AliRoot.git] / HERWIG / src / hwrgau.f


CDECK  ID>, HWRGAU.

*CMZ :-        -19/05/99  11.11.56  by  Mike Seymour

*-- Author :    Mike Seymour

C-----------------------------------------------------------------------

      FUNCTION HWRGAU(J,A,B)

C-----------------------------------------------------------------------

C     Gaussian random number, mean A, standard deviation B.

C     Generates uncorrelated pairs and throws one of them away.

C-----------------------------------------------------------------------

      INCLUDE 'HERWIG61.INC'

      DOUBLE PRECISION HWRGAU,HWR,A,B,X,TRASH

      INTEGER J

      EXTERNAL HWR

 10   X=HWR()

      IF (X.LE.ZERO.OR.X.GT.ONE) GOTO 10

      X=SQRT(-TWO*LOG(X))

      CALL HWRAZM(X,X,TRASH)

      HWRGAU=A+B*X

      END
Commit	Line	Data
3820ca8e	1
	2	CDECK ID>, HWRGAU.
	3
	4	*CMZ :- -19/05/99 11.11.56 by Mike Seymour
	5
	6	*-- Author : Mike Seymour
	7
	8	C-----------------------------------------------------------------------
	9
	10	FUNCTION HWRGAU(J,A,B)
	11
	12	C-----------------------------------------------------------------------
	13
	14	C Gaussian random number, mean A, standard deviation B.
	15
	16	C Generates uncorrelated pairs and throws one of them away.
	17
	18	C-----------------------------------------------------------------------
	19
	20	INCLUDE 'HERWIG61.INC'
	21
	22	DOUBLE PRECISION HWRGAU,HWR,A,B,X,TRASH
	23
	24	INTEGER J
	25
	26	EXTERNAL HWR
	27
	28	10 X=HWR()
	29
	30	IF (X.LE.ZERO.OR.X.GT.ONE) GOTO 10
	31
	32	X=SQRT(-TWO*LOG(X))
	33
	34	CALL HWRAZM(X,X,TRASH)
	35
	36	HWRGAU=A+B*X
	37
	38	END