psi=ii * 20 * u.deg, ) image, sig, bg = mock.make_mock_shower_image( geom, model.pdf, intensity=50, nsb_level_pe=1000, ) clean = image.copy() clean[image <= 3.0 * image.mean()] = 0.0 hillas = hillas_parameters(geom.pix_x, geom.pix_y, clean) disp.image = image disp.add_colorbar(ax=axs[ii]) disp.set_limits_percent(95) disp.overlay_moments(hillas, linewidth=3, color='blue') if __name__ == '__main__': hexgeom = io.CameraGeometry.from_name("hess", 1) recgeom = io.make_rectangular_camera_geometry() draw_several_cams(recgeom) draw_several_cams(hexgeom) plt.tight_layout() plt.show()
model = mock.generate_2d_shower_model(centroid=(0.2 - ii * 0.1, -ii * 0.05), width=0.005 + 0.001 * ii, length=0.1 + 0.05 * ii, psi=ii * 20 * u.deg) image, sig, bg = mock.make_mock_shower_image(geom, model.pdf, intensity=50, nsb_level_pe=1000) clean = image.copy() clean[image <= 3.0 * image.mean()] = 0.0 hillas = hillas_parameters(geom.pix_x.value, geom.pix_y.value, clean) disp.image = image disp.set_limits_percent(95) disp.overlay_moments(hillas, linewidth=3, color='blue') if __name__ == '__main__': hexgeom = io.get_camera_geometry("hess", 1) recgeom = io.make_rectangular_camera_geometry() draw_several_cams(recgeom) draw_several_cams(hexgeom) plt.tight_layout() plt.show()