![]() |
HepPDT Reference DocumentationHepPDT |
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