def calibrate_event(event, params, geom_dict=None):
"""
Generic calibrator for events. Calls the calibrator corresponding to the
source of the event, and stores the dl1 (calibrated_image) information into a
new event container.
Parameters
----------
event : container
A `ctapipe` event container
params : dict
REQUIRED:
params['integrator'] - Integration scheme
params['integration_window'] - Integration window size and shift of
integration window centre
(adapted such that window fits into readout).
OPTIONAL:
params['integration_clip_amp'] - Amplitude in p.e. above which the
signal is clipped.
params['integration_calib_scale'] : Identical to global variable
CALIB_SCALE in reconstruct.c in hessioxxx software package. 0.92 is
the default value (corresponds to HESS). The required value changes
between cameras (GCT = 1.05).
params['integration_sigamp'] - Amplitude in ADC counts above pedestal
at which a signal is considered as significant (separate for
high gain/low gain).
geom_dict : dict
Dict of pixel geometry for each telescope. Leave as None for automatic
calculation when it is required.
dict[(num_pixels, focal_length)] = `ctapipe.io.CameraGeometry`
Returns
-------
calibrated : container
A new `ctapipe` event container containing the dl1 information, and a
reference to all other information contained in the original event
container.
"""
# Obtain relevent calibrator
switch = {'hessio': partial(mc.calibrate_mc, event=event, params=params)}
try:
calibrator = switch[event.meta['source']]
except KeyError:
log.exception("no calibration created for data origin: '{}'".format(
event.meta['source']))
raise
# KPK: should not copy the event here! there is no reason to
# Copying is
# up to the user if they want to do it, not in the algorithms.
# calibrated = copy(event)
# params stored in metadata
event.dl1.meta.update(params)
# Fill dl1
event.dl1.reset()
for telid in event.dl0.tels_with_data:
nchan = event.inst.num_channels[telid]
npix = event.inst.num_pixels[telid]
event.dl1.tel[telid] = CalibratedCameraContainer()
# Get geometry
int_dict, inverted = integrator_dict()
geom = None
# Check if geom is even needed for integrator
if inverted[params['integrator']] in integrators_requiring_geom():
if geom_dict is not None and telid in geom_dict:
geom = geom_dict[telid]
else:
log.debug("[calib] Guessing camera geometry")
geom = CameraGeometry.guess(*event.inst.pixel_pos[telid],
event.inst.optical_foclen[telid])
log.debug("[calib] Camera geometry found")
if geom_dict is not None:
geom_dict[telid] = geom
pe, window, data_ped, peakpos = calibrator(telid=telid, geom=geom)
tel = event.dl1.tel[telid]
tel.calibrated_image = pe
tel.peakpos = peakpos
for chan in range(nchan):
tel.integration_window[chan] = window[chan]
tel.pedestal_subtracted_adc[chan] = data_ped[chan]
return event
def calibrate_event(event, params, geom_dict=None):
"""
Generic calibrator for events. Calls the calibrator corresponding to the
source of the event, and stores the dl1 (pe_charge) information into a
new event container.
Parameters
----------
event : container
A `ctapipe` event container
params : dict
REQUIRED:
params['integrator'] - Integration scheme
params['integration_window'] - Integration window size and shift of
integration window centre
(adapted such that window fits into readout).
OPTIONAL:
params['integration_clip_amp'] - Amplitude in p.e. above which the
signal is clipped.
params['integration_calib_scale'] : Identical to global variable
CALIB_SCALE in reconstruct.c in hessioxxx software package. 0.92 is
the default value (corresponds to HESS). The required value changes
between cameras (GCT = 1.05).
params['integration_sigamp'] - Amplitude in ADC counts above pedestal
at which a signal is considered as significant (separate for
high gain/low gain).
geom_dict : dict
Dict of pixel geometry for each telescope. Leave as None for automatic
calculation when it is required.
dict[(num_pixels, focal_length)] = `ctapipe.io.CameraGeometry`
Returns
-------
calibrated : container
A new `ctapipe` event container containing the dl1 information, and a
reference to all other information contained in the original event
container.
"""
# Obtain relevent calibrator
switch = {
'hessio':
partial(mc.calibrate_mc, event=event, params=params)
}
try:
calibrator = switch[event.meta.source]
except KeyError:
log.exception("no calibration created for data origin: '{}'"
.format(event.meta.source))
raise
calibrated = copy(event)
# Add dl1 to the event container (if it hasn't already been added)
try:
calibrated.add_item("dl1", RawData())
calibrated.dl1.run_id = event.dl0.run_id
calibrated.dl1.event_id = event.dl0.event_id
calibrated.dl1.tels_with_data = event.dl0.tels_with_data
calibrated.dl1.calibration_parameters = params
except AttributeError:
pass
# Fill dl1
calibrated.dl1.tel = dict() # clear the previous telescopes
for telid in event.dl0.tels_with_data:
nchan = event.dl0.tel[telid].num_channels
npix = event.dl0.tel[telid].num_pixels
calibrated.dl1.tel[telid] = CalibratedCameraData(telid)
calibrated.dl1.tel[telid].num_channels = nchan
calibrated.dl1.tel[telid].num_pixels = npix
# Get geometry
int_dict, inverted = integrator_dict()
geom = None
cam_dimensions = (event.dl0.tel[telid].num_pixels,
event.meta.optical_foclen[telid])
# Check if geom is even needed for integrator
if inverted[params['integrator']] in integrators_requiring_geom():
if geom_dict is not None and telid in geom_dict:
geom = geom_dict[telid]
else:
log.debug("[calib] Guessing camera geometry")
geom = CameraGeometry.guess(*event.meta.pixel_pos[telid],
event.meta.optical_foclen[telid])
log.debug("[calib] Camera geometry found")
if geom_dict is not None:
geom_dict[telid] = geom
pe, window, data_ped, peakpos = calibrator(telid=telid, geom=geom)
calibrated.dl1.tel[telid].pe_charge = pe
calibrated.dl1.tel[telid].peakpos = peakpos
for chan in range(nchan):
calibrated.dl1.tel[telid].integration_window[chan] = window[chan]
calibrated.dl1.tel[telid].pedestal_subtracted_adc[chan] = \
data_ped[chan]
return calibrated
def calibrate_event(event, params, geom_dict=None):
"""
Generic calibrator for events. Calls the calibrator corresponding to the
source of the event, and stores the dl1 (pe_charge) information into a
new event container.
Parameters
----------
event : container
A `ctapipe` event container
params : dict
REQUIRED:
params['integrator'] - Integration scheme
params['integration_window'] - Integration window size and shift of
integration window centre
(adapted such that window fits into readout).
OPTIONAL:
params['integration_clip_amp'] - Amplitude in p.e. above which the
signal is clipped.
params['integration_calib_scale'] : Identical to global variable
CALIB_SCALE in reconstruct.c in hessioxxx software package. 0.92 is
the default value (corresponds to HESS). The required value changes
between cameras (GCT = 1.05).
params['integration_sigamp'] - Amplitude in ADC counts above pedestal
at which a signal is considered as significant (separate for
high gain/low gain).
geom_dict : dict
Dict of pixel geometry for each telescope. Leave as None for automatic
calculation when it is required.
dict[(num_pixels, focal_length)] = `ctapipe.io.CameraGeometry`
Returns
-------
calibrated : container
A new `ctapipe` event container containing the dl1 information, and a
reference to all other information contained in the original event
container.
"""
# Obtain relevent calibrator
switch = {'hessio': partial(mc.calibrate_mc, event=event, params=params)}
try:
calibrator = switch[event.meta.source]
except KeyError:
log.exception("no calibration created for data origin: '{}'".format(
event.meta.source))
raise
calibrated = copy(event)
# Add dl1 to the event container (if it hasn't already been added)
try:
calibrated.add_item("dl1", RawData())
calibrated.dl1.run_id = event.dl0.run_id
calibrated.dl1.event_id = event.dl0.event_id
calibrated.dl1.tels_with_data = event.dl0.tels_with_data
calibrated.dl1.calibration_parameters = params
except AttributeError:
pass
# Fill dl1
calibrated.dl1.tel = dict() # clear the previous telescopes
for telid in event.dl0.tels_with_data:
nchan = event.dl0.tel[telid].num_channels
npix = event.dl0.tel[telid].num_pixels
calibrated.dl1.tel[telid] = CalibratedCameraData(telid)
calibrated.dl1.tel[telid].num_channels = nchan
calibrated.dl1.tel[telid].num_pixels = npix
# Get geometry
int_dict, inverted = integrator_dict()
geom = None
cam_dimensions = (event.dl0.tel[telid].num_pixels,
event.meta.optical_foclen[telid])
# Check if geom is even needed for integrator
if inverted[params['integrator']] in integrators_requiring_geom():
if geom_dict is not None and telid in geom_dict:
geom = geom_dict[telid]
else:
log.debug("[calib] Guessing camera geometry")
geom = CameraGeometry.guess(*event.meta.pixel_pos[telid],
event.meta.optical_foclen[telid])
log.debug("[calib] Camera geometry found")
if geom_dict is not None:
geom_dict[telid] = geom
pe, window, data_ped, peakpos = calibrator(telid=telid, geom=geom)
calibrated.dl1.tel[telid].pe_charge = pe
calibrated.dl1.tel[telid].peakpos = peakpos
for chan in range(nchan):
calibrated.dl1.tel[telid].integration_window[chan] = window[chan]
calibrated.dl1.tel[telid].pedestal_subtracted_adc[chan] = \
data_ped[chan]
return calibrated
def calibrate_event(event, params, geom_dict=None):
"""
Generic calibrator for events. Calls the calibrator corresponding to the
source of the event, and stores the dl1 (calibrated_image) information into a
new event container.
Parameters
----------
event : container
A `ctapipe` event container
params : dict
REQUIRED:
params['integrator'] - Integration scheme
params['integration_window'] - Integration window size and shift of
integration window centre
(adapted such that window fits into readout).
OPTIONAL:
params['integration_clip_amp'] - Amplitude in p.e. above which the
signal is clipped.
params['integration_calib_scale'] : Identical to global variable
CALIB_SCALE in reconstruct.c in hessioxxx software package. 0.92 is
the default value (corresponds to HESS). The required value changes
between cameras (GCT = 1.05).
params['integration_sigamp'] - Amplitude in ADC counts above pedestal
at which a signal is considered as significant (separate for
high gain/low gain).
geom_dict : dict
Dict of pixel geometry for each telescope. Leave as None for automatic
calculation when it is required.
dict[(num_pixels, focal_length)] = `ctapipe.io.CameraGeometry`
Returns
-------
calibrated : container
A new `ctapipe` event container containing the dl1 information, and a
reference to all other information contained in the original event
container.
"""
# Obtain relevent calibrator
switch = {
'hessio':
partial(mc.calibrate_mc, event=event, params=params)
}
try:
calibrator = switch[event.meta['source']]
except KeyError:
log.exception("no calibration created for data origin: '{}'"
.format(event.meta['source']))
raise
# KPK: should not copy the event here! there is no reason to
# Copying is
# up to the user if they want to do it, not in the algorithms.
# calibrated = copy(event)
# params stored in metadata
event.dl1.meta.update(params)
# Fill dl1
event.dl1.reset()
for telid in event.dl0.tels_with_data:
nchan = event.inst.num_channels[telid]
npix = event.inst.num_pixels[telid]
event.dl1.tel[telid] = CalibratedCameraContainer()
# Get geometry
int_dict, inverted = integrator_dict()
geom = None
# Check if geom is even needed for integrator
if inverted[params['integrator']] in integrators_requiring_geom():
if geom_dict is not None and telid in geom_dict:
geom = geom_dict[telid]
else:
log.debug("[calib] Guessing camera geometry")
geom = CameraGeometry.guess(*event.inst.pixel_pos[telid],
event.inst.optical_foclen[telid])
log.debug("[calib] Camera geometry found")
if geom_dict is not None:
geom_dict[telid] = geom
pe, window, data_ped, peakpos = calibrator(telid=telid, geom=geom)
tel = event.dl1.tel[telid]
tel.calibrated_image = pe
tel.peakpos = peakpos
for chan in range(nchan):
tel.integration_window[chan] = window[chan]
tel.pedestal_subtracted_adc[chan] = data_ped[chan]
return event
