ISAJET/doc/ztext.doc

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1996/03/08 17:27:17  mclareni
   6 * Isajet714
   7 *
   8 *
   9
  10                           ISAJET ZEBRA INTERFACE
  11                      F.E. Paige and S.D. Protopopescu
  12                             Physics Department
  13                       Brookhaven National Laboratory
  14                              Upton, NY 11973
  15                               March 30, 1991
  16
  17
  18      [1] Overview
  19      ============
  20
  21       The ZEBRA  version of ISAJET  contains subroutines  which replace the
  22 usual  output package  for ISAJET  with one using  the CERN  package ZEBRA.
  23 This version  provides some  features which  are not part  of the  standard
  24 ISAJET  output. Momenta  and masses  of the  partons are  rescaled from the
  25 final  particles so that  they match  exactly  (this is not  so in ISAJET).
  26 Resonances   decaying  strongly are  not part  of the  particle  list while
  27 particles with heavy quarks decaying  weakly are part of a vertex list with
  28 the corresponding  vertex position  and appropriate  links to the particles
  29 belonging  to that  vertex.  Thus there  is no need  to  generate secondary
  30 vertices for charmed or other heavy quark particles and the event structure
  31 makes it easy to find which particles come from the decay of a heavy quark.
  32 The  association of  particles  with  parent partons  is also  made easy by
  33 reference  links. For simple  calorimeter  studies banks  are also provided
  34 with  energy sums  over  calorimeter  cells and  banks  containing only the
  35 leptons (see  below). If only those  are written  instead of particle banks
  36 the resulting files are more compact and require considerably less CPU time
  37 to analyze (albeit the information is rather limited). In addition jets are
  38 calculated by  using a cone  algorithm on the partons,  these can be useful
  39 for studying  efficiencies of jet finding  algorithms and comparing results
  40 to theoretical jets. The results of parton jets are in Zebra banks PJET.
  41
  42       NOTE: The present  version can only handle  one run at a time, so you
  43 cannot supply a command  file with many runs for  now. This problem will be
  44 fixed in a future release.
  45
  46     [2] Creating and running a job
  47     ==============================
  48
  49     There are command files in BNLCL6::$2$DUA14:[ISAJET.ISALIBRARY] to make
  50 EXE files for event generation and analysis.
  51
  52     To make an EXE generating ISAJET events which are then written to
  53 a file in ZEBRA format users should run command file:
  54
  55 $ @isa$area:MAKE_ISAJET_ZEBRA
  56
  57 this starts from  ISAJET.PAM and  makes  ISAJET_ZEBRA.FOR, ISAJET_ZEBRA.OLB
  58 and  ISAJET_ZEBRA.EXE.  It requires  that  ISA$AREA de  defined as the area
  59 where ISAJET.PAM and command files  reside. It also requires that CERN$AREA
  60 be defined   as the area  for the  CERN  libraries. The above  command file
  61 assumes symbol YPATCHY is defined. If running on a machine other than a VAX
  62 looking at the command file should give a pretty good idea of what needs to
  63 be done.
  64
  65 If you  already have ISAJET_ZEBRA.OLB you can run command file
  66
  67 $ @isa$area:LINK_ISAJET_ZEBRA
  68
  69 to make ISAJET_ZEBRA.EXE, this program runs with an interactive
  70 dialog to generate ISAJET instructions. You can supply your own
  71 ISARUN main program if you wish to bypass the dialog, see command
  72 file MAKE_ISAJET.COM for the simplest example.
  73
  74 User can provide any of 3 logical functions to reduce the
  75 size of output files:
  76 1) EDIT(I)     = if false event will not be written out
  77                  I is event number
  78                  Zebra banks have not yet been generated, so
  79                  user must work from standard ISAJET common blocks
  80 2) ISEDIT()    = same as EDIT but Zebra banks have been filled
  81 3) PEDIT(ID,P) = if false particle will not be included in
  82                  the list of particles written out
  83                  ID= particle ID, P(4)= particle momentum vector
  84
  85
  86 To generate a user program that reads an ISAJET output file in ZEBRA
  87 format users should run
  88
  89 $ @isa$area:MAKE_ISZRUN
  90 this makes a sample user analysis subroutine ISZUSR and an ISZRUN.EXE
  91 program. ISZUSR gives an example of how to use the utility subroutines
  92 mentioned in next section for analysis.
  93
  94 If user already has an ISZUSR subroutine he can generate ISZRUN.EXE by
  95
  96 $ @isa$area:LINK_ISZRUN
  97
  98
  99     [3] Utilities
 100     =============
 101
 102     A number of utilities  are provided for the  user to make it easier for
 103 him to use the  information in ZEBRA form without  him having to learn much
 104 about  ZEBRA. It is   nevertheless  recommended that  the user  familiarize
 105 himself with ZEBRA by  reading the first chapter  of the ZEBRA manual. Also
 106 to make  full  use of  the power  of  ZEBRA users  may   eventually want to
 107 familiarize themselves with some of  the utilities provided with it such as
 108 MZDROP and ZSORT.
 109
 110 o For every ZEBRA bank of name ISxx (most ISAJET ZEBRA banks start with
 111   IS) there exist:
 112     1) GZISxx : an integer function returning pointer to first bank in
 113                 the linear structure
 114     2) GTISxx : a subroutine returning the contents of the bank
 115     3) PRISxx : a subroutine generating a printout of the banks with name ISxx
 116
 117 o To generate a full dump:
 118         CALL PRTEVZ(unit #)
 119   This subroutine calls all the PRISxx subroutines.
 120
 121 o Users can provide a primary vertex for each event by:
 122         CALL SETVTX(XYZ)
 123   where XYZ is a 3-vector with values x,y,z
 124   To supply a function for the vertex user should supply subroutine
 125   INTVTX(XYZ) which returns an XYZ vector with each call. The
 126   default returns 0,0,0.
 127
 128 o The parameters for generating PJET banks can be controlled by
 129       CALL PJPSET(N,ALG,ETCUT,DRCUT,MAXIT,IR,MUON)
 130   where
 131                                                           Default values
 132           N     = NUMBER OF PJET ALGORITHMS                      1
 133           ALG   = ALGORITHM TPYE 'CONE' OR 'ANGLE'            'CONE'
 134           ETCUT = ET CUT FR PJET JETS                           5.0
 135           DRCUT = CONE CUT (IN R OR ANGLE)                      0.45
 136           MAXIT = MAXIMUM NUMBER OF ITERATIONS                   1
 137           IR    = INIT RADIATION FLAG 0=USE IR,   1=NO IR        0
 138           MUON  = MUON FLAG    0=NO MUONS, 1=use MUONS in PJET   0
 139
 140
 141       Other subroutines of interest are:
 142 o ISZJET  find jets using ISCL banks
 143 o ISZUSR  example of a user subroutine for analysis
 144
 145
 146
 147     [4] Other Subroutines
 148     =====================
 149
 150     We list here the subroutines that are part of ISAJET_ZEBRA:
 151
 152 GENVTX  generate secondary vertices (for heavy quark decays)
 153 ISABFL  fill ISAB bank (begin of run bank)
 154 ISACFL  fill ISAC and ISCL banks (calorimeter description and
 155         calorimeter cells with non-zero energy deposited).
 156         Entry point ISACIN calculates all numbers needed for
 157         calorimeter description.
 158 ISAEFL  fill ISAE bank (main header)
 159 ISAFFL  fill ISAF bank (end of run bank)
 160 ISAJFL  fill ISAJ banks (primary jets)
 161 ISALFL  fill ISAL banks (leptons)
 162 ISAPFL  fill ISV1 and ISP1 banks (vertices and particles)
 163 ISAQFL  fill ISAQ banks (secondary partons, both initial and final)
 164 ISCMFL  fill ISCM (copy of command instructions)
 165 ISMEAR  generates smeared calorimeter data (modifies ISCL banks)
 166 PJETFL  fill PJET banks (jets obtained from partons)
 167 QRECAL  to recalculate momenta and masses of partons starting
 168         from the final particles
 169
 170
 171
 172
 173
 174    [5] Brief description of Zebra banks
 175    ====================================
 176
 177       The beginning and end of run  records have Zebra banks named ISAB and
 178 ISAF  while  the  whole  event  hangs  from  Zebra  bank  ISAE. A  complete
 179 description of the banks is given in section [6].
 180
 181 ISAB data describing event generation              (in beginning record)
 182 ISCM copy of command file used for event generation(       "           )
 183 ISAF cumulative results from run (in end record)
 184 ISAE general event information such as weight, type, etc.
 185 ISAJ primary jet momenta and masses
 186 ISAQ final and initial partons (momenta and masses)
 187 ISV1 primary and short lived vertices (heavy quarks), includes
 188      id,mass and 4-momenta of decaying particle
 189 ISP1 particles associated with ISV1 (connected via a structural link)
 190      id,mass and 4-momenta
 191 ISAC simple calorimeter description (no depth)
 192 ISCL energy deposition and cell description in calorimeter cell
 193 ISMR description of smearing parameters
 194 ISAL stable leptons (including neutrinos):id,mass and 4-momenta
 195      with reference links to calorimeter cell and ISP1 if they
 196      exist
 197 ISJT jets found by subroutine ISZJET.
 198 PJHD header for PJET banks
 199 PJET jets found adding over partons
 200 PJPT pointers to parton banks (ISAQ) contributing to a jet in PJET
 201
 202       Whether   banks    ISV1,ISP1,ISAC,ISCL  and  ISAL are  written out is
 203 optional. The ISARUN job  asks user to chose an  option ISAP, ISAC or ISAL.
 204 If ISAP is selected then  ISV1 and ISP1 are  written out, if ISAC then ISAC
 205 and  ISCL are   written out  and  if ISAL  then  ISAL  is  written out. Any
 206 combination can be selected, the  instruction ISAPISACISAL causes all banks
 207 to be written out.
 208
 209       The banks ISJT are created by calling ISZJET, the banks ISAC and ISCL
 210 must be available. One can create  the banks ISAC,ISCL and/or ISAL starting
 211 from a file  containing ISV1 and  ISP1 by calling  ISACFL and/or ISALFL for
 212 each event.
 213
 214
 215
 216    [6] Zebra banks documentation
 217    =============================
 218
 219
 220                      +------------------+
 221                      |      ISAE        |
 222                      | event descriptor |
 223                      |                  |
 224                      +------------------+
 225         struc. link         |||   |  |
 226                             |||   |  +------------------------------+
 227                             |||   +-----------------------------+   |
 228        +-------+------------+|                                  |   |
 229        |       |             |                                  |   |
 230   ______     ______       ______    ______                      |   |
 231  | ISAJ \   | ISAQ \     | ISV1 |  | ISV1 \                     |   |
 232  |  -1   >  |  -2   >    |  -3  |--| (-3)  >                    |   |
 233  |______/   |______/     |______|  |______/                     |   |
 234               :  ^         | :...<..: |                         |   |
 235               :  :         |   [-2]   |                         |   |
 236   [PJET -1].<..  :         |          |                         |   |
 237                  : [-2] ______      ______                      |   |
 238                  :.....| ISP1 \    | ISP1 \                     |   |
 239   [PJET -3] ..<........|  -1   >   |  -1   >                    |   |
 240                        |______/    |______/                     |   |
 241                         : |[-1]      : : :                      |   |
 242                     [-4]: v     <....: : :[-4]                  |   |
 243                         : free  <......: :                      |   |
 244                         v                v                      |   |
 245                       ISV2             ISV2                     |   |
 246                    (end vertex)     (end vertex)                |   |
 247                                                                 |   |
 248                              +------------+---------------------+   |
 249                              |            |                         |
 250                              |            |                         |
 251                          _________      ______                     ______
 252                         |  ISAC   |    | ISAL \                   | PJHD \
 253                         |   -6    |    |  -7   >                  |  -8   >
 254                         |_________|    |______/                   |______/
 255                               |                                     |
 256             +-------+---------+--+                                 ______
 257             |       |            |               [ISAQ -1]..>.....| PJET \
 258          ______     ______       ______                           |  -1   >
 259         | ISCL \   | ISMR \     | ISJT \                          |______/
 260         |  -1   >  |  -2   >    |  -3   >                             |
 261         |______/   |______/     |______/                           ______
 262                                                                   | PJPT |
 263                                                      [ISAQ]...<...|  -1  |
 264                                                                   |______|
 265
 266      +-----------+            +-----------+
 267      |    ISAB   |            |    ISAF   |
 268      | begin run |            |  end  run |
 269      |           |            |           |
 270      +-----------+            +-----------+
 271            |
 272      +----------+
 273      |   ISCM   |
 274      |    -1    |
 275      +----------+
 276
 277
 278
 279 C=======================================================================
 280 C
 281 C  Bank name: ISAC
 282 C  Author   : SDP
 283 C  Date     : June 14,1989
 284 C
 285 C                    Calorimeter descriptor
 286 C
 287 C    LQ     Q/IQ
 288 C ___________________________________________________________________________
 289 C     -3            struct.  to ISJT (jets found using ISZJET)
 290 C     -2            struct.  to smearing parameters (ISMR)
 291 C     -1            struct.  to calorimeter cells (ISCL)
 292 C      0            next     to next ISAC (only one at present)
 293 C     +1            up       to ISAE
 294 C     +2            origin   to ISAE for first and previous ISAC for others
 295 C ............................................................................
 296 C            I-5             calorimeter number
 297 C             -4             bank name,  'ISAC'
 298 C             -3             NL=2
 299 C             -2             NS=0
 300 C             -1             ND=6
 301 C              0             STATUS
 302 C             +1     I       number of phi cells
 303 C              2     I         "       eta     "
 304 C              3     F      phi cell size
 305 C              4     F      eta     "   "
 306 C              5     F      minimum eta
 307 C              6     F      maximum eta
 308 C=======================================================================
 309
 310
 311
 312 C=======================================================================
 313 C
 314 C  Bank name: ISAE
 315 C  Author   : SDP
 316 C  Date     : June 19,1986     DH add ISAM
 317 C
 318 C                Top level bank and event descriptor for ISAZEB
 319 C
 320 C    LQ     Q/IQ
 321 C ___________________________________________________________________________
 322 C     -8            struct.  to PJHD (parton jet structure header)
 323 C     -7            struct.  to lepton bank (ISAL)
 324 C     -6            struct.  to pseudo cal. (ISAC)
 325 C     -5            struct.  to unassociated particles (ISP3)
 326 C     -4            struct.  to long lived vertices (ISV2)
 327 C     -3            struct.  to short  "      "     (ISV1)
 328 C     -2            struct.  to final/initial bank (ISAQ)
 329 C     -1            struct.  to primary jet (ISAJ)
 330 C      0            next     to next ISAE
 331 C     +1            up
 332 C     +2            origin
 333 C ............................................................................
 334 C            I-5             bank number
 335 C             -4             bank name,  'ISAE'
 336 C             -3             NL=7
 337 C             -2             NS=7
 338 C             -1             ND=18
 339 C              0             STATUS
 340 C             +1     I       event id(1)
 341 C              2     I       event id(2)
 342 C              3     I       event number
 343 C              4     I       reaction type
 344 C              5     I       number of primary jet banks
 345 C              6     I         "    of stable parton banks (final+initial)
 346 C              7     I         "    of PJET banks
 347 C              8     I         "    of particle banks
 348 C              9     I         "    of vertex banks
 349 C             10     I         "    of lepton banks
 350 C             11     F       cross section in microbarns
 351 C             12     F       weight
 352 C             13     F       effective q**2
 353 C             14     F       hard scattering invariant s
 354 C             15     F        "      "          "      t
 355 C             16     F        "      "          "      u
 356 C             17     D        Seed (part 1)
 357 C             18     D        Seed (part 2)
 358 C=======================================================================
 359
 360
 361
 362 C=======================================================================
 363 C
 364 C  Bank name: ISAJ
 365 C  Author   : SDP
 366 C  Date     : June 19,1986
 367 C
 368 C                Primary jet information
 369 C
 370 C    LQ     Q/IQ
 371 C ___________________________________________________________________________
 372 C    L-1            ref. to parent primary (ISAJ) (0 if original)
 373 C     +0            next     to next jet (ISAJ)
 374 C     +1            up       to ISAE
 375 C     +2            origin   to ISAE
 376 C ............................................................................
 377 C            I-5             Primary jet number
 378 C             -4             bank name,  'ISAJ'
 379 C             -3             NL=0
 380 C             -2             NS=0
 381 C             -1             ND=9
 382 C              0             STATUS
 383 C             +1     I       jet ID
 384 C              2     F       px
 385 C              3     F       py
 386 C              4     F       pz
 387 C              5     F       p
 388 C              6     F       mass
 389 C              7     F       phi
 390 C              8     F       theta
 391 C              9     F       eta
 392 C=======================================================================
 393
 394
 395
 396 C=======================================================================
 397 C
 398 C  Bank name: ISAL
 399 C  Author   : SDP
 400 C  Date     : June 19,1986
 401 C
 402 C                  Lepton bank (from primary and short lived vertices)
 403 C
 404 C    LQ     Q/IQ
 405 C ___________________________________________________________________________
 406 C    L-4            ref.  to particle bank  (ISP1)
 407 C     -3            ref.  to primary jet    (ISAJ)
 408 C     -2            ref.  to initial parton (ISAQ)
 409 C     -1            struct.  (free)
 410 C      0            next     to next ISAL
 411 C     +1            up       to ISAE
 412 C     +2            origin   to ISAE for first and previous ISAL for others
 413 C ............................................................................
 414 C            I-5             particle number
 415 C             -4             bank name,  'ISAL'
 416 C             -3             NL=4
 417 C             -2             NS=1
 418 C             -1             ND=9
 419 C              0             STATUS
 420 C             +1     I       ISAJET particle ID
 421 C              2     F       px
 422 C              3     F       py
 423 C              4     F       pz
 424 C              5     F       p
 425 C              6     F       mass
 426 C              7     F       phi
 427 C              8     F       theta
 428 C              9     F       eta
 429 C=======================================================================
 430
 431
 432
 433 C=======================================================================
 434 C
 435 C  Bank name: ISAQ
 436 C  Author   : SDP
 437 C  Date     : June 19,1986
 438 C
 439 C            Final and initial partons
 440 C
 441 C    LQ     Q/IQ
 442 C ___________________________________________________________________________
 443 C    L-2            ref.  to PJET (parton structure jets)
 444 C     -1            ref.  to primary jet (ISAJ)
 445 C                            0 for initial parton
 446 C      0            next     to next ISAQ
 447 C     +1            up       to ISAE
 448 C     +2            origin   to ISAE for first and previous ISAQ for others
 449 C ............................................................................
 450 C            I-5             final/initial parton number
 451 C             -4             bank name,  'ISAQ'
 452 C             -3             NL=2
 453 C             -2             NS=0
 454 C             -1             ND=9
 455 C              0             STATUS
 456 C             +1     I       final/initial parton type
 457 C              2     F       px
 458 C              3     F       py
 459 C              4     F       pz
 460 C              5     F       E
 461 C              6     F       mass
 462 C              7     F       phi
 463 C              8     F       theta
 464 C              9     F       eta
 465 C=======================================================================
 466
 467
 468
 469 C=======================================================================
 470 C
 471 C  Bank name: ISCL
 472 C  Author   : SDP
 473 C  Date     : June 19,1986
 474 C
 475 C               Calorimeter cell bank
 476 C
 477 C    LQ     Q/IQ
 478 C ___________________________________________________________________________
 479 C    L-2            ref.  to reconstructed jet (ISJT)
 480 C     -1            ref.  to parent parton (ISAQ)
 481 C      0            next     to next ISCL
 482 C     +1            up       to ISAC
 483 C     +2            origin   to ISAC for first and previous ISAL for others
 484 C ............................................................................
 485 C            I-5             cell number (ordinal)
 486 C             -4             bank name,  'ISCL'
 487 C             -3             NL=2
 488 C             -2             NS=0
 489 C             -1             ND=10
 490 C              0             STATUS
 491 C             +1     I       10000*(phi index) + y index
 492 C              2     I       10000*(# of charged trks) + # of gammas
 493 C              3     F       e.m. energy deposited
 494 C              4     F       had.  "        "
 495 C              5     F       sin(theta) at center of cell
 496 C              6     F       cos(theta)         "
 497 C              7     F       sin(phi)           "
 498 C              8     F       cos(phi)           "
 499 C              9     F       phi                "
 500 C             10     F       eta                "
 501 C===========================================================================
 502
 503
 504
 505
 506 C=======================================================================
 507 C
 508 C  Bank name: ISJT
 509 C  Author   : SDP
 510 C  Date     : June 19,1986
 511 C
 512 C                Jets found by ISZJET
 513 C
 514 C    LQ     Q/IQ
 515 C ___________________________________________________________________________
 516 C     -1            to possible parent parton (ISAQ)
 517 C     +0            next     to next jet (ISJT)
 518 C     +1            up       to ISJT
 519 C     +2            origin   to ISAC
 520 C ............................................................................
 521 C            I-5             jet number
 522 C             -4             bank name,  'ISJT'
 523 C             -3             NL=0
 524 C             -2             NS=0
 525 C             -1             ND=6
 526 C              0             STATUS
 527 C             +1     F       ET (transverse energy)
 528 C              2     F       px
 529 C              3     F       py
 530 C              4     F       pz
 531 C              5     F       E
 532 C              6     F       mass
 533 C              7     F       phi
 534 C              8     F       theta
 535 C              9     F       eta
 536 C=======================================================================
 537
 538
 539
 540
 541 C=======================================================================
 542 C
 543 C  Bank name: ISP1
 544 C  Author   : SDP
 545 C  Date     : June 19,1986
 546 C
 547 C             Particle bank (from primary and short lived vertices)
 548 C
 549 C    LQ     Q/IQ
 550 C ___________________________________________________________________________
 551 C    L-5            ref.  to GCAH
 552 C    L-4            ref.  to secondary vertex (ISV2)
 553 C     -3            ref.  to primary jet    (ISAJ)
 554 C     -2            ref.  to initial parton (ISAQ)
 555 C     -1            struct.  (free)
 556 C      0            next     to next ISP1
 557 C     +1            up       to ISV1
 558 C     +2            origin   to ISV1 for first and previous ISP1 for others
 559 C ............................................................................
 560 C            I-5             particle number
 561 C             -4             bank name,  'ISP1'
 562 C             -3             NL=4
 563 C             -2             NS=1
 564 C             -1             ND=9
 565 C              0             STATUS
 566 C             +1     I       ISAJET particle ID
 567 C              2     F       px
 568 C              3     F       py
 569 C              4     F       pz
 570 C              5     F       E
 571 C              6     F       mass
 572 C              7     F       phi
 573 C              8     F       theta
 574 C              9     F       eta
 575 C=======================================================================
 576
 577
 578
 579 C=======================================================================
 580 C
 581 C  Bank Name : PJET
 582 C  Author    : Chip Stewart
 583 C  Date      :  7-NOV-1989 18:10:09.84
 584 C  Version   : 2.0
 585 C
 586 C  Bank description : PARTON JET BANK
 587 C      PJET is constructed by applying a jet algorithm (CONE or
 588 C      OPENING ANGLE) to either Partons given in the ISAQ bank
 589 C      or Paricles given in the ISP1 bank.
 590 C
 591 C     LQ     Q/IQ
 592 C-----------------------------------------------------------------------
 593 C     +2          Down   link to SPARE
 594 C     -1          Down   link to PJPT
 595 C      0          Next   link to
 596 C     +1          Up     link to PJHD
 597 C     +2          Origin link to PJHD
 598 C.......................................................................
 599 C             -5         Bank number
 600 C             -4         Bank name, 'PJET'
 601 C             -3         NL = 2
 602 C             -2         NS = 2
 603 C             -1         ND = 12
 604 C              0         Status
 605 C             +1     I   NV = 2 Version Number
 606 C              2     F   Et ( =Pt in version NV=1 )
 607 C              3     F   Px
 608 C              4     F   Py
 609 C              5     F   Pz
 610 C              6     F   E
 611 C              7     F   mass
 612 C              8     F   phi
 613 C              9     F   theta
 614 C             10     F   eta
 615 C             11     I   CLASS - (eg. for mapping PJET with RECO )
 616 C             12     I   NEXT
 617 C=======================================================================
 618 C
 619 C  Definitions:
 620 C  -----------
 621 C
 622 C  E  - Sum(Ei) over all the Partons/Particles included in the jet.
 623 C
 624 C  Px - Sum(Pxi).
 625 C  Py - Sum(Pyi).
 626 C  Pz - Sum(Pzi).
 627 C
 628 C  Et - Sum[Ei*sin(Thetai)]
 629 C
 630 C  Comment: In version NV=1 word 2 was Pt = SQRT ( Px**2 + Py**2 )
 631 C
 632 C  Mass - SQRT [ E**2 - ( Px**2 + Py**2 + Pz**2 )]
 633 C
 634 C  For Theta, Phi and Eta the following code was used
 635 C  to calculate them from Px,Py,Pz
 636 C
 637 C  PARAMETER( SMALL = 1.0E-5 )
 638 C  Phi=ATAN2(Py,Px+SMALL)   ;   IF(Phi.LT.0) Phi=Phi+TWOPI
 639 C  Theta=ACOS((Pz+SMALL)/(SQRT(Px*Px+Py*Py+Pz*Pz)+SMALL))
 640 C  Eta=-ALOG(TAN(Theta/2.)+SMALL)
 641 C
 642 C=======================================================================
 643
 644
 645
 646 C=======================================================================
 647 C
 648 C  Bank Name : PJHD
 649 C  Author    : Chip Stewart
 650 C  Date      :  7-NOV-1989 17:57:58.00
 651 C
 652 C  Bank description : Header for PJET structure
 653 C
 654 C     LQ     Q/IQ
 655 C-----------------------------------------------------------------------
 656 C     -2          DOWN   spare
 657 C     -1          DOWN   LINK TO PJET
 658 C      0          Next   link to PJHD
 659 C     +1          Up     link to ISAE
 660 C     +2          Origin link to ISAE
 661 C.......................................................................
 662 C             -5         Bank number
 663 C             -4         Bank name, 'PJHD'
 664 C             -3         NL = 2
 665 C             -2         NS = 2
 666 C             -1         ND = 8
 667 C              0         Status
 668 C              1     I   NV=2 Version number
 669 C              2     I   N_ALG - 1 for DR CONE , 2 for CMS OPEN ANGLE
 670 C              3     I   NJET - NUMBER OF PARTON JETS
 671 C              4     F   algorithm parameters DR_CONE_CUT/OPENING_ANGLE_CUT
 672 C              5     F   algorithm parameters JET_ET_CUT
 673 C              6     I   algorithm parameters MAX ITERATIONS
 674 C              7     I   Initial Radiation switch 0-use IR      1-no IR
 675 C              8     I   MUON switch (IF NV=1 MUON NOT USED)
 676 C                        0-no MUONS    1-use MUONS
 677 C=======================================================================
 678
 679
 680
 681 C=======================================================================
 682 C
 683 C  Bank Name : PJPT
 684 C  Author    : Chip Stewart
 685 C  Date      : 15-Dec-1989
 686 C
 687 C  Bank description : pointers to Partons (ISAQ) in a given jet
 688 C
 689 C     LQ     Q/IQ
 690 C-----------------------------------------------------------------------
 691 C     -N+1        Ref Link to parton bank ISAQ for Nth parton in this PJET
 692 C     -2          Ref Link to parton bank ISAQ for 1st    "    "  "    "
 693 C     -1          not used
 694 C      0          Next   link to none
 695 C     +1          Up     link to PJET
 696 C     +2          Origin link to PJET
 697 C.......................................................................
 698 C             -5         Bank number
 699 C             -4         Bank name, 'PJPT'
 700 C             -3         NL = N+1
 701 C             -2         NS = 0
 702 C             -1         ND = 1 - number of data words in PJPT bank = 1
 703 C              0         Status
 704 C             +1     I   Bank version (=1)
 705 C=======================================================================
 706