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

#include "isajet/pilot.h"
      REAL*8 FUNCTION SZJJ2(P1, P2, P3, P4, P5, IM)
C
C Function generated by Madgraph + hand coding
C Returns amplitude squared summed/ave over colors
C for the point in phase space P1,P2,P3,P4,P5
C and helicity NHEL(1)...
C for the process: g g -> z q(im) qb(im)
C with Madgraph code IM
C
#if defined(CERNLIB_IMPNONE)
      IMPLICIT NONE
#endif
C  
C CONSTANTS
C  
      INTEGER    NEXTERNAL,   NCOMB                     
      PARAMETER (NEXTERNAL=5, NCOMB= 48)
C  
C ARGUMENTS 
C  
      REAL*8 P1(0:3),P2(0:3),P3(0:3),P4(0:3),P5(0:3)
      INTEGER IM                             
C  
C LOCAL VARIABLES 
C  
      INTEGER NHEL(NEXTERNAL,NCOMB),NTRY                                         
      REAL*8 T
      REAL*8 ZJJ2                                                             
      INTEGER IHEL
      LOGICAL GOODHEL(NCOMB)
      DATA GOODHEL/NCOMB*.FALSE./
      DATA NTRY/0/
      DATA (NHEL(IHEL,  1),IHEL=1,5) / -1, -1, -1, -1, -1/
      DATA (NHEL(IHEL,  2),IHEL=1,5) / -1, -1, -1, -1,  1/
      DATA (NHEL(IHEL,  3),IHEL=1,5) / -1, -1, -1,  1, -1/
      DATA (NHEL(IHEL,  4),IHEL=1,5) / -1, -1, -1,  1,  1/
      DATA (NHEL(IHEL,  5),IHEL=1,5) / -1, -1,  0, -1, -1/
      DATA (NHEL(IHEL,  6),IHEL=1,5) / -1, -1,  0, -1,  1/
      DATA (NHEL(IHEL,  7),IHEL=1,5) / -1, -1,  0,  1, -1/
      DATA (NHEL(IHEL,  8),IHEL=1,5) / -1, -1,  0,  1,  1/
      DATA (NHEL(IHEL,  9),IHEL=1,5) / -1, -1,  1, -1, -1/
      DATA (NHEL(IHEL, 10),IHEL=1,5) / -1, -1,  1, -1,  1/
      DATA (NHEL(IHEL, 11),IHEL=1,5) / -1, -1,  1,  1, -1/
      DATA (NHEL(IHEL, 12),IHEL=1,5) / -1, -1,  1,  1,  1/
      DATA (NHEL(IHEL, 13),IHEL=1,5) / -1,  1, -1, -1, -1/
      DATA (NHEL(IHEL, 14),IHEL=1,5) / -1,  1, -1, -1,  1/
      DATA (NHEL(IHEL, 15),IHEL=1,5) / -1,  1, -1,  1, -1/
      DATA (NHEL(IHEL, 16),IHEL=1,5) / -1,  1, -1,  1,  1/
      DATA (NHEL(IHEL, 17),IHEL=1,5) / -1,  1,  0, -1, -1/
      DATA (NHEL(IHEL, 18),IHEL=1,5) / -1,  1,  0, -1,  1/
      DATA (NHEL(IHEL, 19),IHEL=1,5) / -1,  1,  0,  1, -1/
      DATA (NHEL(IHEL, 20),IHEL=1,5) / -1,  1,  0,  1,  1/
      DATA (NHEL(IHEL, 21),IHEL=1,5) / -1,  1,  1, -1, -1/
      DATA (NHEL(IHEL, 22),IHEL=1,5) / -1,  1,  1, -1,  1/
      DATA (NHEL(IHEL, 23),IHEL=1,5) / -1,  1,  1,  1, -1/
      DATA (NHEL(IHEL, 24),IHEL=1,5) / -1,  1,  1,  1,  1/
      DATA (NHEL(IHEL, 25),IHEL=1,5) /  1, -1, -1, -1, -1/
      DATA (NHEL(IHEL, 26),IHEL=1,5) /  1, -1, -1, -1,  1/
      DATA (NHEL(IHEL, 27),IHEL=1,5) /  1, -1, -1,  1, -1/
      DATA (NHEL(IHEL, 28),IHEL=1,5) /  1, -1, -1,  1,  1/
      DATA (NHEL(IHEL, 29),IHEL=1,5) /  1, -1,  0, -1, -1/
      DATA (NHEL(IHEL, 30),IHEL=1,5) /  1, -1,  0, -1,  1/
      DATA (NHEL(IHEL, 31),IHEL=1,5) /  1, -1,  0,  1, -1/
      DATA (NHEL(IHEL, 32),IHEL=1,5) /  1, -1,  0,  1,  1/
      DATA (NHEL(IHEL, 33),IHEL=1,5) /  1, -1,  1, -1, -1/
      DATA (NHEL(IHEL, 34),IHEL=1,5) /  1, -1,  1, -1,  1/
      DATA (NHEL(IHEL, 35),IHEL=1,5) /  1, -1,  1,  1, -1/
      DATA (NHEL(IHEL, 36),IHEL=1,5) /  1, -1,  1,  1,  1/
      DATA (NHEL(IHEL, 37),IHEL=1,5) /  1,  1, -1, -1, -1/
      DATA (NHEL(IHEL, 38),IHEL=1,5) /  1,  1, -1, -1,  1/
      DATA (NHEL(IHEL, 39),IHEL=1,5) /  1,  1, -1,  1, -1/
      DATA (NHEL(IHEL, 40),IHEL=1,5) /  1,  1, -1,  1,  1/
      DATA (NHEL(IHEL, 41),IHEL=1,5) /  1,  1,  0, -1, -1/
      DATA (NHEL(IHEL, 42),IHEL=1,5) /  1,  1,  0, -1,  1/
      DATA (NHEL(IHEL, 43),IHEL=1,5) /  1,  1,  0,  1, -1/
      DATA (NHEL(IHEL, 44),IHEL=1,5) /  1,  1,  0,  1,  1/
      DATA (NHEL(IHEL, 45),IHEL=1,5) /  1,  1,  1, -1, -1/
      DATA (NHEL(IHEL, 46),IHEL=1,5) /  1,  1,  1, -1,  1/
      DATA (NHEL(IHEL, 47),IHEL=1,5) /  1,  1,  1,  1, -1/
      DATA (NHEL(IHEL, 48),IHEL=1,5) /  1,  1,  1,  1,  1/
