HIJING/hijing1_36/hifun.F

   1 * $Id$
   2 C
   3 C
   4 C The next three subroutines are for Monte Carlo generation
   5 C according to a given function FHB. One calls first HIFUN
   6 C with assigned channel number I, low and up limits. Then to
   7 C generate the distribution one can call HIRND(I) which gives
   8 C you a random number generated according to the given function.
   9 C
  10         SUBROUTINE HIFUN(I,XMIN,XMAX,FHB)
  11 #include "hijhb.inc"
  12         EXTERNAL FHB
  13         FNORM=GAUSS1(FHB,XMIN,XMAX,0.001)
  14         DO 100 J=1,201
  15                 XX(I,J)=XMIN+(XMAX-XMIN)*(J-1)/200.0
  16                 XDD=XX(I,J)
  17                 RR(I,J)=GAUSS1(FHB,XMIN,XDD,0.001)/FNORM
  18 100     CONTINUE
  19         RETURN
  20         END