[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / hifun.F

* $Id$
C
C
C The next three subroutines are for Monte Carlo generation 
C according to a given function FHB. One calls first HIFUN 
C with assigned channel number I, low and up limits. Then to 
C generate the distribution one can call HIRND(I) which gives 
C you a random number generated according to the given function.
C 
	SUBROUTINE HIFUN(I,XMIN,XMAX,FHB)
#include "hijhb.inc"
	EXTERNAL FHB
	FNORM=GAUSS1(FHB,XMIN,XMAX,0.001)
	DO 100 J=1,201
		XX(I,J)=XMIN+(XMAX-XMIN)*(J-1)/200.0
		XDD=XX(I,J)
		RR(I,J)=GAUSS1(FHB,XMIN,XDD,0.001)/FNORM
100	CONTINUE
	RETURN
	END
Commit	Line	Data
e74335a4	1	* $Id$
	2	C
	3	C
	4	C The next three subroutines are for Monte Carlo generation
	5	C according to a given function FHB. One calls first HIFUN
	6	C with assigned channel number I, low and up limits. Then to
	7	C generate the distribution one can call HIRND(I) which gives
	8	C you a random number generated according to the given function.
	9	C
	10	SUBROUTINE HIFUN(I,XMIN,XMAX,FHB)
	11	#include "hijhb.inc"
	12	EXTERNAL FHB
	13	FNORM=GAUSS1(FHB,XMIN,XMAX,0.001)
	14	DO 100 J=1,201
	15	XX(I,J)=XMIN+(XMAX-XMIN)*(J-1)/200.0
	16	XDD=XX(I,J)
	17	RR(I,J)=GAUSS1(FHB,XMIN,XDD,0.001)/FNORM
	18	100 CONTINUE
	19	RETURN
	20	END