ISAJET/code/getpt.F

   1 #include "isajet/pilot.h"
   2       SUBROUTINE GETPT(PT0,PTMEAN)
   3 C          GENERATE PT WITH 1/(1+B*PT**2)**4 DISTRIBUTION
   4 C          (APPROXIMATELY AN EXPONENTIAL FOR PT < 2 GEV.)
   5 C          CON1=16/(3*PI)
   6 C          CON2=-1/3
   7 #include "isajet/itapes.inc"
   8       DATA CON1/1.697652726/,CON2/-.3333333333/
   9       R=RANF()
  10       ARG=AMAX1(R**CON2-1.,0.)
  11       PT0=PTMEAN*CON1*SQRT(ARG)
  12       RETURN
  13       END