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

#include "isajet/pilot.h"
      LOGICAL FUNCTION SETTYP(LPRT)
C
C          Set JETTYPE flags and WMODE flags for WPAIR.
C          Set WMODES and ZMODES flags for secondary W+- and Z0.
C          Return .FALSE. if no error, .TRUE. otherwise.
C
C          Ver 7.18: Initialize all GOQ to false (limit = MXGOQ)
C                    Use LISTSS for Higgs if GOMSSM
C          Ver 7.29: SUSY Higgs decays are done in SETHSS and SIGHSS
C                    using LISTSS order, so SUSY list should be used.
C                    I.e., 7.18 fix was wrong.
C
#if defined(CERNLIB_IMPNONE)
      IMPLICIT NONE
#endif
#include "isajet/itapes.inc"
#include "isajet/keys.inc"
#include "isajet/types.inc"
#include "isajet/q1q2.inc"
#include "isajet/xmssm.inc"
C
      INTEGER JET,K,I,IW,LPRT
      INTEGER NLIST
      CHARACTER*8 WORD,BLANK,LIST(30),LISTW(4),LISTXY(4),LISTSS(85)
      DATA BLANK/'        '/
      DATA LIST/'GL','UP','UB','DN','DB','ST','SB','CH','CB','BT','BB',
     $'TP','TB','NUE','ANUE','E-','E+','NUM','ANUM','MU-','MU+',
     $'NUT','ANUT','TAU-','TAU+','GM','W+','W-','Z0','HIGGS'/
      DATA LISTW/'GM','W+','W-','Z0'/
      DATA LISTXY/'Y','YB','X','XB'/
      DATA LISTSS/'GLSS',
     $'UPSSL','UBSSL','DNSSL','DBSSL','STSSL','SBSSL','CHSSL','CBSSL',
     $'BTSS1','BBSS1','TPSS1','TBSS1',
     $'UPSSR','UBSSR','DNSSR','DBSSR','STSSR','SBSSR','CHSSR','CBSSR',
     $'BTSS2','BBSS2','TPSS2','TBSS2',
     $'W1SS+','W1SS-','W2SS+','W2SS-','Z1SS','Z2SS','Z3SS','Z4SS',
     $'NUEL','ANUEL','EL-','EL+','NUML','ANUML','MUL-','MUL+',
     $'NUTL','ANUTL','TAU1-','TAU1+','ER-','ER+','MUR-','MUR+',
     $'TAU2-','TAU2+',
     $'GL','UP','UB','DN','DB','ST','SB','CH','CB','BT','BB',
     $'TP','TB','NUE','ANUE','E-','E+','NUM','ANUM','MU-','MU+',
     $'NUT','ANUT','TAU-','TAU+','GM','W+','W-','Z0',
     $'HL0','HH0','HA0','H+','H-'/
C
      SETTYP=.FALSE.
C       
      IF(KEYS(5)) GO TO 5
      IF(KEYS(2).AND.GOMSSM) GO TO 5
      IF(KEYS(6).OR.KEYS(9)) GO TO 6
      IF(KEYS(7).AND..NOT.GOMSSM) GO TO 7
      IF(KEYS(7).AND.GOMSSM) GO TO 5
      IF(KEYS(10).AND.GOMSSM) GO TO 5
C
C          JETTYPE flags all processes except WPAIR and HIGGS.
C          NJTTYP is set in READIN to number of non-blank values read.
C          Check for legal jet type names and set appropriate flags.
C
      DO 1000 JET=1,MXGOJ
        IF(NJTTYP(JET).EQ.0) GO TO 1000
C          Initialize everything to .FALSE.
        GOALL(JET)=.FALSE.
        DO 1100 K=1,MXGOQ
          GOQ(K,JET)=.FALSE.
1100    CONTINUE
C          Loop over non-blank JETTYPE entries
        DO 1200 I=1,NJTTYP(JET)
          WORD=JETYP(I,JET)
C          Blank
          IF(WORD.EQ.BLANK) THEN
            GO TO 1200
          ENDIF
C          All
          IF(WORD.EQ.'ALL     ') THEN
            GOALL(JET)=.TRUE.
            DO 1210 K=1,MXGOQ
1210        GOQ(K,JET)=.TRUE.
            GO TO 1000
          ENDIF
C          Quarks
          IF(WORD.EQ.'QUARKS  ') THEN
            DO 1220 K=2,13
1220        GOQ(K,JET)=.TRUE.
            GO TO 1200
          ENDIF
C          Charged leptons
          IF(WORD.EQ.'LEPTONS ') THEN
            DO 1230 K=16,24,4
            GOQ(K,JET)=.TRUE.
1230        GOQ(K+1,JET)=.TRUE.
            GO TO 1200
          ENDIF
C          Neutrinos
          IF(WORD.EQ.'NUS     ') THEN
            DO 1240 K=14,22,4
            GOQ(K,JET)=.TRUE.
1240        GOQ(K+1,JET)=.TRUE.
            GO TO 1200
          ENDIF
