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

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