Adding macros to create Calibration objects
[u/mrichter/AliRoot.git] / HIJING / hijing1_36 / g.F
CommitLineData
e74335a4 1* $Id$
2C
3C
4C
5C
6 FUNCTION G(Y1,Y2,PT2)
7 IMPLICIT REAL*8 (A-H,O-Z)
bc676b8e 8#define BLANKET_SAVE
e74335a4 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))
9f320b1f 15 DUMMY = 1.D0-XT**2/X1/X2
16 if (DUMMY.LT.0.D0) DUMMY = 0.D0
17 Z=DSQRT(DUMMY)
e74335a4 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)
24C
25 CALL PARTON(F,X1,X2,PT2)
26C
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)
29C
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)
32C
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)
36C
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)
39C
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)
42C
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)
45C
46 G5=AF*F(1,7)*F(2,7)*SUBCRS5(T,U)
47C
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)
52C
53 G7=F(1,7)*F(2,7)*SUBCRS7(T,U)
54C
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