The Init method of AliITSreconstruction has to be called by the user. This was done...
[u/mrichter/AliRoot.git] / ISAJET / isatape / prtlst.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 SUBROUTINE PRTLST(JTLIS,AMY,AMX)
3C
4C List defined particles. AMY, AMX are the masses of the
5C fourth generation quarks. If a negative mass is given,
6C then these are not listed.
7C This must be linked with ISAJET, including ALDATA.
8C
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
19C
20C Initialize SUSY masses to 0. Remember offset of 1 from KL.
21C
22 DO 100 I=22,NQLEP
23 AMLEP(I)=0.
24100 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
32C
33C Loop over IDENT's
34C
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)
52C 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