#pragma rtGlobals=1		// Use modern global access method.
 
//	Example:
//		Make/O/N=5 xWave0=p, yWave0=gnoise(1)
//		Make/O/N=5 xWave1=p, yWave1=gnoise(1)
//		String xWaveList = "xWave0;xWave1;"
//		String yWaveList = "yWave0;yWave1;"
//		FitAndGraphAllXYPairs(xWaveList, yWaveList)
 
Function FitAndGraphXYPair(xWave, yWave)
	Wave xWave
	Wave yWave
 
	Display yWave vs xWave
	ModifyGraph mode=3,marker=19,rgb=(0,0,65535)		// Round blue markers
	DoUpdate
 
	CurveFit /TBOX=256 line yWave /X=xWave /D
 
	String textboxName = "CF_" + NameOfWave(yWave)
	String text
 

// Adds a legend and a tag for each trace to the specified graph.
// Example:
//	Make/O/N=10 wave0=p, wave1=p+1, wave2=p+2
//	Display wave0, wave1, wave2
//	AutoAddAnnotations("")
Function AutoAddAnnotations(graphName)
	String graphName			// "" for top graph
 
	Legend/C/N=legend0		// Add legend
 
	String traceList = TraceNameList(graphName, ";", 1)
 
	Variable index = 0
	do
		String traceName = StringFromList(index, traceList)
		if (strlen(traceName) == 0)
			break				// All done.
		endif
 
		Wave w = TraceNameToWaveRef(graphName, traceName)
 

Function/S GetAxisLabel(gname, axisname)
	String gname, axisname
 
	String grecreation = WinRecreation(gname, 0)
	String oneLine
	String labelStr = ""
	String part1, part2
 
	Variable i=0
	do
		oneLine = StringFromList(i, grecreation, "\r")
		if (stringmatch(oneLine, "\tLabel*"))
			String regexp = "^\\tLabel "+axisname+" \\\"(.*)\\\""
			SplitString/E=regexp oneLine, labelStr
			break;
		endif
		i += 1

//  The following function finds A, Tau, and B from an exponential decay by
//	(1) taking the FFT, and 
//	(2) using the complicated function for alpha found in Rev. Sci. Instrum. 79, 023108 (2008)
Function FFTFit(RDTWave, A, Tau, B)
	wave RDTWave
	Variable &A, &Tau, &B
 
	FFT/OUT=1/DEST=RDTWave_FFT RDTWave
	variable tm = pnt2x(RDTWave, (numpnts(RDTWave)-1)) - pnt2x(RDTWave, 0)
	variable omega = pnt2x(RDTWave_FFT,1)*2*pi
	variable alpha = ( 1/ deltax(RDTWave) )* ln(cos(2*pi/numpnts(RDTWave))+(real(RDTWave_FFT[1])/imag(RDTWave_FFT[1]))*sin(2*pi/numpnts(RDTWave) )  )