def display_event(event):
"""an extremely inefficient display. It creates new instances of
CameraDisplay for every event and every camera, and also new axes
for each event. It's hacked, but it works
"""
print("Displaying... please wait (this is an inefficient implementation)")
global fig
ntels = len(event.r0.tels_with_data)
fig.clear()
plt.suptitle("EVENT {}".format(event.r0.event_id))
disps = []
for ii, tel_id in enumerate(event.r0.tels_with_data):
print("\t draw cam {}...".format(tel_id))
nn = int(ceil(sqrt(ntels)))
ax = plt.subplot(nn, nn, ii + 1)
geom = event.inst.subarray.tel[tel_id].camera
disp = CameraDisplay(geom, ax=ax, title="CT{0}".format(tel_id))
disp.pixels.set_antialiaseds(False)
disp.autoupdate = False
disp.cmap = random.choice(cmaps)
chan = 0
signals = event.r0.tel[tel_id].adc_sums[chan].astype(float)
signals -= signals.mean()
disp.image = signals
disp.set_limits_percent(95)
disp.add_colorbar()
disps.append(disp)
return disps
def display_event(event, geoms):
"""an extremely inefficient display. It creates new instances of
CameraDisplay for every event and every camera, and also new axes
for each event. It's hacked, but it works
"""
print("Displaying... please wait (this is an inefficient implementation)")
global fig
ntels = len(event.r0.tels_with_data)
fig.clear()
plt.suptitle("EVENT {}".format(event.r0.event_id))
disps = []
for ii, tel_id in enumerate(event.r0.tels_with_data):
print("\t draw cam {}...".format(tel_id))
nn = int(ceil(sqrt(ntels)))
ax = plt.subplot(nn, nn, ii + 1)
x, y = event.inst.pixel_pos[tel_id]
geom = geoms[tel_id]
disp = CameraDisplay(geom, ax=ax, title="CT{0}".format(tel_id))
disp.pixels.set_antialiaseds(False)
disp.autoupdate = False
disp.cmap = 'afmhot'
chan = 0
signals = event.r0.tel[tel_id].adc_sums[chan].astype(float)
signals -= signals.mean()
disp.image = signals
disp.set_limits_percent(95)
disp.add_colorbar()
disps.append(disp)
return disps
def plot_event(event, reco, pdf):
cams = [
event.inst.subarray.tels[i].camera for i in event.r0.tels_with_data
]
cams = [c for c in cams if c.cam_id in allowed_cameras]
n_tels = len(cams)
p = 1
params = {}
pointing_azimuth = {}
pointing_altitude = {}
for telescope_id, dl1 in event.dl1.tel.items():
camera = event.inst.subarray.tels[telescope_id].camera
if camera.cam_id not in allowed_cameras:
continue
nn = int(np.ceil(np.sqrt(n_tels)))
ax = plt.subplot(nn, nn, p)
p += 1
boundary_thresh, picture_thresh = cleaning_level[camera.cam_id]
mask = tailcuts_clean(camera,
dl1.image[0],
boundary_thresh=boundary_thresh,
picture_thresh=picture_thresh,
min_number_picture_neighbors=1)
#
if mask.sum() < 3: # only two pixel remaining. No luck anyways.
continue
h = hillas_parameters(
camera[mask],
dl1.image[0, mask],
)
disp = CameraDisplay(camera, ax=ax, title="CT{0}".format(telescope_id))
disp.pixels.set_antialiaseds(False)
disp.autoupdate = False
disp.add_colorbar()
# Show the camera image and overlay Hillas ellipse and clean pixels
disp.image = dl1.image[0]
disp.cmap = 'viridis'
disp.highlight_pixels(mask, color='white')
disp.overlay_moments(h, color='red', linewidth=5)
pointing_azimuth[
telescope_id] = event.mc.tel[telescope_id].azimuth_raw * u.rad
pointing_altitude[
telescope_id] = event.mc.tel[telescope_id].altitude_raw * u.rad
params[telescope_id] = h
return reco.predict(params, event.inst, pointing_altitude,
pointing_azimuth)