Merging the VirtualMC branch to the main development branch (HEAD)
[u/mrichter/AliRoot.git] / ISAJET / doc / output.doc
CommitLineData
0795afa3 1\newpage
2\section{Output\label{OUTPUT}}
3
4 The output tape or file contains three types of records. A
5beginning record is written by a call to ISAWBG before generating a set
6of events; an event record is written by a call to ISAWEV for each
7event; and an end record is written for each run by a call to ISAWND.
8These subroutines load the common blocks described below into a single
9\begin{verbatim}
10COMMON/ZEVEL/ZEVEL(1024)
11\end{verbatim}
12and write it out when it is full. A subroutine RDTAPE, described in
13the next section, inverts this process so that the user can analyze
14the event.
15
16 ZEVEL is written out to TAPEj by a call to BUFOUT. For the CDC
17version IF = PAIRPAK is selected; BUFOUT first packs two words from
18ZEVEL into one word in
19\begin{verbatim}
20COMMON/ZVOUT/ZVOUT(512)
21\end{verbatim}
22using subroutine PAIRPAK and then does a buffer out of ZVOUT to TAPEj.
23Typically at least two records are written per event. For all other
24computers IF=STDIO is selected, and ZEVEL is written out with a
25standard FORTRAN unformatted write.
26
27\subsection{Beginning Record}
28
29 At the start of each run ISAWBG is called. It writes out the
30following common blocks:
31\begin{verbatim}
32#include "dylim.inc"
33\end{verbatim}
34\begin{tabular}{lcl}
35QMIN,QMAX &=& $W$ mass limits\\
36QTMIN,QTMAX &=& $W$ $q_t$ limits\\
37YWMIN,YWMAX &=& $W$ $\eta$ rapidity limits\\
38XWMIN,XWMAX &=& $W$ $x_F$ limits\\
39THWMIN,THWMAX &=& $W$ $\theta$ limits\\
40PHWMIN,PHWMAX &=& $W$ $\phi$ limits\\
41\end{tabular}
42
43\begin{verbatim}
44#include "idrun.inc"
45\end{verbatim}
46\begin{tabular}{lcl}
47IDVER &=& program version\\
48IDG(1) &=& run date (10000$\times$month+100$\times$day+year)\\
49IDG(2) &=& run time (10000$\times$hour+100$\times$minute+second)\\
50IEVT &=& event number\\
51\end{tabular}
52
53\begin{verbatim}
54#include "jetlim.inc"
55\end{verbatim}
56\begin{tabular}{lcl}
57PMIN,PMAX &=& jet momentum limits\\
58PTMIN,PTMAX &=& jet $p_t$ limits\\
59YJMIN,YJMAX &=& jet $\eta$ rapidity limits\\
60PHIMIN,PHIMAX &=& jet $\phi$ limits\\
61THMIN,THMAX &=& jet $\theta$ limits\\
62\end{tabular}
63
64\begin{verbatim}
65#include "keys.inc"
66\end{verbatim}
67\begin{tabular}{lcl}
68KEYON &=& normally TRUE, FALSE if no good reaction\\
69KEYS &=& TRUE if reaction I is chosen\\
70 && 1 for TWOJET\\
71 && 2 for E+E-\\
72 && 3 for DRELLYAN\\
73 && 4 for MINBIAS\\
74 && 5 for SUPERSYM\\
75 && 6 for WPAIR\\
76REAC &=& character reaction code\\
77\end{tabular}
78
79\begin{verbatim}
80#include "primar.inc"
81\end{verbatim}
82\begin{tabular}{lcl}
83NJET &=& number of jets per event\\
84SCM &=& square of com energy\\
85HALFE &=& beam energy\\
86ECM &=& com energy\\
87IDIN &=& ident code for initial beams\\
88NEVENT &=& number of events to be generated\\
89NTRIES &=& maximum number of tries for good jet parameters\\
90NSIGMA &=& number of extra events to determine SIGF\\
91\end{tabular}
92
93\begin{verbatim}
94#include "q1q2.inc"
95\end{verbatim}
96\begin{tabular}{lcl}
97GOQ(I,K) &=& TRUE if quark type I allowed for jet k\\
98 && I = 1 2 3 4 5 6 7 8 9 10 11 12 13\\
99 && \ \ $\Rightarrow$ $g$ $u$ $\bar u$ $d$ $\bar d$ $s$
100 $\bar s$ $c$ $\bar c$ $b$ $\bar b$ $t$ $\bar t$\\
101 && I = 14 15 16 17 18 19 20 21 22 23 24 25\\
102 && \ \ $\Rightarrow$ $\nu_e$ $\bar\nu_e$ $e^-$ $e^+$
103 $\nu_\mu$ $\bar\nu_\mu$ $\mu^-$ $\mu^+$ $\nu_\tau$
104 $\bar\nu_\tau$ $\tau^-$ $\tau^+$\\
105GOALL(K) &=& TRUE if all jet types allowed\\
106GODY(I) &=& TRUE if $W$ type I is allowed.\\
107 I= 1 2 3 4\\
108 GM W+ W- Z0\\
109STDDY &=& TRUE if standard DRELLYAN\\
110GOWW(I,K) &=& TRUE if I is allowed in the decay of K for WPAIR.\\
111ALLWW(K) &=& TRUE if all allowed in the decay of K for WPAIR.\\
112\end{tabular}
113
114\begin{verbatim}
115#include "qcdpar.inc"
116\end{verbatim}
117\begin{tabular}{lcl}
118ALAM &=& QCD scale $\Lambda$\\
119ALAM2 &=& QCD scale $\Lambda^2$\\
120CUTJET &=& cutoff for generating secondary partons\\
121ISTRUC &=& 3 for Eichten (EHLQ), \\
122 &=& 4 for Duke (DO) \\
123 &=& 5 for CTEQ 2L\\
124 &=& 6 for CTEQ 3L\\
125 &=& $-999$ for PDFLIB\\
126\end{tabular}
127
128\begin{verbatim}
129#include "qlmass.inc"
130\end{verbatim}
131\begin{tabular}{lcl}
132AMLEP(6:8) &=& $t$,$y$,$x$ masses, only elements written\\
133\end{tabular}
134
135\subsection{Event Record}
136
137 For each event ISAWEV is called. It writes out the following
138common blocks:
139\begin{verbatim}
140#include "final.inc"
141\end{verbatim}
142\begin{tabular}{lcl}
143SIGF &=& integrated cross section, only element written\\
144\end{tabular}
145
146\begin{verbatim}
147#include "idrun.inc"
148\end{verbatim}
149\begin{tabular}{lcl}
150IDVER &=& program version\\
151IDG &=& run identification\\
152IEVT &=& event number\\
153\end{tabular}
154
155\begin{verbatim}
156#include "jetpar.inc"
157\end{verbatim}
158\begin{tabular}{lcl}
159P &=& jet momentum $\vert\vec p\vert$\\
160PT &=& jet $p_t$\\
161YJ &=& jet $\eta$ rapidity\\
162PHI &=& jet $\phi$\\
163XJ &=& jet $x_F$\\
164TH &=& jet $\theta$\\
165CTH &=& jet $\cos(\theta)$\\
166STH &=& jet $\sin(\theta)$\\
167JETTYP &=& jet type. The code is listed under /Q1Q2/ above\\
168 && {\it continued\dots}\\
169\end{tabular}
170
171\begin{tabular}{lcl}
172SHAT,THAT,UHAT &=& hard scattering $\hat s$, $\hat t$, $\hat u$\\
173QSQ &=& effective $Q^2$\\
174X1,X2 &=& initial parton $x_F$\\
175PBEAM &=& remaining beam momentum\\
176QMW &=& $W$ mass\\
177QW &=& $W$ momentum\\
178QTW &=& $W$ transverse momentum\\
179YW &=& $W$ rapidity\\
180XW &=& $W$ $x_F$\\
181THW &=& $W$ $\theta$\\
182QTMW &=& $\sqrt{q_{t,W}^2+Q^2}$\\
183PHIW &=& $W$ $\phi$\\
184SHAT1,THAT1,UHAT1 &=& invariants for $W$ decay\\
185JWTYP &=& $W$ type. The code is listed under /Q1Q2/ above.\\
186ALFQSQ &=& QCD coupling $\alpha_s(Q^2)$\\
187CTHW &=& $W$ $\cos(\theta)$\\
188STHW &=& $W$ $\sin(\theta)$\\
189Q0W &=& $W$ energy\\
190\end{tabular}
191
192\begin{verbatim}
193#include "jetset.inc"
194\end{verbatim}
195\begin{tabular}{lcl}
196NJSET &=& number of partons\\
197PJSET(1,I) &=& $p_x$ of parton I\\
198PJSET(2,I) &=& $p_y$ of parton I\\
199PJSET(3,I) &=& $p_z$ of parton I\\
200PJSET(4,I) &=& $p_0$ of parton I\\
201PJSET(5,I) &=& mass of parton I\\
202JORIG(I) &=& JPACK*JET+K if I is a decay product of K.\\
203 && IF K=0 then I is a primary parton.\\
204 && (JET = 1,2,3 for final jets.)\\
205 && (JET = 11,12 for initial jets.)\\
206JTYPE(I) &=& IDENT code for parton I\\
207JDCAY(I) &=& JPACK*K1+K2 if K1 and K2 are decay products of I.\\
208 && If JDCAY(I)=0 then I is a final parton\\
209MXJSET &=& dimension for /JETSET/ arrays.\\
210JPACK &=& packing integer for /JETSET/ arrays.\\
211\end{tabular}
212
213\begin{verbatim}
214#include "jetsig.inc"
215\end{verbatim}
216\begin{tabular}{lcl}
217SIGMA &=& cross section summed over types\\
218SIGS(I) &=& cross section for reaction I (not written)\\
219NSIGS &=& number of nonzero cross sections (not written)\\
220INOUT(I) &=& packed partons for process I (not written)\\
221MXSIGS &=& dimension for JETSIG arrays (not written)\\
222SIGEVT &=& partial cross section for selected channel\\
223\end{tabular}
224
225\begin{verbatim}
226#include "partcl.inc"
227\end{verbatim}
228\begin{tabular}{lcl}
229NPTCL &=& number of particles\\
230PPTCL(1,I) &=& $p_x$ for particle I\\
231PPTCL(2,I) &=& $p_y$ for particle I\\
232PPTCL(3,I) &=& $p_z$ for particle I\\
233PPTCL(4,I) &=& $p_0$ for particle I\\
234PPTCL(5,I) &=& mass for particle I\\
235IORIG(I) &=& IPACK*JET+K if I is a decay product of K.\\
236 &=& -(IPACK*JET+K) if I is a primary particle from\\
237 && parton K in /JETSET/.\\
238 &=& 0 if I is a primary beam particle.\\
239 && (JET = 1,2,3 for final jets.)\\
240 && (JET = 11,12 for initial jets.)\\
241IDENT(I) &=& IDENT code for particle I\\
242IDCAY(I) &=& IPACK*K1+K2 if decay products are K1-K2 inclusive.\\
243 && If IDCAY(I)=0 then particle I is stable.\\
244MXPTCL &=& dimension for /PARTCL/ arrays.\\
245IPACK &=& packing integer for /PARTCL/ arrays.\\
246\end{tabular}
247
248\begin{verbatim}
249#include "pinits.inc"
250\end{verbatim}
251\begin{tabular}{lcl}
252PINITS(1,I) &=& $p_x$ for initial parton I\\
253PINITS(2,I) &=& $p_y$ for initial parton I\\
254PINITS(3,I) &=& $p_z$ for initial parton I\\
255PINITS(4,I) &=& $p_0$ for initial parton I\\
256PINITS(5,I) &=& mass for initial parton I\\
257IDINIT(I) &=& IDENT for initial parton I\\
258\end{tabular}
259
260\begin{verbatim}
261#include "pjets.inc"
262\end{verbatim}
263\begin{tabular}{lcl}
264PJETS(1,I) &=& $p_x$ for jet I\\
265PJETS(2,I) &=& $p_y$ for jet I\\
266PJETS(3,I) &=& $p_z$ for jet I\\
267PJETS(4,I) &=& $p_0$ for jet I\\
268PJETS(5,I) &=& mass for jet I\\
269IDJETS(I) &=& IDENT code for jet I\\
270QWJET(1) &=& $p_x$ for $W$\\
271QWJET(2) &=& $p_y$ for $W$\\
272QWJET(3) &=& $p_z$ for $W$\\
273QWJET(4) &=& $p_0$ for $W$\\
274QWJET(5) &=& mass for $W$\\
275IDENTW &=& IDENT CODE for $W$\\
276PPAIR(1,I) &=& $p_x$ for WPAIR decay product I\\
277PPAIR(2,I) &=& $p_y$ for WPAIR decay product I\\
278PPAIR(3,I) &=& $p_z$ for WPAIR decay product I\\
279PPAIR(4,I) &=& $p_0$ for WPAIR decay product I\\
280PPAIR(5,I) &=& mass for WPAIR decay product I\\
281IDPAIR(I) &=& IDENT code for WPAIR product I\\
282JPAIR(I) &=& JETTYPE code for WPAIR product I\\
283NPAIR &=& 2 for $W^\pm\gamma$ events, 4 for $WW$ events\\
284\end{tabular}
285
286\begin{verbatim}
287#include "totals.inc"
288\end{verbatim}
289\begin{tabular}{lcl}
290NKINPT &=& number of kinematic points generated.\\
291NWGEN &=& number of W+jet events accepted.\\
292NKEEP &=& number of events kept.\\
293SUMWT &=& sum of weighted cross sections.\\
294WT &=& current weight. (SIGMA$\times$WT = event weight.)\\
295\end{tabular}
296
297\begin{verbatim}
298#include "wsig.inc"
299\end{verbatim}
300\begin{tabular}{lcl}
301SIGLLQ &=& cross section for $W$ decay.\\
302\end{tabular}
303
304 Of course irrelevant common blocks such as /WSIG/ for TWOJET
305events are not written out.
306
307\subsection{End Record}
308
309 At the end of a set ISAWND is called. It writes out the
310following common block:
311\begin{verbatim}
312#include "final.inc"
313\end{verbatim}
314\begin{tabular}{lcl}
315NKINF &=& number of points generated to calculate SIGF\\
316SIGF &=& integrated cross section for this run\\
317ALUM &=& equivalent luminosity for this run\\
318ACCEPT &=& ratio of events kept over events generated\\
319NRECS &=& number of physical records for this run\\
320\end{tabular}
321
322 Events within a given run have uniform weight. Separate runs can
323be combined together using SIGF/NEVENT as the weight per event. This
324gives a true cross section in mb units.
325
326 The user can replace subroutines ISAWBG, ISAWEV, and ISAWND to
327write out the events in a different format or to update histograms
328using HBOOK or any similar package.