Context Navigation

V_IQ_Utilities.ipf @ 1064

Last change on this file since 1064 was 1064, checked in by srkline, 5 years ago

many minor changes after real VSANS data collected.

additional procedures added to allow easy correction of the incorrect header information from NICE.

Most notable addition is the pinhole resolution added to the calculation and the I(q) output. White beam is also treated (incorrectly) as a gaussian distrivution, but the results of smeared fitting look to be quite good.

Trimming and sorting routines are now (pinhole) resolution aware.

File identification routines have been updated to use the proper definitions of "purpose" and "intent". Both fields are now in the catalog, to allow for better sorting.

File size: 34.4 KB

Rev	Line
[984]	1	#pragma TextEncoding = "MacRoman" // For details execute DisplayHelpTopic "The TextEncoding Pragma"
	2	#pragma rtGlobals=3 // Use modern global access method and strict wave access.
	3
[1023]	4
[984]	5	//
[1023]	6	// Operation does no scaling, only the basic (default) trim of the ends, concatenate, sort, and save
[1030]	7	// -- if data has been converted to WORK and the solid angle correction was done, then the data
[1023]	8	// is per unit solid angle, and matches up - at least the simulated data does...
[1030]	9	// It should match up in real VSANS data since the flux conditions are identical for
	10	// all panels, only the geometry is different.
[984]	11	//
[999]	12	//
[1023]	13	// V_DataPlotting.ipf is where the I(q) panel is drawn and the binning is set
	14	//
[999]	15	// see the VCALC BinAllMiddlePanels() for an example of this
	16	// see the binning routines in VC_DetectorBinning_Utils.ipf for the details
	17	//
[1023]	18
[999]	19	// TODO
[1022]	20	//
	21	// -- verify the binning for slit mode. Looks correct, but verify
	22	// -- DOCUMENT
	23	//
[999]	24	// x- detector "B" is currently skipped since the calibration waves are not faked
	25	// when the raw data is loaded. Then the qxqyqz waves are not generated.
	26	//
[1022]	27	// x- REDO the logic here. It's a mess, and will get the calculation wrong
[999]	28	//
[1022]	29	// x- figure out the binning type (where is it set for VSANS?)
	30	// x- don't know, so currently VSANS binning type is HARD-WIRED
	31	// x- figure out when this needs to be called to (force) re-calculate I vs Q
[999]	32	//
[1022]	33
	34
[1036]	35
[1038]	36	//
[1035]	37	// NOTE
	38	// this is the master conversion function
	39	// ***Use no others
	40	// *** When other bin types are developed, DO NOT reassign these numbers.
	41	// instead, skip the old numbers and assign new ones.
	42	// old modes can be removed from the string constant ksBinTypeStr (above), but the
	43	// mode numbers are what many different binning, plotting, and reduction functions are
	44	// switching on. In the future, it may be necessary to change the key (everywhere) to a string
	45	// switch, but for now, stick with the numbers.
	46	Function V_BinTypeStr2Num(binStr)
	47	String binStr
	48
	49	Variable binType
	50	strswitch(binStr) // string switch
	51	case "One":
	52	binType = 1
	53	break // exit from switch
	54	case "Two":
	55	binType = 2
	56	break // exit from switch
	57	case "Four":
	58	binType = 3
	59	break // exit from switch
	60	case "Slit Mode":
	61	binType = 4
	62	break // exit from switch
	63
	64	default: // optional default expression executed
	65	binType = 0
	66	Abort "Binning mode not found"// when no case matches
	67	endswitch
	68	return(binType)
	69	end
	70
[1044]	71	//
	72	// TODO -- binType == 4 (slit mode) should never end up here
	73	// -- new logic in calling routines to dispatch to proper routine
	74	// -- AND need to write the routine for binning_SlitMode
	75	//
	76	Function V_QBinAllPanels_Circular(folderStr,binType)
[999]	77	String folderStr
[1027]	78	Variable binType
[999]	79
	80	// do the back, middle, and front separately
	81
	82	// figure out the binning type (where is it set?)
[1027]	83	Variable ii,delQ
[999]	84	String detStr
	85
[1027]	86	// binType = V_GetBinningPopMode()
[1022]	87
[999]	88	//// TODO:
[1027]	89	//
[1035]	90	// Back detector is handled spearately since there is nothing to combine
	91	//
[999]	92	delQ = SetDeltaQ(folderStr,"B")
	93
	94	// dispatch based on binning type
[1027]	95	if(binType == 1 \|\| binType == 2 \|\| binType == 3)
[999]	96	VC_fDoBinning_QxQy2D(folderStr, "B") //normal binning, nothing to combine
	97	endif
	98
	99	// TODO -- this is only a temporary fix for slit mode
	100	if(binType == 4)
	101	/// this is for a tall, narrow slit mode
	102	VC_fBinDetector_byRows(folderStr,"B")
	103	endif
	104
	105
	106
	107	// these are the binning types where detectors are not combined
	108	// other combined binning is below the loop
	109	for(ii=0;ii<ItemsInList(ksDetectorListNoB);ii+=1)
	110	detStr = StringFromList(ii, ksDetectorListNoB, ";")
	111
	112	// set delta Q for binning
	113	delQ = SetDeltaQ(folderStr,detStr)
	114
	115	// dispatch based on binning type
	116	if(binType==1)
	117	VC_fDoBinning_QxQy2D(folderStr,detStr)
	118	endif
	119
	120	// TODO -- this is only a temporary fix for slit mode
	121	if(binType == 4)
	122	/// this is for a tall, narrow slit mode
	123	VC_fBinDetector_byRows(folderStr,detStr)
	124	endif
	125
	126	endfor
	127
	128	// bin in pairs
	129	if(binType == 2)
	130	VC_fDoBinning_QxQy2D(folderStr,"MLR")
	131	VC_fDoBinning_QxQy2D(folderStr,"MTB")
	132	VC_fDoBinning_QxQy2D(folderStr,"FLR")
	133	VC_fDoBinning_QxQy2D(folderStr,"FTB")
	134	endif
	135
	136	// bin everything on front or middle together
	137	if(binType == 3)
	138	VC_fDoBinning_QxQy2D(folderStr,"MLRTB")
	139	VC_fDoBinning_QxQy2D(folderStr,"FLRTB")
	140	endif
	141
	142	return(0)
	143	End
	144
[1027]	145	// concatenates and sorts the 1D data in "type" WORK folder
	146	// uses the current display if type==""
	147	//
[1039]	148	Function V_ConcatenateForSave(pathStr,type,tagStr,binType)
	149	String pathStr,type,tagStr
[1027]	150	Variable binType
	151
	152	// get the current display type, if null string passed in
	153	SVAR curtype = root:Packages:NIST:VSANS:Globals:gCurDispType
	154
	155	if(strlen(type)==0)
	156	type = curType
	157	endif
[999]	158
[984]	159	// trim the data if needed
[1002]	160	// remove the q=0 point from the back detector, if it's there
[1027]	161	// does not need to know binType
[1002]	162	V_RemoveQ0_B(type)
[984]	163
	164	// concatenate the data sets
[1023]	165	// TODO x- figure out which binning was used (this is done in V_1DConcatenate())
[1002]	166	// clear the old tmp waves first, if they still exist
[1039]	167	// SetDataFolder $("root:Packages:NIST:VSANS:"+type)
	168	SetDataFolder $(pathStr+type)
[1064]	169	Killwaves/Z tmp_q,tmp_i,tmp_s,tmp_sq,tmp_qb,tmp_fs
[1002]	170	setDataFolder root:
[1039]	171	V_1DConcatenate(pathStr,type,tagStr,binType)
[984]	172
	173	// sort the data set
