hirnd.f

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

C
C
C
  FUNCTION hirnd(I)
  common/hijhb/rr(10,201),xx(10,201)
  SAVE  /hijhb/
  common/ranseed/nseed
  SAVE  /ranseed/
  rx=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
  hirnd=(xx(i,j)+xx(i,j+1))/2.0
  RETURN
  END