hirnd2.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file hirnd2.f

C
C
C
C
C This generate random number between XMIN and XMAX
  FUNCTION hirnd2(I,XMIN,XMAX)
  common/hijhb/rr(10,201),xx(10,201)
  SAVE  /hijhb/
  common/ranseed/nseed
  SAVE  /ranseed/
  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))*atl_ran(nseed)
  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