[1039]	174	V_TmpSort1D(pathStr,type)
[984]	175
[1027]	176	return(0)
	177	End
	178
	179	//
[1057]	180	// this is only called from the button on the data panel (not anymore)
[1027]	181	// so the type is the currently displayed type, and the binning is from the panel
	182	//
[1039]	183	Function V_SimpleSave1DData(pathStr,type,tagStr,saveName)
	184	String pathStr,type,tagStr,saveName
[1027]	185
	186	//
	187	// get the current display type, if null string passed in
	188	SVAR curtype = root:Packages:NIST:VSANS:Globals:gCurDispType
	189	Variable binType = V_GetBinningPopMode()
	190
[1039]	191	V_ConcatenateForSave(pathStr,curType,tagStr,binType)
[1027]	192
[984]	193	// write out the data set to a file
[1027]	194	if(strlen(saveName)==0)
	195	Execute "V_GetNameForSave()"
	196	SVAR newName = root:saveName
	197	saveName = newName
	198	endif
	199
[1039]	200	V_Write1DData(pathStr,curtype,saveName)
[984]	201
[999]	202	End
[984]	203
[1027]	204
[999]	205	Proc V_GetNameForSave(str)
	206	String str
[1027]	207	String/G root:saveName=str
[984]	208	End
	209
	210
[1002]	211	// blindly assumes that there is only one zero at the top of the wave
	212	// could be more sophisticated in the future...
	213	Function V_RemoveQ0_B(type)
	214	String type
	215
	216	SetDataFolder $("root:Packages:NIST:VSANS:"+type)
	217
	218	WAVE/Z qBin = qBin_qxqy_B
	219	WAVE/Z iBin = iBin_qxqy_B
	220	WAVE/Z eBin = eBin_qxqy_B
	221	WAVE/Z nBin = nBin_qxqy_B
	222	WAVE/Z iBin2 = iBin2_qxqy_B
	223
[1064]	224	// resolution waves
	225	Wave/Z sigQ = sigmaQ_qxqy_B
	226	Wave/Z qBar = qBar_qxqy_B
	227	Wave/Z fSubS = fSubS_qxqy_B
	228
[1002]	229	if(qBin[0] == 0)
[1064]	230	DeletePoints 0, 1, qBin,iBin,eBin,nBin,iBin2,sigQ,qBar,fSubS
[1002]	231	endif
	232
	233	SetDataFolder root:
	234	return(0)
	235	end
	236
	237
[984]	238	// concatentate data in folderStr
	239	//
	240	// TODO:
[1051]	241	// -- !!! Resolution waves are currently skipped - these must be added
	242	//
[1023]	243	// x- this currently ignores the binning type (one, two, etc. )
[1036]	244	// x- change the Concatenate call to use the waveList, to eliminate the need to declare all of the waves
[984]	245	// -- this currently assumes that all of the waves exist
	246	// -- need robust error checking for wave existence
	247	// -- wave names are hard-wired and their name and location may be different in the future
[1023]	248	// x- if different averaging options were chosen (bin type of 2, 4 etc) then
[984]	249	// although waves may exist, they may not be the right ones to use. There
	250	// will be a somewhat complex selection process
[1002]	251	// x- detector B is currently skipped
[984]	252	//
[1023]	253	// this seems like a lot of extra work to do something so simple...but it's better than a loop
[984]	254	//
	255	// root:Packages:NIST:VSANS:RAW:iBin_qxqy_FB
[1023]	256	//
	257	// binType = 1 = one
	258	// binType = 2 = two
	259	// binType = 3 = four
	260	// binType = 4 = Slit Mode
	261	//
[1027]	262	// if binType is passed in as -9999, get the binning mode from the popup
	263	// otherwise the value is assumed good (from a protocol)
	264	//
[1039]	265	// pathStr must have the trailing colon
	266	// tagStr is normally null, but is "_trim" for data to be trimmed
[1036]	267	//
[1039]	268	Function V_1DConcatenate(pathStr,folderStr,tagStr,binType)
	269	String pathStr,folderStr,tagStr
[1027]	270	Variable binType
[984]	271
[1027]	272	if(binType==-9999)
	273	binType = V_GetBinningPopMode()
	274	endif
[1023]	275
[1039]	276	// SetDataFolder $("root:Packages:NIST:VSANS:"+folderStr)
	277	SetDataFolder $(pathStr+folderStr)
[1023]	278
[1031]	279	//kill these waves before starting, or the new concatenation will be added to the old
[1064]	280	KillWaves/Z tmp_q,tmp_i,tmp_s,tmp_qb,tmp_sq,tmp_fs
[1051]	281
	282	NVAR gIgnoreDetB = root:Packages:NIST:VSANS:Globals:gIgnoreDetB
[1031]	283
[1036]	284	String waveListStr=""
	285	if(binType == 1)
	286	// q-values
[1051]	287	waveListStr=""
	288	if(!gIgnoreDetB)
	289	waveListStr = "qBin_qxqy_B" + tagStr + ";"
	290	endif
[1039]	291	waveListStr += "qBin_qxqy_MB" + tagStr + ";"
	292	waveListStr += "qBin_qxqy_MT" + tagStr + ";"
	293	waveListStr += "qBin_qxqy_ML" + tagStr + ";"
	294	waveListStr += "qBin_qxqy_MR" + tagStr + ";"
	295	waveListStr += "qBin_qxqy_FB" + tagStr + ";"
	296	waveListStr += "qBin_qxqy_FT" + tagStr + ";"
	297	waveListStr += "qBin_qxqy_FL" + tagStr + ";"
	298	waveListStr += "qBin_qxqy_FR" + tagStr + ";"
[1036]	299
	300	Concatenate/NP/O waveListStr, tmp_q
	301
	302	//intensity
[1051]	303	waveListStr=""
	304	if(!gIgnoreDetB)
	305	waveListStr = "iBin_qxqy_B" + tagStr + ";"
	306	endif
[1039]	307	waveListStr += "iBin_qxqy_MB" + tagStr + ";"
	308	waveListStr += "iBin_qxqy_MT" + tagStr + ";"
	309	waveListStr += "iBin_qxqy_ML" + tagStr + ";"
	310	waveListStr += "iBin_qxqy_MR" + tagStr + ";"
	311	waveListStr += "iBin_qxqy_FB" + tagStr + ";"
	312	waveListStr += "iBin_qxqy_FT" + tagStr + ";"
	313	waveListStr += "iBin_qxqy_FL" + tagStr + ";"
	314	waveListStr += "iBin_qxqy_FR" + tagStr + ";"
	315	// waveListStr = "iBin_qxqy_B;iBin_qxqy_MB;iBin_qxqy_MT;iBin_qxqy_ML;iBin_qxqy_MR;"
	316	// waveListStr += "iBin_qxqy_FB;iBin_qxqy_FT;iBin_qxqy_FL;iBin_qxqy_FR;"
[1036]	317
	318	Concatenate/NP/O waveListStr, tmp_i
	319
	320	//error
[1051]	321	waveListStr=""
	322	if(!gIgnoreDetB)
	323	waveListStr = "eBin_qxqy_B" + tagStr + ";"
	324	endif
