source: sans/Dev/trunk/NCNR_User_Procedures/Reduction/VSANS/V_Test_RAW_Panel.ipf @ 1025

Last change on this file since 1025 was 1025, checked in by srkline, 6 years ago

many changes to get the basics of a reduction protocol working

File size: 36.5 KB
Line 
1#pragma rtGlobals=3             // Use modern global access method and strict wave access.
2#pragma version=1.0
3#pragma IgorVersion=6.1
4
5//
6// this will become the equivalent of "RawWindowHook"
7//
8// Procedures to:
9// display the detector data
10// visualization tools
11// mouse interaction
12// status information necessary to understand the data
13// buttons to more functionality to process data
14
15
16// TODO
17//
18// -- have the status automatically fill in when a new file is loaded, rather than needing a click of the "status" button
19// -- need a place somewhere to show the currently displayed folder
20// -- checkboxes for "active" corrections?
21// -- display of Q, counts, QxQy, X and Y
22// -- do I need a color bar? or is this not used at all? I like it to be there, or the colors are a waste of information
23//              (then... where do I put the color bar?)
24// -- define the "hook" function, and attach it to the panel (or the main detector subwindow?)
25//
26
27
28//
29// call this after loading data to either draw the data display panel or to update the contents
30//
31// TODO
32// -- make sure that the "type" input is correctly used for the updating of the data, values, etc.
33// -- add a procedure to define the global variables for pos, counts, QxQy, etc.
34//
35Function V_UpdateDisplayInformation(type)
36        String type
37       
38        DoWindow/F VSANS_Data
39//      Print V_flag
40        if(V_flag==0)
41       
42                VSANSDataPanelGlobals()
43               
44                Execute "VSANS_DataPanel()"             //draws the panel
45
46        endif
47       
48        // TODO: update the information here  - in either case
49        // what isn't automatically picked up? What is "stale" on the display?
50        String/G root:Packages:NIST:VSANS:Globals:gCurDispType = type
51       
52        // fake a click on all three tabs - to populate the data
53        V_FakeTabClick(2)
54        V_FakeTabClick(1)
55        V_FakeTabClick(0)
56
57        V_LogLinDisplayPref()           // resets the display depending on the preference state
58       
59        // either clear the status, or fake a click
60//      root:Packages:NIST:VSANS:Globals:gStatusText = "status info box"
61
62        V_FakeStatusButtonClick()
63
64        V_FakeRestorePanelsButtonClick()                //so the panels display correctly
65
66       
67//      DoWindow/T VSANS_Data,type + " VSANS_Data"
68
69        String newTitle = "WORK_"+type
70        DoWindow/F VSANS_Data
71        DoWindow/T VSANS_Data, newTitle
72        KillStrings/Z newTitle
73       
74end
75
76Function V_LogLinDisplayPref()
77
78        // get the state of the log/lin button, and make sure preferences are obeyed
79        // log/lin current state is in the S_UserData string (0=linear, 1=log)
80        ControlInfo/W=VSANS_Data button_log
81        Variable curState
82        curState = str2num(S_UserData)
83
84        NVAR gLogScalingAsDefault = root:Packages:NIST:VSANS:Globals:gLogScalingAsDefault
85        if(curState != gLogScalingAsDefault)
86                STRUCT WMButtonAction ba
87                ba.eventCode = 2                //fake mouse click
88                V_LogLinButtonProc(ba)
89        endif
90       
91        return(0)
92end
93
94
95//
96// creates/initializes the globals for display of the data panel
97//
98Function VSANSDataPanelGlobals()
99
100        SetDataFolder root:Packages:NIST:VSANS:Globals
101       
102        Variable/G gXPos=0
103        Variable/G gYPos=0
104        Variable/G gQX=0
105        Variable/G gQY=0
106        Variable/G gQQ=0
107        Variable/G gNCounts=0
108        String/G gCurDispFile = "default string"
109        String/G gCurTitle = ""
110        String/G gCurDispType = ""
111        String/G gStatusText = "status"
112        String/G gLastLoadedFile=""
113       
114        SetDataFolder root:
115End
116
117
118// TODO
119//
120// -- now that the sliders work, label them and move them to a better location
121// -- logical location for all of the buttons
122// -- add raw data load button, load/draw mask button
123// x- fill in the proper window title in the DoWindow/T command
124// -- add help text for all of the controls
125// -- tab order? can I set this?
126//
127Window VSANS_DataPanel() : Panel
128        PauseUpdate; Silent 1           // building window...
129        NewPanel /W=(37,45,1038,719) /N=VSANS_Data
130//      ShowTools/A
131        ModifyPanel cbRGB=(65535,60076,49151)
132
133        String curFolder = root:Packages:NIST:VSANS:Globals:gCurDispType
134        DoWindow/T VSANS_Data,curFolder + " VSANS_Data"
135        SetWindow VSANS_Data,hook(dataHook)=VSANSDataHook,hookevents=2
136       
137        SetDrawLayer UserBack
138        SetDrawEnv linethick= 2,dash= 1,fillpat= 0
139        DrawRect 200,70,310,160
140        SetDrawEnv linethick= 2,dash= 1,fillpat= 0
141        DrawRect 320,70,430,160
142        SetDrawEnv linethick= 2,dash= 1,fillpat= 0
143        DrawRect 440,70,550,160
144       
145        SetDrawEnv fsize= 18
146        DrawText 230,115,"Front"
147        SetDrawEnv fsize= 18
148        DrawText 348,115,"Middle"
149        SetDrawEnv fsize= 18
150        DrawText 476,115,"Back"
151       
152        ToolsGrid visible=1
153
154
155        TabControl tab0,pos={13,41},size={572,617},proc=V_DataTabProc,tabLabel(0)="Front"
156        TabControl tab0,tabLabel(1)="Middle",tabLabel(2)="Back",value= 2,focusRing=0
157
158// on the side 
159        Button button_status,pos={607,146},size={70,20},proc=V_StatusButtonProc,title="Status"
160        Button button_IvsQ,pos={689,113},size={70,20},proc=V_IvsQPanelButtonProc,title="I vs. Q"
161        Button button_file_m,pos={619,55},size={50,20},proc=V_File_minus_ButtonProc,title="File <"
162        Button button_file_p,pos={679,55},size={50,20},proc=V_File_plus_ButtonProc,title="File >"
163        Button button_log,pos={689,146},size={70,20},proc=V_LogLinButtonProc,title="isLin",userData="0"
164        Button button_tab_p,pos={648,81},size={50,20},proc=V_Tab_p_ButtonProc,title="Tab >"
165        Button button_isolate,pos={606,114},size={70,20},proc=V_IsolateButtonProc,title="Isolate"
166        Button button_toWork,pos={770,113},size={70,20},proc=V_ToWorkFileButtonProc,title="to WORK"
167
168        TitleBox title_file,pos={606,178},size={76,20},variable= root:Packages:NIST:VSANS:Globals:gLastLoadedFile
169        TitleBox title_dataPresent,pos={606,210},size={76,20},variable= root:Packages:NIST:VSANS:Globals:gCurDispFile
170        TitleBox title_status,pos={606,240},size={200,200},variable= root:Packages:NIST:VSANS:Globals:gStatusText
171       
172        Button button_tagFile,pos={603,412},size={70,20},proc=V_TagFileButtonProc,title="Tag File"
173        Button button_tagFile,disable=2
174        Button button_saveIQ,pos={720,412},size={70,20},proc=V_SaveIQ_ButtonProc,title="Save I(Q)"
175        Button button_BeamCtr,pos={603,450},size={70,20},proc=V_BeamCtrButtonProc,title="Beam Ctr"
176        Button button_SpreadPanels,pos={603,488},size={100,20},proc=V_SpreadPanelButtonProc,title="Spread Panels"
177        Button button_RestorePanels,pos={603,526},size={100,20},proc=V_RestorePanelButtonProc,title="Restore Panels"
178
179
180// on the tabs, always visible
181        TitleBox title_xy,pos={24,71},size={76,20},variable= root:Packages:NIST:VSANS:Globals:gLastLoadedFile
182        Slider slider_hi,pos={558,224},size={16,80},proc=V_HiMapSliderProc
183        Slider slider_hi,limits={0,1,0},value= 1,ticks= 0
184        Slider slider_lo,pos={558,315},size={16,80},proc=V_LowMapSliderProc
185        Slider slider_lo,limits={0,1,0},value= 0,ticks= 0
186
187        SetVariable xpos,pos={22,97},size={50,17},title="X "
188        SetVariable xpos,limits={-Inf,Inf,0},value= root:Packages:NIST:VSANS:Globals:gXPos
189        SetVariable xpos,help={"x-position on the detector"},frame=0,noedit=1
190        SetVariable ypos,pos={22,121},size={50,17},title="Y "
191        SetVariable ypos,limits={-Inf,Inf,0},value= root:Packages:NIST:VSANS:Globals:gYPos
192        SetVariable ypos,help={"y-position on the detector"},frame=0,noedit=1
193        SetVariable counts,pos={22,151},size={150,17},title="Counts "
194        SetVariable counts,limits={-Inf,Inf,0},value= root:Packages:NIST:VSANS:Globals:gNCounts
195        SetVariable counts,help={"Neutron counts"},frame=0,noedit=1
196        SetVariable qxval,pos={83,94},size={85,17},title="qX"
197        SetVariable qxval,help={"q value in the x-direction on the detector"},frame=0,noedit=1
198        SetVariable qxval,format="%+7.5f",limits={-Inf,Inf,0},value= root:Packages:NIST:VSANS:Globals:gQX
199        SetVariable qyval,pos={83,113},size={85,17},title="qY"
200        SetVariable qyval,help={"q value in the y-direction on the detector"},frame=0,noedit=1
201        SetVariable qyval,format="%+7.5f",limits={-Inf,Inf,0},value= root:Packages:NIST:VSANS:Globals:gQY
202        SetVariable q_pos,pos={83,132},size={85,17},title="q "
203        SetVariable q_pos,help={"q-value on the detector at (x,y)"},format="%+7.5f"
204        SetVariable q_pos,limits={-Inf,Inf,0},value= root:Packages:NIST:VSANS:Globals:gQQ,frame=0,noedit=1
205       
206        Make/O/D tmp_asdf
207        // for back panels (in pixels?)
208//      Display/W=(50,185,545,620)/HOST=# tmp_asdf
209        Display/W=(50,185,517,620)/HOST=# tmp_asdf
210        RenameWindow #,det_panelsB
211        ModifyGraph mode=2              // mode = 2 = dots
212        ModifyGraph marker=19
213        ModifyGraph rgb=(0,0,0)
214        ModifyGraph tick=2,mirror=1
215        Label left "pixels"
216        Label bottom "pixels"   
217        SetActiveSubwindow ##
218       
219        // for middle panels (in pixels?)       
220        Display/W=(50,185,517,620)/HOST=# tmp_asdf
221        RenameWindow #,det_panelsM
222        ModifyGraph mode=2              // mode = 2 = dots
223        ModifyGraph marker=19
224        ModifyGraph rgb=(0,0,0)
225        ModifyGraph tick=2,mirror=1
226        Label left "pixels"
227        Label bottom "pixels"   
228        SetActiveSubwindow ##
229       
230        // for front panels (in pixels?)       
231        Display/W=(50,185,517,620)/HOST=# tmp_asdf
232        RenameWindow #,det_panelsF
233        ModifyGraph mode=2              // mode = 2 = dots
234        ModifyGraph marker=19
235        ModifyGraph rgb=(0,0,0)
236        ModifyGraph tick=2,mirror=1
237        Label left "pixels"
238        Label bottom "pixels"   
239        SetActiveSubwindow ##
240       
241EndMacro
242
243//
244// event code 4: mouse moved
245// event code 11: keyboard events
246//
247// mouse moved is the only event that I really care about for the data display.
248//
249// TODO:
250// -- figure out how to respond only to events in the main window
251// -- figure out which is the correct image to respond "from"
252// -- More complete documentation of how the hook is identifying what graph is "under" the mouse
253//    AND what assumptions are behind this identification
254//
255//
256Function VSANSDataHook(s)
257        STRUCT WMWinHookStruct &s
258
259        Variable hookResult = 0
260
261        switch(s.eventCode)
262                case 0:                         // Activate
263                        // Handle activate
264//                      Print "Activate"
265                        break
266
267                case 1:                         // Deactivate
268                        // Handle deactivate
269                        break
270                       
271                case 3:         //mouse down
272//                      Print "mouse down"
273//
274// TODO (Way in the future -- I could make the small graphs into "buttons" by responding to a "mouse up" (not down)
275//    that hits in one of the small graph areas, and treat that as a click on that tab
276// -- instead of this, I can try to get the focus rectangle to work more easily to move with the TAB,
277//    if I can intercept the keystroke (event 11), see below.
278//
279                        break
280                       
281                case 4:         // mouse moved
282                        NVAR xloc = root:Packages:NIST:VSANS:Globals:gXPos
283                        NVAR yloc = root:Packages:NIST:VSANS:Globals:gYPos
284                        NVAR gQX = root:Packages:NIST:VSANS:Globals:gQX
285                        NVAR gQY = root:Packages:NIST:VSANS:Globals:gQY
286                        NVAR gQQ = root:Packages:NIST:VSANS:Globals:gQQ
287                        NVAR gNCounts = root:Packages:NIST:VSANS:Globals:gNCounts
288                        SVAR gCurDispFile = root:Packages:NIST:VSANS:Globals:gCurDispFile
289                        SVAR gCurDispType = root:Packages:NIST:VSANS:Globals:gCurDispType               //the current folder
290                        Variable xaxval,yaxval,tab
291                       
292                        // is the mouse location within the "main" display window?
293                        // if so, do something, if not, do nothing?
294                        // right now, the "main" display is at (50,185,545,620). its name depends on the active tab
295                       
296//                              xloc = s.mouseLoc.h
297//                              yloc = s.mouseLoc.v
298
299//                      if out of bounds, exit now
300//              TODO - currently the values are hard-wired. eliminate this later if the size of the graph changes
301                        if(s.mouseLoc.h < 50 || s.mouseLoc.h > 545 || s.mouseLoc.v < 185 || s.mouseLoc.v > 620)
302                                break
303                        endif   
304                       
305//                      if(in bounds)
306//                              get the point location
307//                              update the globals --
308//                              but which data instance am I pointing to?
309//                              deduce the carriage and panel, and calculate Q
310//                      endif
311
312                        GetWindow $s.winName activeSW
313                        String activeSubwindow = S_value                // returns something like: "VSANS_Data#det_panelsF"
314                               
315                        xaxval= AxisValFromPixel("","bottom",s.mouseLoc.h)
316                        yaxval= AxisValFromPixel("","left",s.mouseLoc.v)
317                        xloc = round(xaxval)
318                        yloc = round(yaxval)
319                       
320                        // which tab is selected? -this is the main graph panel (subwindow may not be the active one!)
321                        ControlInfo/W=VSANS_Data tab0
322                        tab = V_Value
323                        if(tab == 0)
324                                activeSubwindow = "VSANS_Data#det_panelsF"
325                        elseif (tab == 1)
326                                activeSubwindow = "VSANS_Data#det_panelsM"
327                        else
328                                activeSubwindow = "VSANS_Data#det_panelsB"
329                        endif
330                       
331                        // which images are here?
332                        String detStr="",imStr,carriageStr
333                        String currentImageRef
334                        String imageList = ImageNameList(activeSubwindow,";")
335                        Variable ii,nIm,testX,testY,xctr,yctr,sdd,lam,pixSizeX,pixSizeY
336                        nIm = ItemsInList(imageList,";")
337                        gCurDispFile = imageList
338                        if(nIm==0)
339                                break           //problem, get out
340                        endif
341
342                        // images were added in the order TBLR, so look back through in the order RLBT, checking each to see if
343                        // the xy value is found on that (scaled) array
344                                               
345                        // loop backwards through the list of panels (may only be one if on the back)
346                        for(ii=nIm-1;ii>=0;ii-=1)
347                                Wave w = ImageNameToWaveRef(activeSubwindow,StringFromList(ii, imageList,";"))
348                               
349                                // which, if any image is the mouse xy location on?
350                                // use a multidemensional equivalent to x2pnt: (ScaledDimPos - DimOffset(waveName, dim))/DimDelta(waveName,dim)
351                                testX = trunc( (xloc - DimOffset(w,0))/DimDelta(w,0) )
352                                testY = trunc( (yloc - DimOffset(w,1))/DimDelta(w,1) )
353                               
354                                if( (testX > 0 && testX < DimSize(w,0)) && (testY > 0 && testY < DimSize(w,1)) )
355                                        // we're in-bounds on this wave
356                                       
357                                        // count value to the global
358                                        gNCounts = w[testX][testY]
359                                       
360                                        // deduce the detector panel
361                                        currentImageRef = StringFromList(ii, imageList,";")     //the image instance ##
362                                        // string is "data", or "data#2" etc. - so this returns "", "1", "2", or "3"
363                                        imStr = StringFromList(1, currentImageRef,"#")         
364                                        carriageStr = activeSubWindow[strlen(activeSubWindow)-1]
365                                       
366                                        if(cmpstr(carriageStr,"B")==0)
367                                                detStr = carriageStr
368                                        else
369                                                if(strlen(imStr)==0)
370                                                        imStr = "9"                     // a dummy value so I can replace it later
371                                                endif
372                                                detStr = carriageStr+imStr              // "F2" or something similar
373                                                detStr = ReplaceString("9", detStr, "T")        // ASSUMPTION :::: instances 0123 correspond to TBLR
374                                                detStr = ReplaceString("1", detStr, "B")        // ASSUMPTION :::: this is the order that the panels
375                                                detStr = ReplaceString("2", detStr, "L")        // ASSUMPTION :::: are ALWAYS added to the graph
376                                                detStr = ReplaceString("3", detStr, "R")        // ASSUMPTION ::::
377                                        endif
378                                        gCurDispFile = detStr
379
380                                        // now figure out q
381                                        // calculate the q-values, will be different depending on which panel is up (pixel size, geometry, etc.)
382                                        // TODO: !!!! get rid of the hard-wired values
383                                        // TODO: be sure that the units from HDF are what I expect
384                                        // TODO: beam center XY are pixels in the file, expected in the function, but are better suited for mm or cm
385                                        // TODO: units of xy pixel size are likely wrong
386//                                      xctr = V_getDet_beam_center_x(gCurDispType,detStr)              //written in pixels
387//                                      yctr = V_getDet_beam_center_y(gCurDispType,detStr)
388                                        xctr = V_getDet_beam_center_x_mm(gCurDispType,detStr)           //written in mm
389                                        yctr = V_getDet_beam_center_y_mm(gCurDispType,detStr)   
390                                       
391                                        sdd = V_getDet_ActualDistance(gCurDispType,detStr)      / 100   //written in cm, pass in meters
392                                        lam = V_getVSWavelength(gCurDispType)           //A
393//                                      pixSizeX = V_getDet_x_pixel_size(gCurDispType,detStr)           // written mm? need mm
394//                                      pixSizeY = V_getDet_y_pixel_size(gCurDispType,detStr)           // written mm? need mm
395//
396// TODO: these q-values ignore the non-linear corrections!!!
397// -- What can I do about this?
398                                        String destPath = "root:Packages:NIST:VSANS:"+gCurDispType
399                                        Wave data_realDistX = $(destPath + ":entry:instrument:detector_"+detStr+":data_realDistX")
400                                        Wave data_realDistY = $(destPath + ":entry:instrument:detector_"+detStr+":data_realDistY")     
401                                       
402// TODO: figure out what coordinates I need to pass -- xloc, yloc, textX, testY, (+1 on any?)                           
403                                        //gQQ = VC_CalcQval(testX,testY,xctr,yctr,sdd,lam,pixSizeX,pixSizeY)
404                                        //gQX = VC_CalcQX(testX,testY,xctr,yctr,sdd,lam,pixSizeX,pixSizeY)
405                                        //gQY = VC_CalcQY(testX,testY,xctr,yctr,sdd,lam,pixSizeX,pixSizeY)
406                                        gQQ = V_CalcQval(testX,testY,xctr,yctr,sdd,lam,data_realDistX,data_realDistY)
407                                        gQX = V_CalcQX(testX,testY,xctr,yctr,sdd,lam,data_realDistX,data_realDistY)
408                                        gQY = V_CalcQY(testX,testY,xctr,yctr,sdd,lam,data_realDistX,data_realDistY)
409
410                                        ii = -1         //look no further, set ii to bad value to exit the for loop
411                                endif   //end if(mouse is over a detector panel)
412                        endfor          // end loop over list of displayed images
413               
414                        break
415                       
416                        case 11: // keyboard event
417                                // TODO -- figure out why I'm not getting the TAB keystroke
418                                //  -- I want to be able to use the tab to change the focus only between File <.> and Tab > buttons, not everything
419                                // see the help file section "Keyboard Events" for an example and "WMWinHookStruct"
420                               
421                                //Print "key code = ",s.specialKeyCode
422                                //hookresult = 1                //if non-zero, we handled it and Igor will ignore it
423                                break
424                // And so on . . .
425        endswitch
426
427        return hookResult               // 0 if nothing done, else 1
428End
429
430
431// ********
432//
433// this procedure does most of the work for drawing the panels, setting the proper log/lin
434// scaling, the color scale, and the location based on the active tab
435//
436//lots to to here:
437//
438// - 1 - display the appropriate controls for each tab, and hide the others
439// - 2 - display the correct detector data for each tab, and remove the others from the graph
440// -----?? can I draw 3 graphs, and just put the right one on top?? move the other two to the side?
441//
442//
443// TODO
444//  -- add all of the controls of the VCALC panel (log scaling, adjusting the axes, etc.)
445//  x- get the panel to be correctly populated first, rather than needing to click everywhere to fill in
446//  x- remove the dependency on VCALC being initialized first, and using dummy waves from there...
447//
448//
449Function V_DataTabProc(tca) : TabControl
450        STRUCT WMTabControlAction &tca
451
452        switch( tca.eventCode )
453                case 2: // mouse up
454                        Variable tab = tca.tab
455               
456//                      SetDataFolder root:Packages:NIST:VSANS:VCALC
457                        SetDataFolder root:
458                        RemoveFromGraph/Z /W=VSANS_Data#det_panelsB tmp_asdf
459                        RemoveFromGraph/Z /W=VSANS_Data#det_panelsM tmp_asdf
460                        RemoveFromGraph/Z /W=VSANS_Data#det_panelsF tmp_asdf
461                        SetDataFolder root:
462                       
463                        SVAR dataType = root:Packages:NIST:VSANS:Globals:gCurDispType
464                       
465// make sure log scaling is correct
466                        NVAR state = root:Packages:NIST:VSANS:Globals:gIsLogScale
467                        if(State == 0)
468                                // lookup wave
469                                Wave LookupWave = root:Packages:NIST:VSANS:Globals:linearLookupWave
470                        else
471                                // lookup wave - the linear version
472                                Wave LookupWave = root:Packages:NIST:VSANS:Globals:logLookupWave
473                        endif
474                                               
475                       
476                        //************
477                        // -- can I use "ReplaceWave/W=VSANS_Data#det_panelsB allinCDF" to do this?
478                        // -- only works for "B", since for M and F panels, all 4 data sets are named "data"
479                        // in their respective folders...
480                       
481                       
482                        // get the slider values for the color mapping
483                        Variable lo,hi,lo_B,hi_B
484                        Variable lo_MT,lo_MB,lo_MR,lo_ML
485                        Variable lo_FT,lo_FB,lo_FR,lo_FL
486                        Variable hi_MT,hi_MB,hi_MR,hi_ML
487                        Variable hi_FT,hi_FB,hi_FR,hi_FL
488                        ControlInfo slider_lo
489                        lo = V_Value
490                        ControlInfo slider_hi
491                        hi = V_Value
492                       
493                       
494                        String tmpStr
495                        Variable ii
496                        if(tab==2)
497                                tmpStr = ImageNameList("VSANS_Data#det_panelsB",";")
498                                // for some odd reason, it appears that I need to work from the back of the list
499                                // since the traces get "renumbered" as I take them off !! A do loop may be a better choice
500                                if(ItemsInList(tmpStr) > 0)
501                                        do
502                                                RemoveImage /W=VSANS_Data#det_panelsB $(StringFromList(0,tmpStr,";"))           //get 1st item
503                                                tmpStr = ImageNameList("VSANS_Data#det_panelsB",";")                                                            //refresh list
504                                        while(ItemsInList(tmpStr) > 0)
505                                endif
506                               
507                                SetDataFolder $("root:Packages:NIST:VSANS:"+dataType+":entry:instrument:detector_B")
508                                Wave det_B=data
509                               
510                                CheckDisplayed /W=VSANS_Data#det_panelsB det_B 
511                                if(V_flag == 0)         // 0 == data is not displayed, so append it
512                                        AppendImage/W=VSANS_Data#det_panelsB det_B
513                                        lo_B = lo*(WaveMax(det_B) - WaveMin(det_B)) + WaveMin(det_B)
514                                        hi_B = hi*(WaveMax(det_B) - WaveMin(det_B)) + WaveMin(det_B)
515                                        ModifyImage/W=VSANS_Data#det_panelsB ''#0 ctab= {lo_B,hi_B,ColdWarm,0}          // don't autoscale {*,*,ColdWarm,0}
516                                endif
517                                MoveSubWindow/W=VSANS_Data#det_panelsB fnum=(50,185,517,620)
518                                MoveSubWindow/W=VSANS_Data#det_panelsM fnum=(320,70,430,160)
519                                MoveSubWindow/W=VSANS_Data#det_panelsF fnum=(200,70,310,160)
520
521                                ModifyImage/W=VSANS_Data#det_panelsB ''#0 ctabAutoscale=0,lookup= LookupWave
522                               
523                                // make the plot square
524                                ModifyGraph/W=VSANS_Data#det_panelsB width={Aspect,1}
525                               
526                                SetActiveSubWindow VSANS_Data#det_panelsB
527                                SetDataFolder root:
528                        endif
529       
530                        if(tab==1)
531                                tmpStr = ImageNameList("VSANS_Data#det_panelsM",";")
532                                // for some odd reason, it appears that I need to work from the back of the list
533                                // since the traces get "renumbered" as I take them off !! A do loop may be a better choice
534                                if(ItemsInList(tmpStr) > 0)
535                                        do
536                                                RemoveImage /W=VSANS_Data#det_panelsM $(StringFromList(0,tmpStr,";"))           //get 1st item
537                                                tmpStr = ImageNameList("VSANS_Data#det_panelsM",";")                                                            //refresh list
538                                        while(ItemsInList(tmpStr) > 0)
539                                endif
540
541                                SetDataFolder $("root:Packages:NIST:VSANS:"+dataType+":entry:instrument:detector_ML")
542                                Wave det_ML=data                               
543                                SetDataFolder $("root:Packages:NIST:VSANS:"+dataType+":entry:instrument:detector_MR")
544                                Wave det_MR=data
545                                SetDataFolder $("root:Packages:NIST:VSANS:"+dataType+":entry:instrument:detector_MT")
546                                Wave det_MT=data
547                                SetDataFolder $("root:Packages:NIST:VSANS:"+dataType+":entry:instrument:detector_MB")
548                                Wave det_MB=data
549
550                                CheckDisplayed /W=VSANS_Data#det_panelsM det_MR
551                                if(V_flag == 0)
552                                        AppendImage/W=VSANS_Data#det_panelsM det_MT             //order is important here to get LR on "top" of display
553                                        AppendImage/W=VSANS_Data#det_panelsM det_MB
554                                        AppendImage/W=VSANS_Data#det_panelsM det_ML
555                                        AppendImage/W=VSANS_Data#det_panelsM det_MR
556                                        lo_MT = lo*(WaveMax(det_MT) - WaveMin(det_MT)) + WaveMin(det_MT)
557                                        hi_MT = hi*(WaveMax(det_MT) - WaveMin(det_MT)) + WaveMin(det_MT)
558                                        lo_MB = lo*(WaveMax(det_MB) - WaveMin(det_MB)) + WaveMin(det_MB)
559                                        hi_MB = hi*(WaveMax(det_MB) - WaveMin(det_MB)) + WaveMin(det_MB)
560                                        lo_ML = lo*(WaveMax(det_ML) - WaveMin(det_ML)) + WaveMin(det_ML)
561                                        hi_ML = hi*(WaveMax(det_ML) - WaveMin(det_ML)) + WaveMin(det_ML)
562                                        lo_MR = lo*(WaveMax(det_MR) - WaveMin(det_MR)) + WaveMin(det_MR)
563                                        hi_MR = hi*(WaveMax(det_MR) - WaveMin(det_MR)) + WaveMin(det_MR)
564                                       
565                                        ModifyImage/W=VSANS_Data#det_panelsM ''#0 ctab= {lo_MT,hi_MT,ColdWarm,0}                // ''#n means act on the nth image (there are 4)
566                                        ModifyImage/W=VSANS_Data#det_panelsM ''#1 ctab= {lo_MB,hi_MB,ColdWarm,0}
567                                        ModifyImage/W=VSANS_Data#det_panelsM ''#2 ctab= {lo_ML,hi_ML,ColdWarm,0}
568                                        ModifyImage/W=VSANS_Data#det_panelsM ''#3 ctab= {lo_MR,hi_MR,ColdWarm,0}
569                                endif
570                                MoveSubWindow/W=VSANS_Data#det_panelsM fnum=(50,185,517,620)
571                                MoveSubWindow/W=VSANS_Data#det_panelsB fnum=(440,70,550,160)
572                                MoveSubWindow/W=VSANS_Data#det_panelsF fnum=(200,70,310,160)
573
574                                ModifyImage/W=VSANS_Data#det_panelsM ''#0 ctabAutoscale=0,lookup= LookupWave
575                                ModifyImage/W=VSANS_Data#det_panelsM ''#1 ctabAutoscale=0,lookup= LookupWave
576                                ModifyImage/W=VSANS_Data#det_panelsM ''#2 ctabAutoscale=0,lookup= LookupWave
577                                ModifyImage/W=VSANS_Data#det_panelsM ''#3 ctabAutoscale=0,lookup= LookupWave
578                               
579                                // make the plot square
580                                ModifyGraph/W=VSANS_Data#det_panelsM width={Aspect,1}
581                                                       
582                                SetActiveSubWindow VSANS_Data#det_panelsM
583                                SetDataFolder root:
584                        endif
585
586                        if(tab==0)
587                                tmpStr = ImageNameList("VSANS_Data#det_panelsF",";")
588                                // for some odd reason, it appears that I need to work from the back of the list
589                                // since the traces get "renumbered" as I take them off !! A do loop may be a better choice
590                                if(ItemsInList(tmpStr) > 0)
591                                        do
592                                                RemoveImage /W=VSANS_Data#det_panelsF $(StringFromList(0,tmpStr,";"))           //get 1st item
593                                                tmpStr = ImageNameList("VSANS_Data#det_panelsF",";")                                                            //refresh list
594                                        while(ItemsInList(tmpStr) > 0)
595                                endif
596
597                                SetDataFolder $("root:Packages:NIST:VSANS:"+dataType+":entry:instrument:detector_FL")
598                                Wave det_FL=data
599                                SetDataFolder $("root:Packages:NIST:VSANS:"+dataType+":entry:instrument:detector_FR")
600                                Wave det_FR=data
601                                SetDataFolder $("root:Packages:NIST:VSANS:"+dataType+":entry:instrument:detector_FT")
602                                Wave det_FT=data                               
603                                SetDataFolder $("root:Packages:NIST:VSANS:"+dataType+":entry:instrument:detector_FB")
604                                Wave det_FB=data
605                                                               
606                                CheckDisplayed /W=VSANS_Data#det_panelsF det_FL
607                                if(V_flag == 0)
608                                        AppendImage/W=VSANS_Data#det_panelsF det_FT
609                                        AppendImage/W=VSANS_Data#det_panelsF det_FB
610                                        AppendImage/W=VSANS_Data#det_panelsF det_FL
611                                        AppendImage/W=VSANS_Data#det_panelsF det_FR
612                                        lo_FT = lo*(WaveMax(det_FT) - WaveMin(det_FT)) + WaveMin(det_FT)
613                                        hi_FT = hi*(WaveMax(det_FT) - WaveMin(det_FT)) + WaveMin(det_FT)
614                                        lo_FB = lo*(WaveMax(det_FB) - WaveMin(det_FB)) + WaveMin(det_FB)
615                                        hi_FB = hi*(WaveMax(det_FB) - WaveMin(det_FB)) + WaveMin(det_FB)
616                                        lo_FL = lo*(WaveMax(det_FL) - WaveMin(det_FL)) + WaveMin(det_FL)
617                                        hi_FL = hi*(WaveMax(det_FL) - WaveMin(det_FL)) + WaveMin(det_FL)
618                                        lo_FR = lo*(WaveMax(det_FR) - WaveMin(det_FR)) + WaveMin(det_FR)
619                                        hi_FR = hi*(WaveMax(det_FR) - WaveMin(det_FR)) + WaveMin(det_FR)
620                                       
621                                        ModifyImage/W=VSANS_Data#det_panelsF ''#0 ctab= {lo_FT,hi_FT,ColdWarm,0}
622                                        ModifyImage/W=VSANS_Data#det_panelsF ''#1 ctab= {lo_FB,hi_FB,ColdWarm,0}
623                                        ModifyImage/W=VSANS_Data#det_panelsF ''#2 ctab= {lo_FL,hi_FL,ColdWarm,0}
624                                        ModifyImage/W=VSANS_Data#det_panelsF ''#3 ctab= {lo_FR,hi_FR,ColdWarm,0}
625                                endif
626                                MoveSubWindow/W=VSANS_Data#det_panelsF fnum=(50,185,517,620)
627                                MoveSubWindow/W=VSANS_Data#det_panelsB fnum=(440,70,550,160)
628                                MoveSubWindow/W=VSANS_Data#det_panelsM fnum=(320,70,430,160)
629                               
630                                ModifyImage/W=VSANS_Data#det_panelsF ''#0 ctabAutoscale=0,lookup= LookupWave
631                                ModifyImage/W=VSANS_Data#det_panelsF ''#1 ctabAutoscale=0,lookup= LookupWave
632                                ModifyImage/W=VSANS_Data#det_panelsF ''#2 ctabAutoscale=0,lookup= LookupWave
633                                ModifyImage/W=VSANS_Data#det_panelsF ''#3 ctabAutoscale=0,lookup= LookupWave
634
635                                // make the plot square
636                                ModifyGraph/W=VSANS_Data#det_panelsF width={Aspect,1}                           
637       
638                                SetActiveSubWindow VSANS_Data#det_panelsF
639                                SetDataFolder root:
640                        endif
641
642                                               
643                        break
644                case -1: // control being killed
645                        break
646        endswitch
647
648        return 0
649End
650
651// fake restore panels button click
652Function V_FakeRestorePanelsButtonClick()
653
654        STRUCT WMButtonAction ba
655        ba.eventCode = 2
656        V_RestorePanelButtonProc(ba)
657       
658        return(0)
659End
660
661
662// fake status button click
663Function V_FakeStatusButtonClick()
664
665        STRUCT WMButtonAction ba
666        ba.eventCode = 2
667        V_StatusButtonProc(ba)
668       
669        return(0)
670End
671
672// fake click on each tab to populate the data
673Function V_FakeTabClick(tab)
674        Variable tab
675       
676        STRUCT WMTabControlAction tca
677
678        tca.eventCode = 2               //fake mouse up
679        tca.tab = tab
680        V_DataTabProc(tca)
681       
682        TabControl tab0,win=VSANS_Data,value= tab               //select the proper tab
683        return(0)
684End
685
686// TODO
687//
688// move one file number back
689//
690Function V_File_minus_ButtonProc(ba) : ButtonControl
691        STRUCT WMButtonAction &ba
692
693        switch( ba.eventCode )
694                case 2: // mouse up
695                        // click code here
696                        V_LoadPlotAndDisplayRAW(-1)
697                       
698                        break
699                case -1: // control being killed
700                        break
701        endswitch
702
703        return 0
704End
705
706// TODO
707//
708// move one file number forward
709//
710Function V_File_plus_ButtonProc(ba) : ButtonControl
711        STRUCT WMButtonAction &ba
712
713        switch( ba.eventCode )
714                case 2: // mouse up
715                        // click code here
716                        V_LoadPlotAndDisplayRAW(1)
717
718                        break
719                case -1: // control being killed
720                        break
721        endswitch
722
723        return 0
724End
725
726//
727// button that mimics a click on the tab, cycling through the tabs 0->1->2->0 etc.
728// only goes one direction
729//
730Function V_Tab_p_ButtonProc(ba) : ButtonControl
731        STRUCT WMButtonAction &ba
732
733        switch( ba.eventCode )
734                case 2: // mouse up
735                        // click code here
736                        ControlInfo/W=VSANS_Data tab0
737                       
738                        V_Value += 1
739                        if(V_Value == 3)
740                                V_Value = 0             //reset to 0
741                        endif
742                        V_FakeTabClick(V_Value)
743
744// update the status when the tab is clicked                   
745                        STRUCT WMButtonAction sa
746                        sa.eventCode = 2
747                        V_StatusButtonProc(sa)
748
749                        break
750                case -1: // control being killed
751                        break
752        endswitch
753
754        return 0
755End
756
757
758// See V_Detector_Isolate.ipf
759// isolates a single panel to allow a better view of the details
760// useful for T/B panels which are partially blocked from view
761//
762// will open a separate panel to display the selected detector
763// (more to do here, depending what is necessary for instrument troubleshooting)
764// - like being able to turn corrections on/off and view with different axes (pix, mm, Q)
765//
766Function V_IsolateButtonProc(ba) : ButtonControl
767        STRUCT WMButtonAction &ba
768
769        switch( ba.eventCode )
770                case 2: // mouse up
771                        // click code here
772                        V_DetectorIsolate()
773                        break
774                case -1: // control being killed
775                        break
776        endswitch
777
778        return 0
779End
780
781//
782// if the data display is RAW, convert to the specified WORK data type
783//
784// TODO
785// -- better error checking
786// -- if the data type is not RAW, can I Copy Folder instead?
787//
788Function V_ToWorkFileButtonProc(ba) : ButtonControl
789        STRUCT WMButtonAction &ba
790
791        switch( ba.eventCode )
792                case 2: // mouse up
793                        // click code here
794                        //Convert_to_Workfile(newtype, doadd) // a proc
795                        Execute "V_Convert_to_Workfile()"
796
797                        break
798                case -1: // control being killed
799                        break
800        endswitch
801
802        return 0
803End
804
805
806// TODO
807//
808// opens a separate panel with the I(q) representation of the data
809// ? controls here to select how the data is processed/grouped/saved, etc.
810//
811// -- currently just the graph, no controls
812// -- this re-bins the data every time by calling V_QBinAllPanels(folderStr)
813Function V_IvsQPanelButtonProc(ba) : ButtonControl
814        STRUCT WMButtonAction &ba
815
816        switch( ba.eventCode )
817                case 2: // mouse up
818                        // click code here
819                       
820                        V_PlotData_Panel()
821                       
822                        break
823                case -1: // control being killed
824                        break
825        endswitch
826
827        return 0
828End
829
830// TODO
831//
832// gets the status of the currently displayed file and dumps it to the panel (not the cmd window)
833// - lots to decide here about what is the important stuff to display. There's a lot more information for VSANS
834//
835Function V_StatusButtonProc(ba) : ButtonControl
836        STRUCT WMButtonAction &ba
837
838        switch( ba.eventCode )
839                case 2: // mouse up
840                        // click code here
841                       
842                        // figure out wether to display per carraige, or the whole file
843                        SVAR str = root:Packages:NIST:VSANS:Globals:gStatusText
844                        SVAR type = root:Packages:NIST:VSANS:Globals:gCurDispType               //what folder
845                        // which tab active
846                        ControlInfo/W=VSANS_Data tab0
847                        Variable curTab = V_value
848                       
849                        str = "Current data is from "+ type + "\r"
850                        str += "Description = "+V_getSampleDescription(type) + "\r"
851                        str += "Wavelength is "+num2str(V_getWavelength(type)) + " A \r"
852                        if(curTab == 2)
853                                str += "SDD B = "+num2str(V_getDet_ActualDistance(type,"B")) + " cm \r"         //V_getDet_distance(fname,detStr)
854                        endif
855                        if(curTab == 1)
856                                str += "SDD ML = "+num2str(V_getDet_ActualDistance(type,"ML")) + " cm \r"
857                                str += "SDD MR = "+num2str(V_getDet_ActualDistance(type,"MR")) + " cm \r"
858                                str += "SDD MT = "+num2str(V_getDet_ActualDistance(type,"MT")) + " cm \r"
859                                str += "SDD MB = "+num2str(V_getDet_ActualDistance(type,"MB")) + " cm \r"
860                        endif
861                        if(curTab == 0)
862                                str += "SDD FL = "+num2str(V_getDet_ActualDistance(type,"FL")) + " cm \r"
863                                str += "SDD FR = "+num2str(V_getDet_ActualDistance(type,"FR")) + " cm \r"
864                                str += "SDD FT = "+num2str(V_getDet_ActualDistance(type,"FT")) + " cm \r"
865                                str += "SDD FB = "+num2str(V_getDet_ActualDistance(type,"FB")) + " cm \r"
866                        endif
867                       
868                       
869                       
870                       
871                       
872                        break
873                case -1: // control being killed
874                        break
875        endswitch
876
877        return 0
878End
879
880
881// TODO:
882// x- link this to the preferences for the display. this is done in UpdateDisplayInformation (the main call) so that
883//     the panels are rescaled only once, rather than toggled three times (F, M, B) if I call from the tabProc
884// -- come up with a better definition of the log lookup wave (> 1000 pts, what is the first point)
885// -- make an equivalent linear wave
886// -- hard wire it in so it is created at initialization and stored someplace safe
887// -- catch the error if it doesn't exist (re-make the wave as needed)
888//
889// Using the ModifyImage log=1 keyword fails for values of zero in the data, which is a common
890// occurrence with count data. the display just goes all gray, so that's not an option. Use the lookup wave instead
891//
892// toggle the (z) value of the display log/lin
893//
894Function V_LogLinButtonProc(ba) : ButtonControl
895        STRUCT WMButtonAction &ba
896
897        switch( ba.eventCode )
898                case 2: // mouse up
899                        // click code here
900                       
901                        // which tab is active? does it matter - or do I log-scale everything?
902                        // log/lin current state is in the S_UserData string (0=linear, 1=log)
903                        ControlInfo/W=VSANS_Data button_log
904                        Variable curState,newState
905                        String newStateStr,newTitleStr
906                       
907                       
908                        curState = str2num(S_UserData)
909                       
910                        if(curState == 0)
911                                newState = 1
912                                newStateStr="1"
913                                newTitleStr = "isLog"
914                                // lookup wave
915                                Wave LookupWave = root:Packages:NIST:VSANS:Globals:logLookupWave
916                        else
917                                newState = 0
918                                newStateStr="0"
919                                newTitleStr = "isLin"
920                                // lookup wave - the linear version
921                                Wave LookupWave = root:Packages:NIST:VSANS:Globals:linearLookupWave
922                        endif
923                       
924                        // update the button and the global value
925                        Button button_log,userData=newStateStr,title=newTitleStr
926                        NVAR state = root:Packages:NIST:VSANS:Globals:gIsLogScale
927                        state = newState
928                       
929                        // on the front:                       
930                        ModifyImage/W=VSANS_Data#det_panelsF ''#0 ctabAutoscale=0,lookup= LookupWave
931                        ModifyImage/W=VSANS_Data#det_panelsF ''#1 ctabAutoscale=0,lookup= LookupWave
932                        ModifyImage/W=VSANS_Data#det_panelsF ''#2 ctabAutoscale=0,lookup= LookupWave
933                        ModifyImage/W=VSANS_Data#det_panelsF ''#3 ctabAutoscale=0,lookup= LookupWave
934                        //on the middle:
935//                      ModifyImage/W=VSANS_Data#det_panelsM ''#0 log=newState
936//                      ModifyImage/W=VSANS_Data#det_panelsM ''#1 log=newState
937//                      ModifyImage/W=VSANS_Data#det_panelsM ''#2 log=newState
938//                      ModifyImage/W=VSANS_Data#det_panelsM ''#3 log=newState
939                               
940                        ModifyImage/W=VSANS_Data#det_panelsM ''#0 ctabAutoscale=0,lookup= LookupWave
941                        ModifyImage/W=VSANS_Data#det_panelsM ''#1 ctabAutoscale=0,lookup= LookupWave
942                        ModifyImage/W=VSANS_Data#det_panelsM ''#2 ctabAutoscale=0,lookup= LookupWave
943                        ModifyImage/W=VSANS_Data#det_panelsM ''#3 ctabAutoscale=0,lookup= LookupWave
944
945
946                        // on the back:
947                        ModifyImage/W=VSANS_Data#det_panelsB ''#0 ctabAutoscale=0,lookup= LookupWave
948
949                        break
950                case -1: // control being killed
951                        break
952        endswitch
953
954        return 0
955End
956
957// TODO
958// possibly function to "tag" files right here in the disaply with things
959// like their intent, or other values that reduction will need, kind of like a "quick patch"
960// with limited functionality (since full function would be a nightmare!)
961Function V_TagFileButtonProc(ba) : ButtonControl
962        STRUCT WMButtonAction &ba
963
964        switch( ba.eventCode )
965                case 2: // mouse up
966                        // click code here
967                        DoAlert 0, "TagFileButtonProc(ba) unfinished - this may be used to 'tag' a file as scatter, trans, emp, bkg, etc."
968                       
969                               
970                        break
971                case -1: // control being killed
972                        break
973        endswitch
974
975        return 0
976End
977
978// TODO
979// -- fill in more functionality
980// -- currently a straight concatentation of all data, no options
981// -- maybe allow save of single panels?
982// -- any other options?
983Function V_SaveIQ_ButtonProc(ba) : ButtonControl
984        STRUCT WMButtonAction &ba
985
986        switch( ba.eventCode )
987                case 2: // mouse up
988                        // click code here
989
990                        Execute "V_Combine1DData()"                     
991                               
992                        break
993                case -1: // control being killed
994                        break
995        endswitch
996
997        return 0
998End
999
1000//TODO
1001//
1002//link this to the beam center finding panel
1003//
1004Function V_BeamCtrButtonProc(ba) : ButtonControl
1005        STRUCT WMButtonAction &ba
1006
1007        switch( ba.eventCode )
1008                case 2: // mouse up
1009                        // click code here
1010                        V_FindBeamCenter()
1011                        break
1012                case -1: // control being killed
1013                        break
1014        endswitch
1015
1016        return 0
1017End
1018
1019//TODO
1020//
1021// this "spreads" the display of panels to a nominal separation for easier viewing
1022//
1023Function V_SpreadPanelButtonProc(ba) : ButtonControl
1024        STRUCT WMButtonAction &ba
1025
1026        switch( ba.eventCode )
1027                case 2: // mouse up
1028                        // click code here
1029                        V_SpreadOutPanels()
1030                       
1031                        break
1032                case -1: // control being killed
1033                        break
1034        endswitch
1035
1036        return 0
1037End
1038
1039//TODO
1040//
1041// this "restores" the display of panels to their actual position based on the apparent beam center
1042//
1043Function V_RestorePanelButtonProc(ba) : ButtonControl
1044        STRUCT WMButtonAction &ba
1045
1046        switch( ba.eventCode )
1047                case 2: // mouse up
1048                        // click code here
1049                        V_RestorePanels()
1050                               
1051                        break
1052                case -1: // control being killed
1053                        break
1054        endswitch
1055
1056        return 0
1057End
1058
1059// TODO
1060//
1061// link this slider to the "high" end of the color mapping for whatever is currently displayed
1062//
1063// -- see Buttons.ipf for the old SANS implementation
1064//
1065Function V_HiMapSliderProc(sa) : SliderControl
1066        STRUCT WMSliderAction &sa
1067
1068        switch( sa.eventCode )
1069                case -1: // control being killed
1070                        break
1071                default:
1072                        if( sa.eventCode & 1 ) // value set
1073                                Variable curval = sa.curval
1074                                ControlInfo tab0
1075                                V_FakeTabClick(V_Value)
1076                               
1077//                              ControlInfo slider_lo
1078//                              V_MakeImageLookupTables(10000,V_Value,curval)
1079                        endif
1080                        break
1081        endswitch
1082
1083        return 0
1084End
1085
1086// TODO
1087//
1088// link this slider to the "low" end of the color mapping for whatever is currently displayed
1089//
1090// -- see Buttons.ipf for the old SANS implementation
1091//
1092Function V_LowMapSliderProc(sa) : SliderControl
1093        STRUCT WMSliderAction &sa
1094
1095        switch( sa.eventCode )
1096                case -1: // control being killed
1097                        break
1098                default:
1099                        if( sa.eventCode & 1 ) // value set
1100                                Variable curval = sa.curval
1101                                ControlInfo tab0
1102                                V_FakeTabClick(V_Value)
1103
1104//                              ControlInfo slider_hi
1105//                              V_MakeImageLookupTables(10000,curval,V_Value)
1106                        endif
1107                        break
1108        endswitch
1109
1110        return 0
1111End
1112
1113
Note: See TracBrowser for help on using the repository browser.