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

#include "isajet/pilot.h"
      FUNCTION WZSU(T,U,T1,U1,T3,U3,P1,P2)
C          DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980).
C          INTERFERENCE OF S AND U GRAPHS.
#include "isajet/itapes.inc"
#include "isajet/wwpar.inc"
      DIMENSION P1(4),P2(4)
#if defined(CERNLIB_DOUBLE)
      DOUBLE PRECISION WZSU
      DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
      DOUBLE PRECISION WM4,ZM4,WZM2,CSXCU,EPF
#endif
      WM4=WM2**2
      ZM4=ZM2**2
      WZM2=WM2*ZM2
      CSXCU=CS*CU/U
      WZSU=
     1 +CSXCU*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-128.*ZM2*U1**2*U3-64.*ZM2*T3*WM4-6
     1 4.*ZM2*S13*WM4+128.*T1*U1*T3*U3-64.*T1*T3*WZM2-64.*T1*S13*WZM2-3
     1 2.*T1**2*WZM2+64.*U1*U3*WZM2-64.*U1*S13*WZM2-128.*U1**2*U3**2-32
     1 .*U1**2*WZM2-64.*U1**2*ZM4-128.*T3*S13*WZM2-64.*T3**2*WZM2-64.*S
     1 13**2*WZM2-64.*WM4*ZM4)
      WZSU=WZSU
     1 +CSXCU*CV3*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-64.*U1*T3*U3+64.*U1*U3*S13+32.*U1**2*U3+64.*T3*WZM2+64.*S
     1 13*WZM2)
     1 +CSXCU*CV3*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+64.*U1*U3+32.*U1*S13+64.*T3*S13+64
     1 .*S13**2-32.*WZM2)
      WZSU=WZSU
     1 +CSXCU*CV3*T**2*(32.*U1*U3-32.*WZM2)
     1 +CSXCU*CV3*T**2*U*(-32.*S13)
      WZSU=WZSU
     1 +CSXCU*CV3*U*(64.*WM2*T1*T3-32.*WM2*U1*T3-64.*WM2*U1*U3+64.*WM2*
     1 T3*S13+64.*WM2*T3**2+32.*WM2*ZM4+64.*ZM2*T1*T3+32.*ZM2*T1*S13+32
     1 .*ZM2*T1**2-64.*ZM2*U1*T3-64.*ZM2*U1*U3-64.*ZM2*U1*S13+64.*ZM2*U
     1 1**2+32.*ZM2*WM4-32.*T1*U1*T3-64.*T1*U1*U3+64.*T1*T3*S13+64.*T1*
     1 T3**2+64.*T1*WZM2+32.*T1*ZM4+32.*T1**2*T3-128.*U1*T3*U3-128.*U1*
     1 U3*S13+64.*U1**2*U3+64.*T3*WZM2+32.*T3*WM4+32.*S13*WZM2)
      WZSU=WZSU
     1 +CSXCU*CV3*U**2*(64.*U1*T3+64.*U1*S13)
     1 +CSXCU*CA3*(32.*WM2*ZM2*WM4+32.*WM2*ZM2*ZM4+64.*WM2*T1*ZM4+64.*W
     1 M2*T3*ZM4+64.*WM2*S13*ZM4+64.*ZM2*T1*U1*T3-64.*ZM2*T1*U1*U3+64.*
     1 ZM2*T1*WM4-128.*ZM2*U1**2*U3+64.*ZM2*T3*WM4+64.*ZM2*S13*WM4+64.*
     1 T1*T3*WZM2+64.*T1*S13*WZM2+32.*T1**2*WZM2-64.*U1*U3*WZM2-64.*U1*
     1 S13*WZM2-32.*U1**2*WZM2-64.*U1**2*ZM4+64.*WM4*ZM4)
      WZSU=WZSU
     1 +CSXCU*CA3*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+32.*U1**2*U3-64.*T3*WZM2-64.*S13*WZM2)
     1 +CSXCU*CA3*T*U*(32.*WM2*T3-32.*WM2*S13+32.*ZM2*T1-32.*ZM2*S13+32
     1 .*T1*T3-32.*T1*S13+64.*U1*T3-64.*U1*U3+32.*U1*S13+32.*WZM2)
     1 +CSXCU*CA3*T**2*(-32.*U1*U3+32.*WZM2)
      WZSU=WZSU
     1 +CSXCU*CA3*T**2*U*(32.*S13)
     1 +CSXCU*CA3*U*(-64.*WM2*T1*T3-32.*WM2*U1*T3+64.*WM2*U1*U3-32.*WM2
     1 *ZM4-64.*ZM2*T1*T3-32.*ZM2*T1*S13-32.*ZM2*T1**2-64.*ZM2*U1*T3+64
     1 .*ZM2*U1*U3-64.*ZM2*U1*S13+64.*ZM2*U1**2-32.*ZM2*WM4-32.*T1*U1*T
     1 3+64.*T1*U1*U3-64.*T1*WZM2-32.*T1*ZM4-32.*T1**2*T3+64.*U1**2*U3-
     1 64.*T3*WZM2-32.*T3*WM4-32.*S13*WZM2)
      WZSU=WZSU
     1 +CSXCU*CA3*U**2*(64.*U1*T3+64.*U1*S13)
     1 +EPF(P1,P2,P3,Q1)*CSXCU*CV3*(32.*ZM2*T1+32.*ZM2*U1+32.*WZM2)
     1 +EPF(P1,P2,P3,Q1)*CSXCU*CV3*T*(-16.*ZM2)
     1 +EPF(P1,P2,P3,Q1)*CSXCU*CV3*U*(-16.*ZM2)
     1 +EPF(P1,P2,P3,Q3)*CSXCU*CV3*(64.*WM2*T1+64.*WM2*U1+64.*WM2*S13+6
     1 4.*T1*U1+64.*T1*S13+32.*T1**2+64.*U1*S13+32.*U1**2+32.*WM4)
      WZSU=WZSU
     1 +EPF(P1,P2,P3,Q3)*CSXCU*CV3*T*(-16.*WM2-16.*T1-16.*U1-32.*S13)
     1 +EPF(P1,P2,P3,Q3)*CSXCU*CV3*U*(-16.*WM2-16.*T1-16.*U1-32.*S13)
     1 +EPF(P1,P2,Q1,Q3)*CSXCU*CV3*(64.*WM2*T3+64.*WM2*U3+64.*ZM2*T1+64
     1 .*ZM2*U1+64.*T1*T3+64.*T1*U3+64.*U1*T3+64.*U1*U3+64.*WZM2)
     1 +EPF(P1,P2,Q1,Q3)*CSXCU*CV3*T*(-32.*WM2-32.*ZM2-32.*T1-32.*U1-32
     1 .*T3-32.*U3)
      WZSU=WZSU
     1 +32.*EPF(P1,P2,Q1,Q3)*CSXCU*CV3*T*U
     1 +16.*EPF(P1,P2,Q1,Q3)*CSXCU*CV3*T**2
     1 +EPF(P1,P2,Q1,Q3)*CSXCU*CV3*U*(-32.*WM2-32.*ZM2-32.*T1-32.*U1-32
     1 .*T3-32.*U3)
     1 +16.*EPF(P1,P2,Q1,Q3)*CSXCU*CV3*U**2
     1 +EPF(P1,P3,Q1,Q3)*CSXCU*CV3*(-64.*WM2*T3-32.*ZM2*T1-32.*ZM2*U1-6
     1 4.*T1*T3-64.*U1*T3-32.*WZM2)
      WZSU=WZSU
     1 +EPF(P1,P3,Q1,Q3)*CSXCU*CV3*T*(32.*WM2+16.*ZM2+32.*T1+32.*U1+32.
     1 *T3)
     1 -16.*EPF(P1,P3,Q1,Q3)*CSXCU*CV3*T*U
     1 -16.*EPF(P1,P3,Q1,Q3)*CSXCU*CV3*T**2
     1 +EPF(P1,P3,Q1,Q3)*CSXCU*CV3*U*(16.*ZM2+32.*T3)
     1 +EPF(P2,P3,Q1,Q3)*CSXCU*CV3*(64.*WM2*U3+32.*ZM2*T1+32.*ZM2*U1+64
     1 .*T1*U3+64.*U1*U3+32.*WZM2)
      WZSU=WZSU
     1 +EPF(P2,P3,Q1,Q3)*CSXCU*CV3*T*(-16.*ZM2-32.*U3)
     1 +16.*EPF(P2,P3,Q1,Q3)*CSXCU*CV3*T*U
     1 +EPF(P2,P3,Q1,Q3)*CSXCU*CV3*U*(-32.*WM2-16.*ZM2-32.*T1-32.*U1-32
     1 .*U3)
     1 +16.*EPF(P2,P3,Q1,Q3)*CSXCU*CV3*U**2
      RETURN
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	FUNCTION WZSU(T,U,T1,U1,T3,U3,P1,P2)
	3	C DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980).
	4	C INTERFERENCE OF S AND U 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 WZSU
	10	DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
	11	DOUBLE PRECISION WM4,ZM4,WZM2,CSXCU,EPF
	12	#endif
	13	WM4=WM2**2
	14	ZM4=ZM2**2
	15	WZM2=WM2*ZM2
	16	CSXCU=CS*CU/U
	17	WZSU=
	18	1 +CSXCUCV3(-32.WM2ZM2WM4-32.WM2ZM2ZM4-64.WM2T1*ZM4+128.
	19	1 WM2U1T3U3-64.WM2T3ZM4-64.WM2S13ZM4+64.ZM2T1U1T3+64
	20	1 .ZM2T1U1U3-64.ZM2T1WM4-128.ZM2U12U3-64.ZM2T3*WM4-6
	21	1 4.ZM2S13WM4+128.T1U1T3U3-64.T1T3WZM2-64.T1S13*WZM2-3
	22	1 2.T12WZM2+64.U1U3WZM2-64.U1S13WZM2-128.U12U3**2-32
	23	1 .U12WZM2-64.U12ZM4-128.T3S13WZM2-64.T3*2WZM2-64.*S
	24	1 13*2WZM2-64.WM4ZM4)
	25	WZSU=WZSU
	26	1 +CSXCUCV3T(-32.WM2U1U3+64.WM2ZM4-32.ZM2U1T3-32.ZM2*U
	27	1 1U3+32.ZM2U1S13+32.ZM2U1*2+64.ZM2WM4-32.T1U1U3+64.*T
	28	1 1WZM2-64.U1T3U3+64.U1U3S13+32.U1*2U3+64.T3WZM2+64.*S
	29	1 13*WZM2)
	30	1 +CSXCUCV3TU(-32.WM2T3+32.WM2S13-32.ZM2T1+32.ZM2S13-3
	31	1 2.T1T3+32.T1S13+64.U1T3+64.U1U3+32.U1S13+64.T3S13+64
	32	1 .S132-32.WZM2)
	33	WZSU=WZSU
	34	1 +CSXCUCV3T*2(32.U1U3-32.*WZM2)
	35	1 +CSXCUCV3T*2U(-32.S13)
	36	WZSU=WZSU
	37	1 +CSXCUCV3U(64.WM2T1T3-32.WM2U1T3-64.WM2U1U3+64.WM2
	38	1 T3S13+64.WM2T32+32.WM2ZM4+64.ZM2T1T3+32.ZM2T1*S13+32
	39	1 .ZM2T1*2-64.ZM2U1T3-64.ZM2U1U3-64.ZM2U1S13+64.ZM2U
	40	1 1*2+32.ZM2WM4-32.T1U1T3-64.T1U1U3+64.T1T3S13+64.T1
	41	1 T3*2+64.T1WZM2+32.T1ZM4+32.T1*2T3-128.U1T3U3-128.U1*
	42	1 U3S13+64.U1*2U3+64.T3WZM2+32.T3WM4+32.S13WZM2)
	43	WZSU=WZSU
	44	1 +CSXCUCV3U*2(64.U1T3+64.U1S13)
	45	1 +CSXCUCA3(32.WM2ZM2WM4+32.WM2ZM2ZM4+64.WM2T1ZM4+64.W
	46	1 M2T3ZM4+64.WM2S13ZM4+64.ZM2T1U1T3-64.ZM2T1U1U3+64.
	47	1 ZM2T1WM4-128.ZM2U1*2U3+64.ZM2T3WM4+64.ZM2S13WM4+64.*
	48	1 T1T3WZM2+64.T1S13WZM2+32.T1*2WZM2-64.U1U3WZM2-64.U1*
	49	1 S13WZM2-32.U1*2WZM2-64.U12ZM4+64.WM4ZM4)
	50	WZSU=WZSU
	51	1 +CSXCUCA3T(32.WM2U1U3-64.WM2ZM4-32.ZM2U1T3+32.ZM2*U1
	52	1 U3+32.ZM2U1S13+32.ZM2U1*2-64.ZM2WM4+32.T1U1U3-64.*T1
	53	1 WZM2+32.U1*2U3-64.T3WZM2-64.S13WZM2)
	54	1 +CSXCUCA3TU(32.WM2T3-32.WM2S13+32.ZM2T1-32.ZM2S13+32
	55	1 .T1T3-32.T1S13+64.U1T3-64.U1U3+32.U1S13+32.*WZM2)
	56	1 +CSXCUCA3T*2(-32.U1U3+32.*WZM2)
	57	WZSU=WZSU
	58	1 +CSXCUCA3T*2U(32.S13)
	59	1 +CSXCUCA3U(-64.WM2T1T3-32.WM2U1T3+64.WM2U1U3-32.*WM2
	60	1 ZM4-64.ZM2T1T3-32.ZM2T1S13-32.ZM2T12-64.ZM2U1T3+64
	61	1 .ZM2U1U3-64.ZM2U1S13+64.ZM2U1*2-32.ZM2WM4-32.T1U1T
	62	1 3+64.T1U1U3-64.T1WZM2-32.T1ZM4-32.T1*2T3+64.U12U3-
	63	1 64.T3WZM2-32.T3WM4-32.S13WZM2)
	64	WZSU=WZSU
