function [width, minAbsorption, detunings, absorptions] = loadSimulation(dataFile)
	data = load(dataFile);
	detunings = data(:,1);
	absorptions = data(:,2);
	[width, minAbsorption, transmission, features] = getPeak(detunings, absorptions);
	figure
	hold all
	plot(detunings, absorptions)
	plot(features(:,1), features(:,2), 'r*')

end