From 918bd1e3c7349a50c6ea99e9810c1b0be5367f05 Mon Sep 17 00:00:00 2001 From: josh0-jrg Date: Thu, 7 Nov 2024 06:34:21 -0800 Subject: [PATCH] Fix: expand FZ to it's isolated value in dt TODO: make FV depend on time-stamp --- flamedisx/lz/WS2024_cuts_and_acceptances.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/flamedisx/lz/WS2024_cuts_and_acceptances.py b/flamedisx/lz/WS2024_cuts_and_acceptances.py index e672856c..131e27cd 100644 --- a/flamedisx/lz/WS2024_cuts_and_acceptances.py +++ b/flamedisx/lz/WS2024_cuts_and_acceptances.py @@ -109,7 +109,7 @@ def WS2024_fiducial_volume_cut(x,y,dt): ~((phi <= phi_slices[i]) & (phi > phi_slices[i + 1])) for i in range(n_phi_slices)] # Calculate dR_phi by replacing relevant points in loops over phi slices, then return - dR_phi = np.zeros(len(R)) + dR_phi = np.zeros_like(x) for i, p in enumerate(phi_cuts): dR_phi[p] = R[p] - phi_walls[i](dt[p]) @@ -121,7 +121,7 @@ def WS2024_fiducial_volume_cut(x,y,dt): mask = ((dR_phi < (contour + expansion)) & expandable) | ((dR_phi < contour) & ~expandable) #cut the drift time - dt_cut = (dt > 71) & (dt < 1030) + dt_cut = (dt > 71) & (dt < 1034) return dt_cut&mask&(dR_phi<=0)