[1039]	325	waveListStr += "eBin_qxqy_MB" + tagStr + ";"
	326	waveListStr += "eBin_qxqy_MT" + tagStr + ";"
	327	waveListStr += "eBin_qxqy_ML" + tagStr + ";"
	328	waveListStr += "eBin_qxqy_MR" + tagStr + ";"
	329	waveListStr += "eBin_qxqy_FB" + tagStr + ";"
	330	waveListStr += "eBin_qxqy_FT" + tagStr + ";"
	331	waveListStr += "eBin_qxqy_FL" + tagStr + ";"
	332	waveListStr += "eBin_qxqy_FR" + tagStr + ";"
	333	// waveListStr = "eBin_qxqy_B;eBin_qxqy_MB;eBin_qxqy_MT;eBin_qxqy_ML;eBin_qxqy_MR;"
	334	// waveListStr += "eBin_qxqy_FB;eBin_qxqy_FT;eBin_qxqy_FL;eBin_qxqy_FR;"
[1036]	335
	336	Concatenate/NP/O waveListStr, tmp_s
[1064]	337
	338	//sigma Q
	339	waveListStr=""
	340	if(!gIgnoreDetB)
	341	waveListStr = "sigmaQ_qxqy_B" + tagStr + ";"
	342	endif
	343	waveListStr += "sigmaQ_qxqy_MB" + tagStr + ";"
	344	waveListStr += "sigmaQ_qxqy_MT" + tagStr + ";"
	345	waveListStr += "sigmaQ_qxqy_ML" + tagStr + ";"
	346	waveListStr += "sigmaQ_qxqy_MR" + tagStr + ";"
	347	waveListStr += "sigmaQ_qxqy_FB" + tagStr + ";"
	348	waveListStr += "sigmaQ_qxqy_FT" + tagStr + ";"
	349	waveListStr += "sigmaQ_qxqy_FL" + tagStr + ";"
	350	waveListStr += "sigmaQ_qxqy_FR" + tagStr + ";"
	351	// waveListStr = "sigmaQ_qxqy_B;sigmaQ_qxqy_MB;sigmaQ_qxqy_MT;sigmaQ_qxqy_ML;sigmaQ_qxqy_MR;"
	352	// waveListStr += "sigmaQ_qxqy_FB;sigmaQ_qxqy_FT;sigmaQ_qxqy_FL;sigmaQ_qxqy_FR;"
	353
	354	Concatenate/NP/O waveListStr, tmp_sq
	355
	356	//Q bar
	357	waveListStr=""
	358	if(!gIgnoreDetB)
	359	waveListStr = "qBar_qxqy_B" + tagStr + ";"
	360	endif
	361	waveListStr += "qBar_qxqy_MB" + tagStr + ";"
	362	waveListStr += "qBar_qxqy_MT" + tagStr + ";"
	363	waveListStr += "qBar_qxqy_ML" + tagStr + ";"
	364	waveListStr += "qBar_qxqy_MR" + tagStr + ";"
	365	waveListStr += "qBar_qxqy_FB" + tagStr + ";"
	366	waveListStr += "qBar_qxqy_FT" + tagStr + ";"
	367	waveListStr += "qBar_qxqy_FL" + tagStr + ";"
	368	waveListStr += "qBar_qxqy_FR" + tagStr + ";"
	369	// waveListStr = "qBar_qxqy_B;qBar_qxqy_MB;qBar_qxqy_MT;qBar_qxqy_ML;qBar_qxqy_MR;"
	370	// waveListStr += "qBar_qxqy_FB;qBar_qxqy_FT;qBar_qxqy_FL;qBar_qxqy_FR;"
	371
	372	Concatenate/NP/O waveListStr, tmp_qb
	373
	374	//shadow fs
	375	waveListStr=""
	376	if(!gIgnoreDetB)
	377	waveListStr = "fSubS_qxqy_B" + tagStr + ";"
	378	endif
	379	waveListStr += "fSubS_qxqy_MB" + tagStr + ";"
	380	waveListStr += "fSubS_qxqy_MT" + tagStr + ";"
	381	waveListStr += "fSubS_qxqy_ML" + tagStr + ";"
	382	waveListStr += "fSubS_qxqy_MR" + tagStr + ";"
	383	waveListStr += "fSubS_qxqy_FB" + tagStr + ";"
	384	waveListStr += "fSubS_qxqy_FT" + tagStr + ";"
	385	waveListStr += "fSubS_qxqy_FL" + tagStr + ";"
	386	waveListStr += "fSubS_qxqy_FR" + tagStr + ";"
	387	// waveListStr = "fSubS_qxqy_B;fSubS_qxqy_MB;fSubS_qxqy_MT;fSubS_qxqy_ML;fSubS_qxqy_MR;"
	388	// waveListStr += "fSubS_qxqy_FB;fSubS_qxqy_FT;fSubS_qxqy_FL;fSubS_qxqy_FR;"
	389
	390	Concatenate/NP/O waveListStr, tmp_fs
	391
[1036]	392	endif
	393
	394	if(binType == 2)
	395	// q-values
[1051]	396	waveListStr=""
	397	if(!gIgnoreDetB)
	398	waveListStr = "qBin_qxqy_B" + tagStr + ";"
	399	endif
[1039]	400	waveListStr += "qBin_qxqy_MTB" + tagStr + ";"
	401	waveListStr += "qBin_qxqy_MLR" + tagStr + ";"
	402	waveListStr += "qBin_qxqy_FTB" + tagStr + ";"
	403	waveListStr += "qBin_qxqy_FLR" + tagStr + ";"
[1036]	404
[1039]	405	// waveListStr = "qBin_qxqy_B;qBin_qxqy_MTB;qBin_qxqy_MLR;"
	406	// waveListStr += "qBin_qxqy_FTB;qBin_qxqy_FLR;"
	407
[1036]	408	Concatenate/NP/O waveListStr, tmp_q
	409
	410	//intensity
[1051]	411	waveListStr=""
	412	if(!gIgnoreDetB)
	413	waveListStr = "iBin_qxqy_B" + tagStr + ";"
	414	endif
[1039]	415	waveListStr += "iBin_qxqy_MTB" + tagStr + ";"
	416	waveListStr += "iBin_qxqy_MLR" + tagStr + ";"
	417	waveListStr += "iBin_qxqy_FTB" + tagStr + ";"
	418	waveListStr += "iBin_qxqy_FLR" + tagStr + ";"
[1036]	419
[1039]	420	// waveListStr = "iBin_qxqy_B;iBin_qxqy_MTB;iBin_qxqy_MLR;"
	421	// waveListStr += "iBin_qxqy_FTB;iBin_qxqy_FLR;"
	422
[1036]	423	Concatenate/NP/O waveListStr, tmp_i
	424
	425	//error
[1051]	426	waveListStr=""
	427	if(!gIgnoreDetB)
	428	waveListStr = "eBin_qxqy_B" + tagStr + ";"
	429	endif
[1039]	430	waveListStr += "eBin_qxqy_MTB" + tagStr + ";"
	431	waveListStr += "eBin_qxqy_MLR" + tagStr + ";"
	432	waveListStr += "eBin_qxqy_FTB" + tagStr + ";"
	433	waveListStr += "eBin_qxqy_FLR" + tagStr + ";"
	434
	435	// waveListStr = "eBin_qxqy_B;eBin_qxqy_MTB;eBin_qxqy_MLR;"
	436	// waveListStr += "eBin_qxqy_FTB;eBin_qxqy_FLR;"
[1036]	437
	438	Concatenate/NP/O waveListStr, tmp_s
