check_difc_diff.py
· 1.3 KiB · Python
原始文件
# import mantid algorithms, numpy and matplotlib
from mantid.simpleapi import *
import matplotlib.pyplot as plt
import numpy as np
LoadDiffCal(
InstrumentName="PG3",
Filename="/SNS/PG3/shared/CALIBRATION/2026-1_11A_CAL/2026A_PAC_d62140_2026-03-05.h5",
WorkspaceName="zyp"
)
LoadDiffCal(
InstrumentName="PG3",
Filename="/SNS/PG3/shared/CALIBRATION/2026-1_11A_CAL/2026A_PAC_d62140_2026-03-05_correct.h5",
WorkspaceName="zq"
)
all_difc_zyp = []
all_difc_zq = []
all_difc_diff = []
all_mask_zyp = []
all_mask_zq = []
all_mask_diff = []
for i in range(43120):
difc_zyp = mtd["zyp_cal"].row(i)["difc"]
difc_zq = mtd["zq_cal"].row(i)["difc"]
all_difc_zyp.append(difc_zyp)
all_difc_zq.append(difc_zq)
all_difc_diff.append(difc_zyp - difc_zq)
all_mask_zyp.append(int(mtd["zyp_mask"].readY(i)[0]))
all_mask_zq.append(int(mtd["zq_mask"].readY(i)[0]))
all_mask_diff.append(int(mtd["zyp_mask"].readY(i)[0]) - int(mtd["zq_mask"].readY(i)[0]))
with open("/SNS/PG3/shared/Yuanpeng/PG3_calib_test/PAC_difc_compare.dat", "w") as f:
for i in range(43120):
f.write(
"{0:10.2f}{1:10.2f}{2:10.2f}{3:10d}{4:10d}{5:10d}\n".format(
all_difc_zyp[i], all_difc_zq[i], all_difc_diff[i],
all_mask_zyp[i], all_mask_zq[i], all_mask_diff[i]
)
)
| 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 | ) |
| 45 |