From 5f673e3d6f42c2081eef3b5e2f0db91cf7c82af1 Mon Sep 17 00:00:00 2001 From: Jaafar Chakrani Date: Tue, 17 Feb 2026 08:29:31 -0800 Subject: [PATCH 1/2] Skip the periodic reset during the hit cycle --- larndsim/fee.py | 9 --------- 1 file changed, 9 deletions(-) diff --git a/larndsim/fee.py b/larndsim/fee.py index 095769bb..1fcd2ecb 100644 --- a/larndsim/fee.py +++ b/larndsim/fee.py @@ -603,15 +603,6 @@ def get_adc_values(pixels_signals, break q = 0 - - if apply_periodic_reset and ic % detector.PERIODIC_RESET_CYCLES == periodic_reset_phase: - q_sum = xoroshiro128p_normal_float32(rng_states, ip) * detector.RESET_NOISE_CHARGE * e - true_q = 0 - last_periodic_reset = ic - for itrk in range(current_fractions.shape[2]): - current_fractions[ip][iadc][itrk] = 0 - ic += 1 - continue if detector.BUFFER_RISETIME > 0: conv_start = max(last_reset, floor(ic - 10*detector.BUFFER_RISETIME/detector.TIME_SAMPLING)) From 9faffbc56d7acf7aa5d80b68dfb5a2c775798f17 Mon Sep 17 00:00:00 2001 From: Jaafar Chakrani Date: Tue, 17 Feb 2026 08:35:46 -0800 Subject: [PATCH 2/2] Skip the periodic reset during the hit cycle --- larndsim/fee.py | 18 +++++++++--------- 1 file changed, 9 insertions(+), 9 deletions(-) diff --git a/larndsim/fee.py b/larndsim/fee.py index 1fcd2ecb..045d6797 100644 --- a/larndsim/fee.py +++ b/larndsim/fee.py @@ -603,6 +603,15 @@ def get_adc_values(pixels_signals, break q = 0 + + if apply_periodic_reset and ic % detector.PERIODIC_RESET_CYCLES == periodic_reset_phase: + q_sum = xoroshiro128p_normal_float32(rng_states, ip) * detector.RESET_NOISE_CHARGE * e + true_q = 0 + last_periodic_reset = ic + for itrk in range(current_fractions.shape[2]): + current_fractions[ip][iadc][itrk] = 0 + ic += 1 + continue if detector.BUFFER_RISETIME > 0: conv_start = max(last_reset, floor(ic - 10*detector.BUFFER_RISETIME/detector.TIME_SAMPLING)) @@ -638,15 +647,6 @@ def get_adc_values(pixels_signals, while ic <= integrate_end: q = 0 - if apply_periodic_reset and ic % detector.PERIODIC_RESET_CYCLES == periodic_reset_phase: - q_sum = xoroshiro128p_normal_float32(rng_states, ip) * detector.RESET_NOISE_CHARGE * e - true_q = 0 - last_periodic_reset = ic - for itrk in range(current_fractions.shape[2]): - current_fractions[ip][iadc][itrk] = 0 - ic += 1 - continue - if detector.BUFFER_RISETIME > 0: conv_start = max(last_reset, floor(ic - 10*detector.BUFFER_RISETIME/detector.TIME_SAMPLING)) for jc in range(conv_start, min(ic+1, curre.shape[0])):