minor fix in decay corr

AMYPAD · Jan 10, 2024 · 62db412 · 62db412
1 parent e46b35b
commit 62db412
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/amypet/align.py b/amypet/align.py
@@ -930,7 +930,7 @@ def align_break(
         lmbd = np.log(2) / nimpa.resources.riLUT[istp]['thalf']
 
         # > decay correction factor
-        dcycrr = 1 / np.exp(-lmbd * td)
+        dcycrr = np.exp(lmbd * td)
     else:
         dcycrr = 1.
 

diff --git a/amypet/align_brkdyn_ct.py b/amypet/align_brkdyn_ct.py
@@ -353,7 +353,8 @@ def align_break_petct(niidat, cts, Cnt, qcpth=None, refpetidx=None, use_stored=F
     outdct['faffines'] = faffines
     outdct['fnii_com'] = algn_frm[0]['fnii_com'] + algn_frm[1]['fnii_com']
 
-    np.save(fout, outdct)
+    if use_stored:
+        np.save(fout, outdct)
 
     return outdct