Context Navigation

Auto_Fit.ipf @ 920

Last change on this file since 920 was 920, checked in by srkline, 10 years ago
Adding to batch fitting. Now automatically saves the fit XY waves and coefficients as separate text files that can be loaded and appended later, or imported into a different program.
File size: 62.9 KB

Line
1	#pragma rtGlobals=1 // Use modern global access method.
2
3
4	// This was originally written 2001-2003 ish. This works very differently than SANSview, and I haven't
5	// mentioned that this exists, so that their ideas will be new. Still, some of this may be
6	// serviceable for use within the Igor package, with some cleanup of the interface and functionality.
7	//
8	//
9	//
10	// to do - July 2010
11	//
12	// X- smeared models don't work at all (now they work)
13	// x- the report is good, but the parsing is not great - the format
14	// could be more convenient to read back in to a table (text col, num cols (matrix that can be plotted?)
15	// x- Maybe the report could be compiled as waves. Then they would be numeric. Nothing to plot them against
16	// other than file names, but with a table of values and names, one could manually enter a meaningful
17	// column of values to plot against.
18	// x- Make_HoldConstraintEps_Waves() is currently disabled, both on the macros menu and the checkbox on the panel
19	// Constraints and holds are part of the matrix already, so this would only be for epsilon
20	// x- Many of the items on the Macros menu are just junk. Remove them.
21	//
22	// x- savePath is obsolete. Remove the button. Report saving is sent to the procedure defined in Wrapper.ipf, and uses the home path
23	// x- add a checkbox for "all" to select all data files (function is already set -- CheckAllFiles())
24	//
25	// x- make it easy for the user to save the results in a whole table.
26	// x- remind to rename the results (or move to a folder?) before another batch is run
27	// x- what to plot against?
28	//
29	// x- add a little text box of explanation on each of the tabs at the bottom - to give simple
30	// instructions of what to do on each tab. This will need to change with each tab.
31	//
32	// x- need a help button to link to the help file
33	// -- Comprehensive instructions (after features are finalized)
34	//
35	// -- add a way to save the setup and guesses. May be tough if the file list changes...
36
37
38
39	Menu "Macros"
40	SubMenu "Auto Fit"
41	"Initialize AutoFit Panel"
42	"-"
43	// "Generate_Data_Checklist"
44	// "Make_HoldConstraintEps_Waves" //this is currently disabled
45	// "-"
46	// "CompileAndLoadReports"
47	"Compile_Reports"
48	"Load Compiled Reports"
49	// "Compile_GlobalFit_Reports"
50	"-"
51	"Close Open Notebooks"
52	"Print Open Notebooks"
53	"-"
54	"Generate_PNG_list"
55	"Layout_PNGs"
56
57	End
58	End
59
60
61	//***********
62	//Automatic fitting routines for SANS or NSE data
63	//requires model functions and their associated files
64	//
65	//**************
66
67
68
69	//path to the data
70	// ? have this automatically do the listing?
71	Function PDPButton(ctrlName) : ButtonControl
72	String ctrlName
73
74	//set the global string to the selected pathname
75	NewPath/O/M="pick the data folder" dataPath
76	PathInfo/S dataPath
77	Return(0)
78	End
79
80	Function SavePathButtonProc(ctrlName) : ButtonControl
81	String ctrlName
82
83	//set the global string to the selected pathname
84	NewPath/O/M="pick the save folder" savePath
85	PathInfo/S savePath
86	Return(0)
87	End
88
89	//
90	// this re-sizes AND re-initializes the guessMatrix
91	// and is only done when the files are re-listed, using the list button - since the
92	// new list will screw up the correspondence with the guess matrix
93	//
94	Function FileListButtonProc(ctrlName) : ButtonControl
95	String ctrlName
96
97	DoAlert 1,"This will re-initialize the file list and ALL of the entries in the guess matrix. Do you want to proceed?"
98	if(V_flag!=1) //if the answer is not yes, get out
99	return(0)
100	endif
101
102	String fileList=""
103	fileList = IndexedFile(dataPath,-1,"????")
104	fileList = RemoveFromList(".DS_Store", fileList, ";",0)
105	List2TextWave(Filelist,";","root:AutoFit:fileListWave")
106	Wave/T lw = $"root:AutoFit:fileListWave"
107	Redimension/N=(-1,2) lw
108	lw[][1] = ""
109	Variable num = DimSize(lw,0)
110
111	Wave sw = $"root:AutoFit:fileSelWave"
112	Redimension/B/N=(num,2) sw
113	sw[][1] = (sw[p][1] == 0) ? 2^5 : sw[p][1] //if the element is zero, make it 32, otherwise leave it alone
114
115	// force a re-size (and re-initialization) of the guess-hold-constr-range matrices
116	ToMatrixButtonProc("")
117
118	return(0)
119	End
120
121
122	//currently set up to toggle state, based on the first row
123	//
124	Function CheckAllFiles(ctrlName,checked) : CheckBoxControl
125	String ctrlName
126	Variable checked
127
128	WAVE sel = root:AutoFit:fileSelWave
129	if( (sel[0][1] & 0x10) == 16 ) //first box is checked
130	sel[][1] = 2^5 // uncheck them all
131	else
132	sel[][1] = 2^5 + 2^4 //bit 5== checkbox, bit 4== checked
133	endif
134	return(0)
135	End
136
137
138	//
139	// this re-sizes AND re-initializes the guessMatrix
140	// and is only done when the files are re-listed, using the list button - since the
141	// new list will screw up the correspondence with the guess matrix
142	//
143	Function ToMatrixButtonProc(ctrlName) : ButtonControl
144	String ctrlName
145
146	SetDataFolder root:AutoFit
147	WAVE sel = fileSelWave
148	NVAR numPar=numPar
149	Variable jj,numSelFiles=0
150
151	//update the number of steps for the progress bar
152	// commented out - SRK - 28MAR05
153	// NVAR cur=curProgress
154	// NVAR totSteps=endProgress
155	// cur=0
156	// totSteps=numSelFiles*3 //3 fits per file
157	// ValDisplay progressValdisp,limits={0,totSteps,0}
158	// ControlUpdate/W=AutoFitPanel progressValdisp
159
160	//redimension the matrices and lists
161	Redimension/N=(numPar,DimSize(sel,0)) guessMatrix,HoldMatrix,constrMatrix,rangeMatrix
162	Redimension/N=(numPar) guessList,holdList,constrList,rangeList,guessSel,holdSel,constrSel,rangeSel
163
164	//Set the matrices to zero
165	WAVE/T aa = guessMatrix
166	WAVE/T bb = holdMatrix
167	WAVE/T cc = constrMatrix
168	WAVE/T rr = rangeMatrix
169	aa="1"
170	bb="0"
171	cc="" //set to no constraints
172	rr="0" // zero in points 0 and 1 force a fit of all of the data (default)
173	//set the Lists to "0"
174	WAVE/T dd=guessList
175	WAVE/T ee=holdList
176	WAVE/T ff=constrList
177	WAVE/T rrl=rangeList
178	dd="0"
179	ee="0"
180	ff="0"
181	rrl="0"
182	//Set the sel waves to 2 (editable)
183	WAVE gg=guessSel
184	WAVE hh=holdSel
185	WAVE ii=constrSel
186	WAVE rrs=rangeSel
187	gg=2
188	hh=2
189	ii=2
190	rrs=2
191
192	SetDataFolder root:
193
194	DisplayGuess()
195	tabProc("",0) //force the Guess tab to top
196	TabControl tabC,value=0
197
198	return(0)
199	End
200
201
202	//updates the lists (text waves) to display, depending on the current highlighted file
203	Function DisplayGuess()
204
205	//find the selected data file -[][0] is the file, not the checkbox
206	Variable row
207	row = FindSelectedRowInLB(0) // send 0 to look at the column with filenames
208
209	//update the text waves based on the values in the matrix
210	WAVE/T gm = root:AutoFit:guessMatrix
211	WAVE/T gl = root:AutoFit:guessList
212
213	// the hold matrix
214	WAVE/T hm=root:AutoFit:holdMatrix
215	WAVE/T hl=root:AutoFit:holdList
216
217	// the constraint matrix
218	WAVE/T cm=root:AutoFit:constrMatrix
219	WAVE/T cl=root:AutoFit:constrList
220
221	// the range matrix
222	WAVE/T rm=root:AutoFit:rangeMatrix
223	WAVE/T rl=root:AutoFit:rangeList
224
225	WAVE sel = root:AutoFit:fileSelWave
226	if( (sel[row][1] & 0x10) == 16 ) //box is checked
227
228	gl = gm[p][row]
229	hl = hm[p][row]
230	cl = cm[p][row]
231	rl = rm[p][row]
232	else
233	//file is not checked, don't display guess
234	gl = "Not checked"
235	hl = "Not checked"
236	cl = "Not checked"
237	rl = "Not checked"
238	endif
239
240	return(0)
241	End
242
243
244	//column 0 is the filenames, 1 is the checkboxes
245	//returns the row that is selected, -1 if not found
246	Function FindSelectedRowInLB(col)
247	Variable col
248
249	WAVE sel = root:AutoFit:fileSelWave
250	Variable jj=0,row=-1
251	if(col==0)
252	do
253	if( (sel[jj][0] & 0x01) == 1 ) //returns 1 (bit 0) if selected, 0 otherwise
254	row=jj
255	//print "Row = ",jj
256	break
257	endif
258	jj+=1
259	while(jj<DimSize(sel,0))
260	else
261	do
262	if( (sel[jj][1] & 0x01) == 1 ) //the checkbox is selected, but not necessarily checked!
263	row=jj
264	//print "Row = ",jj
265	break
266	endif
267	jj+=1
268	while(jj<DimSize(sel,0))
269	endif
270
271	return(row)
272	end
273
274
275	Proc InitializeAutoFitPanel()
276	if (wintype("AutoFitPanel") == 0)
277	InitializeAutoFit()
278	AutoFitPanel()
279	tabProc("",0) //force the tab control to display correctly
280	else
281	DoWindow/F AutoFitPanel
282	endif
283	// RangeCheckProc("",1) //check the box to use all of the data
284	End
285
286	Function InitializeAutoFit()
287	//set up all of the globals for the lists, globals, etc...
288	NewDataFolder/O/S root:AutoFit
289	Make/O/T/N=(1,2) fileListWave //list boxes need text waves
290	Make/O/B/N=(1,2) fileSelWave //for the checkboxes
291	Make/O/T/N=1 guessList,holdList,ConstrList,rangeList //list boxes need text waves
292	Make/O/T/N=(1,1) guessMatrix,holdMatrix,rangeMatrix //guesses, etc... need real values (later)
293	Make/O/T/N=(1,1) constrMatrix //constraint matrix is text
294	Make/O/B/N=1 guessSel,holdSel,ConstrSel,rangeSel //SelWave is byte data
295	Variable/G numPar=0,ptLow=0,ptHigh=0,fitTol=1e-3,curProgress=0,endProgress=10
296	Variable/G startTicks=0
297
298	String/G gStatFormat="\K(65535,0,0)\f01"
299	String/G gStatus=gStatFormat+"the fit engine is currently idle"
300	string/g gExt="ABC"
301
302	String/G guessStr = "Guess is based on\rthe unsmeared\rcoefficients on the\rCurve Fit Setup\rpanel"
303	String/G holdStr = "Enter 0 to fit\rEnter 1 to hold"
304	String/G constStr = "Enter low and\rhigh range as:\rlow;high\rfor each parameter\rLeave blank if\rno constraints"
305	String/G rangeStr = "Enter the point\rrange as:\rfirst pt in 1st row\rlast pt in 2nd row\rLeave zero to\rfit all data"
306
307
308	fileListWave = ""
309
310	//make the second column of the file sel wave a checkbox (set bit 5)
311	fileSelWave[][1] = 2^5
312
313	SetDataFolder root:
314	End
315
316	Window AutoFitPanel()
317	PauseUpdate; Silent 1 // building window...
318	NewPanel /W=(434,44,1063,422) /K=1
319	DoWindow/C AutoFitPanel
320	SetDrawLayer UserBack
321	/// DrawText 11,418,"Point Range"
322	/// DrawText 5,167,"Suffix"
323	// DrawLine 8,362,287,362
324	PopupMenu popup0,pos={13,6},size={175,20},title="Pick a Function"
325	PopupMenu popup0,mode=1,value= #"User_FunctionPopupList()",proc=AF_FuncPopupProc
326
327	ListBox lb pos={11,90},size={280,230},proc=FileListboxProc
328	ListBox lb listWave=root:AutoFit:fileListWave,selWave=root:AutoFit:fileSelWave
329	ListBox lb editStyle=1,mode=5
330	ListBox lb userColumnResize=1,widths={200,80}
331
332	/// Button DelButton,pos={245,61},size={40,20},proc=DelButtonProc,title="Del",disable=2
333	Button PathButton,pos={13,61},size={50,20},proc=PDPButton,title="Path..."
334	Button FileListButton,pos={82,61},size={50,20},proc=FileListButtonProc,title="List"
335	CheckBox AllFilesCheck,pos={196,64},size={32,14},proc=CheckAllFiles,title="Select All Files"
336	CheckBox AllFilesCheck,value= 0
337	ListBox guessBox,pos={345,102},size={145,208},disable=1,proc=UpdateGuessMatrixProc
338	ListBox guessBox,frame=2,listWave=root:AutoFit:guessList
339	ListBox guessBox,selWave=root:AutoFit:guessSel,mode= 2,selRow= 0
340	Button FillAllGuessButton,pos={517,110},size={50,20},disable=1,proc=FillAllGuessButtonProc,title="Fill All"
341	Button FillAllHoldButton,pos={517,110},size={50,20},disable=1,proc=FillAllHoldButtonProc,title="Fill All"
342	Button FillAllConstrButton,pos={517,110},size={50,20},proc=FillAllConstrButtonProc,title="Fill All"
343	Button FillAllRangeB,pos={517,110},size={50,20},disable=1,proc=FillAllRangeButtonProc,title="Fill All"
344	SetVariable NumParams,pos={13,35},size={150,15},proc=SetNumParamProc,title="Number of Parameters"
345	SetVariable NumParams,limits={2,Inf,0},value= root:AutoFit:numPar
346	/// CheckBox typeCheck,pos={207,31},size={32,14},title="NSE Data?",value=0
347	/// SetVariable fitTol,pos={80,208},size={80,15},title="Fit Tol"
348	/// SetVariable fitTol,limits={0.0001,0.1,0},value= root:AutoFit:fitTol
349	CheckBox epsilonCheck,pos={156,340},size={32,14},value=0,title="Use Epsilon Wave?"
350
351	TabControl tabC,pos={334,75},size={273,244},proc=tabProc,tabLabel(0)="Guess"
352	TabControl tabC,tabLabel(1)="Hold",tabLabel(2)="Constraint",tabLabel(3)="Range",value= 0
353	/// CheckBox rangeCheck,pos={92,404},size={32,14},proc=RangeCheckProc,title="All"
354	/// CheckBox rangeCheck,value= 1,disable=2
355
356	//// SetVariable lowPt,pos={136,404},size={60,15},title="low"
357	//// SetVariable lowPt,limits={0,Inf,0},value= root:AutoFit:ptLow,noedit=1,disable=1
358	//// SetVariable highPt,pos={201,404},size={60,15},title=" to "
359	//// SetVariable highPt,limits={0,Inf,0},value= root:AutoFit:ptHigh,noedit=1,disable=1
360	ListBox holdBox,pos={345,102},size={145,208},disable=1,proc=UpdateHoldMatrixProc
361	ListBox holdBox,frame=2,listWave=root:AutoFit:holdList
362	ListBox holdBox,selWave=root:AutoFit:holdSel,mode= 2,selRow= 2
363	ListBox ConstrBox,pos={345,102},size={145,208},proc=UpdateConstrMatrixProc
364	ListBox ConstrBox,frame=2,listWave=root:AutoFit:ConstrList
365	ListBox ConstrBox,selWave=root:AutoFit:ConstrSel,mode= 2,selRow= 2
366	ListBox RangeBox,pos={345,102},size={145,208},proc=UpdateRangeMatrixProc
367	ListBox RangeBox,frame=2,listWave=root:AutoFit:rangeList
368	ListBox RangeBox,selWave=root:AutoFit:RangeSel,mode= 2,selRow= 2
369	/// Button MatrixButton,pos={12,205},size={60,20},proc=ToMatrixButtonProc,title="Matrix",disable=2
370	Button DoItButton,pos={339,336},size={80,20},proc=DoTheFitsButtonProc,title="Do the fits"
371	/// Button savePathButton,pos={82,61},size={80,20},proc=SavePathButtonProc,title="Save Path..."
372	TitleBox tb1,pos={479,330},size={128,12},anchor=MC,variable=root:AutoFit:gStatus,frame=0
373	Button button0,pos={14,336},size={60,20},title="Plot",proc=LoadForGuessProc
374	/// SetVariable extStr,pos={4,170},size={40,15},title=" ",value= root:AutoFit:gExt
375
376	Button GuessCoefB,pos={519,144},size={50,20},title="Guess",proc=UseCoefAsGuess
377	Button GuessHoldB,pos={519,144},size={50,20},title="Guess",disable=1,proc=UseHoldAsGuess
378	/// Button GuessConstrB,pos={198,440},size={50,20},title="Guess",disable=1,proc=UseConstraintsAsGuess
379
380	TitleBox infoTitleBox pos={501,174},fixedSize=1,size={96,130},variable=root:AutoFit:guessStr
381
382	ValDisplay progressValdisp,pos={446,349},size={161,7},title="00:00"
383	ValDisplay progressValdisp,limits={0,root:AutoFit:endProgress,0},barmisc={0,0},value= root:AutoFit:curProgress
384
385	Button AF_DoneButton,pos={541,8},size={70,20},proc=AF_DoneButtonProc,title="Done"
386	Button AF_HelpButton,pos={496,8},size={30,20},proc=AF_HelpButtonProc,title="?"
387
388	EndMacro
389
390	// close the panel
391	// + other windows?
392	//
393	Proc AF_DoneButtonProc(ctrlName): ButtonControl
394	String ctrlName
395	DoWindow/K AutoFitPanel
396	end
397
398	//open the Help file for the AutoFit Panel
399	Function AF_HelpButtonProc(ba) : ButtonControl
400	STRUCT WMButtonAction &ba
401
402	switch( ba.eventCode )
403	case 2: // mouse up
404	// click code here
405	DisplayHelpTopic/Z/K=1 "Auto Fit (or Batch fitting)"
406	if(V_flag !=0)
407	DoAlert 0,"The Auto Fit Help file could not be found"
408	endif
409	break
410	endswitch
411
412	return 0
413	End
414
415
416	Function AF_FuncPopupProc(pa) : PopupMenuControl
417	STRUCT WMPopupAction &pa
418
419	switch( pa.eventCode )
420	case 2: // mouse up
421	Variable popNum = pa.popNum
422	String popStr = pa.popStr
423	String name=pa.ctrlName
424	Variable num=0,nc
425
426	String coefStr = getFunctionCoef(popStr)
427	String str
428
429	NVAR pts1 = root:AutoFit:numPar
430	SVAR suff = root:AutoFit:gExt
431
432	Wave/Z coef=$("root:"+coefStr)
433	if(WaveExists(coef))
434	num=numpnts(coef)
435	suff = getModelSuffix(popStr) //get the suffix from the unsmeared model function
436	else
437	//try getting rid of "Smeared" from the beginning of the function string
438	nc = strlen(popStr)
439	str = getFunctionCoef(popStr[7,nc-1])
440	Wave/Z cw = $("root:"+str)
441	if(WaveExists(cw))
442	num=numpnts(cw)
443	endif
444	suff = getModelSuffix(popStr[7,nc-1]) //get the suffix from the unsmeared model function
445	endif
446
447	// set the globals
448	pts1 = num
449
450
451
452	// set the size of the matrix of initial guesses
453	SetNumParamProc("",num,"","")
454
455	break
456	endswitch
457
458	return 0
459	End
460
461
462	Function tabProc(name,tab)
463	String name
464	Variable tab
465
466	String infoStr="",guessStr,holdStr,constStr,rangeStr
467	SVAR str0 = root:AutoFit:guessStr
468	SVAR str1 = root:AutoFit:holdStr
469	SVAR str2 = root:AutoFit:constStr
470	SVAR str3 = root:AutoFit:rangeStr
471
472
473
474	Button FillAllGuessButton disable= (tab!=0)
475	Button GuessCoefB disable=(tab!=0)
476	ListBox guessBox disable= (tab!=0)
477
478
479	Button FillAllHoldButton disable= (tab!=1)
480	Button GuessHoldB disable=(tab!=1)
481	ListBox holdBox disable=(tab!=1)
482
483	Button FillAllConstrButton disable= (tab!=2)
484	// Button GuessConstrB disable=(tab!=2)
485	ListBox ConstrBox disable=(tab!=2)
486
487	//no buttons on the range tab
488	ListBox RangeBox disable=(tab!=3)
489	Button FillAllRangeB disable= (tab!=3)
490
491	if(tab==0)
492	TitleBox infoTitleBox variable=str0
493	endif
494	if(tab==1)
495	TitleBox infoTitleBox variable=str1
496	endif
497	if(tab==2)
498	TitleBox infoTitleBox variable=str2
499	endif
500	if(tab==3)
501	TitleBox infoTitleBox variable=str3
502	endif
503
504
505	return(0)
506	End
507
508
509	Function FillAllGuessButtonProc(ctrlName) : ButtonControl
510	String ctrlName
511
512	//fills all of the guesses with the "current" list of guesses
513	WAVE/T gl=root:AutoFit:guessList
514	WAVE/T gm=root:AutoFit:guessMatrix
515	gm[][] = gl[p]
516	return(0)
517	End
518
519	Function FillAllHoldButtonProc(ctrlName) : ButtonControl
520	String ctrlName
521
522	//fills all of the guesses with the "current" list of guesses
523	WAVE/T hl=root:AutoFit:holdList
524	WAVE/T hm=root:AutoFit:holdMatrix
525	hm[][] = hl[p]
526	return(0)
527	End
528
529	Function FillAllConstrButtonProc(ctrlName) : ButtonControl
530	String ctrlName
531
532	//fills all of the guesses with the "current" list of guesses
533	WAVE/T cl=root:AutoFit:constrList
534	WAVE/T cm=root:AutoFit:ConstrMatrix
535	cm[][] = cl[p]
536	return(0)
537	End
538
539	Function FillAllRangeButtonProc(ctrlName) : ButtonControl
540	String ctrlName
541
542	//fills all of the range values with the top values
543	WAVE/T rl=root:AutoFit:rangeList
544	WAVE/T rm=root:AutoFit:rangeMatrix
545	rm[][] = rl[p]
546	return(0)
547	End
548
549	Function SetNumParamProc(ctrlName,varNum,varStr,varName) : SetVariableControl
550	String ctrlName
551	Variable varNum
552	String varStr
553	String varName
554
555	//set the number of parameters for the function
556	Variable/G numPar=varNum
557
558	// adjust the number of parameters in the guesses
559	SetMatrixRows(varNum)
560
561	return(0)
562	End
563
564	//expand or collapse the matrices to the number of parameters
565	Function SetMatrixRows(newRows)
566	Variable newRows
567
568	SetDataFolder root:AutoFit
569
570	WAVE/T aa = guessMatrix
571	WAVE/T bb = holdMatrix
572	WAVE/T cc = constrMatrix
573	WAVE/T rr = rangeMatrix
574
575	Redimension/N=(newRows,-1) aa,bb,cc,rr
576
577	//redimension the lists too
578	WAVE/T dd=guessList
579	WAVE/T ee=holdList
580	WAVE/T ff=constrList
581	WAVE/T rrl=rangeList
582	Redimension/N=(newRows) dd,ee,ff,rrl
583
584	WAVE gg=guessSel
585	WAVE hh=holdSel
586	WAVE ii=constrSel
587	WAVE rrs=rangeSel
588	//set the selection wave for the individual lists to 2 == editable
589	Redimension/N=(newRows) gg,hh,ii,rrs
590	gg=2
591	hh=2
592	ii=2
593	rrs=2
594
595	SetDataFolder root:
596
597	return 0
598	end
599
600
601	Function FileListboxProc(ctrlName,row,col,event)
602	String ctrlName // name of this control
603	Variable row // row if click in interior, -1 if click in title
604	Variable col // column number
605	Variable event // event code
606
607	if(event==4 && col==0) //selection of a file, not a checkbox
608	//then update the listbox display of guesses for the new file
609	// Print "file selected"
610	DisplayGuess()
611	Endif
612
613	return 0 // other return values reserved
614	End
615
616	Function UpdateGuessMatrixProc(ctrlName,row,col,event)
617	String ctrlName // name of this control
618	Variable row // row if click in interior, -1 if click in title
619	Variable col // column number
620	Variable event // event code
621
622	if(event==7) //finished edit of a cell, update the matrix of guesses
623	//Print "updating"
624	WAVE/T gl=root:AutoFit:guessList
625	WAVE/T gm=root:AutoFit:guessMatrix
626	//put new value in guessMatrix
627	Variable datafile
628	datafile = FindSelectedRowInLB(0) //selected row in file list
629	gm[row][datafile] = gl[row]
630	//print "row= ",row
631	Endif
632	return 0 // other return values reserved
633	End
634
635	Function UpdateHoldMatrixProc(ctrlName,row,col,event)
636	String ctrlName // name of this control
637	Variable row // row if click in interior, -1 if click in title
638	Variable col // column number
639	Variable event // event code
640
641	if(event==7) //finished edit of a cell, update the matrix of guesses
642	//Print "updating"
643	WAVE/T hl=root:AutoFit:holdList
644	WAVE/T hm=root:AutoFit:holdMatrix
645	// put new value in holdMatrix
646	Variable datafile
647	datafile = FindSelectedRowInLB(0)
648	hm[row][datafile] = hl[row]
649	//print "row= ",row
650	Endif
651	return 0 // other return values reserved
652	End
653
654	Function UpdateConstrMatrixProc(ctrlName,row,col,event)
655	String ctrlName // name of this control
656	Variable row // row if click in interior, -1 if click in title
657	Variable col // column number
658	Variable event // event code
659
660	if(event==7) //finished edit of a cell, update the matrix of guesses
661	//Print "updating"
662	WAVE/T cl=root:AutoFit:ConstrList
663	WAVE/T cm=root:AutoFit:constrMatrix
664	//put new value in constrMatrix
665	Variable datafile
666	datafile = FindSelectedRowInLB(0)
667	cm[row][datafile] = cl[row]
668	//print "row= ",row
669	Endif
670	return 0 // other return values reserved
671	End
672
673	Function UpdateRangeMatrixProc(ctrlName,row,col,event)
674	String ctrlName // name of this control
675	Variable row // row if click in interior, -1 if click in title
676	Variable col // column number
677	Variable event // event code
678
679	if(event==7) //finished edit of a cell, update the matrix of guesses
680	//Print "updating"
681	WAVE/T cl=root:AutoFit:RangeList
682	WAVE/T cm=root:AutoFit:RangeMatrix
683	//put new value in constrMatrix
684	Variable datafile
685	datafile = FindSelectedRowInLB(0)
686	cm[row][datafile] = cl[row]
687	//print "row= ",row
688	Endif
689	return 0 // other return values reserved
690	End
691
692	//returns the total number of files checked
693	Function numChecked()
694
695	WAVE sel = root:AutoFit:fileSelWave
696
697	Variable ii,num=0
698	for(ii=0;ii<DimSize(sel,0);ii+=1)
699	if( (sel[ii][1] & 0x10) == 16 )
700	num += 1
701	endif
702	endfor
703	return num
704
705	End
706
707
708	//Function RangeCheckProc(ctrlName,checked) : CheckBoxControl
709	// String ctrlName
710	// Variable checked
711	//
712	// SetVariable lowPt disable=(checked)
713	// SetVariable highPt disable=(checked)
714	//End
715
716	Function DoTheFitsButtonProc(ctrlName) : ButtonControl
717	String ctrlName
718
719	//read the values from the panel
720	ControlInfo/W=AutoFitPanel popup0
721	String funcStr=S_Value
722	String fileName
723	NVAR numPar=root:AutoFit:numPar
724	NVAR startTicks=root:AutoFit:startTicks
725	NVAR curProgress=root:AutoFit:curProgress
726	NVAR endProgress=root:AutoFit:endProgress
727
728	WAVE/T fileWave=root:AutoFit:fileListWave
729	WAVE sel = root:AutoFit:fileSelWave
730
731	WAVE/T guessMatrix=root:AutoFit:guessMatrix
732	WAVE/T holdMatrix=root:AutoFit:holdMatrix
733	WAVE/T constrMatrix=root:AutoFit:constrMatrix
734	WAVE/T rangeMatrix=root:AutoFit:rangeMatrix
735
736	Variable ii=0,numFiles=DimSize(fileWave,0),t1=ticks,t2,matrixindex,count=0
737	//print numFiles
738	//for the timer-progress bar - timer shows time, bar shows number of fits remaining
739	curProgress = 0
740	endProgress = 3*numChecked() //3 fits per file that is fit
741	ValDisplay progressValdisp,limits={0,endProgress,0}
742	ControlUpdate/W=AutoFitPanel progressValdisp
743
744	//print "numChecked() = ",numChecked()
745	startTicks=t1 //starting now
746	for(ii=0;ii<numFiles;ii+=1) //loop over everything in the file list, dispatching only checked items
747	if( (sel[ii][1] & 0x10) == 16 ) //the box is checked, then go fit it
748	fileName = fileWave[ii][0]
749	DoOneFit(funcStr,fileName,numPar,fileWave,ii,guessMatrix,holdMatrix,constrMatrix,rangeMatrix)
750	count+=1
751	endif
752	endfor
753	UpdateStatusString("All files processed")
754	t2=ticks
755	t1=(t2-t1)/60
756	Printf "All %d files processed in %g seconds =(%g minutes) =(%g hours)\r",count,t1,t1/60,t1/60/60
757	return(0)
758	End
759
760
761	Function DoOneFit(funcStr,fileName,numPar,fileWave,ii,guessMatrix,holdMatrix,constrMatrix,rangeMatrix)
762	String funcStr,fileName
763	Variable numPar
764	WAVE/T fileWave //length <= ii
765	Variable ii //index in the full list
766	WAVE/T guessMatrix,holdMatrix,constrMatrix,rangeMatrix
767
768	String list="",fullName="",dataname="",str=""
769	//be sure all is clear (no data, window... etc)
770	DoWindow/K AutoGraph //no error is reported
771	DoWindow/K AutoGraph_NSE
772	list=WaveList("__",";","") //not case sensitive?
773	DoKillList(list) //get rid of all of the old data
774	list=WaveList("wave*",";","") //"wave(n)" can be generated if not all NSE times are collected
775	DoKillList(list)
776
777	//SANS or NSE data?
778	ControlInfo/W=AutoFitPanel typeCheck
779	Variable isNSE=v_value
780	//load in the data
781	PathInfo dataPath
782	fullName = S_path
783	//fullName += fileWave[ii]
784	fullName += fileName
785	if(isNSE)
786	//load, don't plot
787	str = "A_LoadNSEDataWithName(\""+fullName+"\",0)"
788	Execute str
789	// A_LoadNSEDataWithName(fullName,0) //load, don't plot
790	else
791	//load, default name, no plot, overwrite data if needed
792	str = "A_LoadOneDDataToName(\""+fullName+"\",\"\",0,1)"
793	Execute str
794	// A_LoadOneDDataToName(fullName,"",0,1) //load, default name, no plot, overwrite data if needed
795	endif
796	dataName=GetFileNameFromPathNoSemi(fullName)
797	dataname=CleanupName(dataname,0)
798
799	//graph it - creates window named "AutoGraph"
800	if(isNSE)
801	Execute "DoGraph_NSE(\""+dataname+"\")"
802	else
803	Execute "DoGraph(\""+dataname+"\")"
804	endif
805
806	//generate the needed strings to add to the execute stmt
807	// rangeStr, guessStr,constrStr,holdStr
808	// and waves - the coefficient wave and paramter wave
809	Variable useRes,useEps,useCursors,useConstr
810	Variable val,pt1,pt2
811	useEps = 0
812	useConstr = 0
813	useCursors = 0
814	useRes = 0
815
816	String holdStr="",rangeStr="",cStr="",epsilonStr=""
817	String DF="root:"+dataname+":"
818
819
820	Make/O/D/N=(numPar) myCoef //always named "myCoef"
821	Wave cw = root:myCoef
822
823	//range string
824	useCursors = ParseRangeString(ii,dataname,rangeMatrix,pt1,pt2) //read range from the matrix, return values and set flag
825
826
827	//params to hold
828	holdStr = ParseHoldString(ii,numpar,holdMatrix)
829	//initial guesses
830	ParseGuessString(ii,numpar,guessMatrix,"myCoef")
831	//constraints
832	cStr = ParseConstraints(ii,numpar,constrMatrix,"myConstraints")
833	WAVE/Z constr = $"myConstraints"
834	if(strlen(cStr)!=0)
835	useConstr = 1
836	endif
837	//epsilon wave
838	epsilonStr = GetEpsilonWave()
839	if(strlen(epsilonStr) > 0)
840	Wave eps=$epsilonStr //in the root folder
841	useEps = 1
842	endif
843
844	//fit it
845	//these two global variables must be created for IGOR to set them, so that
846	//any errors in the fitting can be diagnosed
847	Variable/G V_FitError=0 //0=no err, 1=error,(2^1+2^0)=3=singular matrix
848	Variable/G V_FitQuitReason=0 //0=ok,1=maxiter,2=user stop,3=no chisq decrease
849	NVAR tol=root:AutoFit:fitTol
850	Variable/G V_FitTol=tol //default is 0.0001
851	Variable/G V_fitOptions=4 //suppress the dialog during fitting (we're not waiting for OK, anyways)
852
853
854	Wave yw = $(DF+dataName+"_i")
855	Wave xw = $(DF+dataName+"_q")
856	Wave sw = $(DF+dataName+"_s")
857
858	if(stringmatch(funcStr, "Smear*")) // if it's a smeared function, need a struct
859	useRes=1
860	endif
861
862	Struct ResSmearAAOStruct fs
863	WAVE/Z resW = $(DF+dataName+"_res") //these may not exist, if 3-column data is used
864	WAVE/Z fs.resW = resW
865	WAVE yw=$(DF+dataName+"_i")
866	WAVE xw=$(DF+dataName+"_q")
867	WAVE sw=$(DF+dataName+"_s")
868	Wave fs.coefW = cw
869	Wave fs.yW = yw
870	Wave fs.xW = xw
871
872	Duplicate/O yw $(DF+"FitYw")
873	WAVE fitYw = $(DF+"FitYw")
874	fitYw = NaN
875
876	NVAR useGenCurveFit = root:Packages:NIST:gUseGenCurveFit
877	///// SEE FitWrapper() for more details.
878
879	Variable nPass = 1
880
881	do // outer do is the loop over passes
882	do // this inner loop is abig switch to select the correct FutFunc to dispatch
883	if(useGenCurveFit)
884	#if !(exists("GenCurveFit"))
885	// XOP not available
886	useGenCurveFit = 0
887	Abort "Genetic Optimiztion XOP not available. Reverting to normal optimization."
888	#endif
889	//send everything to a function, to reduce the clutter
890	// useEps and useConstr are not needed
891	// pass the structure to get the current waves, including the trimmed USANS matrix
892	Variable chi,pt
893
894	chi = DoGenCurveFit(useRes,useCursors,sw,fitYw,fs,funcStr,ParseHoldString(ii,numpar,holdMatrix),val,lolim,hilim,pt1,pt2)
895	pt = val
896
897	break
898
899	endif
900
901
902	if(useRes && useEps && useCursors && useConstr) //do it all
903	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw[pt1,pt2] /X=xw /W=sw /I=1 /E=eps /D=fitYw /C=constr /STRC=fs
904	break
905	endif
906
907	if(useRes && useEps && useCursors) //no constr
908	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw[pt1,pt2] /X=xw /W=sw /I=1 /E=eps /D=fitYw /STRC=fs
909	break
910	endif
911
912	if(useRes && useEps && useConstr) //no crsr
913	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /E=eps /D=fitYw /C=constr /STRC=fs
914	break
915	endif
916
917	if(useRes && useCursors && useConstr) //no eps
918	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw[pt1,pt2] /X=xw /W=sw /I=1 /D=fitYw /C=constr /STRC=fs
919	break
920	endif
921
922	if(useRes && useCursors) //no eps, no constr
923	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw[pt1,pt2] /X=xw /W=sw /I=1 /D=fitYw /STRC=fs
924	break
925	endif
926
927	if(useRes && useEps) //no crsr, no constr
928	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /E=eps /D=fitYw /STRC=fs
929	break
930	endif
931
932	if(useRes && useConstr) //no crsr, no eps
933	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /D=fitYw /C=constr /STRC=fs
934	break
935	endif
936
937	if(useRes) //just res
938	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /D=fitYw /STRC=fs
939	break
940	endif
941
942	///// same as above, but all without useRes (no /STRC flag)
943	if(useEps && useCursors && useConstr) //do it all
944	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw[pt1,pt2] /X=xw /W=sw /I=1 /E=eps /D=fitYw /C=constr
945	break
946	endif
947
948	if(useEps && useCursors) //no constr
949	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw[pt1,pt2] /X=xw /W=sw /I=1 /E=eps /D=fitYw
950	break
951	endif
952
953	if(useEps && useConstr) //no crsr
954	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /E=eps /D=fitYw /C=constr
955	break
956	endif
957
958	if(useCursors && useConstr) //no eps
959	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw[pt1,pt2] /X=xw /W=sw /I=1 /D=fitYw /C=constr
960	break
961	endif
962
963	if(useCursors) //no eps, no constr
964	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw[pt1,pt2] /X=xw /W=sw /I=1 /D=fitYw
965	break
966	endif
967
968	if(useEps) //no crsr, no constr
969	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /E=eps /D=fitYw
970	break
971	endif
972
973	if(useConstr) //no crsr, no eps
974	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /D=fitYw /C=constr
975	break
976	endif
977
978	//just a plain vanilla fit
979	FuncFit/H=ParseHoldString(ii,numpar,holdMatrix) /NTHR=0 $funcStr cw, yw /X=xw /W=sw /I=1 /D=fitYw
980
981	while(0) //always exit the inner do to select the FutFunc syntax. The break will exit this loop
982
983	if(nPass == 1)
984	UpdateStatusString(dataname +" Fit #1")
985	Variable/G gChiSq1 = V_chisq
986	endif
987
988	if(nPass == 2)
989	UpdateStatusString(dataname +" Fit #2")
990	Variable/G gChiSq2 = V_chisq
991	endif
992
993	if(nPass == 3)
994	UpdateStatusString(dataname +" Fit #3")
995	endif
996
997	nPass += 1
998	while(nPass < 4)
999
1000	// append the fit
1001	// need to manage duplicate copies
1002	// Don't plot the full curve if cursors were used (set fitYw to NaN on entry...)
1003	String traces=TraceNameList("", ";", 1 ) //"" as first parameter == look on the target graph
1004	if(strsearch(traces,"FitYw",0) == -1)
1005	if(useGenCurveFit && useCursors)
1006	WAVE trimX = trimX
1007	AppendtoGraph fitYw vs trimX
1008	else
1009	AppendToGraph FitYw vs xw
1010	endif
1011	else
1012	RemoveFromGraph FitYw
1013	if(useGenCurveFit && useCursors)
1014	WAVE trimX = trimX
1015	AppendtoGraph fitYw vs trimX
1016	else
1017	AppendToGraph FitYw vs xw
1018	endif
1019	endif
1020	ModifyGraph lsize(FitYw)=2,rgb(FitYw)=(52428,1,1)
1021
1022
1023	//do the report calling the function, not the proc
1024	String topGraph= WinName(0,1) //this is the topmost graph
1025	SVAR suffix = root:AutoFit:gExt //don't try to read as getModelSuffix(funcStr), since the smeared() may not be plotted
1026	String parStr=GetWavesDataFolder(cw,1)+ WaveList("param"+"_"+suffix, "", "TEXT:1," ) // this is hopefully one wave
1027
1028	if(isNSE)
1029	GenerateReport_NSE(funcStr,dataname,"par","myCoef",1)
1030	else
1031	W_GenerateReport(funcStr,dataname,$parStr,$"myCoef",1,V_chisq,W_sigma,V_npnts,V_FitError,V_FitQuitReason,V_startRow,V_endRow,topGraph)
1032
1033	// AutoFit_GenerateReport(funcStr,dataname,"par","myCoef",1)
1034	endif
1035
1036
1037	// if checked, also write out the fitted waves and coefficients as a text file, maybe an Igor Text format, or just plain text.
1038	if(1)
1039	AF_SaveFitWaves(dataName,$"myCoef",W_sigma,fitYw,xw)
1040	endif
1041
1042	return(0)
1043	End
1044
1045	// if checked, also write out the fitted waves and coefficients as a text file, maybe an Igor Text format, or just plain text.
1046	//
1047	Function AF_SaveFitWaves(dataName,cw,cw_err,yw,xw)
1048	String dataName
1049	Wave cw,cw_err,yw,xw
1050
1051	String saveName,formatStr
1052	Variable refNum
1053
1054	PathInfo dataPath
1055	saveName = S_path+dataName+"_fit.txt"
1056	formatStr="%15.4g\t%15.4g"
1057	formatStr += "\r"
1058
1059	// Print saveName
1060
1061	Open refNum as saveName
1062
1063	fprintf refnum,"Model data created %s\r\n",(date()+" "+time())
1064	wfprintf refnum,formatStr,xw,yw
1065
1066	close refnum
1067
1068	// then the coefficients
1069	saveName = S_path+dataName+"_coef.txt"
1070	Open refNum as saveName
1071
1072	fprintf refnum,"Model coefficients and One Std. Deviation %s\r\n",(date()+" "+time())
1073	wfprintf refnum,formatStr,cw,cw_err
1074
1075	close refnum
1076
1077
1078	return(0)
1079	End
1080
1081
1082
1083	//// to be able to simply fill in the pt range of the yw[lo,hi]
1084	//Function LowRange()
1085	// NVAR val = root:AutoFit:ptLow
1086	// return(val)
1087	//End
1088	//
1089	//Function HighRange()
1090	// NVAR val = root:AutoFit:ptHigh
1091	// return(val)
1092	//End
1093
1094	//updates the status string, forcing a redraw twice, to overcome what seems to be a bug in Igor4
1095	//also updates the progress bar and the countdown
1096	Function UpdateStatusString(newStr)
1097	String newStr
1098
1099	SVAR gStatus=root:AutoFit:gStatus
1100	SVAR gStatFormat=root:AutoFit:gStatFormat
1101	gStatus=""
1102	ControlUpdate/W=AutoFitPanel tb1 //clears the line
1103	gStatus=gStatFormat+newStr
1104	ControlUpdate/W=AutoFitPanel tb1
1105
1106	//update the progress bar and timer
1107	NVAR cur=root:AutoFit:curProgress
1108	NVAR endProgress=root:AutoFit:endProgress
1109	NVAR startTicks=root:AutoFit:startTicks
1110	cur+=1
1111
1112	// Print "****"
1113	// Print "cur = ",cur
1114
1115	Variable t2=ticks,projTot,projRemain
1116	String titleStr=""
1117
1118	projTot = (t2-startTicks)/(cur-1)*endProgress
1119	projRemain = projTot - (t2-startTicks)
1120	projRemain /= 60.15 //to seconds
1121	DoWindow/F AutoFitPanel
1122	ValDisplay progressValdisp,title=Secs2time(projRemain,5)
1123	ControlUpdate/W=AutoFitPanel progressValdisp
1124
1125	// Print "End Progress = ",endProgress
1126	// Print "ProjTot = ",projTot/60.15
1127	// Print "ProjRemain = ",projRemain
1128
1129	return(0)
1130	End
1131
1132	// returns a string with the name of the epsilon wave, or null if it doesn't exist
1133	Function/S GetEpsilonWave()
1134
1135	String retStr="",testStr
1136	SVAR extStr = root:AutoFit:gExt
1137	//do you want to use epsilon wave?
1138	ControlInfo/W=AutoFitPanel epsilonCheck
1139	if(V_Value == 0) //not checked
1140	return(retStr)
1141	endif
1142
1143	testStr = "Epsilon_"+extStr
1144	if(waveexists($testStr) != 0)
1145	// return(" /E="+testStr)
1146	return(testStr)
1147	endif
1148	return(retStr)
1149	End
1150
1151	//returns the y data to fit - "_i" suffix, and range if selected
1152	// enter the range as the first element[0] as low point, second element[1] as high point
1153	// leave the other elements blank
1154	//
1155	// if either point value is not entered, the whole data set is used
1156	//
1157	// -- if both of the elements are between 0 and 1, assume they are q-values (check for fractional part)
1158	// and find the appropriate points...
1159	//
1160	Function ParseRangeString(set,dataStr,rMat,lo,hi)
1161	Variable set
1162	String dataStr
1163	WAVE/T rMat
1164	Variable &lo,&hi
1165
1166	Variable useCursors
1167
1168	lo = str2num(rMat[0][set])
1169	hi = str2num(rMat[1][set])
1170
1171	if(lo == 0 && hi == 0)
1172	return(0)
1173	endif
1174
1175	if(numtype(lo) != 0 \|\| numtype(hi) != 0) //not a normal number
1176	return(0) //don't use cursor values, use all of the dataset
1177	endif
1178
1179
1180	if(lo>=hi) //error in specification, or all data desired (both zero or undefined)
1181	return(0) //use all of the dataset
1182	endif
1183
1184
1185	if(trunc(lo) != lo \|\| trunc(hi) != hi) //values have fractional parts
1186	WAVE xW = $("root:"+dataStr+":"+dataStr + "_q")
1187	FindLevel/P/Q xW, lo
1188	lo = trunc(V_levelX) //find the corresponding point values
1189	FindLevel/P/Q xW, hi
1190	hi = trunc(V_levelX)
1191	endif
1192
1193	Variable/G root:AutoFit:ptLow = lo
1194	Variable/G root:AutoFit:ptHigh = hi
1195
1196	return(1)
1197	End
1198
1199	//returns the holdStr, i.e. "1010101"
1200	Function/S ParseHoldString(set,npar,hMat)
1201	Variable set,npar
1202	WAVE/T hMat
1203
1204	Variable ii
1205	String retStr=""
1206
1207	for(ii=0;ii<npar;ii+=1)
1208	if(strlen(hMat[ii][set]) == 0)
1209	retStr += "0"
1210	else
1211	retStr += hMat[ii][set]
1212	endif
1213	endfor
1214
1215	// Print retStr
1216
1217	return(retStr)
1218	End
1219
1220	//returns only error - wave must already exist
1221	Function ParseGuessString(set,npar,gMat,outWStr)
1222	Variable set,npar
1223	WAVE/T gMat
1224	String outWStr
1225
1226	Variable ii
1227	String retStr=outWStr
1228
1229	WAVE outWave=$outWStr
1230	for(ii=0;ii<npar;ii+=1)
1231	outWave[ii] = str2num(gMat[ii][set])
1232	endfor
1233	return(0)
1234	End
1235
1236	//the constraint wave is created here
1237	Function/S ParseConstraints(set,npar,cMat,outWStr)
1238	Variable set,npar
1239	WAVE/T cMat
1240	String outWStr
1241
1242	String retStr="",loStr="",hiStr=""
1243	Variable ii=0,jj=0
1244	Make/O/T/N=0 $outWStr
1245	WAVE/T cw = $outWStr
1246	//return "/C="+outWStr if parsed ok, or null if not
1247
1248	String listStr=""
1249	for(ii=0;ii<npar;ii+=1)
1250	listStr=cMat[ii][set]
1251	loStr=StringFromList(0, listStr ,";")
1252	hiStr=StringFromList(1, listStr ,";")
1253	if(cmpstr(loStr,"")!=0) //something there...
1254	InsertPoints jj,1,cw //insert before element jj
1255	cw[jj] = "K"+num2Str(ii)+">"+loStr
1256	jj+=1
1257	endif
1258	if(cmpstr(hiStr,"")!=0) //something there...
1259	InsertPoints jj,1,cw //insert before element jj
1260	cw[jj] = "K"+num2Str(ii)+"<"+hiStr
1261	jj+=1
1262	endif
1263	endfor
1264	if(numpnts(cw) > 0 )
1265	retStr = "/C=" + outWStr
1266	Endif
1267	// print "strlen, conStr= ",strlen(retStr),retstr
1268
1269	return(retStr) //null string will be returned if no constraints
1270	End
1271
1272
1273	//*************************
1274	Proc DoGraph(dataname)
1275	String dataname
1276
1277	SetDataFolder $dataname
1278
1279	PauseUpdate; Silent 1 // building window...
1280	Display /W=(5,42,301,313)/K=1 $(dataname+"_i") vs $(dataname+"_q")
1281	DoWindow/C AutoGraph
1282	ModifyGraph mode($(dataname+"_i"))=3
1283	ModifyGraph marker($(dataname+"_i"))=19
1284	ModifyGraph rgb($(dataname+"_i"))=(1,4,52428)
1285	ModifyGraph msize($(dataname+"_i"))=2
1286	ModifyGraph grid=1
1287	ModifyGraph log=1
1288	ModifyGraph tickUnit=1
1289	//
1290	ModifyGraph mirror=2
1291	//
1292	ErrorBars $(dataname+"_i") Y,wave=($(dataname+"_s"),$(dataname+"_s"))
1293	// Legend/A=LB
1294	// Legend/A=LB/X=0.00/Y=-32 //defaults to the top right
1295	Legend/A=MT/E // /X=-16.16/Y=5.17/E
1296
1297	SetDataFolder root:
1298
1299	EndMacro
1300
1301	Proc DoGraph_NSE(dataname)
1302	String dataname
1303
1304	SetDataFolder $dataname
1305
1306	PauseUpdate; Silent 1 // building window...
1307	Display /W=(5,42,301,313)/K=1 $(dataname+"_i") vs $(dataname+"_q")
1308	DoWindow/C AutoGraph_NSE
1309	ModifyGraph mode($(dataname+"_i"))=3
1310	ModifyGraph marker($(dataname+"_i"))=19
1311	ModifyGraph rgb($(dataname+"_i"))=(1,4,52428)
1312	ModifyGraph msize($(dataname+"_i"))=3
1313	ModifyGraph zColor($(dataname+"_i"))={$(dataname+"__Q"),$(dataname+"__Q")[0]0.99,,Rainbow}
1314	ModifyGraph grid=1
1315	ModifyGraph mirror=2
1316	ModifyGraph log(left)=1
1317	ErrorBars $(dataname+"_i") Y,wave=($(dataname+"_s"),$(dataname+"_s"))
1318	Legend/A=MB
1319	SetAxis/A/E=1 bottom
1320	ModifyGraph standoff=0
1321	SetAxis left 0.1,1
1322	//ModifyGraph tickUnit(left)=1
1323	EndMacro
1324
1325	//*************************
1326
1327
1328	//**************************
1329	// Generate the report
1330	//**************************
1331
1332	////must have AutoGraph as the name of the graph window (any size)
1333	//// func is the name of the function (for print only)
1334	////par and coef are the exact names of the waves
1335	////yesSave==1 will save the file(name=func+time)
1336	////
1337	////
1338	//// general report function from the wrapper is used instead
1339	//// W_GenerateReport(funcStr,folderStr,$parStr,cw,yesSave,V_chisq,W_sigma,V_npnts,V_FitError,V_FitQuitReason,V_startRow,V_endRow,topGraph)
1340	////
1341	//Function AutoFit_GenerateReport(func,dataname,par,myCoef,yesSave)
1342	// String func,dataname,par,myCoef
1343	// Variable yesSave
1344	//
1345	// String str,pictStr="P_"
1346	// Wave sigWave=$"W_sigma"
1347	// Wave ans=$myCoef
1348	//// Wave/T param=$par
1349	// SVAR ext=root:AutoFit:gExt
1350	// Wave/T param=$"parameters_"+ext
1351	//
1352	// NVAR V_chisq = V_chisq
1353	// NVAR V_npnts = V_npnts
1354	// NVAR V_FitError = V_FitError
1355	// NVAR V_FitQuitReason = V_FitQuitReason
1356	// NVAR chi1=gChiSq1
1357	// NVAR chi2=gChiSq2
1358	// NVAR V_startRow = V_startRow
1359	// NVAR V_endRow = V_endRow
1360	//
1361	// // bring report up
1362	// DoWindow/F Report
1363	// if (V_flag == 0) // Report notebook doesn't exist ?
1364	// NewNotebook/W=(10,45,550,620)/F=1/N=Report as "Report"
1365	// endif
1366	//
1367	// // delete old stuff
1368	// Notebook Report selection={startOfFile, endOfFile}, text="\r", selection={startOfFile, startOfFile}
1369	//
1370	// // insert title
1371	// Notebook Report newRuler=Title, justification=1, rulerDefaults={"Times", 16, 1, (0, 0, 0)}
1372	// sprintf str, "Fit to %s, %s, %s\r\r", func,Secs2Date(datetime, 0), time()
1373	// Notebook Report ruler=Title, text=str
1374	//
1375	// // insert fit results
1376	// Variable num=numpnts(ans),ii=0
1377	// Notebook Report ruler=Normal; Notebook Report margins={18,18,504}, tabs={63 + 3*8192}
1378	// str = "Data file: " + dataname + "\r\r"
1379	// Notebook Report text=str
1380	// do
1381	// sprintf str, "%s = %g±%g\r", param[ii],ans[ii],sigwave[ii]
1382	// Notebook Report text=str
1383	// ii+=1
1384	// while(ii<num)
1385	//
1386	// //
1387	// Wave dataXw = $(dataname+"_q")
1388	// Variable nData=numpnts(dataXw)
1389	//
1390	// //
1391	// sprintf str,"chisq = %g , %g , %g\r",chi1,chi2,V_chisq
1392	// Notebook Report textRGB=(65000,0,0),fstyle=1,text=str
1393	// sprintf str,"Npnts = %g, sqrt(X^2/N) = %g\r",V_npnts,sqrt(V_chisq/V_npnts)
1394	// Notebook Report textRGB=(0,0,0),fstyle=0, text=str
1395	// sprintf str "Fitted range = [%d,%d] = %g < Q < %g\r",V_startRow,V_endRow,dataXw(V_startRow),dataXw(V_endRow)
1396	// Notebook Report textRGB=(0,0,0),fstyle=0, text=str
1397	// sprintf str,"V_FitError = %g\t\tV_FitQuitReason = %g\r",V_FitError,V_FitQuitReason
1398	// Notebook Report textRGB=(65000,0,0),fstyle=1,text=str
1399	// Notebook Report ruler=Normal
1400	//
1401	// // insert graphs
1402	// Notebook Report picture={AutoGraph(0, 0, 400, 300), 0, 1}, text="\r"
1403	//
1404	//
1405	// //Notebook Report picture={Table1, 0, 0}, text="\r"
1406	//
1407	// // show the top of the report
1408	// Notebook Report selection= {startOfFile, startOfFile}, findText={"", 1}
1409	//
1410	// //save the notebook
1411	// if(yesSave)
1412	// String nameStr=CleanupName(func,0)
1413	// nameStr = nameStr[0,8] //shorten the name
1414	// nameStr += "_"+dataname
1415	// //make sure the name is no more than 31 characters
1416	// namestr = namestr[0,30] //if shorter than 31, this will NOT pad to 31 characters
1417	// Print "file saved as ",nameStr
1418	// SaveNotebook /O/P=savePath/S=2 Report as nameStr
1419	// //save the graph separately as a PICT file, 2x screen
1420	// pictStr += nameStr
1421	// pictStr = pictStr[0,28] //need a shorter name - why?
1422	// DoWindow/F AutoGraph
1423	// // E=-5 is png @screen resolution
1424	// // E=2 is PICT @2x screen resolution
1425	// SavePICT /E=2/O/I/P=savePath /W=(0,0,3,3) as pictStr
1426	// Endif
1427	//
1428	// // ???maybe print the notebook too?
1429	//End
1430
1431
1432
1433	//
1434	//Proc DoGenerateReport(fname,dataname,param,coef,yesno)
1435	// String fname,dataname="",param,coef,yesno="No"
1436	// Prompt fname,"function name",popup,FunctionList("!Sme*",";","KIND:15,")
1437	// Prompt dataname,"name of the dataset"
1438	// Prompt param,"parameter names",popup,WaveList("par*",";","")
1439	// Prompt coef,"coefficient wave",popup,WaveList("coe*",";","")
1440	// Prompt yesno,"Save the report?",popup,"Yes;No;"
1441	//
1442	// //Print "fname = ",fname
1443	// //Print "param = ",param
1444	// //Print "coef = ",coef
1445	// //Print "yesno = ",yesno
1446	//
1447	// Variable doSave=0
1448	// if(cmpstr(yesno,"Yes")==0)
1449	// doSave=1
1450	// Endif
1451	// GenerateReport(fname,dataname,param,coef,dosave)
1452	//
1453	//End
1454
1455
1456
1457	//must have AutoGraph_NSE as the name of the graph window (any size)
1458	// func is the name of the function (for print only)
1459	//par and coef are the exact names of the waves
1460	//yesSave==1 will save the file(name=func+time)
1461	//
1462	Function GenerateReport_NSE(func,dataname,par,coef,yesSave)
1463	String func,dataname,par,coef
1464	Variable yesSave
1465
1466	String str
1467	Wave sigWave=$"W_sigma"
1468	Wave ans=$coef
1469	Wave/T param=$par
1470
1471	NVAR V_chisq = V_chisq
1472	NVAR V_npnts = V_npnts
1473	NVAR V_FitError = V_FitError
1474	NVAR V_FitQuitReason = V_FitQuitReason
1475
1476	// bring report up
1477	DoWindow/F Report
1478	if (V_flag == 0) // Report notebook doesn't exist ?
1479	NewNotebook/W=(10,45,550,620)/F=1/N=Report as "Report"
1480	endif
1481
1482	// delete old stuff
1483	Notebook Report selection={startOfFile, endOfFile}, text="\r", selection={startOfFile, startOfFile}
1484
1485	// insert title
1486	Notebook Report newRuler=Title, justification=1, rulerDefaults={"Times", 16, 1, (0, 0, 0)}
1487	sprintf str, "Fit to %s, %s, %s\r\r", func,Secs2Date(datetime, 0), time()
1488	Notebook Report ruler=Title, text=str
1489
1490	// insert fit results
1491	Variable num=numpnts(ans),ii=0
1492	Notebook Report ruler=Normal; Notebook Report margins={18,18,504}, tabs={63 + 3*8192}
1493	str = "Data file: " + dataname + "\r\r"
1494	Notebook Report text=str
1495	do
1496	sprintf str, "%s = %g±%g\r", param[ii],ans[ii],sigwave[ii]
1497	Notebook Report text=str
1498	ii+=1
1499	while(ii<num)
1500
1501	sprintf str,"chisq = %g\r",V_chisq
1502	Notebook Report textRGB=(65000,0,0),fstyle=1,text=str
1503	sprintf str,"Npnts = %g\r",V_npnts
1504	Notebook Report textRGB=(0,0,0),fstyle=0, text=str
1505	sprintf str,"V_FitError = %g\rV_FitQuitReason = %g\r",V_FitError,V_FitQuitReason
1506	Notebook Report textRGB=(65000,0,0),fstyle=1,text=str
1507	Notebook Report ruler=Normal
1508
1509	// insert graphs
1510	Notebook Report picture={AutoGraph_NSE(0, 0, 400, 300), 0, 1}, text="\r"
1511	//Notebook Report picture={Table1, 0, 0}, text="\r"
1512
1513	// show the top of the report
1514	Notebook Report selection= {startOfFile, startOfFile}, findText={"", 1}
1515
1516	//save the notebook
1517	if(yesSave)
1518	String nameStr=CleanupName(func,0)
1519	nameStr = nameStr[0,8] //shorten the name
1520	nameStr += "_"+dataname
1521	//make sure the name is no more than 31 characters
1522	namestr = namestr[0,30] //if shorter than 31, this will NOT pad to 31 characters
1523	Print "file saved as ",nameStr
1524	// SaveNotebook /O/P=savePath/S=2 Report as nameStr
1525	SaveNotebook /O/P=home/S=2 Report as nameStr
1526	Endif
1527
1528	//save a pict file, just of the graph
1529	namestr = "PNG_" + namestr
1530	namestr = namestr[0,28] //PICT names in IGOR must be shorter, to allow auto-naming?
1531	DoWindow/F AutoGraph_NSE
1532	// SavePICT /E=2/O/I/P=savePath /W=(0,0,3,3) as nameStr
1533	SavePICT /E=2/O/I/P=home /W=(0,0,3,3) as nameStr
1534	End
1535
1536
1537	//*********************
1538	//
1539	//*********************
1540	// List utilities
1541	//*********************
1542	//
1543	// List2TextWave is also in SANS_Utilities.ipf, but to keep from having big lists of dependencies,
1544	// have definitions here that are local and Static
1545
1546	//
1547	Static Function List2TextWave(list,sep,waveStr)
1548	String list,sep,waveStr
1549
1550	Variable n= ItemsInList(list,sep)
1551	Make/O/T/N=(n) $waveStr= StringFromList(p,list,sep)
1552	End
1553
1554	Function List2NumWave(list,sep,waveStr)
1555	String list,sep,waveStr
1556
1557	Variable n= ItemsInList(list,sep)
1558	Make/O/D/N=(n) $waveStr= str2num( StringFromList(p,list,sep) )
1559	End
1560
1561	Function/S TextWave2List(w,sep)
1562	Wave/T w
1563	String sep
1564
1565	String newList=""
1566	Variable n=numpnts(w),ii=0
1567	do
1568	newList += w[ii] + sep
1569	ii+=1
1570	while(ii<n)
1571	return(newList)
1572	End
1573
1574	//for single precision numerical waves
1575	Function/S NumWave2List(w,sep)
1576	Wave w
1577	String sep
1578
1579	String newList="",temp=""
1580	Variable n=numpnts(w),ii=0,val
1581	do
1582	val=w[ii]
1583	temp=""
1584	sprintf temp,"%g",val
1585	newList += temp
1586	newList += sep
1587	ii+=1
1588	while(ii<n)
1589	return(newList)
1590	End
1591
1592	Function DoKillList(list)
1593	String list
1594
1595	String item=""
1596	do
1597	item = StringFromList(0, list ,";" )
1598	KillWaves/Z $item
1599	list = RemoveFromList(item, list, ";")
1600	while(ItemsInList(list, ";")>0)
1601	End
1602
1603	//*********************
1604
1605
1606
1607	//*******************
1608	//Notebook Utils
1609	//*******************
1610
1611	Proc PrintOpenNotebooks()
1612
1613	String list=WinList("*", ";", "WIN:16")
1614	String item = ""
1615	do
1616	item=StringFromList(0,list,";")
1617	Print item
1618	PrintNotebook $item
1619	list = RemoveFromList(item, list, ";")
1620	while(ItemsInList(list, ";")>0)
1621	End
1622
1623	Proc CloseOpenNotebooks()
1624
1625	String list=WinList("*", ";", "WIN:16")
1626	String item = ""
1627	do
1628	item=StringFromList(0,list,";")
1629	DoWindow/K $item
1630	list = RemoveFromList(item, list, ";")
1631	while(ItemsInList(list, ";")>0)
1632	End
1633
1634	Proc Generate_PNG_list()
1635	fGenerate_PNG_list()
1636	End
1637
1638	//spits up a list of PNGS
1639	Function fGenerate_PNG_list()
1640
1641	// String List=IndexedFile(savePath, -1, ".png")
1642	String List=IndexedFile(home, -1, ".png")
1643	List2TextWave(List,";","PNG_files")
1644	WAVE/T PNGs=$"PNG_files"
1645	Edit/K=1 PNGs
1646	End
1647
1648	Proc Layout_PNGS(pStr)
1649	String pStr=""
1650	Prompt pStr,"wave of PNGs to use",popup,WaveList("PNG*", ";", "")
1651
1652	String List=TextWave2List($pStr,";")
1653	String item = ""
1654	//kill old picts from memory
1655	KillPicts/A/Z
1656	//make a new layout
1657	Layout/C=1 as "PNG_Layout"
1658	DoWindow/C PNGLayout
1659	do
1660	item=StringFromList(0,List,";")
1661	//load each PNG, and append it to the layout
1662	Print "load item = ",item
1663	// LoadPICT /O/Q/P=savePath item
1664	LoadPICT /O/Q/P=home item
1665	DoWindow/F PNGLayout //be sure layout is on top
1666	AppendLayoutObject /F=1/W=PNGLayout picture $CleanupName(item,0)
1667	//AppendToLayout $item
1668	//Print item
1669	List = RemoveFromList(item, List, ";")
1670	while(ItemsInList(List, ";")>0)
1671	//tile the PNGs in the layout
1672	Tile/O=8
1673	End
1674
1675	//spit up a list of everything in the folder that is a text file
1676	Proc Generate_Data_Checklist()
1677	fGenerate_Data_list()
1678	End
1679	Function fGenerate_Data_list()
1680
1681	String List=IndexedFile(dataPath, -1, "TEXT")
1682	List2TextWave(List,";","data_files")
1683	WAVE/T files=$"data_files"
1684	Duplicate/O/T files ModelToUse,Done
1685	ModelToUse = ""
1686	Done = ""
1687	Edit/K=1 files,ModelToUse,Done
1688	End
1689
1690	//***************
1691
1692	// with all of the notebooks open for a given model (with N parameters)
1693	// this will comile the reports, save the text notebook, then load it back in
1694	// and put it in a table for easy printing
1695	//
1696	Proc CompileAndLoadReports()
1697	Compile_Reports()
1698	// SaveNotebook /O/P=savePath/S=2 Compilation
1699	SaveNotebook /O/P=home/S=2 Compilation
1700	// Print "Saved as: ",S_Path
1701	DoWindow/K Compilation
1702	LoadCompiledReports(S_Path)
1703	End
1704
1705	Proc LoadCompiledReports(str)
1706	String str
1707	LoadWave/J/D/W/E=1/K=0/V={"\t"," $",0,0} str
1708	End
1709
1710	Proc Compile_Reports(Number_of_Parameters,tagStr)
1711	Variable Number_of_Parameters=(root:AutoFit:numPar)
1712	String tagStr="v1"
1713
1714	fCompile_ReportsTable(Number_of_Parameters,tagStr)
1715
1716	fCompile_Reports(Number_of_Parameters)
1717
1718	DoWindow/F AutoFitResults
1719
1720	End
1721
1722	Proc Compile_GlobalFit_Reports(WhichParameter,FileStr,ParamStr,ParamErrStr)
1723	Variable WhichParameter=7
1724	String FileStr="FileName_",ParamStr="par_",ParamErrStr="parErr_"
1725	fCompile_GlobalFit_Reports(WhichParameter,FileStr,ParamStr,ParamErrStr)
1726	End
1727
1728	////
1729	Function fCompile_Reports(nPar)
1730	Variable nPar
1731
1732	String list=WinList("*", ";", "WIN:16")
1733	String item = "",textStr="",newStr=""
1734	Variable sstop,dum,chi
1735
1736	Variable ii
1737	NewNotebook/F=0 /N=Compilation as "Compilation"
1738	Notebook Compilation text="Variable Name \tValues\tErrors\r"
1739	do
1740	item=StringFromList(0,list,";")
1741	DoWindow/F $item
1742
1743	Notebook $item selection={(2,0), (3,0)} //paragraph 3 (starts from 0) = filename
1744	GetSelection notebook,$item,2
1745	textStr=S_Selection
1746	textStr=textStr[0,strlen(textStr)-2] //remove CR
1747	textStr=textStr[10,strlen(textStr)-1] //remove "DATA FILE:
1748	Notebook Compilation text="File:\t"+textStr+"\t"
1749	Notebook Compilation text=textStr+"_err\r"
1750
1751	// printf "%s %s\r",textStr,textStr+"_err"
1752	//
1753	// results are written as: "%s = \t%g\t±\t%g\r"
1754	//
1755	for(ii=0;ii<nPar;ii+=1) //gather the parameters
1756
1757	Notebook $item selection={(ii+4,0), (ii+5,0)} //paragraph 5 = parameter 0
1758	GetSelection notebook,$item,2
1759	textStr=S_Selection
1760	textStr=textStr[0,strlen(textStr)-2]
1761	//textStr=textStr[13,strlen(textStr)-1] // remove "parameter"
1762	textStr = ReplaceString("=",textStr,"")
1763	textStr = ReplaceString("±\t",textStr,"")
1764	Notebook Compilation text=textStr+"\r"
1765	//printf "%s\r",textStr
1766	//Print strlen(textStr)
1767	endfor
1768	//get the chi-squared/N value
1769	ii+=1
1770	Notebook $item selection={(ii+4,0), (ii+5,0)} //
1771	GetSelection notebook,$item,2
1772	textStr=S_Selection
1773
1774	// Print textStr
1775
1776	sscanf textStr,"Npnts = %g\t\tSqrt(X^2/N) = %g\r",dum,chi
1777	sprintf textStr,"Sqrt(X^2/N)\t%s\t\r",num2str(chi)
1778	Notebook Compilation text=textStr
1779
1780	Notebook Compilation text="\r" //separate items
1781
1782	list = RemoveFromList(item, list, ";")
1783	while(ItemsInList(list, ";")>0)
1784
1785	DoWindow/F Compilation
1786	Notebook Compilation selection={startOfFile,startOfFile}
1787	End
1788
1789	////
1790	// compiles the results into a table.
1791	// add tagStr to the end to (try) to make sure that the compilations are unique
1792	Function fCompile_ReportsTable(nPar,tagStr)
1793	Variable nPar
1794	String tagStr
1795
1796	String list=WinList("*", ";", "WIN:16")
1797	String item = "",textStr="",newStr="",str
1798	Variable sstop,dum,chi
1799
1800	Variable ii,numRep,jj,val1,val2,pt
1801	numRep=ItemsInList(list,";")
1802
1803	Make/O/T/N=(numRep) $("fittedFiles"+"_"+tagStr)
1804	Make/O/D/N=(numRep) $("chiSQ"+"_"+tagStr)
1805	Wave/T fittedFiles = $("fittedFiles"+"_"+tagStr)
1806	Wave chiSq = $("chiSQ"+"_"+tagStr)
1807	Edit/N=AutoFitResults fittedFiles,chiSQ
1808
1809	for(ii=0;ii<nPar;ii+=1) //waves for the parameters
1810	Make/O/D/N=(numRep) $("par"+num2str(ii)+"_"+tagStr),$("par"+num2str(ii)+"_err"+"_"+tagStr)
1811	AppendToTable $("par"+num2str(ii)+"_"+tagStr),$("par"+num2str(ii)+"_err"+"_"+tagStr)
1812	endfor
1813
1814	for(jj=0;jj<numRep;jj+=1)
1815
1816	item=StringFromList(jj,list,";")
1817	DoWindow/F $item
1818	Notebook $item selection={(2,0), (3,0)} //paragraph 3 (starts from 0) = filename
1819	GetSelection notebook,$item,2
1820	textStr=S_Selection
1821	textStr=textStr[0,strlen(textStr)-2] //remove CR
1822	textStr=textStr[10,strlen(textStr)-1] //remove "DATA FILE:
1823
1824	fittedFiles[jj] = textStr
1825	//
1826	// results are written as: "%s = \t%g\t±\t%g\r"
1827	//
1828	for(ii=0;ii<nPar;ii+=1) //gather the parameters
1829	Notebook $item selection={(ii+4,0), (ii+5,0)} //paragraph 5 = parameter 0
1830	GetSelection notebook,$item,2
1831	textStr=S_Selection
1832	textStr=textStr[0,strlen(textStr)-2]
1833	// find the "="
1834	pt = strsearch(textStr,"=",0,0)
1835	textStr = textStr[pt,strlen(textStr)-1]
1836	sscanf textStr,"= \t%g\t ± \t%g\r",val1,val2
1837
1838	// Print textStr
1839
1840	Wave w1 = $("par"+num2str(ii)+"_"+tagStr)
1841	Wave w2 = $("par"+num2str(ii)+"_err"+"_"+tagStr)
1842
1843	w1[jj] = val1
1844	w2[jj] = val2
1845
1846	//textStr = ReplaceString("=",textStr,"")
1847	//textStr = ReplaceString("±\t",textStr,"")
1848	//Notebook Compilation text=textStr+"\r"
1849	//printf "%s\r",textStr
1850	//Print strlen(textStr)
1851	endfor
1852	//get the chi-squared/N value
1853	ii += 1
1854	Notebook $item selection={(ii+4,0), (ii+5,0)} //
1855	GetSelection notebook,$item,2
1856	textStr=S_Selection
1857
1858	// Print textStr
1859
1860	sscanf textStr,"Npnts = %g\t\tSqrt(X^2/N) = %g\r",dum,chi
1861
1862	chiSQ[jj] = chi
1863	endfor
1864
1865
1866	return(0)
1867	End
1868
1869
1870	////
1871	Function fCompile_GlobalFit_Reports(nPar,Files,Param,Param_err)
1872	Variable nPar
1873	String Files,Param,Param_Err
1874
1875	String list=SortList(WinList("*", ";", "WIN:16"))
1876	String item = "",textStr="",fileStr,funcStr,dumStr
1877	Variable sstop,dum,chi,nFiles,val1,val2
1878
1879	Variable ii,jj
1880
1881	nFiles = ItemsInList(list)
1882	Make/O/D/N=(nFiles) $Param,$Param_err
1883	Make/O/T/N=(nFiles) $Files
1884	WAVE pw = $Param
1885	WAVE pw_err = $Param_err
1886	WAVE/T fw = $Files
1887
1888	for(jj=0;jj<nFiles;jj+=1)
1889
1890	item=StringFromList(jj,list,";")
1891	DoWindow/F $item
1892
1893	Notebook $item selection={(6,0), (7,0)} //paragraph 7 (starts from 0) = filename
1894	GetSelection notebook,$item,2
1895	textStr=S_Selection
1896	sscanf textStr,"Data Set:%s ; Function:%s\r",fileStr,funcStr
1897	fw[jj] = fileStr
1898	// Notebook Compilation text=fileStr+"\r"
1899	// Notebook Compilation text=funcStr+"\r"
1900
1901
1902	// printf "%s %s\r",textStr,textStr+"_err"
1903
1904	for(ii=nPar;ii<=nPar;ii+=1) //gather the parameters
1905
1906	Notebook $item selection={(ii+7,0), (ii+8,0)} //paragraph 8 = parameter 0
1907	GetSelection notebook,$item,2
1908	textStr=S_Selection
1909	// Print textStr
1910	sscanf textStr,"%s %g +- %g\r",dumStr,val1,val2
1911	pw[jj] = val1
1912	pw_err[jj] = val2
1913	endfor
1914
1915	endfor
1916
1917	End
1918
1919
1920
1921	//DOES NOT graph the data, does not create a weighting wave
1922	//NSE data assumed to be 4-column
1923	//q (all the same) - t(ns) - I(q,t) - dI(q,t)
1924	Static Function LoadNSEDataWithName(fileStr)
1925	String fileStr
1926	//Load the waves, using default waveX names
1927	//if no path or file is specified for LoadWave, the default Mac open dialog will appear
1928	LoadWave/G/D/A/Q fileStr
1929	String fileName = S_fileName
1930
1931
1932	String w0,w1,w2,n0,n1,n2,wt,w3,n3
1933	Variable rr,gg,bb
1934
1935	// put the names of the four loaded waves into local names
1936	n0 = StringFromList(0, S_waveNames ,";" )
1937	n1 = StringFromList(1, S_waveNames ,";" )
1938	n2 = StringFromList(2, S_waveNames ,";" )
1939	n3 = StringFromList(3, S_waveNames ,";" )
1940
1941
1942	//remove the semicolon AND period from files from the VAX
1943	w0 = CleanupName((S_fileName+"__Q"),0)
1944	w1 = CleanupName((S_fileName+"__X"),0)
1945	w2 = CleanupName((S_fileName+"__Y"),0)
1946	w3 = CleanupName((S_fileName+"__S"),0)
1947
1948	if(exists(w0) !=0)
1949	//waves already exist
1950	KillWaves $n0,$n1,$n2,$n3 // kill the default waveX that were loaded
1951	return(0)
1952	endif
1953
1954	// Rename to give nice names
1955	Rename $n0, $w0
1956	Rename $n1, $w1
1957	Rename $n2, $w2
1958	Rename $n3, $w3
1959
1960	End
1961
1962	/////
1963	//wStr is the coef_ wave with the initial guess
1964	Function FillCoefGuess(wStr)
1965	String wStr
1966
1967	Wave w=$wStr
1968	if (! WaveExists(w) )
1969	DoAlert 0, "The coefficient wave "+wStr+" does not exist. Either create the wave or enter the correct suffix"
1970	return(0)
1971	endif
1972	Wave/T guess=root:AutoFit:guessMatrix
1973	Variable ii=0,num=numpnts(w)
1974	Variable tmp
1975
1976	Variable selRow,matrixRow
1977	selRow = FindSelectedRowInLB(0) // send 0 to look at the column with filenames
1978
1979	do
1980	tmp=w[ii]
1981	guess[ii][selRow] = num2str(tmp)
1982	ii+=1
1983	while(ii<num)
1984	End
1985
1986	/////
1987	//wStr is the hold_ wave with the initial guess
1988	Function FillHoldGuess(wStr)
1989	String wStr
1990
1991	Wave w=$wStr
1992	if (! WaveExists(w) )
1993	DoAlert 0, "The hold wave "+wStr+" does not exist. Either create the wave or enter the correct suffix"
1994	return(0)
1995	endif
1996	Wave/T hold=root:AutoFit:holdMatrix
1997	Variable ii=0,num=numpnts(w)
1998	Variable tmp
1999
2000	Variable selRow,matrixRow
2001	selRow = FindSelectedRowInLB(0) // send 0 to look at the column with filenames
2002
2003	do
2004	tmp=w[ii]
2005	hold[ii][selRow] = num2str(tmp)
2006	ii+=1
2007	while(ii<num)
2008	End
2009
2010	/////
2011	//wStr is the Constr_ wave with the initial guess
2012	// --- this is a Text wave
2013	Function FillConstrGuess(wStr)
2014	String wStr
2015
2016	Wave/T w=$wStr
2017	if (! WaveExists(w) )
2018	DoAlert 0, "The Constraint wave "+wStr+" does not exist. Either create the (TEXT) wave or enter the correct suffix"
2019	return(0)
2020	endif
2021	Wave/T constr=root:AutoFit:constrMatrix
2022	Variable ii=0,num=numpnts(w)
2023	String tmp
2024
2025	Variable selRow,matrixRow
2026	selRow = FindSelectedRowInLB(0) // send 0 to look at the column with filenames
2027
2028	do
2029	tmp=w[ii]
2030	constr[ii][selRow] = tmp
2031	ii+=1
2032	while(ii<num)
2033	End
2034
2035	////
2036	//need the extension of the model wave to append
2037	//
2038	// it will plot the UN-smeared model for generating the initial guesses
2039	// -- fitting will use Smeared or unsmeared, depending on the popup selection
2040	//
2041	//
2042	Function fLoadSelected(extStr)
2043	string extStr
2044
2045	Variable datafile,ii
2046	String fullName,dataName,str
2047
2048	WAVE/T fileWave=root:AutoFit:fileListWave
2049	Wave sel = root:AutoFit:fileSelWave
2050	ii=0
2051	do
2052	if((sel[ii][0] & 0x01) == 1) // the file is selected in the list box
2053	datafile = ii
2054	break
2055	endif
2056	ii+=1
2057	while(ii<numpnts(sel))
2058
2059	//load in the data
2060	PathInfo dataPath
2061	fullName = S_path
2062	fullName += fileWave[datafile][0]
2063
2064	str = "A_LoadOneDDataToName(\""+fullName+"\",\"\",0,1)" //don't plot, force overwrite
2065	Execute str
2066
2067	dataName=GetFileNameFromPathNoSemi(fullName)
2068	dataname=CleanupName(dataname,0)
2069
2070	//graph it - creates window named "AutoGraph"
2071	DoWindow/K tmpGraph
2072	Execute "DoGraph(\""+dataname+"\")"
2073	DoWindow/C tmpGraph
2074	// Wave/Z yw = $("ywave_"+extStr)
2075	// Wave/Z xw = $("xwave_"+extStr)
2076
2077	if(exists("ywave_"+extStr) && exists("xwave_"+extStr))
2078	AppendtoGraph $("ywave_"+extStr) vs $("xwave_"+extStr)
2079	else
2080	DoAlert 0,"Model data with the suffix \""+extStr+"\" does not exist. Either plot the model or enter the correct suffix"
2081	endif
2082
2083	// ControlInfo/W=AutoFitPanel popup0
2084	// if(cmpstr(S_Value[0,4], "Smear" ) == 0 )
2085	// SetDataFolder $dataname
2086	// if(exists("smeared_"+extStr) && exists("smeared_qvals"))
2087	//
2088	// AppendtoGraph $("smeared_"+extStr) vs $("smeared_qvals")
2089	// else
2090	// SetDataFolder root:
2091	// DoAlert 0,"Smeared Model data with the suffix \""+extStr+"\" does not exist. Either plot the model or enter the correct suffix"
2092	// endif
2093	//
2094	// else
2095	// if(exists("ywave_"+extStr) && exists("xwave_"+extStr))
2096	// AppendtoGraph $("ywave_"+extStr) vs $("xwave_"+extStr)
2097	// else
2098	// DoAlert 0,"Model data with the suffix \""+extStr+"\" does not exist. Either plot the model or enter the correct suffix"
2099	// endif
2100	// endif
2101
2102
2103	SetDataFolder root:
2104	return(0)
2105	End
2106
2107	//loads the selected file in a temporary graph
2108	Function LoadForGuessProc(ctrlName) : ButtonControl
2109	String ctrlName
2110
2111	SVAR ext=root:AutoFit:gExt
2112	fLoadSelected(ext)
2113	return(0)
2114	End
2115
2116	//makes the holdWave and Constraint wave that are needed to set the values in the list boxes (and the matrix)
2117	// these are not generated when plotting a model.
2118	Function Make_HoldConstraintEps_Waves()
2119
2120	SVAR ext=root:AutoFit:gExt
2121	Variable num=numpnts($("coef_"+ext))
2122	if(WaveExists($("coef_"+ext)) == 0 )
2123	DoAlert 0,"coef_"+ext+" does not exist. You must set the prefix correctly and plot the model first"
2124	else
2125	Make/O/D/N=(num) $("Hold_"+ext)
2126	Make/O/T/N=(num) $("Constr_"+ext)
2127	Make/O/D/N=(num) $("Epsilon_"+ext)=0.0003
2128
2129	String str=WinList("*", ";","WIN:2"),sel
2130	Prompt sel, "Select Table to Append Waves", popup str
2131	DoPrompt "asdf",sel
2132	if (V_Flag)
2133	return 0 // user canceled
2134	endif
2135	DoWindow/F $sel
2136	AppendToTable $("Hold_"+ext), $("Constr_"+ext), $("Epsilon_"+ext)
2137	endif
2138	End
2139	//needs the name (as a string) of the coefficient wave
2140	//it assumes "coef_" plus an extension
2141	//
2142	Function UseCoefAsGuess(ctrlName) : ButtonControl
2143	String ctrlName
2144
2145	String wStr = "coef_"
2146	SVAR ext=root:AutoFit:gExt
2147	FillCoefGuess(wStr+ext)
2148	DisplayGuess()
2149	return(0)
2150	End
2151
2152	//needs the name (as a string) of the hold string wave
2153	//it assumes "hold_" plus an extension
2154	//
2155	Function UseHoldAsGuess(ctrlName) : ButtonControl
2156	String ctrlName
2157
2158	String wStr = "hold_"
2159	SVAR ext=root:AutoFit:gExt
2160	FillHoldGuess(wStr+ext)
2161	DisplayGuess()
2162	return(0)
2163	End
2164
2165	//needs the name (as a string) of the coefficient wave
2166	//it assumes "constr" plus an extension
2167	//
2168	Function UseConstraintsAsGuess(ctrlName) : ButtonControl
2169	String ctrlName
2170
2171	// String wStr = "Constr_"
2172	// SVAR ext=root:AutoFit:gExt
2173	// FillConstrGuess(wStr+ext)
2174	// DisplayGuess()
2175	return(0)
2176	End
2177

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

Download in other formats:

Original Format