[u/mrichter/AliRoot.git] / ISAJET / code / getpt.F

#include "isajet/pilot.h"
      SUBROUTINE GETPT(PT0,PTMEAN)
C          GENERATE PT WITH 1/(1+B*PT**2)**4 DISTRIBUTION
C          (APPROXIMATELY AN EXPONENTIAL FOR PT < 2 GEV.)
C          CON1=16/(3*PI)
C          CON2=-1/3
#include "isajet/itapes.inc"
      DATA CON1/1.697652726/,CON2/-.3333333333/
      R=RANF()
      ARG=AMAX1(R**CON2-1.,0.)
      PT0=PTMEAN*CON1*SQRT(ARG)
      RETURN
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	SUBROUTINE GETPT(PT0,PTMEAN)
	3	C GENERATE PT WITH 1/(1+BPT2)*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=PTMEANCON1SQRT(ARG)
	12	RETURN
	13	END