[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)
	IF (XT.GT.1.0) XT=1.0 
	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)
c87d167b	13	IF (XT.GT.1.0) XT=1.0
e74335a4	14	YMX1=DLOG(1.0/XT+DSQRT(1.0/XT**2-1.0))
	15	Y1=2.0YMX1X(2)-YMX1
	16	YMX2=DLOG(2.0/XT-DEXP(Y1))
	17	YMN2=DLOG(2.0/XT-DEXP(-Y1))
	18	Y2=(YMX2+YMN2)*X(3)-YMN2
	19	FJET=2.0YMX1(YMX2+YMN2)(HINT1(1)2/4.0-HIPR1(8)*2)
	20	& *G(Y1,Y2,PT2)/2.0
	21	RETURN
	22	END