More volume overlaps corrected
[u/mrichter/AliRoot.git] / ISAJET / code / wwst.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 FUNCTION WWST(T,U,T1,U1,T3,U3,P1,P2)
3C DECAY DISTRIBUTION FOR W+ W- PAIRS FROM SCHOONSCHIP(1980).
4C 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