Changeset 862 for sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS

Timestamp:

Aug 13, 2012 12:30:32 PM (11 years ago)

Author:

srkline

Message:

Added some testing procedures to TISANE.ipf to double-check how the words are being interpreted, and so that I fully understand what is going on. No functionality has been changed or added yet.

File:

• sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/TISANE.ipf (modified) (1 diff)

sans/Dev/trunk/NCNR_User_Procedures/Reduction/SANS/TISANE.ipf

@@ -346 +346 @@
         SetDataFolder root:
 End 
+
+
+////////// SRK - junk I'm testing
+
+//
+Function PreLoadEvents()
+        
+        NVAR time_msw = root:Packages:NIST:gTISANE_time_msw
+        NVAR time_lsw = root:Packages:NIST:gTISANE_time_lsw
+        NVAR t_longest = root:Packages:NIST:gTISANE_t_longest
+        
+        SVAR filepathstr = root:Packages:NIST:gTISANE_logfile
+        SetDataFolder root:Packages:NIST:TISANE
+
+        Variable fileref
+        String buffer
+        Variable dataval,timeval,type,numLines,verbose
+        Variable xval,yval,rollBit,nRoll,roll_time
+        Variable Xmax, yMax
+        xMax = 127              // number the detector from 0->127 
+        yMax = 127
+        
+        verbose = 0
+        numLines = 0
+
+// this gets me the number of lines. not terribly useful        
+//      Open/R fileref as filepathstr
+//      do
+//              numLines += 1
+//              FReadLine fileref, buffer
+//              if (strlen(buffer) == 0)
+//                      numLines -= 1                   //last FReadLine wasn't really a line
+//                      break
+//              endif
+//      while(1)
+//      Close fileref
+        
+        // what I really need is the number of XY events
+        Variable numXYevents,num1,num2,num3,num0
+        numXYevents = 0
+        num0 = 0
+        num1 = 0
+        num2 = 0
+        num3 = 0
+        
+        Open/R fileref as filepathstr
+        do
+                do
+                        FReadLine fileref, buffer                       //skip the "blank" lines that have one character
+                while(strlen(buffer) == 1)              
+
+                if (strlen(buffer) == 0)
+                        break
+                endif
+                
+                sscanf buffer,"%x",dataval
+                
+                // two most sig bits (31-30)
+                type = (dataval & 0xC0000000)/1073741824                //right shift by 2^30
+                
+                if(type == 0)
+                        num0 += 1
+                        numXYevents += 1
+                endif
+                if(type == 2)
+                        num2 += 1
+                        numXYevents += 1
+                endif
+                if(type == 1)
+                        num1 += 1
+                endif
+                if(type == 3)
+                        num3 += 1
+                endif   
+                
+        while(1)
+        Close fileref
+//              done counting the number of XY events
+
+
+//      
+        Printf "numXYevents = %d\r",numXYevents
+        Printf "XY = num0 = %d\r",num0
+        Printf "XY time = num2 = %d\r",num2
+        Printf "time MSW = num1 = %d\r",num1
+        Printf "Rollover = num3 = %d\r",num3
+
+        Printf "num0 + num2 = %d\r",num0+num2
+        
+        Make/O/U/N=(numXYevents) xLoc,yLoc
+        Make/O/D/N=(numXYevents) timePt
+        xLoc=0
+        yLoc=0
+        timePt=0
+
+        
+        Variable ii=0
+        nRoll = 0               //number of rollover events
+        roll_time = 2^26                //units of 10-7 sec
+        
+        time_msw=0
+        
+        Open/R fileref as filepathstr
+        do
+                do
+                        FReadLine fileref, buffer                       //skip the "blank" lines that have one character
+                while(strlen(buffer) == 1)              
+
+                if (strlen(buffer) == 0)
+                        break
+                endif
+                
+                sscanf buffer,"%x",dataval
+                
+                //This is hideous, but that is what you get for IGOR not doing structs like C
+                //
+                //Constant ATXY = 0
+                //Constant ATXYM = 2
+                //Constant ATMIR = 1
+                //Constant ATMAR = 3
+                //
+//              type = (dataval & ~(2^32 - 2^30 -1))/2^30
+
+                // two most sig bits (31-30)
+                type = (dataval & 0xC0000000)/1073741824                //right shift by 2^30
+                
+//              if(ii<10)
+//                      Print "buffer=", buffer//, type
+//                      Print "strlen(buffer) =",strlen(buffer)
+//                      printf "%b      %u\r",dataval,type
+//              endif
+                
+                
+// just to print out some of the diagnostic information 
+//              if(ii < 7)
+//                      verbose = 1
+//              else
+//                      verbose = 0
+//              endif
+                
+//      
+// not sure that the XY positions are correctly interpreted             
+// X seems OK, but Y doesn't match teabag output (it does if the (yMax - ) is dropped, and just the 8 bits are used)
+//              
+                switch(type)
+                        case ATXY:
+                                if(verbose)             
+                                        printf "XY : "          
+                                endif
+//                              xval = ~(dataval & ~(2^32 - 2^8)) & 127
+//                              yval = ((dataval & ~(2^32 - 2^16 ))/2^8) & 127
+//                              time_lsw = (dataval & ~(2^32 - 2^29))/2^16
+
+                                xval = xMax - (dataval & 255)                                           //last 8 bits (7-0)
+                                yval = (dataval & 65280)/256                                            //bits 15-8, right shift by 2^8
+                                time_lsw = (dataval & 536805376)/65536                  //13 bits, 28-16, right shift by 2^16
+
+                                timeval = trunc( nRoll*roll_time + (time_msw * (8192)) + time_lsw )             //left shift msw by 2^13, then add in lsw, as an integer
+                                if (timeval > t_longest) 
+                                        t_longest = timeval
+                                endif
+                                xLoc[ii] = xval
+                                yLoc[ii] = yval
+                                timePt[ii] = timeval
+
+                                if(verbose)             
+//                                      printf "%u : %u : %u : %u\r",dataval,time_lsw,time_msw,timeval
+                                        printf "%u : %u : %u : %u\r",dataval,timeval,xval,yval
+                                endif                           
+
+//                              b = FindBin(timeval,nslices)
+//                              slicedData[xval][yval][b] += 1
+
+                                ii+=1
+                                
+                                break
+                        case ATXYM:
+                                if(verbose)
+                                        printf "XYM : "
+                                endif
+//                              xval = ~(dataval & ~(2^32 - 2^8)) & 127
+//                              yval = ((dataval & ~(2^32 - 2^16 ))/2^8) & 127
+//                              time_lsw =  (dataval & ~(2^32 - 2^29 ))/2^16
+
+                                xval = xMax - (dataval & 255)                                           //last 8 bits (7-0)
+                                yval = (dataval & 65280)/256                                            //bits 15-8, right shift by 2^8
+
+                                time_lsw = (dataval & 536805376)/65536                  //13 bits, 28-16, right shift by 2^16
+
+                                if(verbose)
+//                                      printf "%u : %u : %u : %u\r",dataval,time_lsw,time_msw,timeval
+                                        printf "%u : %u : %u : %u\r",dataval,timeval,xval,yval
+                                endif
+                                
+                                xLoc[ii] = xval
+                                yLoc[ii] = yval
+                                
+                                // don't fill in the time yet, or increment the index ii
+                                // the next event MUST be ATMIR with the MSW time bits
+                                //
+                                break
+                        case ATMIR:
+                                if(verbose)
+                                        printf "MIR : "
+                                endif
+//                              time_msw =  (dataval & ~(2^32 - 2^29 ))/2^16
+                                time_msw =  (dataval & 536805376)/65536                 //13 bits, 28-16, right shift by 2^16
+                                timeval = trunc( nRoll*roll_time + (time_msw * (8192)) + time_lsw )
+                                if (timeval > t_longest) 
+                                        t_longest = timeval
+                                endif
+                                if(verbose)
+//                                      printf "%u : %u : %u : %u\r",dataval,time_lsw,time_msw,timeval
+                                        printf "%u : %u : %u : %u : %u\r",dataval,time_lsw,time_msw,timeval,t_longest
+                                endif
+                                
+                                // the XY position was in the previous event ATXYM
+                                timePt[ii] = timeval
+                                
+                                ii+=1
+                                
+//                              b = FindBin(timeval,nslices)
+//                              slicedData[xval][yval][b] += 1
+
+                                break
+                        case ATMAR:
+                                if(verbose)
+                                        printf "MAR : "
+                                endif
+                                
+                                // do something with the rollover event?
+                                
+                                // check bit 29
+                                rollBit = (dataval & 0x20000000)/536870912              //bit 29 only , shift by 2^29
+                                nRoll += 1
+                                
+                                if(verbose)
+                                        printf "%u : %u \r",rollBit,nRoll
+                                endif
+                                
+                                break
+                endswitch
+                                
+        while(1)
+        
+        
+        Close fileref
+        
+        SetDataFolder root:
+End 
+
+///