[1064]	439
	440	// sigma Q
	441	waveListStr=""
	442	if(!gIgnoreDetB)
	443	waveListStr = "sigmaQ_qxqy_B" + tagStr + ";"
	444	endif
	445	waveListStr += "sigmaQ_qxqy_MTB" + tagStr + ";"
	446	waveListStr += "sigmaQ_qxqy_MLR" + tagStr + ";"
	447	waveListStr += "sigmaQ_qxqy_FTB" + tagStr + ";"
	448	waveListStr += "sigmaQ_qxqy_FLR" + tagStr + ";"
	449
	450	// waveListStr = "sigmaQ_qxqy_B;sigmaQ_qxqy_MTB;sigmaQ_qxqy_MLR;"
	451	// waveListStr += "sigmaQ_qxqy_FTB;sigmaQ_qxqy_FLR;"
	452
	453	Concatenate/NP/O waveListStr, tmp_sq
	454
	455	// Q bar
	456	waveListStr=""
	457	if(!gIgnoreDetB)
	458	waveListStr = "qBar_qxqy_B" + tagStr + ";"
	459	endif
	460	waveListStr += "qBar_qxqy_MTB" + tagStr + ";"
	461	waveListStr += "qBar_qxqy_MLR" + tagStr + ";"
	462	waveListStr += "qBar_qxqy_FTB" + tagStr + ";"
	463	waveListStr += "qBar_qxqy_FLR" + tagStr + ";"
	464
	465	// waveListStr = "qBar_qxqy_B;qBar_qxqy_MTB;qBar_qxqy_MLR;"
	466	// waveListStr += "qBar_qxqy_FTB;qBar_qxqy_FLR;"
	467
	468	Concatenate/NP/O waveListStr, tmp_qb
	469
	470	// shadow fs
	471	waveListStr=""
	472	if(!gIgnoreDetB)
	473	waveListStr = "fSubS_qxqy_B" + tagStr + ";"
	474	endif
	475	waveListStr += "fSubS_qxqy_MTB" + tagStr + ";"
	476	waveListStr += "fSubS_qxqy_MLR" + tagStr + ";"
	477	waveListStr += "fSubS_qxqy_FTB" + tagStr + ";"
	478	waveListStr += "fSubS_qxqy_FLR" + tagStr + ";"
	479
	480	// waveListStr = "fSubS_qxqy_B;fSubS_qxqy_MTB;fSubS_qxqy_MLR;"
	481	// waveListStr += "fSubS_qxqy_FTB;fSubS_qxqy_FLR;"
	482
	483	Concatenate/NP/O waveListStr, tmp_fs
	484
[1036]	485	endif
	486
	487	if(binType == 3)
	488	// q-values
[1051]	489	waveListStr=""
	490	if(!gIgnoreDetB)
	491	waveListStr = "qBin_qxqy_B" + tagStr + ";"
	492	endif
[1039]	493	waveListStr += "qBin_qxqy_MLRTB" + tagStr + ";"
	494	waveListStr += "qBin_qxqy_FLRTB" + tagStr + ";"
	495
	496	// waveListStr = "qBin_qxqy_B;qBin_qxqy_MLRTB;qBin_qxqy_FLRTB;"
[1036]	497
	498	Concatenate/NP/O waveListStr, tmp_q
	499
	500	//intensity
[1051]	501	waveListStr=""
	502	if(!gIgnoreDetB)
	503	waveListStr = "iBin_qxqy_B" + tagStr + ";"
	504	endif
[1039]	505	waveListStr += "iBin_qxqy_MLRTB" + tagStr + ";"
	506	waveListStr += "iBin_qxqy_FLRTB" + tagStr + ";"
[1036]	507
[1039]	508	// waveListStr = "iBin_qxqy_B;iBin_qxqy_MLRTB;iBin_qxqy_FLRTB;"
	509
[1036]	510	Concatenate/NP/O waveListStr, tmp_i
	511
	512	//error
[1051]	513	waveListStr=""
	514	if(!gIgnoreDetB)
	515	waveListStr = "eBin_qxqy_B" + tagStr + ";"
	516	endif
[1039]	517	waveListStr += "eBin_qxqy_MLRTB" + tagStr + ";"
	518	waveListStr += "eBin_qxqy_FLRTB" + tagStr + ";"
	519
	520	// waveListStr = "eBin_qxqy_B;eBin_qxqy_MLRTB;eBin_qxqy_FLRTB;"
[1036]	521
	522	Concatenate/NP/O waveListStr, tmp_s
[1064]	523
	524	// sigma Q
	525	waveListStr=""
	526	if(!gIgnoreDetB)
	527	waveListStr = "sigmaQ_qxqy_B" + tagStr + ";"
	528	endif
	529	waveListStr += "sigmaQ_qxqy_MLRTB" + tagStr + ";"
	530	waveListStr += "sigmaQ_qxqy_FLRTB" + tagStr + ";"
	531
	532	// waveListStr = "sigmaQ_qxqy_B;sigmaQ_qxqy_MLRTB;sigmaQ_qxqy_FLRTB;"
	533
	534	Concatenate/NP/O waveListStr, tmp_sq
	535
	536	// Q bar
	537	waveListStr=""
	538	if(!gIgnoreDetB)
	539	waveListStr = "qBar_qxqy_B" + tagStr + ";"
	540	endif
	541	waveListStr += "qBar_qxqy_MLRTB" + tagStr + ";"
	542	waveListStr += "qBar_qxqy_FLRTB" + tagStr + ";"
	543
	544	// waveListStr = "qBar_qxqy_B;qBar_qxqy_MLRTB;qBar_qxqy_FLRTB;"
	545
	546	Concatenate/NP/O waveListStr, tmp_qb
	547
	548	// shadow fs
	549	waveListStr=""
	550	if(!gIgnoreDetB)
	551	waveListStr = "fSubS_qxqy_B" + tagStr + ";"
	552	endif
	553	waveListStr += "fSubS_qxqy_MLRTB" + tagStr + ";"
	554	waveListStr += "fSubS_qxqy_FLRTB" + tagStr + ";"
	555
	556	// waveListStr = "fSubS_qxqy_B;fSubS_qxqy_MLRTB;fSubS_qxqy_FLRTB;"
	557
	558	Concatenate/NP/O waveListStr, tmp_fs
[1036]	559	endif
	560
	561	// TODO - This is the identical set of waves as for the case of binType = 1.
	562	// they have the same names, but are averaged differently since it's slit mode.
	563	// I have separated this, since in practice the TB panels are probably best to ignore
	564	// and NOT include in the averaging since the Qy range is so limited.
	565	if(binType == 4)
	566	// q-values
[1051]	567	waveListStr=""
	568	if(!gIgnoreDetB)
	569	waveListStr = "qBin_qxqy_B" + tagStr + ";"
	570	endif
	571	// waveListStr += "qBin_qxqy_MB" + tagStr + ";"
	572	// waveListStr += "qBin_qxqy_MT" + tagStr + ";"
[1039]	573	waveListStr += "qBin_qxqy_ML" + tagStr + ";"
	574	waveListStr += "qBin_qxqy_MR" + tagStr + ";"
[1051]	575	// waveListStr += "qBin_qxqy_FB" + tagStr + ";"
	576	// waveListStr += "qBin_qxqy_FT" + tagStr + ";"
[1039]	577	waveListStr += "qBin_qxqy_FL" + tagStr + ";"
	578	waveListStr += "qBin_qxqy_FR" + tagStr + ";"
	579	// waveListStr = "qBin_qxqy_B;qBin_qxqy_MB;qBin_qxqy_MT;qBin_qxqy_ML;qBin_qxqy_MR;"
	580	// waveListStr += "qBin_qxqy_FB;qBin_qxqy_FT;qBin_qxqy_FL;qBin_qxqy_FR;"
