Graduation-Project/script/roi.py

28 lines
725 B
Python

import numpy as np
from matplotlib import pyplot as plt
or1 = [10, 11, 12, 13, 14]
or2 = [12.14, 13.35, 14.57, 15.79, 17]
energy = [600, 650, 700, 750, 800, 850, 900, 950]
markers = ["^", "+", "s", "D", "."]
data = np.loadtxt("test/eff.txt")
data = data.reshape(10, 8)
roi = np.loadtxt("test/roi.txt")
roi = roi.reshape(10, 8)
roi = 1 - roi / data
plt.figure(dpi=400)
for k in range(5):
plt.scatter(energy, roi[k], label="{:d}cm".format(or1[k]), marker=markers[k])
plt.xlabel("Energy [keV]", fontsize=15)
plt.ylabel(r"Energy Leakage Rate $\epsilon$", fontsize=15)
plt.legend(loc="lower left")
plt.tight_layout()
plt.subplots_adjust()
# plt.show()
plt.savefig("/home/fox/Documents/Thesis/figures/chap3/Fig7(a).png")