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