]>
Commit | Line | Data |
---|---|---|
0795afa3 | 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. |