C ----------
C BEGIN CODE
C ----------
      SZJJ2 = 0D0
      NTRY=NTRY+1
      DO IHEL=1,NCOMB
          IF (GOODHEL(IHEL) .OR. NTRY .LT. 10) THEN
             T=ZJJ2(P1, P2, P3, P4, P5,NHEL(1,IHEL),IM) 
             SZJJ2 = SZJJ2 + T
              IF (T .GT. 0D0 .AND. .NOT. GOODHEL(IHEL)) THEN
                GOODHEL(IHEL)=.TRUE.
              ENDIF
          ENDIF
      ENDDO
      SZJJ2 = SZJJ2 /  4D0 
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	REAL*8 FUNCTION SZJJ2(P1, P2, P3, P4, P5, IM)
	3	C
	4	C Function generated by Madgraph + hand coding
	5	C Returns amplitude squared summed/ave over colors
	6	C for the point in phase space P1,P2,P3,P4,P5
	7	C and helicity NHEL(1)...
	8	C for the process: g g -> z q(im) qb(im)
	9	C with Madgraph code IM
	10	C
	11	#if defined(CERNLIB_IMPNONE)
	12	IMPLICIT NONE
	13	#endif
	14	C
	15	C CONSTANTS
	16	C
	17	INTEGER NEXTERNAL, NCOMB
	18	PARAMETER (NEXTERNAL=5, NCOMB= 48)
	19	C
	20	C ARGUMENTS
	21	C
	22	REAL*8 P1(0:3),P2(0:3),P3(0:3),P4(0:3),P5(0:3)
	23	INTEGER IM
	24	C
	25	C LOCAL VARIABLES
	26	C
	27	INTEGER NHEL(NEXTERNAL,NCOMB),NTRY
	28	REAL*8 T
	29	REAL*8 ZJJ2
	30	INTEGER IHEL
	31	LOGICAL GOODHEL(NCOMB)
	32	DATA GOODHEL/NCOMB*.FALSE./
	33	DATA NTRY/0/
	34	DATA (NHEL(IHEL, 1),IHEL=1,5) / -1, -1, -1, -1, -1/
	35	DATA (NHEL(IHEL, 2),IHEL=1,5) / -1, -1, -1, -1, 1/
	36	DATA (NHEL(IHEL, 3),IHEL=1,5) / -1, -1, -1, 1, -1/
	37	DATA (NHEL(IHEL, 4),IHEL=1,5) / -1, -1, -1, 1, 1/
	38	DATA (NHEL(IHEL, 5),IHEL=1,5) / -1, -1, 0, -1, -1/
	39	DATA (NHEL(IHEL, 6),IHEL=1,5) / -1, -1, 0, -1, 1/
	40	DATA (NHEL(IHEL, 7),IHEL=1,5) / -1, -1, 0, 1, -1/
	41	DATA (NHEL(IHEL, 8),IHEL=1,5) / -1, -1, 0, 1, 1/
	42	DATA (NHEL(IHEL, 9),IHEL=1,5) / -1, -1, 1, -1, -1/
	43	DATA (NHEL(IHEL, 10),IHEL=1,5) / -1, -1, 1, -1, 1/
	44	DATA (NHEL(IHEL, 11),IHEL=1,5) / -1, -1, 1, 1, -1/
	45	DATA (NHEL(IHEL, 12),IHEL=1,5) / -1, -1, 1, 1, 1/
	46	DATA (NHEL(IHEL, 13),IHEL=1,5) / -1, 1, -1, -1, -1/
	47	DATA (NHEL(IHEL, 14),IHEL=1,5) / -1, 1, -1, -1, 1/
	48	DATA (NHEL(IHEL, 15),IHEL=1,5) / -1, 1, -1, 1, -1/
	49	DATA (NHEL(IHEL, 16),IHEL=1,5) / -1, 1, -1, 1, 1/
	50	DATA (NHEL(IHEL, 17),IHEL=1,5) / -1, 1, 0, -1, -1/
	51	DATA (NHEL(IHEL, 18),IHEL=1,5) / -1, 1, 0, -1, 1/
	52	DATA (NHEL(IHEL, 19),IHEL=1,5) / -1, 1, 0, 1, -1/
	53	DATA (NHEL(IHEL, 20),IHEL=1,5) / -1, 1, 0, 1, 1/
	54	DATA (NHEL(IHEL, 21),IHEL=1,5) / -1, 1, 1, -1, -1/
	55	DATA (NHEL(IHEL, 22),IHEL=1,5) / -1, 1, 1, -1, 1/
	56	DATA (NHEL(IHEL, 23),IHEL=1,5) / -1, 1, 1, 1, -1/
	57	DATA (NHEL(IHEL, 24),IHEL=1,5) / -1, 1, 1, 1, 1/
	58	DATA (NHEL(IHEL, 25),IHEL=1,5) / 1, -1, -1, -1, -1/
	59	DATA (NHEL(IHEL, 26),IHEL=1,5) / 1, -1, -1, -1, 1/
	60	DATA (NHEL(IHEL, 27),IHEL=1,5) / 1, -1, -1, 1, -1/
	61	DATA (NHEL(IHEL, 28),IHEL=1,5) / 1, -1, -1, 1, 1/
	62	DATA (NHEL(IHEL, 29),IHEL=1,5) / 1, -1, 0, -1, -1/
	63	DATA (NHEL(IHEL, 30),IHEL=1,5) / 1, -1, 0, -1, 1/
	64	DATA (NHEL(IHEL, 31),IHEL=1,5) / 1, -1, 0, 1, -1/
