]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | FUNCTION WZTU(T,U,T1,U1,T3,U3,P1,P2) | |
3 | C DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980). | |
4 | C 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 |