Changeset 1050 for sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_EventModeProcessing.ipf

Timestamp:

Jun 30, 2017 12:28:55 PM (6 years ago)

Author:

srkline

Message:

cleaning up a lot of the TODO items from the code.

File:

• sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_EventModeProcessing.ipf (modified) (24 diffs)

sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_EventModeProcessing.ipf

@@ -49 +49 @@
 
 
-// TODO 
-// x- these dimensions are hard-wired (OK)
+//  
+// x- these dimensions are hard-wired (OK for now, will need to be additional definitions for Back detector)
+//
+// XBINS is for an individual panel
+// NTUBES is the total number of tubes = (4)(48)=192
 //
 Constant XBINS=48
@@ -63 +66 @@
 
 
+// Initialization of the VSANS event mode panel
 Proc V_Show_Event_Panel()
         DoWindow/F VSANS_EventModePanel
@@ -71 +75 @@
 End
 
-// TODO:
+// 
 //  x- need an index table with the tube <-> panel correspondence
 //    NO, per Phil, the tube numbering is sequential:
@@ -127 +131 @@
 
 
-// for decimation
-        Variable/G root:Packages:NIST:VSANS:Event:gEventFileTooLarge = 150              // 150 MB considered too large
+// for decimation (not used for VSANS - may be added back in the future
+        Variable/G root:Packages:NIST:VSANS:Event:gEventFileTooLarge = 1500             // 1500 MB considered too large
         Variable/G root:Packages:NIST:VSANS:Event:gDecimation = 100
         Variable/G root:Packages:NIST:VSANS:Event:gEvent_t_longest_decimated = 0
 
-// for large file splitting
+// for large file splitting (unused)
         String/G root:Packages:NIST:VSANS:Event:gSplitFileList = ""             // a list of the file names as split
         
-// for editing
+// for editing (unused)
         Variable/G root:Packages:NIST:VSANS:Event:gStepTolerance = 5            // 5 = # of standard deviations from mean. See PutCursorsAtStep()
         
@@ -142 +146 @@
 
 
-
+//
+// the main panel for VSANS Event mode
+// -- a duplicate of the SANS panel, with the functions I'm not using disabled.
+//  could be added back in the future
+//
 Proc VSANS_EventModePanel()
         PauseUpdate; Silent 1           // building window...
@@ -191 +199 @@
         Button button1,pos={389,103},size={120,20},fSize=12,proc=V_ProcessEventLog_Button,title="Bin Event Data"
 
-
+// NEW FOR VSANS
         Button button21,pos={580,70},size={120,20},proc=V_SplitToPanels_Button,title="Split to Panels"
         Button button22,pos={580,90},size={120,20},proc=V_GraphPanels_Button,title="Show Panels"
@@ -240 +248 @@
 EndMacro
 
-
+//
+// takes the data that is loaded and binned as a fake
+// (192 x 128) panel to the four LRTB panels, each with 48 tubes
+//
 Function V_SplitToPanels_Button(ba) : ButtonControl
         STRUCT WMButtonAction &ba
@@ -257 +268 @@
 End
 
+
+//
+// after splitting the data into the 4 panels, draws a simple graph to display the panels
+//
 Function V_GraphPanels_Button(ba) : ButtonControl
         STRUCT WMButtonAction &ba
@@ -272 +287 @@
         return 0
 End
+
 
 // mode selector
@@ -477 +493 @@
 end
 
+//
 // for oscillatory mode
+//
+// Allocate the space for the data as if it was a single 192 x 128 pixel array.
+// The (4) 48-tube panels will be split out later. this means that as defined,
+// XBINS is for an individual panel
+// NTUBES is the total number of tubes = (4)(48)=192
 //
 Function V_Osc_ProcessEventLog(ctrlName)
@@ -615 +637 @@
 // - since I don't know if I've sorted or un-sorted. Osc mode always forces a re-sort and a re-index
 //
+//
+// Allocate the space for the data as if it was a single 192 x 128 pixel array.
+// The (4) 48-tube panels will be split out later. this means that as defined,
+// XBINS is for an individual panel
+// NTUBES is the total number of tubes = (4)(48)=192
+//
 Function V_Stream_ProcessEventLog(ctrlName)
         String ctrlName
