#pragma rtGlobals=1		// Use modern global access method.

Function SwWave(Foldername,PosWavePrefix,NegWavePrefix,Temperature,posthres1,posthres2,negthres1,negthres2) //thresholds need to be indexed as well
	Wave Temperature
	Variable Foldername,PosWavePrefix,NegWavePrefix,posthres1,posthres2,negthres1,negthres2
	Variable NumT,h,i,j,k,l, NumPointsPosSwWave,NumPointsNegSwWave //might have to define PosFilename PosWavename etc. I don't think so, but maybe
	SVAR PosFilename, NegFilename, PosLowWave, PosMediumWave, PosLargeWave, NegLowWave, NegMediumWave, NegLargeWave
	Wave PositiveWave,NegativeWave
	h=0
	For (h=4; h<NumT;h+=1) //will go through each temperature run
		Sprintf PosFilename "%s\\%s%-5.3f", Foldername, PosWavePrefix, Temperature[h] //need to sprintf Temperature too, currently this will only work for one experiment
		LoadWave /H PosFilename
		Sprintf PosWavename "%s%-5.3f", PosWavePrefix, Temperature[h]
		Sprintf NegFilename "%s\\%s%-5.3f", Foldername, NegWavePrefix, Temperature[h] //need to sprintf Temperature too
		LoadWave /H NegFilename
		Sprintf NegWavename "%s%-5.3f", NegWavePrefix, Temperature[h]
		//Now the specific file should be loaded and reference-able (Pos/NegWavename)
		
		
		NumT=numpnts(Temperature) //should give 32, to access each temperature sweep
		i=0 //counting variables
		j=0
		k=0
		l=0
		make /O/N=0 PosLowSwitchWave, PosMedSwitchWave, PosLargeSwitchWave //make the three different waves to have the positive histogram sorted into
		make /O/N=0 NegLowSwitchWave, NegMedSwitchWave, NegLargeSwitchWave //make the three different waves to have the negative histogram sorted into
		PositiveWave=$PosWavename
		NegativeWave=$NegWavename
		NumPointsPosSwWave=numpnts(PositiveWave) //should give NumPointsPosSwWave=1000
		NumPointsNegSwWave=numpnts(NegativeWave)//should give NumPointsNegSwWave=1000
		do 
			if ($PosWavename[i]<posthres1[h]) //Low Current Waves
				InsertPoints j,1,$PosLowWave
				$PosLowWave[j]=$PosWavename[i]
				j+=1
			elseif ($PosWavename[i]>posthres2[h]) //Large Current Waves
				InsertPoints k,1,$PosLargeWave
				$PosLargeWave[k]=$PosWavename[i]
				k+=1
			else
				InsertPoints l,1,$PosMediumWave //Medium Current Waves
				$PosMediumWave[l]=$PosWavename[i]
				l+=1
			endif
			i+=1
		while (i<NumPointsPosSwWave)
		Sprintf PosLowWave "%s\\%-5.3f", PosLowSwitchWave, Temperature[h]
		Save $PosLowWave
		Sprintf PosMediumWave "%s\\%-5.3f", PosMedSwitchWave, Temperature[h]
		Save $PosMediumWave
		Sprintf PosLargeWave "%s\\%-5.3f", PosLargeSwitchWave, Temperature[h]
		Save $PosLargeWave
		end
		do
			if ($NegWavename[i]<negthres1[h]) //Low Current Waves
				InsertPoints j,1,$NegLowWave
				$NegLowWave[j]=$NegWavename[i]
				j+=1
			elseif ($NegWavename[i]>negthres2[h]) //Large Current Waves
				InsertPoints k,1,$NegLargeWave
				$NegLargeWave[k]=$NegWavename[i]
				k+=1
			else
				InsertPoints l,1,$NegMediumWave //Medium Current Waves
				$NegMediumWave[l]=$NegWavename[i]
				l+=1
			endif
			i+=1
		while (i<NumPointsNegSwWave)
		Sprintf NegLowWave "%s\\%-5.3f", NegLowSwitchWave, Temperature[h]
		Save $NegLowWave
		Sprintf NegMediumWave "%s\\%-5.3f", NegMedSwitchWave, Temperature[h]
		Save $NegMediumWave
		Sprintf NegLargeWave "%s\\%-5.3f", NegLargeSwitchWave, Temperature[h]
		Save $NegLargeWave
		end
	endfor
end