The source code for estimating the fly's orientation, in matlab (2011a)