fix pmcx dettpsf overwrites input bug

fangq · Feb 23, 2024 · 90e2419 · 90e2419
1 parent e28fd4c
commit 90e2419
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/pmcx/pmcx/utils.py b/pmcx/pmcx/utils.py
@@ -182,6 +182,8 @@ def dettpsf(detp, detnum, prop, time):
     output:
         tpsf: caculated temporal point spread function curve of the specified detector
     """
+    detp = detp.copy()
+
     # select the photon data of the specified detector
     detp["ppath"] = detp["ppath"][detp["detid"] == detnum, :]
     detp["detid"] = detp["detid"][detp["detid"] == detnum]
@@ -370,7 +372,9 @@ def detphoton(detp, medianum, savedetflag, issaveref=None, srcnum=None):
             newdetp["detid"] = detp[0, :].astype(int).transpose()
             if np.any(newdetp["detid"] > 65535):
                 newdetp["srcid"] = np.right_shift(newdetp["detid"], 16, dtype=np.int32)
-                newdetp["detid"] = np.bitwise_and(newdetp["detid"], 0xff, dtype=np.int32)
+                newdetp["detid"] = np.bitwise_and(
+                    newdetp["detid"], 0xFF, dtype=np.int32
+                )
         c0 = 1
 
     length = medianum