def test_neighbor_peak_window_sum(camera_waveforms): waveforms, subarray = camera_waveforms extractor = NeighborPeakWindowSum(subarray=subarray) charge, pulse_time = extractor(waveforms, telid=1) assert_allclose(charge[0], 94.671, rtol=1e-3) assert_allclose(pulse_time[0], 54.116092, rtol=1e-3) extractor.lwt = 4 charge, pulse_time = extractor(waveforms, telid=1) assert_allclose(charge[0], 220.418657, rtol=1e-3) assert_allclose(pulse_time[0], 48.717848, rtol=1e-3)
def test_neighbor_peak_window_sum(camera_waveforms): waveforms, camera = camera_waveforms nei = camera.neighbor_matrix_where extractor = NeighborPeakWindowSum() extractor.neighbors = nei charge, pulse_time = extractor(waveforms) assert_allclose(charge[0], 94.671, rtol=1e-3) assert_allclose(pulse_time[0], 54.116092, rtol=1e-3) extractor.lwt = 4 charge, pulse_time = extractor(waveforms) assert_allclose(charge[0], 220.418657, rtol=1e-3) assert_allclose(pulse_time[0], 48.717848, rtol=1e-3)
def test_neighbor_peak_window_sum(camera_waveforms): waveforms, camera = camera_waveforms nei = camera.neighbor_matrix_where extractor = NeighborPeakWindowSum() extractor.neighbors = nei charge, pulse_time = extractor(waveforms) assert_allclose(charge[0][0], 94.671, rtol=1e-3) assert_allclose(charge[1][0], 426.887, rtol=1e-3) assert_allclose(pulse_time[0][0], 54.116092, rtol=1e-3) assert_allclose(pulse_time[1][0], 62.038344, rtol=1e-3) extractor.lwt = 4 charge, pulse_time = extractor(waveforms) assert_allclose(charge[0][0], 220.418657, rtol=1e-3) assert_allclose(charge[1][0], 426.887, rtol=1e-3) assert_allclose(pulse_time[0][0], 48.717848, rtol=1e-3) assert_allclose(pulse_time[1][0], 62.038344, rtol=1e-3)