FILHIS COMMENTS
*-- Author :    Stephen J. Maxfield   27/02/93
      SUBROUTINE FILHIS
**-------------------------------------------------------------                                               
*                                                                                                             
*     Fill some LOOK histograms for monitoring                                                                
*     Forward Tracker Pattern Recognition.                                                                    
*                                                                                                             
*--------------------------------------------------------------                                               
*                                                                                                             
*KEEP,BCS.                                                                                                    
*KEEP,BOSMDL.                                                                                                 
C     ------BOSMDL                                                                                            
C     ------                                                                                                  
*KEEP,H1EVDT.                                                                                                 
*                                                                                                             
*  IDATA  type of information (HEAD bank word 6) :                                                            
*                                                                                                             
*                       0 - real data H1                                                                      
*                       1 - MC data H1SIM                                                                     
*                       2 - real data CERN tests                                                              
*                       3 - MC data ARCET                                                                     
*                                                                                                             
*  MONTE = .TRUE.   if IDATA=1                                                                                
*  KEVENT = event processed counter for H1REC                                                                 
*                                                                                                             
*KEEP,FDIFLG.                                                                                                 
*KEEP,FPTFLG.                                                                                                 
*KEND.                                                                                                        
*------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.                                                                                                        
*                                                                                                             
*                                                                                                             
*                                                                                                             
*-----------------------------------------------------------------                                            
*       Fill Histograms.                                                                                      
*                                                                                                             
*     Zero all Work Bank Indices...                                                                           
*                                                                                                             
*      'Pointering' bank FPUR...                                                                              
        IND = NLINK('FPUR',NBN)
        CALL BKTOW(IW,'FPUR',NBN,IW,INDPUR,*900)
*                                                                                                             
*       Number of tracks found...                                                                             
        CALL SHS(101, 0, FLOAT(NTRK))                                                                  
*                                                                                                             
*                                                                                                             
*                                                                                                             
*       Level '0' track parameters...                                                                         
        INDTUR = NLINK('FTUR',NBN)
*                                                                                                             
*       List of radial hits...                                                                                
        IND = NLINK('FRUX',NBN)
        CALL BKTOW(IW,'FRUX',NBN,IW,INDRUX,*900)
*                                                                                                             
*       List of planar hits...                                                                                
        IND = NLINK('FPUX',NBN)
        CALL BKTOW(IW,'FPUX',NBN,IW,INDPUX,*900)
*                                                                                                             
*       Radial hit data (drifts etc.)...                                                                      
        INDRLC = NLINK('FRLC',NBN)
*                                                                                                             
*       Planar hit data (drifts etc.)...                                                                      
        INDPLC = NLINK('FPLC',NBN)
*                                                                                                             
*       Loop over pattern recognised tracks                                                                   
*                                                                                                             
*                                                                                                             
          CALL SHS(102, 0, PHI  )                                                                      
          CALL SHS(103, 0, THETA)                                                                      
          CALL SHS(104, 0, PMEAS)                                                                      
            CALL SHS(105, 0, 1./PMEAS)                                                                 
            CALL SHS(106, 0, LOG(1./PMEAS))                                                            
*                                                                                                             
*                                                                                                             
*         Get hit data for each radial hit on the track...                                                    
*           Drift from FRLC bank...                                                                           
            CALL SHS( 110,0, DRIFT)                                                                    
            CALL SHS( 111,0, RADIUS)                                                                   
            CALL SHD( 112,0, DRIFT, RADIUS)                                                            
*         Get hit data for each planar hit on the track...                                                    
*           Drift from FPLC bank...                                                                           
            CALL SHS( 113,0, DRIFT)                                                                    
*                  ...Monte Carlo                                                                             
*                                                                                                             
*               ...Loop over tracks                                                                           
*                                                                                                             
*----------- Done ----------------------------------------                                                    
*     Drop all Work banks (unused ones should have index 0 so O.K.)                                           
      CALL WDROP(IW,INDPUR)                                                                            
      CALL WDROP(IW,INDRUX)                                                                            
      CALL WDROP(IW,INDPUX)                                                                            
*                                                                                                             
*