def test_nb_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 geom = example_event.inst.subarray.tel[telid].camera nei = geom.neighbor_matrix_where integrator = NeighbourPeakIntegrator() integrator.neighbours = nei integration, peakpos, window = integrator.extract_charge(data_ped)
def test_nb_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 geom = event.inst.subarray.tel[telid].camera nei = geom.neighbor_matrix_where integrator = NeighbourPeakIntegrator(None, None) integrator.neighbours = nei integration, peakpos, window = integrator.extract_charge(data_ped) assert_almost_equal(integration[0][0], -64, 0) assert_almost_equal(integration[1][0], -64, 0) assert peakpos[0][0] == 20 assert peakpos[1][0] == 20
def test_nb_peak_integration(test_event): telid = 11 event = deepcopy(test_event) data = event.r0.tel[telid].waveform 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 geom = event.inst.subarray.tel[telid].camera nei = geom.neighbor_matrix_where integrator = NeighbourPeakIntegrator() integrator.neighbours = nei integration, peakpos, window = integrator.extract_charge(data_ped) assert_almost_equal(integration[0][0], -64, 0) assert_almost_equal(integration[1][0], -64, 0) assert peakpos[0][0] == 20 assert peakpos[1][0] == 20
def setup(self): kwargs = dict(config=self.config, tool=self) self.dl0 = CameraDL0Reducer(**kwargs) self.dl1 = CameraDL1Calibrator(**kwargs) self.cal = CameraCalibrator(r1_product=self.calibrator) self.cross = CrossCorrelation() self.glob_peak = GlobalPeakIntegrator() self.local_peak = LocalPeakIntegrator() self.neighbour = NeighbourPeakIntegrator() self.aver = AverageWfPeakIntegrator()
def test_nb_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 geom = CameraGeometry.guess(*event.inst.pixel_pos[telid], event.inst.optical_foclen[telid]) nei = geom.neighbors integrator = NeighbourPeakIntegrator(None, None) integrator.neighbours = nei integration, peakpos, window = integrator.extract_charge(data_ped) assert_almost_equal(integration[0][0], -64, 0) assert_almost_equal(integration[1][0], -64, 0) assert peakpos[0][0] == 20 assert peakpos[1][0] == 20
def test_nb_peak_integration(camera_waveforms): waveforms, camera = camera_waveforms nei = camera.neighbor_matrix_where integrator = NeighbourPeakIntegrator() integrator.neighbours = nei charge, _, _ = integrator.extract_charge(waveforms) assert_allclose(charge[0][0], 94.671, rtol=1e-3) assert_allclose(charge[1][0], 426.887, rtol=1e-3) integrator.lwt = 4 integration, _, _ = integrator.extract_charge(waveforms) assert_allclose(integration[0][0], 220.418657, rtol=1e-3) assert_allclose(integration[1][0], 426.887, rtol=1e-3)