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()