if __name__ == "__main__": # get file paths # data obtained from http://www.physics.emory.edu/faculty/weeks/data/arches/ abspath = os.path.abspath(__file__) dname = os.path.dirname(abspath) os.chdir(dname) mypath = "../data" oldfiles = [join(mypath, f) for f in listdir(mypath) if isfile(join(mypath, f)) and re.search(r'.*yale\.txt', f)] newfiles = [join(mypath, f) for f in listdir(mypath) if isfile(join(mypath, f)) and re.search(r'^sav.*\.txt', f)] # processing data old_interface = OldInterface(oldfiles) old_interface.cal_angles() concave_list = old_interface.pick_concave() plotter = Plotter.PlotCircle(concave_list) #Plotter.AngleDistribution_shaded(old_interface.combine_angle_with_width()) new_interface = NewInterface(newfiles) new_interface.cal_angles() angle_data = new_interface.combine_angle_with_width() Plotter.AngleDistribution(angle_data) Plotter.AngleDistribution_shaded(angle_data) Plotter.NDistribution(new_interface.combine_N_with_width()) Plotter.AngleDistributionVSN(new_interface.combine_angle_with_N()) Plotter.AngleDistributionVSN_shaded(new_interface.combine_angle_with_N()) Plotter.AngleGravityN_shaded(new_interface.combine_angle_with_N_and_g())