This code was developed in Summer 2011 by 2 summer students (Philip Versteeg, UU, Redmer Bertens, UU) under the supervision of Michele Floris. It can serve as a starting point for AOD analysis.