*-- Author : S.Burke SUBROUTINE FFHUNT(JDIGP,JDIGR,JMAX) ********************************************************************** * * * Hunt for the best candidate for the 'true' track, given a list * * of digis. * * * * JDIGP and JDIGR are pointers to the first digis on the track. * * JMAX is returned as the STR index of the 'true' track. * * * * Note that the FRPX, FRRX, FPUX and FRUX banks are assumed to have * * been copied to work banks, with pointers in /FFWBI/. * * * ********************************************************************** *KEEP,FKNPL. * * Per-track values can go in H1WORK; note that LTRUE and LFIRST must * be set at least per event. * * This is about 36k words long; the remaining common blocks are * about 3.6k in total. Some of this could be in /H1WORK/, but the * blocks would have to be reorganised. * * /FKPROJ/ * /FKFILT/ * /FKSMTH/ * /FKINT/ * /FKRSID/ * /FKTRUE/ * /FKDBG/ *KEEP,FFWBI. * Work bank indices (note that INDKTR is *NOT* a work bank index!) *KEEP,BCS. *KEEP,STFUNCT. * index of element before row number IROW * index of L'th element of row number IROW * L'th integer element of the IROW'th row of bank with index IND * L'th real element of the IROW'th row of bank with index IND *KEND. ********************************************************************** * Allow for hit banks to be nonexistant CALL VZERO(ITRUE,2000) * Loop over planar digis looking for a corresponding hit ... IF (IDIG.EQ.JDIG) CALL FFNUM(1,JHIT,ITRUE) * ... and the same for radials IF (IDIG.EQ.JDIG) CALL FFNUM(2,JHIT,ITRUE) * ITRUE now contains the number of hits for each track ... *