[1036]	581
	582	Concatenate/NP/O waveListStr, tmp_q
	583
	584	//intensity
[1051]	585	waveListStr=""
	586	if(!gIgnoreDetB)
	587	waveListStr = "iBin_qxqy_B" + tagStr + ";"
	588	endif
	589	// waveListStr += "iBin_qxqy_MB" + tagStr + ";"
	590	// waveListStr += "iBin_qxqy_MT" + tagStr + ";"
[1039]	591	waveListStr += "iBin_qxqy_ML" + tagStr + ";"
	592	waveListStr += "iBin_qxqy_MR" + tagStr + ";"
[1051]	593	// waveListStr += "iBin_qxqy_FB" + tagStr + ";"
	594	// waveListStr += "iBin_qxqy_FT" + tagStr + ";"
[1039]	595	waveListStr += "iBin_qxqy_FL" + tagStr + ";"
	596	waveListStr += "iBin_qxqy_FR" + tagStr + ";"
	597	// waveListStr = "iBin_qxqy_B;iBin_qxqy_MB;iBin_qxqy_MT;iBin_qxqy_ML;iBin_qxqy_MR;"
	598	// waveListStr += "iBin_qxqy_FB;iBin_qxqy_FT;iBin_qxqy_FL;iBin_qxqy_FR;"
[1036]	599
	600	Concatenate/NP/O waveListStr, tmp_i
	601
	602	//error
[1051]	603	waveListStr=""
	604	if(!gIgnoreDetB)
	605	waveListStr = "eBin_qxqy_B" + tagStr + ";"
	606	endif
	607	// waveListStr += "eBin_qxqy_MB" + tagStr + ";"
	608	// waveListStr += "eBin_qxqy_MT" + tagStr + ";"
[1039]	609	waveListStr += "eBin_qxqy_ML" + tagStr + ";"
	610	waveListStr += "eBin_qxqy_MR" + tagStr + ";"
[1051]	611	// waveListStr += "eBin_qxqy_FB" + tagStr + ";"
	612	// waveListStr += "eBin_qxqy_FT" + tagStr + ";"
[1039]	613	waveListStr += "eBin_qxqy_FL" + tagStr + ";"
	614	waveListStr += "eBin_qxqy_FR" + tagStr + ";"
	615	// waveListStr = "eBin_qxqy_B;eBin_qxqy_MB;eBin_qxqy_MT;eBin_qxqy_ML;eBin_qxqy_MR;"
	616	// waveListStr += "eBin_qxqy_FB;eBin_qxqy_FT;eBin_qxqy_FL;eBin_qxqy_FR;"
[1036]	617
	618	Concatenate/NP/O waveListStr, tmp_s
[1064]	619
	620	//sigma Q
	621	waveListStr=""
	622	if(!gIgnoreDetB)
	623	waveListStr = "sigmaQ_qxqy_B" + tagStr + ";"
	624	endif
	625	// waveListStr += "sigmaQ_qxqy_MB" + tagStr + ";"
	626	// waveListStr += "sigmaQ_qxqy_MT" + tagStr + ";"
	627	waveListStr += "sigmaQ_qxqy_ML" + tagStr + ";"
	628	waveListStr += "sigmaQ_qxqy_MR" + tagStr + ";"
	629	// waveListStr += "sigmaQ_qxqy_FB" + tagStr + ";"
	630	// waveListStr += "sigmaQ_qxqy_FT" + tagStr + ";"
	631	waveListStr += "sigmaQ_qxqy_FL" + tagStr + ";"
	632	waveListStr += "sigmaQ_qxqy_FR" + tagStr + ";"
	633	// waveListStr = "sigmaQ_qxqy_B;sigmaQ_qxqy_MB;sigmaQ_qxqy_MT;sigmaQ_qxqy_ML;sigmaQ_qxqy_MR;"
	634	// waveListStr += "sigmaQ_qxqy_FB;sigmaQ_qxqy_FT;sigmaQ_qxqy_FL;sigmaQ_qxqy_FR;"
	635
	636	Concatenate/NP/O waveListStr, tmp_sq
	637
	638	//Q bar
	639	waveListStr=""
	640	if(!gIgnoreDetB)
	641	waveListStr = "qBar_qxqy_B" + tagStr + ";"
	642	endif
	643	// waveListStr += "qBar_qxqy_MB" + tagStr + ";"
	644	// waveListStr += "qBar_qxqy_MT" + tagStr + ";"
	645	waveListStr += "qBar_qxqy_ML" + tagStr + ";"
	646	waveListStr += "qBar_qxqy_MR" + tagStr + ";"
	647	// waveListStr += "qBar_qxqy_FB" + tagStr + ";"
	648	// waveListStr += "qBar_qxqy_FT" + tagStr + ";"
	649	waveListStr += "qBar_qxqy_FL" + tagStr + ";"
	650	waveListStr += "qBar_qxqy_FR" + tagStr + ";"
	651	// waveListStr = "qBar_qxqy_B;qBar_qxqy_MB;qBar_qxqy_MT;qBar_qxqy_ML;qBar_qxqy_MR;"
	652	// waveListStr += "qBar_qxqy_FB;qBar_qxqy_FT;qBar_qxqy_FL;qBar_qxqy_FR;"
	653
	654	Concatenate/NP/O waveListStr, tmp_qb
	655
	656	//shadow fs
	657	waveListStr=""
	658	if(!gIgnoreDetB)
	659	waveListStr = "fSubS_qxqy_B" + tagStr + ";"
	660	endif
	661	// waveListStr += "fSubS_qxqy_MB" + tagStr + ";"
	662	// waveListStr += "fSubS_qxqy_MT" + tagStr + ";"
	663	waveListStr += "fSubS_qxqy_ML" + tagStr + ";"
	664	waveListStr += "fSubS_qxqy_MR" + tagStr + ";"
	665	// waveListStr += "fSubS_qxqy_FB" + tagStr + ";"
	666	// waveListStr += "fSubS_qxqy_FT" + tagStr + ";"
	667	waveListStr += "fSubS_qxqy_FL" + tagStr + ";"
	668	waveListStr += "fSubS_qxqy_FR" + tagStr + ";"
	669	// waveListStr = "fSubS_qxqy_B;fSubS_qxqy_MB;fSubS_qxqy_MT;fSubS_qxqy_ML;fSubS_qxqy_MR;"
	670	// waveListStr += "fSubS_qxqy_FB;fSubS_qxqy_FT;fSubS_qxqy_FL;fSubS_qxqy_FR;"
	671
	672	Concatenate/NP/O waveListStr, tmp_fs
	673
[1036]	674	endif
	675
	676	// Can't kill here, since they are still needed to sort and write out!
	677	// KillWaves/Z tmp_q,tmp_i,tmp_s,tmp_res0,tmp_res1,tmp_res2,tmp_res3
	678
	679	SetDataFolder root:
	680
	681	return(0)
	682	End
	683
	684	// TODO:
	685	// -- resolution waves are ignored, since they don't exist (yet)
	686	// -- only a sort is done, no rescaling of data sets
	687	// (it's too late now anyways, since the data was concatenated)
	688	//
	689	// see Auto_Sort() in the SANS Automation ipf for the rest of the details of
	690	// how to combine the resolution waves (they also need to be concatenated, which is currently not done)
	691	//
[1039]	692	Function V_TmpSort1D(pathStr,folderStr)
	693	String pathStr,folderStr
[1036]	694
[1039]	695	SetDataFolder $(pathStr+folderStr)
[1036]	696
	697	Wave qw = tmp_q
	698	Wave iw = tmp_i
	699	Wave sw = tmp_s
