def make_plot(): fluence_set = [(0, 25), (25, 35), (35, 45), (45, 75), (75, 300)] spectra = abs_get_processed_data.get_damage_specs() f, axs = plt.subplots(1, 2, sharex=True, sharey='row', figsize=(5.5, 3)) for ind, spec in enumerate(spectra): fluences = fluence_set[ind] label = str(round(int( FLUENCE_CONVERSION * fluences[0]))) + ' to ' + str( round(int(FLUENCE_CONVERSION * fluences[1]))) + ' mJ/cm$^2$' axs[0].plot(spec['minus']['phot'], spec['minus']['spec'], label=label) axs[1].plot(spec['plus']['phot'], spec['plus']['spec']) plt.gcf().legend(loc=7, title='Maximum Previous\nFluence', frameon=True) _format_damage(axs) fig_path = os.path.join(os.path.dirname(__file__), 'figure_s3.eps') plt.savefig(fig_path, dpi=600) fig_path = os.path.join(os.path.dirname(__file__), 'figure_s3.jpg') plt.savefig(fig_path, dpi=600)
def get_source_data(): specs = abs_get_processed_data.get_damage_specs() source_data = { '0 to 11 mJ/cm^2 Photon Energy': specs[0]['minus']['phot'], '0 to 11 mJ/cm^2 Plus': specs[0]['plus']['spec'], '0 to 11 mJ/cm^2 Minus': specs[0]['minus']['spec'], '11 to 15 mJ/cm^2 Photon Energy': specs[1]['minus']['phot'], '11 to 15 mJ/cm^2 Plus': specs[1]['plus']['spec'], '11 to 15 mJ/cm^2 Minus': specs[1]['minus']['spec'], '15 to 19 mJ/cm^2 Photon Energy': specs[2]['minus']['phot'], '15 to 19 mJ/cm^2 Plus': specs[2]['plus']['spec'], '15 to 19 mJ/cm^2 Minus': specs[2]['minus']['spec'], '19 to 33 mJ/cm^2 Photon Energy': specs[3]['minus']['phot'], '19 to 33 mJ/cm^2 Plus': specs[3]['plus']['spec'], '19 to 33 mJ/cm^2 Minus': specs[3]['minus']['spec'], '33 to 132 mJ/cm^2 Photon Energy': specs[4]['minus']['phot'], '33 to 132 mJ/cm^2 Plus': specs[4]['plus']['spec'], '33 to 132 mJ/cm^2 Minus': specs[4]['minus']['spec'] } return source_data
def sanity_plot(): specs = abs_get_processed_data.get_damage_specs() f, axs = plt.subplots(1, 2, sharex=True) for spec in specs: axs[0].plot(spec['minus']['phot'], spec['minus']['spec']) axs[1].plot(spec['minus']['phot'], spec['plus']['spec'])