[u/mrichter/AliRoot.git] / ISAJET / code / wztu.F

#include "isajet/pilot.h"
      FUNCTION WZTU(T,U,T1,U1,T3,U3,P1,P2)
C          DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980).
C          ALL T AND U GRAPH TERMS.
#include "isajet/itapes.inc"
#include "isajet/wwpar.inc"
      DIMENSION P1(4),P2(4)
#if defined(CERNLIB_DOUBLE)
      DOUBLE PRECISION WZTU
      DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
      DOUBLE PRECISION WM4,ZM4,WZM2,CTXCT,CTXCU,CUXCU,EPF
#endif
      WM4=WM2**2
      ZM4=ZM2**2
      WZM2=WM2*ZM2
      CTXCT=CT**2/T**2
      CTXCU=CT*CU/(T*U)
      CUXCU=CU**2/U**2
      WZTU=
     1 +CTXCT*CV3*(-128.*WM2*T1*T3**2-64.*WM2*T1*ZM4-64.*ZM2*T1**2*T3-6
     1 4.*ZM2*T3*WM4-128.*T1*T3*WZM2-64.*T1**2*T3**2-32.*T1**2*ZM4-64.*
     1 T3**2*WM4-32.*WM4*ZM4)
      WZTU=WZTU
     1 +CTXCT*CV3*T*(64.*WM2*T1*T3-32.*WM2*U1*T3-64.*WM2*T3*S13+64.*WM2
     1 *T3**2+32.*WM2*ZM4+64.*ZM2*T1*T3-32.*ZM2*T1*S13+32.*ZM2*T1**2+32
     1 .*ZM2*WM4-32.*T1*U1*T3-64.*T1*T3*S13+64.*T1*T3**2+64.*T1*WZM2+32
     1 .*T1*ZM4+32.*T1**2*T3+64.*T3*WZM2+32.*T3*WM4-32.*S13*WZM2)
     1 +CTXCT*CV3*T**2*(-32.*WM2*T3+32.*WM2*S13-32.*ZM2*T1+32.*ZM2*S13-
     1 32.*T1*T3+32.*T1*S13+32.*U1*T3+64.*T3*S13-32.*WZM2)
      WZTU=WZTU
     1 +CTXCT*CV3*T**3*(-32.*S13)
     1 +CTXCT*CA3*(64.*WM2*T1*ZM4+64.*ZM2*T1**2*T3+64.*ZM2*T3*WM4+128.*
     1 T1*T3*WZM2+32.*T1**2*ZM4+32.*WM4*ZM4)
     1 +CTXCT*CA3*T*(-64.*WM2*T1*T3-32.*WM2*U1*T3-32.*WM2*ZM4-64.*ZM2*T
     1 1*T3+32.*ZM2*T1*S13-32.*ZM2*T1**2-32.*ZM2*WM4-32.*T1*U1*T3-64.*T
     1 1*WZM2-32.*T1*ZM4-32.*T1**2*T3-64.*T3*WZM2-32.*T3*WM4+32.*S13*WZ
     1 M2)
      WZTU=WZTU
     1 +CTXCT*CA3*T**2*(32.*WM2*T3-32.*WM2*S13+32.*ZM2*T1-32.*ZM2*S13+3
     1 2.*T1*T3-32.*T1*S13+32.*U1*T3+32.*WZM2)
     1 +CTXCT*CA3*T**3*(32.*S13)
      WZTU=WZTU
     1 +CTXCU*CV3*(32.*WM2*ZM2*WM4+32.*WM2*ZM2*ZM4+64.*WM2*T1*ZM4-128.*
     1 WM2*U1*T3*U3+64.*WM2*T3*ZM4+64.*WM2*S13*ZM4-64.*ZM2*T1*U1*T3-64.
     1 *ZM2*T1*U1*U3+64.*ZM2*T1*WM4+64.*ZM2*T3*WM4+64.*ZM2*S13*WM4-128.
     1 *T1*U1*T3*U3+64.*T1*T3*WZM2+64.*T1*S13*WZM2+32.*T1**2*WZM2-64.*U
     1 1*U3*WZM2+64.*U1*S13*WZM2+32.*U1**2*WZM2+128.*T3*S13*WZM2+64.*T3
     1 **2*WZM2+64.*S13**2*WZM2+64.*WM4*ZM4)
      WZTU=WZTU
     1 +CTXCU*CV3*T*(32.*WM2*U1*U3-64.*WM2*ZM4+32.*ZM2*U1*T3+32.*ZM2*U1
     1 *U3-32.*ZM2*U1*S13-32.*ZM2*U1**2-64.*ZM2*WM4+32.*T1*U1*U3-64.*T1
     1 *WZM2+64.*U1*T3*U3-64.*U1*U3*S13-32.*U1**2*U3-64.*T3*WZM2-64.*S1
     1 3*WZM2)
     1 +CTXCU*CV3*T*U*(32.*WM2*T3-32.*WM2*S13+32.*ZM2*T1-32.*ZM2*S13+32
     1 .*T1*T3-32.*T1*S13-64.*U1*T3-32.*U1*S13-64.*T3*S13-64.*S13**2+32
     1 .*WZM2)
      WZTU=WZTU
     1 +CTXCU*CV3*T**2*(-32.*U1*U3+32.*WZM2)
     1 +CTXCU*CV3*T**2*U*(32.*S13)
     1 +CTXCU*CV3*U*(-64.*WM2*T1*T3+32.*WM2*U1*T3-64.*WM2*T3*S13-64.*WM
     1 2*T3**2-32.*WM2*ZM4-64.*ZM2*T1*T3-32.*ZM2*T1*S13-32.*ZM2*T1**2-3
     1 2.*ZM2*WM4+32.*T1*U1*T3-64.*T1*T3*S13-64.*T1*T3**2-64.*T1*WZM2-3
     1 2.*T1*ZM4-32.*T1**2*T3-64.*T3*WZM2-32.*T3*WM4-32.*S13*WZM2)
      WZTU=WZTU
     1 +CTXCU*CA3*(-32.*WM2*ZM2*WM4-32.*WM2*ZM2*ZM4-64.*WM2*T1*ZM4-64.*
     1 WM2*T3*ZM4-64.*WM2*S13*ZM4-64.*ZM2*T1*U1*T3+64.*ZM2*T1*U1*U3-64.
     1 *ZM2*T1*WM4-64.*ZM2*T3*WM4-64.*ZM2*S13*WM4-64.*T1*T3*WZM2-64.*T1
     1 *S13*WZM2-32.*T1**2*WZM2+64.*U1*U3*WZM2+64.*U1*S13*WZM2+32.*U1**
     1 2*WZM2-64.*WM4*ZM4)
      WZTU=WZTU
     1 +CTXCU*CA3*T*(-32.*WM2*U1*U3+64.*WM2*ZM4+32.*ZM2*U1*T3-32.*ZM2*U
     1 1*U3-32.*ZM2*U1*S13-32.*ZM2*U1**2+64.*ZM2*WM4-32.*T1*U1*U3+64.*T
     1 1*WZM2-32.*U1**2*U3+64.*T3*WZM2+64.*S13*WZM2)
     1 +CTXCU*CA3*T*U*(-32.*WM2*T3+32.*WM2*S13-32.*ZM2*T1+32.*ZM2*S13-3
     1 2.*T1*T3+32.*T1*S13-64.*U1*T3-32.*U1*S13-32.*WZM2)
     1 +CTXCU*CA3*T**2*(32.*U1*U3-32.*WZM2)
      WZTU=WZTU
     1 +CTXCU*CA3*T**2*U*(-32.*S13)
     1 +CTXCU*CA3*U*(64.*WM2*T1*T3+32.*WM2*U1*T3+32.*WM2*ZM4+64.*ZM2*T1
     1 *T3+32.*ZM2*T1*S13+32.*ZM2*T1**2+32.*ZM2*WM4+32.*T1*U1*T3+64.*T1
     1 *WZM2+32.*T1*ZM4+32.*T1**2*T3+64.*T3*WZM2+32.*T3*WM4+32.*S13*WZM
     1 2)
     1 +CUXCU*CV3*(-64.*ZM2*U1**2*U3-64.*U1**2*U3**2-32.*U1**2*ZM4)
      WZTU=WZTU
     1 +CUXCU*CV3*T*U*(32.*U1*U3)
     1 +CUXCU*CV3*U*(-32.*WM2*U1*U3-32.*ZM2*U1*T3-32.*ZM2*U1*U3-32.*ZM2
     1 *U1*S13+32.*ZM2*U1**2-32.*T1*U1*U3-64.*U1*T3*U3-64.*U1*U3*S13+32
     1 .*U1**2*U3)
     1 +CUXCU*CV3*U**2*(32.*U1*T3+32.*U1*S13)
     1 +CUXCU*CA3*(-64.*ZM2*U1**2*U3-32.*U1**2*ZM4)
      WZTU=WZTU
     1 +CUXCU*CA3*T*U*(-32.*U1*U3)
     1 +CUXCU*CA3*U*(32.*WM2*U1*U3-32.*ZM2*U1*T3+32.*ZM2*U1*U3-32.*ZM2*
     1 U1*S13+32.*ZM2*U1**2+32.*T1*U1*U3+32.*U1**2*U3)
     1 +CUXCU*CA3*U**2*(32.*U1*T3+32.*U1*S13)
     1 +EPF(P1,P2,P3,Q1)*CTXCU*CV3*(-32.*ZM2*T1-16.*WZM2)
     1 +EPF(P1,P2,P3,Q1)*CTXCU*CV3*T*(16.*ZM2)
      WZTU=WZTU
     1 +EPF(P1,P2,P3,Q3)*CTXCU*CV3*(-48.*WM2*T1-16.*WM2*U1-32.*WM2*S13-
     1 32.*T1*U1-64.*T1*S13-32.*T1**2-16.*WM4)
     1 +EPF(P1,P2,P3,Q3)*CTXCU*CV3*T*(16.*WM2+16.*T1+16.*U1+32.*S13)
     1 +EPF(P1,P2,Q1,Q3)*CTXCU*CV3*(-32.*WM2*T3-32.*WM2*U3-64.*ZM2*T1-6
     1 4.*T1*T3-64.*T1*U3-32.*WZM2)
     1 +EPF(P1,P2,Q1,Q3)*CTXCU*CV3*T*(16.*WM2+32.*ZM2+32.*T1+32.*T3+32.
     1 *U3)
      WZTU=WZTU
     1 -16.*EPF(P1,P2,Q1,Q3)*CTXCU*CV3*T*U
     1 -16.*EPF(P1,P2,Q1,Q3)*CTXCU*CV3*T**2
     1 +EPF(P1,P2,Q1,Q3)*CTXCU*CV3*U*(16.*WM2+32.*T1)
     1 +EPF(P1,P3,Q1,Q3)*CTXCU*CV3*(32.*WM2*T3+32.*ZM2*T1+64.*T1*T3+16.
     1 *WZM2)
     1 +EPF(P1,P3,Q1,Q3)*CTXCU*CV3*T*(-16.*WM2-16.*ZM2-32.*T1-32.*T3)
      WZTU=WZTU
     1 +16.*EPF(P1,P3,Q1,Q3)*CTXCU*CV3*T**2
     1 +EPF(P2,P3,Q1,Q3)*CTXCU*CV3*(-32.*WM2*U3-32.*ZM2*T1-64.*T1*U3-16
     1 .*WZM2)
     1 +EPF(P2,P3,Q1,Q3)*CTXCU*CV3*T*(16.*ZM2+32.*U3)
     1 -16.*EPF(P2,P3,Q1,Q3)*CTXCU*CV3*T*U
     1 +EPF(P2,P3,Q1,Q3)*CTXCU*CV3*U*(16.*WM2+32.*T1)+0.
      RETURN
      END
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=CT2/T2
	17	CTXCU=CTCU/(TU)
	18	CUXCU=CU2/U2
	19	WZTU=
	20	1 +CTXCTCV3(-128.WM2T1T32-64.WM2T1ZM4-64.ZM2T1*2T3-6
	21	1 4.ZM2T3WM4-128.T1T3WZM2-64.T12T3*2-32.T1*2ZM4-64.*
	22	1 T3*2WM4-32.WM4ZM4)
	23	WZTU=WZTU
	24	1 +CTXCTCV3T(64.WM2T1T3-32.WM2U1T3-64.WM2T3S13+64.*WM2
	25	1 T32+32.WM2ZM4+64.ZM2T1T3-32.ZM2T1S13+32.ZM2T1*2+32
	26	1 .ZM2WM4-32.T1U1T3-64.T1T3S13+64.T1T3*2+64.T1*WZM2+32
	27	1 .T1ZM4+32.T12T3+64.T3WZM2+32.T3WM4-32.S13WZM2)
	28	1 +CTXCTCV3T*2(-32.WM2T3+32.WM2S13-32.ZM2T1+32.ZM2S13-
	29	1 32.T1T3+32.T1S13+32.U1T3+64.T3S13-32.*WZM2)
	30	WZTU=WZTU
	31	1 +CTXCTCV3T*3(-32.*S13)
	32	1 +CTXCTCA3(64.WM2T1ZM4+64.ZM2T12T3+64.ZM2T3WM4+128.
	33	1 T1T3WZM2+32.T12ZM4+32.WM4ZM4)
	34	1 +CTXCTCA3T(-64.WM2T1T3-32.WM2U1T3-32.WM2ZM4-64.ZM2*T
	35	1 1T3+32.ZM2T1S13-32.ZM2T1*2-32.ZM2WM4-32.T1U1T3-64.*T
	36	1 1WZM2-32.T1ZM4-32.T1*2T3-64.T3WZM2-32.T3WM4+32.S13WZ
	37	1 M2)
	38	WZTU=WZTU
	39	1 +CTXCTCA3T*2(32.WM2T3-32.WM2S13+32.ZM2T1-32.ZM2S13+3
	40	1 2.T1T3-32.T1S13+32.U1T3+32.*WZM2)
	41	1 +CTXCTCA3T*3(32.*S13)
	42	WZTU=WZTU
	43	1 +CTXCUCV3(32.WM2ZM2WM4+32.WM2ZM2ZM4+64.WM2T1ZM4-128.
	44	1 WM2U1T3U3+64.WM2T3ZM4+64.WM2S13ZM4-64.ZM2T1U1*T3-64.
	45	1 ZM2T1U1U3+64.ZM2T1WM4+64.ZM2T3WM4+64.ZM2S13*WM4-128.
	46	1 T1U1T3U3+64.T1T3WZM2+64.T1S13WZM2+32.T12WZM2-64.*U
	47	1 1U3WZM2+64.U1S13WZM2+32.U1*2WZM2+128.T3S13WZM2+64.T3
	48	1 *2WZM2+64.S132WZM2+64.WM4ZM4)
	49	WZTU=WZTU
	50	1 +CTXCUCV3T(32.WM2U1U3-64.WM2ZM4+32.ZM2U1T3+32.ZM2*U1
	51	1 U3-32.ZM2U1S13-32.ZM2U1*2-64.ZM2WM4+32.T1U1U3-64.*T1
	52	1 WZM2+64.U1T3U3-64.U1U3S13-32.U1*2U3-64.T3WZM2-64.*S1
	53	1 3*WZM2)
	54	1 +CTXCUCV3TU(32.WM2T3-32.WM2S13+32.ZM2T1-32.ZM2S13+32
	55	1 .T1T3-32.T1S13-64.U1T3-32.U1S13-64.T3S13-64.S13*2+32
	56	1 .*WZM2)
	57	WZTU=WZTU
	58	1 +CTXCUCV3T*2(-32.U1U3+32.*WZM2)
	59	1 +CTXCUCV3T*2U(32.S13)
	60	1 +CTXCUCV3U(-64.WM2T1T3+32.WM2U1T3-64.WM2T3S13-64.*WM
	61	1 2T32-32.WM2ZM4-64.ZM2T1T3-32.ZM2T1S13-32.ZM2T1*2-3
	62	1 2.ZM2WM4+32.T1U1T3-64.T1T3S13-64.T1T3*2-64.T1*WZM2-3
	63	1 2.T1ZM4-32.T12T3-64.T3WZM2-32.T3WM4-32.S13WZM2)
	64	WZTU=WZTU
