GEANT321/gphys/glandg.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:25  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.21  by  S.Giani
  11 *-- Author :
  12       SUBROUTINE GLANDG(YRAN)
  13 C.
  14 C.    ******************************************************************
  15 C.    *                                                                *
  16 C.    *  Copy of the CERN library routine GENLAN                       *
  17 C.    *  Generation of LANDAU-distributed random numbers by 4-point    *
  18 C.    *  interpolation in the previously-tabulated inverse cumulative  *
  19 C.    *  distribution                                                  *
  20 C.    *                                                                *
  21 C.    *    ==>Called by : GLANDO                                       *
  22 C.    *                                                                *
  23 C.    ******************************************************************
  24 C.
  25       PARAMETER (RANGE1=0.807069,RANGE2=0.994869)
  26       PARAMETER (X1BOT =0.0     ,X2BOT =0.791240)
  27       PARAMETER (GP1INV=118.836 ,GP2INV=486.178 )
  28 C
  29       DIMENSION XCUM1(100) , XCUM2(100)
  30       DIMENSION RNDM(1)
  31 C
  32 C====>  1ST TABLE OF INVERSE CUMULATIVE LANDAU POINTS. ( 0.0<P<0.833 )
  33       DATA XCUM1(  1) / -2.57000 /   ,   XCUM1(  2) / -2.15412 /
  34       DATA XCUM1(  3) / -1.94167 /   ,   XCUM1(  4) / -1.79583 /
  35       DATA XCUM1(  5) / -1.67975 /   ,   XCUM1(  6) / -1.58090 /
  36       DATA XCUM1(  7) / -1.49341 /   ,   XCUM1(  8) / -1.41397 /
  37       DATA XCUM1(  9) / -1.34057 /   ,   XCUM1( 10) / -1.27185 /
  38       DATA XCUM1( 11) / -1.20686 /   ,   XCUM1( 12) / -1.14490 /
  39       DATA XCUM1( 13) / -1.08545 /   ,   XCUM1( 14) / -1.02809 /
  40       DATA XCUM1( 15) /  -.97249 /   ,   XCUM1( 16) /  -.91839 /
  41       DATA XCUM1( 17) /  -.86557 /   ,   XCUM1( 18) /  -.81386 /
  42       DATA XCUM1( 19) /  -.76308 /   ,   XCUM1( 20) /  -.71311 /
  43       DATA XCUM1( 21) /  -.66384 /   ,   XCUM1( 22) /  -.61515 /
  44       DATA XCUM1( 23) /  -.56696 /   ,   XCUM1( 24) /  -.51919 /
  45       DATA XCUM1( 25) /  -.47176 /   ,   XCUM1( 26) /  -.42462 /
  46       DATA XCUM1( 27) /  -.37769 /   ,   XCUM1( 28) /  -.33092 /
  47       DATA XCUM1( 29) /  -.28426 /   ,   XCUM1( 30) /  -.23767 /
  48       DATA XCUM1( 31) /  -.19109 /   ,   XCUM1( 32) /  -.14447 /
  49       DATA XCUM1( 33) /  -.09779 /   ,   XCUM1( 34) /  -.05100 /
  50       DATA XCUM1( 35) /  -.00405 /   ,   XCUM1( 36) /   .04308 /
  51       DATA XCUM1( 37) /   .09044 /   ,   XCUM1( 38) /   .13806 /
  52       DATA XCUM1( 39) /   .18598 /   ,   XCUM1( 40) /   .23423 /
  53       DATA XCUM1( 41) /   .28285 /   ,   XCUM1( 42) /   .33187 /
  54       DATA XCUM1( 43) /   .38134 /   ,   XCUM1( 44) /   .43128 /
  55       DATA XCUM1( 45) /   .48174 /   ,   XCUM1( 46) /   .53275 /
  56       DATA XCUM1( 47) /   .58435 /   ,   XCUM1( 48) /   .63658 /
  57       DATA XCUM1( 49) /   .68948 /   ,   XCUM1( 50) /   .74309 /
  58       DATA XCUM1( 51) /   .79746 /   ,   XCUM1( 52) /   .85263 /
  59       DATA XCUM1( 53) /   .90865 /   ,   XCUM1( 54) /   .96556 /
  60       DATA XCUM1( 55) /  1.02342 /   ,   XCUM1( 56) /  1.08228 /
  61       DATA XCUM1( 57) /  1.14219 /   ,   XCUM1( 58) /  1.20322 /
  62       DATA XCUM1( 59) /  1.26542 /   ,   XCUM1( 60) /  1.32887 /
  63       DATA XCUM1( 61) /  1.39362 /   ,   XCUM1( 62) /  1.45976 /
  64       DATA XCUM1( 63) /  1.52736 /   ,   XCUM1( 64) /  1.59650 /
  65       DATA XCUM1( 65) /  1.66727 /   ,   XCUM1( 66) /  1.73976 /
  66       DATA XCUM1( 67) /  1.81407 /   ,   XCUM1( 68) /  1.89032 /
  67       DATA XCUM1( 69) /  1.96860 /   ,   XCUM1( 70) /  2.04905 /
  68       DATA XCUM1( 71) /  2.13180 /   ,   XCUM1( 72) /  2.21699 /
  69       DATA XCUM1( 73) /  2.30477 /   ,   XCUM1( 74) /  2.39531 /
  70       DATA XCUM1( 75) /  2.48878 /   ,   XCUM1( 76) /  2.58540 /
  71       DATA XCUM1( 77) /  2.68535 /   ,   XCUM1( 78) /  2.78889 /
  72       DATA XCUM1( 79) /  2.89626 /   ,   XCUM1( 80) /  3.00775 /
  73       DATA XCUM1( 81) /  3.12364 /   ,   XCUM1( 82) /  3.24429 /
  74       DATA XCUM1( 83) /  3.37005 /   ,   XCUM1( 84) /  3.50136 /
  75       DATA XCUM1( 85) /  3.63866 /   ,   XCUM1( 86) /  3.78246 /
  76       DATA XCUM1( 87) /  3.93334 /   ,   XCUM1( 88) /  4.09194 /
  77       DATA XCUM1( 89) /  4.25900 /   ,   XCUM1( 90) /  4.43533 /
  78       DATA XCUM1( 91) /  4.62186 /   ,   XCUM1( 92) /  4.81960 /
  79       DATA XCUM1( 93) /  5.02974 /   ,   XCUM1( 94) /  5.25368 /
  80       DATA XCUM1( 95) /  5.49312 /   ,   XCUM1( 96) /  5.74987 /
  81       DATA XCUM1( 97) /  6.02605 /   ,   XCUM1( 98) /  6.32428 /
  82       DATA XCUM1( 99) /  6.64773 /   ,   XCUM1(100) /  7.00000 /
  83 C
  84 C====>  2ND TABLE OF INVERSE CUMULATIVE LANDAU POINTS. ( 0.791<P<0.995 )
  85       DATA XCUM2(  1) /   5.50000 /   ,   XCUM2(  2) /   5.56120 /
  86       DATA XCUM2(  3) /   5.62347 /   ,   XCUM2(  4) /   5.68684 /
  87       DATA XCUM2(  5) /   5.75133 /   ,   XCUM2(  6) /   5.81699 /
  88       DATA XCUM2(  7) /   5.88383 /   ,   XCUM2(  8) /   5.95191 /
  89       DATA XCUM2(  9) /   6.02125 /   ,   XCUM2( 10) /   6.09190 /
  90       DATA XCUM2( 11) /   6.16391 /   ,   XCUM2( 12) /   6.23732 /
  91       DATA XCUM2( 13) /   6.31219 /   ,   XCUM2( 14) /   6.38855 /
  92       DATA XCUM2( 15) /   6.46646 /   ,   XCUM2( 16) /   6.54597 /
  93       DATA XCUM2( 17) /   6.62712 /   ,   XCUM2( 18) /   6.70998 /
  94       DATA XCUM2( 19) /   6.79460 /   ,   XCUM2( 20) /   6.88103 /
  95       DATA XCUM2( 21) /   6.96935 /   ,   XCUM2( 22) /   7.05962 /
  96       DATA XCUM2( 23) /   7.15192 /   ,   XCUM2( 24) /   7.24633 /
  97       DATA XCUM2( 25) /   7.34294 /   ,   XCUM2( 26) /   7.44182 /
  98       DATA XCUM2( 27) /   7.54306 /   ,   XCUM2( 28) /   7.64676 /
  99       DATA XCUM2( 29) /   7.75300 /   ,   XCUM2( 30) /   7.86188 /
 100       DATA XCUM2( 31) /   7.97351 /   ,   XCUM2( 32) /   8.08800 /
 101       DATA XCUM2( 33) /   8.20548 /   ,   XCUM2( 34) /   8.32610 /
 102       DATA XCUM2( 35) /   8.44997 /   ,   XCUM2( 36) /   8.57725 /
 103       DATA XCUM2( 37) /   8.70808 /   ,   XCUM2( 38) /   8.84262 /
 104       DATA XCUM2( 39) /   8.98103 /   ,   XCUM2( 40) /   9.12349 /
 105       DATA XCUM2( 41) /   9.27021 /   ,   XCUM2( 42) /   9.42141 /
 106       DATA XCUM2( 43) /   9.57730 /   ,   XCUM2( 44) /   9.73812 /
 107       DATA XCUM2( 45) /   9.90410 /   ,   XCUM2( 46) /  10.07552 /
 108       DATA XCUM2( 47) /  10.25265 /   ,   XCUM2( 48) /  10.43584 /
 109       DATA XCUM2( 49) /  10.62540 /   ,   XCUM2( 50) /  10.82169 /
 110       DATA XCUM2( 51) /  11.02508 /   ,   XCUM2( 52) /  11.23598 /
 111       DATA XCUM2( 53) /  11.45487 /   ,   XCUM2( 54) /  11.68222 /
 112       DATA XCUM2( 55) /  11.91855 /   ,   XCUM2( 56) /  12.16441 /
 113       DATA XCUM2( 57) /  12.42045 /   ,   XCUM2( 58) /  12.68734 /
 114       DATA XCUM2( 59) /  12.96580 /   ,   XCUM2( 60) /  13.25663 /
 115       DATA XCUM2( 61) /  13.56075 /   ,   XCUM2( 62) /  13.87912 /
 116       DATA XCUM2( 63) /  14.21280 /   ,   XCUM2( 64) /  14.56296 /
 117       DATA XCUM2( 65) /  14.93088 /   ,   XCUM2( 66) /  15.31799 /
 118       DATA XCUM2( 67) /  15.72593 /   ,   XCUM2( 68) /  16.15650 /
 119       DATA XCUM2( 69) /  16.61170 /   ,   XCUM2( 70) /  17.09378 /
 120       DATA XCUM2( 71) /  17.60523 /   ,   XCUM2( 72) /  18.14884 /
 121       DATA XCUM2( 73) /  18.72792 /   ,   XCUM2( 74) /  19.34622 /
 122       DATA XCUM2( 75) /  20.00796 /   ,   XCUM2( 76) /  20.71792 /
 123       DATA XCUM2( 77) /  21.48176 /   ,   XCUM2( 78) /  22.30618 /
 124       DATA XCUM2( 79) /  23.19862 /   ,   XCUM2( 80) /  24.16809 /
 125       DATA XCUM2( 81) /  25.22547 /   ,   XCUM2( 82) /  26.38320 /
 126       DATA XCUM2( 83) /  27.65658 /   ,   XCUM2( 84) /  29.06490 /
 127       DATA XCUM2( 85) /  30.63048 /   ,   XCUM2( 86) /  32.38211 /
 128       DATA XCUM2( 87) /  34.35555 /   ,   XCUM2( 88) /  36.59607 /
 129       DATA XCUM2( 89) /  39.16212 /   ,   XCUM2( 90) /  42.13195 /
 130       DATA XCUM2( 91) /  45.61204 /   ,   XCUM2( 92) /  49.74758 /
 131       DATA XCUM2( 93) /  54.74189 /   ,   XCUM2( 94) /  60.89935 /
 132       DATA XCUM2( 95) /  68.67370 /   ,   XCUM2( 96) /  78.81480 /
 133       DATA XCUM2( 97) /  92.61047 /   ,   XCUM2( 98) / 112.50807 /
 134       DATA XCUM2( 99) / 143.78539 /   ,   XCUM2(100) / 200.00000 /
 135 C.
 136 C.    ------------------------------------------------------------------
 137 C.
 138       CALL GRNDM(RNDM,1)
 139       X = RNDM(1)
 140       IF(X .LT. 0.004) THEN
 141 *
 142 * *** Extreme left-hand tail
 143          YRAN = -SQRT(ABS(LOG(X)))
 144 *
 145       ELSEIF(X.LE.RANGE1) THEN
 146 *
 147 * *** 4-point interpolation in the first cumulative table
 148          TABPO1 = (X-X1BOT)*GP1INV
 149          J = TABPO1 + 1
 150          J = MAX(J,2)
 151          J = MIN(J,98)
 152          P = TABPO1 - (J-1)
 153          A = (P+1.0) * XCUM1(J+2) - (P-2.0) * XCUM1(J-1)
 154          B = (P-1.0) * XCUM1(J)   - P       * XCUM1(J+1)
 155          YRAN = A*P*(P-1.0)*0.16666667+B*(P+1.0)*(P-2.0)*0.5
 156 *
 157       ELSEIF(X.LE.RANGE2) THEN
 158 *
 159 * *** 4-point interpolation in the first cumulative table
 160          TABPO2 = (X-X2BOT)*GP2INV
 161          J = TABPO2 + 1
 162          J = MAX(J,2)
 163          J = MIN(J,98)
 164          P = TABPO2 - (J-1)
 165          A = (P+1.0) * XCUM2(J+2) - (P-2.0) * XCUM2(J-1)
 166          B = (P-1.0) * XCUM2(J)   - P       * XCUM2(J+1)
 167          YRAN = A*P*(P-1.0)*0.16666667+B*(P+1.0)*(P-2.0)*0.5
 168 *
 169       ELSE
 170 *
 171 * *** 1/x**2 sampling for extreme Landau tail
 172          CALL GRNDM(RNDM,1)
 173          YRAN=200./RNDM(1)
 174 *
 175       ENDIF
 176       END