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