HIJING/hijing1_36/g.F

   1 * $Id$
   2 C
   3 C
   4 C
   5 C
   6         FUNCTION G(Y1,Y2,PT2)
   7         IMPLICIT REAL*8 (A-H,O-Z)
   8 #define BLANKET_SAVE
   9 #include "hiparnt.inc"
  10         DIMENSION F(2,7)
  11         SAVE
  12         XT=2.0*DSQRT(PT2)/HINT1(1)
  13         X1=0.50*XT*(DEXP(Y1)+DEXP(Y2))
  14         X2=0.50*XT*(DEXP(-Y1)+DEXP(-Y2))
  15         DUMMY = 1.D0-XT**2/X1/X2
  16         if (DUMMY.LT.0.D0) DUMMY = 0.D0
  17         Z=DSQRT(DUMMY)
  18         SS=X1*X2*HINT1(1)**2
  19         T=-(1.0-Z)/2.0
  20         U=-(1.0+Z)/2.0
  21         AF=3.0
  22         DLAM=HIPR1(15)
  23         APH=12.0*3.1415926/(33.0-2.0*AF)/DLOG(PT2/DLAM**2)
  24 C
  25         CALL PARTON(F,X1,X2,PT2)
  26 C
  27         G11=( (F(1,1)+F(1,2))*(F(2,3)+F(2,4)+F(2,5)+F(2,6))
  28      &      +(F(1,3)+F(1,4))*(F(2,5)+F(2,6)) )*SUBCRS1(T,U)
  29 C
  30         G12=( (F(2,1)+F(2,2))*(F(1,3)+F(1,4)+F(1,5)+F(1,6))
  31      &      +(F(2,3)+F(2,4))*(F(1,5)+F(1,6)) )*SUBCRS1(U,T)
  32 C
  33         G13=(F(1,1)*F(2,1)+F(1,2)*F(2,2)+F(1,3)*F(2,3)+F(1,4)*F(2,4)
  34      &      +F(1,5)*F(2,5)+F(1,6)*F(2,6))*(SUBCRS1(U,T)
  35      &      +SUBCRS1(T,U)-8.D0/T/U/27.D0)
  36 C
  37         G2=(AF-1)*(F(1,1)*F(2,2)+F(2,1)*F(1,2)+F(1,3)*F(2,4)
  38      &     +F(2,3)*F(1,4)+F(1,5)*F(2,6)+F(2,5)*F(1,6))*SUBCRS2(T,U)
  39 C
  40         G31=(F(1,1)*F(2,2)+F(1,3)*F(2,4)+F(1,5)*F(2,6))*SUBCRS3(T,U)
  41         G32=(F(2,1)*F(1,2)+F(2,3)*F(1,4)+F(2,5)*F(1,6))*SUBCRS3(U,T)
  42 C
  43         G4=(F(1,1)*F(2,2)+F(2,1)*F(1,2)+F(1,3)*F(2,4)+F(2,3)*F(1,4)+
  44      1  F(1,5)*F(2,6)+F(2,5)*F(1,6))*SUBCRS4(T,U)
  45 C
  46         G5=AF*F(1,7)*F(2,7)*SUBCRS5(T,U)
  47 C
  48         G61=F(1,7)*(F(2,1)+F(2,2)+F(2,3)+F(2,4)+F(2,5)
  49      &      +F(2,6))*SUBCRS6(T,U)
  50         G62=F(2,7)*(F(1,1)+F(1,2)+F(1,3)+F(1,4)+F(1,5)
  51      &      +F(1,6))*SUBCRS6(U,T)
  52 C
  53         G7=F(1,7)*F(2,7)*SUBCRS7(T,U)
  54 C
  55         G=(G11+G12+G13+G2+G31+G32+G4+G5+G61+G62+G7)*HIPR1(17)*
  56      1  3.14159D0*APH**2/SS**2
  57         RETURN
  58         END