]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | FUNCTION WWSS(T,U,T1,U1,T3,U3) | |
3 | C DECAY DISTRIBUTION FOR W+ W- PAIRS FROM SCHOONSCHIP(1980). | |
4 | C SQUARE OF S GRAPH. | |
5 | #include "isajet/itapes.inc" | |
6 | #include "isajet/wwpar.inc" | |
7 | #if defined(CERNLIB_DOUBLE) | |
8 | DOUBLE PRECISION WWSS | |
9 | DOUBLE PRECISION T,U,T1,U1,T3,U3 | |
10 | DOUBLE PRECISION CV2A2 | |
11 | #endif | |
12 | CV2A2=CV**2+CA**2 | |
13 | WWSS= | |
14 | 1 +CV*CA*T*(-64.*T1*U1*T3+64.*T1*U1*U3+64.*T1*T3*U3+64.*T1*T3**2-6 | |
15 | 1 4.*T1**2*T3-64.*U1*T3*U3-64.*U1*U3**2+64.*U1**2*U3) | |
16 | 1 +CV*CA*T*U*(-128.*T1*U3-64.*T1*S13+128.*U1*T3+64.*U1*S13+64.*T3* | |
17 | 1 S13-64.*U3*S13) | |
18 | 1 +CV*CA*T*WM2*(128.*T1*U3+64.*T1*S13-64.*T1**2-128.*U1*T3+64.*U1* | |
19 | 1 S13+64.*U1**2-64.*T3*S13+64.*T3**2-64.*U3*S13-64.*U3**2) | |
20 | 1 +CV*CA*T**2*(-64.*T1*U3-64.*T1*S13+64.*U1*T3+64.*T3*S13) | |
21 | 1 +CV*CA*U*(-64.*T1*U1*T3+64.*T1*U1*U3+64.*T1*T3*U3+64.*T1*T3**2-6 | |
22 | 1 4.*T1**2*T3-64.*U1*T3*U3-64.*U1*U3**2+64.*U1**2*U3) | |
23 | WWSS=WWSS | |
24 | 1 +CV*CA*U*WM2*(128.*T1*U3-64.*T1*S13-64.*T1**2-128.*U1*T3-64.*U1* | |
25 | 1 S13+64.*U1**2+64.*T3*S13+64.*T3**2+64.*U3*S13-64.*U3**2) | |
26 | 1 +CV*CA*U**2*(-64.*T1*U3+64.*U1*T3+64.*U1*S13-64.*U3*S13) | |
27 | 1 +CV*CA*WM2*(128.*T1*U1*T3-128.*T1*U1*U3-128.*T1*T3*U3-128.*T1*T3 | |
28 | 1 **2+128.*T1**2*T3+128.*U1*T3*U3+128.*U1*U3**2-128.*U1**2*U3) | |
29 | 1 +CV*CA*WM2**2*(128.*T1*S13+128.*T1**2-128.*U1*S13-128.*U1**2-128 | |
30 | 1 .*T3*S13-128.*T3**2+128.*U3*S13+128.*U3**2) | |
31 | 1 +CV2A2*(128.*T1*U1*T3*U3-64.*T1**2*T3**2-64.*U1**2*U3**2) | |
32 | WWSS=WWSS | |
33 | 1 +CV2A2*T*(-32.*T1*U1*T3-32.*T1*U1*U3-32.*T1*T3*U3-64.*T1*T3*S13+ | |
34 | 1 32.*T1*T3**2+32.*T1**2*T3-32.*U1*T3*U3+64.*U1*U3*S13+32.*U1*U3** | |
35 | 1 2+32.*U1**2*U3) | |
36 | 1 +CV2A2*T*U*(64.*T1*U3+32.*T1*S13+64.*U1*T3+32.*U1*S13+32.*T3*S13 | |
37 | 1 +32.*U3*S13+64.*S13**2) | |
38 | 1 +CV2A2*T*WM2*(-64.*T1*U3-32.*T1*S13+32.*T1**2-64.*U1*T3+32.*U1*S | |
39 | 1 13+32.*U1**2-32.*T3*S13+32.*T3**2+32.*U3*S13+32.*U3**2) | |
40 | 1 +CV2A2*T**2*(32.*T1*U3+32.*T1*S13+32.*U1*T3+32.*T3*S13) | |
41 | WWSS=WWSS | |
42 | 1 +CV2A2*U*(-32.*T1*U1*T3-32.*T1*U1*U3-32.*T1*T3*U3+64.*T1*T3*S13+ | |
43 | 1 32.*T1*T3**2+32.*T1**2*T3-32.*U1*T3*U3-64.*U1*U3*S13+32.*U1*U3** | |
44 | 1 2+32.*U1**2*U3) | |
45 | 1 +CV2A2*U*WM2*(-64.*T1*U3+32.*T1*S13+32.*T1**2-64.*U1*T3-32.*U1*S | |
46 | 1 13+32.*U1**2+32.*T3*S13+32.*T3**2-32.*U3*S13+32.*U3**2) | |
47 | 1 +CV2A2*U**2*(32.*T1*U3+32.*U1*T3+32.*U1*S13+32.*U3*S13) | |
48 | 1 +CV2A2*WM2*(64.*T1*U1*T3+64.*T1*U1*U3+64.*T1*T3*U3-64.*T1*T3**2- | |
49 | 1 64.*T1**2*T3+64.*U1*T3*U3-64.*U1*U3**2-64.*U1**2*U3) | |
50 | WWSS=WWSS | |
51 | 1 +CV2A2*WM2**2*(-64.*T1*S13-64.*T1**2-64.*U1*S13-64.*U1**2-64.*T3 | |
52 | 1 *S13-64.*T3**2-64.*U3*S13-64.*U3**2-64.*S13**2) | |
53 | WWSS=2.*WWSS | |
54 | RETURN | |
55 | END |