ISAJET/doc/sample.doc

   1 \newpage
   2 \section{Sample Jobs\label{SAMPLE}}
   3
   4       The simplest ISAJET job reads a user-supplied parameter file and
   5 writes a data file and a listing file. The following is an example of
   6 a parameter file which generates each type of event:
   7 \begin{verbatim}
   8 SAMPLE TWOJET JOB
   9 800.,100,2,50/
  10 TWOJET
  11 PT
  12 50,100,50,100/
  13 JETTYPE1
  14 'GL'/
  15 JETTYPE2
  16 'UP','UB','DN','DB','ST','SB'/
  17 END
  18 SAMPLE DRELLYAN JOB
  19 800.,100,2,50/
  20 DRELLYAN
  21 QMW
  22 80,100/
  23 WTYPE
  24 'W+','W-'/
  25 END
  26 SAMPLE MINBIAS JOB
  27 800.,100,2,50/
  28 MINBIAS
  29 END
  30 SAMPLE WPAIR JOB
  31 800.,100,2,50/
  32 WPAIR
  33 PT
  34 50,100,50,100/
  35 JETTYPE1
  36 'W+','W-','Z0'/
  37 JETTYPE2
  38 'W+','W-','Z0'/
  39 WMODE1
  40 'E+','E-','NUS'/
  41 WMODE2
  42 'QUARKS'/
  43 END
  44 SAMPLE HIGGS JOB FOR SSC
  45 40000,100,1,1/
  46 HIGGS
  47 QMH
  48 400,1600/
  49 HMASS
  50 800/
  51 JETTYPE1
  52 'Z0'/
  53 JETTYPE2
  54 'Z0'/
  55 WMODE1
  56 'MU+','MU-'/
  57 WMODE2
  58 'E+','E-'/
  59 PT
  60 50,20000,50,20000/
  61 END
  62 SAMPLE SUSY JOB
  63 1800,100,1,10/
  64 SUPERSYM
  65 PT
  66 50,100,50,100/
  67 JETTYPE1
  68 'GLSS','SQUARKS'/
  69 JETTYPE2
  70 'GLSS','SQUARKS'/
  71 GAUGINO
  72 60,1,40,40/
  73 SQUARK
  74 80.3,80.3,80.5,81.6,85,110/
  75 FORCE
  76 29,30,1,-1/
  77 FORCE
  78 21,29,1/
  79 FORCE
  80 22,29,2/
  81 FORCE
  82 23,29,3/
  83 FORCE
  84 24,29,4/
  85 FORCE
  86 25,29,5/
  87 FORCE
  88 26,29,6/
  89 END
  90 SAMPLE MSSM JOB FOR TEVATRON
  91 1800.,100,1,1/
  92 SUPERSYM
  93 BEAMS
  94 'P','AP'/
  95 MSSMA
  96 200,-200,500,2/
  97 MSSMB
  98 200,200,200,200,200/
  99 MSSMC
 100 200,200,200,200,200,0,0,0/
 101 JETTYPE1
 102 'GLSS'/
 103 JETTYPE2
 104 'SQUARKS'/
 105 PT
 106 100,300,100,300/
 107 END
 108 SAMPLE MSSM SUGRA JOB FOR LHC
 109 14000,100,1,10/
 110 SUPERSYM
 111 PT
 112 50,500,50,500/
 113 SUGRA
 114 247,302,-617.5,10,-1/
 115 TMASS
 116 175/
 117 END
 118 SAMPLE SUGRA HIGGS JOB USING DEFAULT QMH RANGE
 119 14000,100,20,50/
 120 HIGGS
 121 SUGRA
 122 200,200,0,2,+1/
 123 HTYPE
 124 'HA0'/
 125 JETTYPE1
 126 'GAUGINOS','SLEPTONS'/
 127 JETTYPE2
 128 'GAUGINOS','SLEPTONS'/
 129 END
 130 SAMPLE E+E- TO SUGRA JOB WITH POLARIZED BEAMS AND BREM/BEAMSTRAHLUNG
 131 500.,100,1,1/
 132 E+E-
 133 SUGRA
 134 125,125,0,3,1/
 135 TMASS
 136 175,-1,1/
 137 EPOL
 138 -.9,0./
 139 EEBEAM
 140 200.,500.,.1072,.12/
 141 JETTYPE1
 142 'ALL'/
 143 JETTYPE2
 144 'ALL'/
 145 NTRIES
 146 10000/
 147 END
 148 SAMPLE WH JOB
 149 2000,100,0,0/
 150 WHIGGS
 151 BEAMS
 152 'P','AP'/
 153 HMASS
 154 100./
 155 JETTYPE1
 156 'W+','W-','HIGGS'/
 157 JETTYPE2
 158 'W+','W-','HIGGS'/
 159 WMODE1
 160 'ALL'/
 161 WMODE2
 162 'ALL'/
 163 PT
 164 10,300,10,300/
 165 END
 166 SAMPLE EXTRA DIMENSIONS JOB
 167 14000,100,1,100/
 168 EXTRADIM
 169 QMW
 170 5,1000/
 171 QTW
 172 500,1000/
 173 EXTRAD
 174 2,1000,.FALSE./
 175 END
 176 SAMPLE ZJJ JOB AT LHC
 177 14000,100,1,100/
 178 ZJJ
 179 PT
 180 20,7000,20,7000,20,7000/
 181 MIJLIM
 182 0,0,20,7000/
 183 MTOT
 184 100,500/
 185 NSIGMA
 186 200/
 187 NTRIES
 188 10000/
 189 END
 190 STOP
 191 \end{verbatim}
 192 \noindent See Section~\ref{INPUT} of this manual for a complete list
 193 of the possible commands in a parameter file. Note that all input to
 194 ISAJET must be in {\it UPPER} case only.
 195
 196       Subroutine RDTAPE is supplied to read events from an ISAJET data
 197 file, which is a machine-dependent binary file. It restores the event
 198 data to the FORTRAN common blocks described in Section~\ref{OUTPUT}.
 199 The skeleton of an analysis job using HBOOK and PAW from the CERN
 200 Program Library is provided in patch ISAPLT but is not otherwise
 201 supported. A Zebra output format based on code from the D0
 202 Collaboration is also provided in patch ISAZEB; see the separate
 203 documentation in patch ISZTEXT.
 204
 205 \subsection{DEC VMS}
 206
 207       On a VAX or ALPHA running VMS, ISAJET can be compiled by
 208 executing the .COM file contained in P=ISAUTIL,D=MAKEVAX. Extract this
 209 deck as ISAMAKE.COM and type
 210 \begin{verbatim}
 211 @ISAMAKE
 212 \end{verbatim}
 213 This will run YPATCHY with the pilot patches described in
 214 Section~\ref{PATCHY} and the VAX flag to extract the source code,
 215 decay table, and documentation. The source code is compiled and made
 216 into a library, which is linked with the following main program,
 217 \begin{verbatim}
 218       PROGRAM ISARUN
 219 C          MAIN PROGRAM FOR ISAJET
 220       OPEN(UNIT=1,STATUS='OLD',FORM='FORMATTED',READONLY)
 221       OPEN(UNIT=2,STATUS='NEW',FORM='UNFORMATTED')
 222       OPEN(UNIT=3,STATUS='OLD',FORM='FORMATTED')
 223       OPEN(UNIT=4,STATUS='NEW',FORM='FORMATTED')
 224       CALL ISAJET(-1,2,3,4)
 225       STOP
 226       END
 227 \end{verbatim}
 228 to produce ISAJET.EXE. Two other executables, ISASUSY.EXE and
 229 ISASUGRA.EXE, will also be produced to calculate SUSY masses and decay
 230 modes without generating events. Temporary files can be removed by
 231 typing
 232 \begin{verbatim}
 233 @ISAMAKE CLEAN
 234 \end{verbatim}
 235
 236       Create an input file \verb|JOBNAME.PAR| following the examples
 237 above or the instructions in Section~\ref{INPUT} and run ISAJET with
 238 the command
 239 \begin{verbatim}
 240 @ISAJET JOBNAME
 241 \end{verbatim}
 242 using the ISAJET.COM file contained P=ISAUTIL,D=RUNVAX. This will
 243 create a binary output file \verb|JOBNAME.DAT| and a listing file
 244 \verb|JOBNAME.LIS|. Analyze the output data using the commands
 245 described in Section~\ref{TAPE}.
 246
 247       There is also an simple interactive interface to ISAJET which
 248 will prompt the user for commands, write a parameter file, and
 249 optionally execute it.
 250
 251 \subsection{IBM VM/CMS}
 252
 253       On an IBM mainframe running VM/CMS, run YPATCHY with the pilot
 254 patches described in Section~\ref{PATCHY} and the IBM flag to extract
 255 the source code, decay table, and documentation. Compile the source
 256 code and link it with the main program
 257 \begin{verbatim}
 258       PROGRAM ISARUN
 259 C          MAIN PROGRAM FOR ISAJET
 260       OPEN(UNIT=1,STATUS='OLD',FORM='FORMATTED')
 261       OPEN(UNIT=2,STATUS='NEW',FORM='UNFORMATTED')
 262       OPEN(UNIT=3,STATUS='OLD',FORM='FORMATTED')
 263       OPEN(UNIT=4,STATUS='NEW',FORM='FORMATTED')
 264       CALL ISAJET(-1,2,3,4)
 265       STOP
 266       END
 267 \end{verbatim}
 268 to make ISAJET MODULE.
 269
 270       Create a file called \verb|JOBNAME INPUT| containing ISAJET
 271 input commands following the examples above or the instructions in
 272 Section~\ref{INPUT}. Then run ISAJET using ISAJET EXEC, which is
 273 contained in P=ISAUTIL,D=RUNIBM.  The events will be produced on
 274 \verb|JOBNAME DATA A| and the listing on \verb|JOBNAME OUTPUT A|.
 275
 276 \subsection{Unix}
 277
 278       The Makefile contained in P=ISAUTIL,D=MAKEUNIX has been tested
 279 on DEC Ultrix, Hewlett Packard HP-UX, IBM RS/6000 AIX, Linux, Silicon
 280 Graphics IRIX, Sun SunOS, and Sun Solaris. It should work with minor
 281 modifications on almost any Unix system with /bin/csh, \verb|ypatchy|
 282 or \verb|nypatchy|, and a reasonable Fortran 77 compiler. Extract the
 283 Makefile and edit it, changing the installation parameters to reflect
 284 your system. Note in particular that CERNlib is usually compiled with
 285 underscores postpended to all external names; you must choose the
 286 appropriate compiler option if you intend to link with it. Then type
 287 \begin{verbatim}
 288 make
 289 \end{verbatim}
 290 This should produce an executable \verb|isajet.x| for the event
 291 generator, which links the code with the following main program:
 292 \begin{verbatim}
 293       PROGRAM RUNJET
 294       CHARACTER*60 FNAME
 295       READ 1000, FNAME
 296 1000  FORMAT(A)
 297       PRINT 1020, FNAME
 298 1020  FORMAT(1X,'Data file      = ',A)
 299       OPEN(2,FILE=FNAME,STATUS='NEW',FORM='UNFORMATTED')
 300       READ 1000, FNAME
 301       PRINT 1030, FNAME
 302 1030  FORMAT(1X,'Parameter file = ',A)
 303       OPEN(3,FILE=FNAME,STATUS='OLD',FORM='FORMATTED')
 304       READ 1000, FNAME
 305       PRINT 1040, FNAME
 306 1040  FORMAT(1X,'Listing file   = ',A)
 307       OPEN(4,FILE=FNAME,STATUS='NEW',FORM='FORMATTED')
 308       READ 1000, FNAME
 309       OPEN(1,FILE=FNAME,STATUS='OLD',FORM='FORMATTED')
 310       CALL ISAJET(-1,2,3,4)
 311       STOP
 312       END
 313 \end{verbatim}
 314 Two other executables, \verb|isasusy.x| and \verb|isasugra.x|, will
 315 also be produced to calculate SUSY masses and decay modes without
 316 generating events. Type
 317 \begin{verbatim}
 318 make clean
 319 \end{verbatim}
 320 to delete the temporary files.
 321
 322       Most Unix systems do not allow two jobs to read the same decay
 323 table file at the same time. The shell script in P=ISAUTIL,D=RUNUNIX
 324 copies the decay table to a temporary file to avoid this problem.
 325 Extract this file as \verb|isajet|. Create an input file
 326 \verb|jobname.par| following the examples above or the instructions in
 327 Section~\ref{INPUT} and run ISAJET with the command
 328 \begin{verbatim}
 329 isajet jobname
 330 \end{verbatim}
 331 This will create a binary output file \verb|jobname.dat| and a listing
 332 file \verb|jobname.lis|. Analyze the output data using the commands
 333 described in Section~\ref{TAPE}.
 334
 335       This section only describes running ISAJET as a standalone
 336 program and generating output in machine-dependent binary form. The
 337 user may elect to analyze events as they are generated; this is
 338 discussed in Section~\ref{MAIN} of this manual.