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

#include "isajet/pilot.h"
      FUNCTION ZZALL(T,U,T1,U1,T3,U3,P1,P2)
C          DECAY DISTRIBUTION FOR Z0 Z0 PAIRS FROM SCHOONSCHIP(1980).
C          INCLUDES TT, TU, AND UU TERMS.
#include "isajet/itapes.inc"
#include "isajet/wwpar.inc"
      DIMENSION P1(4),P2(4)
#if defined(CERNLIB_DOUBLE)
      DOUBLE PRECISION ZZALL
      DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
      DOUBLE PRECISION CVAVA,DVAVA,CVVVV,CVVAA,TT,TU,UU,EPF
#endif
      TTUU=(T*U)**2
      CVAVA=CV*CA*(CV1*CA3+CA1*CV3)/TTUU
      DVAVA=CV*CA*(CV1*CA3-CA1*CV3)/TTUU
      CVVVV=(CV**2+CA**2)*CV1*CV3/TTUU
      CVVAA=(CV**2+CA**2)*CA1*CA3/TTUU
      TT=T**2
      TU=T*U
      UU=U**2
      ZZALL=
     1 +CVAVA*T*U*(-16.*ZM2*T1*U1*T3+16.*ZM2*T1*U1*U3+16.*ZM2*T1*T3*U3-
     1 16.*ZM2*U1*T3*U3-16.*ZM2**2*T1*S13-8.*ZM2**2*T1**2+16.*ZM2**2*U1
     1 *S13+8.*ZM2**2*U1**2+16.*ZM2**2*T3*S13+8.*ZM2**2*T3**2-16.*ZM2**
     1 2*U3*S13-8.*ZM2**2*U3**2)
     1 +CVAVA*T*U**2*(16.*ZM2*T1*S13-16.*ZM2*T3*S13)
     1 +CVAVA*T**2*(16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-8.*ZM2**2*U1**2+8
     1 .*ZM2**2*U3**2)
     1 +CVAVA*T**2*U*(-16.*ZM2*U1*S13+16.*ZM2*U3*S13)
     1 +CVAVA*U**2*(-16.*ZM2*T1*T3**2+16.*ZM2*T1**2*T3+8.*ZM2**2*T1**2-
     1 8.*ZM2**2*T3**2)
      ZZALL=ZZALL
     1 +DVAVA*T*U*(-16.*ZM2*T1*U1*T3+16.*ZM2*T1*U1*U3-16.*ZM2*T1*T3*U3+
     1 16.*ZM2*U1*T3*U3-32.*ZM2**2*T1*T3-16.*ZM2**2*T1*S13-8.*ZM2**2*T1
     1 **2+32.*ZM2**2*U1*U3+16.*ZM2**2*U1*S13+8.*ZM2**2*U1**2-16.*ZM2**
     1 2*T3*S13-8.*ZM2**2*T3**2+16.*ZM2**2*U3*S13+8.*ZM2**2*U3**2-32.*Z
     1 M2**3*T1+32.*ZM2**3*U1-32.*ZM2**3*T3+32.*ZM2**3*U3)
     1 +DVAVA*T*U**2*(16.*ZM2*T1*S13+16.*ZM2*T3*S13-16.*ZM2**2*U1-16.*Z
     1 M2**2*U3-32.*ZM2**3)
     1 +DVAVA*T*U**3*(8.*ZM2**2)
      ZZALL=ZZALL
     1 +DVAVA*T**2*(-16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-32.*ZM2**2*U1*U3
     1 -8.*ZM2**2*U1**2-8.*ZM2**2*U3**2-16.*ZM2**3*U1-16.*ZM2**3*U3-8.*
     1 ZM2**4)
     1 +DVAVA*T**2*U*(-16.*ZM2*U1*S13-16.*ZM2*U3*S13+16.*ZM2**2*T1+16.*
     1 ZM2**2*T3+32.*ZM2**3)
     1 +DVAVA*T**3*U*(-8.*ZM2**2)
     1 +DVAVA*U**2*(16.*ZM2*T1*T3**2+16.*ZM2*T1**2*T3+32.*ZM2**2*T1*T3+
     1 8.*ZM2**2*T1**2+8.*ZM2**2*T3**2+16.*ZM2**3*T1+16.*ZM2**3*T3+8.*Z
     1 M2**4)
      ZZALL=ZZALL
     1 +CVVVV*T*U*(-16.*ZM2*T1*U1*T3-16.*ZM2*T1*U1*U3-16.*ZM2*T1*T3*U3-
     1 16.*ZM2*U1*T3*U3+16.*ZM2**2*T1*S13+8.*ZM2**2*T1**2+16.*ZM2**2*U1
     1 *S13+8.*ZM2**2*U1**2+16.*ZM2**2*T3*S13+8.*ZM2**2*T3**2+16.*ZM2**
     1 2*U3*S13+8.*ZM2**2*U3**2+16.*ZM2**2*S13**2+16.*ZM2**3*T1+16.*ZM2
     1 **3*U1+16.*ZM2**3*T3+16.*ZM2**3*U3+32.*ZM2**3*S13+32.*ZM2**4-32.
     1 *T1*U1*T3*U3)
     1 +CVVVV*T*U**2*(-16.*ZM2*T1*S13-16.*ZM2*T3*S13-8.*ZM2**2*U1-8.*ZM
     1 2**2*U3-16.*ZM2**2*S13-16.*ZM2**3-32.*T1*T3*S13)
      ZZALL=ZZALL
     1 +CVVVV*T*U**3*(4.*ZM2**2)
     1 +CVVVV*T**2*(-16.*ZM2*U1*U3**2-16.*ZM2*U1**2*U3-16.*ZM2**2*U1*U3
     1 -8.*ZM2**2*U1**2-8.*ZM2**2*U3**2-8.*ZM2**3*U1-8.*ZM2**3*U3-4.*ZM
     1 2**4-16.*U1**2*U3**2)
     1 +CVVVV*T**2*U*(-16.*ZM2*U1*S13-16.*ZM2*U3*S13-8.*ZM2**2*T1-8.*ZM
     1 2**2*T3-16.*ZM2**2*S13-16.*ZM2**3-32.*U1*U3*S13)
     1 +CVVVV*T**2*U**2*(-16.*S13**2)
     1 +CVVVV*T**3*U*(4.*ZM2**2)
      ZZALL=ZZALL
     1 +CVVVV*U**2*(-16.*ZM2*T1*T3**2-16.*ZM2*T1**2*T3-16.*ZM2**2*T1*T3
     1 -8.*ZM2**2*T1**2-8.*ZM2**2*T3**2-8.*ZM2**3*T1-8.*ZM2**3*T3-4.*ZM
     1 2**4-16.*T1**2*T3**2)
     1 +CVVAA*T*U*(-16.*ZM2**3*T1-16.*ZM2**3*U1-16.*ZM2**3*T3-16.*ZM2**
     1 3*U3-32.*ZM2**3*S13-32.*ZM2**4)
     1 +CVVAA*T*U**2*(8.*ZM2**2*U1+8.*ZM2**2*U3+16.*ZM2**2*S13+16.*ZM2*
     1 *3)
     1 +CVVAA*T*U**3*(-4.*ZM2**2)
      ZZALL=ZZALL
     1 +CVVAA*T**2*(16.*ZM2**2*U1*U3+8.*ZM2**3*U1+8.*ZM2**3*U3+4.*ZM2**
     1 4)
     1 +CVVAA*T**2*U*(8.*ZM2**2*T1+8.*ZM2**2*T3+16.*ZM2**2*S13+16.*ZM2*
     1 *3)
     1 +CVVAA*T**3*U*(-4.*ZM2**2)
     1 +CVVAA*U**2*(16.*ZM2**2*T1*T3+8.*ZM2**3*T1+8.*ZM2**3*T3+4.*ZM2**
     1 4)
