From Philippe & Laurent: new variant of MUON visualization.

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

* $Id$
C
C
C
C
C	This generate random number between XMIN and XMAX
	FUNCTION HIRND2(I,XMIN,XMAX)
#define BLANKET_SAVE
#include "hijhb.inc"
	SAVE
	IF(XMIN.LT.XX(I,1)) XMIN=XX(I,1)
	IF(XMAX.GT.XX(I,201)) XMAX=XX(I,201)
	JMIN=1+200*(XMIN-XX(I,1))/(XX(I,201)-XX(I,1))
	JMAX=1+200*(XMAX-XX(I,1))/(XX(I,201)-XX(I,1))
	RX=RR(I,JMIN)+(RR(I,JMAX)-RR(I,JMIN))*RLU_HIJING(0)
	JL=0
	JU=202
10	IF(JU-JL.GT.1) THEN
	   JM=(JU+JL)/2
	   IF((RR(I,201).GT.RR(I,1)).EQV.(RX.GT.RR(I,JM))) THEN
	      JL=JM
	   ELSE
	      JU=JM
	   ENDIF
	GO TO 10
	ENDIF
	J=JL
	IF(J.LT.1) J=1
	IF(J.GE.201) J=200
	HIRND2=(XX(I,J)+XX(I,J+1))/2.0
	RETURN
	END