@@ -803 +831 @@
         return(0)       
 End
+
 
 // TODO
@@ -877 +906 @@
 
 
-// TODO:
-//
+// 
+// -- The "file too large" check has currently been set to 1.5 GB
 //
 //      
@@ -887 +916 @@
         Variable err=0
         Variable fileref,totBytes
-        NVAR fileTooLarge = root:Packages:NIST:VSANS:Event:gEventFileTooLarge           //limit load to 150MB
+        NVAR fileTooLarge = root:Packages:NIST:VSANS:Event:gEventFileTooLarge           //limit load to 1500MB
 
         SVAR filename = root:Packages:NIST:VSANS:Event:gEvent_logfile
@@ -910 +939 @@
         endif
 
-// TODO - decide if I want (or need) to keep this       
+//  keep this   , but set to 1.5 GB
+// since I'm now in 64-bit space
 /// Abort if the files are too large
-//      Open/R fileref as fileName
-//              FStatus fileref
-//      Close fileref
-//
-//      totBytes = V_logEOF/1e6         //in MB
-//      if(totBytes > fileTooLarge)
-//              sprintf abortStr,"File is %g MB, larger than the limit of %g MB. Split and Decimate.",totBytes,fileTooLarge
-//              Abort abortStr
-//      endif
-//      
-//      Print "TotalBytes = ",totBytes
+        Open/R fileref as fileName
+                FStatus fileref
+        Close fileref
+//
+        totBytes = V_logEOF/1e6         //in MB
+        if(totBytes > fileTooLarge)
+                sprintf abortStr,"File is %g MB, larger than the limit of %g MB. Split and Decimate.",totBytes,fileTooLarge
+                Abort abortStr
+        endif
+//      
+        Print "TotalBytes (MB) = ",totBytes
         
 
@@ -936 +966 @@
 // assign to the proper panels
 
-// TODO:
-// x- (YES - this is MUCH faster) what if I do the JointHistogram first, then break out the blocks of the 
-//  3D sliced data into the individual panels. Then the sort operation could be skipped,
-//  since it would implicitly be done during the histogram operation
+// 
+// x- (YES - this is MUCH faster)  if I do the JointHistogram first, then break out the blocks of the 
+//  3D sliced data into the individual panels. Then the sort operation can be skipped,
+//  since it is implicitly be done during the histogram operation
 // x- go back and redimension as needed to get the 128 x 192 histogram to work
 // x- MatrixOp or a wave assignemt should be able to break up the 3D
@@ -947 +977 @@
                 Duplicate/O eventTime timePt
 
-// TODO:
+// 
 // x- for processing, initially treat all of the tubes along x, and 128 pixels along y
 //   panels can be transposed later as needed to get the orientation correct
@@ -962 +992 @@
 //s_toc()
 
-
 //
 // switch the "active" panel to the selected group (1-4) (5 concatenates them all together)
@@ -969 +998 @@
 // and redimension them to be sure that they are double precision
 //
-
+//
 //      V_SwitchTubeGroup(1)
-
-//      
+//
+//      
+
 //tic()
         Wave timePt=timePt
@@ -981 +1011 @@
         
         NVAR gResol = root:Packages:NIST:VSANS:Event:gResol             //timeStep in clock frequency (Hz)
+        Printf "Time Step = 1/Frequency (ns) = %g\r",(1/gResol)*1e9
+        
 /////
 // now do a little processing of the times based on the type of data
@@ -1004 +1036 @@
 // MODE_TISANE
 
+
 // MODE_TOF
         if(mode == MODE_TOF)            // TOF mode - don't adjust the times, we get periodic t0 to reset t=0
@@ -1028 +1061 @@
 // -- MUCH faster to count the number of lines to remove, then delete (N)
 // rather then delete them one-by-one in the do-loop
+//
 Function V_CleanupTimes(xLoc,yLoc,timePt)
         Wave xLoc,yLoc,timePt
@@ -1076 +1110 @@
 
 
-// TODO (DONE)
+// (DONE)
 // this "fails" for data sets that have 3 or 4 slices, as the ModifyImage command
 // interprets the data as being RGB - and so does nothing.
@@ -1137 +1171 @@
 //  for 64-bit values:
 // http://calc.penjee.com
+//
+//
 //
 //              K0 = 536870912