fjet.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file fjet.f

C
C
C
        FUNCTION fjet(X,WGT)
        IMPLICIT REAL*8(a-h,o-z)
        REAL hipr1(100),hint1(100)
        common/hiparnt/hipr1,ihpr2(50),hint1,ihnt2(50)
        SAVE  /hiparnt/
        dimension x(10)
        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