hifun.f

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

C
C
C The next three subroutines are for Monte Carlo generation 
C according to a given function FHB. One calls first HIFUN 
C with assigned channel number I, low and up limits. Then to 
C generate the distribution one can call HIRND(I) which gives 
C you a random number generated according to the given function.
C 
        SUBROUTINE hifun(I,XMIN,XMAX,FHB)
        common/hijhb/rr(10,201),xx(10,201)
        SAVE  /hijhb/
        EXTERNAL fhb
        fnorm=gauss1(fhb,xmin,xmax,0.001)
        DO 100 j=1,201
                xx(i,j)=xmin+(xmax-xmin)*(j-1)/200.0
                xdd=xx(i,j)
                rr(i,j)=gauss1(fhb,xmin,xdd,0.001)/fnorm
100     CONTINUE
        RETURN
        END