*-- Author : Stephen J. Maxfield SUBROUTINE FPOKER **: FPOKER 40000 SM. New routine for calibration checking. **---------------------------------------------------------------------- *KEEP,BCS. *KEEP,BOSMDL. C ------BOSMDL C ------ *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,FDIFLG. *KEND. * COMMONs for planar found tracks * Pointers to radials associated with planar tracks * Tan of 1/2 wedge angle... * Binning parameters for histograms... * Statement functions for TABLE access... *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 *KEND. * Book LOOK histograms * Histograms for t-to-d... CALL BHS(KHIS1, 0, 200, -DPMAX, DPMAX) CALL BHS(KHIS2, 0, 20, DLO, DHI) * Histograms for Lorentz angle... CALL BHS(KHIS3, 0, 100, -50., 50.) CALL BHS(KHIS4, 0, 20, DLO, DHI) * NPP is number of planar-based tracks... * Hit data... INFRLC = MLINK(IW,'FRLC',0) * Auxiliary hit data... INFAUX = MLINK(IW,'FAUX',0) * Loop over planar-based tracks... * Which PLANAR supermodules have the hits on the track. Actually * always have a full segments worth (9-12 hits) or none at all * so this is overkill! * Loop over the radial hits on this track. Only mods 0 and 1 * Radial hit on this track? * Accept 'sandwich' configurations only... * Get Phi and R at this wire plane from Phi-z R-z fit parameters. * Limit radius range. * Predicted drift...corrected for stagger. * Predicted radius along wire direction (if no Lorenz angle) * Get max allowed drift (Position of cathode plane less a * 3mm tolerance)... * ...and cut out region near cathode. * Drift time, corrected for T0 and radius and pre-scaled * by approxfrom F0R8. * ...and radius of hit at wire from Charge division * (Predicted radius at wire if no Lorentz angle) - (measured * radius) * ...slope of DRR vs. predicted drift is tan(alpha). * Fill histograms of slices in predicted drift. CALL SHS(2000+KBIN, 0, DTSCA) CALL SHS(3000+KBIN, 0, DRP) CALL SHS(10000+KBINL, 0, DRR) CALL SHS(11000+KBINL, 0, DRP) *--------------- * ...end loop over radial hits * ...end loop over planar based tracks. * Write(6,*) ' Fpoker hits', ITOTHT *