1 #include "isajet/pilot.h"
4 C Print initial conditions and limits for generating jets
6 #if defined(CERNLIB_IMPNONE)
9 #include "isajet/itapes.inc"
10 #include "isajet/force.inc"
11 #include "isajet/mbgen.inc"
12 #include "isajet/qcdpar.inc"
13 #include "isajet/qlmass.inc"
14 #include "isajet/wcon.inc"
15 #include "isajet/jetpar.inc"
16 #include "isajet/keys.inc"
17 #include "isajet/kkgrav.inc"
18 #include "isajet/frgpar.inc"
19 #include "isajet/nodcay.inc"
20 #include "isajet/prtout.inc"
21 #include "isajet/seed.inc"
22 #include "isajet/types.inc"
23 #include "isajet/q1q2.inc"
24 #include "isajet/jetlim.inc"
25 #include "isajet/primar.inc"
26 #include "isajet/ptpar.inc"
27 #include "isajet/idrun.inc"
28 #include "isajet/dylim.inc"
29 #include "isajet/hcon.inc"
30 #include "isajet/isloop.inc"
31 #include "isajet/xmssm.inc"
33 INTEGER I,II,K,NPRT,I1,I2,I3,J1,INDEX,IQ,KK,KKK,NN,N0J
35 CHARACTER*8 LSTRUC(6),LMODE(5),STUF(6),IDFMT(2)
45 CHARACTER*8 LPRSS(NPRSS)
46 DATA LSTRUC/'OWENS','BAIER','EICHTEN','DUKE','CTEQ2L','CTEQ3L'/
47 DATA WTITL/'GM','W+','W-','Z0'/
49 DATA IDPRSS/21,22,23,24,25,26,41,42,43,44,45,46,
50 $31,32,33,34,35,36,52,54,56,
51 $29,30,40,50,60,39,49,82,83,84,86/
56 1000 FORMAT('1',44('*')/' *',42X,'*'/
58 C ' *',42X,'*'/' ',44('*')/)
60 C Print title, reaction, energy, number, run id
61 WRITE(ITLIS,1010) TITLE
62 1010 FORMAT(/11X,10A8)
63 WRITE(ITLIS,1020) NJET
64 1020 FORMAT(/2X,'NUMBER OF JETS TO BE GENERATED PER EVENT',I3)
66 IDFMT(I)=LABEL(IDIN(I))
68 WRITE(ITLIS,1030) IDFMT(1),IDFMT(2),ECM
69 1030 FORMAT(/2X,A8,' ON ',A8,' AT COM ENERGY',E15.4)
70 WRITE(ITLIS,1040) REAC,NEVENT
71 1040 FORMAT(/2X,'NUMBER OF ',A8,' EVENTS TO BE GENERATED',I10)
72 IF(NEVPRT.GT.0) WRITE(ITLIS,1050) NEVPRT,NJUMP
73 1050 FORMAT(/2X,'PRINT A MAXIMUM OF ',I6,
74 C ' EVENTS SKIPPING ',I6,' EVENTS AT A TIME')
76 1060 FORMAT(/2X,'RUN ID ',2I10.6)
78 C Print W/Higgs parameters
80 IF(KEYS(3).OR.KEYS(7).OR.KEYS(11)) THEN
84 IF(.NOT.GODY(I)) GOTO 200
94 IF(IHTYPE.EQ.82) STUF(II)='HL0'
95 IF(IHTYPE.EQ.83) STUF(II)='HH0'
96 IF(IHTYPE.EQ.84) STUF(II)='HA0'
98 WRITE(ITLIS,2000) (STUF(K),K=1,II)
99 2000 FORMAT(//10X,'PARAMETERS FOR',4(2X,A8))
100 WRITE(ITLIS,2010) QMIN,QMAX
101 2010 FORMAT(' MASS LIMITS',15X,2E15.4)
102 WRITE(ITLIS,3020) QTMIN,QTMAX
103 WRITE(ITLIS,3030) THWMIN,THWMAX
104 WRITE(ITLIS,3040) PHWMIN,PHWMAX
105 WRITE(ITLIS,3050) YWMIN,YWMAX
106 WRITE(ITLIS,3060) XWMIN,XWMAX
128 IF(NPRT.NE.0) WRITE(ITLIS,2070) (STUF(K),K=1,NPRT)
131 C Print jet parameters
143 3000 FORMAT(//10X,'JET NO.',I3,/)
144 WRITE(ITLIS,3010) PMIN(I),PMAX(I)
145 3010 FORMAT(' MOMENTUM LIMITS',11X,2E15.4)
146 WRITE(ITLIS,3020) PTMIN(I),PTMAX(I)
147 3020 FORMAT(' TRANSVERSE MOMENTUM LIMITS',2E15.4)
148 WRITE(ITLIS,3030) THMIN(I),THMAX(I)
149 3030 FORMAT(' THETA LIMITS',14X,2E15.4)
150 WRITE(ITLIS,3040) PHIMIN(I),PHIMAX(I)
151 3040 FORMAT(' PHI LIMITS',16X,2E15.4)
152 WRITE(ITLIS,3050) YJMIN(I),YJMAX(I)
153 3050 FORMAT(' RAPIDITY (Y) LIMITS',7X,2E15.4)
154 WRITE(ITLIS,3060) XJMIN(I),XJMAX(I)
155 3060 FORMAT(' FEYNMAN X LIMITS',10X,2E15.4)
156 IF(.NOT.GOALL(I)) WRITE(ITLIS,3070) (JETYP(K,I),K=1,NJTTYP(I))
157 3070 FORMAT(' JET TYPE',23X,A8,1X,A8,1X,A8,1X,A8,1X,A8)
158 IF((KEYS(6).OR.KEYS(7)).AND..NOT.ALLWW(I))
159 $ WRITE(ITLIS,3080) (WWTYP(K,I),K=1,NWWTYP(I))
160 3080 FORMAT(' DECAY MODES',20X,A8,1X,A8,1X,A8,1X,A8,1X,A8)
183 IF(NPRT.EQ.0) GOTO 310
184 WRITE(ITLIS,2070) (STUF(K),K=1,NPRT)
185 2070 FORMAT(/5X,'FOLLOWING PARAMETERS HAVE BEEN FIXED ',6A8)
188 C Print structure functions, QCD parameters, W parameters,
189 C and other parameters changed from their default values.
191 IF(KEYS(1).OR.KEYS(3).OR.KEYS(5).OR.KEYS(6).OR.KEYS(7)
192 $.OR.KEYS(8).OR.KEYS(9).OR.KEYS(11).OR.KEYS(12)) THEN
194 WRITE(ITLIS,4000) LSTRUC(ISTRUC)
195 4000 FORMAT(//1X,A8,' STRUCTURE FUNCTIONS')
196 ELSEIF(ISTRUC.EQ.-999) THEN
198 4001 FORMAT(//1X,' PDFLIB STRUCTURE FUNCTIONS')
201 WRITE(ITLIS,4010) ALAM,CUTJET
202 4010 FORMAT(//' QCD LAMBDA =',F10.4,10X,'JET CUTOFF MASS =',F10.3)
203 WRITE(ITLIS,4020) AMLEP(6),AMLEP(7),AMLEP(8)
204 4020 FORMAT(/' HEAVY QUARK MASSES =',3F8.2)
205 IF(LOC(36).NE.0) THEN
206 CALL FLAVOR(80,I1,I2,I3,J1,INDEX)
207 WRITE(ITLIS,4030) (AMLEP(INDEX+K),K=1,9)
208 4030 FORMAT(/' HIGGS MASSES =',6F8.2/15X,3F8.2)
211 IF(KEYS(5).AND..NOT.GOMSSM) THEN
216 WRITE(ITLIS,4040) (AM(KK),KK=1,6)
217 4040 FORMAT(//' SQUARK MASSES ',7F9.2)
218 WRITE(ITLIS,4050) (AML(KK),KK=1,6)
219 4050 FORMAT(' SLEPTON MASSES ',7F9.2)
224 WRITE(ITLIS,4060) (AM(KK),KK=1,4)
225 4060 FORMAT(' GAUGINO MASSES (WITH SIGNS) ',7F9.2)
226 IF(LOC(44).NE.0) WRITE(ITLIS,4070) XGENSS(9),
227 $ (XGENSS(KK),KK=1,8)
228 4070 FORMAT(/' FRAG. PARAM. XGENSS = ',9F8.3)
233 LPRSS(I)=LABEL(IDPRSS(I))
234 AMPRSS(I)=AMASS(IDPRSS(I))
236 WRITE(ITLIS,4100) (LPRSS(I),AMPRSS(I),I=1,NPRSS)
237 4100 FORMAT(/' MSSM MASSES (WITHOUT SIGNS):'/
238 $ 10(' M(',A5,') = ',F10.3,5X,'M(',A5,') = ',F10.3,
239 $ 5X,'M(',A5,') = ',F10.3/),
240 $ ' M(',A5,') = ',F10.3,5X,'M(',A5,') = ',F10.3)
242 WRITE(ITLIS,4110) XTBSS,XMUSS,XATSS
243 4110 FORMAT(/' OTHER MSSM PARAMETERS:'/
244 $ ' TAN(BETA) = ',F10.3,5X,'MU = ',F10.3,5X,'A_t = ',F10.3)
247 WRITE(ITLIS,4120) XM0SU,XMHSU,XA0SU,XTGBSU,XSMUSU
248 4120 FORMAT(/' DETERMINED FROM SUGRA INPUT:'/
249 $ ' M_0 =',F10.3,' M_(1/2) =',F10.3,' A_0 =',F10.3/
250 $ ' TAN(BETA) =',F10.4,' SGN(MU) =',F6.1)
256 IF(LOC(30).NE.0) WRITE(ITLIS,4200) XGEN
257 4200 FORMAT(/' FRAGMENTATION PARAMETER XGEN =',6F8.3)
258 IF(LOC(31).NE.0) WRITE(ITLIS,4210) SIGQT
259 4210 FORMAT(/' FRAGMENTATION PARAMETER SIGQT =',F8.3)
260 IF(KEYS(2).OR.KEYS(3).OR.KEYS(6).OR.KEYS(7)) THEN
261 WRITE(ITLIS,4220) SIN2W,WMASS(2),WMASS(4),WGAM(2),WGAM(4)
262 4220 FORMAT(//' WEINBERG MODEL',5X,'SIN**2(THETA-W)=',F8.4/
263 $ ' MASSES = ',F8.2,',',F8.2,' WIDTHS = ',F8.3,',',F8.3)
266 IF(.NOT.STDDY) WRITE(ITLIS,4230) CUTOFF,CUTPOW
267 4230 FORMAT(/' CUTOFF FUNCTION IS QT**2=',E11.4,'*Q**',E11.4)
268 IF(LOC(50).NE.0) WRITE(ITLIS,4240) WFUDGE
269 4240 FORMAT(/' W fudge factor (WFUDGE) = ',F8.3)
272 WRITE(ITLIS,4250) HMASS,HGAM
273 4250 FORMAT(/' HIGGS MASS = ',F8.2,' WIDTH = ',F8.3)
275 WRITE(ITLIS,4260) XSEED
276 4260 FORMAT(/' SEED FOR RANDOM NUMBER GENERATOR = ',A24)
277 IF(LOC(13).NE.0) WRITE(ITLIS,4270) FRPAR
278 4270 FORMAT(//' FRAGMENTATION PARAMETERS ',8E11.3,2(/,26X,8E11.3))
279 IF(LOC(34).NE.0) WRITE(ITLIS,4280) MNPOM,MXPOM
280 4280 FORMAT(//' NUMBER OF POMERONS =',I4,' TO',I4)
282 C Print KKG parameters
284 WRITE(ITLIS,4291) NEXTRAD
285 WRITE(ITLIS,4292) MASSD
286 WRITE(ITLIS,4296) SURFD
287 WRITE(ITLIS,4295) KKGSD
288 WRITE(ITLIS,4297) UVCUT
289 4291 FORMAT(//' NB EXTRA-DIMENSIONS',7X,I4)
290 4292 FORMAT(' SCALE M_D',15X,E15.4)
291 4295 FORMAT(' KKGSD FACTOR',12X,E15.4)
292 4296 FORMAT(' SD SURFACE',14X,E15.4)
293 4297 FORMAT(' UV CUTOFF',17X,I4)
298 IF(LOC(46).NE.0) THEN
300 4300 FORMAT(//' SECONDARY W DECAY MODES:')
302 WRITE(ITLIS,4310) (WMODES(KKK,1),KKK=1,NN)
303 4310 FORMAT(' W+ --> '
304 $ ,A8,1X,A8,1X,A8,1X,A8,1X,A8,1X,A8,1X,A8,1X,A8)
306 WRITE(ITLIS,4320) (WMODES(KKK,2),KKK=1,NN)
307 4320 FORMAT(' W- --> '
308 $ ,A8,1X,A8,1X,A8,1X,A8,1X,A8,1X,A8,1X,A8,1X,A8)
309 WRITE(ITLIS,4330) (WMODES(KKK,3),KKK=1,NN)
310 4330 FORMAT(' Z0 --> '
311 $ ,A8,1X,A8,1X,A8,1X,A8,1X,A8,1X,A8,1X,A8,1X,A8)
314 IF(NODCAY) WRITE(ITLIS,4400)
315 4400 FORMAT(//' NO DECAYS OF HADRONS WILL BE GENERATED')
316 IF(NOETA) WRITE(ITLIS,4410)
317 4410 FORMAT(//' NO DECAYS OF ETAS WILL BE GENERATED')
318 IF(NOPI0) WRITE(ITLIS,4420)
319 4420 FORMAT(//' NO DECAYS OF PI0S WILL BE GENERATED')
320 IF(NONUNU) WRITE(ITLIS,4430)
321 4430 FORMAT(//' NO DECAYS Z0---->NU+NU WILL BE GENERATED')
322 IF(NOEVOL) WRITE(ITLIS,4440)
323 4440 FORMAT(//' NO QCD JET EVOLUTION WILL BE DONE')
324 IF(NOHADR) WRITE(ITLIS,4450)
325 4450 FORMAT(//' NO JET HADRONIZATION WILL BE DONE')
326 IF(GOGMSB.AND.NOGRAV) WRITE(ITLIS,4460)
327 4460 FORMAT(//' NO GRAVITINO DECAYS WILL BE GENERATED')
329 C Print forced decay modes and M.E. flag
333 4500 FORMAT(//7X,'FORCED DECAY MODES (MATRIX ELEMENT FLAGS)'/
334 $ 2X,'PART',6X,'DECAY MODE')
336 IF(IFORCE(I).EQ.0) GOTO 450
340 IF(MFORCE(K,I).EQ.0) GO TO 451
341 LMODE(K)=LABEL(MFORCE(K,I))
343 WRITE(ITLIS,4510) L0,(LMODE(K),K=1,5),MEFORC(I)
344 4510 FORMAT(2X,6A10,'(M.E. =',I5,')')
348 C Print multiple evolution/fragmentation information
350 IF(NEVOLV.NE.1.OR.NFRGMN.NE.1) THEN
351 WRITE(ITLIS,4600) NEVOLV,NFRGMN
353 $ ' MULTIPLE EVOLUTION AND FRAGMENTATION VERSION'/
354 $ ' EVENTS WILL BE EVOLVED',I6,' TIMES'/
355 $ ' AND FRAGMENTED ',I6,' TIMES'/)