3 C*********************************************************************
5 SUBROUTINE LUFOWO_HIJING(H10,H20,H30,H40)
7 C...Purpose: to calculate the first few Fox-Wolfram moments.
8 #include "lujets_hijing.inc"
9 #include "ludat1_hijing.inc"
10 #include "ludat2_hijing.inc"
12 C...Copy momenta for particles and calculate H0.
17 IF(K(I,1).LE.0.OR.K(I,1).GT.10) GOTO 110
18 IF(MSTU(41).GE.2) THEN
19 KC=LUCOMP_HIJING(K(I,2))
20 IF(KC.EQ.0.OR.KC.EQ.12.OR.KC.EQ.14.OR.KC.EQ.16.OR.
22 IF(MSTU(41).GE.3.AND.KCHG(KC,2).EQ.0.AND.LUCHGE_HIJING(K(I,2))
25 IF(N+NP.GE.MSTU(4)-MSTU(32)-5) THEN
27 $ ,'(LUFOWO_HIJING:) no more memory left in LUJETS_HIJING')
37 P(N+NP,4)=SQRT(P(I,1)**2+P(I,2)**2+P(I,3)**2)
43 C...Very low multiplicities (0 or 1) not considered.
46 $ ,'(LUFOWO_HIJING:) too few particles for analysis')
54 C...Calculate H1 - H4.
61 CTHE=(P(I1,1)*P(I2,1)+P(I1,2)*P(I2,2)+P(I1,3)*P(I2,3))/
63 H10=H10+P(I1,4)*P(I2,4)*CTHE
64 H20=H20+P(I1,4)*P(I2,4)*(1.5*CTHE**2-0.5)
65 H30=H30+P(I1,4)*P(I2,4)*(2.5*CTHE**3-1.5*CTHE)
66 H40=H40+P(I1,4)*P(I2,4)*(4.375*CTHE**4-3.75*CTHE**2+0.375)
69 C...Calculate H1/H0 - H4/H0. Output.