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

#include "isajet/pilot.h"
      FUNCTION WWTT(T,U,T1,U1,T3,U3)
C          DECAY DISTRIBUTION FOR W+ W- PAIRS FROM SCHOONSCHIP(1980).
C          SQUARE OF T GRAPH.
#include "isajet/itapes.inc"
#include "isajet/wwpar.inc"
#if defined(CERNLIB_DOUBLE)
      DOUBLE PRECISION WWTT
      DOUBLE PRECISION T,U,T1,U1,T3,U3
#endif
      WWTT=
     1(+CQ**2*(-512.*T1**2*T3**2)
     1 +CQ**2*T*(-512.*T1*U1*T3-512.*T1*T3*S13+512.*T1*T3**2)
     1 +CQ**2*T*WM2*(-512.*U1*T3-512.*T3*S13+512.*T3**2)
     1 +CQ**2*T**2*(512.*U1*T3+512.*T3*S13)
     1 +CQ**2*WM2*(-1024.*T1*T3**2)
     1 +CQ**2*WM2**2*(-512.*T3**2))/T**2
      WWTT=2.*WWTT
      RETURN
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	FUNCTION WWTT(T,U,T1,U1,T3,U3)
	3	C DECAY DISTRIBUTION FOR W+ W- PAIRS FROM SCHOONSCHIP(1980).
	4	C SQUARE OF T GRAPH.
	5	#include "isajet/itapes.inc"
	6	#include "isajet/wwpar.inc"
	7	#if defined(CERNLIB_DOUBLE)
	8	DOUBLE PRECISION WWTT
	9	DOUBLE PRECISION T,U,T1,U1,T3,U3
	10	#endif
	11	WWTT=
	12	1(+CQ*2(-512.T12T3**2)
	13	1 +CQ*2T(-512.T1U1T3-512.T1T3S13+512.T1T3*2)
	14	1 +CQ*2TWM2(-512.U1T3-512.T3S13+512.T3*2)
	15	1 +CQ*2T*2(512.U1T3+512.T3S13)
	16	1 +CQ*2WM2(-1024.T1T3*2)
	17	1 +CQ*2WM2*2(-512.T32))/T*2
	18	WWTT=2.*WWTT
	19	RETURN
	20	END