HepPDT Reference Documentation

HepPDT

examListIsajetInterface.F

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

Generated on Fri Dec 4 14:05:23 2009 for HepPDT by  doxygen 1.4.7