C          Explicit types
C          E+E- now also contains W+, W-, Z0
          IF(KEYS(10).OR.KEYS(11).OR.KEYS(12)) THEN
            NLIST=30
          ELSE IF(KEYS(2)) THEN
            NLIST=29
          ELSE
            NLIST=25
          ENDIF
          DO 1250 K=1,NLIST
            IF(WORD.EQ.LIST(K)) THEN
              GOQ(K,JET)=.TRUE.
              GO TO 1200
            ENDIF
1250      CONTINUE
C          Special types for TWOJET
          DO 1270 K=1,4
            IF(KEYS(1).AND.WORD.EQ.LISTXY(K)) THEN
              GOQ(13+K,JET)=.TRUE.
              GO TO 1200
            ENDIF
1270      CONTINUE
C          Special type for PHOTON
          IF(KEYS(8).AND.WORD.EQ.LISTW(1)) THEN
            GOQ(26,JET)=.TRUE.
            GO TO 1200
          ENDIF
C          Error
          WRITE(ITLIS,1300) WORD,JET
1300      FORMAT(1X,A8,' IS NOT RECOGNIZABLE FOR JETTYPE',I1)
          SETTYP=.TRUE.
1200    CONTINUE
1000  CONTINUE
      GO TO 4000
C
C          JETTYPE flags for SUSY
C
5     DO 5000 JET=1,2
        IF(NJTTYP(JET).EQ.0) GO TO 5000
        GOALL(JET)=.FALSE.
        DO 5100 K=1,MXGOQ
5100    GOQ(K,JET)=.FALSE.
        DO 5200 I=1,NJTTYP(JET)
          WORD=JETYP(I,JET)
C          Blank
          IF(WORD.EQ.BLANK) THEN
            GO TO 5200
          ENDIF
C          All
          IF(WORD.EQ.'ALL     ') THEN
            GOALL(JET)=.TRUE.
            DO 5210 K=1,85
5210        GOQ(K,JET)=.TRUE.
            GO TO 5000
          ENDIF
C          Squarks
          IF(WORD.EQ.'SQUARKS ') THEN
            DO 5220 K=2,25
5220        GOQ(K,JET)=.TRUE.
            GO TO 5200
          ENDIF
C           Gauginos
          IF(WORD.EQ.'GAUGINOS') THEN
            DO 5230 K=26,33
5230        GOQ(K,JET)=.TRUE.
            GO TO 5200
          ENDIF
C           Sleptons
          IF(WORD.EQ.'SLEPTONS') THEN
            DO 5240 K=34,51
5240        GOQ(K,JET)=.TRUE.
            GO TO 5200
          ENDIF
C          Explicit susy types
          DO 5300 K=1,85
            IF(WORD.EQ.LISTSS(K)) THEN
              GOQ(K,JET)=.TRUE.
              GO TO 5200
            ENDIF
5300      CONTINUE
5200    CONTINUE
5000  CONTINUE
      GO TO 4000
C
C          JETTYPE and WMODE flags for WPAIR
C          NJTTYP and NWWTYP are the number of non-blank values.
C
6     DO 2000 JET=1,2
        IF(NJTTYP(JET).EQ.0) GO TO 2300
C          Initialize to FALSE
        GOALL(JET)=.FALSE.
        DO 2100 K=1,4
2100    GOQ(K,JET)=.FALSE.
C
C          Loop over non-blank JETTYPE flags
C
        DO 2200 I=1,NJTTYP(JET)
          WORD=JETYP(I,JET)
C          Blank
          IF(WORD.EQ.BLANK) THEN
            GO TO 2200
          ENDIF
C          All
          IF(WORD.EQ.'ALL     ') THEN
            GOALL(JET)=.TRUE.
            DO 2210 K=1,4
2210        GOQ(K,JET)=.TRUE.
            GO TO 2300
          ENDIF
C          Explicit types
          DO 2220 K=1,4
            IF(WORD.EQ.LISTW(K)) THEN
              GOQ(K,JET)=.TRUE.
              GO TO 2200
            ENDIF
2220      CONTINUE
C          Error
          WRITE(ITLIS,1300) WORD,JET
          SETTYP=.TRUE.
2200    CONTINUE
C
C          Loop over nonblank WMODE flags
C
2300    IF(NWWTYP(JET).EQ.0) GO TO 2000
        ALLWW(JET)=.FALSE.
C         Initialize everything to FALSE
        DO 2350 K=1,25
2350    GOWW(K,JET)=.FALSE.
C
        DO 2400 I=1,NWWTYP(JET)
          WORD=WWTYP(I,JET)
          IF(WORD.NE.BLANK) NWWTYP(JET)=I
C          Blank
          IF(WORD.EQ.BLANK) THEN
            GO TO 2400
          ENDIF
C          All
          IF(WORD.EQ.'ALL     ') THEN
            ALLWW(JET)=.TRUE.
            DO 2410 K=1,25
2410        GOWW(K,JET)=.TRUE.
            GO TO 2000
          ENDIF
