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


Function SongStim_IClamp()
	
	DoOpenMultiFileDialog()
	
	NVAR nlist = root:numFilesSelected
	WAVE/T stimlist = root:stimlist
	
	variable inval = nlist
	
	variable basedur, nreps
	Prompt nreps, "number reps song list"
	Prompt basedur, "baseline time in sec"
	DoPrompt "birdsong trial", nreps, basedur
	
	variable rep
	for(rep=0;rep<nreps;rep+=1)
		permute(inval)
		WAVE outseq=root:outseq
		
		variable item, pick
		for(item=0;item<nlist;item+=1)	// Initialize variables;continue test
			variable  sampfreq, songsample, tdur, normval, scaley, sampdelta, basesample, datasample, delta, beginwav, endwav
			string pathfromopen, sampfreqst, tit
		
			pick=outseq[item]
			pathfromopen = stimlist[item]		//get the stim file path returned in outputPaths list of DoOpenMultiFileDialog() function using the randomized index generated by purmute(inval)			
	
			SndLoadWave/M=1/O anywav, (pathfromopen)  //load the .WAV file from the stim file path determined by outputPaths list of DoOpenMultiFileDialog()
			wavestats/Q anywav
			sampfreqst=stringbykey("RATE", S_info)  //gets the sample frequency of the original .WAV file from the info strings returned by SndLoadWave
			//sampfreq=str2num(sampfreqst)
			sampfreq=22050
			tit=S_fileName
		
			songsample = V_npnts
			basesample =(basedur)*(sampfreq)
		
			tdur=songsample/sampfreq
			delta=1/sampfreq
		
			Make/D/N=(songsample) datawav, playwav
			WAVE anywav=root:anywav
			playwav=anywav
			InsertPoints 0,(basesample), playwav      //insert x zeros to the beginning of the .WAV to be played so that there is a baseline measurement of the cell's activity
			InsertPoints 0,(basesample), datawav
			//I may want to figure out how to add points after as well

			setscale/P x, 0, (delta), datawav, playwav
			normval = max(V_max, abs(V_min))
			playwav /=normval
			scaley = 6
			playwav *= scaley
			
			DAQmx_Scan/DEV="NI6259"/TRIG={"/NI6259/ao/starttrigger"}/BKG WAVES= "datawav, 1"
			DAQmx_WaveformGen/DEV="NI6259"/NPRD=1 "playwav, 0;"
			//fDAQmx_WF_WaitUntilFinished("NI6259",-1)
			fDAQmx_ScanWait("NI6259")
			//fDAQmx_ScanStop("NI6259")
			//fDAQmx_WaveformStop("NI6259")
		
			string datawavst, basewavst
			datawavst= (tit)+(num2str(rep))
			
			rename datawav, $(datawavst)
			rename playwav, $(datawavst) + "play"
			
			KillWaves anywav
			//, playwav
		endfor												
	endfor
 End
 
 Function/S DoOpenMultiFileDialog()
	Variable refNum
	String message = "Select one or more files"
	String outputPaths
	String fileFilters = "Data Files (*.txt,*.dat,*.csv,*.wav):.txt,.dat,.csv,.wav;"
	fileFilters += "All Files:.*;"

	Open/D/R/MULT=1/F=fileFilters/M=message refNum
	outputPaths = S_fileName
	
	if (strlen(outputPaths) == 0)
		Print "Cancelled"
	else
		Variable/G numFilesSelected = ItemsInList(outputPaths, "\r")
		Make/T/N=(numFilesSelected)/O stimlist
		Variable i
		for(i=0; i<numFilesSelected; i+=1)
			String path = StringFromList(i, outputPaths, "\r")
			Printf "%d: %s\r", i, path	
			stimlist[i]=path
		endfor
	endif
	
	return outputPaths		// Will be empty if user canceled
End

Function/WAVE permute(inval)
	variable inval
	Make/O/N=(inval) outseq
	variable i, j, val, fail, test, emax=(inval)/2
	for (i=1; i<=inval; i+=1)
		fail=0
		do
		val=  floor(emax + enoise(emax))
		fail=0
			for(j=0; j<(i-1); j+=1)
				if (val==outseq[j]) 
					fail+=1
				endif
			print j, fail
			endfor	
		
		
		while(fail!=0)
	 		 	
	 	outseq[i-1]=val
	 print outseq
	endfor
End