Event topology routines and ENFLW

To use the event topology routines described in Chapter 10 with these energy-flow objects, use option 'EF' with subroutine QJOPTR (see 10.1):

      CALL QJOPTR('EF',' ')

Example:

The following code calculates the total energy energy of an event and finds the thrust using energy flow objects.

      E=0.
      DO 10 I = KFEFT, KLEFT
         E=E + QE(I)
   10 CONTINUE
C---    Find thrust
      CALL QJOPTR('EF', ' ')
      CALL QJTHRU(THRU, 'THRU', KRECO)

Jets based on energy flow objects using QJMMCL with YCUT = 0.003 (see Sec. 10.11.1) are stored in the EJET bank. If the EFLJ card is used instead of the EFLW card, the EFT section will be filled as described above, and these jets will be stored in the JET section. The jets may be accessed with DO loops (KFJET, KLJET, KNJET) or with the particle name 'EJET' using the functions KPDIR and KFOLLO. The energy flow objects making up these jets can be found with XSAME as described in Sec. 8.1.7. To save time, these jets may be used as input for jet-finding with a higher YCUT (see 10.11.1) by calling QJOPTR with the option EJ:

     CALL QJOPTR('EJ',' ').

XSAME may be used to find the original energy flow objects (in the EFT section) making up the final jets.