ISAJET/doc/output.doc

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