GEANT321/neutron/rnmaxf.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:21:58  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/04 23/02/95  14.46.01  by  S.Giani
  11 *-- Author :
  12 C*********************************************************************
  13       FUNCTION RNMAXF(T)
  14 C T := most probable value of distribution
  15 C*********************************************************************
  16       DATA FF/0./
  17       SAVE FF,R1SQ,W,U
  18       U=EXPRNF(U)
  19       IF(FF) 30 ,10 ,30
  20    10 R1=FLTRNF(R1)
  21       R2=FLTRNF(R2)
  22       R1SQ=R1*R1
  23       R2SQ=R2*R2
  24       RSQ=R1SQ+R2SQ
  25       IF(RSQ-1.) 20 ,20 ,10
  26    20 W=EXPRNF(W)/RSQ
  27       FF=1.
  28       RNMAXF=(R2SQ*W+U)*T
  29       GO TO 40
  30    30 FF=0.
  31       RNMAXF=(R1SQ*W+U)*T
  32    40 RETURN
  33       END