Context Navigation

DataSetHandling.ipf @ 1239

Last change on this file since 1239 was 820, checked in by ajj, 12 years ago
Modified ReSortDataSet?() to sort resolution data, handle 6 vs 3 column data and ignore USANS data Added button to dataset management panel to do resort
File size: 58.5 KB

Line
1	#pragma rtGlobals=1 // Use modern global access method.
2
3
4	///////// SRK - VERY SIMPLE batch converter has been added: see
5	//
6	// Function batchXML26ColConvert()
7	//
8	// Function batchGrasp26ColConvert()
9	//
10	// these both need a real interface, and a way to better define the name of the
11	// converted output file. And some retention of the header would be nice too...
12	//
13
14
15	// Functions and interfaces to manage datasets now that they are in data folders
16	// Planned interface
17	// - Panel to select/load data and then select operation
18	// Planned functions:
19	// - Rename a data set - AJJ Done Nov 2009
20	// - Duplicate a data set - AJJ Done Nov 2009
21	// - Subtract one data set from another
22	// - Add one data set to another
23	// - Divide data sets
24	// - Multiply data sets
25	// - Save a data folder to file
26
27	/////////////////// Data Management Panel ///////////////////////////////////////////////////////
28
29	//
30	Proc MakeDMPanel()
31	DoWindow/F DataManagementPanel
32	if(V_flag==0)
33	fMakeDMPanel()
34	endif
35	End
36
37	Proc fMakeDMPanel()
38	PauseUpdate; Silent 1 // building window...
39	NewPanel /W=(459,44,959,460)/N=DataManagementPanel/K=2 as "Data Set Management"
40	ModifyPanel fixedSize=1,cbRGB=(30000,60000,60000)
41
42	//Main bit of panel
43	GroupBox grpBox_0,pos={20,10},size={460,150}
44	GroupBox grpBox_1,pos={20,180},size={460,70}
45	GroupBox grpBox_2,pos={20,270},size={460,40}
46
47	GroupBox grpBox_3,pos={20,330},size={460,40}
48
49	Button DS_button,title="Load 1D Data Set",pos={300,20},size={150,20}
50	Button DS_button,proc=DM_LoadDataSetProc
51	Button Unload_button,title="Unload 1D Data Set",pos={300,50},size={150,20}
52	Button Unload_button,proc=DM_UnloadProc
53	Button Save_button,title="Save 1D Data Set",pos={300,80},size={150,20}
54	Button Save_button,proc=DM_SaveProc
55	Button ReSort_button,title="Re-Sort 1D Data Set",pos={300,130},size={150,20}
56	Button ReSort_button,proc=DM_ReSortProc
57	PopupMenu DS_popup,pos={30,40},size={318,20},title="Data Set ",proc=DM_PopupProc
58	PopupMenu DS_popup,mode=1,value= #"DM_DataSetPopupList()"
59
60	CheckBox XMLStateCtrl,pos={30,82},size={124,14},title="XML Output Enabled (change in Preferences)"
61	CheckBox XMLStateCtrl,help={"Default output format is canSAS XML rather than NIST 6 column"}
62	CheckBox XMLStateCtrl,value= root:Packages:NIST:gXML_Write,disable=2
63
64	Button Rename_button,title="Rename",pos={75,220},size={150,20}
65	Button Rename_button,proc=DM_RenameProc
66	Button Duplicate_button,title="Duplicate",pos={275,220},size={150,20}
67	Button Duplicate_button,proc=DM_DuplicateProc
68
69	SetVariable NewName_setvar,title="New Name (max 25 characters)",pos={50,190},size={400,20}
70	SetVariable NewName_setvar,fsize=12,value=_STR:"",proc=DMNameSetvarproc,live=1
71
72	Button SaveAsXML_button,title="Save as canSAS XML",pos={75,280},size={150,20}
73	Button SaveAsXML_button,proc=DMSaveAsXMLproc
74
75	Button SaveAs6col_button,title="Save as NIST 6 column",pos={275,280},size={160,20}
76	Button SaveAs6col_button,proc=DMSaveAs6colproc
77
78	Button BatchConvertData_button,title="Batch Convert Format of 1D Data Files",pos={75,340},size={350,20}
79	Button BatchConvertData_button,proc=DMBatchConvertProc
80
81	Button DMDone_button,title="Done",pos={360,370},size={60,20}
82	Button DMDone_button,proc=DMDoneButtonProc
83	Button DMHelp_button,title="?",pos={440,370},size={30,20}
84	Button DMHelp_button,proc=DMHelpButtonProc
85
86
87	ControlInfo/W=DataManagementPanel DS_popup
88	if (cmpstr(S_Value,"No data loaded") == 0)
89	SetVariable NewName_setvar,value=_STR:"dataset_copy"
90	else
91	//fake call to popup
92	//// -- can't use STRUCT in a Proc, only a function
93	// STRUCT WMPopupAction pa
94	// pa.win = "DataManagementPanel"
95	// pa.ctrlName = "DS_popup"
96	// pa.eventCode = 2
97	// DM_PopupProc(pa)
98	endif
99
100	End
101
102	Function DMSaveAs6colproc(ba) : ButtonControl
103	STRUCT WMButtonAction &ba
104
105	switch( ba.eventCode)
106	case 2:
107	ControlInfo/W=DataManagementPanel DS_popup
108	String folderName = S_Value
109	fReWrite1DData(folderName,"tab","CRLF")
110	break
111	endswitch
112
113	return 0
114	End
115
116	Function DMSaveAsXMLproc(ba) : ButtonControl
117	STRUCT WMButtonAction &ba
118
119	switch( ba.eventCode)
120	case 2:
121	ControlInfo/W=DataManagementPanel DS_popup
122	String folderName = S_Value
123	ReWrite1DXMLData(folderName)
124	break
125	endswitch
126
127	return 0
128	End
129
130
131	//Function SaveDataSetToFile(folderName)
132	// String folderName
133	//
134	// String protoStr = ""
135	// //Check for history string in folder
136	// //If it doesn't exist then
137	//
138	// //Do saving of data file.
139	//
140	// NVAR gXML_Write = root:Packages:NIST:gXML_Write
141	//
142	// if (gXML_Write == 1)
143	// ReWrite1DXMLData(folderName)
144	// else
145	// fReWrite1DData(folderName,"tab","CRLF")
146	// endif
147	//
148	// //Include history string to record what was done?
149	//
150	//End
151
152
153
154
155
156
157	Function DMBatchConvertProc(ba) : ButtonControl
158	STRUCT WMButtonAction &ba
159
160	switch( ba.eventCode)
161	case 2:
162	Execute "MakeNCNRBatchConvertPanel()"
163	break
164	endswitch
165
166	return 0
167	End
168
169
170	Function DMNameSetvarproc(sva) : SetVariableControl
171	STRUCT WMSetVariableAction &sva
172
173	switch( sva.eventCode )
174	case 1: // mouse up
175	case 2: // Enter key
176	case 3: // Live update
177	String sv = sva.sval
178	if( strlen(sv) > 25 )
179	sv= sv[0,24]
180	SetVariable $(sva.ctrlName),win=$(sva.win),value=_STR:sv
181	ControlUpdate /W=$(sva.win) $(sva.ctrlName)
182	Beep
183	endif
184	break
185	endswitch
186	return 0
187	End
188
189	Function DM_RenameProc(ba) : ButtonControl
190	STRUCT WMButtonAction &ba
191
192	String DS,NewName
193
194	ControlInfo/W=$(ba.win) DS_popup
195	DS = S_Value
196
197	ControlInfo/W=$(ba.win) NewName_setvar
198	NewName = CleanupName(S_Value, 0 ) //clean up any bad characters, and put the cleaned string back
199	SetVariable NewName_setvar,value=_STR:NewName
200
201	switch (ba.eventcode)
202	case 2: // mouse up
203	RenameDataSet(DS,NewName)
204	ControlUpdate /W=$(ba.win) DS_Popup
205	break
206	endswitch
207
208
209
210	return 0
211	end
212
213	Function DM_DuplicateProc(ba) : ButtonControl
214	STRUCT WMButtonAction &ba
215
216	String DS,NewName
217
218	ControlInfo/W=$(ba.win) DS_popup
219	DS = S_Value
220
221	ControlInfo/W=$(ba.win) NewName_setvar
222	NewName = CleanupName(S_Value, 0 ) //clean up any bad characters, and put the cleaned string back
223	SetVariable NewName_setvar,value=_STR:NewName
224
225	switch (ba.eventcode)
226	case 2: // mouse up
227	DuplicateDataSet(DS,NewName,0)
228	ControlUpdate /W=$(ba.win) DS_Popup
229	break
230	endswitch
231
232	return 0
233	end
234
235	Function DM_SaveProc(ba) : ButtonControl
236	STRUCT WMButtonAction &ba
237
238	switch(ba.eventCode)
239	case 2:
240	ControlInfo/W=$(ba.win) DS_popup
241	SaveDataSetToFile(S_Value)
242	break
243	endswitch
244
245	return 0
246	end
247
248	Function DM_UnloadProc(ba) : ButtonControl
249	STRUCT WMButtonAction &ba
250
251	switch (ba.eventcode)
252	case 2: // mouse up
253
254	String savDF=GetDataFolder(1)
255	String DF
256	ControlInfo /W=$(ba.win) DS_Popup
257	DF = S_Value
258
259	print DF
260	//check for horrific null output from control
261	if (cmpstr(DF,"") != 0)
262
263	SetDataFolder DF
264	KillVariables/A
265	SetDataFolder savDF
266
267	KillDataFolder/Z $DF
268	ControlUpdate /W=$(ba.win) DS_Popup
269
270	ControlInfo/W=DataManagementPanel DS_popup
271	if (cmpstr(S_Value,"No data loaded") == 0)
272	SetVariable NewName_setvar,value=_STR:"dataset_copy"
273	else
274	//fake call to popup
275	STRUCT WMPopupAction pa
276	pa.win = "DataManagementPanel"
277	pa.ctrlName = "DS_popup"
278	pa.eventCode = 2
279	DM_PopupProc(pa)
280	endif
281	endif
282	break
283	endswitch
284
285	return 0
286	end
287
288
289	Function DM_PopupProc(pa) : PopupMenuControl
290	STRUCT WMPopupAction &pa
291
292	String resultName
293
294	switch( pa.eventCode)
295	case 2:
296	//print "Called by "+pa.ctrlname+" with value "+pa.popStr
297	ControlInfo/W=$(pa.win) $(pa.ctrlName)
298	String popStr = S_Value
299	if (stringmatch(pa.ctrlname,"DS") == 1)
300	resultName = stringfromlist(0,popStr,"_")+"_copy"
301
302	SetVariable NewName_setvar win=$(pa.win), value=_STR:resultName
303	endif
304	break
305	endswitch
306
307
308	End
309
310
311	//Must follow naming scheme to match buttons to popups
312	//"Name_button" goes with "Name_popup"
313	Function DM_LoadDataSetProc(ba) : ButtonControl
314	STRUCT WMButtonAction &ba
315
316
317	switch( ba.eventCode )
318	case 2: // mouse up
319	// click code here
320	String cmd = "A_LoadOneDDataWithName(\"\","+num2str(0)+")"
321	Execute cmd
322
323	SVAR gLastFileName = root:packages:NIST:gLastFileName
324
325	String windowName = ba.win
326	String popupName = StringFromList(0,ba.ctrlName,"_")+"_popup"
327
328	ControlUpdate/W=$(windowName) $(popupName)
329	//instead of a simple controlUpdate, better to pop the menu to make sure that the other menus follow
330	// convoluted method to find the right item and pop the menu.
331
332	String list,folderStr
333	Variable num
334	folderStr = CleanupName(gLastFileName,0)
335	list = DM_DataSetPopupList()
336	num=WhichListItem(folderStr,list,";",0,0)
337	if(num != -1)
338	PopupMenu $(popupName),mode=num+1,win=$(windowName)
339	ControlUpdate/W=$(windowName) $(popupName)
340
341	if (cmpstr(popupName,"DS_popup") == 0)
342	//send fake mouse action to popup to update old name if
343	Struct WMPopupAction pa
344	pa.eventCode = 2 //fake mouse up
345	pa.win = windowName
346	pa.ctrlName = "DS_popup"
347	DM_PopupProc(pa)
348	endif
349	endif
350	break
351	endswitch
352
353	return 0
354	End
355
356
357	Function/S DMGetDSName(dsNum)
358	Variable dsNum
359
360	String ctrlName
361	if (dsNum == 1)
362	ctrlName = "DS1_popup"
363	elseif (dsNum == 2)
364	ctrlName = "DS2_popup"
365	endif
366
367	ControlInfo/W=DataManagementPanel $(ctrlName)
368
369	Return S_Value
370
371	End
372
373
374	Function DMDoneButtonProc(ba) : ButtonControl
375	STRUCT WMButtonAction &ba
376
377	String win = ba.win
378
379	switch (ba.eventCode)
380	case 2:
381	DoWindow/K DataManagementPanel
382	break
383	endswitch
384
385	return 0
386	End
387
388	Function DMHelpButtonProc(ba) : ButtonControl
389	STRUCT WMButtonAction &ba
390
391	String win = ba.win
392
393	switch (ba.eventCode)
394	case 2:
395	// click code here
396	DisplayHelpTopic/Z/K=1 "Data Set Management"
397	if(V_flag !=0)
398	DoAlert 0,"The Data Set Management Help file could not be found"
399	endif
400	break
401	endswitch
402
403	return 0
404	End
405
406	Function DM_ReSortProc(ba) : ButtonControl
407	STRUCT WMButtonAction &ba
408
409	switch(ba.eventCode)
410	case 2:
411	ControlInfo/W=$(ba.win) DS_popup
412	ReSortDataSet(S_Value)
413	break
414	endswitch
415
416	return 0
417	end
418
419	/////////////////////// Batch Data Conversion Panel ////////////////////////////////////
420	//
421	//
422
423	Proc MakeNCNRBatchConvertPanel()
424	NCNRInitBatchConvert()
425	DoWindow/F NCNRBatchConvertPanel
426	if(V_flag==0)
427	fMakeNCNRBatchConvertPanel()
428	endif
429	End
430
431
432	Function NCNRInitBatchConvert()
433	NewDataFolder/O/S root:Packages:NIST:BatchConvert
434	Make/O/T/N=1 filewave=""
435	Make/O/N=1 selWave=0
436	Variable/G ind=0,gRadioVal=1
437	SetDataFolder root:
438	End
439
440	Proc fMakeNCNRBatchConvertPanel()
441	PauseUpdate; Silent 1 // building window...
442	NewPanel /W=(658,347,1018,737)/N=NCNRBatchConvertPanel/K=2 as "Batch Convert 1D Data"
443	// NewPanel /W=(658,347,1018,737)/N=NCNRBatchConvertPanel as "Batch Convert 1D Data"
444	ModifyPanel cbRGB=(40000,50000,32896)
445	ModifyPanel fixedSize=1
446
447	ListBox fileList,pos={13,11},size={206,179}
448	ListBox fileList,listWave=root:Packages:NIST:BatchConvert:fileWave
449	ListBox fileList,selWave=root:Packages:NIST:BatchConvert:selWave,mode= 4
450
451	Button button7,pos={238,20},size={100,20},proc=NCNRBatchConvertNewFolder,title="New Folder"
452	Button button7,help={"Select a new data folder"}
453	TitleBox msg0,pos={238,160},size={100,30},title="\JCShift-click to\rselect multiple files"
454	TitleBox msg0,frame=0,fixedSize=1
455
456	GroupBox filterGroup,pos={13,200},size={206,60},title="Filter list by input file type"
457	CheckBox filterCheck_1,pos={24,220},size={36,14},title="XML",value= 1,mode=1, proc=BC_filterCheckProc
458	CheckBox filterCheck_2,pos={24,239},size={69,14},title="ABS or AVE",value= 0,mode=1, proc=BC_filterCheckProc
459	CheckBox filterCheck_3,pos={100,220},size={69,14},title="none",value= 0,mode=1, proc=BC_filterCheckProc
460
461	Button button8,pos={238,75},size={100,20},proc=NCNRBatchConvertHelpProc,title="Help"
462	Button button9,pos={238,47},size={100,20},proc=NCNRBatchConvertRefresh,title="Refresh List"
463	Button button10,pos={238,102},size={100,20},proc=NCNRBatchConvertSelectAll,title="Select All"
464	Button button0,pos={238,130},size={100,20},proc=NCNRBatchConvertDone,title="Done"
465
466	GroupBox outputGroup,pos={13,270},size={206,60},title="Output File Type"
467	CheckBox outputCheck_1,pos={24,289},size={36,14},title="XML",value= 0,mode=1, proc=BC_outputCheckProc
468	CheckBox outputCheck_2,pos={24,309},size={69,14},title="ABS or AVE",value= 1,mode=1, proc=BC_outputCheckProc
469
470	Button button6,pos={13,350},size={206,20},proc=NCNRBatchConvertFiles,title="Convert File(s)"
471	Button button6,help={"Converts the files to the format selected"}
472
473
474
475	End
476
477	Function BC_filterCheckProc(ctrlName,checked) : CheckBoxControl
478	String ctrlName
479	Variable checked
480
481	NVAR gRadioVal= root:Packages:NIST:BatchConvert:gRadioVal
482
483	strswitch (ctrlName)
484	case "filterCheck_1":
485	gRadioVal= 1
486	break
487	case "filterCheck_2":
488	gRadioVal= 2
489	break
490	case "filterCheck_3":
491	gRadioVal= 3
492	break
493	endswitch
494	CheckBox filterCheck_1,value= gRadioVal==1
495	CheckBox filterCheck_2,value= gRadioVal==2
496	CheckBox filterCheck_3,value= gRadioVal==3
497
498	NCNRBatchConvertGetList()
499
500	return(0)
501	End
502
503	Function BC_outputCheckProc(ctrlName,checked) : CheckBoxControl
504	String ctrlName
505	Variable checked
506
507	if(cmpstr("outputCheck_1",ctrlName)==0)
508	CheckBox outputCheck_1,value=checked
509	CheckBox outputCheck_2,value=!checked
510	else
511	CheckBox outputCheck_1,value=!checked
512	CheckBox outputCheck_2,value=checked
513	endif
514
515	return(0)
516	End
517
518	Function NCNRBatchConvertFiles(ba) : ButtonControl
519	STRUCT WMButtonAction &ba
520
521
522	switch (ba.eventCode)
523	case 2:
524
525	//check the input/output as best I can (none may be the input filter)
526	Variable inputType,outputType=1
527	NVAR gRadioVal= root:Packages:NIST:BatchConvert:gRadioVal
528	inputType = gRadioVal
529	ControlInfo outputCheck_1
530	if(V_value==1)
531	outputType = 1 //xml
532	else
533	outputType = 2 //6-col
534	endif
535
536	if(inputType==outputType)
537	DoAlert 0,"Input and output types are the same. Nothing will be converted"
538	return(0)
539	endif
540
541
542	// input and output are different, proceed
543
544	Wave/T fileWave=$"root:Packages:NIST:BatchConvert:fileWave"
545	Wave sel=$"root:Packages:NIST:BatchConvert:selWave"
546
547	String fname="",pathStr="",newFileName=""
548	Variable ii,num
549	PathInfo catPathName //this is where the files are
550	pathStr=S_path
551
552	// process the selected items
553	num=numpnts(sel)
554	ii=0
555	do
556	if(sel[ii] == 1)
557	fname=pathStr + fileWave[ii]
558
559	if(outputType == 1)
560	convertNISTtoNISTXML(fname)
561	endif
562
563	if(outputType == 2)
564	convertNISTXMLtoNIST6Col(fname)
565	endif
566	endif
567	ii+=1
568	while(ii<num)
569
570	break
571	endswitch
572
573	return 0
574	End
575
576	Function NCNRBatchConvertSelectAll(ba) : ButtonControl
577	STRUCT WMButtonAction &ba
578
579	switch (ba.eventcode)
580	case 2:
581	Wave sel=$"root:Packages:NIST:BatchConvert:selWave"
582	sel = 1
583	break
584	endswitch
585
586	End
587
588	Function NCNRBatchConvertNewFolder(ba) : ButtonControl
589	STRUCT WMButtonAction &ba
590
591	switch (ba.eventcode)
592	case 2:
593	A_PickPath()
594	NCNRBatchConvertGetList()
595	break
596	endswitch
597
598	End
599
600
601	// filter is a bit harsh - will need to soften this by presenting an option to enter the suffix
602	//
603	Function NCNRBatchConvertGetList()
604
605	//make sure that path exists
606	PathInfo catPathName
607	if (V_flag == 0)
608	Abort "Folder path does not exist - use \"New Folder\" button"
609	Endif
610
611	String newList = A_ReducedDataFileList(""),tmpList=""
612	Variable num
613
614	NVAR gRadioVal= root:Packages:NIST:BatchConvert:gRadioVal
615	ControlInfo filterCheck_1
616	if(gRadioVal == 1)
617	//keep XML data
618	tmpList = ListMatch(newList, "*.ABSx" ,";") //note that ListMatch is not case-sensitive
619	tmpList += ListMatch(newList, "*.AVEx" ,";")
620	tmpList += ListMatch(newList, "*.CORx" ,";")
621	tmpList += ListMatch(newList, "*.xml" ,";")
622	else
623	if(gRadioVal ==2)
624	//keep ave, abs data
625	tmpList = ListMatch(newList, "*.ABS" ,";")
626	tmpList += ListMatch(newList, "*.AVE" ,";")
627	else
628	//return everything
629	tmpList = newList
630	endif
631	endif
632	newList = tmpList
633
634	num=ItemsInList(newlist,";")
635	WAVE/T fileWave=$"root:Packages:NIST:BatchConvert:fileWave"
636	WAVE selWave=$"root:Packages:NIST:BatchConvert:selWave"
637	Redimension/N=(num) fileWave
638	Redimension/N=(num) selWave
639	fileWave = StringFromList(p,newlist,";")
640	Sort filewave,filewave
641
642	return 0
643	End
644
645	Function NCNRBatchConvertRefresh(ba) : ButtonControl
646	STRUCT WMButtonAction &ba
647
648	switch (ba.eventCode)
649	case 2:
650	NCNRBatchConvertGetList()
651	break
652	endswitch
653
654	return 0
655	End
656
657
658	Function NCNRBatchConvertDone(ba) : ButtonControl
659	STRUCT WMButtonAction &ba
660
661	switch (ba.eventCode)
662	case 2:
663	DoWindow/K NCNRBatchConvertPanel
664	break
665	endswitch
666
667	return 0
668	End
669
670	Function NCNRBatchConvertHelpProc(ba) : ButtonControl
671	STRUCT WMButtonAction &ba
672
673	String win = ba.win
674
675	switch (ba.eventCode)
676	case 2:
677	// click code here
678	DisplayHelpTopic/Z/K=1 "Batch Data Conversion"
679	if(V_flag !=0)
680	DoAlert 0,"The Batch Data Conversion Help file could not be found"
681	endif
682	break
683	endswitch
684
685	return 0
686	End
687
688	/////////////////////// Data Arithmetic Panel /////////////////////////////////////////
689
690	//
691	Function MakeDAPanel()
692	DoWindow/F DataArithmeticPanel
693	if(V_flag==0)
694	fMakeDAPanel()
695	else
696	DoWindow/F DAPlotPanel
697	endif
698
699	return(0)
700	End
701
702	Function fMakeDAPanel()
703	PauseUpdate; Silent 1 // building window...
704	DoWindow/K DataArithmeticPanel
705	NewPanel /W=(459,44,959,404)/N=DataArithmeticPanel/K=2 as "Data Set Arithmetic"
706	ModifyPanel fixedSize=1
707
708	//Main bit of panel
709	GroupBox grpBox_0,pos={20,10},size={460,105}
710
711	Button DS1_button,pos={300,20},size={150,20},proc=DA_LoadDataSetProc,title="Load 1D Data Set 1"
712	Button DS1_button,valueColor=(65535,0,0),userdata="DS1"
713	PopupMenu DS1_popup,pos={30,21},size={318,20},title="Data Set 1"
714	PopupMenu DS1_popup,mode=1,value= #"DM_DataSetPopupList()"
715	PopupMenu DS1_popup,proc=DA_PopupProc
716	PopupMenu DS1_popup,fsize=12,fcolor=(65535,0,0),valueColor=(65535,0,0)
717
718	Button DS2_button,pos={300,50},size={150,20},proc=DA_LoadDataSetProc,title="Load 1D Data Set 2"
719	Button DS2_button,valueColor=(0,0,65535),userdata="DS2"
720	PopupMenu DS2_popup,pos={30,51},size={318,20},title="Data Set 2"
721	PopupMenu DS2_popup,mode=1,value= #"DM_DataSetPopupList()"
722	PopupMenu DS2_popup,proc=DA_PopupProc
723	PopupMenu DS2_popup,fsize=12,fcolor=(0,0,65535),valueColor=(0,0,65535)
724
725	Button DAPlot_button,title="Plot",pos={100,85},size={150,20}
726	Button DAPlot_button,proc=DAPlotButtonProc
727	Button DADone_button,title="Done",pos={360,85},size={60,20}
728	Button DADone_button,proc=DADoneButtonProc
729	Button DAHelp_button,title="?",pos={440,85},size={30,20}
730	Button DAHelp_button,proc=DAHelpButtonProc
731
732
733	//Tabs
734	TabControl DATabs,pos={20,120},size={460,220},tabLabel(0)="Subtract", proc=DATabsProc
735	TabControl DATabs,tablabel(1)="Add",tablabel(2)="Multiply",tablabel(3)="Divide"
736	TabControl DATabs,value=0
737
738	Button DACalculate_button,title="Calculate",pos={50,310},size={150,20}
739	Button DACalculate_button,proc=DACalculateProc
740	Button DASave_button,title="Save Result",pos={300,310},size={150,20}
741	Button DASave_button,proc=DASaveProc,disable=2
742	Button DACursors_button,title="Get Matching Range",pos={175,250},size={150,20}
743	Button DACursors_button,proc=DACursorButtonProc
744
745	SetVariable DAResultName_sv,title="Result Name (max 25 characters)",pos={50,280},size={400,20}
746	SetVariable DAResultName_Sv,fsize=12,proc=DANameSetvarproc,live=1
747	//Update the result name
748	ControlInfo/W=DataArithmeticPanel DS1_popup
749	if (cmpstr(S_Value,"No data loaded") == 0)
750	SetVariable DAResultName_sv,value=_STR:"SubtractionResult"
751	else
752	//fake call to popup
753	STRUCT WMPopupAction pa
754	pa.win = "DataArithmeticPanel"
755	pa.ctrlName = "DS1_popup"
756	pa.eventCode = 2
757	DA_PopupProc(pa)
758	endif
759
760	CheckBox DANoDS2_cb,title="Data Set 2 = 1?",pos={300,180}
761	CheckBox DANoDS2_cb,proc=DANoDS2Proc
762
763	ValDisplay DARangeStar_vd,title="Start",pos={40,220},size={100,20},fsize=12,value=_NUM:0
764	ValDisplay DARangeEnd_vd,title="End ",pos={160,220},size={100,20},fsize=12,value=_NUM:0
765
766	SetVariable DAScale_sv,title="Scale Factor (f)",pos={280,220},size={180,20},fsize=12,value=_NUM:1
767
768	GroupBox grpBox_1,pos={30,210},size={440,70}
769
770	NewPanel/HOST=DataArithmeticPanel/N=arithDisplay/W=(50,150,170,190)
771
772	//Update the result name
773	ControlInfo/W=DataArithmeticPanel DS1_popup
774
775
776
777	arithDisplayProc(0)
778
779	End
780
781
782	Function MakeDAPlotPanel()
783	PauseUpdate; Silent 1 // building window...
784	DoWindow/K DAPlotPanel
785	NewPanel /W=(14,44,454,484)/N=DAPlotPanel/K=1 as "Data Set Arithmetic"
786	ModifyPanel fixedSize=1
787
788	Display/HOST=DAPlotPanel/N=DAPlot/W=(0,0,440,400)
789	Legend
790	ShowInfo
791	SetActiveSubWindow DAPlotPanel
792	Checkbox DAPlot_log_cb, title="Log I(q)", pos={20,410},value=0
793	Checkbox DAPlot_log_cb, proc=DALogLinIProc
794	Checkbox DAPlot_lin_cb, title="High Q Linear", pos={100,410},value=0
795	Checkbox DAPlot_lin_cb, proc=DAHighQLinProc
796
797	End
798
799	Function AddDAPlot(dataset)
800	Variable dataset
801
802	String win = "DataArithmeticPanel"
803	String DS1name,DS2name,ResultName
804
805	switch(dataset)
806	case 1:
807	ControlInfo/W=$(win) DS1_popup
808	DS1name = S_Value
809	Wave qWave = $("root:"+DS1name+":"+DS1name+"_q")
810	Wave iWave = $("root:"+DS1name+":"+DS1name+"_i")
811	Wave errWave = $("root:"+DS1name+":"+DS1name+"_s")
812	AppendToGraph/W=DAPlotPanel#DAPlot iWave vs Qwave
813	ErrorBars/W=DAPlotPanel#DAPlot /T=0 $(DS1name+"_i"), Y wave=(errWave,errWave)
814	ModifyGraph/W=DAPlotPanel#DAPlot rgb($(DS1name+"_i"))=(65535,0,0)
815	ControlInfo/W=$(win) DANoDS2_cb
816	// if (V_Value == 1)
817	// Cursor/W=DAPlotPanel#DAPlot A, $(DS1name+"_i"), leftx(iWave)
818	// Cursor/W=DAPlotPanel#DAPlot/A=0 B, $(DS1name+"_i"), rightx(iWave)
819	// endif
820	break
821	case 2:
822	ControlInfo/W=$(win) DANoDS2_cb
823	if (V_Value == 0)
824	ControlInfo/W=$(win) DS2_popup
825	DS2name = S_Value
826	if(cmpstr(DS2name,"No data loaded")==0)
827	break //in case someone loads set 1, but not set two, then plots
828	endif
829	Wave qWave = $("root:"+DS2name+":"+DS2name+"_q")
830	Wave iWave = $("root:"+DS2name+":"+DS2name+"_i")
831	Wave errWave = $("root:"+DS2name+":"+DS2name+"_s")
832	AppendToGraph/W=DAPlotPanel#DAPlot iWave vs Qwave
833	ErrorBars/W=DAPlotPanel#DAPlot /T=0 $(DS2name+"_i"), Y wave=(errWave,errWave)
834	ModifyGraph/W=DAPlotPanel#DAPlot rgb($(DS2name+"_i"))=(0,0,65535)
835	Cursor/W=DAPlotPanel#DAPlot A, $(DS2name+"_i"), leftx(iWave)
836	Cursor/W=DAPlotPanel#DAPlot/A=0 B, $(DS2name+"_i"), rightx(iWave)
837	else
838	ControlInfo/W=$(win) DS1_popup
839	DS1name = S_Value
840	DuplicateDataSet("root:"+DS1name,"NullSolvent",1)
841	Wave qWave =root:NullSolvent:NullSolvent_q
842	Wave iWave = root:NullSolvent:NullSolvent_i
843	Wave errWave = root:NullSolvent:NullSolvent_s
844	Wave iWaveDS1 = $("root:"+DS1name+":"+DS1name+"_i")
845	iWave = 1
846	errWave = 0
847	AppendToGraph/W=DAPlotPanel#DAPlot iWave vs Qwave
848	ErrorBars/W=DAPlotPanel#DAPlot /T=0 NullSolvent_i, Y wave=(errWave,errWave)
849	ModifyGraph/W=DAPlotPanel#DAPlot rgb(NullSolvent_i)=(0,0,65535)
850	//Cursor/W=DAPlotPanel#DAPlot A, NullSolvent_i, leftx(iWave)
851	//Cursor/W=DAPlotPanel#DAPlot/A=0 B, NullSolvent_i, rightx(iWave)
852	if(strlen(CsrInfo(A,"DAPlotPanel#DAPlot")) == 0) //cursors not already on the graph
853	Cursor/W=DAPlotPanel#DAPlot A, $(DS1Name+"_i"), leftx(iWaveDS1)
854	Cursor/W=DAPlotPanel#DAPlot/A=0 B, $(DS1Name+"_i"), rightx(iWaveDS1)
855	endif
856	endif
857	break
858	case 3:
859	ControlInfo/W=$(win) DAResultName_sv
860	ResultName = S_Value
861	Wave qWave = $("root:"+ResultName+":"+ResultName+"_q")
862	Wave iWave = $("root:"+ResultName+":"+ResultName+"_i")
863	Wave errWave = $("root:"+ResultName+":"+ResultName+"_s")
864	AppendToGraph/W=DAPlotPanel#DAPlot iWave vs Qwave
865	ErrorBars/W=DAPlotPanel#DAPlot /T=0 $(ResultName+"_i"), Y wave=(errWave,errWave)
866	ModifyGraph/W=DAPlotPanel#DAPlot rgb($(ResultName+"_i"))=(0,65535,0)
867	break
868	endswitch
869
870	ControlInfo/W=DAPlotPanel DAPlot_log_cb
871	ModifyGraph/W=DAPlotPanel#DAPlot mode=3, msize=2, marker=19, mirror=1, tick=2, log(bottom)=1,log(left)=V_Value,tickUnit=1
872	End
873
874	Function arithDisplayProc(s)
875	Variable s
876
877	SetActiveSubWindow DataArithmeticPanel#arithDisplay
878
879	switch (s)
880	case 0:
881	//Subtract
882	DrawAction/L=progFront delete
883	DrawRect 0,0,120,40
884	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,65535,0)
885	DrawText 10,32,"I"
886	SetDrawEnv fname="Times", fsize=22, fstyle= 1,textrgb= (0,0,0)
887	DrawText 20,30,"="
888	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (65535,0,0)
889	DrawText 35,32,"I"
890	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,65535)
891	DrawText 95,32,"I"
892	SetDrawEnv fname="Symbol", fsize=22, fstyle= 1,textrgb= (0,0,0)
893	DrawText 52,32,"-"
894	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,0)
895	DrawText 75,30,"f"
896	SetDrawEnv fname="Times", fsize=22, fstyle= 1,textrgb= (0,0,0)
897	DrawText 85,32,"*"
898	break
899	case 1:
900	//Add
901	DrawAction/L=progFront delete
902	DrawRect 0,0,120,40
903	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,65535,0)
904	DrawText 10,32,"I"
905	SetDrawEnv fname="Times", fsize=22, fstyle= 1,textrgb= (0,0,0)
906	DrawText 20,30,"="
907	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (65535,0,0)
908	DrawText 35,32,"I"
909	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,65535)
910	DrawText 95,32,"I"
911	SetDrawEnv fname="Symbol", fsize=22, fstyle= 1,textrgb= (0,0,0)
912	DrawText 52,32,"+"
913	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,0)
914	DrawText 75,30,"f"
915	SetDrawEnv fname="Times", fsize=22, fstyle= 1,textrgb= (0,0,0)
916	DrawText 85,32,"*"
917	break
918	case 2:
919	//Multiply
920	DrawAction/L=progFront delete
921	DrawRect 0,0,120,40
922	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,65535,0)
923	DrawText 10,32,"I"
924	SetDrawEnv fname="Times", fsize=22, fstyle= 1,textrgb= (0,0,0)
925	DrawText 20,30,"="
926	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (65535,0,0)
927	DrawText 35,32,"I"
928	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,65535)
929	DrawText 95,32,"I"
930	SetDrawEnv fname="Symbol", fsize=22, fstyle= 1,textrgb= (0,0,0)
931	DrawText 52,32,"*"
932	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,0)
933	DrawText 66,30,"("
934	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,0)
935	DrawText 75,30,"f"
936	SetDrawEnv fname="Times", fsize=22, fstyle= 1,textrgb= (0,0,0)
937	DrawText 85,32,"*"
938	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,0)
939	DrawText 105,30,")"
940	break
941	case 3:
942	//Divide
943	DrawAction/L=progFront delete
944	DrawRect 0,0,120,40
945	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,65535,0)
946	DrawText 10,32,"I"
947	SetDrawEnv fname="Times", fsize=22, fstyle= 1,textrgb= (0,0,0)
948	DrawText 20,30,"="
949	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (65535,0,0)
950	DrawText 35,32,"I"
951	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,65535)
952	DrawText 95,32,"I"
953	SetDrawEnv fname="Symbol", fsize=22, fstyle= 1,textrgb= (0,0,0)
954	DrawText 52,32,"/"
955	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,0)
956	DrawText 66,30,"("
957	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,0)
958	DrawText 75,30,"f"
959	SetDrawEnv fname="Times", fsize=22, fstyle= 1,textrgb= (0,0,0)
960	DrawText 85,32,"*"
961	SetDrawEnv fname="Times", fsize=22, fstyle= 2,textrgb= (0,0,0)
962	DrawText 105,30,")"
963	break
964	endswitch
965
966	SetActiveSubWindow DataArithmeticPanel
967
968	return 0
969	End
970
971	//Must follow naming scheme to match buttons to popups
972	//"Name_button" goes with "Name_popup"
973	Function DA_LoadDataSetProc(ba) : ButtonControl
974	STRUCT WMButtonAction &ba
975
976	switch( ba.eventCode )
977	case 2: // mouse up
978	// click code here
979	String cmd = "A_LoadOneDDataWithName(\"\","+num2str(0)+")"
980	Execute cmd
981
982	SVAR gLastFileName = root:packages:NIST:gLastFileName
983
984	String windowName = ba.win
985	String popupName = StringFromList(0,ba.ctrlName,"_")+"_popup"
986
987	ControlUpdate/W=$(windowName) $(popupName)
988	//instead of a simple controlUpdate, better to pop the menu to make sure that the other menus follow
989	// convoluted method to find the right item and pop the menu.
990
991	String list,folderStr
992	Variable num
993	folderStr = CleanupName(gLastFileName,0)
994	list = DM_DataSetPopupList()
995	num=WhichListItem(folderStr,list,";",0,0)
996	if(num != -1)
997	PopupMenu $(popupName),mode=num+1,win=$(windowName)
998	ControlUpdate/W=$(windowName) $(popupName)
999
1000	endif
1001	//fake call to popup
1002	STRUCT WMPopupAction pa
1003	pa.win = ba.win
1004	pa.ctrlName = "DS1_popup"
1005	pa.eventCode = 2
1006	DA_PopupProc(pa)
1007	break
1008	endswitch
1009
1010	return 0
1011	End
1012
1013	Function DA_PopupProc(pa) : PopupMenuControl
1014	STRUCT WMPopupAction &pa
1015
1016	String resultName
1017
1018	switch( pa.eventCode)
1019	case 2:
1020	//print "Called by "+pa.ctrlname+" with value "+pa.popStr
1021	ControlInfo/W=$(pa.win) $(pa.ctrlName)
1022	String popStr = S_Value
1023	if (stringmatch(pa.ctrlname,"DS1") == 1)
1024	resultName = stringfromlist(0,popStr,"_")+"_mod"
1025
1026	SetVariable DAResultName_sv win=$(pa.win), value=_STR:resultName
1027	endif
1028	break
1029	endswitch
1030
1031
1032	End
1033
1034	// function to control the drawing of buttons in the TabControl on the main panel
1035	// Naming scheme for the buttons MUST be strictly adhered to... else buttons will
1036	// appear in odd places...
1037	// all buttons are named MainButton_NA where N is the tab number and A is the letter denoting
1038	// the button's position on that particular tab.
1039	// in this way, buttons will always be drawn correctly..
1040	//
1041	Function DATabsProc(tca) : TabControl
1042	STRUCT WMTabControlAction &tca
1043
1044	switch (tca.eventCode)
1045	case 2:
1046	// Print "name,number",name,tab
1047	String ctrlList = ControlNameList("",";"),item="",nameStr=""
1048	Variable num = ItemsinList(ctrlList,";"),ii,onTab
1049	for(ii=0;ii<num;ii+=1)
1050	//items all start w/"DSTabItem_"
1051	item=StringFromList(ii, ctrlList ,";")
1052	nameStr=item[0,9]
1053	if(cmpstr(nameStr,"DATabItem_")==0)
1054	onTab = str2num(item[10])
1055	ControlInfo $item
1056	switch (V_flag)
1057	case 1:
1058	Button $item,disable=(tca.tab!=onTab)
1059	break
1060	case 2:
1061	CheckBox $item,disable=(tca.tab!=onTab)
1062	break
1063	case 3:
1064	PopUpMenu $item,disable=(tca.tab!=onTab)
1065	break
1066	case 5:
1067	SetVariable $item,disable=(tca.tab!=onTab)
1068	break
1069	endswitch
1070	endif
1071	endfor
1072
1073	arithDisplayProc(tca.tab)
1074	break
1075	endswitch
1076	End
1077
1078
1079	Function DACalculateProc(ba) : ButtonControl
1080	STRUCT WMButtonAction &ba
1081
1082	String DS1,DS2,Resultname
1083
1084	switch(ba.eventCode)
1085	case 2:
1086	//Which tab?
1087	ControlInfo/W=$(ba.win) DATabs
1088	Variable tabNum = V_value
1089	//Print "Tab number "+num2str(tabNum)
1090
1091	ControlInfo/W=$(ba.win) DS1_popup
1092	DS1 = S_Value
1093	ControlInfo/W=$(ba.win) DANoDS2_cb
1094	if (V_Value == 0)
1095	ControlInfo/W=$(ba.win) DS2_popup
1096	DS2 = S_Value
1097	else
1098	DS2 = "NullSolvent"
1099	endif
1100	ControlInfo/W=$(ba.win) DAResultName_sv
1101	Resultname = CleanupName(S_Value, 0 ) //clean up any bad characters, and put the cleaned string back
1102	SetVariable DAResultName_sv,value=_STR:ResultName
1103
1104	ControlInfo/W=$(ba.win) DAScale_sv
1105	Variable Scalefactor = V_Value
1106
1107	switch(tabNum)
1108	case 0:
1109	//do subtraction
1110	//print "Subtraction of "+DS2+" from "+DS1+" with sf "+num2str(scalefactor)+ " into "+Resultname
1111	SubtractDataSets(DS1,DS2,Scalefactor,Resultname)
1112	break
1113	case 1:
1114	//do addition
1115	AddDataSets(DS1,DS2,Scalefactor,Resultname)
1116	break
1117	case 2:
1118	//do multiplication
1119	MultiplyDataSets(DS1,DS2,Scalefactor,Resultname)
1120	break
1121	case 3:
1122	//do division
1123	DivideDataSets(DS1,DS2,Scalefactor,Resultname)
1124	break
1125	endswitch
1126
1127	//Sort out plot
1128	//Fake button press to DAPlotButtonProc
1129	STRUCT WMButtonAction ba2
1130	ba2.win = ba.win
1131	ba2.ctrlName = "DAPlot_button"
1132	ba2.eventCode = 2
1133
1134	// I've commented this out - the cursors get reset to the ends since this removes all sets from the graph, and
1135	// then replots them. What is the real purpose of this call? To clear the old result off before adding the
1136	// new one?
1137	// DAPlotButtonProc(ba2)
1138	ba2.userData = ResultName
1139	DAPlotRemoveResult(ba2)
1140
1141
1142	AddDAPlot(3)
1143	DoWindow/F DataArithmeticPanel
1144
1145	//Enable save button now that we have a result to save
1146	Button DASave_Button win=$(ba.win),disable=0
1147
1148	// SetActiveSubWindow DAPlotPanel
1149	endswitch
1150
1151	End
1152
1153	// remove what is not the
1154	//
1155	Function DAPlotRemoveResult(ba) : ButtonControl
1156	STRUCT WMButtonAction &ba
1157
1158	String win = ba.win
1159	String ResultName = ba.userData
1160	String item="",traceList=""
1161	Variable ii=0,num
1162
1163	switch (ba.eventCode)
1164	case 2: //mouse up
1165	//data set 1
1166	ControlInfo/W=$(win) DS1_popup
1167	String DS1 = S_Value
1168
1169	//Get folder for DS2
1170	ControlInfo/W=$(win) DS2_popup
1171	String DS2 = S_Value
1172
1173	// state of the checkbox
1174	ControlInfo/W=$(win) DANoDS2_cb
1175	if(V_Value)
1176	DS2 = "NullSolvent"
1177	endif
1178
1179	DoWindow DAPlotPanel
1180	if (V_Flag == 0)
1181	MakeDAPlotPanel()
1182	else
1183	DoWindow/HIDE=0/F DAPlotPanel
1184	traceList = TraceNameList("DAPlotPanel#DAPlot",";",1)
1185	num=ItemsInList(traceList)
1186	ii=0
1187	do
1188	item = StringFromList(ii,traceList,";")
1189	if (stringmatch(item,ResultName+"*")==1) //it it's the specific trace I've asked to remove
1190	RemoveFromGraph/W=DAPlotPanel#DAPlot $item
1191	elseif (stringmatch(item,DS1+"")==0 && stringmatch(item,DS2+"")==0) //if it's not set1 & not set2
1192	RemoveFromGraph/W=DAPlotPanel#DAPlot $item
1193	endif
1194
1195	ii+=1
1196	while(ii<num)
1197	endif
1198
1199	break
1200	endswitch
1201
1202	return 0
1203	End
1204
1205
1206	Function DAPlotButtonProc(ba) : ButtonControl
1207	STRUCT WMButtonAction &ba
1208
1209	String win = ba.win
1210
1211	switch (ba.eventCode)
1212	case 2:
1213	//mouse up
1214	//Get folder for DS1
1215	DoWindow DAPlotPanel
1216	if (V_Flag == 0)
1217	MakeDAPlotPanel()
1218	else
1219	DoWindow/HIDE=0/F DAPlotPanel
1220	do
1221	String tracename = StringFromList(0,TraceNameList("DAPlotPanel#DAPlot",";",1),";")
1222	if (cmpstr(tracename,"")==0)
1223	break
1224	else
1225	RemoveFromGraph/W=DAPlotPanel#DAPlot $tracename
1226	endif
1227	while(1)
1228	endif
1229
1230	ControlInfo/W=$(win) DS1_popup
1231	String DS1 = S_Value
1232	if (cmpstr(DS1,"") != 0 )
1233	AddDAPlot(1)
1234	endif
1235	//Get folder for DS2
1236	ControlInfo/W=$(win) DS2_popup
1237	String DS2 = S_Value
1238	if (cmpstr(DS2,"") != 0)
1239	AddDAPlot(2)
1240	endif
1241	break
1242	endswitch
1243
1244	return 0
1245	End
1246
1247	Function DADoneButtonProc(ba) : ButtonControl
1248	STRUCT WMButtonAction &ba
1249
1250	String win = ba.win
1251
1252	switch (ba.eventCode)
1253	case 2:
1254	DoWindow/K DAPlotPanel
1255	DoWindow/K DataArithmeticPanel
1256	break
1257	endswitch
1258
1259	return 0
1260	End
1261
1262	Function DAHelpButtonProc(ba) : ButtonControl
1263	STRUCT WMButtonAction &ba
1264
1265	String win = ba.win
1266
1267	switch (ba.eventCode)
1268	case 2:
1269	// click code here
1270	DisplayHelpTopic/Z/K=1 "Data Set Arithmetic"
1271	if(V_flag !=0)
1272	DoAlert 0,"The Data Set Arithmetic Help file could not be found"
1273	endif
1274	break
1275	endswitch
1276
1277	return 0
1278	End
1279
1280	Function DALogLinIProc(cba) : CheckBoxControl
1281	STRUCT WMCheckBoxAction &cba
1282
1283	switch(cba.eventcode)
1284	case 2:
1285
1286	ModifyGraph/W=DAPlotPanel#DAPlot log(left)=cba.checked
1287	ModifyGraph/W=DAPlotPanel#DAPlot log(bottom)=cba.checked
1288	ModifyGraph/W=DAPlotPanel#DAPlot zero(left)=0
1289	SetAxis/A/W=DAPlotPanel#DAPlot
1290
1291	if(cba.checked)
1292	Checkbox DAPlot_lin_cb,value=0 //uncheck lin
1293	endif
1294	endswitch
1295
1296
1297	End
1298
1299	Function DAHighQLinProc(cba) : CheckBoxControl
1300	STRUCT WMCheckBoxAction &cba
1301
1302	switch(cba.eventcode)
1303	case 2:
1304	if(cba.checked)
1305	ModifyGraph/W=DAPlotPanel#DAPlot log=0,zero(left)=1
1306	SetAxis/W=DAPlotPanel#DAPlot left -0.1,0.1
1307	SetAxis/W=DAPlotPanel#DAPlot bottom 0.1,*
1308	SetAxis/W=DAPlotPanel#DAPlot left -0.02,0.02
1309
1310	Checkbox DAPlot_log_cb,value=0 //uncheck the log
1311	endif
1312	endswitch
1313
1314
1315	End
1316
1317
1318	Function DACursorButtonProc(ba) : ButtonControl
1319	STRUCT WMButtonAction &ba
1320
1321	String DS1,DS2
1322
1323	switch(ba.eventCode)
1324	case 2:
1325
1326	ControlInfo/W=$(ba.win) DS1_popup
1327	DS1 = S_Value
1328	ControlInfo/W=$(ba.win) DANoDS2_cb
1329	Variable NoDS2 = V_Value
1330	if (NoDS2 == 0)
1331	ControlInfo/W=$(ba.win) DS2_popup
1332	DS2 = S_Value
1333	else
1334	DS2 = "NullSolvent"
1335	endif
1336
1337	//AJJ Nov 2009 - UGLY - Will have to revisit this when we deal with hierarchical folders
1338	Wave set1_i = $("root:"+DS1+":"+DS1+"_i")
1339	Wave set1_q = $("root:"+DS1+":"+DS1+"_q")
1340	Wave set2_i = $("root:"+DS2+":"+DS2+"_i")
1341	Wave set2_q = $("root:"+DS2+":"+DS2+"_q")
1342	Duplicate/O set1_i tmp_i
1343	Duplicate/O set1_q tmp_q
1344	tmp_i = set1_i / interp(set1_q[p],set2_q,set2_i)
1345
1346	//Get cursors
1347	Variable q1,q2
1348
1349	DoWindow/F DAPlotPanel
1350	SetActiveSubWindow DAPlotPanel#DAPlot
1351
1352	q1 = CsrXWaveRef(A)[pcsr(A)]
1353	q2 = CsrXWaveRef(B)[pcsr(B)]
1354
1355	//Update value display
1356	ValDisplay DARangeStar_vd,value=_NUM:q1, win=$(ba.win)
1357	ValDisplay DARangeEnd_vd,value=_NUM:q2, win=$(ba.win)
1358
1359	//Calculate scalefactor
1360
1361	if (NoDS2 == 1)
1362	Wave avgWave = set1_i
1363	else
1364	Wave avgWave = tmp_i
1365	endif
1366
1367	Variable p1 = BinarySearch(tmp_q,q1)
1368	Variable p2 = BinarySearch(tmp_q,q2)
1369
1370	//print avgWave
1371
1372	WaveStats/Q/R=[p1,p2] avgWave
1373	//print V_avg
1374	//Update sv control
1375	SetVariable DAScale_sv, value=_NUM:V_avg, win=$(ba.win)
1376
1377	KillWaves/Z tmp_i,tmp_q
1378	DoWindow/F DataArithmeticPanel
1379	endswitch
1380
1381	End
1382
1383
1384	Function DANoDS2Proc(cba) : CheckBoxControl
1385	STRUCT WMCheckBoxAction &cba
1386
1387
1388	switch(cba.eventCode)
1389	case 2:
1390	if (cba.checked == 1)
1391	//Disable DS2 popup etc
1392	PopupMenu DS2_popup win=$(cba.win), disable=2
1393	Button DS2_button win=$(cba.win), disable=2
1394	else
1395	//Enable DS2 popup etc
1396	PopupMenu DS2_popup win=$(cba.win), disable=0
1397	Button DS2_button win=$(cba.win), disable=0
1398	endif
1399	//Sort out plot
1400	//Fake button press to DAPlotButtonProc
1401	STRUCT WMButtonAction ba2
1402	ba2.win = cba.win
1403	ba2.ctrlName = "DAPlot_button"
1404	ba2.eventCode = 2
1405	DAPlotButtonProc(ba2)
1406	SetActiveSubWindow DAPlotPanel
1407	DoWindow/F DataArithmeticPanel
1408	endswitch
1409
1410	End
1411
1412	Function DASaveProc(ba) : ButtonControl
1413	STRUCT WMButtonAction &ba
1414
1415	switch(ba.eventCode)
1416	case 2:
1417	ControlInfo/W=$(ba.win) DAResultName_sv
1418	SaveDataSetToFile(S_Value)
1419	break
1420	endswitch
1421
1422
1423	End
1424
1425	/////////////////////// Common Panel Functions ///////////////////////////////////////
1426
1427
1428	// is there a simpler way to do this? I don't think so.
1429	Function/S DM_DataSetPopupList()
1430
1431	String str=GetAList(4)
1432
1433	if(strlen(str)==0)
1434	str = "No data loaded"
1435	endif
1436	str = SortList(str)
1437
1438	return(str)
1439	End
1440
1441
1442	Function DANameSetvarproc(sva) : SetVariableControl
1443	STRUCT WMSetVariableAction &sva
1444
1445	switch( sva.eventCode )
1446	case 1: // mouse up
1447	case 2: // Enter key
1448	case 3: // Live update
1449	String sv = sva.sval
1450	if( strlen(sv) > 25 )
1451	sv= sv[0,24]
1452	SetVariable $(sva.ctrlName),win=$(sva.win),value=_STR:sv
1453	ControlUpdate /W=$(sva.win) $(sva.ctrlName)
1454	Beep
1455	endif
1456	Button DASave_Button win=$(sva.win), disable=2
1457	break
1458	endswitch
1459	return 0
1460	End
1461
1462
1463	////////////////////// Functions to do manipulations ///////////////////////////////////
1464
1465	Function RenameDataSet(dataSetFolder, newName)
1466	String dataSetFolder
1467	String newName
1468
1469	String dataSetFolderParent,basestr,objName
1470	Variable index = 0
1471
1472	//Abuse ParseFilePath to get path without folder name
1473	dataSetFolderParent = ParseFilePath(1,dataSetFolder,":",1,0)
1474	//Abuse ParseFilePath to get basestr
1475	basestr = ParseFilePath(0,dataSetFolder,":",1,0)
1476
1477	// try
1478	RenameDataFolder $(dataSetFolder) $(newName)//; AbortOnRTE
1479
1480
1481	SetDataFolder $(dataSetFolderParent+newName)//; AbortOnRTE
1482	do
1483	objName = GetIndexedObjName("",1,index)
1484	if (strlen(objName) == 0)
1485	break
1486	endif
1487	Rename $(objName) $(ReplaceString(basestr,objName,newName))
1488	index+=1
1489	while(1)
1490	SetDataFolder root:
1491	// catch
1492	// Print "Aborted: " + num2str(V_AbortCode)
1493	// SetDataFolder root:
1494	// endtry
1495	End
1496
1497
1498	Function DuplicateDataSet(dataSetFolder, newName, forceoverwrite)
1499	String dataSetFolder
1500	String newName
1501	Variable forceoverwrite
1502
1503	String dataSetFolderParent,basestr,objName
1504	Variable index = 0
1505
1506	//Abuse ParseFilePath to get path without folder name
1507	dataSetFolderParent = ParseFilePath(1,dataSetFolder,":",1,0)
1508	//Abuse ParseFilePath to get basestr
1509	basestr = ParseFilePath(0,dataSetFolder,":",1,0)
1510
1511	print "Duplicating "+dataSetFolder+" as "+newName
1512
1513	SetDataFolder $(dataSetFolderParent)
1514
1515	if (!DataFolderExists(newName))
1516	NewDataFolder $(newName)
1517	else
1518	if (!forceoverwrite)
1519	DoAlert 1, "A dataset with the name "+newName+" already exists. Overwrite?"
1520	if (V_flag == 2)
1521	return 1
1522	endif
1523	endif
1524	endif
1525
1526	//If we are here, the folder (now) exists and the user has agreed to overwrite
1527	//either in the function call or from the alert.
1528
1529	// here, GetIndexedObjectName copies all of the waves
1530	index = 0
1531	do
1532	objName = GetIndexedObjName(basestr,1,index)
1533	if (strlen(objName) == 0)
1534	break
1535	endif
1536	objname = ":"+basestr+":"+objname
1537	Duplicate/O $(objName) $(ReplaceString(basestr,objName,newName))
1538	index+=1
1539	while(1)
1540
1541	// -- for USANS data, we need the slit height. copy all of the "USANS_*" variables
1542	// may need to augment this for other situations
1543	index = 0
1544	do
1545	objName = GetIndexedObjName(basestr,2,index)
1546	if (strlen(objName) == 0)
1547	break
1548	endif
1549	if(stringmatch(objName,"USANS*") == 1)
1550	objname = ":"+basestr+":"+objname
1551	NVAR tmp = $objName
1552	Variable/G $(ReplaceString(basestr,objName,newName))= tmp
1553	endif
1554	index+=1
1555	while(1)
1556
1557
1558	SetDataFolder root:
1559	return 0
1560	End
1561
1562
1563	// Subtract Set2 From Set1
1564	// Use Result_I = Set1_I - f*Set2_I
1565	Function SubtractDataSets(set1Name,set2Name,set2Scale,resultName)
1566	String set1Name
1567	String set2Name
1568	Variable set2Scale
1569	String resultName
1570
1571	String set1Path = "root:"+set1Name+":"
1572	String set2Path = "root:"+set2Name+":"
1573	String resultPath = "root:"+resultName+":"
1574
1575	SetDataFolder root:
1576	//Create folder for result
1577	//UnloadDataSet(resultName)
1578	//Make the folder
1579	if (DuplicateDataSet(set1Path,resultName,0))
1580	return 1
1581	else
1582	//Do subtraction of I waves - including interpolation if necessary.
1583	Wave result_i = $(resultPath+resultName+"_i")
1584	Wave result_s = $(resultPath+resultName+"_s")
1585	Wave set1_i = $(set1Path+set1Name+"_i")
1586	Wave set1_q = $(set1Path+set1Name+"_q")
1587	Wave set1_s = $(set1Path+set1Name+"_s")
1588	Wave set2_i = $(set2Path+set2Name+"_i")
1589	Wave set2_q = $(set2Path+set2Name+"_q")
1590	Wave set2_s = $(set2Path+set2Name+"_s")
1591
1592	result_i = set1_i - (set2Scale*interp(set1_q[p],set2_q,set2_i))
1593	result_s = sqrt(set1_s^2 + (set2Scale*interp(set1_q[p],set2_q,set2_s))^2 )
1594	//Calculate result error wave - can we produce corrected Q error?
1595
1596	//Generate history string to record what was done?
1597	return 0
1598	endif
1599	End
1600
1601	// Add Set2 to Set1
1602	// Use Result_I = Set1_I + f*Set2_I
1603	Function AddDataSets(set1Name,set2Name,set2Scale,resultName)
1604	String set1Name
1605	String set2Name
1606	Variable set2Scale
1607	String resultName
1608
1609	String set1Path = "root:"+set1Name+":"
1610	String set2Path = "root:"+set2Name+":"
1611	String resultPath = "root:"+resultName+":"
1612
1613	SetDataFolder root:
1614	//Create folder for result
1615	if(DuplicateDataSet(set1Path,resultName,0))
1616	//User said no overwrite
1617	return 1
1618	else
1619	//Do addition of I waves - including interpolation if necessary.
1620	Wave result_i = $(resultPath+resultName+"_i")
1621	Wave result_s = $(resultPath+resultName+"_s")
1622	Wave set1_i = $(set1Path+set1Name+"_i")
1623	Wave set1_q = $(set1Path+set1Name+"_q")
1624	Wave set1_s = $(set1Path+set1Name+"_s")
1625	Wave set2_i = $(set2Path+set2Name+"_i")
1626	Wave set2_q = $(set2Path+set2Name+"_q")
1627	Wave set2_s = $(set2Path+set2Name+"_s")
1628
1629	result_i = set1_i + set2Scale*interp(set1_q[p],set2_q,set2_i)
1630	//Calculate result error wave (note that this is identical to subtraction)
1631	result_s = sqrt(set1_s^2 + (set2Scale*interp(set1_q[p],set2_q,set2_s))^2 )
1632
1633	// - can we produce corrected Q error?
1634	//Generate history string to record what was done?
1635	return 0
1636	endif
1637	End
1638
1639	// Multiply Set1 by Set2
1640	// Use Result_I = Set1_I * (f*Set2_I)
1641	Function MultiplyDataSets(set1Name, set2Name, set2Scale, resultName)
1642	String set1Name
1643	String set2Name
1644	Variable set2Scale
1645	String resultName
1646
1647	String set1Path = "root:"+set1Name+":"
1648	String set2Path = "root:"+set2Name+":"
1649	String resultPath = "root:"+resultName+":"
1650
1651	SetDataFolder root:
1652	//Create folder for result
1653	//Make the folder
1654	if(DuplicateDataSet(set1Path,resultName,0))
1655	//User said no overwrite
1656	return 1
1657	else
1658	//Do multiplcation of I waves - including interpolation if necessary.
1659	Wave result_i = $(resultPath+resultName+"_i")
1660	Wave result_s = $(resultPath+resultName+"_s")
1661	Wave set1_i = $(set1Path+set1Name+"_i")
1662	Wave set1_q = $(set1Path+set1Name+"_q")
1663	Wave set1_s = $(set1Path+set1Name+"_s")
1664	Wave set2_i = $(set2Path+set2Name+"_i")
1665	Wave set2_q = $(set2Path+set2Name+"_q")
1666	Wave set2_s = $(set2Path+set2Name+"_s")
1667
1668	result_i = set1_iset2Scaleinterp(set1_q[p],set2_q,set2_i)
1669	//Calculate result error wave
1670	// sum each of the relative errors, interpolating set 2 intensity and error as needed
1671	// then sqrt
1672	result_s = (set2Scaleinterp(set1_q[p],set2_q,set2_i)set1_s)^2
1673	result_s += (set2Scaleset1_iinterp(set1_q[p],set2_q,set2_s))^2
1674	result_s = sqrt(result_s)
1675
1676	// - can we produce corrected Q error?
1677
1678	//Generate history string to record what was done?
1679	return 0
1680	endif
1681	End
1682
1683	// Divide Set1 by Set2
1684	// Use Result_I = Set1_I / (f*Set2_I)
1685	Function DivideDataSets(set1Name, set2Name, set2Scale, resultName)
1686	String set1Name
1687	String set2Name
1688	Variable set2Scale
1689	String resultName
1690
1691	String set1Path = "root:"+set1Name+":"
1692	String set2Path = "root:"+set2Name+":"
1693	String resultPath = "root:"+resultName+":"
1694
1695	SetDataFolder root:
1696	//Create folder for result
1697	//Make the folder
1698	if(DuplicateDataSet(set1Path,resultName,0))
1699	//User said no overwrite
1700	return 1
1701	else
1702	//Do division of I waves - including interpolation if necessary.
1703	Wave result_i = $(resultPath+resultName+"_i")
1704	Wave result_s = $(resultPath+resultName+"_s")
1705	Wave set1_i = $(set1Path+set1Name+"_i")
1706	Wave set1_q = $(set1Path+set1Name+"_q")
1707	Wave set1_s = $(set1Path+set1Name+"_s")
1708	Wave set2_i = $(set2Path+set2Name+"_i")
1709	Wave set2_q = $(set2Path+set2Name+"_q")
1710	Wave set2_s = $(set2Path+set2Name+"_s")
1711
1712	result_i = set1_i/(set2Scale*interp(set1_q[p],set2_q,set2_i) )
1713	//Calculate result error wave
1714	// sum each of the relative errors, interpolating set 2 intensity and error as needed
1715	// then sqrt
1716	result_s = (set1_s/set2Scale/interp(set1_q[p],set2_q,set2_i))^2
1717	result_s += (interp(set1_q[p],set2_q,set2_s)*set1_i/set2Scale/interp(set1_q[p],set2_q,set2_i)^2)^2
1718	result_s = sqrt(result_s)
1719
1720	// - can we produce corrected Q error?
1721
1722	//Generate history string to record what was done?
1723	return 0
1724	endif
1725	End
1726
1727	Function ReSortDataSet(set1name)
1728	String set1name
1729
1730	String set1Path = "root:"+set1Name+":"
1731	String curPath = GetDataFolder(1)
1732
1733	SetDataFolder set1Path
1734
1735
1736	//Check for resolution wave
1737	if (exists(set1name+"_res"))
1738	Wave reswave = $(set1name+"_res")
1739
1740	//Check for USANS data - we won't resort these for the moment
1741	if (dimsize(reswave, 1) > 4 )
1742	//USANS data, bail out
1743	print "Can't re-sort USANS data yet!"
1744	return 1
1745	endif
1746
1747	//Break out resolution wave into separate waves
1748	Make/O/N=(numpnts($(set1name+"_q"))) res0 = reswave[p][0]
1749	Make/O/N=(numpnts($(set1name+"_q"))) res1 = reswave[p][1]
1750	Make/O/N=(numpnts($(set1name+"_q"))) res2 = reswave[p][2]
1751	Make/O/N=(numpnts($(set1name+"_q"))) res3 = reswave[p][3]
1752
1753	//sort
1754	print "Re-Sorting 4 or 6 Column Data Set with resolution information: "+set1Name
1755	sort $(set1name+"_q"),$(set1name+"_q"),$(set1name+"_i"),$(set1name+"_s"), res0, res1, res2, res3
1756
1757	//Put resolution contents back
1758	reswave[][0] = res0[p]
1759	reswave[][1]= res1[p]
1760	reswave[][2] = res2[p]
1761	reswave[][3] = res3[p]
1762
1763	//cleanup
1764	Killwaves/Z res0, res1, res2, res3
1765	else
1766	//3 Column only
1767	//sort
1768	print "Re-Sorting 3 Column Data Set: "+set1Name
1769	sort $(set1name+"_q"),$(set1name+"_q"),$(set1name+"_i"),$(set1name+"_s")
1770	endif
1771
1772	SetDataFolder curPath
1773
1774	return 0
1775	End
1776
1777
1778	///////////////////////////Other Utility functions ////////////////////////////
1779
1780	Function SaveDataSetToFile(folderName)
1781	String folderName
1782
1783	String protoStr = ""
1784	//Check for history string in folder
1785	//If it doesn't exist then
1786
1787	//Do saving of data file.
1788
1789	NVAR gXML_Write = root:Packages:NIST:gXML_Write
1790
1791	if (gXML_Write == 1)
1792	ReWrite1DXMLData(folderName)
1793	else
1794	fReWrite1DData(folderName,"tab","CRLF")
1795	endif
1796
1797	//Include history string to record what was done?
1798
1799	End
1800
1801
1802
1803	// still need to get the header information, and possibly the SASprocessnote from the XML load into the 6-column header
1804	//
1805	// start by looking in:
1806	// String xmlReaderFolder = "root:Packages:CS_XMLreader:"
1807	// for Title and Title_folder strings -> then the metadata (but the processnote is still not there
1808	//
1809	// may need to get it directly using the filename
1810	Function convertNISTXMLtoNIST6Col(fname)
1811	String fname
1812
1813	String list, item,path
1814	Variable num,ii
1815
1816	//load the XML
1817
1818	LoadNISTXMLData(fname,"",0,0) //no plot, no force overwrite
1819	// Execute "A_LoadOneDDataWithName(\""+fname+"\",0)" //won't plot
1820
1821	// then rewrite what is in the data folder that was just loaded
1822	String basestr = ParseFilePath(0, fname, ":", 1, 0)
1823	baseStr = CleanupName(baseStr,0)
1824	print fname
1825	print basestr
1826
1827	fReWrite1DData_noPrompt(baseStr,"tab","CR")
1828
1829	return(0)
1830	End
1831
1832
1833	///////// SRK - VERY SIMPLE batch converter
1834	// no header information is preserved
1835	// file names are partially preserved
1836	//
1837
1838	/// to use this:
1839	// -open the Plot Manager and set the path
1840	// -run this function
1841	//
1842	// it doesn't matter if the XML ouput flag is set - this overrides.
1843	Function batchXML26ColConvert()
1844
1845	String list, item,path,fname
1846	Variable num,ii
1847
1848	PathInfo CatPathName
1849	path = S_Path
1850
1851	list = A_ReducedDataFileList("")
1852	num = itemsInList(list)
1853	Print num
1854	for(ii=0;ii<num;ii+=1)
1855	item = StringFromList(ii, list ,";")
1856	fname=path + item
1857	Execute "A_LoadOneDDataWithName(\""+fname+"\",0)" //won't plot
1858	// easier to load all, then write out, since the name will be changed
1859	endfor
1860
1861
1862	list = DM_DataSetPopupList()
1863
1864	num = itemsInList(list)
1865	Print num
1866	for(ii=0;ii<num;ii+=1)
1867	item = StringFromList(ii, list ,";")
1868	fReWrite1DData_noPrompt(item,"tab","CR")
1869	endfor
1870
1871	End
1872
1873	///////// SRK - VERY SIMPLE batch converter
1874	// NO header information is preserved
1875	// file names are partially preserved
1876	//
1877
1878	/// to use this:
1879	// -open the Plot Manager and set the path
1880	// -run this function
1881	//
1882	// it doesn't matter if the XML ouput flag is set - this overrides.
1883	//
1884	// The GRASP output data is 5-column Q-I-errI-sigQ-nCells
1885	// which gets read in as q-i-s-ism-fit_ism (as if it was some wierd USANS data format)
1886	// -- so fake the output...
1887	//
1888	Function batchGrasp26ColConvert()
1889
1890	String list, item,path,fname
1891	Variable num,ii,npt
1892
1893	PathInfo CatPathName
1894	path = S_Path
1895
1896	list = A_ReducedDataFileList("")
1897	num = itemsInList(list)
1898	Print num
1899	for(ii=0;ii<num;ii+=1)
1900	item = StringFromList(ii, list ,";")
1901	fname=path + item
1902	Execute "A_LoadOneDDataWithName(\""+fname+"\",0)" //won't plot
1903	// easier to load all, then write out, since the name will be changed
1904	endfor
1905
1906
1907	list = DM_DataSetPopupList()
1908
1909	num = itemsInList(list)
1910
1911	Print num
1912	for(ii=0;ii<num;ii+=1)
1913	item = StringFromList(ii, list ,";")
1914
1915	// fake the 6-column NIST data structure
1916	WAVE qw = $("root:"+item+":"+item+"_q")
1917	npt = numpnts(qw)
1918	Make/O/D/N=(npt,4) $("root:"+item+":"+item+"_res")
1919	WAVE res = $("root:"+item+":"+item+"_res")
1920	WAVE sigQ = $("root:"+item+":"+item+"_ism")
1921	res[][0] = sigQ[p] // sigQ
1922	res[][1] = qw[p] // qBar ~ q
1923	res[][2] = 1 //shadow
1924	res[][3] = qw[p] // q
1925
1926
1927	fReWrite1DData_noPrompt(item,"tab","CR")
1928	endfor
1929
1930	End
1931
1932	// quick version (copied from fReWrite1DData() that NEVER asks for a new fileName
1933	// - and right now, always expect 6-column data, either SANS or USANS (re-writes -dQv)
1934	// - AJJ Nov 2009 : better make sure we always fake 6 columns on reading then....
1935	Function fReWrite1DData_noPrompt(folderStr,delim,term)
1936	String folderStr,delim,term
1937
1938	String formatStr="",fullpath=""
1939	Variable refnum,dialog=1
1940
1941	String dataSetFolderParent,basestr
1942
1943	//setup delimeter and terminator choices
1944	If(cmpstr(delim,"tab")==0)
1945	//tab-delimeted
1946	formatStr="%15.8g\t%15.8g\t%15.8g\t%15.8g\t%15.8g\t%15.8g"
1947	else
1948	//use 3 spaces
1949	formatStr="%15.8g %15.8g %15.8g %15.8g %15.8g %15.8g"
1950	Endif
1951	If(cmpstr(term,"CR")==0)
1952	formatStr += "\r"
1953	Endif
1954	If(cmpstr(term,"LF")==0)
1955	formatStr += "\n"
1956	Endif
1957	If(cmpstr(term,"CRLF")==0)
1958	formatStr += "\r\n"
1959	Endif
1960
1961	//Abuse ParseFilePath to get path without folder name
1962	dataSetFolderParent = ParseFilePath(1,folderStr,":",1,0)
1963	//Abuse ParseFilePath to get basestr
1964	basestr = ParseFilePath(0,folderStr,":",1,0)
1965
1966	//make sure the waves exist
1967	SetDataFolder $(dataSetFolderParent+basestr)
1968	WAVE/Z qw = $(baseStr+"_q")
1969	WAVE/Z iw = $(baseStr+"_i")
1970	WAVE/Z sw = $(baseStr+"_s")
1971	WAVE/Z resw = $(baseStr+"_res")
1972
1973	if(WaveExists(qw) == 0)
1974	Abort "q is missing"
1975	endif
1976	if(WaveExists(iw) == 0)
1977	Abort "i is missing"
1978	endif
1979	if(WaveExists(sw) == 0)
1980	Abort "s is missing"
1981	endif
1982	if(WaveExists(resw) == 0)
1983	Abort "Resolution information is missing."
1984	endif
1985
1986	Duplicate/O qw qbar,sigQ,fs
1987	if(dimsize(resW,1) > 4)
1988	//it's USANS put -dQv back in the last 3 columns
1989	NVAR/Z dQv = USANS_dQv
1990	if(NVAR_Exists(dQv) == 0)
1991	Abort "It's USANS data, and I don't know what the slit height is."
1992	endif
1993	sigQ = -dQv
1994	qbar = -dQv
1995	fs = -dQv
1996	else
1997	//it's SANS
1998	sigQ = resw[p][0]
1999	qbar = resw[p][1]
2000	fs = resw[p][2]
2001	endif
2002
2003	dialog=0
2004	if(dialog)
2005	PathInfo/S catPathName
2006	// fullPath = DoSaveFileDialog("Save data as",fname=baseStr+".txt")
2007	fullPath = DoSaveFileDialog("Save data as",fname=baseStr[0,strlen(BaseStr)-1])
2008	Print fullPath
2009	If(cmpstr(fullPath,"")==0)
2010	//user cancel, don't write out a file
2011	Close/A
2012	Abort "no data file was written"
2013	Endif
2014	//Print "dialog fullpath = ",fullpath
2015	Endif
2016	PathInfo catPathName
2017	fullPath = S_Path + baseStr[0,strlen(BaseStr)-1]
2018
2019	Open refnum as fullpath
2020
2021	fprintf refnum,"Modified data written from folder %s on %s\r\n",baseStr,(date()+" "+time())
2022	wfprintf refnum,formatStr,qw,iw,sw,sigQ,qbar,fs
2023	Close refnum
2024
2025	KillWaves/Z sigQ,qbar,fs
2026
2027	SetDataFolder root:
2028	return(0)
2029	End
2030
2031	///////////end SRK
2032
2033
2034
2035	///// SRK
2036	///// Rebinning
2037
2038
2039	// main entry procedure for subtraction panel
2040	// re-initializes necessary folders and waves
2041	Proc OpenRebin()
2042	DoWindow/F Rebin_Panel
2043	if(V_Flag==0)
2044	InitRebin()
2045	Rebin() //the panel
2046	// Plot_Sub1D() //the graph
2047	endif
2048
2049	End
2050
2051	Function InitRebin()
2052
2053	NewDataFolder/O/S root:Packages:NIST:Rebin
2054	// Variable/G gPtsBeg1=0
2055	// Variable/G gPtsEnd1=0
2056	Variable/G binning=2
2057
2058	SetDataFolder root:
2059	End
2060
2061	Proc Rebin()
2062
2063	PauseUpdate; Silent 1 // building window...
2064	NewPanel /W=(658,347,920,562)/K=1
2065	DoWindow/C Rebin_Panel
2066	ModifyPanel cbRGB=(65535,48662,45086),fixedSize=1
2067	Button button0,pos={161,178},size={50,20},proc=Rebin_Done,title="Done"
2068	PopupMenu popup0,pos={11,50},size={211,20},title="Data in Memory"
2069	PopupMenu popup0,mode=1,value= #"A_OneDDataInMemory()"
2070	Button button1,pos={60,14},size={150,20},proc=Rebin_Load,title="Load Data From File"
2071	Button button2,pos={118,84},size={100,20},proc=Rebin_Append,title="Append Data"
2072	Button button3,pos={11,84},size={80,20},proc=Rebin_newGraph,title="New Graph"
2073	Button button6,pos={11,138},size={70,20},proc=Rebin_by,title="Rebin by"
2074	Button button_8,pos={160,138},size={90,20},proc=SaveResultBin,title="Save Result"
2075	Button button_9,pos={20,178},size={30,20},proc=Rebin_Help,title="?"
2076
2077	SetVariable end_3,pos={97,140},size={40,14},title=" ",Font="Arial",fsize=10
2078	SetVariable end_3,limits={1,10,1},value= root:Packages:NIST:Rebin:binning
2079
2080	EndMacro
2081
2082	Function Rebin_Help(ba) : ButtonControl
2083	STRUCT WMButtonAction &ba
2084
2085	String win = ba.win
2086
2087	switch (ba.eventCode)
2088	case 2:
2089	// click code here
2090	DisplayHelpTopic/Z/K=1 "Re-Bin Data"
2091	if(V_flag !=0)
2092	DoAlert 0,"The ReBin Help file could not be found"
2093	endif
2094	break
2095	endswitch
2096
2097	return 0
2098	End
2099
2100	Function Rebin_Done(ctrlName) : ButtonControl
2101	String ctrlName
2102	DoWindow/K Rebin_Panel
2103	End
2104
2105	Proc Rebin_Load(ctrlName) : ButtonControl
2106	String ctrlName
2107
2108	A_LoadOneDData()
2109	ControlUpdate/W=Rebin_Panel popup0
2110	End
2111
2112	Function Rebin_Append(ctrlName) : ButtonControl
2113	String ctrlName
2114	//appends data set from memory /// joindre
2115	String iStr,qStr,eStr
2116	Variable rr,gg,bb
2117
2118	//get the current selection
2119	ControlInfo popup0
2120	if(strlen(S_Value)==0)
2121	Abort "You must load data from a file into memory before appending the data"
2122	Endif
2123
2124	A_PM_doAppend(S_Value)
2125
2126	DoWindow/F Rebin_Panel
2127	End
2128
2129
2130	Function Rebin_newGraph(ctrlName) : ButtonControl
2131	String ctrlName
2132
2133	//get the current selection
2134	ControlInfo popup0
2135	if(strlen(S_Value)==0 \|\| cmpstr(S_Value,"No data loaded")==0)
2136	Abort "You must load data from a file into memory before plotting the data"
2137	Endif
2138
2139	A_PM_doNewGraph(S_Value)
2140
2141	DoWindow/F Rebin_Panel
2142	End
2143
2144
2145
2146	Function Rebin_by (ctrlName) : ButtonControl
2147	String ctrlName
2148
2149
2150	Variable len
2151	Variable ii,jj, helplp,kk
2152	String iStr,qStr,eStr,sqStr,qbStr,fsStr,wtStr,folderStr
2153	Variable rr,gg,bb
2154	NVAR binpts = root:Packages:NIST:Rebin:binning
2155
2156
2157
2158	ControlInfo popup0
2159	if(strlen(S_Value)==0)
2160	Abort "You must load data from a file into memory before plotting the data"
2161	Endif
2162	folderStr = S_Value
2163
2164	SetDataFolder $("root:"+folderStr)
2165
2166	Wave w0 = $(folderStr+"_q")
2167	Wave w1 = $(folderStr+"_i")
2168	Wave w2 = $(folderStr+"_s")
2169	Wave resW = $(folderStr+"_res")
2170
2171	len = numpnts(w0)
2172	Make/O/D/N=(len) w3,w4,w5
2173	w3 = resW[p][0] //std dev of resolution fn
2174	w4 = resW[p][1] //mean q-value
2175	w5 = resW[p][2] //beamstop shadow factor
2176
2177	Make/O/D/N=(round((len-1)/binpts)) qbin,Ibin,sigbin,sqbin,qbbin,fsbin,helpTemp
2178
2179	qbin = 0
2180	Ibin = 0
2181	sigbin = 0
2182	sqbin = 0
2183	qbbin = 0
2184	fsbin = 0
2185	helptemp = 0
2186
2187	ii=0
2188	do
2189	qBin = qBin + w0[binpts*p+ii]/binpts
2190	iBin = iBin + w1[binpts*p+ii]/binpts
2191	helptemp = helptemp + w2[binpts*p+ii]^2
2192	sqBin = sqBin + w3[binpts*p+ii]/binpts
2193	qbBin = qbBin + w4[binpts*p+ii]/binpts
2194	fsBin = fsBin + w5[binpts*p+ii]/binpts
2195
2196
2197	ii+=1
2198	while (ii<binpts)
2199
2200	sigBin = sqrt(helptemp)/binpts
2201
2202	CheckDisplayed Ibin //check top graph only
2203	if(V_flag==0)
2204	AppendToGraph Ibin vs qbin
2205	ModifyGraph log=1,mode(Ibin)=4,marker(Ibin)=29,msize(Ibin)=3,rgb(Ibin)=(65535,0,0)
2206	ErrorBars/T=0 ibin Y,wave=(sigbin,sigbin)
2207	endif
2208
2209	KillWaves/Z helpTemp,w3,w4,w5
2210
2211	SetDataFolder root:
2212	End
2213
2214
2215	// duplicates the binned data to its own data folder, fixing up the names
2216	// as needed, then write out the file from the folder using the standard procedures
2217	// then kill the data folder to clean up.
2218	Function SaveResultBin(ctrlName) : ButtonControl
2219	String ctrlName
2220
2221	String finalname,folderStr
2222
2223	ControlInfo popup0
2224	folderStr=S_Value
2225
2226	SetDataFolder $("root:"+folderStr)
2227	WAVE qbin = qBin
2228	WAVE Ibin = iBin
2229	WAVE sigbin = sigBin
2230	WAVE sqbin = sqBin
2231	WAVE qbbin = qbBin
2232	WAVE fsbin = fsBin
2233
2234	NVAR/Z dQv = USANS_dQv
2235	Variable isSANS = (NVAR_Exists(dQv) == 0)
2236	// Print "isSANS = ",isSANS
2237	if(isSANS)
2238	// SANS data - make a smaller resWave
2239	Variable np=numpnts(qBin)
2240	Make/D/O/N=(np,4) $(folderStr+"_resBin")
2241	Wave res = $(folderStr+"_resBin")
2242
2243	res[][0] = sqBin[p] //sigQ
2244	res[][1] = qbBin[p] //qBar
2245	res[][2] = fsBin[p] //fShad
2246	res[][3] = qBin[p] //Qvalues
2247
2248	endif
2249	finalname = folderStr+"bin"
2250
2251	SetDataFolder root:
2252	DuplicateDataSet(folderStr,finalname,1)
2253
2254	// rename the waves so that the correct ones are written out
2255	SetDataFolder $("root:"+finalname)
2256	Killwaves/Z $(finalname+"_q")
2257	Killwaves/Z $(finalname+"_i")
2258	Killwaves/Z $(finalname+"_s")
2259	if(isSANS) //if USANS, keep the resolution waves, so that the writer will recognize the set as USANS
2260	Killwaves/Z $(finalname+"_res")
2261	endif
2262
2263	WAVE qbin = qBin //reset the wave references to the new folder (bin)
2264	WAVE Ibin = iBin
2265	WAVE sigbin = sigBin
2266	if(isSANS)
2267	Wave res = $(finalName+"_resBin")
2268	endif
2269
2270	Rename qbin $(finalname+"_q")
2271	Rename ibin $(finalname+"_i")
2272	Rename sigbin $(finalname+"_s")
2273	if(isSANS)
2274	Rename res $(finalname+"_res")
2275	endif
2276	SetDataFolder root:
2277
2278	// save out the data in the folder
2279	SaveDataSetToFile(finalName)
2280
2281	KillDataFolder/Z $("root:"+finalname)
2282
2283	SetDataFolder root:
2284	return(0)
2285	End
2286
2287
2288
2289
2290	/// end rebinning procedures

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format