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