[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / fjet.F

* $Id$
C
C
C
	FUNCTION FJET(X,WGT)
	IMPLICIT REAL*8(A-H,O-Z)
#include "hiparnt.inc"
	DIMENSION X(10)
	SAVE
	PT2=(HINT1(1)**2/4.0-HIPR1(8)**2)*X(1)+HIPR1(8)**2
	XT=2.0*DSQRT(PT2)/HINT1(1)
	YMX1=DLOG(1.0/XT+DSQRT(1.0/XT**2-1.0))
	Y1=2.0*YMX1*X(2)-YMX1
	YMX2=DLOG(2.0/XT-DEXP(Y1))
	YMN2=DLOG(2.0/XT-DEXP(-Y1))
	Y2=(YMX2+YMN2)*X(3)-YMN2
	FJET=2.0*YMX1*(YMX2+YMN2)*(HINT1(1)**2/4.0-HIPR1(8)**2)
     &		*G(Y1,Y2,PT2)/2.0
	RETURN
	END
Commit	Line	Data
e74335a4	1	* $Id$
	2	C
	3	C
	4	C
	5	FUNCTION FJET(X,WGT)
	6	IMPLICIT REAL*8(A-H,O-Z)
	7	#include "hiparnt.inc"
	8	DIMENSION X(10)
	9	SAVE
	10	PT2=(HINT1(1)2/4.0-HIPR1(8)2)X(1)+HIPR1(8)*2
	11	XT=2.0*DSQRT(PT2)/HINT1(1)
	12	YMX1=DLOG(1.0/XT+DSQRT(1.0/XT**2-1.0))
	13	Y1=2.0YMX1X(2)-YMX1
	14	YMX2=DLOG(2.0/XT-DEXP(Y1))
	15	YMN2=DLOG(2.0/XT-DEXP(-Y1))
	16	Y2=(YMX2+YMN2)*X(3)-YMN2
	17	FJET=2.0YMX1(YMX2+YMN2)(HINT1(1)2/4.0-HIPR1(8)*2)
	18	& *G(Y1,Y2,PT2)/2.0
	19	RETURN
	20	END