This routine calculates the weight matrix of the track with hist stored in /ZTTRJ/ (filled by ZTFOLW) taking measurements, multiple scattering and geometry into account.
INCR : if TRUE, this call is done in the process of transporting a weight matrix, rather than calculating it from scratch.
from /ZTTRJ/ : R,Z,PHI,BETA : track parameters at intersection of each detector element; NPT : number of intersected elements; IORD : pointer into R etc ordered in increasing track- length; ITYP : cylinder/plane (1/2) and measurement/no measurement (+ve/-ve); JREF : the position of the reference surface intersection in rack-length ordering. if INCR = TRUE : From /MATRIX/ : WGT : weight matrix.
to /MATRIX/ : WGT : weight matrix; IGEOM : geometry type of present point (barrel/forward = 1/2)
(To switch between plane and cylindrical geometry.)
(To do geometrical propagation of errors.)
(To include multiple scattering.)
(To add contribution of measurement.)
SUBROUTINE ZTGEOC(WG,IPT) * CALLED ROUTINE(S): * * none. *