Context Navigation

Wrapper_v40.ipf @ 455

Last change on this file since 455 was 449, checked in by srkline, 14 years ago
Adding the model of a cylinder with spherical end caps. By taking all of the limiting cases, 5 models are added. Since the geometries are sufficiently different, I thought it best to add as separate models. Corresponding changes have been made to the Wrapper and the modelPicker
File size: 38.1 KB

Line
1	#pragma rtGlobals=1 // Use modern global access method.
2	#pragma version=4.00
3	#pragma IgorVersion=6.0
4
5	//
6	// need a way of importing more functions into the experiment
7	// ? call the picker panel from the panel?
8	//
9	//
10
11	//Macro OpenWrapperPanel()
12	// Init_WrapperPanel()
13	//End
14
15	Function Init_WrapperPanel()
16	//make sure that folders exist - this is the first initialization to be called
17	NewDataFolder/O root:Packages
18	NewDataFolder/O root:Packages:NIST
19
20	//Create useful globals
21	Variable/G root:Packages:NIST:SANS_ANA_VERSION=4.00
22	String/G root:Packages:NIST:SANS_ANA_EXTENSION="_v40"
23	//Set this variable to 1 to force use of trapezoidal integration routine for USANS smearing
24	Variable/G root:Packages:NIST:USANSUseTrap = 0
25	Variable/G root:Packages:NIST:USANS_dQv = 0.117
26
27	DoWindow/F WrapperPanel
28	if(V_flag==0)
29	if(exists("root:Packages:NIST:coefKWStr")==0)
30	String/G root:Packages:NIST:coefKWStr=""
31	endif
32	if(exists("root:Packages:NIST:suffixKWStr")==0)
33	String/G root:Packages:NIST:suffixKWStr=""
34	endif
35	Execute "WrapperPanel()"
36	endif
37	End
38
39	////////
40	//
41	// if model is Smeared, search the DF for coefficients
42	// if new DF chosen, need to reset
43	// if new model function, need to reset table
44	//
45	// create hold_mod (0/1), constr_low_mod, constr_hi_mod
46	// in either DF (smeared) or in root: (normal)
47	// and put these in the table as needed
48	//
49	Window WrapperPanel()
50	PauseUpdate; Silent 1 // building window...
51	NewPanel /W=(459,44,1113,499)/N=wrapperPanel/K=1 as "Curve Fit Setup"
52	ModifyPanel fixedSize=1
53
54	GroupBox grpBox_0,pos={18,11},size={390,113}
55	GroupBox grpBox_1,pos={426,10},size={207,113}
56	GroupBox grpBox_2 title="No Fit",pos={10,130},size={0,0},frame=1,fSize=10,fstyle=1,fColor=(39321,1,1)
57	GroupBox grpBox_3 title="",pos={10,150},size={0,0},frame=1,fSize=10,fstyle=1,fColor=(39321,1,1)
58
59	Button button_1,pos={280,57},size={120,20},proc=PlotModelFunction,title="Plot 1D Function"
60	Button button_2,pos={300,93},size={100,20},proc=AppendModelToTarget,title="Append 1D"
61	Button button_3,pos={300,20},size={100,20},proc=W_LoadDataButtonProc,title="Load 1D Data"
62	PopupMenu popup_0,pos={30,21},size={218,20},title="Data Set",proc=DataSet_PopMenuProc
63	PopupMenu popup_0,mode=1,value= #"W_DataSetPopupList()"
64	PopupMenu popup_1,pos={30,57},size={136,20},title="Function"
65	PopupMenu popup_1,mode=1,value= #"W_FunctionPopupList()",proc=Function_PopMenuProc
66	PopupMenu popup_2,pos={30,93},size={123,20},title="Coefficients"
67	PopupMenu popup_2,mode=1,value= #"W_CoefPopupList()",proc=Coef_PopMenuProc
68	CheckBox check_0,pos={440,19},size={79,14},title="Use Cursors?",value= 0
69	CheckBox check_0,proc=UseCursorsWrapperProc
70	CheckBox check_1,pos={440,42},size={74,14},title="Use Epsilon?",value= 0
71	CheckBox check_2,pos={440,65},size={95,14},title="Use Constraints?",value= 0
72	CheckBox check_3,pos={530,18},size={72,14},title="2D Functions?",value= 0
73	CheckBox check_3 proc=Toggle2DControlsCheckProc
74	CheckBox check_4,pos={440,85},size={72,14},title="Report?",value= 0
75	CheckBox check_5,pos={454,103},size={72,14},title="Save it?",value= 0
76	//change draw order to put button over text of checkbox
77	Button button_0,pos={520,93},size={100,20},proc=DoTheFitButton,title="Do 1D Fit"
78	Button button_4,pos={520,126},size={100,20},proc=FeedbackButtonProc,title="Feedback"
79	Button button_5,pos={520,150},size={100,20},proc=FitHelpButtonProc,title="Help"
80
81	Edit/W=(20,174,634,435)/HOST=#
82	ModifyTable width(Point)=0
83	RenameWindow #,T0
84	SetActiveSubwindow ##
85	EndMacro
86
87	//open the Help file for the Fit Manager
88	Function FitHelpButtonProc(ba) : ButtonControl
89	STRUCT WMButtonAction &ba
90
91	switch( ba.eventCode )
92	case 2: // mouse up
93	// click code here
94	DisplayHelpTopic/Z/K=1 "Fit Manager"
95	if(V_flag !=0)
96	DoAlert 0,"The Fit Manager Help file could not be found"
97	endif
98	break
99	endswitch
100
101	return 0
102	End
103
104
105	//open the trac page for feedback
106	Function FeedbackButtonProc(ba) : ButtonControl
107	STRUCT WMButtonAction &ba
108
109	switch( ba.eventCode )
110	case 2: // mouse up
111	// click code here
112	OpenTracTicketPage()
113	break
114	endswitch
115
116	return 0
117	End
118
119
120	//obvious use, now finds the most recent data loaded, finds the folder, and pops the menu with that selection...
121	//
122	Function W_LoadDataButtonProc(ba) : ButtonControl
123	STRUCT WMButtonAction &ba
124
125	switch( ba.eventCode )
126	case 2: // mouse up
127	// click code here
128	String topGraph= WinName(0,1)
129	if(strlen(topGraph) != 0)
130	DoWindow/F $topGraph //so that the panel is not on top
131	endif
132	Execute "A_LoadOneDData()"
133
134	ControlUpdate/W=WrapperPanel popup_0
135	//instead of a simple controlUpdate, better to pop the menu to make sure that the other menus follow
136	// convoluted method to find the right item and pop the menu.
137
138	// data is plotted, so get the "new" top graph
139	topGraph= WinName(0,1) //this is the topmost graph, and should exist, but...
140	if(cmpstr(topGraph,"")==0)
141	return(0)
142	endif
143	String list,folderStr
144	Variable num
145	list = TraceNameList(topGraph,";",1) //want the last item in the list
146	num= ItemsInList(list)
147	FolderStr = StringFromList(num-1,list,";")
148	folderStr = folderStr[0,strlen(folderStr)-3] //remove the "_i" that the loader enforced
149	list = W_DataSetPopupList()
150	num=WhichListItem(folderStr,list,";",0,0)
151	if(num != -1)
152	PopupMenu popup_0,mode=num+1,win=WrapperPanel
153	ControlUpdate/W=WrapperPanel popup_0
154
155	// fake mouse up to pop the menu
156	Struct WMPopupAction ps
157	ps.eventCode = 2 //fake mouse up
158	DataSet_PopMenuProc(ps)
159	endif
160	break
161	endswitch
162
163	return 0
164	End
165
166
167	// is there a simpler way to do this? I don't think so.
168	Function/S W_DataSetPopupList()
169
170	String str=GetAList(4)
171
172	if(strlen(str)==0)
173	str = "No data loaded"
174	endif
175	str = SortList(str)
176
177	return(str)
178	End
179
180
181	// show the available models
182	// not the f*(cw,xw) point calculations
183	// not the *X(cw,xw) XOPS
184	//
185	// KIND:10 should show only user-defined curve fitting functions
186	// - not XOPs
187	// - not other user-defined functions
188	Function/S W_FunctionPopupList()
189	String list,tmp
190	list = FunctionList("*",";","KIND:10") //get every user defined curve fit function
191
192	//now start to remove everything the user doesn't need to see...
193
194	tmp = FunctionList("*_proto",";","KIND:10") //prototypes
195	list = RemoveFromList(tmp, list ,";")
196
197	//prototypes that show up if GF is loaded
198	list = RemoveFromList("GFFitFuncTemplate", list)
199	list = RemoveFromList("GFFitAllAtOnceTemplate", list)
200	list = RemoveFromList("NewGlblFitFunc", list)
201	list = RemoveFromList("NewGlblFitFuncAllAtOnce", list)
202	list = RemoveFromList("GlobalFitFunc", list)
203	list = RemoveFromList("GlobalFitAllAtOnce", list)
204	list = RemoveFromList("GFFitAAOStructTemplate", list)
205	list = RemoveFromList("NewGF_SetXWaveInList", list)
206	list = RemoveFromList("NewGlblFitFuncAAOStruct", list)
207
208	// more to remove as a result of 2D/Gizmo
209	list = RemoveFromList("A_WMRunLessThanDelta", list)
210	list = RemoveFromList("WMFindNaNValue", list)
211	list = RemoveFromList("WM_Make3DBarChartParametricWave", list)
212	list = RemoveFromList("UpdateQxQy2Mat", list)
213	list = RemoveFromList("MakeBSMask", list)
214
215	// MOTOFIT/GenFit bits
216	tmp = "GEN_allatoncefitfunc;GEN_fitfunc;GetCheckBoxesState;MOTO_GFFitAllAtOnceTemplate;MOTO_GFFitFuncTemplate;MOTO_NewGF_SetXWaveInList;MOTO_NewGlblFitFunc;MOTO_NewGlblFitFuncAllAtOnce;"
217	list = RemoveFromList(tmp, list ,";")
218
219	// SANS Reduction bits
220	tmp = "ASStandardFunction;Ann_1D_Graph;Avg_1D_Graph;BStandardFunction;CStandardFunction;Draw_Plot1D;MyMat2XYZ;NewDirection;SANSModelAAO_MCproto;Monte_SANS_Threaded;Monte_SANS_NotThreaded;Monte_SANS_W1;Monte_SANS_W2;Monte_SANS;"
221	list = RemoveFromList(tmp, list ,";")
222
223	// USANS Reduction bits
224	tmp = "DSM_Guinier_Fit;RemoveMaskedPoints;"
225	list = RemoveFromList(tmp, list ,";")
226
227	//more functions from analysis models (2008)
228	tmp = "Barbell_Inner;Barbell_Outer;Barbell_integrand;BCC_Integrand;Integrand_BCC_Inner;Integrand_BCC_Outer;"
229	list = RemoveFromList(tmp, list ,";")
230	tmp = "CapCyl;CapCyl_Inner;CapCyl_Outer;ConvLens;ConvLens_Inner;ConvLens_Outer;"
231	list = RemoveFromList(tmp, list ,";")
232	tmp = "Dumb;Dumb_Inner;Dumb_Outer;FCC_Integrand;Integrand_FCC_Inner;Integrand_FCC_Outer;"
233	list = RemoveFromList(tmp, list ,";")
234	tmp = "Integrand_SC_Inner;Integrand_SC_Outer;SC_Integrand;SphCyl;SphCyl_Inner;SphCyl_Outer;"
235	list = RemoveFromList(tmp, list ,";")
236
237
238	tmp = FunctionList("f*",";","NPARAMS:2") //point calculations
239	list = RemoveFromList(tmp, list ,";")
240
241	tmp = FunctionList("fSmear*",";","NPARAMS:3") //smeared dependency calculations
242	list = RemoveFromList(tmp, list ,";")
243
244	// tmp = FunctionList("*X",";","KIND:4") //XOPs, but these shouldn't show up if KIND:10 is used initially
245	// Print "X* = ",tmp
246	// print " "
247	// list = RemoveFromList(tmp, list ,";")
248
249	//non-fit functions that I can't seem to filter out
250	list = RemoveFromList("BinaryHS_PSF11;BinaryHS_PSF12;BinaryHS_PSF22;EllipCyl_Integrand;PP_Inner;PP_Outer;Phi_EC;TaE_Inner;TaE_Outer;",list,";")
251
252	if(strlen(list)==0)
253	list = "No functions plotted"
254	endif
255
256	list = SortList(list)
257	return(list)
258	End
259
260
261	// show all the appropriate coefficient waves
262	//
263	// also need to search the folder listed in "data set" popup
264	// for smeared coefs
265	//
266	// - or - restrict the coefficient list based on the model function
267	// - new way, filter the possible values based on the data folder and function
268	Function/S W_CoefPopupList()
269
270	String notPlotted="Please plot the function"
271	ControlInfo/W=wrapperpanel popup_1
272	String funcStr=S_Value
273	String coefStr=getFunctionCoef(funcStr)
274
275	if(cmpstr(coefStr,"")==0) //no correspondence in the KW string
276	return(notPlotted)
277	endif
278
279
280	//found a coefficient wave - only two places to look
281	// is it in the root folder?
282	if(exists("root:"+coefStr) != 0)
283	return(coefStr)
284	endif
285
286	// is it in the data folder?
287	ControlInfo/W=wrapperpanel popup_0
288	String folderStr=S_Value
289	if(exists("root:"+folderStr+":"+coefStr) != 0)
290	return(coefStr)
291	endif
292
293	return(notPlotted)
294	End
295
296	// show all the appropriate coefficient waves
297	//
298	// also need to search the folder listed in "data set" popup
299	// for smeared coefs
300	//
301	// - or - restrict the coefficient list based on the model function
302	//
303	// --old way
304	//Function/S W_CoefPopupList()
305	// String list
306	// setDataFolder root:
307	// list = WaveList("coef*",";","")
308	//
309	// ControlInfo/W=wrapperpanel popup_0
310	// if(V_Value != 0) //0== no items in menu
311	// if(DataFolderExists("root:"+S_Value))
312	// SetDataFolder $("root:"+S_Value)
313	// list += WaveList("coef",";","")
314	// endif
315	// endif
316	//
317	// // tmp coefficients that aren't being cleaned up from somewhere...
318	// list = RemoveFromList("temp_coef_1;temp_coef_2;", list ,";")
319	//
320	// if(strlen(list)==0)
321	// list = "No functions plotted"
322	// endif
323	// list = SortList(list)
324	//
325	//// Print itemsinlist(list,";")
326	//
327	// setDataFolder root:
328	// return(list)
329	//End
330
331	// if the coefficients are changed, then update the table
332	//
333	//update the table
334	// may be easier to just kill the subwindow (table) and create a new one
335	// need to set/reset all of the waves in the table
336	//
337	// !! only respond to mouse up
338	//
339	Function Coef_PopMenuProc(pa) : PopupMenuControl
340	STRUCT WMPopupAction &pa
341
342	switch( pa.eventCode )
343	case 2: // mouse up
344	Variable popNum = pa.popNum
345	String popStr = pa.popStr
346	String suffix = getModelSuffix(popStr)
347	ControlInfo/W=WrapperPanel popup_0
348	String folderStr=S_Value
349
350	if(cmpstr(popStr,"Please plot the function")==0)
351	// Print "function not plotted"
352	return(0)
353	endif
354
355	// this if/else/endif should not ever return an error alert
356	// it should simply set the data folder properly
357	if(DataFolderExists("root:"+folderStr))
358	SetDataFolder $("root:"+folderStr)
359	if(!exists(popStr))
360	// must be unsmeared model, work in the root folder
361	SetDataFolder root:
362	if(!exists(popStr)) //this should be fine if the coef filter is working, but check anyhow
363	DoAlert 0,"the coefficient and data sets do not match (1)"
364	return 0
365	endif
366	endif
367	else
368	// must be unsmeared model, work in the root folder
369	SetDataFolder root:
370	if(!exists(popStr)) //this should be fine if the coef filter is working, but check anyhow
371	DoAlert 0,"the coefficient and data sets do not match (2)"
372	return 0
373	endif
374	endif
375
376	// farm the work out to another function?
377	Variable num=numpnts($popStr)
378	// make the necessary waves if they don't exist already
379	if(exists("Hold_"+suffix) == 0)
380	Make/O/D/N=(num) $("epsilon_"+suffix),$("Hold_"+suffix)
381	Make/O/T/N=(num) $("LoLim_"+suffix),$("HiLim_"+suffix)
382	Wave eps = $("epsilon_"+suffix)
383	Wave coef=$popStr
384	if(eps[0] == 0) //if eps already if filled, don't change it
385	eps = abs(coef*1e-4) + 1e-10 //default eps is proportional to the coefficients
386	endif
387	endif
388	// default epsilon values, sometimes needed for the fit
389
390	WAVE/T LoLim = $("LoLim_"+suffix)
391	WAVE/T HiLim = $("HiLim_"+suffix)
392
393	// clear the table (a subwindow)
394	DoWindow/F WrapperPanel // ?? had to add this in during all of the cursor meddling...
395	KillWindow wrapperPanel#T0
396	Edit/W=(20,174,634,435)/HOST=wrapperPanel
397	RenameWindow #,T0
398	// get them onto the table
399	// how do I get the parameter name?
400	String param = WaveList("parameters_"+suffix, "", "TEXT:1," ) //this is hopefully* one wave
401	AppendtoTable/W=wrapperPanel#T0 $param,$(popStr)
402	AppendToTable/W=wrapperPanel#T0 $("Hold_"+suffix),$("LoLim_"+suffix),$("HiLim_"+suffix),$("epsilon_"+suffix)
403	ModifyTable/W=wrapperPanel#T0 width(Point)=0
404
405	SetDataFolder root:
406	break
407	endswitch
408
409	return 0
410	End
411
412	// if the Function is changed, then update the coef popup (if possible) and then the table (if possible)
413	//
414	// !! only respond to mouse up
415	//
416	Function Function_PopMenuProc(pa) : PopupMenuControl
417	STRUCT WMPopupAction &pa
418
419	switch( pa.eventCode )
420	case 2: // mouse up
421	Variable popNum = pa.popNum
422	String funcStr = pa.popStr
423	String coefStr = W_CoefPopupList()
424
425	// Print "coefStr = ",coefStr
426
427	ControlInfo/W=WrapperPanel popup_0
428	String folderStr=S_Value
429
430	String listStr = W_CoefPopupList()
431	Variable num=WhichListItem(coefStr, listStr, ";")
432	String str=StringFromList(num, listStr ,";")
433	// print "str = ",str
434	//set the item in the coef popup, and pop it
435	PopupMenu popup_2 win=WrapperPanel,mode=(num+1)
436
437	Struct WMPopupAction ps
438	ps.eventCode = 2 //fake mouse up
439	ps.popStr = str
440	Coef_PopMenuProc(ps)
441
442	SetDataFolder root:
443	break
444	endswitch
445
446	return 0
447	End
448
449	// if the Data Set is changed, then update the function (if possible)
450	// and the coef popup (if possible) and then the table (if possible)
451	//
452	// !! only respond to mouse up here, and simply send a fake mouse up
453	// to the function pop, which will do what it can do
454	//
455	Function DataSet_PopMenuProc(pa) : PopupMenuControl
456	STRUCT WMPopupAction &pa
457
458	switch( pa.eventCode )
459	case 2: // mouse up
460	// make sure that the cursors are on/off appropriately
461	// let the cursors checkbox decide what to do, sending the current state
462	ControlInfo/W=WrapperPanel check_0
463	STRUCT WMCheckboxAction cba
464	cba.eventCode = 2
465	cba.checked = V_Value
466	UseCursorsWrapperProc(cba)
467
468	// then cascade the function/coefficient popups
469	Struct WMPopupAction ps
470	ps.eventCode = 2 //fake mouse up
471	Function_PopMenuProc(ps)
472
473	SetDataFolder root:
474	break
475	endswitch
476
477	return 0
478	End
479
480
481
482	// always pass this the coefficient string
483	//
484	// either "coef_"
485	// or "smear_coef_"
486	//
487	Function/S getFunctionCoef(funcStr)
488	String funcStr
489
490	SVAR listStr=root:Packages:NIST:coefKWStr
491	String coefStr = StringByKey(funcStr, listStr ,"=",";",0)
492
493	return(coefStr)
494	End
495
496	// always pass this the coefficient string
497	//
498	// either "coef_"
499	// or "smear_coef_"
500	//
501	// does NOT return the leading "_" as part of the suffix
502	Function/S getModelSuffix(coefStr)
503	String coefStr
504
505	SVAR listStr=root:Packages:NIST:suffixKWStr
506	String suffixStr = StringByKey(coefStr, listStr ,"=",";",0)
507
508	return(suffixStr)
509	// Variable pos,start=0
510
511	// if(stringmatch(coefStr,"smear_*") == 1)
512	// start=7 //look forwards to find "_", skipping "smear_coe" if necessary
513	// endif
514	// pos=Strsearch(coefStr,"_",start,0)
515	//Print start, pos
516	// return(coefStr[pos,strlen(coefStr)-1])
517	End
518
519
520
521	// - this is based on stringent "PlotNNN" naming requirements
522	//
523	// ???
524	// - how to kill the generated table and graph, that are not needed now
525	//
526	Function PlotModelFunction(ba) : ButtonControl
527	STRUCT WMButtonAction &ba
528
529	String folderStr,funcStr,coefStr,cmdStr=""
530	Variable useCursors,useEps,useConstr
531
532	Variable killWhat=0 //kill nothing as default
533
534	switch( ba.eventCode )
535	case 2: // mouse up
536	ControlInfo/W=WrapperPanel popup_0
537	folderStr=S_Value
538
539	ControlInfo/W=WrapperPanel popup_1
540	funcStr=S_Value
541
542	// maybe nothing has been loaded yet...
543	if(cmpstr(funcStr,"No functions plotted") == 0)
544	break
545	endif
546
547	// check for smeared or smeared function
548	if(stringmatch(funcStr, "Smear*" )==1)
549	//it's a smeared model
550	// check for the special case of RPA that has an extra parameter
551	if(strsearch(funcStr, "RPAForm", 0 ,0) == -1)
552	sprintf cmdStr, "Plot%s(\"%s\")",funcStr,folderStr //not RPA
553	else
554	sprintf cmdStr, "Plot%s(\"%s\",)",funcStr,folderStr //yes RPA, leave a comma for input
555	endif
556	killWhat = 1
557	else
558	// it's not, don't kill the graph, just the table
559	sprintf cmdStr, "Plot%s()",funcStr
560	killWhat = 2
561	endif
562
563	//Print cmdStr
564	Execute cmdStr
565
566	//pop the function menu to set the proper coefficients
567	DoWindow/F WrapperPanel
568	STRUCT WMPopupAction pa
569	pa.popStr = funcStr
570	pa.eventcode = 2
571	Function_PopMenuProc(pa)
572
573	KillTopGraphAndTable(killWhat) // crude
574
575	break
576	endswitch
577
578	return 0
579	End
580
581	// passing 0 kills nothing
582	// passing 1 kills the top graph and table
583	// passing 2 kills the top table only
584	//
585	Function KillTopGraphAndTable(killwhat)
586	Variable killWhat
587
588	String topGraph= WinName(0,1) //this is the topmost graph
589	String topTable= WinName(0,2) //this is the topmost table
590
591	if(killWhat == 0)
592	return(0)
593	endif
594
595	if(killWhat == 1)
596	KillWindow $topGraph
597	KillWindow $topTable
598	endif
599
600	if(killWhat == 2)
601	KillWindow $topTable
602	endif
603
604	return(0)
605	End
606
607	// How to bypass the step of plot and append?
608	//
609	// do it in two separate events
610	//
611	Function AppendModelToTarget(ba) : ButtonControl
612	STRUCT WMButtonAction &ba
613
614	String coefStr,suffix,yWStr,xWStr,folderStr
615
616	switch( ba.eventCode )
617	case 2: // mouse up
618	ControlInfo/W=WrapperPanel popup_2
619	coefStr=S_Value
620	suffix = getModelSuffix(coefStr)
621
622	// check for smeared or smeared function
623	if(stringmatch(coefStr, "smear*" )==1)
624	//it's a smeared model
625	ControlInfo/W=WrapperPanel popup_0
626	folderStr=S_Value
627	xWStr = "root:"+folderStr+":smeared_qvals"
628	ywStr = "root:"+folderStr+":smeared_"+suffix
629	else
630	// it's not, so it's in the root folder
631	xWStr = "xwave_"+suffix
632	yWStr = "ywave_"+suffix
633	endif
634
635	Wave/Z yw = $yWStr
636	Wave/Z xw = $xWStr
637	if(WaveExists(yw) && WaveExists(xw))
638	AppendtoGraph yw vs xw
639	else
640	DoAlert 0,"The selected model has not been plotted for the selected data set."
641	endif
642	break
643	endswitch
644
645	return 0
646	End
647
648
649	// this should parse the panel and call the FitWrapper() function
650	Function DoTheFitButton(ba) : ButtonControl
651	STRUCT WMButtonAction &ba
652
653	String folderStr,funcStr,coefStr
654	Variable useCursors,useEps,useConstr
655
656	switch( ba.eventCode )
657	case 2: // mouse up
658	ControlInfo/W=WrapperPanel popup_0
659	folderStr=S_Value
660
661	ControlInfo/W=WrapperPanel popup_1
662	funcStr=S_Value
663
664	ControlInfo/W=WrapperPanel popup_2
665	coefStr=S_Value
666
667	ControlInfo/W=WrapperPanel check_0
668	useCursors=V_Value
669	ControlInfo/W=WrapperPanel check_1
670	useEps=V_Value
671	ControlInfo/W=WrapperPanel check_2
672	useConstr=V_Value
673
674	if(!CheckFunctionAndCoef(funcStr,coefStr))
675	DoAlert 0,"The coefficients and function type do not match. Please correct the selections in the popup menus."
676	break
677	endif
678
679	FitWrapper(folderStr,funcStr,coefStr,useCursors,useEps,useConstr)
680
681	// DoUpdate (does not work!)
682	//?? why do I need to force an update ??
683	if(!exists("root:"+folderStr+":"+coefStr))
684	Wave w=$coefStr
685	else
686	Wave w=$("root:"+folderStr+":"+coefStr) //smeared coefs in data folder
687	endif
688	w[0] += 1e-6
689	w[0] -= 1e-6
690
691	break
692	endswitch
693
694	return 0
695	End
696
697	Function CheckFunctionAndCoef(funcStr,coefStr)
698	String funcStr,coefStr
699
700	SVAR listStr=root:Packages:NIST:coefKWStr
701	String properCoefStr = StringByKey(funcStr, listStr ,"=",";",0)
702	if(cmpstr("",properCoefStr)==0)
703	return(0) //false, no match found, so properCoefStr is returned null
704	endif
705	if(cmpstr(coefStr,properCoefStr)==0)
706	return(1) //true, the coef is the correct match
707	endif
708	return(0) //false, wrong coef
709	End
710
711	/////////////////////////////////
712
713	// wrapper to do the desired fit
714	//
715	// folderStr is the data folder for the desired data set
716	//
717	//
718	Function FitWrapper(folderStr,funcStr,coefStr,useCursors,useEps,useConstr)
719	String folderStr,funcStr,coefStr
720	Variable useCursors,useEps,useConstr
721
722	String suffix=getModelSuffix(coefStr)
723
724	SetDataFolder $("root:"+folderStr)
725	if(!exists(coefStr))
726	// must be unsmeared model, work in the root folder
727	SetDataFolder root:
728	if(!exists(coefStr)) //this should be fine if the coef filter is working, but check anyhow
729	DoAlert 0,"the coefficient and data sets do not match"
730	return 0
731	endif
732	endif
733
734	WAVE cw=$(coefStr)
735	Wave hold=$("Hold_"+suffix)
736	Wave/T lolim=$("LoLim_"+suffix)
737	Wave/T hilim=$("HiLim_"+suffix)
738	Wave eps=$("epsilon_"+suffix)
739
740	// fill a struct instance whether I need one or not
741	String DF="root:"+folderStr+":"
742
743	Struct ResSmearAAOStruct fs
744	WAVE/Z resW = $(DF+folderStr+"_res") //these may not exist, if 3-column data is used
745	WAVE/Z fs.resW = resW
746	WAVE yw=$(DF+folderStr+"_i")
747	WAVE xw=$(DF+folderStr+"_q")
748	WAVE sw=$(DF+folderStr+"_s")
749	Wave fs.coefW = cw
750	Wave fs.yW = yw
751	Wave fs.xW = xw
752
753	Duplicate/O yw $(DF+"FitYw")
754	WAVE fitYw = $(DF+"FitYw")
755	fitYw = NaN
756
757	Variable useRes=0
758	if(stringmatch(funcStr, "Smear*")) // if it's a smeared function, need a struct
759	useRes=1
760	endif
761
762	// do not construct constraints for any of the coefficients that are being held
763	// -- this will generate an "unknown error" from the curve fitting
764	Make/O/T/N=0 constr
765	if(useConstr)
766	String constraintExpression
767	Variable i, nPnts=DimSize(lolim, 0),nextRow=0
768	for (i=0; i < nPnts; i += 1)
769	if (strlen(lolim[i]) > 0 && hold[i] == 0)
770	InsertPoints nextRow, 1, constr
771	sprintf constraintExpression, "K%d > %s", i, lolim[i]
772	constr[nextRow] = constraintExpression
773	nextRow += 1
774	endif
775	if (strlen(hilim[i]) > 0 && hold[i] == 0)
776	InsertPoints nextRow, 1, constr
777	sprintf constraintExpression, "K%d < %s", i, hilim[i]
778	constr[nextRow] = constraintExpression
779	nextRow += 1
780	endif
781	endfor
782	endif
783
784	// 20JUN if useCursors is true, and there are no cursors on the specified data set, uncheck and set to false
785	// this is a last line of defense, and should never actually do anything...
786	if(useCursors)
787	useCursors = AreCursorsCorrect(folderStr)
788	endif
789	//if useCursors, and the data is USANS, need to recalculate the matrix if the range is new
790	Variable pt1,pt2,newN,mPt1,mPt2
791	String noteStr
792	if(useCursors && (dimsize(resW,1) > 4) )
793	//where are the cursors, and what is the status of the current matrix?
794	if(pcsr(A) > pcsr(B))
795	pt1 = pcsr(B)
796	pt2 = pcsr(A)
797	else
798	pt1 = pcsr(A)
799	pt2 = pcsr(B)
800	endif
801
802	noteStr = note(resW)
803	mPt1 = NumberByKey("P1",noteStr,"=",";")
804	mPt2 = NumberByKey("P2",noteStr,"=",";")
805	if((mPt1 != pt1) \|\| (mPt2 != pt2) )
806	// need to recalculate
807	USANS_RE_CalcWeights(folderStr,pt1,pt2)
808	Print "Done recalculating the matrix"
809	endif
810
811	Wave trimResW=$(DF+folderStr+"_res"+"t") //put the trimmed resW in the struct for the fit!
812	Wave fs.resW=trimResW
813
814	endif
815
816	// create these variables so that FuncFit will set them on exit
817	Variable/G V_FitError=0 //0=no err, 1=error,(2^1+2^0)=3=singular matrix
818	Variable/G V_FitQuitReason=0 //0=ok,1=maxiter,2=user stop,3=no chisq decrease
819
820
821	// don't use the auto-destination with no flag, it doesn't appear to work correctly
822	// dispatch the fit
823
824	// currently, none of the fit functions are defined as threadsafe, so I don't think that the /NTHR flag really
825	// does anything. The functions themselves can be threaded since they are AAO, and that is probably enough,
826	// since it doesn't make much sense to thread threads. In addition, there is a little-publicized warning
827	// in the WM help file that /C=texWave cannot be used to specify constraints for threadsafe functions!
828	// The textwave would have to be parsed into a constraint matrix first, then passed as /C={cMat,cVec}.
829	// -- just something to watch out for.
830
831	do
832	// Variable t0 = stopMStimer(-2) // corresponding print is at the end of the do-while loop (outside)
833
834	if(useRes && useEps && useCursors && useConstr) //do it all
835	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw[pcsr(A),pcsr(B)] /X=xw /W=sw /I=1 /E=eps /D=fitYw /C=constr /STRC=fs
836	break
837	endif
838
839	if(useRes && useEps && useCursors) //no constr
840	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw[pcsr(A),pcsr(B)] /X=xw /W=sw /I=1 /E=eps /D=fitYw /STRC=fs
841	break
842	endif
843
844	if(useRes && useEps && useConstr) //no crsr
845	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /E=eps /D=fitYw /C=constr /STRC=fs
846	break
847	endif
848
849	if(useRes && useCursors && useConstr) //no eps
850	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw[pcsr(A),pcsr(B)] /X=xw /W=sw /I=1 /D=fitYw /C=constr /STRC=fs
851	break
852	endif
853
854	if(useRes && useCursors) //no eps, no constr
855	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw[pcsr(A),pcsr(B)] /X=xw /W=sw /I=1 /D=fitYw /STRC=fs
856	break
857	endif
858
859	if(useRes && useEps) //no crsr, no constr
860	// Print "timing test for Cylinder_PolyRadius --- the supposedly threaded version ---"
861	// Variable t0 = stopMStimer(-2)
862	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /E=eps /D=fitYw /STRC=fs
863	// t0 = (stopMSTimer(-2) - t0)*1e-6
864	// Printf "CylPolyRad fit time using res and eps and /NTHR=0 time = %g seconds\r\r",t0
865
866	// cw[0] = .01
867	// cw[1] = 20
868	// cw[2] = 400
869	// cw[3] = 0.2
870	// cw[4] = 3e-6
871	// cw[5] = 0.0
872	//
873	// t0 = stopMSTimer(-2)
874	// FuncFit/H=getHStr(hold) $funcStr cw, yw /X=xw /W=sw /I=1 /E=eps /D=fitYw /STRC=fs
875	// t0 = (stopMSTimer(-2) - t0)*1e-6
876	// Printf "CylPolyRad fit time using res and eps and NO THREADING time = %g seconds\r\r",t0
877	break
878	endif
879
880	if(useRes && useConstr) //no crsr, no eps
881	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /D=fitYw /C=constr /STRC=fs
882	break
883	endif
884
885	if(useRes) //just res
886	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /D=fitYw /STRC=fs
887	break
888	endif
889
890	///// same as above, but all without useRes (no /STRC flag)
891	if(useEps && useCursors && useConstr) //do it all
892	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw[pcsr(A),pcsr(B)] /X=xw /W=sw /I=1 /E=eps /D=fitYw /C=constr
893	break
894	endif
895
896	if(useEps && useCursors) //no constr
897	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw[pcsr(A),pcsr(B)] /X=xw /W=sw /I=1 /E=eps /D=fitYw
898	break
899	endif
900
901	if(useEps && useConstr) //no crsr
902	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /E=eps /D=fitYw /C=constr
903	break
904	endif
905
906	if(useCursors && useConstr) //no eps
907	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw[pcsr(A),pcsr(B)] /X=xw /W=sw /I=1 /D=fitYw /C=constr
908	break
909	endif
910
911	if(useCursors) //no eps, no constr
912	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw[pcsr(A),pcsr(B)] /X=xw /W=sw /I=1 /D=fitYw
913	break
914	endif
915
916	if(useEps) //no crsr, no constr
917	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /E=eps /D=fitYw
918	break
919	endif
920
921	if(useConstr) //no crsr, no eps
922	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /D=fitYw /C=constr
923	break
924	endif
925
926	//just a plain vanilla fit
927	FuncFit/H=getHStr(hold) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /D=fitYw
928
929	while(0)
930
931	// t0 = (stopMSTimer(-2) - t0)*1e-6
932	// Printf "fit time = %g seconds\r\r",t0
933
934	// append the fit
935	// need to manage duplicate copies
936	// Don't plot the full curve if cursors were used (set fitYw to NaN on entry...)
937	String traces=TraceNameList("", ";", 1 ) //"" as first parameter == look on the target graph
938	if(strsearch(traces,"FitYw",0) == -1)
939	AppendToGraph FitYw vs xw
940	else
941	RemoveFromGraph FitYw
942	AppendToGraph FitYw vs xw
943	endif
944	ModifyGraph lsize(FitYw)=2,rgb(FitYw)=(0,0,0)
945
946	DoUpdate //force update of table and graph with fitted values (why doesn't this work? - the table still does not update)
947
948	// report the results (to the panel?)
949	print "V_chisq = ",V_chisq
950	print cw
951	WAVE/Z w_sigma
952	print w_sigma
953	String resultStr=""
954
955	//now re-write the results
956	sprintf resultStr,"Chi^2 = %g Sqrt(X^2/N) = %g",V_chisq,sqrt(V_chisq/V_Npnts)
957	resultStr = PadString(resultStr,63,0x20)
958	GroupBox grpBox_2 title=resultStr
959	ControlUpdate/W=WrapperPanel grpBox_2
960	sprintf resultStr,"FitErr = %s : FitQuit = %s",W_ErrorMessage(V_FitError),W_QuitMessage(V_FitQuitReason)
961	resultStr = PadString(resultStr,63,0x20)
962	GroupBox grpBox_3 title=resultStr
963	ControlUpdate/W=WrapperPanel grpBox_3
964
965	Variable yesSave=0,yesReport=0
966	ControlInfo/W=WrapperPanel check_4
967	yesReport = V_Value
968	ControlInfo/W=WrapperPanel check_5
969	yesSave = V_Value
970
971
972	if(yesReport)
973	String parStr=GetWavesDataFolder(cw,1)+ WaveList("param"+suffix, "", "TEXT:1," ) //this is hopefully one wave
974	String topGraph= WinName(0,1) //this is the topmost graph
975
976	DoUpdate //force an update of the graph before making a copy of it for the report
977
978	W_GenerateReport(funcStr,folderStr,$parStr,cw,yesSave,V_chisq,W_sigma,V_npnts,V_FitError,V_FitQuitReason,V_startRow,V_endRow,topGraph)
979	endif
980
981	SetDataFolder root:
982	return(0)
983	End
984
985	// parse something off of a table, or ?
986	Function/S getHStr(hold)
987	Wave hold
988
989	String str=""
990	Variable ii
991	for(ii=0;ii<numpnts(hold);ii+=1)
992	str += num2str(hold[ii])
993	endfor
994
995	// print str
996	if(strsearch(str, "1", 0) == -1)
997	return ("")
998	else
999	return(str)
1000	endif
1001	End
1002
1003	//taken from SRK Auto_Fit, and modified to work better with FitWrapper
1004	//must have AutoGraph as the name of the graph window (any size)
1005	// func is the name of the function (for print only)
1006	//par and coef are the exact names of the waves
1007	//yesSave==1 will save the file(name=func+time)
1008	//
1009	Function W_GenerateReport(func,dataname,param,ans,yesSave,chiSq,sigWave,npts,fitErr,fitQuit,fitStart,fitEnd,topGraph)
1010	String func,dataname
1011	Wave/T param
1012	Wave ans
1013	Variable yesSave,chiSq
1014	Wave sigWave
1015	Variable npts,fitErr,fitQuit,fitStart,fitEnd
1016	String topGraph
1017
1018	String str,pictStr="P_"
1019	String nb="Report"
1020
1021	// bring report up
1022	DoWindow/F Report
1023	if (V_flag == 0) // Report notebook doesn't exist ?
1024	NewNotebook/W=(10,45,550,620)/F=1/N=Report as "Report"
1025	endif
1026	// delete old stuff
1027	Notebook $nb selection={startOfFile, endOfFile}, text="\r", selection={startOfFile, startOfFile}
1028
1029	//setup
1030	Notebook $nb defaultTab=36, statusWidth=252, pageMargins={72,72,72,72}
1031	Notebook $nb showRuler=1, rulerUnits=1, updating={1, 60}
1032	Notebook $nb newRuler=Normal, justification=0, margins={0,0,468}, spacing={0,0,0}, tabs={}, rulerDefaults={"Geneva",10,0,(0,0,0)}
1033	//
1034	// insert title
1035	Notebook $nb newRuler=Title, justification=1, rulerDefaults={"Times", 16, 1, (0, 0, 0)}
1036	sprintf str, "Fit to %s, %s, %s\r\r", func,Secs2Date(datetime, 0), time()
1037	Notebook $nb ruler=Title, text=str
1038
1039	// insert fit results
1040	Variable num=numpnts(ans),ii=0
1041	Notebook $nb ruler=Normal
1042	Notebook $nb margins={18,18,504}, tabs={63 + 3*8192}
1043	str = "Data file: " + dataname + "\r\r"
1044	Notebook $nb text=str
1045	Notebook $nb ruler=Normal
1046	Notebook $nb margins={18,18,504}, tabs={144,234,252}
1047	do
1048	sprintf str, "%s = \t%g\t±\t%g\r", param[ii],ans[ii],sigwave[ii]
1049	Notebook $nb text=str
1050	ii+=1
1051	while(ii<num)
1052
1053	//
1054	// no "fitted range" for 2D data, so make sure that this exists before using it
1055	Wave/Z dataXw = $("root:"+dataname+":"+dataname+"_q")
1056	//
1057	Notebook $nb ruler=Normal
1058	Notebook $nb margins={18,18,504}, tabs={63+3*8192}, fStyle=1, textRGB=(65000,0,0)
1059
1060	sprintf str,"chisq = %g\r",chisq
1061	Notebook $nb textRGB=(65000,0,0),fstyle=1,text=str
1062	sprintf str,"Npnts = %g \t\t Sqrt(X^2/N) = %g\r",npts,sqrt(chiSq/npts)
1063	Notebook $nb textRGB=(0,0,0),fstyle=0, text=str
1064	if(WaveExists(dataXw))
1065	sprintf str "Fitted range = [%d,%d] = %g < Q < %g\r",fitStart,fitEnd,dataXw(fitStart),dataXw(fitEnd)
1066	Notebook $nb textRGB=(0,0,0),fstyle=0, text=str
1067	endif
1068	sprintf str,"FitError = %s\t\tFitQuitReason = %s\r",W_ErrorMessage(FitErr),W_QuitMessage(FitQuit)
1069	Notebook $nb textRGB=(65000,0,0),fstyle=1,text=str
1070	Notebook $nb ruler=Normal
1071
1072	// insert graphs
1073	if(WaveExists(dataXw))
1074	Notebook $nb picture={$topGraph(0, 0, 400, 300), -5, 1}, text="\r"
1075	//
1076	else //must be 2D Gizmo
1077	Execute "ExportGizmo Clip" //this ALWAYS is a PICT or BMP. Gizmo windows are different...
1078	LoadPict/Q/O "Clipboard",tmp_Gizmo
1079	Notebook $nb picture={tmp_Gizmo(0, 0, 400, 300), 0, 1}, text="\r"
1080	endif
1081	//Notebook Report picture={Table1, 0, 0}, text="\r"
1082
1083	// show the top of the report
1084	Notebook $nb selection= {startOfFile, startOfFile}, findText={"", 1}
1085
1086	//save the notebook and the graphic file
1087	if(yesSave)
1088	String nameStr=CleanupName(func,0)
1089	nameStr = nameStr[0,8] //shorten the name
1090	nameStr += "_"+dataname
1091	//make sure the name is no more than 31 characters
1092	namestr = namestr[0,30] //if shorter than 31, this will NOT pad to 31 characters
1093	Print "file saved as ",nameStr
1094	SaveNotebook /O/P=home/S=2 $nb as nameStr
1095	//save the graph separately as a PICT file, 2x screen
1096	pictStr += nameStr
1097	pictStr = pictStr[0,28] //need a shorter name - why?
1098	// DoWindow/F $topGraph
1099	// E=-5 is png @screen resolution
1100	// E=2 is PICT @2x screen resolution
1101	/// SavePICT /E=-5/O/P=home /I/W=(0,0,3,3) as pictStr
1102	if(WaveExists(dataXw))
1103	SavePICT /E=-5/O/P=home/WIN=$topGraph /W=(0,0,400,300) as pictStr
1104	else
1105	Execute "ExportGizmo /P=home as \""+pictStr+"\""
1106	//SavePICT /E=-5/O/P=home/WIN=$topGraph /W=(0,0,400,300) as pictStr
1107	endif
1108	Endif
1109
1110	// ???maybe print the notebook too?
1111	End
1112
1113	Function/S W_ErrorMessage(code)
1114	Variable code
1115
1116	switch (code)
1117	case 0:
1118	return "No Error"
1119	break
1120	case 3: //2^0 + 2^1
1121	return "Singular Matrix"
1122	break
1123	case 5: //(2^0 + 2^2)
1124	return "Out of Memory"
1125	break
1126	case 9: //(2^0 + 2^3)
1127	return "Func= NaN or Inf"
1128	break
1129	default:
1130	return "Unknown error code "+num2str(code)
1131	endswitch
1132	end
1133
1134	Function/S W_QuitMessage(code)
1135	Variable code
1136
1137	switch (code)
1138	case 0:
1139	return "No Error"
1140	break
1141	case 1:
1142	return "Max iterations - re-run fit"
1143	break
1144	case 2:
1145	return "User terminated fit"
1146	break
1147	case 3:
1148	return "No decrease in chi-squared"
1149	break
1150	default:
1151	return "Unknown Quit code "+num2str(code)
1152	endswitch
1153	end
1154
1155	Function Toggle2DControlsCheckProc(cba) : CheckBoxControl
1156	STRUCT WMCheckboxAction &cba
1157
1158	switch( cba.eventCode )
1159	case 2: // mouse up
1160	Variable checked = cba.checked
1161	if(checked)
1162	//print "change the buttons to 2D"
1163	Button button_0,proc=Do2DFitButtonProc,title="Do 2D Fit"
1164	Button button_1,size={120,20},proc=Plot2DFunctionButtonProc,title="Plot 2D Function"
1165	Button button_2,size={100,20},proc=Append2DModelButtonProc,title="Append 2D"
1166	Button button_3,size={100,20},proc=Load2DDataButtonProc,title="Load 2D Data"
1167
1168	Button button_2D_0,pos={550,60},size={70,20},proc=LogToggle2DButtonProc,title="Log/Lin"
1169	Button button_2D_1,pos={520,37},size={100,20},proc=Plot2DButtonProc,title="Plot 2D Data"
1170
1171	Button button_2D_0,disable=0 //visible again, and enabled
1172	Button button_2D_1,disable=0
1173	else
1174	//print "unchecked, change them back to 1D"
1175	Button button_0,pos={520,93},size={100,20},proc=DoTheFitButton,title="Do 1D Fit"
1176	Button button_1,pos={280,57},size={120,20},proc=PlotModelFunction,title="Plot 1D Function"
1177	Button button_2,pos={300,93},size={100,20},proc=AppendModelToTarget,title="Append 1D"
1178	Button button_3,pos={300,20},size={100,20},proc=W_LoadDataButtonProc,title="Load 1D Data"
1179
1180	Button button_2D_0,disable=3 //hide the extra 2D buttons, and disable
1181	Button button_2D_1,disable=3
1182	endif
1183	break
1184	endswitch
1185
1186	return 0
1187	End
1188
1189	// function to either add or remove the cursors from the topmost graph, as needed
1190
1191	Function UseCursorsWrapperProc(cba) : CheckBoxControl
1192	STRUCT WMCheckboxAction &cba
1193
1194
1195	switch( cba.eventCode )
1196	case 2: // mouse up
1197
1198	// check to make sure there really is a "topmost" graph
1199	String topGraph= WinName(0,1) //this is the topmost graph
1200	if(cmpstr(topGraph,"")==0) //no graphs, uncheck and exit
1201	CheckBox check_0,value=0
1202	return(0)
1203	endif
1204
1205	String ciStr = CsrInfo(A , topGraph)
1206
1207	ControlInfo/W=wrapperpanel popup_0
1208	String folderStr=S_Value
1209	String traceList = TraceNameList(topGraph, ";", 1 )
1210
1211	Variable checked = cba.checked
1212
1213	if(checked)
1214	//print "add the cursors to the topmost graph, if the data set is there, or move them"
1215	if(strlen(ciStr)==0 && strsearch(traceList, folderStr, 0) != -1 ) //no cursors && data is there
1216	ShowInfo
1217	Wave yw=$("root:"+folderStr+":"+folderStr+"_i")
1218	Cursor/P/W=$topGraph A, $(folderStr+"_i"),0
1219	Cursor/P/W=$topGraph/A=0 B, $(folderStr+"_i"),numpnts(yw)-1 //deactivate the one at the high Q end
1220	DoUpdate
1221	elseif (strlen(ciStr)!=0 && strsearch(traceList, folderStr, 0) != -1 ) //cursors present, but on wrong data
1222	Wave yw=$("root:"+folderStr+":"+folderStr+"_i")
1223	Cursor/P/W=$topGraph A, $(folderStr+"_i"),0 //move the cursors
1224	Cursor/P/W=$topGraph/A=0 B, $(folderStr+"_i"),numpnts(yw)-1
1225	DoUpdate
1226	endif
1227
1228	AreCursorsCorrect(folderStr)
1229	else
1230	//print "unchecked, remove the cursors"
1231	// go back to the full matrix for the resolution calculation (not if SANS data...)
1232	if(waveExists($("root:"+folderStr+":weights_save")))
1233	Duplicate/O $("root:"+folderStr+":weights_save"), $("root:"+folderStr+":"+folderStr+"_res"),$("root:"+folderStr+":"+folderStr+"_rest")
1234	endif
1235
1236	HideInfo
1237	Cursor/K A
1238	Cursor/K B
1239	endif
1240	break
1241	endswitch
1242
1243	return 0
1244	End
1245
1246	// returns 1 if the specified data is on the top graph && has cursors on it
1247	// returns 0 and unchecks the box if anything is wrong
1248	//
1249	// only call this function if the cursor box is checked, to uncheck it as needed
1250	// if the box is unchecked, leave it be.
1251	//
1252	Function AreCursorsCorrect(folderStr)
1253	String folderStr
1254
1255	String topGraph= WinName(0,1) //this is the topmost graph
1256	if(cmpstr(topGraph,"")==0) //no graphs, uncheck and exit
1257	CheckBox check_0,win=wrapperpanel,value=0
1258	return(0)
1259	endif
1260
1261	String traceAisOn = CsrWave(A , "", 0)
1262	if( strsearch(traceAisOn, folderStr, 0) == -1) //data and cursors don't match
1263	CheckBox check_0,win=wrapperpanel,value=0
1264	HideInfo
1265	Cursor/K A
1266	Cursor/K B
1267	return(0)
1268	endif
1269
1270	return(1)
1271	End

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

Download in other formats:

Original Format