#pragma rtGlobals=1		// Use modern global access method.
//#include <AreaXYBetweenCursors>

Function twogauss(w,x) : FitFunc
	Wave w
	Variable x
       
	//CurveFitDialog/ These comments were created by the Curve Fitting dialog. Altering them will
	//CurveFitDialog/ make the function less convenient to work with in the Curve Fitting dialog.
	//CurveFitDialog/ Equation:
	//CurveFitDialog/ f(x) = W_0 + W_1*x + W_2*exp(-(x-W_3)^2/W_4^2)+W_5*exp(-(x-W_6)^2/W_7^2)
	//CurveFitDialog/ End of Equation
	//CurveFitDialog/ Independent Variables 1
	//CurveFitDialog/ x
	//CurveFitDialog/ Coefficients 8
	//CurveFitDialog/ w[0] = W_0
	//CurveFitDialog/ w[1] = W_1
	//CurveFitDialog/ w[2] = W_2
	//CurveFitDialog/ w[3] = W_3
	//CurveFitDialog/ w[4] = W_4
	//CurveFitDialog/ w[5] = W_5
	//CurveFitDialog/ w[6] = W_6
	//CurveFitDialog/ w[7] = W_7

	return w[0] + w[1]*x + w[2]*exp(-(x-w[3])^2/w[4]^2)+w[5]*exp(-(x-w[6])^2/w[7]^2)
End
		//return w[0] + w[1]*x + w[2]*x^2 + w[3]*x^3 + w[4]*exp(-(x-w[5])^2/w[6]^2)+w[7]*exp(-(x-w[8])^2/w[9]^2)
		//return W[0]+W[1]*exp(-(x-W[2])^2/W[3]^2)+W[4]+W[5]*exp(-(x-W[6])^2/W[7]^2)
//End
Macro load()
//string filename
Killwaves/A/Z
NewPath/O path "C:\Program Files\Andor Newton\data"
//filename=".asc" //specify file location of data
LoadWave/G/A/P=path //filename
Silent 1 ; PauseUpdate

Display/N = Spectrum wave1 vs wave0
ModifyGraph tick=2,mirror=1
Label left "Intensity (a.u.)"
Label bottom "Raman Shift (cm\\S-1\\M)"
showinfo 
 end macro
 //Macro basefitting (numspec)
 //variable numspec, i
 //i=1
 //do
 
 //while (i<numspec+1)
// Display/N = Spectrum1 wave1 vs wave0
//ModifyGraph tick=2,mirror=1
//Label left "Intensity (a.u.)"
//Label bottom "Raman Shift (cm\\S-1\\M)"
//showinfo 
 //end macro

Macro gaussianfit(C1_S, C1_E, nopix, numspec) 
variable nopix, numspec
variable C1_S, C1_E,base,amp,freq1,freq2,width
prompt C1_S, "input the start point of  C=C mode "
prompt C1_E, "input the end point of  C=C mode "
prompt nopix,"enter the dimentions"
prompt numspec,"enter the frames"
//prompt base,"enter the initial guess for base"
//prompt amp,"enter the initial guess for amplitude"
//prompt freq1,"enter theinitial guess for lower frequency"
//prompt freq2,"enter theinitial guess for higher freqency"
//prompt width,"enter the initial guess for width"
variable i

make/o/d/n= (nopix,nopix) freqFIT_1
make/o/d/n= (nopix,nopix) freqFIT_2
//make/o/d/n= (nopix,nopix) width  
make/o/d/n= (nopix,nopix) peakareaFIT_1
make/o/d/n= (nopix,nopix) peakareaFIT_2  
make/o/d/n= (nopix,nopix) widthFIT_1
make/o/d/n= (nopix,nopix) widthFIT_2 
//make/o/d/n= (numspec) wavec
//Make/D/O W_Coef={base,amp,freq1,width,base,amp,freq2,width}
//Change the coefs when you load data set, and the coefs will aodpt the real fit numbers as next initial guess
//Make/T/O TwoGaussConstraint = {"13<K4<20 ", "13<K7 < 20"}
//FuncFit TwoGauss /C=TwoGaussConstraint
i=1
do
      //appendtograph $("wave"+ num2str(i)) vs wave0
      //cursor A,$("wave"+ num2str(i)),C1_S
      //cursor B,$("wave"+ num2str(i)), C1_E
   //Smooth/E=2/S=2 45, $("wave"+ num2str(i))
   // appendtograph $("wave"+ num2str(i)) vs wave0
    //wavec= $("wave"+ num2str(i))
    
      FuncFit/NTHR=1 twogauss W_coef  $("wave"+ num2str(i))[C1_S,C1_E] /X=wave0 /D 
     //twogauss(x) = w[0] + w[1]*x + w[2]*exp(-(x-w[3])^2/w[4]^2)+w[5]*exp(-(x-w[6])^2/w[7]^2)
     freqFIT_1(i-1)=W_coef[3] 
     freqFIT_2(i-1)=W_coef[6]
     //FWHM=2*sqrt (ln2)*width
     widthFIT_1(i-1)=abs(W_coef[4]*1.665)
     widthFIT_2(i-1)=abs(W_coef[7]*1.665)
     //width(i-1)=W[3]
     //peakareaFIT_1(i-1)=W_coef[1]
     //peakareaFIT_2(i-1)=W_coef[5]
     peakareaFIT_1(i-1)=abs(W_coef[2]*W_coef[4]*sqrt(pi))
     peakareaFIT_2(i-1)=abs(W_coef[5]*W_coef[7]*sqrt(pi))
     //removefromgraph fit_$("wave"+ num2str(i))
     //removefromgraph $("wave"+ num2str(i))
         i += 1
		while (i <=numspec)
Newimage /N =fre1 freqFIT_1
ModifyImage freqFIT_1 ctab= {*,*,Rainbow,0}
Newimage /N =fre2 freqFIT_2
ModifyImage freqFIT_2 ctab= {*,*,Rainbow,0}
Newimage /N =inten1 peakareaFIT_1
ModifyImage peakareaFIT_1 ctab= {*,*,Rainbow,0}
Newimage /N =inten2 peakareaFIT_2
ModifyImage peakareaFIT_2 ctab= {*,*,Rainbow,0}
Newimage /N =width1 widthFIT_1
ModifyImage widthFIT_1 ctab= {*,*,Rainbow,0}
Newimage /N =width2 widthFIT_2
ModifyImage widthFIT_2 ctab= {*,*,Rainbow,0}
	 end