65	DATA (NHEL(IHEL, 32),IHEL=1,5) / 1, -1, 0, 1, 1/
66	DATA (NHEL(IHEL, 33),IHEL=1,5) / 1, -1, 1, -1, -1/
67	DATA (NHEL(IHEL, 34),IHEL=1,5) / 1, -1, 1, -1, 1/
68	DATA (NHEL(IHEL, 35),IHEL=1,5) / 1, -1, 1, 1, -1/
69	DATA (NHEL(IHEL, 36),IHEL=1,5) / 1, -1, 1, 1, 1/
70	DATA (NHEL(IHEL, 37),IHEL=1,5) / 1, 1, -1, -1, -1/
71	DATA (NHEL(IHEL, 38),IHEL=1,5) / 1, 1, -1, -1, 1/
72	DATA (NHEL(IHEL, 39),IHEL=1,5) / 1, 1, -1, 1, -1/
73	DATA (NHEL(IHEL, 40),IHEL=1,5) / 1, 1, -1, 1, 1/
74	DATA (NHEL(IHEL, 41),IHEL=1,5) / 1, 1, 0, -1, -1/
75	DATA (NHEL(IHEL, 42),IHEL=1,5) / 1, 1, 0, -1, 1/
76	DATA (NHEL(IHEL, 43),IHEL=1,5) / 1, 1, 0, 1, -1/
77	DATA (NHEL(IHEL, 44),IHEL=1,5) / 1, 1, 0, 1, 1/
78	DATA (NHEL(IHEL, 45),IHEL=1,5) / 1, 1, 1, -1, -1/
79	DATA (NHEL(IHEL, 46),IHEL=1,5) / 1, 1, 1, -1, 1/
80	DATA (NHEL(IHEL, 47),IHEL=1,5) / 1, 1, 1, 1, -1/
81	DATA (NHEL(IHEL, 48),IHEL=1,5) / 1, 1, 1, 1, 1/
82	C ----------
83	C BEGIN CODE
84	C ----------
85	SZJJ2 = 0D0
86	NTRY=NTRY+1
87	DO IHEL=1,NCOMB
88	IF (GOODHEL(IHEL) .OR. NTRY .LT. 10) THEN
89	T=ZJJ2(P1, P2, P3, P4, P5,NHEL(1,IHEL),IM)
90	SZJJ2 = SZJJ2 + T
91	IF (T .GT. 0D0 .AND. .NOT. GOODHEL(IHEL)) THEN
92	GOODHEL(IHEL)=.TRUE.
93	ENDIF
94	ENDIF
95	ENDDO
96	SZJJ2 = SZJJ2 / 4D0
97	END