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'])