aminoAcidMassMapFile = './data/integer_mass_table.txt' # set current directory os.chdir(curDir) # read input with open(inputFile) as f: spectrumLine = f.readline().strip() spectrum = [int(n) for n in spectrumLine.split(" ")] # load SeqUtil methods sys.path.append('code') import SeqUtil importlib.reload(SeqUtil) ### original work #### # load amion acid mass AAMassMap = SeqUtil.loadAminoAcidMass(aminoAcidMassMapFile) # generate sequence from spectrum conv = SeqUtil.spectralConvolution(spectrum) conv = sorted(conv) # output with open(outputFile, "w") as f: f.writelines(" ".join([str(m) for m in conv]))