Skip to content

Commit 516677d

Browse files
committed
fix compton error
1 parent afe3dc4 commit 516677d

File tree

1 file changed

+2
-9
lines changed

1 file changed

+2
-9
lines changed

src/pygama/pargen/survival_fractions.py

Lines changed: 2 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -550,15 +550,8 @@ def compton_sf(
550550
else:
551551
raise ValueError("mode not recognised")
552552

553-
ct_n = len(cut_param[~mask])
554-
ct_err = np.sqrt(len(cut_param[~mask]))
555-
surv_n = len(cut_param[mask])
556-
surv_err = np.sqrt(len(cut_param[mask]))
557-
558-
pc_n = ct_n + surv_n
559-
560-
sf = surv_n / pc_n
561-
err = 100 * sf * (1 - sf) * np.sqrt((ct_err / ct_n) ** 2 + (surv_err / surv_n) ** 2)
553+
sf = len(cut_param[mask]) / len(cut_param)
554+
err = 100 * np.sqrt(sf * (1 - sf)) / len(cut_param)
562555
sf *= 100
563556

564557
return {

0 commit comments

Comments
 (0)