1 #include "isajet/pilot.h"
4 C Copy event information into ZEVEL and call BUFOUT.
5 C If number of words required exceeds MAXLEN-8, the number
6 C of records written=no. of words/(MAXLEN-8)+1
8 #if defined(CERNLIB_IMPNONE)
11 #include "isajet/itapes.inc"
12 #include "isajet/mbgen.inc"
13 #include "isajet/keys.inc"
14 #include "isajet/idrun.inc"
15 #include "isajet/jetpar.inc"
16 #include "isajet/jetset.inc"
17 #include "isajet/jetsig.inc"
18 #include "isajet/partcl.inc"
19 #include "isajet/pjets.inc"
20 #include "isajet/pinits.inc"
21 #include "isajet/primar.inc"
22 #include "isajet/zevel.inc"
23 #include "isajet/final.inc"
24 #include "isajet/totals.inc"
25 #include "isajet/wsig.inc"
27 INTEGER I2,I1,JET,K,IEX,IL,ITA,I,NPSTA
33 CALL MOVLEI(IDVER,IZEVEL(IL),4)
35 CALL MOVLEL(KEYS(1),LZEVEL(IL),MXKEYS)
39 CALL MOVLEV(P(1),ZEVEL(IL),59)
41 CALL MOVLEV(SIGF,ZEVEL(IL),1)
49 C IF ITEVT.LT.0 WRITE ONLY STABLE PARTICLES AND FLAG
50 C BY NPTCL=-(NO. OF STABLE PARTICLES)
56 990 IF(IDCAY(I).EQ.0) NPSTA=NPSTA+1
62 CALL MOVLEV(PJETS(1,1),ZEVEL(IL),IEX)
64 CALL MOVLEI(IDJETS(1),IZEVEL(IL),NJET)
67 IF(KEYS(3).OR.KEYS(7).OR.KEYS(11)) THEN
68 CALL MOVLEV(QWJET(1),ZEVEL(IL),6)
70 CALL MOVLEV(QMW,ZEVEL(IL),16)
77 IF(KEYS(6).OR.KEYS(7)) THEN
81 CALL MOVLEV(PPAIR(1,1),ZEVEL(IL),5*NPAIR)
83 CALL MOVLEI(IDPAIR(1),IZEVEL(IL),NPAIR)
85 CALL MOVLEI(JPAIR(1),IZEVEL(IL),NPAIR)
91 CALL MOVLEI(NKINPT,IZEVEL(IL),5)
93 CALL MOVLEI(NPOM,IZEVEL(IL),1)
96 C /JETSET/ COMMON BLOCK
97 IF(NJSET.LT.1) GOTO 12
99 CALL MOVLEV(PJSET(1,I),ZEVEL(IL),5)
102 IZEVEL(IL+1)=JTYPE(I)
103 IZEVEL(IL+2)=JDCAY(I)
105 IF(IL.LE.MAXLEN-9) GO TO 50
106 IZEVEL(1)=IZEVEL(1)+1
108 IF(I.EQ.NJSET) GO TO 12
111 C /PARTCL/ COMMON BLOCK
112 C IF ITEVT.LT.0, WRITE OUT ONLY STABLE PARTICLES
113 C FLAG BY NPTCL=-(NO. OF STABLE PARTICLES)
114 C SUPPRESS ORIGIN AND DECAY INFORMATION
115 12 IF(NPTCL.EQ.0) GOTO 999
116 IF(ITEVT.GT.0) GOTO 997
117 C ONLY STABLE PARTICLES
119 IF(IDCAY(K).NE.0) GOTO 992
120 JET=IABS(IORIG(K))/1000
121 CALL MOVLEV(PPTCL(1,K),ZEVEL(IL),5)
122 IZEVEL(IL+5)=(JET*10000+IABS(IDENT(K)))*ISIGN(1,IDENT(K))
124 IF(IL.LE.MAXLEN-6) GOTO 992
125 IZEVEL(1)=IZEVEL(1)+1
127 IF(K.EQ.NPTCL) RETURN
132 C NOTE IDCAY CAN EXCEED 2**24 LIMIT OF PAIRPAK
134 CALL MOVLEV(PPTCL(1,K),ZEVEL(IL),5)
135 IZEVEL(IL+5)=IORIG(K)
136 IZEVEL(IL+6)=IDENT(K)
137 IZEVEL(IL+7)=IDCAY(K)/IPACK
138 IZEVEL(IL+8)=MOD(IDCAY(K),IPACK)
140 IF(IL.LE.MAXLEN-9) GOTO 998
141 IZEVEL(1)=IZEVEL(1)+1
143 IF(K.EQ.NPTCL) RETURN
146 IZEVEL(1)=IZEVEL(1)+1