]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | FUNCTION ZZALL(T,U,T1,U1,T3,U3,P1,P2) | |
3 | C DECAY DISTRIBUTION FOR Z0 Z0 PAIRS FROM SCHOONSCHIP(1980). | |
4 | C INCLUDES TT, TU, AND UU 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 ZZALL | |
10 | DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2 | |
11 | DOUBLE PRECISION CVAVA,DVAVA,CVVVV,CVVAA,TT,TU,UU,EPF | |
12 | #endif | |
13 | TTUU=(T*U)**2 | |
14 | CVAVA=CV*CA*(CV1*CA3+CA1*CV3)/TTUU | |
15 | DVAVA=CV*CA*(CV1*CA3-CA1*CV3)/TTUU | |
16 | CVVVV=(CV**2+CA**2)*CV1*CV3/TTUU | |
17 | CVVAA=(CV**2+CA**2)*CA1*CA3/TTUU | |
18 | TT=T**2 | |
19 | TU=T*U | |
20 | UU=U**2 | |
21 | ZZALL= | |
22 | 1 +CVAVA*T*U*(-16.*ZM2*T1*U1*T3+16.*ZM2*T1*U1*U3+16.*ZM2*T1*T3*U3- | |
23 | 1 16.*ZM2*U1*T3*U3-16.*ZM2**2*T1*S13-8.*ZM2**2*T1**2+16.*ZM2**2*U1 | |
24 | 1 *S13+8.*ZM2**2*U1**2+16.*ZM2**2*T3*S13+8.*ZM2**2*T3**2-16.*ZM2** | |
25 | 1 2*U3*S13-8.*ZM2**2*U3**2) | |
26 | 1 +CVAVA*T*U**2*(16.*ZM2*T1*S13-16.*ZM2*T3*S13) | |
27 | 1 +CVAVA*T**2*(16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-8.*ZM2**2*U1**2+8 | |
28 | 1 .*ZM2**2*U3**2) | |
29 | 1 +CVAVA*T**2*U*(-16.*ZM2*U1*S13+16.*ZM2*U3*S13) | |
30 | 1 +CVAVA*U**2*(-16.*ZM2*T1*T3**2+16.*ZM2*T1**2*T3+8.*ZM2**2*T1**2- | |
31 | 1 8.*ZM2**2*T3**2) | |
32 | ZZALL=ZZALL | |
33 | 1 +DVAVA*T*U*(-16.*ZM2*T1*U1*T3+16.*ZM2*T1*U1*U3-16.*ZM2*T1*T3*U3+ | |
34 | 1 16.*ZM2*U1*T3*U3-32.*ZM2**2*T1*T3-16.*ZM2**2*T1*S13-8.*ZM2**2*T1 | |
35 | 1 **2+32.*ZM2**2*U1*U3+16.*ZM2**2*U1*S13+8.*ZM2**2*U1**2-16.*ZM2** | |
36 | 1 2*T3*S13-8.*ZM2**2*T3**2+16.*ZM2**2*U3*S13+8.*ZM2**2*U3**2-32.*Z | |
37 | 1 M2**3*T1+32.*ZM2**3*U1-32.*ZM2**3*T3+32.*ZM2**3*U3) | |
38 | 1 +DVAVA*T*U**2*(16.*ZM2*T1*S13+16.*ZM2*T3*S13-16.*ZM2**2*U1-16.*Z | |
39 | 1 M2**2*U3-32.*ZM2**3) | |
40 | 1 +DVAVA*T*U**3*(8.*ZM2**2) | |
41 | ZZALL=ZZALL | |
42 | 1 +DVAVA*T**2*(-16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-32.*ZM2**2*U1*U3 | |
43 | 1 -8.*ZM2**2*U1**2-8.*ZM2**2*U3**2-16.*ZM2**3*U1-16.*ZM2**3*U3-8.* | |
44 | 1 ZM2**4) | |
45 | 1 +DVAVA*T**2*U*(-16.*ZM2*U1*S13-16.*ZM2*U3*S13+16.*ZM2**2*T1+16.* | |
46 | 1 ZM2**2*T3+32.*ZM2**3) | |
47 | 1 +DVAVA*T**3*U*(-8.*ZM2**2) | |
48 | 1 +DVAVA*U**2*(16.*ZM2*T1*T3**2+16.*ZM2*T1**2*T3+32.*ZM2**2*T1*T3+ | |
49 | 1 8.*ZM2**2*T1**2+8.*ZM2**2*T3**2+16.*ZM2**3*T1+16.*ZM2**3*T3+8.*Z | |
50 | 1 M2**4) | |
51 | ZZALL=ZZALL | |
52 | 1 +CVVVV*T*U*(-16.*ZM2*T1*U1*T3-16.*ZM2*T1*U1*U3-16.*ZM2*T1*T3*U3- | |
53 | 1 16.*ZM2*U1*T3*U3+16.*ZM2**2*T1*S13+8.*ZM2**2*T1**2+16.*ZM2**2*U1 | |
54 | 1 *S13+8.*ZM2**2*U1**2+16.*ZM2**2*T3*S13+8.*ZM2**2*T3**2+16.*ZM2** | |
55 | 1 2*U3*S13+8.*ZM2**2*U3**2+16.*ZM2**2*S13**2+16.*ZM2**3*T1+16.*ZM2 | |
56 | 1 **3*U1+16.*ZM2**3*T3+16.*ZM2**3*U3+32.*ZM2**3*S13+32.*ZM2**4-32. | |
57 | 1 *T1*U1*T3*U3) | |
58 | 1 +CVVVV*T*U**2*(-16.*ZM2*T1*S13-16.*ZM2*T3*S13-8.*ZM2**2*U1-8.*ZM | |
59 | 1 2**2*U3-16.*ZM2**2*S13-16.*ZM2**3-32.*T1*T3*S13) | |
60 | ZZALL=ZZALL | |
61 | 1 +CVVVV*T*U**3*(4.*ZM2**2) | |
62 | 1 +CVVVV*T**2*(-16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-16.*ZM2**2*U1*U3 | |
63 | 1 -8.*ZM2**2*U1**2-8.*ZM2**2*U3**2-8.*ZM2**3*U1-8.*ZM2**3*U3-4.*ZM | |
64 | 1 2**4-16.*U1**2*U3**2) | |
65 | 1 +CVVVV*T**2*U*(-16.*ZM2*U1*S13-16.*ZM2*U3*S13-8.*ZM2**2*T1-8.*ZM | |
66 | 1 2**2*T3-16.*ZM2**2*S13-16.*ZM2**3-32.*U1*U3*S13) | |
67 | 1 +CVVVV*T**2*U**2*(-16.*S13**2) | |
68 | 1 +CVVVV*T**3*U*(4.*ZM2**2) | |
69 | ZZALL=ZZALL | |
70 | 1 +CVVVV*U**2*(-16.*ZM2*T1*T3**2-16.*ZM2*T1**2*T3-16.*ZM2**2*T1*T3 | |
71 | 1 -8.*ZM2**2*T1**2-8.*ZM2**2*T3**2-8.*ZM2**3*T1-8.*ZM2**3*T3-4.*ZM | |
72 | 1 2**4-16.*T1**2*T3**2) | |
73 | 1 +CVVAA*T*U*(-16.*ZM2**3*T1-16.*ZM2**3*U1-16.*ZM2**3*T3-16.*ZM2** | |
74 | 1 3*U3-32.*ZM2**3*S13-32.*ZM2**4) | |
75 | 1 +CVVAA*T*U**2*(8.*ZM2**2*U1+8.*ZM2**2*U3+16.*ZM2**2*S13+16.*ZM2* | |
76 | 1 *3) | |
77 | 1 +CVVAA*T*U**3*(-4.*ZM2**2) | |
78 | ZZALL=ZZALL | |
79 | 1 +CVVAA*T**2*(16.*ZM2**2*U1*U3+8.*ZM2**3*U1+8.*ZM2**3*U3+4.*ZM2** | |
80 | 1 4) | |
81 | 1 +CVVAA*T**2*U*(8.*ZM2**2*T1+8.*ZM2**2*T3+16.*ZM2**2*S13+16.*ZM2* | |
82 | 1 *3) | |
83 | 1 +CVVAA*T**3*U*(-4.*ZM2**2) | |
84 | 1 +CVVAA*U**2*(16.*ZM2**2*T1*T3+8.*ZM2**3*T1+8.*ZM2**3*T3+4.*ZM2** | |
85 | 1 4) | |
86 | C THE EPF TERMS FROM SCHOONSCHIP DID NOT CONTAIN CVAVA, ETC., | |
87 | C BUT DID CONTAIN AN EXPLICIT T*U. THIS WAS REPLACED BY HAND | |
88 | C BY 1./TU. | |
89 | ZZALL=ZZALL | |
90 | 1 +EPF(P1,P2,P3,Q3)/TU*(-32.*CV*CA*CV1*CV3*T1*S13) | |
91 | 1 +EPF(P1,P2,Q1,Q3)/TU*(-32.*ZM2*CV*CA*CV1*CV3*T1-32.*CV*CA*CV1*C | |
92 | 1 V3*T1*T3-32.*CV*CA*CV1*CV3*T1*U3) | |
93 | 1 +EPF(P1,P3,Q1,Q3)/TU*(32.*CV*CA*CV1*CV3*T1*T3) | |
94 | 1 +EPF(P2,P3,Q1,Q3)/TU*(-32.*CV*CA*CV1*CV3*T1*U3)+0. | |
95 | ZZALL=2.*ZZALL | |
96 | RETURN | |
97 | END |