C          THE EPF TERMS FROM SCHOONSCHIP DID NOT CONTAIN CVAVA, ETC.,
C          BUT DID CONTAIN AN EXPLICIT T*U. THIS WAS REPLACED BY HAND
C          BY 1./TU.
      ZZALL=ZZALL
     1 +EPF(P1,P2,P3,Q3)/TU*(-32.*CV*CA*CV1*CV3*T1*S13)
     1 +EPF(P1,P2,Q1,Q3)/TU*(-32.*ZM2*CV*CA*CV1*CV3*T1-32.*CV*CA*CV1*C
     1 V3*T1*T3-32.*CV*CA*CV1*CV3*T1*U3)
     1 +EPF(P1,P3,Q1,Q3)/TU*(32.*CV*CA*CV1*CV3*T1*T3)
     1 +EPF(P2,P3,Q1,Q3)/TU*(-32.*CV*CA*CV1*CV3*T1*U3)+0.
      ZZALL=2.*ZZALL
      RETURN
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	FUNCTION ZZALL(T,U,T1,U1,T3,U3,P1,P2)
	3	C DECAY DISTRIBUTION FOR Z0 Z0 PAIRS FROM SCHOONSCHIP(1980).
	4	C INCLUDES TT, TU, AND UU 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 ZZALL
	10	DOUBLE PRECISION T,U,T1,U1,T3,U3,P1,P2
	11	DOUBLE PRECISION CVAVA,DVAVA,CVVVV,CVVAA,TT,TU,UU,EPF
	12	#endif
	13	TTUU=(TU)*2
	14	CVAVA=CVCA(CV1CA3+CA1CV3)/TTUU
	15	DVAVA=CVCA(CV1CA3-CA1CV3)/TTUU
	16	CVVVV=(CV2+CA2)CV1CV3/TTUU
	17	CVVAA=(CV2+CA2)CA1CA3/TTUU
	18	TT=T**2
	19	TU=T*U
	20	UU=U**2
	21	ZZALL=
	22	1 +CVAVATU(-16.ZM2T1U1T3+16.ZM2T1U1U3+16.ZM2T1T3*U3-
	23	1 16.ZM2U1T3U3-16.ZM22T1S13-8.ZM2*2T1*2+16.ZM2*2U1
	24	1 S13+8.ZM2*2U1*2+16.ZM2*2T3S13+8.ZM2*2T3*2-16.ZM2**
	25	1 2U3S13-8.ZM22U3**2)
	26	1 +CVAVATU*2(16.ZM2T1S13-16.ZM2T3S13)
	27	1 +CVAVAT2(16.ZM2U1U32-16.ZM2U12U3-8.ZM22U1**2+8
	28	1 .ZM22U3**2)
	29	1 +CVAVAT2U(-16.ZM2U1S13+16.ZM2U3*S13)
	30	1 +CVAVAU2(-16.ZM2T1T32+16.ZM2T12T3+8.ZM22T1**2-
	31	1 8.ZM22T3**2)
	32	ZZALL=ZZALL
	33	1 +DVAVATU(-16.ZM2T1U1T3+16.ZM2T1U1U3-16.ZM2T1T3*U3+
	34	1 16.ZM2U1T3U3-32.ZM22T1T3-16.ZM2*2T1S13-8.ZM2*2T1
	35	1 *2+32.ZM2*2U1U3+16.ZM2*2U1S13+8.ZM2*2U1*2-16.ZM2**
	36	1 2T3S13-8.ZM22T3*2+16.ZM2*2U3S13+8.ZM2*2U3*2-32.Z
	37	1 M2*3T1+32.ZM23U1-32.ZM23T3+32.ZM23U3)
	38	1 +DVAVATU*2(16.ZM2T1S13+16.ZM2T3S13-16.ZM22U1-16.*Z
	39	1 M2*2U3-32.ZM2*3)
	40	1 +DVAVATU*3(8.ZM2*2)
	41	ZZALL=ZZALL
	42	1 +DVAVAT2(-16.ZM2U1U32-16.ZM2U12U3-32.ZM22U1*U3
	43	1 -8.ZM22U1*2-8.ZM2*2U3*2-16.ZM2*3U1-16.ZM23U3-8.*
	44	1 ZM2**4)
	45	1 +DVAVAT2U(-16.ZM2U1S13-16.ZM2U3S13+16.ZM2*2T1+16.*
	46	1 ZM2*2T3+32.ZM2*3)
	47	1 +DVAVAT3U(-8.ZM2**2)
	48	1 +DVAVAU2(16.ZM2T1T32+16.ZM2T12T3+32.ZM22T1*T3+
	49	1 8.ZM22T1*2+8.ZM2*2T3*2+16.ZM2*3T1+16.ZM23T3+8.*Z
	50	1 M2**4)
	51	ZZALL=ZZALL
	52	1 +CVVVVTU(-16.ZM2T1U1T3-16.ZM2T1U1U3-16.ZM2T1T3*U3-
	53	1 16.ZM2U1T3U3+16.ZM22T1S13+8.ZM2*2T1*2+16.ZM2*2U1
	54	1 S13+8.ZM2*2U1*2+16.ZM2*2T3S13+8.ZM2*2T3*2+16.ZM2**
	55	1 2U3S13+8.ZM22U3*2+16.ZM2*2S13*2+16.ZM2*3T1+16.*ZM2
	56	1 *3U1+16.ZM23T3+16.ZM23U3+32.ZM23S13+32.ZM2*4-32.
	57	1 T1U1T3U3)
	58	1 +CVVVVTU*2(-16.ZM2T1S13-16.ZM2T3S13-8.ZM22U1-8.*ZM
	59	1 2*2U3-16.ZM22S13-16.ZM23-32.T1T3S13)
	60	ZZALL=ZZALL
	61	1 +CVVVVTU*3(4.ZM2*2)
	62	1 +CVVVVT2(-16.ZM2U1U32-16.ZM2U12U3-16.ZM22U1*U3
	63	1 -8.ZM22U1*2-8.ZM2*2U3*2-8.ZM2*3U1-8.ZM23U3-4.*ZM
	64	1 2*4-16.U1*2U3**2)
