function [data] = loadSimulations(directory, decay_bc, detuning, dephase_bc)
	dataFiles = dir([directory, decay_bc, '*_*_', detuning, '*_', dephase_bc, '*.dat']);
	widths = [];
	minAbsorptions = [];
	contrasts = [];
	drives = [];
	deltas = [];
	for i = 1:length(dataFiles)
		dataFile = [directory, dataFiles(i).name];
		data = load(dataFile);
		detunings = data(:,1);
		absorptions = data(:,2);
		[width, absorption, contrast_left, contrast_right] = getPeak(detunings, absorptions);
		contrast = [contrast_left, contrast_right];
		widths = [widths; width];
		minAbsorptions = [minAbsorptions; absorption];
		contrasts = [contrasts; contrast];
		sections = strfind(dataFile, '_');
		drive = dataFile(sections(1)+1:sections(2)-1);
		drives = [drives; str2num(drive)];
		delta = dataFile(sections(2)+1:sections(3)-1);
		deltas = [deltas; str2num(delta)];
	end
	data = [drives, widths, minAbsorptions, contrasts];
	data = sortrows(data, 1);
	drives = data(:,1); widths = data(:,2); minAbsorptions = data(:,3); contrasts = data(:,4:5);
	loglog(drives, widths);
	hold all
	loglog(drives, minAbsorptions.*10^5);
end