More volume overlaps corrected
[u/mrichter/AliRoot.git] / ISAJET / code / getpt.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 SUBROUTINE GETPT(PT0,PTMEAN)
3C GENERATE PT WITH 1/(1+B*PT**2)**4 DISTRIBUTION
4C (APPROXIMATELY AN EXPONENTIAL FOR PT < 2 GEV.)
5C CON1=16/(3*PI)
6C 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