If ECAL and HCAL objects are topologically associated to each other, the loops described above give access to composite calorimeter objects rather than to each contributing ECAL and HCAL object separately. It is also possible to get access to all ECAL and HCAL objects, regardless of whether or not they are associated to other reconstructed objects. (The loops described below are equivalent to looping through the PECO and PHCO banks.)
The following statements perform a loop over all ECAL objects; see 7.4. (DO loops cannot be used because the objects are not stored in consecutive locations.)
IOBJ = KPDIR ('ECAL', KRECO)
10 IF (IOBJ .EQ. 0) GO TO 999
C... Analysis of the ECAL object IOBJ ...
IOBJ = KFOLLO (IOBJ)
GO TO 10
(The functions KPDIR and KFOLLO are described in
7.4.)
The corresponding loop for HCAL objects is:
IOBJ = KPDIR ('HCAL', KRECO)
10 etc ...