More volume overlaps corrected
[u/mrichter/AliRoot.git] / ISAJET / code / wztu.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 FUNCTION WZTU(T,U,T1,U1,T3,U3,P1,P2)
3C DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980).
4C ALL T AND U GRAPH TERMS.
5#include "isajet/itapes.inc"
6#include "isajet/wwpar.inc"
7 DIMENSION P1(4),P2(4)
8#if defined(CERNLIB_DOUBLE)
9 DOUBLE PRECISION WZTU
10 DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
11 DOUBLE PRECISION WM4,ZM4,WZM2,CTXCT,CTXCU,CUXCU,EPF
12#endif
13 WM4=WM2**2
14 ZM4=ZM2**2
15 WZM2=WM2*ZM2
16 CTXCT=CT**2/T**2
17 CTXCU=CT*CU/(T*U)
18 CUXCU=CU**2/U**2
19 WZTU=
20 1 +CTXCT*CV3*(-128.*WM2*T1*T3**2-64.*WM2*T1*ZM4-64.*ZM2*T1**2*T3-6
21 1 4.*ZM2*T3*WM4-128.*T1*T3*WZM2-64.*T1**2*T3**2-32.*T1**2*ZM4-64.*
22 1 T3**2*WM4-32.*WM4*ZM4)
23 WZTU=WZTU
24 1 +CTXCT*CV3*T*(64.*WM2*T1*T3-32.*WM2*U1*T3-64.*WM2*T3*S13+64.*WM2
25 1 *T3**2+32.*WM2*ZM4+64.*ZM2*T1*T3-32.*ZM2*T1*S13+32.*ZM2*T1**2+32
26 1 .*ZM2*WM4-32.*T1*U1*T3-64.*T1*T3*S13+64.*T1*T3**2+64.*T1*WZM2+32
27 1 .*T1*ZM4+32.*T1**2*T3+64.*T3*WZM2+32.*T3*WM4-32.*S13*WZM2)
28 1 +CTXCT*CV3*T**2*(-32.*WM2*T3+32.*WM2*S13-32.*ZM2*T1+32.*ZM2*S13-
29 1 32.*T1*T3+32.*T1*S13+32.*U1*T3+64.*T3*S13-32.*WZM2)
30 WZTU=WZTU
31 1 +CTXCT*CV3*T**3*(-32.*S13)
32 1 +CTXCT*CA3*(64.*WM2*T1*ZM4+64.*ZM2*T1**2*T3+64.*ZM2*T3*WM4+128.*
33 1 T1*T3*WZM2+32.*T1**2*ZM4+32.*WM4*ZM4)
34 1 +CTXCT*CA3*T*(-64.*WM2*T1*T3-32.*WM2*U1*T3-32.*WM2*ZM4-64.*ZM2*T
35 1 1*T3+32.*ZM2*T1*S13-32.*ZM2*T1**2-32.*ZM2*WM4-32.*T1*U1*T3-64.*T
36 1 1*WZM2-32.*T1*ZM4-32.*T1**2*T3-64.*T3*WZM2-32.*T3*WM4+32.*S13*WZ
37 1 M2)
38 WZTU=WZTU
39 1 +CTXCT*CA3*T**2*(32.*WM2*T3-32.*WM2*S13+32.*ZM2*T1-32.*ZM2*S13+3
40 1 2.*T1*T3-32.*T1*S13+32.*U1*T3+32.*WZM2)
41 1 +CTXCT*CA3*T**3*(32.*S13)
42 WZTU=WZTU
43 1 +CTXCU*CV3*(32.*WM2*ZM2*WM4+32.*WM2*ZM2*ZM4+64.*WM2*T1*ZM4-128.*
44 1 WM2*U1*T3*U3+64.*WM2*T3*ZM4+64.*WM2*S13*ZM4-64.*ZM2*T1*U1*T3-64.
45 1 *ZM2*T1*U1*U3+64.*ZM2*T1*WM4+64.*ZM2*T3*WM4+64.*ZM2*S13*WM4-128.
46 1 *T1*U1*T3*U3+64.*T1*T3*WZM2+64.*T1*S13*WZM2+32.*T1**2*WZM2-64.*U
47 1 1*U3*WZM2+64.*U1*S13*WZM2+32.*U1**2*WZM2+128.*T3*S13*WZM2+64.*T3
48 1 **2*WZM2+64.*S13**2*WZM2+64.*WM4*ZM4)
49 WZTU=WZTU
50 1 +CTXCU*CV3*T*(32.*WM2*U1*U3-64.*WM2*ZM4+32.*ZM2*U1*T3+32.*ZM2*U1
51 1 *U3-32.*ZM2*U1*S13-32.*ZM2*U1**2-64.*ZM2*WM4+32.*T1*U1*U3-64.*T1
52 1 *WZM2+64.*U1*T3*U3-64.*U1*U3*S13-32.*U1**2*U3-64.*T3*WZM2-64.*S1
53 1 3*WZM2)
54 1 +CTXCU*CV3*T*U*(32.*WM2*T3-32.*WM2*S13+32.*ZM2*T1-32.*ZM2*S13+32
55 1 .*T1*T3-32.*T1*S13-64.*U1*T3-32.*U1*S13-64.*T3*S13-64.*S13**2+32
56 1 .*WZM2)
57 WZTU=WZTU
58 1 +CTXCU*CV3*T**2*(-32.*U1*U3+32.*WZM2)
59 1 +CTXCU*CV3*T**2*U*(32.*S13)
60 1 +CTXCU*CV3*U*(-64.*WM2*T1*T3+32.*WM2*U1*T3-64.*WM2*T3*S13-64.*WM
61 1 2*T3**2-32.*WM2*ZM4-64.*ZM2*T1*T3-32.*ZM2*T1*S13-32.*ZM2*T1**2-3
62 1 2.*ZM2*WM4+32.*T1*U1*T3-64.*T1*T3*S13-64.*T1*T3**2-64.*T1*WZM2-3
63 1 2.*T1*ZM4-32.*T1**2*T3-64.*T3*WZM2-32.*T3*WM4-32.*S13*WZM2)
64 WZTU=WZTU
65 1 +CTXCU*CA3*(-32.*WM2*ZM2*WM4-32.*WM2*ZM2*ZM4-64.*WM2*T1*ZM4-64.*
66 1 WM2*T3*ZM4-64.*WM2*S13*ZM4-64.*ZM2*T1*U1*T3+64.*ZM2*T1*U1*U3-64.
67 1 *ZM2*T1*WM4-64.*ZM2*T3*WM4-64.*ZM2*S13*WM4-64.*T1*T3*WZM2-64.*T1
68 1 *S13*WZM2-32.*T1**2*WZM2+64.*U1*U3*WZM2+64.*U1*S13*WZM2+32.*U1**
69 1 2*WZM2-64.*WM4*ZM4)
70 WZTU=WZTU
71 1 +CTXCU*CA3*T*(-32.*WM2*U1*U3+64.*WM2*ZM4+32.*ZM2*U1*T3-32.*ZM2*U
72 1 1*U3-32.*ZM2*U1*S13-32.*ZM2*U1**2+64.*ZM2*WM4-32.*T1*U1*U3+64.*T
73 1 1*WZM2-32.*U1**2*U3+64.*T3*WZM2+64.*S13*WZM2)
74 1 +CTXCU*CA3*T*U*(-32.*WM2*T3+32.*WM2*S13-32.*ZM2*T1+32.*ZM2*S13-3
75 1 2.*T1*T3+32.*T1*S13-64.*U1*T3-32.*U1*S13-32.*WZM2)
76 1 +CTXCU*CA3*T**2*(32.*U1*U3-32.*WZM2)
77 WZTU=WZTU
78 1 +CTXCU*CA3*T**2*U*(-32.*S13)
79 1 +CTXCU*CA3*U*(64.*WM2*T1*T3+32.*WM2*U1*T3+32.*WM2*ZM4+64.*ZM2*T1
80 1 *T3+32.*ZM2*T1*S13+32.*ZM2*T1**2+32.*ZM2*WM4+32.*T1*U1*T3+64.*T1
81 1 *WZM2+32.*T1*ZM4+32.*T1**2*T3+64.*T3*WZM2+32.*T3*WM4+32.*S13*WZM
82 1 2)
83 1 +CUXCU*CV3*(-64.*ZM2*U1**2*U3-64.*U1**2*U3**2-32.*U1**2*ZM4)
84 WZTU=WZTU
85 1 +CUXCU*CV3*T*U*(32.*U1*U3)
86 1 +CUXCU*CV3*U*(-32.*WM2*U1*U3-32.*ZM2*U1*T3-32.*ZM2*U1*U3-32.*ZM2
87 1 *U1*S13+32.*ZM2*U1**2-32.*T1*U1*U3-64.*U1*T3*U3-64.*U1*U3*S13+32
88 1 .*U1**2*U3)
89 1 +CUXCU*CV3*U**2*(32.*U1*T3+32.*U1*S13)
90 1 +CUXCU*CA3*(-64.*ZM2*U1**2*U3-32.*U1**2*ZM4)
91 WZTU=WZTU
92 1 +CUXCU*CA3*T*U*(-32.*U1*U3)
93 1 +CUXCU*CA3*U*(32.*WM2*U1*U3-32.*ZM2*U1*T3+32.*ZM2*U1*U3-32.*ZM2*
94 1 U1*S13+32.*ZM2*U1**2+32.*T1*U1*U3+32.*U1**2*U3)
95 1 +CUXCU*CA3*U**2*(32.*U1*T3+32.*U1*S13)
96 1 +EPF(P1,P2,P3,Q1)*CTXCU*CV3*(-32.*ZM2*T1-16.*WZM2)
97 1 +EPF(P1,P2,P3,Q1)*CTXCU*CV3*T*(16.*ZM2)
98 WZTU=WZTU
99 1 +EPF(P1,P2,P3,Q3)*CTXCU*CV3*(-48.*WM2*T1-16.*WM2*U1-32.*WM2*S13-
100 1 32.*T1*U1-64.*T1*S13-32.*T1**2-16.*WM4)
101 1 +EPF(P1,P2,P3,Q3)*CTXCU*CV3*T*(16.*WM2+16.*T1+16.*U1+32.*S13)
102 1 +EPF(P1,P2,Q1,Q3)*CTXCU*CV3*(-32.*WM2*T3-32.*WM2*U3-64.*ZM2*T1-6
103 1 4.*T1*T3-64.*T1*U3-32.*WZM2)
104 1 +EPF(P1,P2,Q1,Q3)*CTXCU*CV3*T*(16.*WM2+32.*ZM2+32.*T1+32.*T3+32.
105 1 *U3)
106 WZTU=WZTU
107 1 -16.*EPF(P1,P2,Q1,Q3)*CTXCU*CV3*T*U
108 1 -16.*EPF(P1,P2,Q1,Q3)*CTXCU*CV3*T**2
109 1 +EPF(P1,P2,Q1,Q3)*CTXCU*CV3*U*(16.*WM2+32.*T1)
110 1 +EPF(P1,P3,Q1,Q3)*CTXCU*CV3*(32.*WM2*T3+32.*ZM2*T1+64.*T1*T3+16.
111 1 *WZM2)
112 1 +EPF(P1,P3,Q1,Q3)*CTXCU*CV3*T*(-16.*WM2-16.*ZM2-32.*T1-32.*T3)
113 WZTU=WZTU
114 1 +16.*EPF(P1,P3,Q1,Q3)*CTXCU*CV3*T**2
115 1 +EPF(P2,P3,Q1,Q3)*CTXCU*CV3*(-32.*WM2*U3-32.*ZM2*T1-64.*T1*U3-16
116 1 .*WZM2)
117 1 +EPF(P2,P3,Q1,Q3)*CTXCU*CV3*T*(16.*ZM2+32.*U3)
118 1 -16.*EPF(P2,P3,Q1,Q3)*CTXCU*CV3*T*U
119 1 +EPF(P2,P3,Q1,Q3)*CTXCU*CV3*U*(16.*WM2+32.*T1)+0.
120 RETURN
121 END