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