Änderungen von check_difc_diff.py

1 file changed, 44 insertions

check_difc_diff.py(Datei erstellt)

		@@ -0,0 +1,44 @@
1	+	# import mantid algorithms, numpy and matplotlib
2	+	from mantid.simpleapi import *
3	+	import matplotlib.pyplot as plt
4	+	import numpy as np
5	+
6	+	LoadDiffCal(
7	+	InstrumentName="PG3",
8	+	Filename="/SNS/PG3/shared/CALIBRATION/2026-1_11A_CAL/2026A_PAC_d62140_2026-03-05.h5",
9	+	WorkspaceName="zyp"
10	+	)
11	+	LoadDiffCal(
12	+	InstrumentName="PG3",
13	+	Filename="/SNS/PG3/shared/CALIBRATION/2026-1_11A_CAL/2026A_PAC_d62140_2026-03-05_correct.h5",
14	+	WorkspaceName="zq"
15	+	)
16	+
17	+	all_difc_zyp = []
18	+	all_difc_zq = []
19	+	all_difc_diff = []
20	+
21	+	all_mask_zyp = []
22	+	all_mask_zq = []
23	+	all_mask_diff = []
24	+
25	+	for i in range(43120):
26	+	difc_zyp = mtd["zyp_cal"].row(i)["difc"]
27	+	difc_zq = mtd["zq_cal"].row(i)["difc"]
28	+
29	+	all_difc_zyp.append(difc_zyp)
30	+	all_difc_zq.append(difc_zq)
31	+	all_difc_diff.append(difc_zyp - difc_zq)
32	+
33	+	all_mask_zyp.append(int(mtd["zyp_mask"].readY(i)[0]))
34	+	all_mask_zq.append(int(mtd["zq_mask"].readY(i)[0]))
35	+	all_mask_diff.append(int(mtd["zyp_mask"].readY(i)[0]) - int(mtd["zq_mask"].readY(i)[0]))
36	+
37	+	with open("/SNS/PG3/shared/Yuanpeng/PG3_calib_test/PAC_difc_compare.dat", "w") as f:
38	+	for i in range(43120):
39	+	f.write(
40	+	"{0:10.2f}{1:10.2f}{2:10.2f}{3:10d}{4:10d}{5:10d}\n".format(
41	+	all_difc_zyp[i], all_difc_zq[i], all_difc_diff[i],
42	+	all_mask_zyp[i], all_mask_zq[i], all_mask_diff[i]
43	+	)
44	+	)

Neuer Älter