First commit.
[u/mrichter/AliRoot.git] / ISAJET / doc / tape.doc
1 \newpage
2 \section{File Reading\label{TAPE}}
3
4       The FORTRAN instruction
5 \begin{verbatim}
6       CALL RDTAPE(IDEV,IFL)
7 \end{verbatim}
8 will read a beginning record, an end record or an event (which can be
9 more than one record). IDEV is the tape number and
10 \begin{verbatim}
11       IFL=0  for a good read,
12       IFL=-1 for an end of file.
13 \end{verbatim}
14 The information is restored to the common blocks described above. The
15 type of record is contained in
16 \begin{verbatim}
17 #include "rectp.inc"
18 \end{verbatim}
19 \begin{tabular}{lcl}
20 IRECTP            &=& 100 for an event record\\
21 IRECTP            &=& 200 for a beginning record\\
22 IRECTP            &=& 300 for an end record\\
23 IREC              &=& no. of physical records in event record, 0 
24                       otherwise\\
25 \end{tabular}
26
27       The parton momenta from the primary hard scattering are
28 contained in /PJETS/. The parton momenta generated by the QCD cascade
29 are contained in /JETSET/. The hadron momenta both from the QCD jets
30 and from the beam jets are contained in /PARTCL/. The final hadron
31 momenta and the associated pointers should be used to calculate the
32 jet momenta, since they are changed both by the QCD cascade and by
33 hadronization. Particles with IDCAY=0 are stable, while the others are
34 resonances.
35
36       The weight per event needed to produce a weighted histogram in
37 millibarn units is SIGF/NEVENT. The integrated cross section SIGF is
38 calculated by Monte Carlo integration during the run for the given
39 kinematic limits and JETTYPE, WTYPE, and WMODE selections. Any of three
40 methods can be used to find the value of SIGF:
41
42       (1) The current value, which is written out with each event, can
43 be used. To prevent enormous fluctuations at the beginning of a run,
44 NSIGMA extra primary parton events are generated first. The default
45 value, NSIGMA = 20, gives negligible overhead but may not be large
46 enough for good accuracy.
47
48       (2) The value SIGF calculated with the full statistics of the run
49 can be obtained by reading through the tape until an end record
50 (IRECTP=300) is found. After SIGF is saved with a different name, the
51 first event record for the run can be found by backspacing the tape
52 NRECS times.
53
54       (3) Unweighted histograms can be made for the run and the weight
55 added after the end record is found. An implementation of this using
56 special features of HBOOK is contained in ISAPLT.
57
58       The functions AMASS(IDENT), CHARGE(IDENT), and LABEL(IDENT) are
59 available to determine the mass, charge, and character label in A8
60 format. Subroutine FLAVOR returns the quark content of any hadron and
61 may be useful to convert IDENT codes to other schemes. CALL PRTEVT(0)
62 prints an event.