C          Quarks
          IF(WORD.EQ.'QUARKS  ') THEN
            DO 2420 K=2,13
2420        GOWW(K,JET)=.TRUE.
            GO TO 2400
          ENDIF
C          Charged leptons
          IF(WORD.EQ.'LEPTONS ') THEN
            DO 2430 K=16,24,4
            GOWW(K,JET)=.TRUE.
2430        GOWW(K+1,JET)=.TRUE.
            GO TO 2400
          ENDIF
C          Neutrinos
          IF(WORD.EQ.'NUS     ') THEN
            DO 2440 K=14,22,4
            GOWW(K,JET)=.TRUE.
2440        GOWW(K+1,JET)=.TRUE.
            GO TO 2400
          ENDIF
C          Explicit types
          DO 2450 K=1,25
            IF(WORD.EQ.LIST(K)) THEN
              GOWW(K,JET)=.TRUE.
              GO TO 2400
            ENDIF
2450      CONTINUE
C          Error
          WRITE(ITLIS,2500) WORD,JET
2500      FORMAT(1X,A8,' IS NOT A VALID CODE FOR WMODE',I1)
          SETTYP=.TRUE.
2400    CONTINUE
2000  CONTINUE
      GO TO 4000
C
C          JETTYPE and WMODE flags for HIGGS
C          SUSY HIGGS uses LISTSS order and hence SUSY part
C
7     DO 3000 JET=1,2
        IF(NJTTYP(JET).EQ.0) GO TO 3300
C          Initialize to FALSE
        GOALL(JET)=.FALSE.
        DO 3100 K=1,MXGOQ
3100    GOQ(K,JET)=.FALSE.
C
C          Loop over non-blank JETTYPE flags
C
        DO 3200 I=1,NJTTYP(JET)
          WORD=JETYP(I,JET)
C          Blank
          IF(WORD.EQ.BLANK) THEN
            GO TO 3200
          ENDIF
C          All
          IF(WORD.EQ.'ALL     ') THEN
            GOALL(JET)=.TRUE.
            DO 3210 K=1,MXGOQ
3210        GOQ(K,JET)=.TRUE.
            GO TO 3300
          ENDIF
C          Quarks
          IF(WORD.EQ.'QUARKS  ') THEN
            DO 3220 K=2,13
3220        GOQ(K,JET)=.TRUE.
            GO TO 3200
          ENDIF
C          Charged leptons
          IF(WORD.EQ.'LEPTONS ') THEN
            DO 3240 K=16,24,4
            GOQ(K,JET)=.TRUE.
3240        GOQ(K+1,JET)=.TRUE.
            GO TO 3200
          ENDIF
          DO 3250 K=1,85
            IF(WORD.EQ.LIST(K)) THEN
              GOQ(K,JET)=.TRUE.
              GO TO 3200
            ENDIF
3250      CONTINUE
C          Error
          WRITE(ITLIS,1300) WORD,JET
          SETTYP=.TRUE.
3200    CONTINUE
C
C          Loop over nonblank WMODE flags
C
3300    IF(NWWTYP(JET).EQ.0) GO TO 3000
        ALLWW(JET)=.FALSE.
C         Initialize everything to FALSE
        DO 3350 K=1,25
3350    GOWW(K,JET)=.FALSE.
C
        DO 3400 I=1,NWWTYP(JET)
          WORD=WWTYP(I,JET)
          IF(WORD.NE.BLANK) NWWTYP(JET)=I
C          Blank
          IF(WORD.EQ.BLANK) THEN
            GO TO 3400
          ENDIF
C          All
          IF(WORD.EQ.'ALL     ') THEN
            ALLWW(JET)=.TRUE.
            DO 3410 K=1,25
3410        GOWW(K,JET)=.TRUE.
            GO TO 3000
          ENDIF
C          Quarks
          IF(WORD.EQ.'QUARKS  ') THEN
            DO 3420 K=2,13
3420        GOWW(K,JET)=.TRUE.
            GO TO 3400
          ENDIF
C          Charged leptons
          IF(WORD.EQ.'LEPTONS ') THEN
            DO 3430 K=16,24,4
            GOWW(K,JET)=.TRUE.
3430        GOWW(K+1,JET)=.TRUE.
            GO TO 3400
          ENDIF
C          Neutrinos
          IF(WORD.EQ.'NUS     ') THEN
            DO 3440 K=14,22,4
            GOWW(K,JET)=.TRUE.
3440        GOWW(K+1,JET)=.TRUE.
            GO TO 3400
          ENDIF
C          Explicit types
          DO 3450 K=1,25
            IF(WORD.EQ.LIST(K)) THEN
              GOWW(K,JET)=.TRUE.
              GO TO 3400
            ENDIF
3450      CONTINUE
C          Error
          WRITE(ITLIS,2500) WORD,JET
3500      FORMAT(1X,A8,' IS NOT A VALID CODE FOR WMODE',I1)
          SETTYP=.TRUE.
