DOUBLE PRECISION FUNCTION KSGMOM(PW, VP, MVP)
#endif
#if !defined(CLEO_DOUBLE)
      REAL FUNCTION KSGMOM(PW, VP, MVP)
#endif
 
C  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
C  Returns the error in momentum for a track in the W parameterization.
C
C  PW        real array (read)
C            Track parameters in W form (at least 3-vector part)
C
C  VP        real array (read)
C            Covariance matrix
C
C  MVP       integer variable (read)
C            First dimension of VP (must be at least 3)
C
C
C   The momentum can be written as
C
C         P = sqrt(Px**2 + Py**2 + Pz**2)
C
C   so    dP = (Px*dPx + Py*dPy + Pz*dPz) / P
C
C   The covariance matrix of the momentum is then
C
C       sig**2 = A*Vx0*A(t)
C
C    where A = (Px, Py, Pz, 0) / P
C  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>