Context Navigation

U_CALC.ipf @ 1004

Last change on this file since 1004 was 939, checked in by srkline, 9 years ago
correcting typos in U_CALC and SumSANSModels
File size: 58.9 KB

Line
1	#pragma rtGlobals=1 // Use modern global access method.
2	#pragma IgorVersion=6.1
3
4
5	// USANS version of SASCALC
6
7	// to simulate the intensity from a USANS experiment for planning
8	// see John's instrument paper:
9	//
10	// J. Appl. Cryst. (2005) 38 1004-1011.
11	//
12	// SRK JUL 2009
13
14
15	// ideas:
16	//
17	// - more presets
18	//
19	// - NEED -
20	// - printable output that makes sense to an instrument scientist at least
21	// for instrument setup
22	//
23	// - make sure that the # points <-> increment relation meshes with ICP
24
25	//
26	// X plot as countrate, not absolute scale
27	// X 3e-5 cutoff
28	// X ? don't plot lowest angle range (but needs to be in the count time)
29	// - need empty beam and empty cell count rate vs. aperture (Cd vs. Gd?)
30	//
31
32
33	// need to add in empty and background corrections to see "reduced" data
34	// or at least compare to what the empty cell would give in the same count time
35	//
36	// -? model the direct beam?? currently the "red" region from -1 to 0.6 is almost entirely
37	// the primary beam, so it's a bit artificial (qmin is really ~ 3e-5)
38	//
39
40	//
41	// Need T_wide, T_rock, I peak for proper absolute scaling, but I don't know if it's really important
42	// to be able to simlulate to this extent
43	//
44
45
46	//#include "MultScatter_MonteCarlo_2D"
47
48
49
50	// Bring the UCALC panel to the front
51	// ALWAYS initializes the folders and variables
52	// then draws the panel if necessary
53	Proc Show_UCALC()
54
55	DoWindow/F UCALC
56	if(V_Flag==0)
57	Init_UCALC()
58	UCALC_Panel()
59	CalcTotalCountTime()
60	ControlUpdate/W=UCALC U_popup0 //force a pop of the function list
61	Endif
62
63	End
64
65
66	//set up the global values for the angles, # points, and count time
67	Proc Init_UCALC()
68
69	//
70	NewDataFolder/O root:Simulation // required to calculate the RandomDeviate
71	NewDataFolder/O root:Packages:NIST:SAS // required to calculate the RandomDeviate
72
73	//
74	NewDataFolder/O root:Packages:NIST:USANS:SIM //for the fake raw data
75	NewDataFolder/O/S root:Packages:NIST:USANS:Globals:U_Sim //for constants, panel, etc
76
77	Variable/G gAngLow1 = -1
78	Variable/G gAngHigh1 = 0.6
79	Variable/G gNumPts1 = 33
80	Variable/G gCtTime1 = 20
81	Variable/G gIncr1 = 0.05
82
83	Variable/G gAngLow2 = 0.7
84	Variable/G gAngHigh2 = 1.9
85	Variable/G gNumPts2 = 13
86	Variable/G gCtTime2 = 50
87	Variable/G gIncr2 = 0.1
88
89	Variable/G gAngLow3 = 2
90	Variable/G gAngHigh3 = 4.8
91	Variable/G gNumPts3 = 15
92	Variable/G gCtTime3 = 100
93	Variable/G gIncr3 = 0.2
94
95	Variable/G gAngLow4 = 5
96	Variable/G gAngHigh4 = 9.5
97	Variable/G gNumPts4 = 10
98	Variable/G gCtTime4 = 200
99	Variable/G gIncr4 = 0.5
100
101	Variable/G gAngLow5 = 10
102	Variable/G gAngHigh5 = 19
103	Variable/G gNumPts5 = 10
104	Variable/G gCtTime5 = 500
105	Variable/G gIncr5 = 1
106
107	Variable/G gAngLow6 = 20
108	Variable/G gAngHigh6 = 48
109	Variable/G gNumPts6 = 15
110	Variable/G gCtTime6 = 1000
111	Variable/G gIncr6 = 2
112
113	Variable/G gAngLow7 = 50
114	Variable/G gAngHigh7 = 95
115	Variable/G gNumPts7 = 10
116	Variable/G gCtTime7 = 0
117	Variable/G gIncr7 = 5
118
119	// results, setup values
120	String/G gFuncStr=""
121	String/G gTotTimeStr=""
122	Variable/G gAnalyzerOmega = 7.1e-7 //solid angle of the analyzer, in steradians
123	// Variable/G gBeamCurrent=25000 // (?? not used) beam current Ed*I (n/s) for 5/8" diam = 25000 n/s
124	Variable/G gThick=0.1 //sample thickness (cm)
125	Variable/G gSamTrans=0.8
126	Variable/G g_1D_DoABS = 0 //=1 for abs scale, 0 for just counts
127	Variable/G g_1D_PlotCR = 1 //=1 to plot countrate
128	Variable/G g_1D_AddNoise = 1 // add in appropriate noise to simulation
129
130	// a box for the results
131	// ?? maybe not useful to report
132	Variable/G g_1DEstDetCR = 0 // estimated detector count rate
133	Variable/G g_1DTotCts = 0
134	Variable/G g_1DFracScatt= 0 // ??
135	Variable/G g_1DEstTrans = 0 // ? can I calculate this?
136	Variable/G g_MultScattFraction = 0
137
138
139	// not on panel yet
140	Variable/G g_Empirical_EMP = 0 // use an emperical model for empty cell subtraction
141	Variable/G g_EmptyLevel = 0.7
142	Variable/G g_BkgLevel = 0.6
143
144	SetDataFolder root:
145
146	End
147
148	//// make the display panel a graph with a control bar just as in SASCALC
149	//// so that the subwindow syntax doesn't break all of the other functionality
150	////
151	//Window UCALC_Panel() : Graph
152	// PauseUpdate; Silent 1 // building window...
153	// Display /W=(55,44,670,850) /K=1
154	// ModifyGraph cbRGB=(36929,50412,31845)
155	// DoWindow/C UCALC
156	// DoWindow/T UCALC,"USANS Simulation"
157	// ControlBar 320
158	//
159	// GroupBox group0,pos={5,0},size={577,159},title="Instrument Setup"
160	// GroupBox group1,pos={5,165},size={240,147},title="Sample Setup"
161	// GroupBox group2,pos={327,165},size={259,147},title="Results"
162	//
163	// PopupMenu popup0,pos={17,18},size={165,20},title="Sample Aperture Diam (in)"
164	// PopupMenu popup0,mode=3,popvalue="0.625",value="0.25;0.50;0.625;0.75;1.0;1.75;2.0;"
165	// PopupMenu popup2,pos={220,18},size={165,20},title="Presets"
166	// PopupMenu popup2,mode=3,popvalue="Long Count",value="Short Count;Medium Count;Long Count;"
167	// PopupMenu popup2,proc=UCALC_PresetPopup
168	//
169	// SetDataFolder root:Packages:NIST:USANS:Globals:U_Sim
170	//
171	// Variable top=44,pt=0,inc=18
172	// SetVariable setvar1a,pos={12,top},size={100,15},title="theta min",value= gAngLow1
173	// SetVariable setvar1b,pos={119,top},size={100,15},title="theta max",value= gAngHigh1
174	// SetVariable setvar1c,pos={227,top},size={100,15},title="increm",value= gIncr1
175	// SetVariable setvar1d,pos={335,top},size={100,15},title="# points",value= gNumPts1
176	// SetVariable setvar1e,pos={443,top},size={100,15},title="count (s)",value= gCtTime1
177	// SetVariable setvar1a,labelBack=(65535,32768,32768)
178	//
179	// pt += inc
180	// SetVariable setvar2a,pos={12,top+pt},size={100,15},title="theta min",value= gAngLow2
181	// SetVariable setvar2b,pos={119,top+pt},size={100,15},title="theta max",value= gAngHigh2
182	// SetVariable setvar2c,pos={227,top+pt},size={100,15},title="increm",value= gIncr2
183	// SetVariable setvar2d,pos={335,top+pt},size={100,15},title="# points",value= gNumPts2
184	// SetVariable setvar2e,pos={443,top+pt},size={100,15},title="count (s)",value= gCtTime2
185	// SetVariable setvar2a labelBack=(65535,65533,32768)
186	//
187	// pt += inc
188	// SetVariable setvar3a,pos={12,top+pt},size={100,15},title="theta min",value= gAngLow3
189	// SetVariable setvar3b,pos={119,top+pt},size={100,15},title="theta max",value= gAngHigh3
190	// SetVariable setvar3c,pos={227,top+pt},size={100,15},title="increm",value= gIncr3
191	// SetVariable setvar3d,pos={335,top+pt},size={100,15},title="# points",value= gNumPts3
192	// SetVariable setvar3e,pos={443,top+pt},size={100,15},title="count (s)",value= gCtTime3
193	// SetVariable setvar3a labelBack=(32769,65535,32768)
194	//
195	// pt += inc
196	// SetVariable setvar4a,pos={12,top+pt},size={100,15},title="theta min",value= gAngLow4
197	// SetVariable setvar4b,pos={119,top+pt},size={100,15},title="theta max",value= gAngHigh4
198	// SetVariable setvar4c,pos={227,top+pt},size={100,15},title="increm",value= gIncr4
199	// SetVariable setvar4d,pos={335,top+pt},size={100,15},title="# points",value= gNumPts4
200	// SetVariable setvar4e,pos={443,top+pt},size={100,15},title="count (s)",value= gCtTime4
201	// SetVariable setvar4a labelBack=(32768,65535,65535)
202	//
203	// pt += inc
204	// SetVariable setvar5a,pos={12,top+pt},size={100,15},title="theta min",value= gAngLow5
205	// SetVariable setvar5b,pos={119,top+pt},size={100,15},title="theta max",value= gAngHigh5
206	// SetVariable setvar5c,pos={227,top+pt},size={100,15},title="increm",value= gIncr5
207	// SetVariable setvar5d,pos={335,top+pt},size={100,15},title="# points",value= gNumPts5
208	// SetVariable setvar5e,pos={443,top+pt},size={100,15},title="count (s)",value= gCtTime5
209	// SetVariable setvar5a labelBack=(32768,54615,65535)
210	//
211	// pt += inc
212	// SetVariable setvar6a,pos={12,top+pt},size={100,15},title="theta min",value= gAngLow6
213	// SetVariable setvar6b,pos={119,top+pt},size={100,15},title="theta max",value= gAngHigh6
214	// SetVariable setvar6c,pos={227,top+pt},size={100,15},title="increm",value= gIncr6
215	// SetVariable setvar6d,pos={335,top+pt},size={100,15},title="# points",value= gNumPts6
216	// SetVariable setvar6e,pos={443,top+pt},size={100,15},title="count (s)",value= gCtTime6
217	// SetVariable setvar6a labelBack=(44253,29492,58982)
218	//
219	//// the action procedures and limits/increments
220	// SetVariable setvar1a proc=ThetaMinSetVarProc //,limits={-2,0,0.1}
221	// SetVariable setvar2a proc=ThetaMinSetVarProc
222	// SetVariable setvar3a proc=ThetaMinSetVarProc
223	// SetVariable setvar4a proc=ThetaMinSetVarProc
224	// SetVariable setvar5a proc=ThetaMinSetVarProc
225	// SetVariable setvar6a proc=ThetaMinSetVarProc
226	//
227	////
228	// SetVariable setvar1b proc=ThetaMaxSetVarProc //,limits={0.4,1,0.1}
229	// SetVariable setvar2b proc=ThetaMaxSetVarProc
230	// SetVariable setvar3b proc=ThetaMaxSetVarProc
231	// SetVariable setvar4b proc=ThetaMaxSetVarProc
232	// SetVariable setvar5b proc=ThetaMaxSetVarProc
233	// SetVariable setvar6b proc=ThetaMaxSetVarProc
234	////
235	// SetVariable setvar1c proc=IncrSetVarProc,limits={0.01,0.1,0.01}
236	// SetVariable setvar2c proc=IncrSetVarProc,limits={0.02,0.2,0.02}
237	// SetVariable setvar3c proc=IncrSetVarProc,limits={0.05,0.4,0.05}
238	// SetVariable setvar4c proc=IncrSetVarProc,limits={0.1,1,0.1}
239	// SetVariable setvar5c proc=IncrSetVarProc,limits={0.5,5,1}
240	// SetVariable setvar6c proc=IncrSetVarProc,limits={1,10,2}
241	////
242	// SetVariable setvar1d proc=NumPtsSetVarProc,limits={2,50,1}
243	// SetVariable setvar2d proc=NumPtsSetVarProc,limits={2,50,1}
244	// SetVariable setvar3d proc=NumPtsSetVarProc,limits={2,50,1}
245	// SetVariable setvar4d proc=NumPtsSetVarProc,limits={2,50,1}
246	// SetVariable setvar5d proc=NumPtsSetVarProc,limits={2,50,1}
247	// SetVariable setvar6d proc=NumPtsSetVarProc,limits={2,50,1}
248	////
249	// SetVariable setvar1e proc=CtTimeSetVarProc,limits={-1,50000,1}
250	// SetVariable setvar2e proc=CtTimeSetVarProc,limits={-1,50000,10}
251	// SetVariable setvar3e proc=CtTimeSetVarProc,limits={-1,50000,10}
252	// SetVariable setvar4e proc=CtTimeSetVarProc,limits={-1,50000,30}
253	// SetVariable setvar5e proc=CtTimeSetVarProc,limits={-1,50000,100}
254	// SetVariable setvar6e proc=CtTimeSetVarProc,limits={-1,50000,100}
255	//
256	// Button button0,pos={255,180},size={60,20},fColor=(65535,65535,0),proc=U_SimPlotButtonProc,title="Plot"
257	// Button button1,pos={260,286},size={50,20},proc=U_SaveButtonProc,title="Save"
258	//
259	////checkbox for "easy" mode
260	// CheckBox check0 title="Simple mode?",pos={400,19},proc=EnterModeCheckProc,value=1
261	// ThetaEditMode(2) //checked on startup
262	//
263	//// instrument setup
264	// SetVariable U_setvar0_1,pos={20,211},size={160,15},title="Thickness (cm)"
265	// SetVariable U_setvar0_1,limits={0,inf,0.1},value= root:Packages:NIST:USANS:Globals:U_Sim:gThick
266	// SetVariable U_setvar0_3,pos={20,235},size={160,15},title="Sample Transmission"
267	// SetVariable U_setvar0_3,limits={0,1,0.01},value= root:Packages:NIST:USANS:Globals:U_Sim:gSamTrans
268	// PopupMenu U_popup0,pos={20,185},size={165,20},proc=Sim_USANS_ModelPopMenuProc,title="Model"
269	// PopupMenu U_popup0,mode=1,value= #"U_FunctionPopupList()"
270	// SetVariable setvar0,pos={20,259},size={120,15},title="Empty Level"
271	// SetVariable setvar0,limits={0,10,0.01},value= root:Packages:NIST:USANS:Globals:U_Sim:g_EmptyLevel
272	// SetVariable setvar0_1,pos={20,284},size={120,15},title="Bkg Level"
273	// SetVariable setvar0_1,limits={0,10,0.01},value= root:Packages:NIST:USANS:Globals:U_Sim:g_BkgLevel
274	//
275	// CheckBox check0_4 title="Show EMP?",pos={160,260},proc=ShowEMPCheckProc,value=0
276	//
277	// CheckBox check0_2,pos={253,239},size={60,14},title="CountRate?",variable= root:Packages:NIST:USANS:Globals:U_Sim:g_1D_PlotCR
278	// CheckBox check0_3,pos={262,264},size={60,14},title="Noise?",variable= root:Packages:NIST:USANS:Globals:U_Sim:g_1D_AddNoise
279	//
280	//// a box for the results
281	// SetVariable totalTime,pos={338,185},size={150,15},title="Count time (h:m)",value= gTotTimeStr
282	//// ValDisplay valdisp0,pos={338,210},size={220,13},title="Total detector counts"
283	//// ValDisplay valdisp0,limits={0,0,0},barmisc={0,1000},value= root:Packages:NIST:USANS:Globals:U_Sim:g_1DTotCts
284	// ValDisplay valdisp0_2,pos={338,234},size={220,13},title="Fraction of beam scattered"
285	// ValDisplay valdisp0_2,limits={0,0,0},barmisc={0,1000},value= root:Packages:NIST:USANS:Globals:U_Sim:g_1DFracScatt
286	// ValDisplay valdisp0_3,pos={338,259},size={220,13},title="Estimated transmission"
287	// ValDisplay valdisp0_3,limits={0,0,0},barmisc={0,1000},value=root:Packages:NIST:USANS:Globals:U_Sim:g_1DEstTrans
288	//
289	//
290	// SetDataFolder root:
291	//
292	//EndMacro
293
294	// ??make the USANS simulation results into a graph and a separate panel.
295	// the control bar at the top makes the whole thing too large, and control bars are limited
296	// to 500 pix wide, which is really tight
297	// left of 1055 is good for Mac, 955 is better for Win
298	//
299	Window UCALC_Panel()
300	PauseUpdate; Silent 1 // building window...
301	String platform=UpperStr(IgorInfo(2))
302	Variable pos=strsearch(platform,"WINDOWS",0)
303	if(pos >= 0) //windows
304	Display /W=(55,44,855,450) /K=1
305	else //mac
306	Display /W=(55,44,1055,544) /K=1
307	endif
308
309	ModifyGraph cbRGB=(36929,50412,31845)
310	DoWindow/C UCALC
311	DoWindow/T UCALC,"USANS Simulation"
312	ControlBar/L 500
313
314	GroupBox group0,pos={5,1},size={493,177},title="Instrument Setup"
315	GroupBox group1,pos={5,183},size={240,110},title="Sample Setup"
316	GroupBox group2,pos={5,310},size={240,170},title="Results"
317
318	PopupMenu popup0,pos={17,19},size={165,20},title="Sample Aperture Diam (in)"
319	PopupMenu popup0,mode=3,popvalue="0.625",value="0.25;0.50;0.625;0.75;1.0;1.75;2.0;"
320	PopupMenu popup0,proc=UCALC_SampleAperturePopup
321	PopupMenu popup2,pos={220,19},size={165,20},title="Presets"
322	PopupMenu popup2,mode=2,popvalue="Medium Count",value="Short Count;Medium Count;Long Count;"
323	PopupMenu popup2,proc=UCALC_PresetPopup
324
325	SetDataFolder root:Packages:NIST:USANS:Globals:U_Sim
326
327	Variable top=46,pt=0,inc=18,left=0//left=533
328	SetVariable setvar1a,pos={left+17,top},size={90,15},title="theta min",value= gAngLow1
329	SetVariable setvar1b,pos={left+113,top},size={89,15},title="theta max",value= gAngHigh1
330	SetVariable setvar1c,pos={left+209,top},size={89,15},title="increm",value= gIncr1
331	SetVariable setvar1d,pos={left+299,top},size={100,15},title="# points",value= gNumPts1
332	SetVariable setvar1e,pos={left+399,top},size={93,15},title="count (s)",value= gCtTime1
333	// SetVariable setvar1a,labelBack=(65535,32768,32768) //old rainbow
334	SetVariable setvar1a,labelBack=(49858,65535,65535)
335
336	pt += inc
337	SetVariable setvar2a,pos={left+17,top+pt},size={90,15},title="theta min",value= gAngLow2
338	SetVariable setvar2b,pos={left+113,top+pt},size={89,15},title="theta max",value= gAngHigh2
339	SetVariable setvar2c,pos={left+209,top+pt},size={89,15},title="increm",value= gIncr2
340	SetVariable setvar2d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts2
341	SetVariable setvar2e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime2
342	// SetVariable setvar2a labelBack=(65535,65533,32768) //old rainbow
343	SetVariable setvar2a labelBack=(21074,8995,21074)
344
345	pt += inc
346	SetVariable setvar3a,pos={left+17,top+pt},size={90,15},title="theta min",value= gAngLow3
347	SetVariable setvar3b,pos={left+113,top+pt},size={89,15},title="theta max",value= gAngHigh3
348	SetVariable setvar3c,pos={left+209,top+pt},size={89,15},title="increm",value= gIncr3
349	SetVariable setvar3d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts3
350	SetVariable setvar3e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime3
351	// SetVariable setvar3a labelBack=(32769,65535,32768) //old rainbow
352	SetVariable setvar3a labelBack=(0,60652,60652)
353
354	pt += inc
355	SetVariable setvar4a,pos={left+17,top+pt},size={90,15},title="theta min",value= gAngLow4
356	SetVariable setvar4b,pos={left+113,top+pt},size={89,15},title="theta max",value= gAngHigh4
357	SetVariable setvar4c,pos={left+209,top+pt},size={89,15},title="increm",value= gIncr4
358	SetVariable setvar4d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts4
359	SetVariable setvar4e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime4
360	// SetVariable setvar4a labelBack=(32768,65535,65535) //old rainbow
361	SetVariable setvar4a labelBack=(0,51400,0)
362
363	pt += inc
364	SetVariable setvar5a,pos={left+17,top+pt},size={90,15},title="theta min",value= gAngLow5
365	SetVariable setvar5b,pos={left+113,top+pt},size={89,15},title="theta max",value= gAngHigh5
366	SetVariable setvar5c,pos={left+209,top+pt},size={89,15},title="increm",value= gIncr5
367	SetVariable setvar5d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts5
368	SetVariable setvar5e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime5
369	// SetVariable setvar5a labelBack=(32768,54615,65535) //old rainbow
370	SetVariable setvar5a labelBack=(59367,49344,0)
371
372	pt += inc
373	SetVariable setvar6a,pos={left+17,top+pt},size={90,15},title="theta min",value= gAngLow6
374	SetVariable setvar6b,pos={left+113,top+pt},size={89,15},title="theta max",value= gAngHigh6
375	SetVariable setvar6c,pos={left+209,top+pt},size={89,15},title="increm",value= gIncr6
376	SetVariable setvar6d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts6
377	SetVariable setvar6e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime6
378	// SetVariable setvar6a labelBack=(44253,29492,58982) //old rainbow
379	SetVariable setvar6a labelBack=(54998,0,0)
380
381	pt += inc
382	SetVariable setvar7a,pos={left+17,top+pt},size={90,15},title="theta min",value= gAngLow7
383	SetVariable setvar7b,pos={left+113,top+pt},size={89,15},title="theta max",value= gAngHigh7
384	SetVariable setvar7c,pos={left+209,top+pt},size={89,15},title="increm",value= gIncr7
385	SetVariable setvar7d,pos={left+299,top+pt},size={100,15},title="# points",value= gNumPts7
386	SetVariable setvar7e,pos={left+399,top+pt},size={93,15},title="count (s)",value= gCtTime7
387	// SetVariable setvar7a labelBack=(44253,29492,58982) //old rainbow
388	SetVariable setvar7a labelBack=(39321,21845,51657)
389
390	// the action procedures and limits/increments
391	SetVariable setvar1a proc=ThetaMinSetVarProc //,limits={-2,0,0.1}
392	SetVariable setvar2a proc=ThetaMinSetVarProc
393	SetVariable setvar3a proc=ThetaMinSetVarProc
394	SetVariable setvar4a proc=ThetaMinSetVarProc
395	SetVariable setvar5a proc=ThetaMinSetVarProc
396	SetVariable setvar6a proc=ThetaMinSetVarProc
397	SetVariable setvar7a proc=ThetaMinSetVarProc
398
399	//
400	SetVariable setvar1b proc=ThetaMaxSetVarProc //,limits={0.4,1,0.1}
401	SetVariable setvar2b proc=ThetaMaxSetVarProc
402	SetVariable setvar3b proc=ThetaMaxSetVarProc
403	SetVariable setvar4b proc=ThetaMaxSetVarProc
404	SetVariable setvar5b proc=ThetaMaxSetVarProc
405	SetVariable setvar6b proc=ThetaMaxSetVarProc
406	SetVariable setvar7b proc=ThetaMaxSetVarProc
407	//
408	SetVariable setvar1c proc=IncrSetVarProc,limits={0.01,0.1,0.01}
409	SetVariable setvar2c proc=IncrSetVarProc,limits={0.02,0.2,0.02}
410	SetVariable setvar3c proc=IncrSetVarProc,limits={0.05,0.4,0.05}
411	SetVariable setvar4c proc=IncrSetVarProc,limits={0.1,1,0.1}
412	SetVariable setvar5c proc=IncrSetVarProc,limits={0.5,5,1}
413	SetVariable setvar6c proc=IncrSetVarProc,limits={1,10,2}
414	SetVariable setvar7c proc=IncrSetVarProc,limits={1,20,2}
415	//
416	SetVariable setvar1d proc=NumPtsSetVarProc,limits={2,50,1}
417	SetVariable setvar2d proc=NumPtsSetVarProc,limits={2,50,1}
418	SetVariable setvar3d proc=NumPtsSetVarProc,limits={2,50,1}
419	SetVariable setvar4d proc=NumPtsSetVarProc,limits={2,50,1}
420	SetVariable setvar5d proc=NumPtsSetVarProc,limits={2,50,1}
421	SetVariable setvar6d proc=NumPtsSetVarProc,limits={2,50,1}
422	SetVariable setvar7d proc=NumPtsSetVarProc,limits={2,50,1}
423	//
424	SetVariable setvar1e proc=CtTimeSetVarProc,limits={-1,50000,1}
425	SetVariable setvar2e proc=CtTimeSetVarProc,limits={-1,50000,10}
426	SetVariable setvar3e proc=CtTimeSetVarProc,limits={-1,50000,10}
427	SetVariable setvar4e proc=CtTimeSetVarProc,limits={-1,50000,30}
428	SetVariable setvar5e proc=CtTimeSetVarProc,limits={-1,50000,100}
429	SetVariable setvar6e proc=CtTimeSetVarProc,limits={-1,50000,100}
430	SetVariable setvar7e proc=CtTimeSetVarProc,limits={-1,50000,100}
431
432	Button button0,pos={left+280,200},size={130,20},fColor=(65535,65535,0),proc=U_SimPlotButtonProc,title="Simulate USANS"
433	CheckBox check0_2,pos={left+280,250},size={60,14},title="CountRate?",variable= root:Packages:NIST:USANS:Globals:U_Sim:g_1D_PlotCR
434	CheckBox check0_3,pos={left+280,270},size={60,14},title="Noise?",variable= root:Packages:NIST:USANS:Globals:U_Sim:g_1D_AddNoise
435	CheckBox check0_4 title="Show EMP?",pos={left+280,290},proc=ShowEMPCheckProc,value=0
436
437
438	//checkbox for "easy" mode
439	CheckBox check0 title="Simple mode?",pos={left+400,19},proc=EnterModeCheckProc,value=1
440	ThetaEditMode(2) //checked on startup
441
442	// sample setup
443	SetVariable U_setvar0_1,pos={left+20,231},size={160,15},title="Thickness (cm)"
444	SetVariable U_setvar0_1,limits={0,inf,0.1},value= root:Packages:NIST:USANS:Globals:U_Sim:gThick
445	SetVariable U_setvar0_3,pos={left+20,255},size={160,15},title="Sample Transmission"
446	SetVariable U_setvar0_3,limits={0,1,0.01},value= root:Packages:NIST:USANS:Globals:U_Sim:gSamTrans
447	PopupMenu U_popup0,pos={left+20,205},size={165,20},proc=Sim_USANS_ModelPopMenuProc,title="Model"
448	PopupMenu U_popup0,mode=1,value= #"U_FunctionPopupList()"
449
450	//unused as of yet, not sure I'll ever need these
451	// SetVariable setvar0,pos={left+20,279},size={120,15},title="Empty Level",disable=2
452	// SetVariable setvar0,limits={0,10,0.01},value= root:Packages:NIST:USANS:Globals:U_Sim:g_EmptyLevel
453	// SetVariable setvar0_1,pos={left+20,304},size={120,15},title="Bkg Level",disable=2
454	// SetVariable setvar0_1,limits={0,10,0.01},value= root:Packages:NIST:USANS:Globals:U_Sim:g_BkgLevel
455
456
457
458	// a box for the results
459	SetVariable totalTime,pos={left+20,335},size={150,15},title="Count time (h:m)",value= gTotTimeStr
460	SetVariable totalTime,limits={-inf,inf,0},noedit=1
461	SetVariable valdisp0_2,pos={left+20,365},size={200,15},title="Multiple Coherent Scattering"
462	SetVariable valdisp0_2,limits={-inf,inf,0},format="%8f",noedit=1,value=g_MultScattFraction
463	// ValDisplay valdisp0_3,pos={left+20,420},size={220,13},title="Estimated transmission"
464	// ValDisplay valdisp0_3,limits={0,0,0},barmisc={0,1000},value=root:Packages:NIST:USANS:Globals:U_Sim:g_1DEstTrans
465	// ValDisplay valdisp0_3,disable=1
466	Button button1,pos={left+20,390},size={150,20},proc=U_SavePanelProc,title="Save PNG"
467	Button button2,pos={left+20,420},size={150,20},proc=U_ConfigTextProc,title="Config Text"
468	Button button3,pos={left+20,450},size={150,20},proc=U_SaveButtonProc,title="Save Simulated Data"
469
470	// help, done buttons
471	Button U_helpButton,pos={300,440},size={25,20},proc=showUCALCHelp,title="?"
472	Button U_helpButton,help={"Show help file for simulation of USANS Data"}
473	Button U_DoneButton,pos={350,440},size={50,20},proc=UCALCDoneButton,title="Done"
474	Button U_DoneButton,help={"This button will close the panel"}
475
476	SetDataFolder root:
477
478	EndMacro
479
480	Proc UCALCDoneButton(ctrlName): ButtonControl
481	String ctrlName
482	DoWindow/K UCALC
483	end
484
485	Proc showUCALCHelp(ctrlName): ButtonControl
486	String ctrlName
487	DisplayHelpTopic/K=1/Z "UCALC"
488	if(V_flag !=0)
489	DoAlert 0,"The USANS Simulation Help file could not be found"
490	endif
491	end
492
493	// changing theta min - hold incr and #, result is new theta max
494	Function ThetaMinSetVarProc(sva) : SetVariableControl
495	STRUCT WMSetVariableAction &sva
496
497	switch( sva.eventCode )
498	case 1: // mouse up
499	case 2: // Enter key
500	case 3: // Live update
501	Variable ThetaMin = sva.dval
502	String ns=CtrlNumber(sva.ctrlName) //control number as a string
503	NVAR incr = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gIncr"+ns)
504	NVAR NumPts = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gNumPts"+ns)
505	NVAR thetaMax = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gAngHigh"+ns)
506
507	thetaMax = ThetaMin + incr*NumPts
508
509	break
510	endswitch
511
512	return 0
513	End
514
515
516	// changing theta max - hold min and incr, result is new # of points
517	// then need to recalculate the total counting time
518	Function ThetaMaxSetVarProc(sva) : SetVariableControl
519	STRUCT WMSetVariableAction &sva
520
521	switch( sva.eventCode )
522	case 1: // mouse up
523	case 2: // Enter key
524	case 3: // Live update
525	Variable ThetaMax = sva.dval
526	String ns=CtrlNumber(sva.ctrlName) //control number as a string
527	NVAR thetaMin = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gAngLow"+ns)
528	NVAR NumPts = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gNumPts"+ns)
529	NVAR Incr = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gIncr"+ns)
530
531	NumPts = trunc( (thetaMax - ThetaMin) / incr ) + 1 //+1 to count both ends of the interval
532
533	CalcTotalCountTime()
534
535	break
536	endswitch
537
538	return 0
539	End
540
541	// changing increment - hold min and #, result is new max
542	Function IncrSetVarProc(sva) : SetVariableControl
543	STRUCT WMSetVariableAction &sva
544
545	switch( sva.eventCode )
546	case 1: // mouse up
547	case 2: // Enter key
548	case 3: // Live update
549	Variable incr = sva.dval
550	String ns=CtrlNumber(sva.ctrlName) //control number as a string
551	NVAR thetaMin = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gAngLow"+ns)
552	NVAR NumPts = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gNumPts"+ns)
553	NVAR thetaMax = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gAngHigh"+ns)
554
555	thetaMax = ThetaMin + incr*NumPts
556
557	break
558	endswitch
559
560	return 0
561	End
562
563	// changing #pts - hold min and incr, result is new max
564	// then recalculate the total count time
565	Function NumPtsSetVarProc(sva) : SetVariableControl
566	STRUCT WMSetVariableAction &sva
567
568	switch( sva.eventCode )
569	case 1: // mouse up
570	case 2: // Enter key
571	case 3: // Live update
572	Variable NumPts = sva.dval
573	String ns=CtrlNumber(sva.ctrlName) //control number as a string
574	NVAR thetaMin = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gAngLow"+ns)
575	NVAR incr = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gIncr"+ns)
576	NVAR thetaMax = $("root:Packages:NIST:USANS:Globals:U_Sim:"+"gAngHigh"+ns)
577
578	// old way, like ICP
579	// thetaMax = ThetaMin + incr*NumPts
580
581	// new way - to spread the points out over the specified angle range
582	incr = (thetaMax - thetaMin) /( numpts -1 ) //-1 since end points are both counted
583
584	CalcTotalCountTime()
585
586	break
587	endswitch
588
589	return 0
590	End
591
592	// changing count time -
593	// then recalculate the total count time
594	Function CtTimeSetVarProc(sva) : SetVariableControl
595	STRUCT WMSetVariableAction &sva
596
597	switch( sva.eventCode )
598	case 1: // mouse up
599	case 2: // Enter key
600	case 3: // Live update
601	CalcTotalCountTime()
602	break
603	endswitch
604
605	return 0
606	End
607
608	// hard-wired for 7 controls
609	// return value is in seconds
610	// global display string is set with hrs:min
611	Function CalcTotalCountTime()
612
613	Variable ii,num,totTime=0
614	String pathStr="root:Packages:NIST:USANS:Globals:U_Sim:"
615	num=7
616
617	for(ii=1;ii<=num;ii+=1)
618	NVAR ctTime = $(pathStr+"gCtTime"+num2str(ii))
619	NVAR numPts = $(pathStr+"gNumPts"+num2str(ii))
620	if(ctTime>0)
621	totTime += ctTime*numPts
622	endif
623	endfor
624	Variable hrs,mins
625	hrs = trunc(totTime/3600)
626	mins = trunc(mod(totTime,3600)/60)
627	// Printf "Counting time (hr:min) = %d:%d\r",hrs,mins
628
629	SVAR str = $(pathStr+"gTotTimeStr")
630	sprintf str,"%d:%d",hrs,mins
631
632	return(totTime)
633	End
634
635	//returns the control number from the name string
636	// all are setvarNa
637	Function/S CtrlNumber(str)
638	String str
639
640	return(str[6])
641	End
642
643	// changes edit mode of the theta min/max boxes and increment for simplified setup
644	// val = 2 = disable
645	// val = 0 = edit enabled
646	Function ThetaEditMode(val)
647	Variable val
648
649	SetVariable setvar1a,win=UCALC,disable=val
650	SetVariable setvar1b,win=UCALC,disable=val
651	SetVariable setvar1c,win=UCALC,disable=val
652
653	SetVariable setvar2a,win=UCALC,disable=val
654	SetVariable setvar2b,win=UCALC,disable=val
655	SetVariable setvar2c,win=UCALC,disable=val
656
657	SetVariable setvar3a,win=UCALC,disable=val
658	SetVariable setvar3b,win=UCALC,disable=val
659	SetVariable setvar3c,win=UCALC,disable=val
660
661	SetVariable setvar4a,win=UCALC,disable=val
662	SetVariable setvar4b,win=UCALC,disable=val
663	SetVariable setvar4c,win=UCALC,disable=val
664
665	SetVariable setvar5a,win=UCALC,disable=val
666	SetVariable setvar5b,win=UCALC,disable=val
667	SetVariable setvar5c,win=UCALC,disable=val
668
669	SetVariable setvar6a,win=UCALC,disable=val
670	SetVariable setvar6b,win=UCALC,disable=val
671	SetVariable setvar6c,win=UCALC,disable=val
672
673	SetVariable setvar7a,win=UCALC,disable=val
674	SetVariable setvar7b,win=UCALC,disable=val
675	SetVariable setvar7c,win=UCALC,disable=val
676	return(0)
677	End
678
679
680
681	Function EnterModeCheckProc(cba) : CheckBoxControl
682	STRUCT WMCheckboxAction &cba
683
684	switch( cba.eventCode )
685	case 2: // mouse up
686	Variable checked = cba.checked
687
688	if(checked)
689	ThetaEditMode(2)
690	else
691	ThetaEditMode(0)
692	endif
693
694	break
695	endswitch
696
697	return 0
698	End
699
700
701	Function ShowEMPCheckProc(cba) : CheckBoxControl
702	STRUCT WMCheckboxAction &cba
703
704	switch( cba.eventCode )
705	case 2: // mouse up
706	Variable checked = cba.checked
707
708	String list,item,popStr,qval,CR
709	Variable OK=1
710
711	if(exists("root:Packages:NIST:USANS:Globals:Q_2p0")==0) //
712	MakeUSANSEmptyWaves()
713	endif
714
715
716	if(checked)
717	// put it on the graph
718	SetDataFolder root:Packages:NIST:USANS:Globals
719	ControlInfo/W=UCALC popup0
720	popStr = S_Value
721	strswitch(popStr) // string switch
722	case "0.25": // execute if case matches expression
723	qval = "Q_0p25"
724	CR = "CR_0p25"
725	break
726	case "0.50":
727	qval = "Q_0p50"
728	CR = "CR_0p50"
729	break
730	case "0.625":
731	qval = "Q_0p625"
732	CR = "CR_0p625"
733	break
734	case "1.75":
735	qval = "Q_1p75"
736	CR = "CR_1p75"
737	break
738	case "2.0":
739	qval = "Q_2p0"
740	CR = "CR_2p0"
741	break
742	default: // optional default expression executed
743	OK=0
744	endswitch
745
746	if(OK)
747	AppendToGraph/W=UCALC $CR vs $qval
748	ModifyGraph marker=19,mode($CR)=4,msize($CR)=3,rgb($CR)=(0,0,0)
749	endif
750
751	else
752	//take it off of the graph
753	SetDataFolder root:Packages:NIST:USANS:Globals
754	list=WaveList("CR*", ";", "WIN:UCALC")
755	item=StringFromList(0, list ,";") //should be one item
756	if(strlen(item) != 0)
757	RemoveFromGraph/W=UCALC $item
758	endif
759	endif
760
761	break
762	endswitch
763
764	SetDataFolder root:
765	return 0
766	End
767
768
769
770	// based on the angle ranges above (with non-zero count times)
771	// plot where the data points would be
772	//
773	Function U_SimPlotButtonProc(ba) : ButtonControl
774	STRUCT WMButtonAction &ba
775
776
777	switch( ba.eventCode )
778	case 2: // mouse up
779	// click code here
780
781	CalcUSANS()
782
783	break
784	endswitch
785
786	return 0
787	End
788
789
790
791	// Fills a fake USANS folder with the concatenated ranges, and converts to Q
792	// root:Packages:NIST:USANS:SIM
793	//
794	// does the work of calculating and smearing the simluated USANS data set
795	Function CalcUSANS()
796
797	Variable num,ii,firstSet=0
798	String pathStr="root:Packages:NIST:USANS:Globals:U_Sim:gCtTime"
799	String fromType="SWAP",toType="SIM"
800	SVAR USANSFolder = root:Packages:NIST:USANS:Globals:gUSANSFolder
801
802	num = 7 //# of angular ranges
803
804	// only try to plot ranges with non-zero count times
805	firstSet=0
806	for(ii=1;ii<=num;ii+=1)
807	NVAR dum = $(pathStr+num2str(ii))
808	if(dum>0)
809	// print "CtTime = ",dum
810	firstSet += 1
811	LoadSimulatedAngleRange(ii,"SWAP") //overwrite what's in the SWAP folder
812
813	// did not do this step
814	//Convert2Countrate("SWAP",1)
815	//
816
817	if(firstSet==1) //first time, overwrite
818	NewDataWaves("SWAP","SIM")
819	//plus my two waves
820	Duplicate/O $(USANSFolder+":"+fromType+":countingTime"),$(USANSFolder+":"+toType+":countingTime")
821	Duplicate/O $(USANSFolder+":"+fromType+":SetNumber"),$(USANSFolder+":"+toType+":SetNumber")
822
823	else //append to waves in "SIM"
824	AppendDataWaves("SWAP","SIM")
825	// and my two waves
826	ConcatenateData( (USANSFolder+":"+toType+":countingTime"),(USANSFolder+":"+fromType+":countingTime") )
827	ConcatenateData( (USANSFolder+":"+toType+":SetNumber"),(USANSFolder+":"+fromType+":SetNumber") )
828
829	endif
830
831	endif
832	endfor
833
834
835	//sort after all loaded - not by angle, but by Q
836	// Get rid of the negative angles - the smearing integration does not like these!
837	// (may add them back in later, but probably not)
838	UDoAngleSort("SIM")
839
840	ConvertAngle2Qvals("SIM",0)
841
842
843
844	//fill the data with something
845	WAVE qvals = root:Packages:NIST:USANS:SIM:qvals
846	WAVE DetCts = root:Packages:NIST:USANS:SIM:DetCts
847
848
849	//find the Trans Cts for T_Wide
850	// FindTWideCts("EMP")
851
852	//generate a "fake" 1d data folder/set named "Sim_USANS" that can be used for smearing
853	// DOES NOT calculate a matrix, but instead fills 4-5-6 columns with -dQv
854	// so that the trapezoid rule is used
855	//
856
857	FakeUSANSDataFolder(qvals,DetCts,0.117,"Sim_USANS")
858
859	// now calculate the smearing... instead of the counts
860	SVAR funcStr = root:Packages:NIST:USANS:Globals:U_Sim:gFuncStr //set by the popup
861
862	Wave inten = root:Sim_USANS:Sim_USANS_i
863	Wave sigave = root:Sim_USANS:Sim_USANS_s
864	Wave countingTime = root:Sim_USANS:countingTime //counting time per point in the set
865
866	Duplicate/O inten root:Sim_USANS:Smeared_inten // a place for the smeared result for probabilities
867	Wave Smeared_inten = root:Sim_USANS:Smeared_inten
868
869	String coefStr=""
870	Variable sig_sas=0,wavelength = 2.4
871	Variable Imon
872
873
874	NVAR omega = root:Packages:NIST:USANS:Globals:U_Sim:gAnalyzerOmega
875
876	if(exists(funcStr) != 0)
877	FUNCREF SANSModelAAO_proto func=$("fSmeared"+funcStr) //a wrapper for the structure version
878	FUNCREF SANSModelAAO_proto funcUnsmeared=$(funcStr) //unsmeared
879	coefStr = MC_getFunctionCoef(funcStr)
880
881	if(!MC_CheckFunctionAndCoef(funcStr,coefStr))
882	Abort "Function and coefficients do not match. You must plot the unsmeared function before simulation."
883	endif
884
885	// do the smearing calculation
886	func($coefStr,Smeared_inten,qvals)
887
888
889	NVAR thick = root:Packages:NIST:USANS:Globals:U_Sim:gThick
890	NVAR trans = root:Packages:NIST:USANS:Globals:U_Sim:gSamTrans
891	NVAR SimDetCts = root:Packages:NIST:USANS:Globals:U_Sim:g_1DTotCts //summed counts (simulated)
892	NVAR estDetCR = root:Packages:NIST:USANS:Globals:U_Sim:g_1DEstDetCR // estimated detector count rate
893	NVAR fracScat = root:Packages:NIST:USANS:Globals:U_Sim:g_1DFracScatt // fraction of beam captured on detector
894	NVAR estTrans = root:Packages:NIST:USANS:Globals:U_Sim:g_1DEstTrans // estimated transmission of sample
895	// NVAR SimCountTime = root:Packages:NIST:USANS:Globals:U_Sim:gCntTime //counting time used for simulation
896
897	Imon = GetUSANSBeamIntensity() //based on the aperture size, select the beam intensity
898	Print "imon=",imon
899
900	// calculate the scattering cross section simply to be able to estimate the transmission
901	// unfortunately, this calculation is useless in the USANS range. So although it's calculated, the
902	// results are never reported.
903	// -- the main issue is that the integration range is optimized for SANS, and is not useful for USANS
904	// -there are only a few points in the USANS range, and extending it lower caused issues before.
905	CalculateRandomDeviate(funcUnsmeared,$coefStr,wavelength,"root:Packages:NIST:SAS:ran_dev",sig_sas)
906
907	if(sig_sas > 100)
908	DoAlert 0,"SAS cross section > 100. Estimates of multiple scattering are unreliable. Choosing a model with a well-defined Rg may help"
909	endif
910	// if(sig_sas > 100)
911	// sprintf abortStr,"sig_sas = %g. Please check that the model coefficients have a zero background, or the low q is well-behaved.",sig_sas
912	// endif
913	// calculate the multiple scattering fraction for display (10/2009)
914	NVAR mScat = root:Packages:NIST:USANS:Globals:U_Sim:g_MultScattFraction
915	Variable nMax=10,tau
916	mScat=0
917	tau = thick*sig_sas
918	// this sums the normalized scattering P', so the result is the fraction of multiply coherently scattered
919	// neutrons out of those that were scattered
920	for(ii=2;ii<nMax;ii+=1)
921	mScat += tau^(ii)/factorial(ii)
922	// print tau^(ii)/factorial(ii)
923	endfor
924	estTrans = exp(-1thicksig_sas) //thickness and sigma both in units of cm
925	mscat *= (estTrans)/(1-estTrans)
926
927	if(mScat > 0.1)
928	SetVariable valdisp0_2 win=UCALC,labelBack=(65535,32768,32768)
929	else
930	SetVariable valdisp0_2 win=UCALC,labelBack=0
931	endif
932
933
934	Print "Sig_sas = ",sig_sas
935	Print "Estimated Transmission = ",estTrans
936
937	Duplicate/O qvals prob_i
938
939	prob_i = transthickomega*Smeared_inten //probability of a neutron in q-bin(i)
940
941	// Variable P_on = sum(prob_i,-inf,inf)
942	// Print "P_on = ",P_on
943	fracScat = 1-estTrans
944
945	inten = (ImoncountingTime)prob_i
946
947	// do I round to an integer?
948	inten = round(inten)
949
950	// SimDetCts = sum(inten,-inf,inf)
951	// estDetCR = SimDetCts/SimCountTime
952
953
954	NVAR doABS = root:Packages:NIST:USANS:Globals:U_Sim:g_1D_DoABS
955	NVAR plotCR = root:Packages:NIST:USANS:Globals:U_Sim:g_1D_PlotCR
956	NVAR addNoise = root:Packages:NIST:USANS:Globals:U_Sim:g_1D_AddNoise
957
958	sigave = sqrt(inten) // assuming that N is large
959
960	// add in random error in aveint based on the sigave
961	if(addNoise)
962	inten += gnoise(sigave)
963
964	//round to an integer again
965	inten = round(inten)
966	endif
967
968	// convert to absolute scale? Maybe not needed
969	// does nothing yet - need Ipeak, Twide
970	// if(doABS)
971	// Variable kappa = thickomegatransiMonctTime
972	// inten /= kappa
973	// inten /= kappa
974	// endif
975
976	// plot as countrate - maybe easier to visualize, and all of the data overlaps
977	if(plotCR)
978	inten /= countingTime
979	sigave /= countingTime
980	endif
981
982	GraphSIM()
983
984	else
985	//no function plotted, no simulation can be done
986	DoAlert 0,"No function is selected or plotted, so no simulation is done. The default power law function is used."
987
988	inten = U_Power_Law_Model(1e-6,3,0,qvals)
989	// inten = U_SphereForm(1,9000,6e-6,0,qvals)
990
991	GraphSIM()
992
993	endif
994
995
996	end
997
998
999	//sort the data in the "type"folder, based on angle
1000	//carry along all associated waves
1001	//
1002	// ---a duplicate of DoAngleSort(), by modified to
1003	// include counting time and setNumber
1004	//
1005	// also trims the beginning of each data set so that it does not include any negative or zero angles
1006	//
1007	Function UDoAngleSort(type)
1008	String type
1009
1010	SVAR USANSFolder = root:Packages:NIST:USANS:Globals:gUSANSFolder
1011
1012	Wave angle = $(USANSFolder+":"+Type+":Angle")
1013	Wave detCts = $(USANSFolder+":"+Type+":DetCts")
1014	Wave ErrdetCts = $(USANSFolder+":"+Type+":ErrDetCts")
1015	Wave MonCts = $(USANSFolder+":"+Type+":MonCts")
1016	Wave TransCts = $(USANSFolder+":"+Type+":TransCts")
1017	Wave countingTime = $(USANSFolder+":"+Type+":countingTime")
1018	Wave SetNumber = $(USANSFolder+":"+Type+":SetNumber")
1019
1020	Sort Angle DetCts,ErrDetCts,MonCts,TransCts,Angle,countingTime,SetNumber
1021
1022	Variable ii,num,numBad,ang,val
1023	num=numpnts(angle)
1024	ii=0
1025	numBad=0
1026	val = 0 //cutoff value
1027	do
1028	ang = angle[ii]
1029	if(ang <= val)
1030	numBad += 1
1031	else //keep the points
1032	Angle[ii-numBad] = ang
1033	DetCts[ii-numBad] = DetCts[ii]
1034	ErrDetCts[ii-numBad] = ErrDetCts[ii]
1035	MonCts[ii-numBad] = MonCts[ii]
1036	TransCts[ii-numBad] = TransCts[ii]
1037	countingTime[ii-numBad] = countingTime[ii]
1038	SetNumber[ii-numBad] = SetNumber[ii]
1039	endif
1040	ii += 1
1041	while(ii<num)
1042	//trim the end of the waves
1043	DeletePoints num-numBad, numBad, DetCts,ErrDetCts,MonCts,TransCts,Angle,countingTime,SetNumber
1044
1045
1046	return(0)
1047	End
1048
1049	// a simple default function
1050	//
1051	Function U_SphereForm(scale,radius,delrho,bkg,x)
1052	Variable scale,radius,delrho,bkg
1053	Variable x
1054
1055	// variables are:
1056	//[0] scale
1057	//[1] radius (A)
1058	//[2] delrho (A-2)
1059	//[3] background (cm-1)
1060
1061
1062	// calculates scale * f^2/Vol where f=Vol3delrho*((sin(qr)-qrcos(qr))/qr^3
1063	// and is rescaled to give [=] cm^-1
1064
1065	Variable bes,f,vol,f2
1066	//
1067	//handle q==0 separately
1068	If(x==0)
1069	f = 4/3piradius^3delrhodelrhoscale1e8 + bkg
1070	return(f)
1071	Endif
1072
1073	bes = 3(sin(xradius)-xradiuscos(x*radius))/x^3/radius^3
1074	vol = 4pi/3radius^3
1075	f = volbesdelrho // [=] A
1076	// normalize to single particle volume, convert to 1/cm
1077	f2 = f * f / vol * 1.0e8 // [=] 1/cm
1078
1079	return (scale*f2+bkg) // Scale, then add in the background
1080
1081	End
1082
1083	// better default function
1084	Function U_Power_Law_Model(A,m,bgd,x) : FitFunc
1085	Variable A, m,bgd,x
1086	// Input (fitting) variables are:
1087	//[0] Coefficient
1088	//[1] (-) Power
1089	//[2] incoherent background
1090
1091	// local variables
1092	Variable inten, qval
1093	// x is the q-value for the calculation
1094	qval = x
1095	// do the calculation and return the function value
1096
1097	inten = A*qval^-m + bgd
1098	Return (inten)
1099	End
1100
1101
1102	// mimics LoadBT5File
1103	// creates two other waves to identify the set and the counting time for that set
1104	//
1105	Function LoadSimulatedAngleRange(set,type)
1106	Variable set
1107	String type
1108
1109	SVAR USANSFolder = root:Packages:NIST:USANS:Globals:gUSANSFolder
1110
1111	String s=num2str(set)
1112	NVAR angLow = $("root:Packages:NIST:USANS:Globals:U_Sim:gAngLow"+s)
1113	NVAR angHigh = $("root:Packages:NIST:USANS:Globals:U_Sim:gAngHigh"+s)
1114	NVAR numPts = $("root:Packages:NIST:USANS:Globals:U_Sim:gNumPts"+s)
1115	NVAR ctTime = $("root:Packages:NIST:USANS:Globals:U_Sim:gCtTime"+s)
1116	NVAR incr = $("root:Packages:NIST:USANS:Globals:U_Sim:gIncr"+s)
1117
1118	Variable ii, err=0
1119
1120	// generate q-points based on angular range from panel
1121
1122	Make/O/D/N=(numPts) $(USANSFolder+":"+type+":Angle")
1123	Make/O/D/N=(numPts) $(USANSFolder+":"+type+":DetCts")
1124	Make/O/D/N=(numPts) $(USANSFolder+":"+type+":ErrDetCts")
1125	Make/O/D/N=(numPts) $(USANSFolder+":"+type+":MonCts")
1126	Make/O/D/N=(numPts) $(USANSFolder+":"+type+":TransCts")
1127	//
1128	Make/O/D/N=(numPts) $(USANSFolder+":"+type+":CountingTime")
1129	Make/O/D/N=(numPts) $(USANSFolder+":"+type+":SetNumber")
1130
1131	Wave Angle = $(USANSFolder+":"+type+":Angle")
1132	Wave DetCts = $(USANSFolder+":"+type+":DetCts")
1133	Wave ErrDetCts = $(USANSFolder+":"+type+":ErrDetCts")
1134	Wave MonCts = $(USANSFolder+":"+type+":MonCts")
1135	Wave TransCts = $(USANSFolder+":"+type+":TransCts")
1136	Wave countingTime = $(USANSFolder+":"+type+":countingTime")
1137	Wave SetNumber = $(USANSFolder+":"+type+":SetNumber")
1138
1139	countingTime = ctTime
1140	SetNumber = set
1141
1142	for(ii=0;ii<numPts;ii+=1)
1143	Angle[ii] = angLow + ii*incr
1144	DetCts[ii] = set
1145	endfor
1146
1147	//set the wave note for the DetCts
1148	String str=""
1149	str = "FILE:Sim "+s+";"
1150	str += "TIMEPT:"+num2str(ctTime)+";"
1151	str += "PEAKANG:0;" //no value yet
1152	str += "STATUS:;" //no value yet
1153	str += "LABEL:SimSet "+s+";"
1154	str += "PEAKVAL:;" //no value yet
1155	str += "TWIDE:0;" //no value yet
1156	Note DetCts,str
1157
1158	return err // Zero signifies no error.
1159	End
1160
1161	// add SIM data to the graph if it exists and is not already on the graph
1162	//
1163	//
1164	// ** currently, I have changed the graph to not display the angle at the top. I have not yet
1165	// found a way to keep the scaling of the two axes in-sync when the empty cell data is added to the
1166	// graph (versus Q)
1167	//
1168	Function GraphSIM()
1169
1170	// SVAR USANSFolder = root:Packages:NIST:USANS:Globals:gUSANSFolder
1171	// SetDataFolder $(USANSFolder+":SIM")
1172
1173	SetDataFolder root:Sim_USANS
1174	//is it already on the graph?
1175	String list=""
1176	list = Wavelist("Sim_USANS*",";","WIN:UCALC")
1177
1178	if(strlen(list)!=0)
1179	// Print "SIM already on graph"
1180	SetDataFolder root:
1181	return(0)
1182	Endif
1183
1184	//append the data if it exists
1185	If(waveExists($"Sim_USANS_i")==1)
1186	DoWindow/F UCALC
1187
1188	// lines for the no-noise result vs angle
1189	// AppendToGraph/T Smeared_inten vs angle
1190	AppendToGraph Smeared_inten vs Sim_USANS_q
1191	ModifyGraph rgb(Smeared_inten)=(1,12815,52428)
1192	ModifyGraph mode(Smeared_inten)=4,marker(Smeared_inten)=19,msize(Smeared_inten)=2
1193	ModifyGraph tickUnit=1
1194
1195	// colored points for the simulation with noise on top
1196	AppendToGraph Sim_USANS_i vs Sim_USANS_q
1197	ModifyGraph mode(Sim_USANS_i)=3,marker(Sim_USANS_i)=19,msize(Sim_USANS_i)=4
1198	//don't reverse the dbZ21 or the highest angle will have "invisible" light blue error bars
1199	ModifyGraph zColor(Sim_USANS_i)={setNumber,1,7,dBZ21,0} //better for 7 colors
1200	// ModifyGraph zColor(Sim_USANS_i)={setNumber,1,7,Rainbow,0} //force the colors from 1->7
1201	ModifyGraph useMrkStrokeRGB(Sim_USANS_i)=1
1202	ErrorBars/T=0 Sim_USANS_i Y,wave=(Sim_USANS_s,Sim_USANS_s)
1203
1204	ModifyGraph log=1
1205	ModifyGraph mirror(left)=1
1206	// if no top axis, then
1207	// ModifyGraph mirror(bottom)=1
1208	ModifyGraph grid=2
1209	ModifyGraph standoff=0
1210
1211	// to make sure that the scales are the same (but fails on zoom)
1212	NewFreeAxis/O/T top_angle
1213	ModifyFreeAxis/Z top_angle,master= bottom,hook= TransformAngleAxisHook
1214	ModifyGraph lblPos(top_angle)=50,freePos(top_angle)=0
1215	ModifyGraph log(top_angle)=1
1216
1217	SetDrawEnv linefgc= (39321,1,1),dash= 3,linethick= 3
1218	SetDrawEnv xcoord= bottom
1219	DrawLine 3e-05,0.01,3e-05,0.99
1220
1221	Label top_angle "Angle"
1222	Label bottom "Q (1/A)"
1223	Label left "Counts or Count Rate"
1224
1225	Legend
1226	endif
1227
1228	SetDataFolder root:
1229	End
1230
1231	Function TransformAngleAxisHook(s)
1232	STRUCT WMAxisHookStruct &s
1233
1234	s.max= s.max/5.55e-5
1235	s.min= s.min/5.55e-5
1236
1237	return 0
1238	End
1239
1240	//fakes a folder with loaded 1-d usans data, no calculation of the matrix
1241	Function FakeUSANSDataFolder(qval,aveint,dqv,dataFolder)
1242	WAVE qval,aveint
1243	Variable dqv
1244	String dataFolder
1245
1246
1247	String baseStr=dataFolder
1248	if(DataFolderExists("root:"+baseStr))
1249	SetDataFolder $("root:"+baseStr)
1250	else
1251	NewDataFolder/S $("root:"+baseStr)
1252	endif
1253
1254	////overwrite the existing data, if it exists
1255	Duplicate/O qval, $(baseStr+"_q")
1256	Duplicate/O aveint, $(baseStr+"_i")
1257
1258	Duplicate/O qval, $(baseStr+"_s")
1259	Wave sigave = $(baseStr+"_s")
1260
1261	sigave = sqrt(aveint)
1262
1263	// make a resolution matrix for SANS data
1264	Variable np=numpnts(qval)
1265	Make/D/O/N=(np,4) $(baseStr+"_res")
1266	Wave res=$(baseStr+"_res")
1267
1268	res[][0] = -dQv //sigQ
1269	res[][1] = -dQv //qBar
1270	res[][2] = -dQv //fShad
1271	res[][3] = qval[p] //Qvalues
1272
1273	// extra waves of set number and counting time for the simulation
1274	WAVE ctW = root:Packages:NIST:USANS:SIM:countingTime
1275	WAVE setW = root:Packages:NIST:USANS:SIM:setNumber
1276	WAVE ang = root:Packages:NIST:USANS:SIM:Angle
1277	Duplicate/O ctW countingTime
1278	Duplicate/O setW setNumber
1279	Duplicate/O ang Angle
1280
1281
1282	//clean up
1283	SetDataFolder root:
1284
1285	End
1286
1287
1288	Function/S U_FunctionPopupList()
1289	String list,tmp
1290	list = User_FunctionPopupList()
1291
1292	//simplify the display, forcing smeared calculations behind the scenes
1293	tmp = FunctionList("Smear*",";","NPARAMS:1") //smeared dependency calculations
1294	list = RemoveFromList(tmp, list ,";")
1295
1296
1297	if(strlen(list)==0)
1298	list = "No functions plotted"
1299	endif
1300
1301	list = SortList(list)
1302
1303	return(list)
1304	End
1305
1306	Function Sim_USANS_ModelPopMenuProc(pa) : PopupMenuControl
1307	STRUCT WMPopupAction &pa
1308
1309	switch( pa.eventCode )
1310	case 2: // mouse up
1311	Variable popNum = pa.popNum
1312	String popStr = pa.popStr
1313	SVAR gStr = root:Packages:NIST:USANS:Globals:U_Sim:gFuncStr
1314	gStr = popStr
1315
1316	break
1317	endswitch
1318
1319	return 0
1320	End
1321
1322	// if the sample aperture is changed, AND the empty data is displayed, change to the proper data
1323	Function UCALC_SampleAperturePopup(pa) : PopupMenuControl
1324	STRUCT WMPopupAction &pa
1325
1326	switch( pa.eventCode )
1327	case 2: // mouse up
1328	Variable popNum = pa.popNum
1329	String popStr = pa.popStr
1330
1331	ControlInfo/W=UCALC check0_4
1332	if(V_Value==1) //currently checked, need to update
1333
1334	STRUCT WMCheckboxAction cba
1335	cba.checked=0 //"un-check"
1336	cba.eventCode=2
1337	ShowEMPCheckProc(cba)
1338
1339	cba.checked=1 //"re-check"
1340	ShowEMPCheckProc(cba)
1341
1342	endif
1343	break
1344
1345	endswitch
1346
1347	return 0
1348	End
1349
1350
1351	Function UCALC_PresetPopup(pa) : PopupMenuControl
1352	STRUCT WMPopupAction &pa
1353
1354	switch( pa.eventCode )
1355	case 2: // mouse up
1356	Variable popNum = pa.popNum
1357	String popStr = pa.popStr
1358
1359	SetDataFolder root:Packages:NIST:USANS:Globals:U_Sim
1360
1361	NVAR gAngLow1 = gAngLow1
1362	NVAR gAngHigh1 = gAngHigh1
1363	NVAR gNumPts1 = gNumPts1
1364	NVAR gCtTime1 = gCtTime1
1365	NVAR gIncr1 = gIncr1
1366
1367	NVAR gAngLow2 = gAngLow2
1368	NVAR gAngHigh2 = gAngHigh2
1369	NVAR gNumPts2 = gNumPts2
1370	NVAR gCtTime2 = gCtTime2
1371	NVAR gIncr2 = gIncr2
1372
1373	NVAR gAngLow3 = gAngLow3
1374	NVAR gAngHigh3 = gAngHigh3
1375	NVAR gNumPts3 = gNumPts3
1376	NVAR gCtTime3 = gCtTime3
1377	NVAR gIncr3 = gIncr3
1378
1379	NVAR gAngLow4 = gAngLow4
1380	NVAR gAngHigh4 = gAngHigh4
1381	NVAR gNumPts4 = gNumPts4
1382	NVAR gCtTime4 = gCtTime4
1383	NVAR gIncr4 = gIncr4
1384
1385	NVAR gAngLow5 = gAngLow5
1386	NVAR gAngHigh5 = gAngHigh5
1387	NVAR gNumPts5 = gNumPts5
1388	NVAR gCtTime5 = gCtTime5
1389	NVAR gIncr5 = gIncr5
1390
1391	NVAR gAngLow6 = gAngLow6
1392	NVAR gAngHigh6 = gAngHigh6
1393	NVAR gNumPts6 = gNumPts6
1394	NVAR gCtTime6 = gCtTime6
1395	NVAR gIncr6 = gIncr6
1396
1397	NVAR gAngLow7 = gAngLow7
1398	NVAR gAngHigh7 = gAngHigh7
1399	NVAR gNumPts7 = gNumPts7
1400	NVAR gCtTime7 = gCtTime7
1401	NVAR gIncr7 = gIncr7
1402
1403	strswitch(popStr) // string switch
1404	case "Short Count": // execute if case matches expression
1405
1406	gAngLow1 = -1
1407	gAngHigh1 = 0.6
1408	gNumPts1 = 33
1409	gCtTime1 = 10
1410	gIncr1 = 0.05
1411
1412	gAngLow2 = 0.7
1413	gAngHigh2 = 1.9
1414	gNumPts2 = 6
1415	gCtTime2 = 60
1416	gIncr2 = 0.2
1417
1418	gAngLow3 = 2
1419	gAngHigh3 = 4.8
1420	gNumPts3 = 7
1421	gCtTime3 = 120
1422	gIncr3 = 0.4
1423
1424	gAngLow4 = 5
1425	gAngHigh4 = 9.5
1426	gNumPts4 = 5
1427	gCtTime4 = 180
1428	gIncr4 = 0.9
1429
1430	gAngLow5 = 10
1431	gAngHigh5 = 19
1432	gNumPts5 = 5
1433	gCtTime5 = 240
1434	gIncr5 = 1.8
1435
1436	gAngLow6 = 20
1437	gAngHigh6 = 48
1438	gNumPts6 = 15
1439	gCtTime6 = 0
1440	gIncr6 = 2
1441
1442	gAngLow7 = 50
1443	gAngHigh7 = 95
1444	gNumPts7 = 10
1445	gCtTime7 = 0
1446	gIncr7 = 5
1447
1448	break
1449	case "Medium Count":
1450
1451	gAngLow1 = -1
1452	gAngHigh1 = 0.6
1453	gNumPts1 = 33
1454	gCtTime1 = 20
1455	gIncr1 = 0.05
1456
1457	gAngLow2 = 0.7
1458	gAngHigh2 = 1.9
1459	gNumPts2 = 13
1460	gCtTime2 = 50
1461	gIncr2 = 0.1
1462
1463	gAngLow3 = 2
1464	gAngHigh3 = 4.8
1465	gNumPts3 = 15
1466	gCtTime3 = 100
1467	gIncr3 = 0.2
1468
1469	gAngLow4 = 5
1470	gAngHigh4 = 9.5
1471	gNumPts4 = 10
1472	gCtTime4 = 200
1473	gIncr4 = 0.5
1474
1475	gAngLow5 = 10
1476	gAngHigh5 = 19
1477	gNumPts5 = 10
1478	gCtTime5 = 500
1479	gIncr5 = 1
1480
1481	gAngLow6 = 20
1482	gAngHigh6 = 48
1483	gNumPts6 = 15
1484	gCtTime6 = 1000
1485	gIncr6 = 2
1486
1487	gAngLow7 = 50
1488	gAngHigh7 = 95
1489	gNumPts7 = 10
1490	gCtTime7 = 0
1491	gIncr7 = 5
1492
1493	break
1494	case "Long Count":
1495
1496	gAngLow1 = -1
1497	gAngHigh1 = 0.6
1498	gNumPts1 = 33
1499	gCtTime1 = 25
1500	gIncr1 = 0.05
1501
1502	gAngLow2 = 0.7
1503	gAngHigh2 = 1.9
1504	gNumPts2 = 13
1505	gCtTime2 = 100
1506	gIncr2 = 0.1
1507
1508	gAngLow3 = 2
1509	gAngHigh3 = 4.8
1510	gNumPts3 = 15
1511	gCtTime3 = 300
1512	gIncr3 = 0.2
1513
1514	gAngLow4 = 5
1515	gAngHigh4 = 9.5
1516	gNumPts4 = 10
1517	gCtTime4 = 600
1518	gIncr4 = 0.5
1519
1520	gAngLow5 = 10
1521	gAngHigh5 = 19
1522	gNumPts5 = 10
1523	gCtTime5 = 1200
1524	gIncr5 = 1
1525
1526	gAngLow6 = 20
1527	gAngHigh6 = 48
1528	gNumPts6 = 15
1529	gCtTime6 = 2000
1530	gIncr6 = 2
1531
1532	gAngLow7 = 50
1533	gAngHigh7 = 95
1534	gNumPts7 = 10
1535	gCtTime7 = 3000
1536	gIncr7 = 5
1537
1538	break
1539	default:
1540	endswitch
1541
1542	break
1543	endswitch
1544
1545	//update the count time
1546	CalcTotalCountTime()
1547
1548	SetDataFolder root:
1549
1550	return 0
1551	End
1552
1553
1554	// return the beam intensity based on the sample aperture diameter
1555	//
1556	// based on the equation in John's instrument paper
1557	//
1558	Function GetUSANSBeamIntensity()
1559
1560	String popStr
1561	Variable flux,diam,rad
1562
1563	ControlInfo/W=UCALC popup0
1564	popStr = S_Value
1565	diam=str2num(popStr) //in inches
1566	rad = diam/2*25.4 //radius in mm
1567
1568	flux = 662radrad-39.9radradrad+0.70radradrad*rad
1569
1570	return(flux)
1571	End
1572
1573	// based on the angle ranges above (with non-zero count times)
1574	// save fake data points into a fake BT5 data file
1575	//
1576	Function U_SaveButtonProc(ba) : ButtonControl
1577	STRUCT WMButtonAction &ba
1578
1579
1580	switch( ba.eventCode )
1581	case 2: // mouse up
1582	// click code here
1583
1584	Variable num,ii,baseNumber=101,firstSet=0
1585	String pathStr="root:Packages:NIST:USANS:Globals:U_Sim:gCtTime"
1586	SVAR USANSFolder = root:Packages:NIST:USANS:Globals:gUSANSFolder
1587	String baseName="SIMUL"
1588
1589	num = 7 //# of angular ranges
1590
1591	// only try to plot ranges with non-zero count times
1592	firstSet=0
1593	for(ii=1;ii<=num;ii+=1)
1594	NVAR dum = $(pathStr+num2str(ii))
1595	if(dum>0)
1596	firstSet += 1
1597	if(firstSet==1) //first time, ask for base name
1598	Prompt baseName,"Enter a base name for the files"
1599	Prompt baseNumber,"Enter the starting index"
1600	DoPrompt "Enter Information for Save",baseName,baseNumber
1601
1602	if(V_Flag==1) //user canceled
1603	return(1)
1604	endif
1605	endif
1606
1607	SaveFakeUSANS(baseName,baseNumber-1,ii)
1608
1609	endif
1610	endfor
1611
1612	// no save the fake "COR" file that is the combination and rescaling of the data sets
1613	// - the combination is in root:Sim_USANS:
1614	NVAR g_1D_PlotCR = root:Packages:NIST:USANS:Globals:U_Sim:g_1D_PlotCR
1615
1616	if(g_1D_PlotCR) // right now, save COR only if the CR is plotted
1617	SaveFakeCOR(baseName)
1618	endif
1619
1620	break
1621	endswitch
1622
1623	return 0
1624	End
1625
1626
1627	Function SaveFakeCOR(nameStr)
1628	String nameStr
1629
1630	NVAR gThick = root:Packages:NIST:USANS:Globals:U_Sim:gThick
1631	NVAR gSamTrans = root:Packages:NIST:USANS:Globals:U_Sim:gSamTrans
1632	NVAR gAnalyzerOmega = root:Packages:NIST:USANS:Globals:U_Sim:gAnalyzerOmega
1633	NVAR g_MultScattFraction = root:Packages:NIST:USANS:Globals:U_Sim:g_MultScattFraction
1634	SVAR funcStr = root:Packages:NIST:USANS:Globals:U_Sim:gFuncStr
1635
1636	String folder = "root:Sim_USANS:",fileStr=""
1637	String termStr="\r\n" //VAX uses only <CR> as terminator, but only CRLF seems to FTP correctly to VAX
1638	String destStr="",formatStr = "%15.6g %15.6g %15.6g %15.6g %15.6g %15.6g"+termStr
1639
1640	Variable refNum,integer,realval
1641	//***these waves MUST EXIST, or IGOR Pro will crash, with a type 2 error**
1642	WAVE qvals =$(folder + "Sim_USANS_q")
1643	WAVE inten=$(folder + "Sim_USANS_i")
1644	WAVE sig=$(folder + "Sim_USANS_s")
1645
1646	//check each wave
1647	If(!(WaveExists(qvals)))
1648	Abort "qvals DNExist in SaveFakeCOR()"
1649	Endif
1650	If(!(WaveExists(inten)))
1651	Abort "inten DNExist in SaveFakeCOR()"
1652	Endif
1653	If(!(WaveExists(sig)))
1654	Abort "sig DNExist in SaveFakeCOR()"
1655	Endif
1656
1657	//dummy wave for the resolution
1658	Duplicate/O qvals,dumWave
1659	SVAR USANSFolder = root:Packages:NIST:USANS:Globals:gUSANSFolder
1660	NVAR DQv=$(USANSFolder+":Globals:MainPanel:gDQv")
1661	dumWave = - DQv
1662
1663	//temp waves for the rescaling
1664	Duplicate/O qvals,tq,ti,te
1665	ti=inten
1666	te=sig
1667	tq=qvals
1668
1669	//rescale
1670	ti *= 1.0e6 / GetUSANSBeamIntensity() //converts to "SAM" scaled data
1671	ti = 1/(gSamTransgThickgAnalyzerOmega1.0e6) //converts to "COR" scaling (approx, since I don't know Twide)
1672	//rescale error appropriately too
1673	te *= 1.0e6 / GetUSANSBeamIntensity() //converts to "SAM" scaled data
1674	te = 1/(gSamTransgThickgAnalyzerOmega1.0e6) //converts to "COR" scaling (approx, since I don't know Twide)
1675
1676	// add fake error to every few data points to mark it as simulated data
1677	te[2,numpnts(te)-1;7] *= 10
1678
1679	// now trim off the q-values less than 3e-5
1680	if(tq[0] < 3e-5)
1681	do
1682	DeletePoints 0, 1, ti,tq,te,dumWave
1683	while(tq[0] < 3e-5)
1684	endif
1685
1686	String samStr="",empStr="",dateStr="",samLabelStr="",paramStr="",empLevStr="",bkgLevStr="",pkStr=""
1687	samStr = "SIM FILES: "+StringByKey("FILE",note(inten),":",";")
1688	empStr = "EMP FILES: "+"none"
1689	empLevStr = "EMP LEVEL: " + "none"
1690	bkgLevStr = "BKG LEVEL: " + "not used"
1691	paramStr = "Ds = "+num2str(gthick)+" cm ; "
1692	paramStr += "Twide = "+"unknown"+" ; "
1693	paramStr += "Trock = "+num2str(gSamTrans)
1694	pkStr += "Multiple coherent scattering fraction: "+num2str(g_MultScattFraction)
1695
1696	//these strings are always the same
1697	samLabelStr ="LABEL: "+"Simulated USANS Data from: "+funcStr
1698	dateStr = date()+" at "+time()// + " "+Secs2Time(DateTime,2)
1699
1700
1701	fileStr=nameStr+".cor"
1702	//actually open the file
1703	PathInfo savePathName
1704	if(V_flag==1)
1705	Open/P=savePathName refNum as fileStr
1706	else
1707	Open refNum as fileStr
1708	endif
1709
1710	fprintf refnum,"%s"+termStr,samStr
1711	fprintf refnum,"%s"+termStr,dateStr
1712	fprintf refnum,"%s"+termStr,samLabelStr
1713	fprintf refnum,"%s"+termStr,empStr
1714	fprintf refnum,"%s"+termStr,paramStr
1715	fprintf refnum,"%s"+termStr,pkStr
1716	fprintf refnum,"%s"+termStr,empLevStr + " ; "+bkglevStr
1717
1718	//
1719	wfprintf refnum, formatStr, tq,ti,te,dumWave,dumWave,dumWave
1720
1721	Close refnum
1722
1723	Killwaves/Z dumWave,ti,te,tq
1724
1725
1726
1727	return(0)
1728	End
1729
1730
1731	//
1732	// duplicates the BT5 file format
1733	//
1734	Function SaveFakeUSANS(nameStr,num,set)
1735	String nameStr
1736	Variable num,set
1737
1738	String folder = "root:Packages:NIST:USANS:SIM:"
1739	String pathStr="root:Packages:NIST:USANS:Globals:U_Sim:"
1740	String termStr="\r\n" //VAX uses only <CR> as terminator, but only CRLF seems to FTP correctly to VAX
1741
1742	// WAVE DetCts = $(folder+"DetCts") //these are only dummy values
1743	WAVE DetCts = root:Sim_USANS:Sim_USANS_i
1744	WAVE angle = $(folder+"Angle")
1745	WAVE SetNumber = $(folder+"SetNumber")
1746	WAVE countingTime = $(folder+"countingTime")
1747
1748	NVAR ang1 = $(pathStr+"gAngLow"+num2str(set))
1749	NVAR ang2 = $(pathStr+"gAngHigh"+num2str(set))
1750	NVAR incr = $(pathStr+"gIncr"+num2str(set))
1751	SVAR ext = root:Packages:NIST:USANS:Globals:MainPanel:gUExt
1752
1753
1754	Variable refNum,ii,wavePts,numPts,first,last,ctTime,monCt,transDet
1755	String str,fileStr,dateStr
1756
1757	wavePts = numpnts(angle)
1758	for(ii=0;ii<wavePts;ii+=1)
1759	if(setNumber[ii] == set)
1760	first = ii
1761	break
1762	endif
1763	endfor
1764
1765	for(ii=wavePts-1;ii>=0;ii-=1)
1766	if(setNumber[ii] == set)
1767	last = ii
1768	break
1769	endif
1770	endfor
1771
1772	// Print "First, last = ",first,last
1773
1774	// set up some of the strings needed
1775	fileStr=nameStr+num2str(num+set)+ext
1776	dateStr = date()
1777	ctTime = countingTime[first]
1778	numPts = last-first+1
1779
1780	MonCt = ctTime*GetUSANSBeamIntensity()
1781
1782	transDet = 999 //bogus and constant value
1783
1784	//actually open the file
1785	PathInfo savePathName
1786	if(V_flag==1)
1787	Open/P=savePathName refNum as fileStr
1788	else
1789	Open refNum as fileStr
1790	endif
1791
1792	sprintf str,"'%s' '%s' 'I' %d 1 'TIME' %d 'RAW'",fileStr,dateStr,ctTime,numpts
1793	fprintf refnum,"%s"+termStr,str
1794	fprintf refnum,"%s"+termStr," Filename Date Scan Mon Prf Base #pts Type"
1795
1796	fprintf refnum,"Simulated USANS data"+termStr
1797
1798	fprintf refnum,"%s"+termStr," 0 0 0 0 0 0 0 0.0000 0.00000 0.00000 0.00000 2"
1799	fprintf refnum,"%s"+termStr," Collimation Mosaic Wavelength T-Start Incr. H-field #Det "
1800	fprintf refnum,"%s"+termStr," 1 0.00000 0.00000 0.00000"
1801	fprintf refnum," 2 %9.5f %9.5f %9.5f"+termStr,ang1,incr,ang2
1802	fprintf refnum,"%s"+termStr," 3 10.00000 0.00000 10.00000"
1803	fprintf refnum,"%s"+termStr," 4 10.00000 0.00000 10.00000"
1804	fprintf refnum,"%s"+termStr," 5 0.00000 0.00000 0.00000"
1805	fprintf refnum,"%s"+termStr," 6 0.00000 0.00000 0.00000"
1806	fprintf refnum,"%s"+termStr," Mot: Start Step End"
1807	fprintf refnum,"%s"+termStr," A2 MIN MONITOR COUNTS EXTRA "
1808
1809	//loop over the waves, picking out the desired set
1810	//write 2 lines each time
1811	for(ii=first;ii<=last;ii+=1)
1812	sprintf str," %6.3f %6.2f %d %d %d",angle[ii],ctTime/60,MonCt,DetCts[ii]*ctTime,transDet
1813	fprintf refnum,"%s"+termStr,str
1814
1815	sprintf str,"%d,%d,0,%d,0,0,0,0",MonCt,DetCts[ii]*ctTime,transDet
1816	fprintf refnum,"%s"+termStr,str
1817	endfor
1818
1819	Close refnum
1820
1821	return(0)
1822	end
1823
1824
1825	// print out the USANS configuration in some reasonable format
1826	Function/S USANSConfigurationText()
1827
1828	String str="",temp
1829
1830	SetDataFolder root:Packages:NIST:USANS:Globals:U_Sim
1831
1832	// results, setup values
1833	SVAR gTotTimeStr=gTotTimeStr
1834
1835	Variable ii,num
1836	String pathStr="root:Packages:NIST:USANS:Globals:U_Sim:"
1837	num=7
1838
1839	str += "USANS Instrument Configuration:\r\r"
1840	str += "Theta Min Theta Max Increment # Points Count Time\r"
1841
1842
1843	for(ii=1;ii<=num;ii+=1)
1844	NVAR ctTime = $(pathStr+"gCtTime"+num2str(ii))
1845	if(ctTime>0)
1846	NVAR angLow = $(pathStr+"gAngLow"+num2str(ii))
1847	NVAR angHigh = $(pathStr+"gAngHigh"+num2str(ii))
1848	NVAR incr = $(pathStr+"gIncr"+num2str(ii))
1849	NVAR numPts = $(pathStr+"gNumPts"+num2str(ii))
1850
1851	sprintf temp,"%9.3f %9.3f %9.3f %9d %9d\r",angLow,angHigh,incr,numPts,ctTime
1852	str += temp
1853	endif
1854	endfor
1855
1856
1857	sprintf temp,"\r\rTotal Counting Time (HR:MIN) = %s\r",gTotTimeStr
1858	str += temp
1859
1860
1861	setDataFolder root:
1862	return str
1863	End
1864
1865	Function DisplayUCALCText()
1866
1867	if(WinType("USANS_Configuration")==0)
1868	NewNotebook/F=0/K=1/N=USANS_Configuration /W=(480,44,880,369)
1869	endif
1870	//replace the text
1871	Notebook USANS_Configuration selection={startOfFile, endOfFile}
1872	Notebook USANS_Configuration font="Monaco",fSize=10,text=USANSConfigurationText()
1873	return(0)
1874	end
1875
1876
1877	//
1878	Function U_ConfigTextProc(ba) : ButtonControl
1879	STRUCT WMButtonAction &ba
1880
1881
1882	switch( ba.eventCode )
1883	case 2: // mouse up
1884	// click code here
1885
1886	DisplayUCALCText()
1887	break
1888	endswitch
1889
1890	return 0
1891	End
1892
1893	// this will save a graphic of the whole panel that then needs to be opened and printed
1894	// must be a PNG @ screen resolution
1895	//
1896	Function U_SavePanelProc(ba) : ButtonControl
1897	STRUCT WMButtonAction &ba
1898
1899
1900	switch( ba.eventCode )
1901	case 2: // mouse up
1902	// click code here
1903
1904	SavePICT/P=home/E=-5/B=72/SNAP=1
1905
1906	// can I reload and print?
1907	// how will the users know where this went and what to do with it?
1908	//
1909	break
1910	endswitch
1911
1912	return 0
1913	End

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

Download in other formats:

Original Format