3400    CONTINUE
3000  CONTINUE
C
C          Set WMODES and ZMODES flags for secondary W+- and Z0
C
4000  DO 4100 IW=1,3
        IF(NWMODE(IW).EQ.0) GO TO 4100
C           Initialize everything to .FALSE.
        DO 4200 K=1,25
4200    GOWMOD(K,IW)=.FALSE.
C          Loop over non-blank WMODE entries
        DO 4300 I=1,NWMODE(IW)
          WORD=WMODES(I,IW)
C          Blank
          IF(WORD.EQ.BLANK) THEN
            GO TO 4300
          ENDIF
C          All
          IF(WORD.EQ.'ALL     ') THEN
            DO 4310 K=1,25
4310        GOWMOD(K,IW)=.TRUE.
            GO TO 4100
          ENDIF
C          Quarks
          IF(WORD.EQ.'QUARKS  ') THEN
            DO 4320 K=2,13
4320        GOWMOD(K,IW)=.TRUE.
            GO TO 4300
          ENDIF
C          Charged leptons
          IF(WORD.EQ.'LEPTONS ') THEN
            DO 4330 K=16,24,4
            GOWMOD(K,IW)=.TRUE.
4330        GOWMOD(K+1,IW)=.TRUE.
            GO TO 4300
          ENDIF
C          Neutrinos
          IF(WORD.EQ.'NUS     ') THEN
            DO 4340 K=14,22,4
            GOWMOD(K,IW)=.TRUE.
4340        GOWMOD(K+1,IW)=.TRUE.
            GO TO 4300
          ENDIF
C          Explicit types
          DO 4350 K=1,25
            IF(WORD.EQ.LIST(K)) THEN
              GOWMOD(K,IW)=.TRUE.
              GO TO 4300
            ENDIF
4350      CONTINUE
C          Error
          WRITE(ITLIS,4380) WORD
4380      FORMAT(1X,A8,' IS NOT RECOGNIZABLE FOR SECONDARY WS')
          SETTYP=.TRUE.
4300    CONTINUE
4100  CONTINUE
C
C          Loop over nonblank WMODE flags FOR WHIGGS
C
      IF (KEYS(10)) THEN
      DO 6000 JET=1,2
6300    IF(NWWTYP(JET).EQ.0) GO TO 6000
        ALLWW(JET)=.FALSE.
C         Initialize everything to FALSE
        DO 6350 K=1,25
6350    GOWW(K,JET)=.FALSE.
C
        DO 6400 I=1,NWWTYP(JET)
          WORD=WWTYP(I,JET)
          IF(WORD.NE.BLANK) NWWTYP(JET)=I
C          Blank
          IF(WORD.EQ.BLANK) THEN
            GO TO 6400
          ENDIF
C          All
          IF(WORD.EQ.'ALL     ') THEN
            ALLWW(JET)=.TRUE.
            DO 6410 K=1,25
6410        GOWW(K,JET)=.TRUE.
            GO TO 6000
          ENDIF
C          Quarks
          IF(WORD.EQ.'QUARKS  ') THEN
            DO 6420 K=2,13
6420        GOWW(K,JET)=.TRUE.
            GO TO 6400
          ENDIF
C          Charged leptons
          IF(WORD.EQ.'LEPTONS ') THEN
            DO 6430 K=16,24,4
            GOWW(K,JET)=.TRUE.
6430        GOWW(K+1,JET)=.TRUE.
            GO TO 6400
          ENDIF
C          Neutrinos
          IF(WORD.EQ.'NUS     ') THEN
            DO 6440 K=14,22,4
            GOWW(K,JET)=.TRUE.
6440        GOWW(K+1,JET)=.TRUE.
            GO TO 6400
          ENDIF
C          Explicit types
          DO 6450 K=1,25
            IF(WORD.EQ.LIST(K)) THEN
              GOWW(K,JET)=.TRUE.
              GO TO 6400
            ENDIF
6450      CONTINUE
C          Error
          WRITE(ITLIS,6500) WORD,JET
6500      FORMAT(1X,A8,' IS NOT A VALID CODE FOR WMODE',I1)
          SETTYP=.TRUE.
6400    CONTINUE
6000  CONTINUE
      END IF
      RETURN
      END
