*-- Author : Stephen J. Maxfield
      SUBROUTINE FPKANL
*D: FPKANL.......SM. Add options for left/right separation                                                    
**: FPKANL 40000 SM. New routine for calibration checking.                                                    
**----------------------------------------------------------------------                                      
*     Does not work in multi-processor environment.                                                           
*                                                                                                             
*                                                                                                             
*KEEP,BCS.                                                                                                    
*KEND.                                                                                                        
*                                                                                                             
*--------------------------------------------------------------------                                         
      CALL SAREA('FPOKE', 0)                                                                           
*     Book histograms...                                                                                      
       CALL BVEC(KSIDE*100, 0, 6)                                                                      
       CALL STEXT(KSIDE*100, 4,'D-time (scaled) vs. Dist(predicted)')                                  
      CALL BVEC( 12000, 0, 6)                                                                          
      CALL STEXT(12000, 4,'DeltaR vs. Drift predicted)')                                               
*     Analyse the data.                                                                                       
*     ------- --- ----                                                                                        
*     Analysis of histogram results.                                                                          
*     Do peakparm analysis of the pred drift histograms...                                                    
*      Get average predicted drift distance in the slice...                                                   
       CALL GHSTAT('HS', JHIS2, 0, NENT, SUMW, RNEFF, XST, YST)                                        
       CALL HPEAK('HS',JHIS1, 0, NPK, PDAT)                                                            
*          peak position and error on...                                                                      
*          comment out next line for 'full width errors'                                                      
*          Hence Drift time vs. predicted  drift distance:-                                                   
           CALL SVEC(KSIDE*100, 0, CVEC)                                                               
*      Now purge figures - no longer needed.                                                                  
       CALL PURGEF(JHIS1)                                                                              
       CALL PURGEF(JHIS2)                                                                              
*     Do peakparm analysis of the Lorentz Angle Histograms...                                                 
*      Get average predicted drift distance in the slice...                                                   
       CALL GHSTAT('HS', JHIS2, 0, NENT, SUMW, RNEFF, XST, YST)                                        
       CALL HPEAK('HS',JHIS1, 0, NPK, PDAT)                                                            
*          peak position and error on...                                                                      
*          comment out next line for 'full width errors'                                                      
*          Hence Delta R vs. predicted drift distance:-                                                       
           CALL SVEC(12000, 0, CVEC)                                                                   
*      Now purge figures - no longer needed.                                                                  
       CALL PURGEF(JHIS1)                                                                              
       CALL PURGEF(JHIS2)                                                                              
*     Extraction of calibration data. Not for online at moment...                                             
*