*-- Author : Stephen Burke SUBROUTINE FFRAD(JDIG,NPLAN) *-----------------------------------------Updates 27/07/93------- **: FFRAD 30907 RP. Farm changes. *-----------------------------------------Updates 03/05/93------- **: FFRAD 30907 SB. Radius ignored if there are planar hits. *-----------------------------------------Updates 06/08/92------- **: FFRAD 30907 SB. Fix bad radius bug. *-----------------------------------------Updates 13/02/92------- **: FFRAD 30205.SB. ERRLOG error numbers changed. *-----------------------------------------Updates 07/02/92------- **: FFRAD 30205.SB. Remove unused FKMEAS sequence. *-----------------------------------------Updates 24/01/92------- **: FFRAD 30205.SB. ERRLOG message format changed. *-----------------------------------------Updates---------------- ********************************************************************** * * * Fill the Kalman filter arrays from the radial DIGI bank * * * ********************************************************************** *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,FFSTEE. *KEEP,FFGEO. *KEEP,FKFLAG. *KEEP,FKMEAS. *KEEP,FKCONS. *KEEP,FFWBI. * Work bank indices (note that INDKTR is *NOT* a work bank index!) *KEEP,BCS. *KEEP,CNSTBF. *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 *KEEP,STFCLW. * statement functions acting on the BOS COMMON LW * index of element before row number LWROW * index of L-th element of row number LWROW * L-th integer element of the LWROW'th row in bank with index LND * L-th real element of the LWROW'th row in bank with index LND * *KEEP,FTFUNCT. * Statement functions for RADIAL Chamber data access. * Using Channel Number J * Module, Wedge-pair and Z-plane numbers... * Statement function for obtaining WEDGE numbers(0-47) of * wires at plus and minus ends of Cell numbers * Statement function for obtaining IOS wire number (1-36) * Statement functions for PLANAR Chamber data access. * Using Channel Number J * Module, orientation, W-cell and Z-plane numbers... * IPSMD in range 0:8 Planar module number. * * IOS wire number (runs from 0 to 36) * SB plane numbers (1-72) from cell number * Module, orientation, wire and (typical) cell number from plane * number in the range 1-72 (planars, radials and combined) *KEND. ********************************************************************** CALL ERRLOG(361,'F:FFRAD: Radial digi on mapped-out plane') * Remember the digi pointer and drift sign * DD = RBTAB(INDLC(2),2,JDIG) * Use assymetric drift * DD = -RBTAB(INDLC(2),2,JDIG) * Use assymetric drift * Measured position and error * WMES(1,JPL) = DD + SBTAB(INDG1(2),3+3*KWED,ICELL+1) * Use geometric stagger (3 & 6 -> 8 & 9) WMES(2,JPL) = RR + FLOREN(RR,ABS(DD),DSIGN) *