1 #include "isajet/pilot.h"
4 C Set up user data with interactive commands from standard
7 C V6.26: Fixed FORCE to allow automatic top decays and to use
8 C correct limit for maximum number.
10 #if defined(CERNLIB_IMPNONE)
13 #include "isajet/itapes.inc"
14 #include "isajet/mbgen.inc"
15 #include "isajet/force.inc"
16 #include "isajet/eepar.inc"
17 #include "isajet/dkytab.inc"
18 #include "isajet/qcdpar.inc"
19 #include "isajet/idrun.inc"
20 #include "isajet/frgpar.inc"
21 #include "isajet/keys.inc"
22 #include "isajet/prtout.inc"
23 #include "isajet/seed.inc"
24 #include "isajet/types.inc"
25 #include "isajet/primar.inc"
26 #include "isajet/jetlim.inc"
27 #include "isajet/nodcay.inc"
28 #include "isajet/wcon.inc"
29 #include "isajet/dylim.inc"
30 #include "isajet/qlmass.inc"
31 #include "isajet/q1q2.inc"
32 #include "isajet/jetpar.inc"
33 #include "isajet/isloop.inc"
34 #include "isajet/tcpar.inc"
35 #include "isajet/xmssm.inc"
36 #include "isajet/sugnu.inc"
37 #if defined(CERNLIB_PDFLIB)
38 #include "isajet/w50510.inc"
41 #if defined(CERNLIB_PDFLIB)
42 #include "isajet/w50517.inc"
45 #include "isajet/hcon.inc"
46 #include "isajet/kkgrav.inc"
47 #include "isajet/mglims.inc"
49 #include "isajet/const.inc"
52 CHARACTER*8 TTL(10),WORD,LSTRUC,BLANK,LIST(29)
53 CHARACTER*8 WTYP(4),RDID(2)
55 INTEGER IFL,I1,I2,I3,J1,I,IKEY,IJ,J,KK,IDXQK
56 INTEGER NSEL,K,KFORCE(5),INDEX,IDG1,IDG2,IDG3,IDG4,IDXQ,IDXLEP
57 #if defined(CERNLIB_SINGLE)
60 #if defined(CERNLIB_DOUBLE)
64 REAL EMIN,YLIM1,YLIM2,ANGL1,ANGL2,XLIM1,XLIM2,AMW,AMZ
67 #if defined(CERNLIB_PDFLIB)
68 CHARACTER*20 PDFPAR(20)
70 #if (defined(CERNLIB_PDFLIB))&&(defined(CERNLIB_SINGLE))
72 REAL DX,DSCALE,DXPDF(-6:6)
74 #if (defined(CERNLIB_PDFLIB))&&(defined(CERNLIB_DOUBLE))
75 DOUBLE PRECISION PDFVAL(20)
76 DOUBLE PRECISION DX,DSCALE,DXPDF(-6:6)
81 C Overlapping variable flags.
82 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,
83 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,
86 C Duplicate data to avoid requiring ALDATA
87 DATA LIST/'GL','UP','UB','DN','DB','ST','SB','CH','CB','BT','BB',
88 $'TP','TB','NUE','ANUE','E-','E+','NUM','ANUM','MU-','MU+',
89 $'NUT','ANUT','TAU-','TAU+','ALL','QUARKS','LEPTONS','NUS'/
97 11 FORMAT(//,10X,' Use CAPITAL LETTERS only for all commands.',//,
98 1 10X,' Give title for this run.')
105 IF(TTL(1).EQ.'STOP ') GOTO 999
107 C Read energy and no. of events
109 PRINT*,' total energy, # of events, # to print, # to skip.'
110 READ*,ECM,NEVENT,NEVPRT,NJUMP
111 WRITE(ITCOM,*)ECM,NEVENT,NEVPRT,NJUMP
113 C Reset all variables if title is not 'SAME'
115 IF(TTL(1).NE.'SAME ') THEN
121 C Read reaction and set keys and NJET
124 22 FORMAT(/,' Choose one of the following reactions:',/,
125 $' TWOJET E+E- DRELLYAN MINBIAS SUSY WPAIR HIGGS',
126 $' PHOTON TCOLOR WHIGGS')
133 IF(REAC.EQ.'TWOJET ') THEN
137 ELSEIF(REAC.EQ.'E+E- ') THEN
141 ELSEIF(REAC.EQ.'DRELLYAN') THEN
145 ELSEIF(REAC.EQ.'MINBIAS ') THEN
149 ELSEIF(REAC.EQ.'SUPERSYM'.OR.REAC.EQ.'SUSY ') THEN
153 ELSEIF(REAC.EQ.'WPAIR ') THEN
157 ELSEIF(REAC.EQ.'HIGGS ') THEN
161 ELSEIF(REAC.EQ.'PHOTON ') THEN
165 ELSEIF(REAC.EQ.'TCOLOR ') THEN
169 ELSEIF(REAC.EQ.'WHIGGS ') THEN
174 PRINT* ,' You must choose a valid reaction, try again.'
192 C Read keyword. For each recognized keyword read corresponding
193 C variables and set LOC flag.
198 $' Choose among the following:',/,
200 $' JETTYPE1 JETTYPE2 JETTYPE3 P PHI PT TH X Y WMODE1 WMODE2',/,
201 $' Drell Yan/Higgs limits:',/,
202 $' HTYPE PHIW QMH QMW QTW THW WTYPE XW YW',/,
204 $' FORCE FORCE1 NODECAY NOETA NOEVOLVE NOFRGMNT NOPI0',/,
205 $' Physics parameters:',/,
206 $' CUTJET CUTOFF FRAGMENT GAUGINO HMASS HMASSES LAMBDA',/,
207 $' MSSMA MSSMB MSSMC MSSMD MSSME SIGQT SIN2W SLEPTON SQUARK',/,
208 $' SUGRA TCMASS TMASS WMASS XGEN',/,
210 $' BEAMS EPOL NPOMERON NSIGMA NTRIES PDFLIB SEED STRUC WFUDGE',/,
211 $' WMMODE WPMODE Z0MODE',/,
212 $' Terminate with END. Use HELP for help.',/)
214 100 PRINT* ,' Give a variable name.'
216 IF(WORD.EQ.'HELP ') GO TO 1000
220 IF(WORD.EQ.'END ') THEN
226 889 FORMAT(//,10X,' NEW RUN',/,' If you are finished answer STOP.',/,
227 1 ' If you answer SAME the parameters from previous run will be',
228 2 ' used,',/,' unless you request otherwise.')
233 IF(WORD.EQ.'P ') THEN
237 READ*, (PMIN(K),PMAX(K),K=1,NJET)
238 WRITE(ITCOM,*)(PMIN(K),PMAX(K),K=1,NJET)
244 IF(WORD.EQ.'Y ') THEN
247 PRINT 563,YLIM1,YLIM2
248 READ*, (YJMIN(K),YJMAX(K),K=1,NJET)
249 WRITE(ITCOM,*)(YJMIN(K),YJMAX(K),K=1,NJET)
255 IF(WORD.EQ.'X ') THEN
258 PRINT 563,XLIM1,XLIM2
259 READ*, (XJMIN(K),XJMAX(K),K=1,NJET)
260 WRITE(ITCOM,*)(XJMIN(K),XJMAX(K),K=1,NJET)
266 IF(WORD.EQ.'PHI ') THEN
269 PRINT 563,ANGL1,ANGL2
270 READ*, (PHIMIN(K),PHIMAX(K),K=1,NJET)
271 WRITE(ITCOM,*)(PHIMIN(K),PHIMAX(K),K=1,NJET)
277 IF(WORD.EQ.'TH '.OR.WORD.EQ.'THETA ') THEN
281 READ*, (THMIN(K),THMAX(K),K=1,NJET)
282 WRITE(ITCOM,*)(THMIN(K),THMAX(K),K=1,NJET)
289 IF(WORD.EQ.'PT '.OR.WORD.EQ.'PPERP ') THEN
293 READ*, (PTMIN(K),PTMAX(K),K=1,NJET)
294 WRITE(ITCOM,*)(PTMIN(K),PTMAX(K),K=1,NJET)
301 IF(WORD.EQ.'NODECAY ') THEN
306 WRITE(ITCOM,571) NODCAY
312 IF(WORD.EQ.'NOETA ') THEN
316 WRITE(ITCOM,571) NOETA
322 IF(WORD.EQ.'NOPI0 ') THEN
326 WRITE(ITCOM,571) NOPI0
332 IF(WORD.EQ.'BEAMS ') THEN
334 58 PRINT*,' Select each beam enclosed in single quotes.'
335 PRINT*,' Allowed names are: P N AP AN'
336 READ *,RDID(1),RDID(2)
340 IF(RDID(K).EQ.'P ') IDIN(K)=+1120
341 IF(RDID(K).EQ.'AP ') IDIN(K)=-1120
342 IF(RDID(K).EQ.'N ') IDIN(K)=+1220
343 IF(RDID(K).EQ.'AN ') IDIN(K)=-1220
345 IF(IDIN(1)*IDIN(2).EQ.0) GOTO 58
347 WRITE(ITCOM,4) RDID(1),RDID(2)
352 IF(WORD.EQ.'FRAGMENT') THEN
362 IF(WORD.EQ.'SEED ') THEN
372 C Keywords JETTYPE1, JETTYPE2, JETTYPE3
373 IF(WORD.EQ.'JETTYPE1'.OR.WORD.EQ.'JETTYPE2'.OR.WORD.EQ.
376 C Print allowed types
380 6001 FORMAT(' Give a list of jet types enclosed in single quotes ',
381 $ 'and separated by commas.',/,' Terminate list with a /')
382 6002 FORMAT(' The recognized W types are: W+, W-, Z0, GM and ALL')
386 6003 FORMAT(' The recognized jet types are:',/,4(1X,8A8,/),1X,5A8)
390 6004 FORMAT(' A fourth generation is also allowed:',/,
391 $ ' X, XB for up, Y, YB for down.',/,
392 $ ' The fourth generation is normally turned off.')
394 IF(WORD.EQ.'JETTYPE1') IJ=1
395 IF(WORD.EQ.'JETTYPE2') IJ=2
396 IF(WORD.EQ.'JETTYPE3') IJ=3
397 READ*,(JETYP(K,IJ),K=1,25)
399 61 IF(JETYP(K,IJ).NE.BLANK) NJTTYP(IJ)=NJTTYP(IJ)+1
400 C Check that only legal jet types are in the list
401 IF(SETTYP(0)) GOTO 60
402 WRITE(ITCOM,4)(JETYP(K,IJ),K=1,NJTTYP(IJ))
404 IF(NJTTYP(IJ).LT.25) WRITE(ITCOM,*) '/'
410 IF(WORD.EQ.'SIN2W ') THEN
420 IF(WORD.EQ.'TMASS ') THEN
423 READ*, AMLEP(6),AMLEP(7),AMLEP(8)
424 WRITE(ITCOM,*) AMLEP(6),AMLEP(7),AMLEP(8)
429 C Keyword QMH (note that it uses same variable as QMW)
430 IF(WORD.EQ.'QMH ') THEN
432 PRINT*,' Give minimum and maximum MASS for Higgs.'
434 WRITE(ITCOM,*)QMIN,QMAX
440 IF(WORD.EQ.'QMW ') THEN
442 PRINT*,' Give minimum and maximum MASS for W(Z0).'
444 WRITE(ITCOM,*)QMIN,QMAX
450 IF(WORD.EQ.'QTW ') THEN
452 PRINT*,' Give minimum and maximum PT for W(Z0).'
454 PRINT*,' Fix QTW to zero for standard Drell-Yan.'
456 WRITE(ITCOM,*) QTMIN,QTMAX
462 IF(WORD.EQ.'YW ') THEN
464 PRINT*,' Give minimum and maximum Y for W(Z0).'
465 PRINT 563,YLIM1,YLIM2
467 WRITE(ITCOM,*)YWMIN,YWMAX
473 IF(WORD.EQ.'XW ') THEN
475 PRINT*,' Give minimum and maximum X for W(Z0).'
476 PRINT 563,XLIM1,XLIM2
478 WRITE(ITCOM,*)XWMIN,XWMAX
484 IF(WORD.EQ.'THW ') THEN
486 PRINT*,' Give minimum and maximum THETA for W(Z0).'
489 WRITE(ITCOM,*)THWMIN,THWMAX
495 IF(WORD.EQ.'PHIW ') THEN
497 PRINT*,' Give minimum and maximum PHI for W(Z0).'
498 PRINT 563,ANGL1,ANGL2
500 WRITE(ITCOM,*)PHWMIN,PHWMAX
506 IF(WORD.EQ.'NONUNU ') THEN
510 WRITE(ITCOM,571) NONUNU
516 IF(WORD.EQ.'WTYPE ') THEN
518 92 PRINT*,' Choose from one set of bosons: W+,W- or Z0,GM'
519 PRINT*,' Enclose each in single quotes and terminate with / '
525 IF(WTYP(K).EQ.'GM ') GODY(1)=.TRUE.
526 IF(WTYP(K).EQ.'W+ ') GODY(2)=.TRUE.
527 IF(WTYP(K).EQ.'W- ') GODY(3)=.TRUE.
528 IF(WTYP(K).EQ.'Z0 ') GODY(4)=.TRUE.
531 IF(GODY(2).OR.GODY(3)) JWTYP=3
533 IF((GODY(3).OR.GODY(2)).AND.(GODY(1).OR.GODY(4)).OR.JWTYP
545 IF(WORD.EQ.'LAMBDA ') THEN
556 IF(WORD.EQ.'NTRIES ') THEN
566 IF(WORD.EQ.'CUTOFF ') THEN
570 WRITE(ITCOM,*)CUTOFF,CUTPOW
576 IF(WORD.EQ.'XGEN ') THEN
586 IF(WORD.EQ.'SIGQT ') THEN
596 IF(WORD.EQ.'CUTJET ') THEN
600 WRITE(ITCOM,*) CUTJET
606 IF(WORD.EQ.'WFUDGE ') THEN
610 WRITE(ITCOM,*) WFUDGE
616 IF(WORD.EQ.'STRUC ') THEN
619 69 PRINT*,' Choose a structure function: OWENS, BAIER or EICHTEN'
620 PRINT*,' Enclose in single quotes'
622 IF(LSTRUC.EQ.'OWENS ') ISTRUC=1
623 IF(LSTRUC.EQ.'BAIER ') ISTRUC=2
624 IF(LSTRUC.EQ.'EICHTEN ') ISTRUC=3
625 IF(LSTRUC.EQ.'DUKE '.OR.LSTRUC.EQ.'DO ') ISTRUC=4
626 IF(LSTRUC.EQ.'CTEQ2L ') ISTRUC=5
627 IF(LSTRUC.EQ.'CTEQ '.OR.LSTRUC.EQ.'CTEQ3L ') ISTRUC=6
628 IF(ISTRUC.EQ.0) GOTO 69
629 WRITE(ITCOM,*) LSTRUC
635 IF(WORD.EQ.'NPOMERON') THEN
637 93 PRINT*,' Give min and max; allowed range is 1 to 20.'
639 IF(MNPOM.LT.1.OR.MNPOM.GT.MXPOM.OR.MXPOM.GT.LIMPOM)
641 WRITE(ITCOM,*)MNPOM,MXPOM
647 IF(WORD.EQ.'FORCE ') THEN
648 PRINT*,' Force decay of a particle and its anti-particle.'
649 PRINT*,' Use FORCE1 if you want to decay only one of them.'
650 PRINT*,' Give a particle IDENT and IDENTs for decay mode.'
651 PRINT*,' Max. number of decay products is 5.'
653 IF(NFTEMP.LE.MXFORC) THEN
656 72 MFORCE(K,NFTEMP)=0
657 READ*,IFORCE(NFTEMP),(MFORCE(K,NFTEMP),K=1,5)
658 WRITE(ITCOM,*)IFORCE(NFTEMP),(MFORCE(K,NFTEMP),K=1,5)
661 PRINT*,' You exceeded the maximum allowed forced decays.'
667 IF(WORD.EQ.'FORCE1 ') THEN
668 PRINT*,' Force decay of a particle.'
669 PRINT*,' Give a particle IDENT and IDENTs for decay mode.'
670 PRINT*,' Max. number of decay products is 5.'
672 IF(NFTEMP.LE.MXFORC) THEN
675 73 MFORCE(K,NFTEMP)=0
676 READ*,IFORCE(NFTEMP),(MFORCE(K,NFTEMP),K=1,5)
677 WRITE(ITCOM,*)IFORCE(NFTEMP),(MFORCE(K,NFTEMP),K=1,5)
680 PRINT*,' You exceeded the maximun allowed forced decays.'
685 C Keyword HMASSES - also see HMASS
686 IF(WORD.EQ.'HMASSES ') THEN
689 $ ' Give Higgs masses (HIGGS,H10,H20,H30,H40,H1+,H2+,H1++,H2++)'
690 CALL FLAVOR(80,I1,I2,I3,J1,INDEX)
691 READ*,(AMLEP(INDEX+K),K=1,9)
692 WRITE(ITCOM,*)(AMLEP(INDEX+K),K=1,9)
697 C Keywords WMODE1,WMODE2
698 IF(WORD.EQ.'WMODE1 '.OR.WORD.EQ.'WMODE2 ') THEN
702 IF(WORD.EQ.'WMODE1 ') IJ=1
703 IF(WORD.EQ.'WMODE2 ') IJ=2
704 READ*,(WWTYP(K,IJ),K=1,25)
706 372 IF(WWTYP(K,IJ).NE.BLANK) NWWTYP(IJ)=NWWTYP(IJ)+1
707 IF(SETTYP(0)) GOTO 95
708 WRITE(ITCOM,4)(WWTYP(K,IJ),K=1,NWWTYP(IJ))
709 IF(NWWTYP(IJ).LT.25) WRITE(ITCOM,*) '/'
715 IF(WORD.EQ.'NOEVOLVE') THEN
719 WRITE(ITCOM,571) NOEVOL
725 IF(WORD.EQ.'NOHADRON') THEN
729 WRITE(ITCOM,571) NOEVOL
735 IF(WORD.EQ.'GAUGINO ') THEN
737 PRINT*,' Give masses (GLSS,GMSS,W+SS,W-SS)'
738 CALL FLAVOR(29,I1,I2,I3,J1,IDG1)
739 CALL FLAVOR(30,I1,I2,I3,J1,IDG2)
740 CALL FLAVOR(39,I1,I2,I3,J1,IDG3)
741 CALL FLAVOR(40,I1,I2,I3,J1,IDG4)
742 READ*, AMLEP(IDG1),AMLEP(IDG2),AMLEP(IDG3),AMLEP(IDG4)
743 WRITE(ITCOM,*) AMLEP(IDG1),AMLEP(IDG2),AMLEP(IDG3),AMLEP(IDG4)
749 IF(WORD.EQ.'SQUARK ') THEN
751 PRINT*,' Give masses (UPSS,DNSS,STSS,CHSS,BTSS,TPSS).'
752 CALL FLAVOR(21,I1,I2,I3,J1,IDXQK)
753 READ*, (AMLEP(IDXQK+K-1),K=1,6)
754 WRITE(ITCOM,*) (AMLEP(IDXQK+K-1),K=1,6)
760 IF(WORD.EQ.'SLEPTON ') THEN
762 PRINT*,' Give masses (NUESS,E-SS,NUMSS,MU-SS,NUTSS,T-SS).'
763 CALL FLAVOR(31,I1,I2,I3,J1,IDXLEP)
764 READ*, (AMLEP(IDXLEP+K-1),K=1,6)
765 WRITE(ITCOM,*) (AMLEP(IDXLEP+K-1),K=1,6)
771 IF(WORD.EQ.'NSIGMA ') THEN
775 WRITE(ITCOM,*) NSIGMA
781 IF(WORD.EQ.'XGENSS ') THEN
784 READ*, XGENSS(9),(XGENSS(KK),KK=1,8)
785 WRITE(ITCOM,*) XGENSS(9),(XGENSS(KK),KK=1,8)
790 C Keyword HMASS - just standard Higgs
791 IF(WORD.EQ.'HMASS ') THEN
794 CALL FLAVOR(81,I1,I2,I3,J1,INDEX)
796 WRITE(ITCOM,*) AMLEP(INDEX)
801 C Keywords WPMODE, WMMODE, Z0MODE
802 IF(WORD.EQ.'WPMODE '.OR.WORD.EQ.'WMMODE '
803 $.OR.WORD.EQ.'Z0MODE ') THEN
804 IF(WORD.EQ.'WPMODE ') IJ=1
805 IF(WORD.EQ.'WMMODE ') IJ=2
806 IF(WORD.EQ.'Z0MODE ') IJ=3
810 READ *, (WMODES(K,IJ),K=1,25)
812 463 IF(WMODES(K,IJ).NE.BLANK) NWMODE(IJ)=NWMODE(IJ)+1
813 WRITE(ITCOM,*) (WMODES(K,IJ),K=1,NWMODE(IJ))
819 IF(WORD.EQ.'WMASS ') THEN
822 WRITE(ITCOM,*) AMW,AMZ
827 CALL FLAVOR(80,I1,I2,I3,J,INDEX)
829 CALL FLAVOR(90,I1,I2,I3,J,INDEX)
836 IF(WORD.EQ.'NEVOLVE ') THEN
840 WRITE(ITCOM,*) NEVOLV
846 IF(WORD.EQ.'NHADRON ') THEN
856 IF(WORD.EQ.'TCMASS ') THEN
860 WRITE(ITCOM,*) TCMRHO,TCGRHO
866 IF(WORD.EQ.'MSSMA ') THEN
869 READ *, XGLSS,XMUSS,XHASS,XTBSS
870 WRITE(ITCOM,*) XGLSS,XMUSS,XHASS,XTBSS
877 IF(WORD.EQ.'MSSMB ') THEN
880 READ*, XQ1SS,XDRSS,XURSS,XL1SS,XERSS
881 WRITE(ITCOM,*) XQ1SS,XDRSS,XURSS,XL1SS,XERSS
888 IF(WORD.EQ.'MSSMC ') THEN
891 READ*, XQ3SS,XBRSS,XTRSS,XL3SS,XTARSS,XATSS,XABSS,XATASS
892 WRITE(ITCOM,*)XQ3SS,XBRSS,XTRSS,XL3SS,XTARSS,XATSS,XABSS,XATASS
898 C Keyword PDFLIB: parameters for PDFLIB
899 #if defined(CERNLIB_PDFLIB)
900 IF(WORD.EQ.'PDFLIB ') THEN
907 READ*, (PDFPAR(I),PDFVAL(I),I=1,20)
909 IF(PDFPAR(I).NE.' ') THEN
910 WRITE(ITCOM,5401) PDFPAR(I)
911 5401 FORMAT(1H',A,1H')
912 WRITE(ITCOM,*) PDFVAL(I)
923 IF(WORD.EQ.'SUGRA ') THEN
926 READ*, XM0SU,XMHSU,XA0SU,XTGBSU,XSMUSU
927 WRITE(ITCOM,*) XM0SU,XMHSU,XA0SU,XTGBSU,XSMUSU
935 IF(WORD.EQ.'HTYPE ') THEN
936 PRINT*,' Enter Higgs type (HL0, HH0, HA0) in single quotes'
939 IF(HTYPE.EQ.'HL0 ') IHTYPE=82
940 IF(HTYPE.EQ.'HH0 ') IHTYPE=83
941 IF(HTYPE.EQ.'HA0 ') IHTYPE=84
943 PRINT*, 'Invalid Higgs type - try again'
953 IF(WORD.EQ.'EPOL ') THEN
957 WRITE(ITCOM,*) PLEM,PLEP
963 IF(WORD.EQ.'MSSMD ') THEN
966 READ*,XQ2SS,XSRSS,XCRSS,XL2SS,XMRSS
967 WRITE(ITCOM,*) XQ2SS,XSRSS,XCRSS,XL2SS,XMRSS
973 IF(WORD.EQ.'MSSME ') THEN
977 WRITE(ITCOM,*) XM1SS,XM2SS
983 IF(WORD.EQ.'GMSB ') THEN
986 READ*, XLAMGM,XMESGM,XN5GM,XTGBSU,XSMUSU,XCMGV
987 WRITE(ITCOM,*) XLAMGM,XMESGM,XN5GM,XTGBSU,XSMUSU,XCMGV
994 C Keyword NUSUG1: optional GUT scale gaugino masses
995 IF(WORD.EQ.'NUSUG1 ') THEN
998 READ*, XNUSUG(1),XNUSUG(2),XNUSUG(3)
999 WRITE(ITCOM,*) XNUSUG(1),XNUSUG(2),XNUSUG(3)
1004 C Keyword NUSUG2: optional GUT scale A terms
1005 IF(WORD.EQ.'NUSUG2 ') THEN
1008 READ*, XNUSUG(6),XNUSUG(5),XNUSUG(4)
1009 WRITE(ITCOM,*) XNUSUG(6),XNUSUG(5),XNUSUG(4)
1014 C Keyword NUSUG3: optional GUT scale Higgs masses
1015 IF(WORD.EQ.'NUSUG3 ') THEN
1018 READ*, XNUSUG(7),XNUSUG(8)
1019 WRITE(ITCOM,*) XNUSUG(7),XNUSUG(8)
1024 C Keyword NUSUG4: optional GUT scale 1st/2nd gen. masses
1025 IF(WORD.EQ.'NUSUG4 ') THEN
1028 READ*, XNUSUG(13),XNUSUG(11),XNUSUG(12),XNUSUG(10)
1030 WRITE(ITCOM,*) XNUSUG(13),XNUSUG(11),XNUSUG(12),XNUSUG(10)
1036 C Keyword NUSUG5: optional GUT scale 3rd gen. masses
1037 IF(WORD.EQ.'NUSUG5 ') THEN
1040 READ*, XNUSUG(18),XNUSUG(16),XNUSUG(17),XNUSUG(15)
1042 WRITE(ITCOM,*) XNUSUG(18),XNUSUG(16),XNUSUG(17),XNUSUG(15)
1048 C Keyword NOGRAV: No gravitino decays
1049 IF(WORD.EQ.'NOGRAV ') THEN
1053 WRITE(ITCOM,571) NOGRAV
1058 C Keyword MGVTNO: Sets the gravitino mass
1059 IF(WORD.EQ.'MGVTNO ') THEN
1063 WRITE(ITCOM,*) XMGVTO
1068 C Keyword AL3UNI: Impose alpha_s unification at M_GUT
1069 IF(WORD.EQ.'AL3UNI ') THEN
1073 WRITE(ITCOM,*) AL3UNI
1078 C Keyword GMSB2: additional GMSB parameters
1079 IF(WORD.EQ.'GMSB2 ') THEN
1082 READ*, XRSLGM,XDHDGM,XDHUGM,XDYGM,XN51GM,XN52GM,XN53GM
1083 WRITE(ITCOM,*) XRSLGM,XDHDGM,XDHUGM,XDYGM,XN51GM,XN52GM,XN53GM
1088 C Keyword EEBREM: do bremsstralung
1089 IF(WORD.EQ.'EEBREM ') THEN
1092 READ*, RSHMIN,RSHMAX
1093 WRITE(ITCOM,*) RSHMIN,RSHMAX
1099 C Keyword EEBEAM: do beamstralung
1100 IF(WORD.EQ.'EEBEAM ') THEN
1103 READ*, RSHMIN,RSHMAX,UPSLON,SIGZ
1104 WRITE(ITCOM,*) RSHMIN,RSHMAX,UPSLON,SIGZ
1111 C Keyword QMKKG (QMW for EXTRADIM)
1112 IF(WORD.EQ.'QMKKG ') THEN
1114 PRINT*,'KK graviton mass limits?'
1116 WRITE(ITCOM,*) QMIN,QMAX
1121 C Keyword QTKKG (QTW for EXTRADIM)
1122 IF(WORD.EQ.'QTKKG ') THEN
1124 PRINT*,'KK graviton qt limits?'
1126 WRITE(ITCOM,*) QTMIN,QTMAX
1131 C Keyword EXTRAD for EXTRADIM
1132 IF(WORD.EQ.'EXTRAD ') THEN
1134 PRINT*,'KK delta n, mass, UVCUT flag?'
1135 READ*, NEXTRAD,MASSD,UVCUT
1136 WRITE(ITCOM,*) NEXTRAD,MASSD,UVCUT
1142 IF(WORD.EQ.'MIJLIM ') THEN
1144 PRINT*,'Dijet pair mass limits i,j,min,max?'
1145 READ*, JLIM1,JLIM2,AMLIM1,AMLIM2
1146 WRITE(ITCOM,*) JLIM1,JLIM2,AMLIM1,AMLIM2
1152 IF(WORD.EQ.'MTOT ') THEN
1154 PRINT*,'Total mass limits?'
1155 READ*, EHMGMN,EHMGMX
1156 WRITE(ITCOM,*) EHMGMN,EHMGMX
1161 C Keyword SUGRHN: enter right-neutrino parameters
1162 IF(WORD.EQ.'SUGRHN ') THEN
1165 READ*, XMN3NR,XMAJNR,XANSS,XNRSS
1166 WRITE(ITCOM,*) XMN3NR,XMAJNR,XANSS,XNRSS
1172 IF(WORD.EQ.'AMSB ') THEN
1175 READ*, XM0SU,XMHSU,XTGBSU,XSMUSU
1176 WRITE(ITCOM,*) XM0SU,XMHSU,XTGBSU,XSMUSU
1185 IF(WORD.EQ.'SSBCSC ') THEN
1189 WRITE(ITCOM,*) XSBCS
1199 562 FORMAT(' Give jet limits: min1,max1,min2,max2, etc.',
1200 1' Terminate list with a /.',/,
1201 2' To fix a variable give only minimum value.',/)
1202 563 FORMAT(' Allowed range is ',F9.2,' - ',F9.2)
1203 572 FORMAT(' T or F?')
1204 584 FORMAT(' Value?')
1205 585 FORMAT(' Give M(gluino), mu, M(ha), tan(beta)')
1206 586 FORMAT(' Give M(ul), M(dr), M(ur), M(el), M(er), ')
1207 587 FORMAT(' Give M(tl), M(br), M(tr), M(taul), M(taur), A(t),',
1209 588 FORMAT(' Give PDFLIB keyword, value, ...')
1210 589 FORMAT(' Give M_0, M_(1/2), A_0, tan(beta), sgn(mu)')
1211 590 FORMAT(' Give polarizations P_L(e-), P_L(e+)')
1212 591 FORMAT(' Give optional M(cl), M(sr), M(cr), M(mul), M(mur), ')
1213 592 FORMAT(' Give optional M_1, M_2')
1214 593 FORMAT(' Give LAM, M_MES, XN5, tan(beta), sgn(mu), C_mgv')
1215 594 FORMAT(' Give optional GUT scale M_1, M_2, M_3')
1216 595 FORMAT(' Give optional GUT scale A_t, A_b, A_tau')
1217 596 FORMAT(' Give optional GUT scale m_Hd, m_Hu')
1218 597 FORMAT(' Give optional M(ul), M(dr), M(ur), M(el), M(er), ')
1219 598 FORMAT(' Give optional M(tl), M(br), M(tr), M(Ll), M(Lr), ')
1220 599 FORMAT(' Enter optional m(gravitino), ')
1221 600 FORMAT(' Enter .TRUE. if alpha_s unification is wanted:, ')
1222 601 FORMAT(' Enter Rsl,dmH_d,dmH_u,RDY,n5_1,n5_2,n5_3: ')
1223 602 FORMAT(' Enter RSHMIN,RSHMAX for bremsstrahlung: ')
1224 603 FORMAT(' Enter RSHMIN,RSHMAX,UPSLON,SIGZ for beamstrahlung: ')
1225 604 FORMAT(' Enter M_N3, M_MAJ, A_N, M_NRSS for RHN model: ')
1226 605 FORMAT(' Give M_0, M_(3/2), tan(beta), sgn(mu)')
1227 606 FORMAT(' Enter Q_max= max scale choice for SUSY BCs')
1229 C CHECK FOR ERRORS AND OVERLAPPING VARIABLES.
1230 C PRINT ERROR MESSAGE OR WARNINGS
1232 IF(LOC(2)*LOC(5).NE.0) PRINT 2001
1233 IF(LOC(2)*LOC(6).NE.0) PRINT 2001
1234 IF(LOC(15).NE.0.OR.LOC(37).NE.0) DUMY=SETTYP(0)
1239 IF(LOC(I1)*LOC(I2)*LOC(I3).NE.0) PRINT 1001
1241 1001 FORMAT(//,2X,'YOU HAVE GIVEN LIMITS FOR AN OVERLAPPING SET OF',
1242 C ' VARIABLES. SET MINIMIZING PPERP INTERVAL WILL BE USED')
1243 2001 FORMAT(//,2X,' YOU CANNOT GIVE LIMITS FOR BOTH THETA AND Y. MAKE',
1244 C ' UP YOUR MIND. JOB TERMINATED')
1245 2003 FORMAT(/,' YOU CANNOT RUN WS AND Z0 OR GAMMAS AT THE SAME TIME')
1246 2004 FORMAT(//,' PARAMETER OUT OF RANGE. ')