HillasReconstructor¶
-
class
ctapipe.reco.hillas_reconstructor.
HillasReconstructor
(**kwargs)[source]¶ Bases:
ctapipe.reco.reco_algorithms.Reconstructor
class that reconstructs the direction of an atmospheric shower using a simple hillas parametrisation of the camera images it provides a direction estimate in two steps and an estimate for the shower’s impact position on the ground.
so far, it does neither provide an energy estimator nor an uncertainty on the reconstructed parameters
Methods Summary
__call__
(event)Perform the full shower geometry reconstruction on the input event.
estimate_core_position
(event, hillas_dict, …)Estimate the core position by intersection the major ellipse lines of each telescope.
estimate_direction
(hillas_planes)calculates the origin of the gamma as the weighted average direction of the intersections of all hillas planes
estimate_h_max
(hillas_planes)Estimate the max height by intersecting the lines of the cog directions of each telescope.
initialize_hillas_planes
(hillas_dict, …)Creates a dictionary of
HillasPlane
from a dictionary of hillas parametersMethods Documentation
-
__call__
(event)[source]¶ Perform the full shower geometry reconstruction on the input event.
- Parameters
- eventcontainer
-
estimate_core_position
(event, hillas_dict, array_pointing, corrected_angle_dict, hillas_planes)[source]¶ Estimate the core position by intersection the major ellipse lines of each telescope.
- Parameters
- hillas_dict: dict[HillasContainer]
dictionary of hillas moments
- array_pointing: SkyCoord[HorizonFrame]
Pointing direction of the array
- Returns
- core_x: u.Quantity
estimated x position of impact
- core_y: u.Quantity
estimated y position of impact
Notes
The part of the algorithm taking into account divergent pointing mode and the usage of a corrected psi angle is explained in [gasparetto] section 7.1.4.
-
estimate_direction
(hillas_planes)[source]¶ calculates the origin of the gamma as the weighted average direction of the intersections of all hillas planes
- Returns
- gammashape (3) numpy array
direction of origin of the reconstructed shower as a 3D vector
- crossingsshape (n,3) list
an error estimate
-
estimate_h_max
(hillas_planes)[source]¶ Estimate the max height by intersecting the lines of the cog directions of each telescope.
- Returns
- astropy.unit.Quantity
the estimated max height
-
initialize_hillas_planes
(hillas_dict, subarray, telescopes_pointings, array_pointing)[source]¶ Creates a dictionary of
HillasPlane
from a dictionary of hillas parameters- Parameters
- hillas_dictdictionary
dictionary of hillas moments
- subarrayctapipe.instrument.SubarrayDescription
subarray information
- telescopes_pointings: dictionary
dictionary of pointing direction per each telescope
- array_pointing: SkyCoord[AltAz]
pointing direction of the array
Notes
The part of the algorithm taking into account divergent pointing mode and the correction to the psi angle is explained in [gasparetto] section 7.1.4.
-