def lorentzian(x, height=1., center=0., width=1.):
	""" defined such that height is the height when x==x0 """
	halfWSquared = (width/2.)**2
	return (height * halfWSquared) / ((x - center)**2 + halfWSquared)