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

#include "isajet/pilot.h"
      FUNCTION WZSS(T,U,T1,U1,T3,U3,P1,P2)
C          DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980).
C          SQUARE OF S GRAPH.
#include "isajet/itapes.inc"
#include "isajet/wwpar.inc"
      DIMENSION P1(4),P2(4)
#if defined(CERNLIB_DOUBLE)
      DOUBLE PRECISION WZSS
      DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
      DOUBLE PRECISION WM4,ZM4,WZM2,CSXCS
#endif
      WM4=WM2**2
      ZM4=ZM2**2
      WZM2=WM2*ZM2
      CSXCS=CS**2
      WZSS=
     1 +CSXCS*CV3*(-32.*WM2*ZM2*WM4-32.*WM2*ZM2*ZM4-128.*WM2*T1*T3**2-1
     1 28.*WM2*T1*ZM4+128.*WM2*U1*T3*U3-64.*WM2*T3*ZM4-64.*WM2*S13*ZM4+
     1 64.*ZM2*T1*U1*T3+64.*ZM2*T1*U1*U3-64.*ZM2*T1*WM4-64.*ZM2*T1**2*T
     1 3-64.*ZM2*U1**2*U3-128.*ZM2*T3*WM4-64.*ZM2*S13*WM4+128.*T1*U1*T3
     1 *U3-192.*T1*T3*WZM2-64.*T1*S13*WZM2-64.*T1**2*T3**2-32.*T1**2*WZ
     1 M2-32.*T1**2*ZM4+64.*U1*U3*WZM2-64.*U1*S13*WZM2-64.*U1**2*U3**2-
     1 32.*U1**2*WZM2-32.*U1**2*ZM4-128.*T3*S13*WZM2-64.*T3**2*WZM2-64.
     1 *T3**2*WM4-64.*S13**2*WZM2-96.*WM4*ZM4)
      WZSS=WZSS
     1 +CSXCS*CV3*T*(64.*WM2*T1*T3-32.*WM2*U1*T3-32.*WM2*U1*U3-64.*WM2*
     1 T3*S13+64.*WM2*T3**2+96.*WM2*ZM4+64.*ZM2*T1*T3-32.*ZM2*T1*S13+32
     1 .*ZM2*T1**2-32.*ZM2*U1*T3-32.*ZM2*U1*U3+32.*ZM2*U1*S13+32.*ZM2*U
     1 1**2+96.*ZM2*WM4-32.*T1*U1*T3-32.*T1*U1*U3-64.*T1*T3*S13+64.*T1*
     1 T3**2+128.*T1*WZM2+32.*T1*ZM4+32.*T1**2*T3-64.*U1*T3*U3+64.*U1*U
     1 3*S13+32.*U1**2*U3+128.*T3*WZM2+32.*T3*WM4+32.*S13*WZM2)

      WZSS=WZSS
     1 +CSXCS*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+32.*U1*U3+32.*U1*S13+64.*T3*S13+64
     1 .*S13**2-32.*WZM2)
     1 +CSXCS*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+32.*U1*U3+64.*T3*S13-64.*WZM2)
     1 +CSXCS*CV3*T**2*U*(-32.*S13)
     1 +CSXCS*CV3*T**3*(-32.*S13)
      WZSS=WZSS
     1 +CSXCS*CV3*U*(64.*WM2*T1*T3-32.*WM2*U1*T3-32.*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-32.*ZM2*U1*T3-32.*ZM2*U1*U3-32.*ZM2*U1*S13+32.*ZM2*U
     1 1**2+32.*ZM2*WM4-32.*T1*U1*T3-32.*T1*U1*U3+64.*T1*T3*S13+64.*T1*
     1 T3**2+64.*T1*WZM2+32.*T1*ZM4+32.*T1**2*T3-64.*U1*T3*U3-64.*U1*U3
     1 *S13+32.*U1**2*U3+64.*T3*WZM2+32.*T3*WM4+32.*S13*WZM2)
     1 +CSXCS*CV3*U**2*(32.*U1*T3+32.*U1*S13)
      WZSS=WZSS
     1 +CSXCS*CA3*(32.*WM2*ZM2*WM4+32.*WM2*ZM2*ZM4+128.*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*T1**2*T3-64.*ZM2*U1**2*U3+128.*ZM2*T3*WM4+64
     1 .*ZM2*S13*WM4+192.*T1*T3*WZM2+64.*T1*S13*WZM2+32.*T1**2*WZM2+32.
     1 *T1**2*ZM4-64.*U1*U3*WZM2-64.*U1*S13*WZM2-32.*U1**2*WZM2-32.*U1*
     1 *2*ZM4+96.*WM4*ZM4)
      WZSS=WZSS
     1 +CSXCS*CA3*T*(-64.*WM2*T1*T3-32.*WM2*U1*T3+32.*WM2*U1*U3-96.*WM2
     1 *ZM4-64.*ZM2*T1*T3+32.*ZM2*T1*S13-32.*ZM2*T1**2-32.*ZM2*U1*T3+32
     1 .*ZM2*U1*U3+32.*ZM2*U1*S13+32.*ZM2*U1**2-96.*ZM2*WM4-32.*T1*U1*T
     1 3+32.*T1*U1*U3-128.*T1*WZM2-32.*T1*ZM4-32.*T1**2*T3+32.*U1**2*U3
     1 -128.*T3*WZM2-32.*T3*WM4-32.*S13*WZM2)
     1 +CSXCS*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-32.*U1*U3+32.*U1*S13+32.*WZM2)
      WZSS=WZSS
     1 +CSXCS*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.*U1*U3+64.*WZM2)
     1 +CSXCS*CA3*T**2*U*(32.*S13)
     1 +CSXCS*CA3*T**3*(32.*S13)
      WZSS=WZSS
     1 +CSXCS*CA3*U*(-64.*WM2*T1*T3-32.*WM2*U1*T3+32.*WM2*U1*U3-32.*WM2
     1 *ZM4-64.*ZM2*T1*T3-32.*ZM2*T1*S13-32.*ZM2*T1**2-32.*ZM2*U1*T3+32
     1 .*ZM2*U1*U3-32.*ZM2*U1*S13+32.*ZM2*U1**2-32.*ZM2*WM4-32.*T1*U1*T
     1 3+32.*T1*U1*U3-64.*T1*WZM2-32.*T1*ZM4-32.*T1**2*T3+32.*U1**2*U3-
     1 64.*T3*WZM2-32.*T3*WM4-32.*S13*WZM2)
     1 +CSXCS*CA3*U**2*(32.*U1*T3+32.*U1*S13)
      RETURN
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	FUNCTION WZSS(T,U,T1,U1,T3,U3,P1,P2)
	3	C DECAY DISTRIBUTION FOR W- Z0 PAIRS FROM SCHOONSCHIP(1980).
	4	C SQUARE OF S GRAPH.
	5	#include "isajet/itapes.inc"
	6	#include "isajet/wwpar.inc"
	7	DIMENSION P1(4),P2(4)
	8	#if defined(CERNLIB_DOUBLE)
	9	DOUBLE PRECISION WZSS
	10	DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
	11	DOUBLE PRECISION WM4,ZM4,WZM2,CSXCS
	12	#endif
	13	WM4=WM2**2
	14	ZM4=ZM2**2
	15	WZM2=WM2*ZM2
	16	CSXCS=CS**2
	17	WZSS=
	18	1 +CSXCSCV3(-32.WM2ZM2WM4-32.WM2ZM2ZM4-128.WM2T1T3*2-1
	19	1 28.WM2T1ZM4+128.WM2U1T3U3-64.WM2T3ZM4-64.WM2S13*ZM4+
	20	1 64.ZM2T1U1T3+64.ZM2T1U1U3-64.ZM2T1WM4-64.ZM2T12T
	21	1 3-64.ZM2U1*2U3-128.ZM2T3WM4-64.ZM2S13WM4+128.T1U1*T3
	22	1 U3-192.T1T3WZM2-64.T1S13WZM2-64.T1*2T3*2-32.T1*2WZ
	23	1 M2-32.T12ZM4+64.U1U3WZM2-64.U1S13WZM2-64.U12U3**2-
	24	1 32.U12WZM2-32.U12ZM4-128.T3S13WZM2-64.T3*2WZM2-64.
	25	1 T32WM4-64.S132WZM2-96.WM4ZM4)
	26	WZSS=WZSS
	27	1 +CSXCSCV3T(64.WM2T1T3-32.WM2U1T3-32.WM2U1U3-64.WM2
	28	1 T3S13+64.WM2T32+96.WM2ZM4+64.ZM2T1T3-32.ZM2T1*S13+32
	29	1 .ZM2T1*2-32.ZM2U1T3-32.ZM2U1U3+32.ZM2U1S13+32.ZM2U
	30	1 1*2+96.ZM2WM4-32.T1U1T3-32.T1U1U3-64.T1T3S13+64.T1
	31	1 T3*2+128.T1WZM2+32.T1ZM4+32.T1*2T3-64.U1T3U3+64.U1*U
	32	1 3S13+32.U1*2U3+128.T3WZM2+32.T3WM4+32.S13WZM2)
	33
	34	WZSS=WZSS
	35	1 +CSXCSCV3TU(-32.WM2T3+32.WM2S13-32.ZM2T1+32.ZM2S13-3
	36	1 2.T1T3+32.T1S13+64.U1T3+32.U1U3+32.U1S13+64.T3S13+64
	37	1 .S132-32.WZM2)
	38	1 +CSXCSCV3T*2(-32.WM2T3+32.WM2S13-32.ZM2T1+32.ZM2S13-
	39	1 32.T1T3+32.T1S13+32.U1T3+32.U1U3+64.T3S13-64.*WZM2)
	40	1 +CSXCSCV3T*2U(-32.S13)
	41	1 +CSXCSCV3T*3(-32.*S13)
	42	WZSS=WZSS
	43	1 +CSXCSCV3U(64.WM2T1T3-32.WM2U1T3-32.WM2U1U3+64.WM2
	44	1 T3S13+64.WM2T32+32.WM2ZM4+64.ZM2T1T3+32.ZM2T1*S13+32
	45	1 .ZM2T1*2-32.ZM2U1T3-32.ZM2U1U3-32.ZM2U1S13+32.ZM2U
	46	1 1*2+32.ZM2WM4-32.T1U1T3-32.T1U1U3+64.T1T3S13+64.T1
	47	1 T3*2+64.T1WZM2+32.T1ZM4+32.T1*2T3-64.U1T3U3-64.U1*U3
	48	1 S13+32.U1*2U3+64.T3WZM2+32.T3WM4+32.S13WZM2)
	49	1 +CSXCSCV3U*2(32.U1T3+32.U1S13)
	50	WZSS=WZSS
	51	1 +CSXCSCA3(32.WM2ZM2WM4+32.WM2ZM2ZM4+128.WM2T1ZM4+64.
	52	1 WM2T3ZM4+64.WM2S13ZM4+64.ZM2T1U1T3-64.ZM2T1U1*U3+64.
	53	1 ZM2T1WM4+64.ZM2T12T3-64.ZM2U1*2U3+128.ZM2T3*WM4+64
	54	1 .ZM2S13WM4+192.T1T3WZM2+64.T1S13WZM2+32.T1*2WZM2+32.
	55	1 T12ZM4-64.U1U3WZM2-64.U1S13WZM2-32.U12WZM2-32.U1
	56	1 2ZM4+96.WM4ZM4)
	57	WZSS=WZSS
	58	1 +CSXCSCA3T(-64.WM2T1T3-32.WM2U1T3+32.WM2U1U3-96.*WM2
	59	1 ZM4-64.ZM2T1T3+32.ZM2T1S13-32.ZM2T12-32.ZM2U1T3+32
	60	1 .ZM2U1U3+32.ZM2U1S13+32.ZM2U1*2-96.ZM2WM4-32.T1U1T
	61	1 3+32.T1U1U3-128.T1WZM2-32.T1ZM4-32.T1*2T3+32.U12U3
	62	1 -128.T3WZM2-32.T3WM4-32.S13WZM2)
	63	1 +CSXCSCA3TU(32.WM2T3-32.WM2S13+32.ZM2T1-32.ZM2S13+32
	64	1 .T1T3-32.T1S13+64.U1T3-32.U1U3+32.U1S13+32.*WZM2)
65	WZSS=WZSS
66	1 +CSXCSCA3T*2(32.WM2T3-32.WM2S13+32.ZM2T1-32.ZM2S13+3
67	1 2.T1T3-32.T1S13+32.U1T3-32.U1U3+64.*WZM2)
68	1 +CSXCSCA3T*2U(32.S13)
69	1 +CSXCSCA3T*3(32.*S13)
70	WZSS=WZSS
71	1 +CSXCSCA3U(-64.WM2T1T3-32.WM2U1T3+32.WM2U1U3-32.*WM2
72	1 ZM4-64.ZM2T1T3-32.ZM2T1S13-32.ZM2T12-32.ZM2U1T3+32
73	1 .ZM2U1U3-32.ZM2U1S13+32.ZM2U1*2-32.ZM2WM4-32.T1U1T
74	1 3+32.T1U1U3-64.T1WZM2-32.T1ZM4-32.T1*2T3+32.U12U3-
75	1 64.T3WZM2-32.T3WM4-32.S13WZM2)
76	1 +CSXCSCA3U*2(32.U1T3+32.U1S13)
77	RETURN
78	END