Merging the VirtualMC branch to the main development branch (HEAD)
[u/mrichter/AliRoot.git] / ISAJET / doc / sample.doc
CommitLineData
0795afa3 1\newpage
2\section{Sample Jobs\label{SAMPLE}}
3
4 The simplest ISAJET job reads a user-supplied parameter file and
5writes a data file and a listing file. The following is an example of
6a parameter file which generates each type of event:
7\begin{verbatim}
8SAMPLE TWOJET JOB
9800.,100,2,50/
10TWOJET
11PT
1250,100,50,100/
13JETTYPE1
14'GL'/
15JETTYPE2
16'UP','UB','DN','DB','ST','SB'/
17END
18SAMPLE DRELLYAN JOB
19800.,100,2,50/
20DRELLYAN
21QMW
2280,100/
23WTYPE
24'W+','W-'/
25END
26SAMPLE MINBIAS JOB
27800.,100,2,50/
28MINBIAS
29END
30SAMPLE WPAIR JOB
31800.,100,2,50/
32WPAIR
33PT
3450,100,50,100/
35JETTYPE1
36'W+','W-','Z0'/
37JETTYPE2
38'W+','W-','Z0'/
39WMODE1
40'E+','E-','NUS'/
41WMODE2
42'QUARKS'/
43END
44SAMPLE HIGGS JOB FOR SSC
4540000,100,1,1/
46HIGGS
47QMH
48400,1600/
49HMASS
50800/
51JETTYPE1
52'Z0'/
53JETTYPE2
54'Z0'/
55WMODE1
56'MU+','MU-'/
57WMODE2
58'E+','E-'/
59PT
6050,20000,50,20000/
61END
62SAMPLE SUSY JOB
631800,100,1,10/
64SUPERSYM
65PT
6650,100,50,100/
67JETTYPE1
68'GLSS','SQUARKS'/
69JETTYPE2
70'GLSS','SQUARKS'/
71GAUGINO
7260,1,40,40/
73SQUARK
7480.3,80.3,80.5,81.6,85,110/
75FORCE
7629,30,1,-1/
77FORCE
7821,29,1/
79FORCE
8022,29,2/
81FORCE
8223,29,3/
83FORCE
8424,29,4/
85FORCE
8625,29,5/
87FORCE
8826,29,6/
89END
90SAMPLE MSSM JOB FOR TEVATRON
911800.,100,1,1/
92SUPERSYM
93BEAMS
94'P','AP'/
95MSSMA
96200,-200,500,2/
97MSSMB
98200,200,200,200,200/
99MSSMC
100200,200,200,200,200,0,0,0/
101JETTYPE1
102'GLSS'/
103JETTYPE2
104'SQUARKS'/
105PT
106100,300,100,300/
107END
108SAMPLE MSSM SUGRA JOB FOR LHC
10914000,100,1,10/
110SUPERSYM
111PT
11250,500,50,500/
113SUGRA
114247,302,-617.5,10,-1/
115TMASS
116175/
117END
118SAMPLE SUGRA HIGGS JOB USING DEFAULT QMH RANGE
11914000,100,20,50/
120HIGGS
121SUGRA
122200,200,0,2,+1/
123HTYPE
124'HA0'/
125JETTYPE1
126'GAUGINOS','SLEPTONS'/
127JETTYPE2
128'GAUGINOS','SLEPTONS'/
129END
130SAMPLE E+E- TO SUGRA JOB WITH POLARIZED BEAMS AND BREM/BEAMSTRAHLUNG
131500.,100,1,1/
132E+E-
133SUGRA
134125,125,0,3,1/
135TMASS
136175,-1,1/
137EPOL
138-.9,0./
139EEBEAM
140200.,500.,.1072,.12/
141JETTYPE1
142'ALL'/
143JETTYPE2
144'ALL'/
145NTRIES
14610000/
147END
148SAMPLE WH JOB
1492000,100,0,0/
150WHIGGS
151BEAMS
152'P','AP'/
153HMASS
154100./
155JETTYPE1
156'W+','W-','HIGGS'/
157JETTYPE2
158'W+','W-','HIGGS'/
159WMODE1
160'ALL'/
161WMODE2
162'ALL'/
163PT
16410,300,10,300/
165END
166SAMPLE EXTRA DIMENSIONS JOB
16714000,100,1,100/
168EXTRADIM
169QMW
1705,1000/
171QTW
172500,1000/
173EXTRAD
1742,1000,.FALSE./
175END
176SAMPLE ZJJ JOB AT LHC
17714000,100,1,100/
178ZJJ
179PT
18020,7000,20,7000,20,7000/
181MIJLIM
1820,0,20,7000/
183MTOT
184100,500/
185NSIGMA
186200/
187NTRIES
18810000/
189END
190STOP
191\end{verbatim}
192\noindent See Section~\ref{INPUT} of this manual for a complete list
193of the possible commands in a parameter file. Note that all input to
194ISAJET must be in {\it UPPER} case only.
195
196 Subroutine RDTAPE is supplied to read events from an ISAJET data
197file, which is a machine-dependent binary file. It restores the event
198data to the FORTRAN common blocks described in Section~\ref{OUTPUT}.
199The skeleton of an analysis job using HBOOK and PAW from the CERN
200Program Library is provided in patch ISAPLT but is not otherwise
201supported. A Zebra output format based on code from the D0
202Collaboration is also provided in patch ISAZEB; see the separate
203documentation in patch ISZTEXT.
204
205\subsection{DEC VMS}
206
207 On a VAX or ALPHA running VMS, ISAJET can be compiled by
208executing the .COM file contained in P=ISAUTIL,D=MAKEVAX. Extract this
209deck as ISAMAKE.COM and type
210\begin{verbatim}
211@ISAMAKE
212\end{verbatim}
213This will run YPATCHY with the pilot patches described in
214Section~\ref{PATCHY} and the VAX flag to extract the source code,
215decay table, and documentation. The source code is compiled and made
216into a library, which is linked with the following main program,
217\begin{verbatim}
218 PROGRAM ISARUN
219C 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}
228to produce ISAJET.EXE. Two other executables, ISASUSY.EXE and
229ISASUGRA.EXE, will also be produced to calculate SUSY masses and decay
230modes without generating events. Temporary files can be removed by
231typing
232\begin{verbatim}
233@ISAMAKE CLEAN
234\end{verbatim}
235
236 Create an input file \verb|JOBNAME.PAR| following the examples
237above or the instructions in Section~\ref{INPUT} and run ISAJET with
238the command
239\begin{verbatim}
240@ISAJET JOBNAME
241\end{verbatim}
242using the ISAJET.COM file contained P=ISAUTIL,D=RUNVAX. This will
243create a binary output file \verb|JOBNAME.DAT| and a listing file
244\verb|JOBNAME.LIS|. Analyze the output data using the commands
245described in Section~\ref{TAPE}.
246
247 There is also an simple interactive interface to ISAJET which
248will prompt the user for commands, write a parameter file, and
249optionally execute it.
250
251\subsection{IBM VM/CMS}
252
253 On an IBM mainframe running VM/CMS, run YPATCHY with the pilot
254patches described in Section~\ref{PATCHY} and the IBM flag to extract
255the source code, decay table, and documentation. Compile the source
256code and link it with the main program
257\begin{verbatim}
258 PROGRAM ISARUN
259C 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}
268to make ISAJET MODULE.
269
270 Create a file called \verb|JOBNAME INPUT| containing ISAJET
271input commands following the examples above or the instructions in
272Section~\ref{INPUT}. Then run ISAJET using ISAJET EXEC, which is
273contained 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
279on DEC Ultrix, Hewlett Packard HP-UX, IBM RS/6000 AIX, Linux, Silicon
280Graphics IRIX, Sun SunOS, and Sun Solaris. It should work with minor
281modifications on almost any Unix system with /bin/csh, \verb|ypatchy|
282or \verb|nypatchy|, and a reasonable Fortran 77 compiler. Extract the
283Makefile and edit it, changing the installation parameters to reflect
284your system. Note in particular that CERNlib is usually compiled with
285underscores postpended to all external names; you must choose the
286appropriate compiler option if you intend to link with it. Then type
287\begin{verbatim}
288make
289\end{verbatim}
290This should produce an executable \verb|isajet.x| for the event
291generator, which links the code with the following main program:
292\begin{verbatim}
293 PROGRAM RUNJET
294 CHARACTER*60 FNAME
295 READ 1000, FNAME
2961000 FORMAT(A)
297 PRINT 1020, FNAME
2981020 FORMAT(1X,'Data file = ',A)
299 OPEN(2,FILE=FNAME,STATUS='NEW',FORM='UNFORMATTED')
300 READ 1000, FNAME
301 PRINT 1030, FNAME
3021030 FORMAT(1X,'Parameter file = ',A)
303 OPEN(3,FILE=FNAME,STATUS='OLD',FORM='FORMATTED')
304 READ 1000, FNAME
305 PRINT 1040, FNAME
3061040 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}
314Two other executables, \verb|isasusy.x| and \verb|isasugra.x|, will
315also be produced to calculate SUSY masses and decay modes without
316generating events. Type
317\begin{verbatim}
318make clean
319\end{verbatim}
320to delete the temporary files.
321
322 Most Unix systems do not allow two jobs to read the same decay
323table file at the same time. The shell script in P=ISAUTIL,D=RUNUNIX
324copies the decay table to a temporary file to avoid this problem.
325Extract this file as \verb|isajet|. Create an input file
326\verb|jobname.par| following the examples above or the instructions in
327Section~\ref{INPUT} and run ISAJET with the command
328\begin{verbatim}
329isajet jobname
330\end{verbatim}
331This will create a binary output file \verb|jobname.dat| and a listing
332file \verb|jobname.lis|. Analyze the output data using the commands
333described in Section~\ref{TAPE}.
334
335 This section only describes running ISAJET as a standalone
336program and generating output in machine-dependent binary form. The
337user may elect to analyze events as they are generated; this is
338discussed in Section~\ref{MAIN} of this manual.