]>
Commit | Line | Data |
---|---|---|
0795afa3 | 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 |