Commit	Line	Data
0795afa3	1	#include "isajet/pilot.h"
	2	LOGICAL FUNCTION SETTYP(LPRT)
	3	C
	4	C Set JETTYPE flags and WMODE flags for WPAIR.
	5	C Set WMODES and ZMODES flags for secondary W+- and Z0.
	6	C Return .FALSE. if no error, .TRUE. otherwise.
	7	C
	8	C Ver 7.18: Initialize all GOQ to false (limit = MXGOQ)
	9	C Use LISTSS for Higgs if GOMSSM
	10	C Ver 7.29: SUSY Higgs decays are done in SETHSS and SIGHSS
	11	C using LISTSS order, so SUSY list should be used.
	12	C I.e., 7.18 fix was wrong.
	13	C
	14	#if defined(CERNLIB_IMPNONE)
	15	IMPLICIT NONE
	16	#endif
	17	#include "isajet/itapes.inc"
	18	#include "isajet/keys.inc"
	19	#include "isajet/types.inc"
	20	#include "isajet/q1q2.inc"
	21	#include "isajet/xmssm.inc"
	22	C
	23	INTEGER JET,K,I,IW,LPRT
	24	INTEGER NLIST
	25	CHARACTER*8 WORD,BLANK,LIST(30),LISTW(4),LISTXY(4),LISTSS(85)
	26	DATA BLANK/' '/
	27	DATA LIST/'GL','UP','UB','DN','DB','ST','SB','CH','CB','BT','BB',
	28	$'TP','TB','NUE','ANUE','E-','E+','NUM','ANUM','MU-','MU+',
	29	$'NUT','ANUT','TAU-','TAU+','GM','W+','W-','Z0','HIGGS'/
	30	DATA LISTW/'GM','W+','W-','Z0'/
	31	DATA LISTXY/'Y','YB','X','XB'/
	32	DATA LISTSS/'GLSS',
	33	$'UPSSL','UBSSL','DNSSL','DBSSL','STSSL','SBSSL','CHSSL','CBSSL',
	34	$'BTSS1','BBSS1','TPSS1','TBSS1',
	35	$'UPSSR','UBSSR','DNSSR','DBSSR','STSSR','SBSSR','CHSSR','CBSSR',
	36	$'BTSS2','BBSS2','TPSS2','TBSS2',
	37	$'W1SS+','W1SS-','W2SS+','W2SS-','Z1SS','Z2SS','Z3SS','Z4SS',
	38	$'NUEL','ANUEL','EL-','EL+','NUML','ANUML','MUL-','MUL+',
	39	$'NUTL','ANUTL','TAU1-','TAU1+','ER-','ER+','MUR-','MUR+',
	40	$'TAU2-','TAU2+',
	41	$'GL','UP','UB','DN','DB','ST','SB','CH','CB','BT','BB',
	42	$'TP','TB','NUE','ANUE','E-','E+','NUM','ANUM','MU-','MU+',
	43	$'NUT','ANUT','TAU-','TAU+','GM','W+','W-','Z0',
	44	$'HL0','HH0','HA0','H+','H-'/
	45	C
	46	SETTYP=.FALSE.
	47	C
	48	IF(KEYS(5)) GO TO 5
	49	IF(KEYS(2).AND.GOMSSM) GO TO 5
	50	IF(KEYS(6).OR.KEYS(9)) GO TO 6
	51	IF(KEYS(7).AND..NOT.GOMSSM) GO TO 7
	52	IF(KEYS(7).AND.GOMSSM) GO TO 5
	53	IF(KEYS(10).AND.GOMSSM) GO TO 5
	54	C
	55	C JETTYPE flags all processes except WPAIR and HIGGS.
	56	C NJTTYP is set in READIN to number of non-blank values read.
	57	C Check for legal jet type names and set appropriate flags.
	58	C
	59	DO 1000 JET=1,MXGOJ
	60	IF(NJTTYP(JET).EQ.0) GO TO 1000
	61	C Initialize everything to .FALSE.
	62	GOALL(JET)=.FALSE.
	63	DO 1100 K=1,MXGOQ
	64	GOQ(K,JET)=.FALSE.
