First commit.
[u/mrichter/AliRoot.git] / 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.