*-- Author : Stephen Burke
SUBROUTINE FFSTART(SSTART,CSTART,ZSTART)
**********************************************************************
* *
* Set up the starting vector for the Kalman filter *
* *
* SSTART and CSTART give the initial state vector, and its *
* covariance, at z = ZSTART. *
* *
**********************************************************************
*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,FKCNTL.
*KEEP,FKFLAG.
*KEEP,FKCONS.
*KEEP,FKPROJ.
*KEEP,FKSMTH.
*KEEP,FKRSID.
*KEEP,FKTRUE.
*KEEP,FFSTEE.
*KEEP,FFDBG.
*KEEP,FKDBG.
*KEEP,FKINT.
*KEEP,FRLORA.
*KEEP,BCS.
*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.
**********************************************************************
* Set up the start and end points
CALL ERRLOG(391,'F:FFSTART: Starting plane mapped out')
* Swim to 1st plane (ignoring MS)
CALL FKTRAN(DZ,ZSTART,SSTART,SPRO(1,JSTART),DTRAN)
CALL FKMUL(CSTART,DTRAN,CPRO(1,1,JSTART))
* Fix for low-momentum tracks
* Starting errors - these are just fixed at the moment
CALL UCOPY(CPRO(1,1,JSTART),CTEMP,50)
CALL VZERO(CPRO(1,1,JSTART),50)
*