65	1 +CSXCUCA3U*2(64.U1T3+64.U1S13)
66	1 +EPF(P1,P2,P3,Q1)CSXCUCV3(32.ZM2T1+32.ZM2U1+32.WZM2)
67	1 +EPF(P1,P2,P3,Q1)CSXCUCV3T(-16.*ZM2)
68	1 +EPF(P1,P2,P3,Q1)CSXCUCV3U(-16.*ZM2)
69	1 +EPF(P1,P2,P3,Q3)CSXCUCV3(64.WM2T1+64.WM2U1+64.WM2*S13+6
70	1 4.T1U1+64.T1S13+32.T12+64.U1S13+32.U1*2+32.WM4)
71	WZSU=WZSU
72	1 +EPF(P1,P2,P3,Q3)CSXCUCV3T(-16.WM2-16.T1-16.U1-32.S13)
73	1 +EPF(P1,P2,P3,Q3)CSXCUCV3U(-16.WM2-16.T1-16.U1-32.S13)
74	1 +EPF(P1,P2,Q1,Q3)CSXCUCV3(64.WM2T3+64.WM2U3+64.ZM2*T1+64
75	1 .ZM2U1+64.T1T3+64.T1U3+64.U1T3+64.U1U3+64.*WZM2)
76	1 +EPF(P1,P2,Q1,Q3)CSXCUCV3T(-32.WM2-32.ZM2-32.T1-32.U1-32
77	1 .T3-32.U3)
78	WZSU=WZSU
79	1 +32.EPF(P1,P2,Q1,Q3)CSXCUCV3T*U
80	1 +16.EPF(P1,P2,Q1,Q3)CSXCUCV3T**2
81	1 +EPF(P1,P2,Q1,Q3)CSXCUCV3U(-32.WM2-32.ZM2-32.T1-32.U1-32
82	1 .T3-32.U3)
83	1 +16.EPF(P1,P2,Q1,Q3)CSXCUCV3U**2
84	1 +EPF(P1,P3,Q1,Q3)CSXCUCV3(-64.WM2T3-32.ZM2T1-32.ZM2*U1-6
85	1 4.T1T3-64.U1T3-32.*WZM2)
86	WZSU=WZSU
87	1 +EPF(P1,P3,Q1,Q3)CSXCUCV3T(32.WM2+16.ZM2+32.T1+32.U1+32.
88	1 *T3)
89	1 -16.EPF(P1,P3,Q1,Q3)CSXCUCV3T*U
90	1 -16.EPF(P1,P3,Q1,Q3)CSXCUCV3T**2
91	1 +EPF(P1,P3,Q1,Q3)CSXCUCV3U(16.ZM2+32.T3)
92	1 +EPF(P2,P3,Q1,Q3)CSXCUCV3(64.WM2U3+32.ZM2T1+32.ZM2*U1+64
93	1 .T1U3+64.U1U3+32.*WZM2)
94	WZSU=WZSU
95	1 +EPF(P2,P3,Q1,Q3)CSXCUCV3T(-16.ZM2-32.U3)
96	1 +16.EPF(P2,P3,Q1,Q3)CSXCUCV3T*U
97	1 +EPF(P2,P3,Q1,Q3)CSXCUCV3U(-32.WM2-16.ZM2-32.T1-32.U1-32
98	1 .*U3)
99	1 +16.EPF(P2,P3,Q1,Q3)CSXCUCV3U**2
100	RETURN
101	END