1 #include "isajet/pilot.h"
4 C Read in user data and execute SETTYP if appropriate values
5 C are set. IFL return values:
6 C IFL = 0 Good parameter set
8 C IFL > 0 Error. Program will continue reading data but
9 C will exit when END or unrecognizable keyword
12 #if defined(CERNLIB_IMPNONE)
15 #include "isajet/itapes.inc"
16 #include "isajet/mbgen.inc"
17 #include "isajet/force.inc"
18 #include "isajet/dkytab.inc"
19 #include "isajet/qcdpar.inc"
20 #include "isajet/eepar.inc"
21 #include "isajet/idrun.inc"
22 #include "isajet/frgpar.inc"
23 #include "isajet/keys.inc"
24 #include "isajet/kkgrav.inc"
25 #include "isajet/prtout.inc"
26 #include "isajet/seed.inc"
27 #include "isajet/types.inc"
28 #include "isajet/primar.inc"
29 #include "isajet/jetlim.inc"
30 #include "isajet/nodcay.inc"
31 #include "isajet/wcon.inc"
32 #include "isajet/dylim.inc"
33 #include "isajet/qlmass.inc"
34 #include "isajet/q1q2.inc"
35 #include "isajet/jetpar.inc"
36 #include "isajet/isloop.inc"
37 #include "isajet/tcpar.inc"
38 #include "isajet/xmssm.inc"
39 #include "isajet/sugnu.inc"
40 #if defined(CERNLIB_PDFLIB)
41 #include "isajet/w50510.inc"
44 #if defined(CERNLIB_PDFLIB)
45 #include "isajet/w50517.inc"
48 #include "isajet/hcon.inc"
49 #include "isajet/mglims.inc"
52 CHARACTER*8 TTL(10),WORD,LSTRUC,BLANK
53 CHARACTER*8 WTYP(4),RDID(2)
57 INTEGER IFL,I1,I2,I3,J1,I,IKEY,IJ,J,KK,IDABS
59 INTEGER NSEL,K,KFORCE(5),INDEX,IDG1,IDG2,IDG3,IDG4,IDXLEP
64 #if defined(CERNLIB_SINGLE)
67 #if defined(CERNLIB_DOUBLE)
70 #if defined(CERNLIB_PDFLIB)
71 CHARACTER*20 PDFPAR(20)
73 #if (defined(CERNLIB_PDFLIB))&&(defined(CERNLIB_SINGLE))
75 REAL DX,DSCALE,DXPDF(-6:6)
77 #if (defined(CERNLIB_PDFLIB))&&(defined(CERNLIB_DOUBLE))
78 DOUBLE PRECISION PDFVAL(20)
79 DOUBLE PRECISION DX,DSCALE,DXPDF(-6:6)
82 C Overlapping variable flags.
83 DATA NLAP/1,2,3, 1,2,7 ,1,2,8, 1,3,5, 1,3,6, 1,3,7, 1,3,8, 1,5,7,
84 X 1,5,8, 1,6,7, 1,6,8, 2,3,7, 2,3,8, 3,5,7, 3,6,7, 3,5,8,
92 10 FORMAT('1',//5X,'***** ',A40,' *****')
94 11 FORMAT(////30X,' COMMANDS READ BY READIN')
102 IF(TTL(1).EQ.'STOP ') THEN
107 C Read energy and no. of events
109 READ(ITCOM,*) ECM,NEVENT,NEVPRT,NJUMP
110 WRITE(ITLIS,*) ECM,NEVENT,NEVPRT,NJUMP
112 C Reset all variables and set process if title is not 'SAME'
114 IF(TTL(1).NE.'SAME ') THEN
128 IF(REAC.EQ.'TWOJET ') THEN
132 ELSEIF(REAC.EQ.'E+E- ') THEN
136 ELSEIF(REAC.EQ.'DRELLYAN') THEN
140 ELSEIF(REAC.EQ.'MINBIAS ') THEN
144 ELSEIF(REAC.EQ.'SUPERSYM'.OR.REAC.EQ.'SUSY ') THEN
148 ELSEIF(REAC.EQ.'WPAIR ') THEN
152 ELSEIF(REAC.EQ.'HIGGS ') THEN
156 ELSEIF(REAC.EQ.'PHOTON ') THEN
160 ELSEIF(REAC.EQ.'TCOLOR ') THEN
164 ELSEIF(REAC.EQ.'WHIGGS ') THEN
168 ELSEIF(REAC.EQ.'EXTRADIM') THEN
172 ELSEIF(REAC.EQ.'ZJJ ') THEN
178 890 WRITE(ITLIS,1999)
188 C Read keyword. For each recognized keyword read corresponding
189 C variables and set LOC flag.
198 IF(WORD.EQ.'END ') THEN
199 C Check for previous error
201 C Check inconsistent limits
202 IF(LOC(2)*LOC(5).NE.0.OR.LOC(2)*LOC(6).NE.0) THEN
206 C Set and check jet types
207 IF(LOC(15).NE.0.OR.LOC(37).NE.0.OR.LOC(46).NE.0) THEN
213 C Check MSSM/SUGRA conflict
214 IF((LOC(51).NE.0.OR.LOC(52).NE.0.OR.LOC(53).NE.0).AND.
219 C Check overlapping limits
224 IF(LOC(I1)*LOC(I2)*LOC(I3).NE.0) WRITE(ITLIS,1001)
227 #if defined(CERNLIB_PDFLIB)
228 IF(ISTRUC.EQ.-999) THEN
231 $ '1********************************'/
233 $ ' * INITIALIZE PDFLIB FOR ISAJET *'/
235 $ ' ********************************'/)
238 CALL PDFSET(PDFPAR,PDFVAL)
239 CALL PFTOPDG(0.5D0,1.0D2,DXPDF)
243 C Check EXTRADIM parameters are set
244 IF(KEYS(11).AND.LOC(72).EQ.0) THEN
245 WRITE(ITLIS,*) 'YOU FORGOT TO SET EXTRAD PARAMETERS'
253 IF(WORD.EQ.'P ') THEN
254 READ(ITCOM,*) (PMIN(K),PMAX(K),K=1,NJET)
255 WRITE(ITLIS,*) (PMIN(K),PMAX(K),K=1,NJET)
261 IF(WORD.EQ.'Y ') THEN
262 READ(ITCOM,*) (YJMIN(K),YJMAX(K),K=1,NJET)
263 WRITE(ITLIS,*) (YJMIN(K),YJMAX(K),K=1,NJET)
269 IF(WORD.EQ.'X ') THEN
270 READ(ITCOM,*) (XJMIN(K),XJMAX(K),K=1,NJET)
271 WRITE(ITLIS,*) (XJMIN(K),XJMAX(K),K=1,NJET)
277 IF(WORD.EQ.'PHI ') THEN
278 READ(ITCOM,*) (PHIMIN(K),PHIMAX(K),K=1,NJET)
279 WRITE(ITLIS,*) (PHIMIN(K),PHIMAX(K),K=1,NJET)
285 IF(WORD.EQ.'TH '.OR.WORD.EQ.'THETA ') THEN
286 READ(ITCOM,*) (THMIN(K),THMAX(K),K=1,NJET)
287 WRITE(ITLIS,*) (THMIN(K),THMAX(K),K=1,NJET)
294 IF(WORD.EQ.'PT '.OR.WORD.EQ.'PPERP ') THEN
295 READ(ITCOM,*) (PTMIN(K),PTMAX(K),K=1,NJET)
296 WRITE(ITLIS,*) (PTMIN(K),PTMAX(K),K=1,NJET)
303 IF(WORD.EQ.'NODECAY ') THEN
304 READ(ITCOM,571) NODCAY
306 WRITE(ITLIS,572) NODCAY
313 IF(WORD.EQ.'NOETA ') THEN
314 READ(ITCOM,571) NOETA
315 WRITE(ITLIS,572) NOETA
321 IF(WORD.EQ.'NOPI0 ') THEN
322 READ(ITCOM,571) NOPI0
323 WRITE(ITLIS,572) NOPI0
329 IF(WORD.EQ.'BEAMS ') THEN
330 READ(ITCOM,*) RDID(1),RDID(2)
331 WRITE(ITLIS,*) RDID(1),RDID(2)
335 IF(RDID(K).EQ.'P ') IDIN(K)=+1120
336 IF(RDID(K).EQ.'AP ') IDIN(K)=-1120
337 IF(RDID(K).EQ.'N ') IDIN(K)=+1220
338 IF(RDID(K).EQ.'AN ') IDIN(K)=-1220
340 IF(IDIN(1)*IDIN(2).EQ.0) THEN
349 IF(WORD.EQ.'FRAGMENT') THEN
357 IF(WORD.EQ.'SEED ') THEN
361 WRITE(XSEED,'(E24.15)') SEED
366 C Keywords JETTYPE1, JETTYPE2, JETTYPE3, ...
367 C (Yes, this is ugly)
368 IF(WORD.EQ.'JETTYPE1'.OR.WORD.EQ.'JETTYPE2'.OR.
369 $WORD.EQ.'JETTYPE3'.OR.WORD.EQ.'JETTYPE4'.OR.
370 $WORD.EQ.'JETTYPE5'.OR.WORD.EQ.'JETTYPE6'.OR.
371 $WORD.EQ.'JETTYPE7'.OR.WORD.EQ.'JETTYPE8') THEN
372 IF(WORD.EQ.'JETTYPE1') IJ=1
373 IF(WORD.EQ.'JETTYPE2') IJ=2
374 IF(WORD.EQ.'JETTYPE3') IJ=3
375 IF(WORD.EQ.'JETTYPE4') IJ=4
376 IF(WORD.EQ.'JETTYPE5') IJ=5
377 IF(WORD.EQ.'JETTYPE6') IJ=6
378 IF(WORD.EQ.'JETTYPE7') IJ=7
379 IF(WORD.EQ.'JETTYPE8') IJ=8
383 READ(ITCOM,*) (JETYP(K,IJ),K=1,30)
385 152 IF(JETYP(K,IJ).NE.BLANK) NJTTYP(IJ)=NJTTYP(IJ)+1
386 WRITE(ITLIS,*) (JETYP(K,IJ),K=1,NJTTYP(IJ))
392 IF(WORD.EQ.'SIN2W ') THEN
400 IF(WORD.EQ.'TMASS ') THEN
401 READ(ITCOM,*) AMLEP(6),AMLEP(7),AMLEP(8)
402 WRITE(ITLIS,*) AMLEP(6),AMLEP(7),AMLEP(8)
408 IF(WORD.EQ.'QMH ') THEN
409 READ(ITCOM,*) QMIN,QMAX
410 WRITE(ITLIS,*) QMIN,QMAX
416 IF(WORD.EQ.'QMW ') THEN
417 READ(ITCOM,*) QMIN,QMAX
418 WRITE(ITLIS,*) QMIN,QMAX
424 IF(WORD.EQ.'QTW ') THEN
425 READ(ITCOM,*) QTMIN,QTMAX
426 WRITE(ITLIS,*) QTMIN,QTMAX
432 IF(WORD.EQ.'YW ') THEN
433 READ(ITCOM,*) YWMIN,YWMAX
434 WRITE(ITLIS,*) YWMIN,YWMAX
440 IF(WORD.EQ.'XW ') THEN
441 READ(ITCOM,*) XWMIN,XWMAX
442 WRITE(ITLIS,*) XWMIN,XWMAX
448 IF(WORD.EQ.'THW ') THEN
449 READ(ITCOM,*) THWMIN,THWMAX
450 WRITE(ITLIS,*) THWMIN,THWMAX
456 IF(WORD.EQ.'PHIW ') THEN
457 READ(ITCOM,*) PHWMIN,PHWMAX
458 WRITE(ITLIS,*) PHWMIN,PHWMAX
464 IF(WORD.EQ.'NONUNU ') THEN
465 READ(ITCOM,571) NONUNU
466 WRITE(ITLIS,572) NONUNU
472 IF(WORD.EQ.'WTYPE ') THEN
480 IF(WTYP(K).EQ.'GM ') GODY(1)=.TRUE.
481 IF(WTYP(K).EQ.'W+ ') GODY(2)=.TRUE.
482 IF(WTYP(K).EQ.'W- ') GODY(3)=.TRUE.
483 IF(WTYP(K).EQ.'Z0 ') GODY(4)=.TRUE.
486 IF(GODY(2).OR.GODY(3)) JWTYP=3
488 IF((GODY(2).OR.GODY(3)).AND.(GODY(1).OR.GODY(4))) THEN
497 IF(WORD.EQ.'LAMBDA ') THEN
506 IF(WORD.EQ.'NTRIES ') THEN
508 WRITE(ITLIS,*) NTRIES
514 IF(WORD.EQ.'CUTOFF ') THEN
515 READ(ITCOM,*) CUTOFF,CUTPOW
516 WRITE(ITLIS,*) CUTOFF,CUTPOW
522 IF(WORD.EQ.'XGEN ') THEN
530 IF(WORD.EQ.'SIGQT ') THEN
538 IF(WORD.EQ.'CUTJET ') THEN
540 WRITE(ITLIS,*) CUTJET
546 IF(WORD.EQ.'WFUDGE ') THEN
548 WRITE(ITLIS,*) WFUDGE
554 IF(WORD.EQ.'STRUC ') THEN
557 WRITE(ITLIS,4) LSTRUC
558 IF(LSTRUC.EQ.'OWENS ') ISTRUC=1
559 IF(LSTRUC.EQ.'BAIER ') ISTRUC=2
560 IF(LSTRUC.EQ.'EICHTEN '.OR.LSTRUC.EQ.'EHLQ ') ISTRUC=3
561 IF(LSTRUC.EQ.'DUKE '.OR.LSTRUC.EQ.'DO ') ISTRUC=4
562 IF(LSTRUC.EQ.'CTEQ2L ') ISTRUC=5
563 IF(LSTRUC.EQ.'CTEQ '.OR.LSTRUC.EQ.'CTEQ3L ') ISTRUC=6
573 IF(WORD.EQ.'NPOMERON') THEN
574 READ(ITCOM,*) MNPOM,MXPOM
575 WRITE(ITLIS,*) MNPOM,MXPOM
576 IF(MNPOM.LT.1.OR.MNPOM.GT.MXPOM.OR.MXPOM.GT.LIMPOM) THEN
585 IF(WORD.EQ.'FORCE ') THEN
587 IF(NFORCE.GT.MXFORC-1) THEN
593 READ(ITCOM,*) IFORCE(NFORCE),(KFORCE(K),K=1,5)
594 CALL ORDER(IFORCE(NFORCE),KFORCE,MFORCE(1,NFORCE),
596 WRITE(ITLIS,*) IFORCE(NFORCE),(MFORCE(K,NFORCE),K=1,5)
607 352 MFORCE(K,NFORCE+1)=IDANTI(MFORCE(K,NFORCE))
615 IF(WORD.EQ.'FORCE1 ') THEN
617 IF(NFORCE.GT.MXFORC) THEN
623 READ(ITCOM,*) IFORCE(NFORCE),(KFORCE(K),K=1,5)
624 CALL ORDER(IFORCE(NFORCE),KFORCE,MFORCE(1,NFORCE),
626 WRITE(ITLIS,*) IFORCE(NFORCE),(MFORCE(K,NFORCE),K=1,5)
627 IF(IABS(IFORCE(NFORCE)).LT.6) THEN
635 C Keyword HMASSES - also see HMASS
636 IF(WORD.EQ.'HMASSES ') THEN
637 CALL FLAVOR(80,I1,I2,I3,J1,INDEX)
638 READ(ITCOM,*) (AMLEP(INDEX+K),K=1,9)
639 WRITE(ITLIS,*) (AMLEP(INDEX+K),K=1,9)
644 C Keywords WMODE1,WMODE2
645 IF(WORD.EQ.'WMODE1 '.OR.WORD.EQ.'WMODE2 ') THEN
646 IF(WORD.EQ.'WMODE1 ') IJ=1
647 IF(WORD.EQ.'WMODE2 ') IJ=2
648 READ(ITCOM,*) (WWTYP(K,IJ),K=1,25)
650 372 IF(WWTYP(K,IJ).NE.BLANK) NWWTYP(IJ)=NWWTYP(IJ)+1
651 WRITE(ITLIS,*) (WWTYP(K,IJ),K=1,NWWTYP(IJ))
657 IF(WORD.EQ.'NOEVOLVE') THEN
658 READ (ITCOM,571) NOEVOL
659 WRITE(ITLIS,572) NOEVOL
665 IF(WORD.EQ.'NOHADRON') THEN
666 READ (ITCOM,571) NOHADR
667 WRITE(ITLIS,572) NOHADR
673 IF(WORD.EQ.'GAUGINO ') THEN
674 CALL FLAVOR(29,I1,I2,I3,J1,IDG1)
675 CALL FLAVOR(30,I1,I2,I3,J1,IDG2)
676 CALL FLAVOR(39,I1,I2,I3,J1,IDG3)
677 CALL FLAVOR(40,I1,I2,I3,J1,IDG4)
678 READ(ITCOM,*) AMLEP(IDG1),AMLEP(IDG2),AMLEP(IDG3),AMLEP(IDG4)
679 WRITE(ITLIS,*) AMLEP(IDG1),AMLEP(IDG2),AMLEP(IDG3),AMLEP(IDG4)
685 IF(WORD.EQ.'SQUARK ') THEN
686 CALL FLAVOR(21,I1,I2,I3,J1,IDXQKL)
687 READ(ITCOM,*) (AMLEP(IDXQKL+K-1),K=1,6)
688 WRITE(ITLIS,*) (AMLEP(IDXQKL+K-1),K=1,6)
689 CALL FLAVOR(41,I1,I2,I3,J1,IDXQKR)
691 AMLEP(IDXQKR+K-1)=AMLEP(IDXQKL+K-1)
698 IF(WORD.EQ.'SLEPTON ') THEN
699 CALL FLAVOR(31,I1,I2,I3,J1,IDXLEP)
700 READ(ITCOM,*) (AMLEP(IDXLEP+K-1),K=1,6)
701 WRITE(ITLIS,*) (AMLEP(IDXLEP+K-1),K=1,6)
707 IF(WORD.EQ.'NSIGMA ') THEN
709 WRITE(ITLIS,*) NSIGMA
715 IF(WORD.EQ.'XGENSS ') THEN
716 READ(ITCOM,*) XGENSS(9),(XGENSS(KK),KK=1,8)
717 WRITE(ITLIS,*) XGENSS(9),(XGENSS(KK),KK=1,8)
722 C Keyword HMASS - just standard Higgs
723 IF(WORD.EQ.'HMASS ') THEN
724 CALL FLAVOR(81,I1,I2,I3,J1,INDEX)
725 READ(ITCOM,*) AMLEP(INDEX)
726 WRITE(ITLIS,*) AMLEP(INDEX)
731 C Keywords WPMODE, WMMODE, Z0MODE
732 IF(WORD.EQ.'WPMODE '.OR.WORD.EQ.'WMMODE '
733 $.OR.WORD.EQ.'Z0MODE ') THEN
734 IF(WORD.EQ.'WPMODE ') IJ=1
735 IF(WORD.EQ.'WMMODE ') IJ=2
736 IF(WORD.EQ.'Z0MODE ') IJ=3
737 READ(ITCOM,*) (WMODES(K,IJ),K=1,25)
739 463 IF(WMODES(K,IJ).NE.BLANK) NWMODE(IJ)=NWMODE(IJ)+1
740 WRITE(ITLIS,*) (WMODES(K,IJ),K=1,NWMODE(IJ))
746 IF(WORD.EQ.'WMASS ') THEN
747 READ(ITCOM,*) AMW,AMZ
748 WRITE(ITLIS,*) AMW,AMZ
753 CALL FLAVOR(80,I1,I2,I3,J,INDEX)
755 CALL FLAVOR(90,I1,I2,I3,J,INDEX)
762 IF(WORD.EQ.'NEVOLVE ') THEN
764 WRITE(ITLIS,*) NEVOLV
770 IF(WORD.EQ.'NHADRON ') THEN
772 WRITE(ITLIS,*) NFRGMN
778 IF(WORD.EQ.'TCMASS ') THEN
779 READ(ITCOM,*) TCMRHO,TCGRHO
780 WRITE(ITLIS,*) TCMRHO,TCGRHO
785 C Keyword MSSMA: gluino, mu, mha, tanb
786 IF(WORD.EQ.'MSSMA ') THEN
787 READ(ITCOM,*) XGLSS,XMUSS,XHASS,XTBSS
788 WRITE(ITLIS,*) XGLSS,XMUSS,XHASS,XTBSS
794 C Keyword MSSMB: 1st generation soft terms
795 IF(WORD.EQ.'MSSMB ') THEN
796 READ(ITCOM,*) XQ1SS,XDRSS,XURSS,XL1SS,XERSS
797 WRITE(ITLIS,*) XQ1SS,XDRSS,XURSS,XL1SS,XERSS
803 C Keyword MSSMC: 3rd generation soft terms
804 IF(WORD.EQ.'MSSMC ') THEN
805 READ(ITCOM,*)XQ3SS,XBRSS,XTRSS,XL3SS,XTARSS,XATSS,XABSS,XATASS
806 WRITE(ITLIS,*)XQ3SS,XBRSS,XTRSS,XL3SS,XTARSS,XATSS,XABSS,XATASS
812 C Keyword PDFLIB: parameters for PDFLIB
813 #if defined(CERNLIB_PDFLIB)
814 IF(WORD.EQ.'PDFLIB ') THEN
819 READ(ITCOM,*) (PDFPAR(I),PDFVAL(I),I=1,20)
821 IF(PDFPAR(I).NE.' ') THEN
822 WRITE(ITLIS,*) PDFPAR(I),PDFVAL(I)
832 IF(WORD.EQ.'SUGRA ') THEN
833 READ(ITCOM,*) XM0SU,XMHSU,XA0SU,XTGBSU,XSMUSU
834 WRITE(ITLIS,*) XM0SU,XMHSU,XA0SU,XTGBSU,XSMUSU
842 IF(WORD.EQ.'HTYPE ') THEN
847 IF(HTYPE.EQ.'HL0 ') IHTYPE=82
848 IF(HTYPE.EQ.'HH0 ') IHTYPE=83
849 IF(HTYPE.EQ.'HA0 ') IHTYPE=84
851 WRITE(ITLIS,2000) HTYPE
858 IF(WORD.EQ.'EPOL ') THEN
859 READ(ITCOM,*) PLEM,PLEP
860 WRITE(ITLIS,*) PLEM,PLEP
865 C Keyword MSSMD: optional 2nd geenration soft terms
866 IF(WORD.EQ.'MSSMD ') THEN
867 READ(ITCOM,*) XQ2SS,XSRSS,XCRSS,XL2SS,XMRSS
868 WRITE(ITLIS,*) XQ2SS,XSRSS,XCRSS,XL2SS,XMRSS
873 C Keyword MSSME: optional U(1) and SU(2) gaugino masses
874 IF(WORD.EQ.'MSSME ') THEN
875 READ(ITCOM,*) XM1SS,XM2SS
876 WRITE(ITLIS,*) XM1SS,XM2SS
881 C Keyword GMSB: gauge-mediated SUSY breaking model
882 IF(WORD.EQ.'GMSB ') THEN
883 READ(ITCOM,*) XLAMGM,XMESGM,XN5GM,XTGBSU,XSMUSU,XCMGV
884 WRITE(ITLIS,*) XLAMGM,XMESGM,XN5GM,XTGBSU,XSMUSU,XCMGV
891 C Keyword NUSUG1: optional GUT scale gaugino masses
892 IF(WORD.EQ.'NUSUG1 ') THEN
893 READ(ITCOM,*) XNUSUG(1),XNUSUG(2),XNUSUG(3)
894 WRITE(ITLIS,*) XNUSUG(1),XNUSUG(2),XNUSUG(3)
899 C Keyword NUSUG2: optional GUT scale A terms
900 IF(WORD.EQ.'NUSUG2 ') THEN
901 READ(ITCOM,*) XNUSUG(6),XNUSUG(5),XNUSUG(4)
902 WRITE(ITLIS,*) XNUSUG(6),XNUSUG(5),XNUSUG(4)
907 C Keyword NUSUG3: optional GUT scale Higgs masses
908 IF(WORD.EQ.'NUSUG3 ') THEN
909 READ(ITCOM,*) XNUSUG(7),XNUSUG(8)
910 WRITE(ITLIS,*) XNUSUG(7),XNUSUG(8)
915 C Keyword NUSUG4: optional GUT scale 1st/2nd gen. masses
916 IF(WORD.EQ.'NUSUG4 ') THEN
917 READ(ITCOM,*) XNUSUG(13),XNUSUG(11),XNUSUG(12),XNUSUG(10)
919 WRITE(ITLIS,*) XNUSUG(13),XNUSUG(11),XNUSUG(12),XNUSUG(10)
925 C Keyword NUSUG5: optional GUT scale 3rd gen. masses
926 IF(WORD.EQ.'NUSUG5 ') THEN
927 READ(ITCOM,*) XNUSUG(18),XNUSUG(16),XNUSUG(17),XNUSUG(15)
929 WRITE(ITLIS,*) XNUSUG(18),XNUSUG(16),XNUSUG(17),XNUSUG(15)
935 C Keyword NOGRAV: No gravitino decays
936 IF(WORD.EQ.'NOGRAV ') THEN
937 READ(ITCOM,571) NOGRAV
938 WRITE(ITLIS,572) NOGRAV
943 C Keyword MGVTNO: Sets the gravitino mass
944 IF(WORD.EQ.'MGVTNO ') THEN
946 WRITE(ITLIS,*) XMGVTO
951 C Keyword AL3UNI: Impose alpha_s unification at M_GUT
952 IF(WORD.EQ.'AL3UNI ') THEN
953 READ(ITCOM,571) AL3UNI
954 WRITE(ITLIS,572) AL3UNI
959 C Keyword GMSB2: additional GMSB parameters
960 IF(WORD.EQ.'GMSB2 ') THEN
961 READ(ITCOM,*) XRSLGM,XDHDGM,XDHUGM,XDYGM,XN51GM,XN52GM,XN53GM
962 WRITE(ITLIS,*) XRSLGM,XDHDGM,XDHUGM,XDYGM,XN51GM,XN52GM,XN53GM
967 C Keyword EEBREM: invoke bremsstrahlung in e+e- reactions
968 IF(WORD.EQ.'EEBREM ') THEN
969 READ(ITCOM,*) RSHMIN,RSHMAX
970 WRITE(ITLIS,*) RSHMIN,RSHMAX
976 C Keyword EEBEAM: invoke beamstrahlung in e+e- reactions
977 IF(WORD.EQ.'EEBEAM ') THEN
978 READ(ITCOM,*) RSHMIN,RSHMAX,UPSLON,SIGZ
979 WRITE(ITLIS,*) RSHMIN,RSHMAX,UPSLON,SIGZ
986 C Keyword QMKKG (QMW for EXTRADIM)
987 IF(WORD.EQ.'QMKKG ') THEN
988 READ(ITCOM,*) QMIN,QMAX
989 WRITE(ITLIS,*) QMIN,QMAX
994 C Keyword QTKKG (QTW for EXTRADIM)
995 IF(WORD.EQ.'QTKKG ') THEN
996 READ(ITCOM,*) QTMIN,QTMAX
997 WRITE(ITLIS,*) QTMIN,QTMAX
1002 C Keyword EXTRAD for EXTRADIM
1003 IF(WORD.EQ.'EXTRAD ') THEN
1004 READ(ITCOM,*) NEXTRAD,MASSD,UVCUT
1005 WRITE(ITLIS,*) NEXTRAD,MASSD,UVCUT
1011 IF(WORD.EQ.'MIJLIM ') THEN
1012 READ(ITCOM,*) JLIM1,JLIM2,AMLIM1,AMLIM2
1013 WRITE(ITLIS,*) JLIM1,JLIM2,AMLIM1,AMLIM2
1014 IF(JLIM1.EQ.0.AND.JLIM2.EQ.0) THEN
1021 ELSEIF(JLIM1.GT.0.AND.JLIM1.LE.NJET.AND.JLIM2.GT.0.AND.
1022 $ JLIM2.LE.NJET) THEN
1023 AMIJMN(JLIM1,JLIM2)=AMLIM1
1024 AMIJMN(JLIM2,JLIM1)=AMLIM1
1025 AMIJMX(JLIM1,JLIM2)=AMLIM2
1026 AMIJMX(JLIM2,JLIM1)=AMLIM2
1036 IF(WORD.EQ.'MTOT ') THEN
1037 READ(ITCOM,*) EHMGMN,EHMGMX
1038 WRITE(ITLIS,*) EHMGMN,EHMGMX
1044 IF(WORD.EQ.'SUGRHN ') THEN
1045 READ(ITCOM,*) XMN3NR,XMAJNR,XANSS,XNRSS
1046 WRITE(ITLIS,*) XMN3NR,XMAJNR,XANSS,XNRSS
1052 IF(WORD.EQ.'AMSB ') THEN
1053 READ(ITCOM,*) XM0SU,XMHSU,XTGBSU,XSMUSU
1054 WRITE(ITLIS,*) XM0SU,XMHSU,XTGBSU,XSMUSU
1063 IF(WORD.EQ.'SSBCSC ') THEN
1065 WRITE(ITLIS,*) XSBCS
1072 WRITE(ITLIS,2000) WORD
1076 C Error message or warnings
1078 1001 FORMAT(//2X,'YOU HAVE GIVEN LIMITS FOR AN OVERLAPPING SET',
1079 $ ' OF VARIABLES. SET MINIMIZING PPERP INTERVAL WILL BE USED.')
1080 1999 FORMAT(//' YOU FORGOT TO SELECT A PROCESS FOR GENERATION.'
1081 $ /' AVAILABLE AT PRESENT ARE ',
1082 $ /' TWOJET E+E- DRELLYAN MINBIAS WPAIR SUPERSYM,'
1083 $ /' HIGGS PHOTON TCOLOR')
1084 2000 FORMAT(//2X,A8,' IS NOT A RECOGNIZABLE PARAMETER. JOB TERMINATED')
1085 2001 FORMAT(//2X,' YOU CANNOT GIVE LIMITS FOR BOTH THETA AND Y.',
1086 $ ' MAKE UP YOUR MIND. JOB TERMINATED.')
1087 2002 FORMAT(/' WHAT IS THAT SUPPOSED TO BE')
1088 2003 FORMAT(/' YOU CANNOT RUN WS AND Z0 OR GAMMAS AT THE',
1089 $ ' SAME TIME. JOB TERMINATED')
1090 2004 FORMAT(//' PARAMETER OUT OF RANGE. JOB TERMINATED.')
1091 2005 FORMAT(//' YOU CANNOT FORCE DECAY OF A QUARK, YOU MUST CHOSE '
1093 2006 FORMAT(//2X,' INVALID JETTYPE VALUES. JOB TERMINATED.')
1094 2007 FORMAT(//2X,'YOU CANNOT USE MSSM AND SUGRA SIMULTANEOUSLY')
1095 2008 FORMAT(//2X,'INVALID JET-JET MASS LIMITS. JOB TERMINATED.')