[1064]	700	Wave sq = tmp_sq
	701	Wave qb = tmp_qb
	702	Wave fs = tmp_fs
[1036]	703
	704
[1064]	705	Sort qw, qw,iw,sw,sq,qb,fs
[1036]	706
	707
	708	SetDataFolder root:
	709	return(0)
	710	End
	711
	712
	713	// TODO
[1039]	714	// (appears to be unused, in favor of the version that uses the global strings)
[1036]	715	// needs:
	716	// -- trim the beamstop out (based on shadow?)
	717	// -- trim out zero q from the file (bad actor in analysis functions)
	718	// -- trim num from the highQ end or lowQ end?
	719	// -- splits the res wave into individual waves in anticipation of concatenation
	720	// -- or -- deal with the res wave after?
	721	//
	722	// -- make a copy of the waves?
	723	// -- then, what is the concatenate function looking for??
	724	//
	725	Function V_Trim1DData(dataFolder,binType,nBeg,nEnd)
	726	String dataFolder
	727	Variable binType,nBeg,nEnd
	728
	729	Variable npt,ii
	730	SetDataFolder $("root:Packages:NIST:VSANS:"+dataFolder)
	731
	732	Printf "%d points removed from beginning, %d points from the end (of each set) before concatenating\r",nbeg,nend
	733
	734	// for each binType block:
	735	// declare the waves
	736	// make a copy of the waves??
	737	// //Break out resolution wave into separate waves
	738	// delete the beginning points from everything
	739	// trim off the last nEnd points from everything
	740	// DeletePoints num-nEnd,nEnd, qw,iw,sw
	741	// // delete all points where the shadow is < 0.98
	742	////Put resolution contents back???
	743
[1023]	744	if(binType == 1)
	745	Wave/Z q_fb = qBin_qxqy_FB
	746	Wave/Z q_ft = qBin_qxqy_FT
	747	Wave/Z q_fl = qBin_qxqy_FL
	748	Wave/Z q_fr = qBin_qxqy_FR
	749	Wave/Z q_mb = qBin_qxqy_MB
	750	Wave/Z q_mt = qBin_qxqy_MT
	751	Wave/Z q_ml = qBin_qxqy_ML
	752	Wave/Z q_mr = qBin_qxqy_MR
	753	Wave/Z q_b = qBin_qxqy_B
[984]	754
[1023]	755	Wave/Z i_fb = iBin_qxqy_FB
	756	Wave/Z i_ft = iBin_qxqy_FT
	757	Wave/Z i_fl = iBin_qxqy_FL
	758	Wave/Z i_fr = iBin_qxqy_FR
	759	Wave/Z i_mb = iBin_qxqy_MB
	760	Wave/Z i_mt = iBin_qxqy_MT
	761	Wave/Z i_ml = iBin_qxqy_ML
	762	Wave/Z i_mr = iBin_qxqy_MR
	763	Wave/Z i_b = iBin_qxqy_B
	764
	765	Wave/Z s_fb = eBin_qxqy_FB
	766	Wave/Z s_ft = eBin_qxqy_FT
	767	Wave/Z s_fl = eBin_qxqy_FL
	768	Wave/Z s_fr = eBin_qxqy_FR
	769	Wave/Z s_mb = eBin_qxqy_MB
	770	Wave/Z s_mt = eBin_qxqy_MT
	771	Wave/Z s_ml = eBin_qxqy_ML
	772	Wave/Z s_mr = eBin_qxqy_MR
	773	Wave/Z s_b = eBin_qxqy_B
[1038]	774
	775
[1036]	776
	777	DeletePoints 0,nBeg, q_fb,q_ft,q_fl,q_fr,q_mb,q_mt,q_ml,q_mr,q_b
	778	DeletePoints 0,nBeg, i_fb,i_ft,i_fl,i_fr,i_mb,i_mt,i_ml,i_mr,i_b
	779	DeletePoints 0,nBeg, s_fb,s_ft,s_fl,s_fr,s_mb,s_mt,s_ml,s_mr,s_b
	780	//since each set may have a different number of points
	781	npt = numpnts(q_fb)
	782	DeletePoints npt-nEnd,nEnd, q_fb,i_fb,s_fb
	783
	784	npt = numpnts(q_ft)
	785	DeletePoints npt-nEnd,nEnd, q_ft,i_ft,s_ft
	786
	787	npt = numpnts(q_fl)
	788	DeletePoints npt-nEnd,nEnd, q_fl,i_fl,s_fl
	789
	790	npt = numpnts(q_fr)
	791	DeletePoints npt-nEnd,nEnd, q_fr,i_fr,s_fr
	792
	793	npt = numpnts(q_mb)
	794	DeletePoints npt-nEnd,nEnd, q_mb,i_mb,s_mb
	795
	796	npt = numpnts(q_mt)
	797	DeletePoints npt-nEnd,nEnd, q_mt,i_mt,s_mt
	798
	799	npt = numpnts(q_ml)
	800	DeletePoints npt-nEnd,nEnd, q_ml,i_ml,s_ml
	801
	802	npt = numpnts(q_mr)
	803	DeletePoints npt-nEnd,nEnd, q_mr,i_mr,s_mr
	804
	805	npt = numpnts(q_b)
	806	DeletePoints npt-nEnd,nEnd, q_b,i_b,s_b
[1023]	807
	808	endif
[984]	809
[1023]	810	if(binType == 2)
	811	Wave/Z q_ftb = qBin_qxqy_FTB
	812	Wave/Z q_flr = qBin_qxqy_FLR
	813	Wave/Z q_mtb = qBin_qxqy_MTB
	814	Wave/Z q_mlr = qBin_qxqy_MLR
	815	Wave/Z q_b = qBin_qxqy_B
	816
	817	Wave/Z i_ftb = iBin_qxqy_FTB
	818	Wave/Z i_flr = iBin_qxqy_FLR
	819	Wave/Z i_mtb = iBin_qxqy_MTB
	820	Wave/Z i_mlr = iBin_qxqy_MLR
	821	Wave/Z i_b = iBin_qxqy_B
[1036]	822
[1023]	823	Wave/Z s_ftb = eBin_qxqy_FTB
	824	Wave/Z s_flr = eBin_qxqy_FLR
	825	Wave/Z s_mtb = eBin_qxqy_MTB
	826	Wave/Z s_mlr = eBin_qxqy_MLR
	827	Wave/Z s_b = eBin_qxqy_B
	828
[1036]	829
	830	DeletePoints 0,nBeg, q_ftb,q_flr,q_mtb,q_mlr,q_b
	831	DeletePoints 0,nBeg, i_ftb,i_flr,i_mtb,i_mlr,i_b
	832	DeletePoints 0,nBeg, s_ftb,s_flr,s_mtb,s_mlr,s_b
	833	//since each set may have a different number of points
	834	npt = numpnts(q_ftb)
	835	DeletePoints npt-nEnd,nEnd, q_ftb,i_ftb,s_ftb
	836
	837	npt = numpnts(q_flr)
	838	DeletePoints npt-nEnd,nEnd, q_flr,i_flr,s_flr
	839
	840	npt = numpnts(q_mtb)
	841	DeletePoints npt-nEnd,nEnd, q_mtb,i_mtb,s_mtb
	842
	843	npt = numpnts(q_mlr)
	844	DeletePoints npt-nEnd,nEnd, q_mlr,i_mlr,s_mlr
	845
	846	npt = numpnts(q_b)
	847	DeletePoints npt-nEnd,nEnd, q_b,i_b,s_b
	848
	849
