//Bob Dougherty November 2, 2010
//Compute CSM of all wave files in a folder
//Also save the narrowband spectra

nMics = 24;
keepCSM = true;
keepSpectra = true;

dir1 = getDirectory("Choose Data Directory");
if(keepCSM){
	outputCSMDirectory = "["+getDirectory("Choose destination directory for CSMs")+"]";
	if(startsWith(outputCSMDirectory,"[")){
		outputCSMDirectory = substring(outputCSMDirectory,1,lengthOf(outputCSMDirectory));
	}
	if(endsWith(outputCSMDirectory,"]")){
		outputCSMDirectory = substring(outputCSMDirectory,0,lengthOf(outputCSMDirectory)-1);
	}
}else{
	outputCSMDirectory = "[]";
}
if(keepSpectra){
	outputSpectraDirectory = "["+getDirectory("Choose destination directory for spectra")+"]";
	if(startsWith(outputSpectraDirectory,"[")){
		outputSpectraDirectory = substring(outputSpectraDirectory,1,lengthOf(outputSpectraDirectory));
	}
	if(endsWith(outputSpectraDirectory,"]")){
		outputSpectraDirectory = substring(outputSpectraDirectory,0,lengthOf(outputSpectraDirectory)-1);
	}
}else{
	outputSpectraDirectory = "[]";
}
list = getFileList(dir1);
for(i=0; i<list.length; i++) {
	if(endsWith(list[i],".wav")){
		run("Audio Reader 3", "data_folder="+dir1+" data_file="+list[i]+" number=512 convert");
		titleData =  getTitle();
		run("CSM IQ FFT"," ");
		titleCSM = getTitle();
		print(outputCSMDirectory+titleCSM);
		if(keepCSM)saveAs("Tiff", outputCSMDirectory+titleCSM);
		if(keepSpectra){
			run("CSM IQ Stats", "first=1 last="+nMics+" leveldelta=0 smoothbins=0 show spectral percent=50 live");
			selectWindow("Results");
			saveAs("Results", outputSpectraDirectory+titleData+"_NBspectra.txt");
			close();
		}
		close();
		close();
	}
}