2 C*********************************************************************
4 SUBROUTINE LUFOWO(H10,H20,H30,H40)
6 C...Purpose: to calculate the first few Fox-Wolfram moments.
7 COMMON/LUJETS/N,K(4000,5),P(4000,5),V(4000,5)
8 COMMON/LUDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
9 COMMON/LUDAT2/KCHG(500,3),PMAS(500,4),PARF(2000),VCKM(4,4)
10 SAVE /LUJETS/,/LUDAT1/,/LUDAT2/
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
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(K(I,2)).EQ.0)
25 IF(N+NP.GE.MSTU(4)-MSTU(32)-5) THEN
26 CALL LUERRM(11,'(LUFOWO:) no more memory left in LUJETS')
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.
45 CALL LUERRM(8,'(LUFOWO:) too few particles for analysis')
53 C...Calculate H1 - H4.
60 CTHE=(P(I1,1)*P(I2,1)+P(I1,2)*P(I2,2)+P(I1,3)*P(I2,3))/
62 H10=H10+P(I1,4)*P(I2,4)*CTHE
63 H20=H20+P(I1,4)*P(I2,4)*(1.5*CTHE**2-0.5)
64 H30=H30+P(I1,4)*P(I2,4)*(2.5*CTHE**3-1.5*CTHE)
65 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.