00001 subroutine list_isajet_init
00002 C
00003 C JTDKY = +/- unit number for decay table file.
00004 C If it is negative, decay table is not printed.
00005 C JTEVT = +/- unit number for output event file.
00006 C If it is negative, only stable particles are written on it.
00007 C JTCOM = unit number for command file.
00008 C JTLIS = unit number for listing.
00009 C
00010 IMPLICIT NONE
00011
00012 C
00013 integer istr,nevt,itotal
00014 INTEGER JTDKY,JTEVT,JTCOM,JTLIS,IFL,ILOOP,IPRT,LOK,ILOOP2
00015 INTEGER INDEC,INDEC2
00016 CHARACTER*132 ISADEC
00017 LOGICAL OK,DONE
00018 SAVE ILOOP,JTDKY,JTEVT,JTCOM,JTLIS
00019
00020 EXTERNAL ALDATA
00021 C
00022 C Initialize ISAJET
00023 C
00024 JTDKY=-1
00025 JTEVT=23
00026 JTCOM=21
00027 JTLIS=22
00028 nevt=1000
00029 itotal=0
00030 C point to standard decay tables
00031 CALL GETENV( 'ISAJET_DIR', ISADEC )
00032 IF ( ISADEC .EQ. ' ' ) THEN
00033 ISADEC = 'isadecay.dat'
00034 ELSE
00035 INDEC = INDEX ( ISADEC , ' ' )
00036 ISADEC(INDEC:INDEC+13) = '/isadecay.dat'
00037 INDEC2 = INDEX ( ISADEC , ' ' ) - 1
00038 c print *,'looking for ',ISADEC(:INDEC2)
00039 ENDIF
00040 OPEN(UNIT=1,FILE=ISADEC,STATUS='OLD')
00041 OPEN(UNIT=JTLIS,FILE='examListIsajet.lpt',STATUS='NEW')
00042 CALL ISAINI(JTDKY,JTEVT,JTCOM,JTLIS)
00043 IPRT=1
00044 C
00045 C print list of defined particles and their translations
00046 C (you have to call ISAINI first)
00047 C use bogus masses for 4th generation quarks so we can see their definitions
00048 C
00049 CALL PRTLST(JTLIS,200.,220.)
00050 C
00051
00052 return
00053 END
00054
00055 subroutine get_label( id, lb)
00056
00057 integer id
00058 character*8 LB,LABEL
00059
00060 lb = ' '
00061
00062 if ( id.ne.0) then
00063 lb = LABEL(id)
00064 endif
00065
00066 return
00067 end
1.5.1-3