Skip to content

Commit

Permalink
[evt.modules.spms] also remove HPGe window masking
Browse files Browse the repository at this point in the history
  • Loading branch information
gipert committed May 10, 2024
1 parent 4355d57 commit 9185cfe
Showing 1 changed file with 8 additions and 22 deletions.
30 changes: 8 additions & 22 deletions src/pygama/evt/modules/spms.py
Original file line number Diff line number Diff line change
Expand Up @@ -331,26 +331,12 @@ def geds_coincidence_classifier(
energy summed over all channels, in p.e.). Derived from forced trigger
data.
"""
# mask for windowing data around the HPGe t0
geds_t0_mask = make_pulse_data_mask(
datainfo,
tcm,
table_names,
a_thr_pe=None,
t_loc_ns=geds_t0_ns,
dt_range_ns=(-1_000, 5_000),
t_loc_default_ns=48_000,
).view_as("ak")

# remove pulses outside the HPGe trigger window
spms_t0 = spms_t0.view_as("ak")[geds_t0_mask]
spms_amp = spms_amp.view_as("ak")[geds_t0_mask]

# get the (HPGe) trigger position
geds_t0_ns = geds_t0_ns.view_as("ak")

ts_data = larveto.l200_combined_test_stat(
spms_amp, spms_amp, geds_t0_ns, ts_bkg_prob, rc_density
return types.Array(
larveto.l200_combined_test_stat(
spms_t0.view_as("ak"),
spms_amp.view_as("ak"),
geds_t0_ns.view_as("ak"),
ts_bkg_prob,
rc_density,
)
)

return types.Array(ts_data)

0 comments on commit 9185cfe

Please sign in to comment.