#files = ['g77-61_blue_sum.fits']
print("Files in directory: ", files)
for filename in files:
print("Current file: ", filename)
#filename = "spSpec-2371-53762-261.fit"
#input_file = pd.read_csv(path + filename, header=None, sep="\s+")
input_file = pd.read_csv(path + filename)
#input_file = fits.open(path + filename)
### Create Spectrum object
spec = Spectrum(input_file, filename, fits=False)
### generate segments
#spec.generate_segments(15)
spec.generate_inflection_segments(sigma=25, cahk=True)
spec.assess_segment_variation()
spec.define_cont_points()
spec.set_segment_continuum()
spec.set_segment_midpoints()
spec.get_continuum_points()
#spec.remove_point(8)
#spec.remove_point([11,17])
####
#spec.remove_point([19,20])
#spec.add_continuum_point((4438.9, 706.02))
#spec.add_continuum_point((4547.34, 648.371))
