check_difc_diff.py(檔案已創建)
| @@ -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 | + | ) | |
上一頁
下一頁