]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | FUNCTION WWST(T,U,T1,U1,T3,U3,P1,P2) | |
3 | C DECAY DISTRIBUTION FOR W+ W- PAIRS FROM SCHOONSCHIP(1980). | |
4 | C INTERFERENCE OF T AND S GRAPHS. | |
5 | #include "isajet/itapes.inc" | |
6 | #include "isajet/wwpar.inc" | |
7 | DIMENSION P1(4),P2(4) | |
8 | #if defined(CERNLIB_DOUBLE) | |
9 | DOUBLE PRECISION WWST | |
10 | DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2 | |
11 | DOUBLE PRECISION CVACQ,EPF | |
12 | #endif | |
13 | CVACQ=(CV+CA)*CQ | |
14 | WWST= | |
15 | 1 +CVACQ*(-256.*T1*U1*T3*U3+256.*T1**2*T3**2) | |
16 | 1 +CVACQ*T*(256.*T1*U1*T3+256.*T1*T3*S13-256.*T1*T3**2+128.*U1*T3* | |
17 | 1 U3-128.*U1*U3*S13-128.*U1**2*U3) | |
18 | 1 +CVACQ*T*U*(-256.*U1*T3-128.*U1*S13-128.*T3*S13-128.*S13**2) | |
19 | 1 +CVACQ*T*WM2*(384.*U1*T3-128.*U1*S13-128.*U1**2+256.*T3*S13-256. | |
20 | 1 *T3**2) | |
21 | 1 +CVACQ*T**2*(-256.*U1*T3-256.*T3*S13) | |
22 | 1 +CVACQ*U*(128.*T1*U1*T3-128.*T1*T3*S13-128.*T1*T3**2) | |
23 | WWST=WWST | |
24 | 1 +CVACQ*U*WM2*(128.*U1*T3-128.*T3*S13-128.*T3**2) | |
25 | 1 +CVACQ*WM2*(-256.*T1*U1*T3+512.*T1*T3**2-256.*U1*T3*U3) | |
26 | 1 +CVACQ*WM2**2*(256.*U1*S13+128.*U1**2+256.*T3*S13+384.*T3**2+128 | |
27 | 1 .*S13**2) | |
28 | 1 +EPF(P1,P2,P3,Q1)*CVACQ*(128.*T3*U3+128.*T3*S13+64.*T3**2+128.*U | |
29 | 1 3*S13+64.*U3**2) | |
30 | 1 +EPF(P1,P2,P3,Q1)*CVACQ*T*(-32.*T3-32.*U3-64.*S13) | |
31 | 1 -32.*EPF(P1,P2,P3,Q1)*CVACQ*T*WM2 | |
32 | WWST=WWST | |
33 | 1 +EPF(P1,P2,P3,Q1)*CVACQ*U*(-32.*T3-32.*U3-64.*S13) | |
34 | 1 -32.*EPF(P1,P2,P3,Q1)*CVACQ*U*WM2 | |
35 | 1 +EPF(P1,P2,P3,Q1)*CVACQ*WM2*(128.*T3+128.*U3+128.*S13) | |
36 | 1 +64.*EPF(P1,P2,P3,Q1)*CVACQ*WM2**2 | |
37 | 1 -32.*EPF(P1,P2,P3,Q3)*CVACQ*T*WM2 | |
38 | 1 -32.*EPF(P1,P2,P3,Q3)*CVACQ*U*WM2 | |
39 | 1 +EPF(P1,P2,P3,Q3)*CVACQ*WM2*(64.*T3+64.*U3) | |
40 | 1 +64.*EPF(P1,P2,P3,Q3)*CVACQ*WM2**2 | |
41 | WWST=WWST | |
42 | 1 +EPF(P1,P3,Q1,Q3)*CVACQ*(128.*U1*T3+128.*U1*U3) | |
43 | 1 +EPF(P1,P3,Q1,Q3)*CVACQ*T*(-64.*U1) | |
44 | 1 +32.*EPF(P1,P3,Q1,Q3)*CVACQ*T*U | |
45 | 1 -32.*EPF(P1,P3,Q1,Q3)*CVACQ*T*WM2 | |
46 | 1 +EPF(P1,P3,Q1,Q3)*CVACQ*U*(-64.*U1-64.*T3-64.*U3) | |
47 | 1 -96.*EPF(P1,P3,Q1,Q3)*CVACQ*U*WM2 | |
48 | 1 +32.*EPF(P1,P3,Q1,Q3)*CVACQ*U**2 | |
49 | 1 +EPF(P1,P3,Q1,Q3)*CVACQ*WM2*(128.*U1+64.*T3+64.*U3) | |
50 | WWST=WWST | |
51 | 1 +64.*EPF(P1,P3,Q1,Q3)*CVACQ*WM2**2 | |
52 | 1 +EPF(P2,P3,Q1,Q3)*CVACQ*(-128.*T1*T3-128.*T1*U3) | |
53 | 1 +EPF(P2,P3,Q1,Q3)*CVACQ*T*(64.*T1+64.*T3+64.*U3) | |
54 | 1 -32.*EPF(P2,P3,Q1,Q3)*CVACQ*T*U | |
55 | 1 +96.*EPF(P2,P3,Q1,Q3)*CVACQ*T*WM2 | |
56 | 1 -32.*EPF(P2,P3,Q1,Q3)*CVACQ*T**2 | |
57 | 1 +EPF(P2,P3,Q1,Q3)*CVACQ*U*(64.*T1) | |
58 | 1 +32.*EPF(P2,P3,Q1,Q3)*CVACQ*U*WM2 | |
59 | 1 +EPF(P2,P3,Q1,Q3)*CVACQ*WM2*(-128.*T1-64.*T3-64.*U3) | |
60 | 1 -64.*EPF(P2,P3,Q1,Q3)*CVACQ*WM2**2 | |
61 | WWST=WWST/T | |
62 | WWST=2.*WWST | |
63 | RETURN | |
64 | END |