The Init method of AliITSreconstruction has to be called by the user. This was done...
[u/mrichter/AliRoot.git] / ISAJET / isatape / prtlst.F
1 #include "isajet/pilot.h"
2       SUBROUTINE PRTLST(JTLIS,AMY,AMX)
3 C
4 C          List defined particles. AMY, AMX are the masses of the 
5 C          fourth generation quarks. If a negative mass is given,
6 C          then these are not listed.
7 C          This must be linked with ISAJET, including ALDATA.
8 C
9 #if defined(CERNLIB_IMPNONE)
10       IMPLICIT NONE
11 #endif
12 #include "isajet/qlmass.inc"
13 #include "isajet/wcon.inc"
14       INTEGER JTLIS
15       REAL AMY,AMX
16       INTEGER IFL1,IFL2,IFL3,JSPIN,INDEX,I,ID
17       REAL AM,CG,AMASS,CHARGE
18       CHARACTER*8 LB,LABEL
19 C
20 C          Initialize SUSY masses to 0. Remember offset of 1 from KL.
21 C
22       DO 100 I=22,NQLEP
23         AMLEP(I)=0.
24 100   CONTINUE
25       AMLEP(7)=AMY
26       AMLEP(8)=AMX
27       CALL FLAVOR(80,IFL1,IFL2,IFL3,JSPIN,INDEX)
28       AMLEP(INDEX)=WMASS(2)
29       CALL FLAVOR(90,IFL1,IFL2,IFL3,JSPIN,INDEX)
30       AMLEP(INDEX)=WMASS(4)
31       WRITE(JTLIS,101) AMY,AMX
32 C
33 C          Loop over IDENT's
34 C
35       DO 200 I=1,40000
36         ID = I
37         CALL FLAVOR(ID,IFL1,IFL2,IFL3,JSPIN,INDEX)
38         IF(AMX.LT.0..OR.AMY.LT.0.) THEN
39           IF(IABS(IFL1).GT.6.OR.IABS(IFL2).GT.6.OR.IABS(IFL3).GT.6)
40      $    GO TO 200
41         ENDIF
42         IF(INDEX.GT.0) THEN
43           LB = LABEL(ID)
44           IF(LB.NE.'ERR') THEN
45             AM = AMASS(ID)
46             CG = CHARGE(ID)
47             WRITE(JTLIS,102) ID,LB,AM,CG,IFL1,IFL2,IFL3,JSPIN,INDEX
48           ENDIF
49         ENDIF
50         ID = -I
51         CALL FLAVOR(ID,IFL1,IFL2,IFL3,JSPIN,INDEX)
52 C          Eliminate bad ID's:
53         IF(INDEX.GT.0) THEN
54           LB = LABEL(ID)
55           IF(LB.NE.'ERR') THEN
56             AM = AMASS(ID)
57             CG = CHARGE(ID)
58             WRITE(JTLIS,102) ID,LB,AM,CG,IFL1,IFL2,IFL3,JSPIN,INDEX
59           ENDIF
60         ENDIF
61  200  CONTINUE
62       RETURN
63  101  FORMAT(10X,'ISAJET PARTICLES, M(Y) =',F10.3,'  M(X) =',F10.3//
64      1 5X,'ID',4X,'PARTICLE',8X,'MASS',4X,'CHARGE',
65      2 4X,'---FLAVOR---',4X,'SPIN',4X,'INDEX')
66  102  FORMAT(1X,I6,4X,A8,F12.6,F10.2,4X,3I4,I8,I9)
67       END