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

#include "isajet/pilot.h"
      FUNCTION WZST(T,U,T1,U1,T3,U3,P1,P2)
C          DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980).
C          INTERFERENCE OF S AND T GRAPHS.
#include "isajet/itapes.inc"
#include "isajet/wwpar.inc"
      DIMENSION P1(4),P2(4)
#if defined(CERNLIB_DOUBLE)
      DOUBLE PRECISION WZST
      DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
      DOUBLE PRECISION WM4,ZM4,WZM2,CSXCT,EPF
#endif
      WM4=WM2**2
      ZM4=ZM2**2
      WZM2=WM2*ZM2
      CSXCT=CS*CT/T
      WZST=
     1 +CSXCT*CV3*(32.*WM2*ZM2*WM4+32.*WM2*ZM2*ZM4+256.*WM2*T1*T3**2+19
     1 2.*WM2*T1*ZM4-128.*WM2*U1*T3*U3+64.*WM2*T3*ZM4+64.*WM2*S13*ZM4-6
     1 4.*ZM2*T1*U1*T3-64.*ZM2*T1*U1*U3+64.*ZM2*T1*WM4+128.*ZM2*T1**2*T
     1 3+192.*ZM2*T3*WM4+64.*ZM2*S13*WM4-128.*T1*U1*T3*U3+320.*T1*T3*WZ
     1 M2+64.*T1*S13*WZM2+128.*T1**2*T3**2+32.*T1**2*WZM2+64.*T1**2*ZM4
     1 -64.*U1*U3*WZM2+64.*U1*S13*WZM2+32.*U1**2*WZM2+128.*T3*S13*WZM2+
     1 64.*T3**2*WZM2+128.*T3**2*WM4+64.*S13**2*WZM2+128.*WM4*ZM4)
      WZST=WZST
     1 +CSXCT*CV3*T*(-128.*WM2*T1*T3+64.*WM2*U1*T3+32.*WM2*U1*U3+128.*W
     1 M2*T3*S13-128.*WM2*T3**2-128.*WM2*ZM4-128.*ZM2*T1*T3+64.*ZM2*T1*
     1 S13-64.*ZM2*T1**2+32.*ZM2*U1*T3+32.*ZM2*U1*U3-32.*ZM2*U1*S13-32.
     1 *ZM2*U1**2-128.*ZM2*WM4+64.*T1*U1*T3+32.*T1*U1*U3+128.*T1*T3*S13
     1 -128.*T1*T3**2-192.*T1*WZM2-64.*T1*ZM4-64.*T1**2*T3+64.*U1*T3*U3
     1 -64.*U1*U3*S13-32.*U1**2*U3-192.*T3*WZM2-64.*T3*WM4)
      WZST=WZST
     1 +CSXCT*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)
     1 +CSXCT*CV3*T**2*(64.*WM2*T3-64.*WM2*S13+64.*ZM2*T1-64.*ZM2*S13+6
     1 4.*T1*T3-64.*T1*S13-64.*U1*T3-32.*U1*U3-128.*T3*S13+96.*WZM2)
     1 +CSXCT*CV3*T**2*U*(32.*S13)
      WZST=WZST
     1 +CSXCT*CV3*T**3*(64.*S13)
     1 +CSXCT*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)
      WZST=WZST
     1 +CSXCT*CA3*(-32.*WM2*ZM2*WM4-32.*WM2*ZM2*ZM4-192.*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-128.*ZM2*T1**2*T3-192.*ZM2*T3*WM4-64.*ZM2*S13*WM4-3
     1 20.*T1*T3*WZM2-64.*T1*S13*WZM2-32.*T1**2*WZM2-64.*T1**2*ZM4+64.*
     1 U1*U3*WZM2+64.*U1*S13*WZM2+32.*U1**2*WZM2-128.*WM4*ZM4)
      WZST=WZST
     1 +CSXCT*CA3*T*(128.*WM2*T1*T3+64.*WM2*U1*T3-32.*WM2*U1*U3+128.*WM
     1 2*ZM4+128.*ZM2*T1*T3-64.*ZM2*T1*S13+64.*ZM2*T1**2+32.*ZM2*U1*T3-
     1 32.*ZM2*U1*U3-32.*ZM2*U1*S13-32.*ZM2*U1**2+128.*ZM2*WM4+64.*T1*U
     1 1*T3-32.*T1*U1*U3+192.*T1*WZM2+64.*T1*ZM4+64.*T1**2*T3-32.*U1**2
     1 *U3+192.*T3*WZM2+64.*T3*WM4)
      WZST=WZST
     1 +CSXCT*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 +CSXCT*CA3*T**2*(-64.*WM2*T3+64.*WM2*S13-64.*ZM2*T1+64.*ZM2*S13-
     1 64.*T1*T3+64.*T1*S13-64.*U1*T3+32.*U1*U3-96.*WZM2)
     1 +CSXCT*CA3*T**2*U*(-32.*S13)
     1 +CSXCT*CA3*T**3*(-64.*S13)
      WZST=WZST
     1 +CSXCT*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 +EPF(P1,P2,P3,Q1)*CSXCT*CV3*(64.*ZM2*T3+64.*ZM2*U3+64.*ZM2*S13+6
     1 4.*T3*U3+64.*T3*S13+32.*T3**2+64.*U3*S13+32.*U3**2+32.*ZM4)
      WZST=WZST
     1 +EPF(P1,P2,P3,Q1)*CSXCT*CV3*T*(-16.*ZM2-16.*T3-16.*U3-32.*S13)
     1 +EPF(P1,P2,P3,Q1)*CSXCT*CV3*U*(-16.*ZM2-16.*T3-16.*U3-32.*S13)
     1 +EPF(P1,P2,P3,Q3)*CSXCT*CV3*(32.*WM2*T3+32.*WM2*U3+32.*WZM2)
     1 +EPF(P1,P2,P3,Q3)*CSXCT*CV3*T*(-16.*WM2)
     1 +EPF(P1,P2,P3,Q3)*CSXCT*CV3*U*(-16.*WM2)
     1 +EPF(P1,P3,Q1,Q3)*CSXCT*CV3*(32.*WM2*T3+32.*WM2*U3+64.*ZM2*U1+64
     1 .*U1*T3+64.*U1*U3+32.*WZM2)
      WZST=WZST
     1 +EPF(P1,P3,Q1,Q3)*CSXCT*CV3*T*(-16.*WM2-32.*U1)
     1 +16.*EPF(P1,P3,Q1,Q3)*CSXCT*CV3*T*U
     1 +EPF(P1,P3,Q1,Q3)*CSXCT*CV3*U*(-16.*WM2-32.*ZM2-32.*U1-32.*T3-32
     1 .*U3)
     1 +16.*EPF(P1,P3,Q1,Q3)*CSXCT*CV3*U**2
     1 +EPF(P2,P3,Q1,Q3)*CSXCT*CV3*(-32.*WM2*T3-32.*WM2*U3-64.*ZM2*T1-6
     1 4.*T1*T3-64.*T1*U3-32.*WZM2)
      WZST=WZST
     1 +EPF(P2,P3,Q1,Q3)*CSXCT*CV3*T*(16.*WM2+32.*ZM2+32.*T1+32.*T3+32.
     1 *U3)
     1 -16.*EPF(P2,P3,Q1,Q3)*CSXCT*CV3*T*U
     1 -16.*EPF(P2,P3,Q1,Q3)*CSXCT*CV3*T**2
     1 +EPF(P2,P3,Q1,Q3)*CSXCT*CV3*U*(16.*WM2+32.*T1)
      RETURN
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	FUNCTION WZST(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 T 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 WZST
	10	DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
	11	DOUBLE PRECISION WM4,ZM4,WZM2,CSXCT,EPF
	12	#endif
	13	WM4=WM2**2
	14	ZM4=ZM2**2
	15	WZM2=WM2*ZM2
	16	CSXCT=CS*CT/T
	17	WZST=
	18	1 +CSXCTCV3(32.WM2ZM2WM4+32.WM2ZM2ZM4+256.WM2T1T3*2+19
	19	1 2.WM2T1ZM4-128.WM2U1T3U3+64.WM2T3ZM4+64.WM2S13*ZM4-6
	20	1 4.ZM2T1U1T3-64.ZM2T1U1U3+64.ZM2T1WM4+128.ZM2T12T
	21	1 3+192.ZM2T3WM4+64.ZM2S13WM4-128.T1U1T3U3+320.T1T3*WZ
	22	1 M2+64.T1S13WZM2+128.T1*2T3*2+32.T1*2WZM2+64.T12ZM4
	23	1 -64.U1U3WZM2+64.U1S13WZM2+32.U12WZM2+128.T3S13*WZM2+
	24	1 64.T32WZM2+128.T32WM4+64.S132WZM2+128.WM4ZM4)
	25	WZST=WZST
	26	1 +CSXCTCV3T(-128.WM2T1T3+64.WM2U1T3+32.WM2U1U3+128.*W
	27	1 M2T3S13-128.WM2T3*2-128.WM2ZM4-128.ZM2T1T3+64.ZM2T1*
	28	1 S13-64.ZM2T1*2+32.ZM2U1T3+32.ZM2U1U3-32.ZM2U1S13-32.
	29	1 ZM2U1*2-128.ZM2WM4+64.T1U1T3+32.T1U1U3+128.T1T3S13
	30	1 -128.T1T3*2-192.T1WZM2-64.T1ZM4-64.T1*2T3+64.U1T3*U3
	31	1 -64.U1U3S13-32.U1*2U3-192.T3WZM2-64.T3WM4)
	32	WZST=WZST
	33	1 +CSXCTCV3TU(32.WM2T3-32.WM2S13+32.ZM2T1-32.ZM2S13+32
	34	1 .T1T3-32.T1S13-64.U1T3-32.U1S13-64.T3S13-64.S13*2+32
	35	1 .*WZM2)
	36	1 +CSXCTCV3T*2(64.WM2T3-64.WM2S13+64.ZM2T1-64.ZM2S13+6
	37	1 4.T1T3-64.T1S13-64.U1T3-32.U1U3-128.T3S13+96.*WZM2)
	38	1 +CSXCTCV3T*2U(32.S13)
	39	WZST=WZST
	40	1 +CSXCTCV3T*3(64.*S13)
	41	1 +CSXCTCV3U(-64.WM2T1T3+32.WM2U1T3-64.WM2T3S13-64.*WM
	42	1 2T32-32.WM2ZM4-64.ZM2T1T3-32.ZM2T1S13-32.ZM2T1*2-3
	43	1 2.ZM2WM4+32.T1U1T3-64.T1T3S13-64.T1T3*2-64.T1*WZM2-3
	44	1 2.T1ZM4-32.T12T3-64.T3WZM2-32.T3WM4-32.S13WZM2)
	45	WZST=WZST
	46	1 +CSXCTCA3(-32.WM2ZM2WM4-32.WM2ZM2ZM4-192.WM2T1*ZM4-64.
	47	1 WM2T3ZM4-64.WM2S13ZM4-64.ZM2T1U1T3+64.ZM2T1U1U3-64
	48	1 .ZM2T1WM4-128.ZM2T12T3-192.ZM2T3WM4-64.ZM2S13WM4-3
	49	1 20.T1T3WZM2-64.T1S13WZM2-32.T12WZM2-64.T12ZM4+64.*
	50	1 U1U3WZM2+64.U1S13WZM2+32.U1*2WZM2-128.WM4ZM4)
	51	WZST=WZST
	52	1 +CSXCTCA3T(128.WM2T1T3+64.WM2U1T3-32.WM2U1U3+128.*WM
	53	1 2ZM4+128.ZM2T1T3-64.ZM2T1S13+64.ZM2T12+32.ZM2U1T3-
	54	1 32.ZM2U1U3-32.ZM2U1S13-32.ZM2U1*2+128.ZM2WM4+64.T1*U
	55	1 1T3-32.T1U1U3+192.T1WZM2+64.T1ZM4+64.T12T3-32.U1*2
	56	1 U3+192.T3WZM2+64.T3*WM4)
	57	WZST=WZST
	58	1 +CSXCTCA3TU(-32.WM2T3+32.WM2S13-32.ZM2T1+32.ZM2S13-3
	59	1 2.T1T3+32.T1S13-64.U1T3-32.U1S13-32.*WZM2)
	60	1 +CSXCTCA3T*2(-64.WM2T3+64.WM2S13-64.ZM2T1+64.ZM2S13-
	61	1 64.T1T3+64.T1S13-64.U1T3+32.U1U3-96.*WZM2)
	62	1 +CSXCTCA3T*2U(-32.S13)
	63	1 +CSXCTCA3T*3(-64.*S13)
	64	WZST=WZST
65	1 +CSXCTCA3U(64.WM2T1T3+32.WM2U1T3+32.WM2ZM4+64.ZM2*T1
66	1 T3+32.ZM2T1S13+32.ZM2T1*2+32.ZM2WM4+32.T1U1T3+64.*T1
67	1 WZM2+32.T1ZM4+32.T1*2T3+64.T3WZM2+32.T3WM4+32.S13WZM
68	1 2)
69	1 +EPF(P1,P2,P3,Q1)CSXCTCV3(64.ZM2T3+64.ZM2U3+64.ZM2*S13+6
70	1 4.T3U3+64.T3S13+32.T32+64.U3S13+32.U3*2+32.ZM4)
71	WZST=WZST
72	1 +EPF(P1,P2,P3,Q1)CSXCTCV3T(-16.ZM2-16.T3-16.U3-32.S13)
73	1 +EPF(P1,P2,P3,Q1)CSXCTCV3U(-16.ZM2-16.T3-16.U3-32.S13)
74	1 +EPF(P1,P2,P3,Q3)CSXCTCV3(32.WM2T3+32.WM2U3+32.WZM2)
75	1 +EPF(P1,P2,P3,Q3)CSXCTCV3T(-16.*WM2)
76	1 +EPF(P1,P2,P3,Q3)CSXCTCV3U(-16.*WM2)
77	1 +EPF(P1,P3,Q1,Q3)CSXCTCV3(32.WM2T3+32.WM2U3+64.ZM2*U1+64
78	1 .U1T3+64.U1U3+32.*WZM2)
79	WZST=WZST
80	1 +EPF(P1,P3,Q1,Q3)CSXCTCV3T(-16.WM2-32.U1)
81	1 +16.EPF(P1,P3,Q1,Q3)CSXCTCV3T*U
82	1 +EPF(P1,P3,Q1,Q3)CSXCTCV3U(-16.WM2-32.ZM2-32.U1-32.T3-32
83	1 .*U3)
84	1 +16.EPF(P1,P3,Q1,Q3)CSXCTCV3U**2
85	1 +EPF(P2,P3,Q1,Q3)CSXCTCV3(-32.WM2T3-32.WM2U3-64.ZM2*T1-6
86	1 4.T1T3-64.T1U3-32.*WZM2)
87	WZST=WZST
88	1 +EPF(P2,P3,Q1,Q3)CSXCTCV3T(16.WM2+32.ZM2+32.T1+32.T3+32.
89	1 *U3)
90	1 -16.EPF(P2,P3,Q1,Q3)CSXCTCV3T*U
91	1 -16.EPF(P2,P3,Q1,Q3)CSXCTCV3T**2
92	1 +EPF(P2,P3,Q1,Q3)CSXCTCV3U(16.WM2+32.T1)
93	RETURN
94	END