1 #include "isajet/pilot.h"
4 C Set the MSSM Higgs parameters in /HCON/.
5 C HMASS = Higgs mass for HTYPE
7 C HGAMSS = Higgs partial widths. Note HGAMSS is not
8 C necessarily diagonal for SUSY decays.
9 C ZSTARS = minimum allowed mass for Z*
11 C Note LISTSS(78) => W+, LISTSS(79) => W-, LISTSS(80) => Z0
13 #if defined(CERNLIB_IMPNONE)
16 #include "isajet/itapes.inc"
17 #include "isajet/hcon.inc"
18 #include "isajet/listss.inc"
19 #include "isajet/q1q2.inc"
20 #include "isajet/ssmode.inc"
21 #include "isajet/sstype.inc"
22 #include "isajet/wcon.inc"
26 INTEGER I,J,N,IQ1,IQ2,IW,K
27 INTEGER LISTJ(25),LISTW(4)
29 DATA LISTJ/9,1,-1,2,-2,3,-3,4,-4,5,-5,6,-6,
30 $11,-11,12,-12,13,-13,14,-14,15,-15,16,-16/
31 DATA LISTW/10,80,-80,90/
36 WRITE(ITLIS,*) ' YOU MUST SELECT AN HTYPE FOR SUSY HIGGS'
37 WRITE(ITLIS,*) ' JOB TERMINATED'
48 C Extract widths from SSMODE common block
49 C Note the only 3-body modes are Zff or Wff
50 C These are added to the ZZ and WW entries in HCONSS,
51 C and the Z* or W* decay is generated later, as for SM Higgs
54 IF(ISSMOD(N).NE.IHTYPE) GO TO 200
56 IF(JSSMOD(3,N).NE.0) THEN
58 IF(IABS(JSSMOD(1,N)).EQ.80) THEN
59 HGAMSS(78,79)=HGAMSS(78,79)+0.5*GSSMOD(N)
60 HGAMSS(79,78)=HGAMSS(79,78)+0.5*GSSMOD(N)
61 ELSEIF(JSSMOD(1,N).EQ.90) THEN
62 HGAMSS(80,80)=HGAMSS(80,80)+GSSMOD(N)
64 WRITE(ITLIS,1000) ISSMOD(N),(JSSMOD(K,N),K=1,5)
65 1000 FORMAT(' SETHSS: UNEXPECTED MODE ',I8,' --> ',5I8)
72 IF(JSSMOD(1,N).NE.LISTSS(I)) GO TO 210
74 IF(JSSMOD(2,N).NE.LISTSS(J)) GO TO 220
75 HGAMSS(I,J)=HGAMSS(I,J)+.5*GSSMOD(N)
76 HGAMSS(J,I)=HGAMSS(J,I)+.5*GSSMOD(N)
81 WRITE(ITLIS,1000) ISSMOD(N),(JSSMOD(K,N),K=1,5)
85 C W* and Z* mass limits
90 ZSTARS(IW,I)=AMASS(LISTW(IW))
93 IF(IQ2.EQ.0) GO TO 320
94 IF(GOWW(IQ1,I).AND.GOWW(IQ2,I)) THEN
95 AM12=AMASS(LISTJ(IQ1))+AMASS(LISTJ(IQ2))+1.0
96 ZSTARS(IW,I)=MIN(ZSTARS(IW,I),AM12)