65	1 +CVVVVT2U(-16.ZM2U1S13-16.ZM2U3S13-8.ZM2*2T1-8.*ZM
66	1 2*2T3-16.ZM22S13-16.ZM23-32.U1U3S13)
67	1 +CVVVVT2U*2(-16.S13*2)
68	1 +CVVVVT3U(4.ZM2**2)
69	ZZALL=ZZALL
70	1 +CVVVVU2(-16.ZM2T1T32-16.ZM2T12T3-16.ZM22T1*T3
71	1 -8.ZM22T1*2-8.ZM2*2T3*2-8.ZM2*3T1-8.ZM23T3-4.*ZM
72	1 2*4-16.T1*2T3**2)
73	1 +CVVAATU(-16.ZM2*3T1-16.ZM23U1-16.ZM23T3-16.ZM2*
74	1 3U3-32.ZM2*3S13-32.ZM2*4)
75	1 +CVVAATU*2(8.ZM22U1+8.ZM22U3+16.ZM22S13+16.ZM2
76	1 *3)
77	1 +CVVAATU*3(-4.ZM2*2)
78	ZZALL=ZZALL
79	1 +CVVAAT2(16.ZM22U1U3+8.ZM2*3U1+8.ZM23U3+4.ZM2*
80	1 4)
81	1 +CVVAAT2U(8.ZM2*2T1+8.ZM22T3+16.ZM22S13+16.ZM2
82	1 *3)
83	1 +CVVAAT3U(-4.ZM2**2)
84	1 +CVVAAU2(16.ZM22T1T3+8.ZM2*3T1+8.ZM23T3+4.ZM2*
85	1 4)
86	C THE EPF TERMS FROM SCHOONSCHIP DID NOT CONTAIN CVAVA, ETC.,
87	C BUT DID CONTAIN AN EXPLICIT T*U. THIS WAS REPLACED BY HAND
88	C BY 1./TU.
89	ZZALL=ZZALL
90	1 +EPF(P1,P2,P3,Q3)/TU(-32.CVCACV1CV3T1*S13)
91	1 +EPF(P1,P2,Q1,Q3)/TU(-32.ZM2CVCACV1CV3T1-32.CVCACV1*C
92	1 V3T1T3-32.CVCACV1CV3T1U3)
93	1 +EPF(P1,P3,Q1,Q3)/TU(32.CVCACV1CV3T1*T3)
94	1 +EPF(P2,P3,Q1,Q3)/TU(-32.CVCACV1CV3T1*U3)+0.
95	ZZALL=2.*ZZALL
96	RETURN
97	END