65	1 +CTXCUCA3(-32.WM2ZM2WM4-32.WM2ZM2ZM4-64.WM2T1ZM4-64.
66	1 WM2T3ZM4-64.WM2S13ZM4-64.ZM2T1U1T3+64.ZM2T1U1*U3-64.
67	1 ZM2T1WM4-64.ZM2T3WM4-64.ZM2S13WM4-64.T1T3WZM2-64.*T1
68	1 S13WZM2-32.T12WZM2+64.U1U3WZM2+64.U1S13WZM2+32.U1*
69	1 2WZM2-64.WM4*ZM4)
70	WZTU=WZTU
71	1 +CTXCUCA3T(-32.WM2U1U3+64.WM2ZM4+32.ZM2U1T3-32.ZM2*U
72	1 1U3-32.ZM2U1S13-32.ZM2U1*2+64.ZM2WM4-32.T1U1U3+64.*T
73	1 1WZM2-32.U1*2U3+64.T3WZM2+64.S13WZM2)
74	1 +CTXCUCA3TU(-32.WM2T3+32.WM2S13-32.ZM2T1+32.ZM2S13-3
75	1 2.T1T3+32.T1S13-64.U1T3-32.U1S13-32.*WZM2)
76	1 +CTXCUCA3T*2(32.U1U3-32.*WZM2)
77	WZTU=WZTU
78	1 +CTXCUCA3T*2U(-32.S13)
79	1 +CTXCUCA3U(64.WM2T1T3+32.WM2U1T3+32.WM2ZM4+64.ZM2*T1
80	1 T3+32.ZM2T1S13+32.ZM2T1*2+32.ZM2WM4+32.T1U1T3+64.*T1
81	1 WZM2+32.T1ZM4+32.T1*2T3+64.T3WZM2+32.T3WM4+32.S13WZM
82	1 2)
83	1 +CUXCUCV3(-64.ZM2U1*2U3-64.U12U3*2-32.U1*2ZM4)
84	WZTU=WZTU
85	1 +CUXCUCV3TU(32.U1U3)
86	1 +CUXCUCV3U(-32.WM2U1U3-32.ZM2U1T3-32.ZM2U1U3-32.*ZM2
87	1 U1S13+32.ZM2U1*2-32.T1U1U3-64.U1T3U3-64.U1U3S13+32
88	1 .U12U3)
89	1 +CUXCUCV3U*2(32.U1T3+32.U1S13)
90	1 +CUXCUCA3(-64.ZM2U1*2U3-32.U12ZM4)
91	WZTU=WZTU
92	1 +CUXCUCA3TU(-32.U1U3)
93	1 +CUXCUCA3U(32.WM2U1U3-32.ZM2U1T3+32.ZM2U1U3-32.ZM2
94	1 U1S13+32.ZM2U12+32.T1U1U3+32.U12U3)
95	1 +CUXCUCA3U*2(32.U1T3+32.U1S13)
96	1 +EPF(P1,P2,P3,Q1)CTXCUCV3(-32.ZM2T1-16.WZM2)
97	1 +EPF(P1,P2,P3,Q1)CTXCUCV3T(16.*ZM2)
98	WZTU=WZTU
99	1 +EPF(P1,P2,P3,Q3)CTXCUCV3(-48.WM2T1-16.WM2U1-32.WM2*S13-
100	1 32.T1U1-64.T1S13-32.T12-16.WM4)
101	1 +EPF(P1,P2,P3,Q3)CTXCUCV3T(16.WM2+16.T1+16.U1+32.S13)
102	1 +EPF(P1,P2,Q1,Q3)CTXCUCV3(-32.WM2T3-32.WM2U3-64.ZM2*T1-6
103	1 4.T1T3-64.T1U3-32.*WZM2)
104	1 +EPF(P1,P2,Q1,Q3)CTXCUCV3T(16.WM2+32.ZM2+32.T1+32.T3+32.
105	1 *U3)
106	WZTU=WZTU
107	1 -16.EPF(P1,P2,Q1,Q3)CTXCUCV3T*U
108	1 -16.EPF(P1,P2,Q1,Q3)CTXCUCV3T**2
109	1 +EPF(P1,P2,Q1,Q3)CTXCUCV3U(16.WM2+32.T1)
110	1 +EPF(P1,P3,Q1,Q3)CTXCUCV3(32.WM2T3+32.ZM2T1+64.T1*T3+16.
111	1 *WZM2)
112	1 +EPF(P1,P3,Q1,Q3)CTXCUCV3T(-16.WM2-16.ZM2-32.T1-32.T3)
113	WZTU=WZTU
114	1 +16.EPF(P1,P3,Q1,Q3)CTXCUCV3T**2
115	1 +EPF(P2,P3,Q1,Q3)CTXCUCV3(-32.WM2U3-32.ZM2T1-64.T1*U3-16
116	1 .*WZM2)
117	1 +EPF(P2,P3,Q1,Q3)CTXCUCV3T(16.ZM2+32.U3)
118	1 -16.EPF(P2,P3,Q1,Q3)CTXCUCV3T*U
119	1 +EPF(P2,P3,Q1,Q3)CTXCUCV3U(16.WM2+32.T1)+0.
120	RETURN
121	END