*-- Author : Stephen J. Maxfield / I.O.Skillicorn SUBROUTINE FPTREZ(KTRK) * * Calculate residuals of Radial hits to PATREC tracks * and make lookagrams of them. * *D Remove vertex shifts in segment residuals. * *KEEP,FRDIMS. *KEEP,FH1WORK. * Planar geometry * * Radial geometry * * Radial data... * * Planar Data * * Pointers into DIGI bank for IOS labelled hits * * Track segment data * * Fit data * * *KEEP,FPTVTX. **the common/VERTEX/ becomes /VERTVV/ (in analogy to /VERTFF/) on the ** 17/6/91, since it is in conflict with the VERTEX module (g.bernardi) ** (note that all these common names should start by F in this deck...) *KEEP,FPTFLG. *KEND. * * Offsets for radials... * Jiggery-pokery to get right wedge number... * Expected drift corrected for shifts... * Measured drifts... CALL SHS(116,0,REZ(I)) CALL SHS(119,0,REZ(I)) CALL SHS(117,0,REZ(I)) CALL SHS(119,0,REZ(I)) CALL SHS(118,0,REZ(I)) CALL SHS(119,0,REZ(I)) * Now look at radial segment data... CALL SHS(241,0,CH) * Track parameters for this track... * IF(IERRF(JP,JWIR) .GT. 1) GOTO 1003 * Z of wire... * Measured drift... * ...and radius (this has been corrected for Lorenz angle). * Expected drift corrected for shifts... * Expected radius corrected for shifts... CALL SHS(209,0,FLOAT(JWIR)) CALL SHS(232+I,0,RESID(KK)) CALL SHS(236 ,0,RESID(KK)) CALL SHD(242 ,0,DRFE(KK),RESID(KK)) CALL SHS(236+I,0,RESR(KK)) CALL SHS(240 ,0,RESR(KK)) CALL SHD(250 ,0,RMEAS(KK),RESR(KK)) * Close to wire stuff... *