[1023]	850	endif
[984]	851
[1023]	852	if(binType == 3)
	853	Wave/Z q_flrtb = qBin_qxqy_FLRTB
	854	Wave/Z q_mlrtb = qBin_qxqy_MLRTB
	855	Wave/Z q_b = qBin_qxqy_B
[984]	856
[1023]	857	Wave/Z i_flrtb = iBin_qxqy_FLRTB
	858	Wave/Z i_mlrtb = iBin_qxqy_MLRTB
[1036]	859	Wave/Z i_b = iBin_qxqy_B
[1023]	860
	861	Wave/Z s_flrtb = eBin_qxqy_FLRTB
	862	Wave/Z s_mlrtb = eBin_qxqy_MLRTB
	863	Wave/Z s_b = eBin_qxqy_B
	864
[1036]	865	DeletePoints 0,nBeg, q_flrtb,q_mlrtb,q_b
	866	DeletePoints 0,nBeg, i_flrtb,i_mlrtb,i_b
	867	DeletePoints 0,nBeg, s_flrtb,s_mlrtb,s_b
	868	//since each set may have a different number of points
	869	npt = numpnts(q_flrtb)
	870	DeletePoints npt-nEnd,nEnd, q_flrtb,i_flrtb,s_flrtb
	871
	872	npt = numpnts(q_mlrtb)
	873	DeletePoints npt-nEnd,nEnd, q_mlrtb,i_mlrtb,s_mlrtb
	874
	875	npt = numpnts(q_b)
	876	DeletePoints npt-nEnd,nEnd, q_b,i_b,s_b
	877
[1023]	878	endif
	879
	880	// TODO - This is the identical set of waves as for the case of binType = 1.
	881	// they have the same names, but are averaged differently since it's slit mode.
	882	// I have separated this, since in practice the TB panels are probably best to ignore
	883	// and NOT include in the averaging since the Qy range is so limited.
	884	if(binType == 4)
	885	Wave/Z q_fb = qBin_qxqy_FB
	886	Wave/Z q_ft = qBin_qxqy_FT
	887	Wave/Z q_fl = qBin_qxqy_FL
	888	Wave/Z q_fr = qBin_qxqy_FR
	889	Wave/Z q_mb = qBin_qxqy_MB
	890	Wave/Z q_mt = qBin_qxqy_MT
	891	Wave/Z q_ml = qBin_qxqy_ML
	892	Wave/Z q_mr = qBin_qxqy_MR
	893	Wave/Z q_b = qBin_qxqy_B
[984]	894
[1023]	895	Wave/Z i_fb = iBin_qxqy_FB
	896	Wave/Z i_ft = iBin_qxqy_FT
	897	Wave/Z i_fl = iBin_qxqy_FL
	898	Wave/Z i_fr = iBin_qxqy_FR
	899	Wave/Z i_mb = iBin_qxqy_MB
	900	Wave/Z i_mt = iBin_qxqy_MT
	901	Wave/Z i_ml = iBin_qxqy_ML
	902	Wave/Z i_mr = iBin_qxqy_MR
	903	Wave/Z i_b = iBin_qxqy_B
	904
	905	Wave/Z s_fb = eBin_qxqy_FB
	906	Wave/Z s_ft = eBin_qxqy_FT
	907	Wave/Z s_fl = eBin_qxqy_FL
	908	Wave/Z s_fr = eBin_qxqy_FR
	909	Wave/Z s_mb = eBin_qxqy_MB
	910	Wave/Z s_mt = eBin_qxqy_MT
	911	Wave/Z s_ml = eBin_qxqy_ML
	912	Wave/Z s_mr = eBin_qxqy_MR
	913	Wave/Z s_b = eBin_qxqy_B
[1036]	914
	915	DeletePoints 0,nBeg, q_fb,q_ft,q_fl,q_fr,q_mb,q_mt,q_ml,q_mr,q_b
	916	DeletePoints 0,nBeg, i_fb,i_ft,i_fl,i_fr,i_mb,i_mt,i_ml,i_mr,i_b
	917	DeletePoints 0,nBeg, s_fb,s_ft,s_fl,s_fr,s_mb,s_mt,s_ml,s_mr,s_b
	918	//since each set may have a different number of points
	919	npt = numpnts(q_fb)
	920	DeletePoints npt-nEnd,nEnd, q_fb,i_fb,s_fb
[984]	921
[1036]	922	npt = numpnts(q_ft)
	923	DeletePoints npt-nEnd,nEnd, q_ft,i_ft,s_ft
[1023]	924
[1036]	925	npt = numpnts(q_fl)
	926	DeletePoints npt-nEnd,nEnd, q_fl,i_fl,s_fl
[1023]	927
[1036]	928	npt = numpnts(q_fr)
	929	DeletePoints npt-nEnd,nEnd, q_fr,i_fr,s_fr
[1023]	930
[1036]	931	npt = numpnts(q_mb)
	932	DeletePoints npt-nEnd,nEnd, q_mb,i_mb,s_mb
[984]	933
[1036]	934	npt = numpnts(q_mt)
	935	DeletePoints npt-nEnd,nEnd, q_mt,i_mt,s_mt
[984]	936
[1036]	937	npt = numpnts(q_ml)
	938	DeletePoints npt-nEnd,nEnd, q_ml,i_ml,s_ml
[984]	939
[1036]	940	npt = numpnts(q_mr)
	941	DeletePoints npt-nEnd,nEnd, q_mr,i_mr,s_mr
[984]	942
[1036]	943	npt = numpnts(q_b)
	944	DeletePoints npt-nEnd,nEnd, q_b,i_b,s_b
	945
[984]	946	endif
	947
	948	SetDataFolder root:
	949	return(0)
	950	end
	951
	952
[999]	953
[1064]	954	/
[1042]	955	Proc V_Load_Data_ITX()
[1038]	956	V_Load_itx("","",0,0)
	957	end
	958