65	1100 CONTINUE
66	C Loop over non-blank JETTYPE entries
67	DO 1200 I=1,NJTTYP(JET)
68	WORD=JETYP(I,JET)
69	C Blank
70	IF(WORD.EQ.BLANK) THEN
71	GO TO 1200
72	ENDIF
73	C All
74	IF(WORD.EQ.'ALL ') THEN
75	GOALL(JET)=.TRUE.
76	DO 1210 K=1,MXGOQ
77	1210 GOQ(K,JET)=.TRUE.
78	GO TO 1000
79	ENDIF
80	C Quarks
81	IF(WORD.EQ.'QUARKS ') THEN
82	DO 1220 K=2,13
83	1220 GOQ(K,JET)=.TRUE.
84	GO TO 1200
85	ENDIF
86	C Charged leptons
87	IF(WORD.EQ.'LEPTONS ') THEN
88	DO 1230 K=16,24,4
89	GOQ(K,JET)=.TRUE.
90	1230 GOQ(K+1,JET)=.TRUE.
91	GO TO 1200
92	ENDIF
93	C Neutrinos
94	IF(WORD.EQ.'NUS ') THEN
95	DO 1240 K=14,22,4
96	GOQ(K,JET)=.TRUE.
97	1240 GOQ(K+1,JET)=.TRUE.
98	GO TO 1200
99	ENDIF
100	C Explicit types
101	C E+E- now also contains W+, W-, Z0
102	IF(KEYS(10).OR.KEYS(11).OR.KEYS(12)) THEN
103	NLIST=30
104	ELSE IF(KEYS(2)) THEN
105	NLIST=29
106	ELSE
107	NLIST=25
108	ENDIF
109	DO 1250 K=1,NLIST
110	IF(WORD.EQ.LIST(K)) THEN
111	GOQ(K,JET)=.TRUE.
112	GO TO 1200
113	ENDIF
114	1250 CONTINUE
115	C Special types for TWOJET
116	DO 1270 K=1,4
117	IF(KEYS(1).AND.WORD.EQ.LISTXY(K)) THEN
118	GOQ(13+K,JET)=.TRUE.
119	GO TO 1200
120	ENDIF
121	1270 CONTINUE
122	C Special type for PHOTON
123	IF(KEYS(8).AND.WORD.EQ.LISTW(1)) THEN
124	GOQ(26,JET)=.TRUE.
125	GO TO 1200
126	ENDIF
127	C Error
128	WRITE(ITLIS,1300) WORD,JET
129	1300 FORMAT(1X,A8,' IS NOT RECOGNIZABLE FOR JETTYPE',I1)
130	SETTYP=.TRUE.
131	1200 CONTINUE
132	1000 CONTINUE
133	GO TO 4000
134	C
135	C JETTYPE flags for SUSY
136	C
137	5 DO 5000 JET=1,2
138	IF(NJTTYP(JET).EQ.0) GO TO 5000
139	GOALL(JET)=.FALSE.
140	DO 5100 K=1,MXGOQ
141	5100 GOQ(K,JET)=.FALSE.
142	DO 5200 I=1,NJTTYP(JET)
143	WORD=JETYP(I,JET)
144	C Blank
145	IF(WORD.EQ.BLANK) THEN
146	GO TO 5200
147	ENDIF
148	C All
149	IF(WORD.EQ.'ALL ') THEN
150	GOALL(JET)=.TRUE.
151	DO 5210 K=1,85
152	5210 GOQ(K,JET)=.TRUE.
153	GO TO 5000
154	ENDIF
155	C Squarks
156	IF(WORD.EQ.'SQUARKS ') THEN
157	DO 5220 K=2,25
158	5220 GOQ(K,JET)=.TRUE.
159	GO TO 5200
160	ENDIF
161	C Gauginos
162	IF(WORD.EQ.'GAUGINOS') THEN
163	DO 5230 K=26,33
164	5230 GOQ(K,JET)=.TRUE.
165	GO TO 5200
166	ENDIF
167	C Sleptons
168	IF(WORD.EQ.'SLEPTONS') THEN
169	DO 5240 K=34,51
170	5240 GOQ(K,JET)=.TRUE.
171	GO TO 5200
172	ENDIF
173	C Explicit susy types
174	DO 5300 K=1,85
175	IF(WORD.EQ.LISTSS(K)) THEN
176	GOQ(K,JET)=.TRUE.
177	GO TO 5200
178	ENDIF
179	5300 CONTINUE
180	5200 CONTINUE
181	5000 CONTINUE
182	GO TO 4000
183	C
184	C JETTYPE and WMODE flags for WPAIR
185	C NJTTYP and NWWTYP are the number of non-blank values.
186	C
187	6 DO 2000 JET=1,2
188	IF(NJTTYP(JET).EQ.0) GO TO 2300
189	C Initialize to FALSE
190	GOALL(JET)=.FALSE.
191	DO 2100 K=1,4
192	2100 GOQ(K,JET)=.FALSE.
193	C
194	C Loop over non-blank JETTYPE flags
195	C
196	DO 2200 I=1,NJTTYP(JET)
197	WORD=JETYP(I,JET)
198	C Blank
199	IF(WORD.EQ.BLANK) THEN
200	GO TO 2200
201	ENDIF
202	C All
203	IF(WORD.EQ.'ALL ') THEN
204	GOALL(JET)=.TRUE.
205	DO 2210 K=1,4
206	2210 GOQ(K,JET)=.TRUE.
207	GO TO 2300
208	ENDIF
209	C Explicit types
210	DO 2220 K=1,4
211	IF(WORD.EQ.LISTW(K)) THEN
212	GOQ(K,JET)=.TRUE.
213	GO TO 2200
214	ENDIF
215	2220 CONTINUE
216	C Error
217	WRITE(ITLIS,1300) WORD,JET
218	SETTYP=.TRUE.
219	2200 CONTINUE
220	C
221	C Loop over nonblank WMODE flags
222	C
223	2300 IF(NWWTYP(JET).EQ.0) GO TO 2000
224	ALLWW(JET)=.FALSE.
225	C Initialize everything to FALSE
226	DO 2350 K=1,25
227	2350 GOWW(K,JET)=.FALSE.
228	C
229	DO 2400 I=1,NWWTYP(JET)
230	WORD=WWTYP(I,JET)
231	IF(WORD.NE.BLANK) NWWTYP(JET)=I
232	C Blank
233	IF(WORD.EQ.BLANK) THEN
234	GO TO 2400
235	ENDIF
236	C All
237	IF(WORD.EQ.'ALL ') THEN
238	ALLWW(JET)=.TRUE.
239	DO 2410 K=1,25
240	2410 GOWW(K,JET)=.TRUE.
241	GO TO 2000
242	ENDIF
243	C Quarks
244	IF(WORD.EQ.'QUARKS ') THEN
245	DO 2420 K=2,13
246	2420 GOWW(K,JET)=.TRUE.
247	GO TO 2400
248	ENDIF
249	C Charged leptons
250	IF(WORD.EQ.'LEPTONS ') THEN
251	DO 2430 K=16,24,4
252	GOWW(K,JET)=.TRUE.
253	2430 GOWW(K+1,JET)=.TRUE.
254	GO TO 2400
255	ENDIF
256	C Neutrinos
257	IF(WORD.EQ.'NUS ') THEN
258	DO 2440 K=14,22,4
259	GOWW(K,JET)=.TRUE.
260	2440 GOWW(K+1,JET)=.TRUE.
261	GO TO 2400
262	ENDIF
263	C Explicit types
264	DO 2450 K=1,25
265	IF(WORD.EQ.LIST(K)) THEN
266	GOWW(K,JET)=.TRUE.
267	GO TO 2400
268	ENDIF
269	2450 CONTINUE
270	C Error
271	WRITE(ITLIS,2500) WORD,JET
272	2500 FORMAT(1X,A8,' IS NOT A VALID CODE FOR WMODE',I1)
273	SETTYP=.TRUE.
274	2400 CONTINUE
275	2000 CONTINUE
276	GO TO 4000
277	C
278	C JETTYPE and WMODE flags for HIGGS
279	C SUSY HIGGS uses LISTSS order and hence SUSY part
280	C
281	7 DO 3000 JET=1,2
282	IF(NJTTYP(JET).EQ.0) GO TO 3300
283	C Initialize to FALSE
284	GOALL(JET)=.FALSE.
285	DO 3100 K=1,MXGOQ
286	3100 GOQ(K,JET)=.FALSE.
287	C
288	C Loop over non-blank JETTYPE flags
289	C
290	DO 3200 I=1,NJTTYP(JET)
291	WORD=JETYP(I,JET)
292	C Blank
293	IF(WORD.EQ.BLANK) THEN
294	GO TO 3200
295	ENDIF
296	C All
297	IF(WORD.EQ.'ALL ') THEN
298	GOALL(JET)=.TRUE.
299	DO 3210 K=1,MXGOQ
300	3210 GOQ(K,JET)=.TRUE.
301	GO TO 3300
302	ENDIF
303	C Quarks
304	IF(WORD.EQ.'QUARKS ') THEN
305	DO 3220 K=2,13
306	3220 GOQ(K,JET)=.TRUE.
307	GO TO 3200
308	ENDIF
309	C Charged leptons
310	IF(WORD.EQ.'LEPTONS ') THEN
311	DO 3240 K=16,24,4
312	GOQ(K,JET)=.TRUE.
313	3240 GOQ(K+1,JET)=.TRUE.
314	GO TO 3200
315	ENDIF
316	DO 3250 K=1,85
317	IF(WORD.EQ.LIST(K)) THEN
318	GOQ(K,JET)=.TRUE.
319	GO TO 3200
320	ENDIF
321	3250 CONTINUE
322	C Error
323	WRITE(ITLIS,1300) WORD,JET
324	SETTYP=.TRUE.
325	3200 CONTINUE
326	C
327	C Loop over nonblank WMODE flags
328	C
329	3300 IF(NWWTYP(JET).EQ.0) GO TO 3000
330	ALLWW(JET)=.FALSE.
331	C Initialize everything to FALSE
332	DO 3350 K=1,25
333	3350 GOWW(K,JET)=.FALSE.
334	C
335	DO 3400 I=1,NWWTYP(JET)
336	WORD=WWTYP(I,JET)
337	IF(WORD.NE.BLANK) NWWTYP(JET)=I
338	C Blank
339	IF(WORD.EQ.BLANK) THEN
340	GO TO 3400
341	ENDIF
342	C All
343	IF(WORD.EQ.'ALL ') THEN
344	ALLWW(JET)=.TRUE.
345	DO 3410 K=1,25
346	3410 GOWW(K,JET)=.TRUE.
347	GO TO 3000
348	ENDIF
349	C Quarks
350	IF(WORD.EQ.'QUARKS ') THEN
351	DO 3420 K=2,13
352	3420 GOWW(K,JET)=.TRUE.
353	GO TO 3400
354	ENDIF
355	C Charged leptons
356	IF(WORD.EQ.'LEPTONS ') THEN
357	DO 3430 K=16,24,4
358	GOWW(K,JET)=.TRUE.
359	3430 GOWW(K+1,JET)=.TRUE.
360	GO TO 3400
361	ENDIF
362	C Neutrinos
363	IF(WORD.EQ.'NUS ') THEN
364	DO 3440 K=14,22,4
365	GOWW(K,JET)=.TRUE.
366	3440 GOWW(K+1,JET)=.TRUE.
367	GO TO 3400
368	ENDIF
369	C Explicit types
370	DO 3450 K=1,25
371	IF(WORD.EQ.LIST(K)) THEN
372	GOWW(K,JET)=.TRUE.
373	GO TO 3400
374	ENDIF
375	3450 CONTINUE
376	C Error
377	WRITE(ITLIS,2500) WORD,JET
378	3500 FORMAT(1X,A8,' IS NOT A VALID CODE FOR WMODE',I1)
379	SETTYP=.TRUE.
380	3400 CONTINUE
381	3000 CONTINUE
382	C
383	C Set WMODES and ZMODES flags for secondary W+- and Z0
384	C
385	4000 DO 4100 IW=1,3
386	IF(NWMODE(IW).EQ.0) GO TO 4100
387	C Initialize everything to .FALSE.
388	DO 4200 K=1,25
389	4200 GOWMOD(K,IW)=.FALSE.
390	C Loop over non-blank WMODE entries
391	DO 4300 I=1,NWMODE(IW)
392	WORD=WMODES(I,IW)
393	C Blank
394	IF(WORD.EQ.BLANK) THEN
395	GO TO 4300
396	ENDIF
397	C All
398	IF(WORD.EQ.'ALL ') THEN
399	DO 4310 K=1,25
400	4310 GOWMOD(K,IW)=.TRUE.
401	GO TO 4100
402	ENDIF
403	C Quarks
404	IF(WORD.EQ.'QUARKS ') THEN
405	DO 4320 K=2,13
406	4320 GOWMOD(K,IW)=.TRUE.
407	GO TO 4300
408	ENDIF
409	C Charged leptons
410	IF(WORD.EQ.'LEPTONS ') THEN
411	DO 4330 K=16,24,4
412	GOWMOD(K,IW)=.TRUE.
413	4330 GOWMOD(K+1,IW)=.TRUE.
414	GO TO 4300
415	ENDIF
416	C Neutrinos
417	IF(WORD.EQ.'NUS ') THEN
418	DO 4340 K=14,22,4
419	GOWMOD(K,IW)=.TRUE.
420	4340 GOWMOD(K+1,IW)=.TRUE.
421	GO TO 4300
422	ENDIF
423	C Explicit types
424	DO 4350 K=1,25
425	IF(WORD.EQ.LIST(K)) THEN
426	GOWMOD(K,IW)=.TRUE.
427	GO TO 4300
428	ENDIF
429	4350 CONTINUE
430	C Error
431	WRITE(ITLIS,4380) WORD
432	4380 FORMAT(1X,A8,' IS NOT RECOGNIZABLE FOR SECONDARY WS')
433	SETTYP=.TRUE.
434	4300 CONTINUE
435	4100 CONTINUE
436	C
437	C Loop over nonblank WMODE flags FOR WHIGGS
438	C
439	IF (KEYS(10)) THEN
440	DO 6000 JET=1,2
441	6300 IF(NWWTYP(JET).EQ.0) GO TO 6000
442	ALLWW(JET)=.FALSE.
443	C Initialize everything to FALSE
444	DO 6350 K=1,25
445	6350 GOWW(K,JET)=.FALSE.
446	C
447	DO 6400 I=1,NWWTYP(JET)
448	WORD=WWTYP(I,JET)
449	IF(WORD.NE.BLANK) NWWTYP(JET)=I
450	C Blank
451	IF(WORD.EQ.BLANK) THEN
452	GO TO 6400
453	ENDIF
454	C All
455	IF(WORD.EQ.'ALL ') THEN
456	ALLWW(JET)=.TRUE.
457	DO 6410 K=1,25
458	6410 GOWW(K,JET)=.TRUE.
459	GO TO 6000
460	ENDIF
461	C Quarks
462	IF(WORD.EQ.'QUARKS ') THEN
463	DO 6420 K=2,13
464	6420 GOWW(K,JET)=.TRUE.
465	GO TO 6400
466	ENDIF
467	C Charged leptons
468	IF(WORD.EQ.'LEPTONS ') THEN
469	DO 6430 K=16,24,4
470	GOWW(K,JET)=.TRUE.
471	6430 GOWW(K+1,JET)=.TRUE.
472	GO TO 6400
473	ENDIF
474	C Neutrinos
475	IF(WORD.EQ.'NUS ') THEN
476	DO 6440 K=14,22,4
477	GOWW(K,JET)=.TRUE.
478	6440 GOWW(K+1,JET)=.TRUE.
479	GO TO 6400
480	ENDIF
481	C Explicit types
482	DO 6450 K=1,25
483	IF(WORD.EQ.LIST(K)) THEN
484	GOWW(K,JET)=.TRUE.
485	GO TO 6400
486	ENDIF
487	6450 CONTINUE
488	C Error
489	WRITE(ITLIS,6500) WORD,JET
490	6500 FORMAT(1X,A8,' IS NOT A VALID CODE FOR WMODE',I1)
491	SETTYP=.TRUE.
492	6400 CONTINUE
493	6000 CONTINUE
494	END IF
495	RETURN
496	END