*-- Author : Stephen Burke SUBROUTINE FFEVT0(DEVT0,ZNOM,VXYZ) ********************************************************************** * * * Calculate the event T0 correction (no e-bunch timing correction * * at the moment) * * * * OUTPUT: * * DEVT0 is the event T0 correction * * ZNOM is the nominal z vertex * * VXYZ is the event vertex position * * * ********************************************************************** * c = 30 cm/nsec *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,FKPROJ. *KEEP,FKSMTH. *KEEP,FKRSID. *KEEP,FKTRUE. *KEEP,FFSTEE. *KEEP,FFDBG. *KEEP,FKDBG. *KEEP,FKINT. *KEEP,FRLORA. *KEEP,BCS. *KEEP,BOSMDL. C ------BOSMDL C ------ *KEEP,CNSTBF. *KEEP,FTANG. * Statement functions for track angle corrections... * (assumes COMMON FRLORA present) * *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. ********************************************************************** * Get the nominal z-vertex if available * The nominal interaction point is deemed to be at x=y=0 CALL VZERO(VXYZ,2) * If the CJC has found a vertex, use it C INVERT = NLINK('CTKV',0) * Lose dependence on CTREC to aid reprocessing * If the error is bigger than 50 cm its not worth doing & CALL UCOPY(IW(INDCR(INVERT,1,JPRIM)),VXYZ,3) * Event T0 correction *