[1037]	959	// TODO
	960	// -- fill in
	961	// -- link somewhere?
	962	//
	963	// a function to load in the individual I(q) sets which were written out to a single
	964	// file, in itx format.
	965	//
	966	// The data, like other 1D data sets, is to be loaded to its own folder under root
	967	//
	968	// Then, the data sets can be plotted as VSANS data sets, depending on which data extensions are present.
	969	// (and color coded)
	970	// (and used for setting the trimming)
	971	// (and...)
	972	//
	973	//
	974	// see A_LoadOneDDataToName(fileStr,outStr,doPlot,forceOverwrite)
	975	//
	976	Function V_Load_itx(fileStr,outStr,doPlot,forceOverwrite)
	977	String fileStr, outstr
	978	Variable doPlot,forceOverwrite
	979
	980	SetDataFolder root: //build sub-folders for each data set under root
	981
	982	// if no fileStr passed in, display dialog now
	983	if (cmpStr(fileStr,"") == 0)
	984	fileStr = DoOpenFileDialog("Select a data file to load")
	985	if (cmpstr(fileStr,"") == 0)
	986	String/G root:Packages:NIST:gLastFileName = ""
	987	return(0) //get out if no file selected
	988	endif
	989	endif
	990
	991	//Load the waves, using default waveX names
	992	//if no path or file is specified for LoadWave, the default Mac open dialog will appear
	993	LoadWave/O/T fileStr
	994	// LoadWave/G/D/A/Q fileStr
	995	String fileNamePath = S_Path+S_fileName
	996	// String basestr = ParseFilePath(3,ParseFilePath(5,fileNamePath,":",0,0),":",0,0)
	997
	998	String basestr
	999	if (!cmpstr(outstr, "")) //Outstr = "", cmpstr returns 0
	1000	// enforce a short enough name here to keep Igor objects < 31 chars
	1001	baseStr = ShortFileNameString(CleanupName(S_fileName,0))
	1002	baseStr = CleanupName(baseStr,0) //in case the user added odd characters
	1003	//baseStr = CleanupName(S_fileName,0)
	1004	else
	1005	baseStr = outstr //for output, hopefully correct length as passed in
	1006	endif
	1007
	1008	// print "basestr :"+basestr
	1009	String fileName = ParseFilePath(0,ParseFilePath(5,filestr,":",0,0),":",1,0)
	1010	// print "filename :"+filename
	1011
	1012	Variable ii,num=ItemsinList(S_waveNames)
	1013
	1014	if(DataFolderExists("root:"+baseStr))
	1015	if (!forceOverwrite)
	1016	DoAlert 1,"The file "+S_filename+" has already been loaded. Do you want to load the new data file, overwriting the data in memory?"
	1017	if(V_flag==2) //user selected No, don't load the data
	1018	SetDataFolder root:
	1019	for(ii=0;ii<num;ii+=1)
	1020	KillWaves $(StringFromList(ii, S_waveNames)) // kill the waves that were loaded
	1021	endfor
	1022	if(DataFolderExists("root:Packages:NIST"))
	1023	String/G root:Packages:NIST:gLastFileName = filename
	1024	endif
	1025	return(0) //quits the macro
	1026	endif
	1027	endif
	1028	SetDataFolder $("root:"+baseStr)
	1029	else
	1030	NewDataFolder/S $("root:"+baseStr)
	1031	endif
	1032
	1033	// ////overwrite the existing data, if it exists
	1034
	1035	// a semicolon-delimited list of wave names loaded
	1036	//S_waveNames
	1037	for(ii=0;ii<num;ii+=1)
	1038	Duplicate/O $("root:"+StringFromList(ii, S_waveNames)), $(StringFromList(ii, S_waveNames))
	1039	endfor
	1040
	1041
	1042	// clean up
	1043	SetDataFolder root:
	1044
	1045	for(ii=0;ii<num;ii+=1)
	1046	KillWaves $(StringFromList(ii, S_waveNames)) // kill the waves that were loaded
	1047	endfor
	1048	// Duplicate/O $("root:"+n0), $w0
	1049	// Duplicate/O $("root:"+n1), $w1
	1050	// Duplicate/O $("root:"+n2), $w2
	1051
	1052	// no resolution matrix to make
	1053
	1054
	1055	return(0)
[1038]	1056	End
	1057
	1058
	1059	// given strings of the number of points to remove, loop over the detectors
[1039]	1060	//
	1061	// TODO
	1062	// -- currently uses global strings or default strings
	1063	// -- if proper strings (non-null) are passed in, they are used, otherwise global, then default
[1038]	1064	Function V_Trim1DDataStr(folderStr,binType,nBegStr,nEndStr)
	1065	String folderStr
	1066	Variable binType
	1067	String nBegStr,nEndStr
	1068
	1069	String detListStr
	1070	if(binType == 1)
	1071	detListStr = ksBinType1
	1072	endif
	1073	if(binType == 2)
	1074	detListStr = ksBinType2
	1075	endif
	1076	if(binType == 3)
	1077	detListStr = ksBinType3
	1078	endif
	1079	if(binType == 4)
	1080	detListStr = ksBinType4
	1081	endif
	1082	if(strlen(detListStr)==0)
	1083	return(0)
	1084	endif
[1039]	1085
[1038]	1086
[1039]	1087	//use global, then default values if null string passed in
[1038]	1088	if(strlen(nBegStr)==0)
[1039]	1089	SVAR/Z gBegPtsStr=root:Packages:NIST:VSANS:Globals:Protocols:gBegPtsStr
	1090	SVAR/Z gEndPtsStr=root:Packages:NIST:VSANS:Globals:Protocols:gEndPtsStr
	1091
	1092	if(!SVAR_exists(gBegPtsStr) \|\| !SVAR_exists(gEndPtsStr) \|\| strlen(gBegPtsStr)==0 \|\| strlen(gEndPtsStr)==0)
	1093	nBegStr = ksBinTrimBegDefault
	1094	nEndStr = ksBinTrimEndDefault
	1095	else
	1096	nBegStr = gBegPtsStr
	1097	nEndStr = gEndPtsStr
	1098	endif
[1038]	1099	endif
	1100
	1101	Variable num, ii,nBeg,nEnd
	1102	String item,detstr
	1103
	1104	num = ItemsInList(detListStr)
	1105	for(ii=0;ii<num;ii+=1)
	1106	detStr = StringFromList(ii, detListStr)
	1107	nBeg = NumberByKey(detStr, nBegStr,"=",";")
	1108	nEnd = NumberByKey(detStr, nEndStr,"=",";")
	1109	V_TrimOneSet(folderStr,detStr,nBeg,nEnd)
	1110	endfor
	1111
	1112	return(0)
	1113	End
	1114
	1115	// TODO
[1064]	1116	// x- make this resolution-aware
[1038]	1117	//
	1118	Function V_TrimOneSet(folderStr,detStr,nBeg,nEnd)
	1119	String folderStr,detStr
	1120	Variable nBeg,nEnd
	1121
	1122	SetDataFolder $("root:Packages:NIST:VSANS:"+folderStr)
	1123
	1124	Printf "%d points removed from beginning, %d points from the end of %s \r",nbeg,nend,detStr
	1125
	1126	// TODO
	1127	// for each binType block:
	1128	// --declare the waves
	1129	// --make a copy of the waves??
	1130	// //--Break out resolution wave into separate waves
	1131	// --delete the beginning points from everything
	1132	// --trim off the last nEnd points from everything
	1133	// --DeletePoints num-nEnd,nEnd, qw,iw,sw
	1134	// // --delete all points where the shadow is < 0.98
	1135	////--Put resolution contents back???
	1136
	1137	Wave/Z qw = $("qBin_qxqy_"+detStr)
	1138	Wave/Z iw = $("iBin_qxqy_"+detStr)
	1139	Wave/Z ew = $("eBin_qxqy_"+detStr)
[1064]	1140	// resolution waves
	1141	Wave/Z sigQ = $("sigmaQ_qxqy_"+detStr)
	1142	Wave/Z qBar = $("qBar_qxqy_"+detStr)
	1143	Wave/Z fSubS = $("fSubS_qxqy_"+detStr)
[1038]	1144
[1064]	1145	DeletePoints 0,nBeg, qw,iw,ew,sigQ,qBar,fSubS
[1038]	1146
	1147	Variable npt
	1148	npt = numpnts(qw)
[1064]	1149	DeletePoints npt-nEnd,nEnd, qw,iw,ew,sigQ,qBar,fSubS
[1038]	1150
	1151	return(0)
[1064]	1152	End
	1153
	1154
	1155	//
	1156	// returns 1 if the val is non-negative, other value
	1157	// indicates that the resoution data is USANS data.
	1158	//
	1159	// TODO:
	1160	// -- this DUPLICATES a same-named SANS procedure, so there could be a clash at some point
	1161	// -- bigger issue - I'll need a better way to identify and load the different resolution
	1162	// conditions with VSANS
	1163	//
	1164	//
	1165	Function isSANSResolution(val)
	1166	Variable val
	1167
	1168	if(val >= 0)
	1169	return(1)
	1170	else
	1171	return(0)
	1172	endif
	1173	End
	1174
	1175

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

Download in other formats:

Original Format