def test_global_peak_integration(camera_waveforms):
waveforms, _ = camera_waveforms
integrator = GlobalPeakIntegrator()
charge, _, _ = integrator.extract_charge(waveforms)
assert_allclose(charge[0][0], 232.559, rtol=1e-3)
assert_allclose(charge[1][0], 425.406, rtol=1e-3)
def test_global_peak_integration(example_event):
telid = list(example_event.r0.tel)[0]
data = example_event.r0.tel[telid].waveform
nsamples = data.shape[2]
ped = example_event.mc.tel[telid].pedestal
data_ped = data - np.atleast_3d(ped / nsamples)
data_ped = np.array([data_ped[0], data_ped[0]]) # Test LG functionality
integrator = GlobalPeakIntegrator()
integration, peakpos, window = integrator.extract_charge(data_ped)
def test_global_peak_integration(example_event):
telid = 11
data = example_event.r0.tel[telid].waveform
nsamples = data.shape[2]
ped = example_event.mc.tel[telid].pedestal
data_ped = data - np.atleast_3d(ped / nsamples)
data_ped = np.array([data_ped[0], data_ped[0]]) # Test LG functionality
integrator = GlobalPeakIntegrator()
integration, peakpos, window = integrator.extract_charge(data_ped)
assert_almost_equal(integration[0][0], 58, 0)
assert_almost_equal(integration[1][0], 58, 0)
assert peakpos[0][0] == 14
assert peakpos[1][0] == 14
def test_global_peak_integration():
telid = 11
event = get_test_event()
data = event.r0.tel[telid].adc_samples
nsamples = data.shape[2]
ped = event.mc.tel[telid].pedestal
data_ped = data - np.atleast_3d(ped/nsamples)
data_ped = np.array([data_ped[0], data_ped[0]]) # Test LG functionality
integrator = GlobalPeakIntegrator(None, None)
integration, peakpos, window = integrator.extract_charge(data_ped)
assert_almost_equal(integration[0][0], 58, 0)
assert_almost_equal(integration[1][0], 58, 0)
assert peakpos[0][0] == 14
assert peakpos[1][0] == 14