High Level Analysis (high_level)

Introduction

Library containing the functions and classes to perform LST high level analysis.

lstchain.high_level.hdu_table Module

Module containing functions to create DL3 event lists and DL3 index files, following recommendations from GADF v0.3

For reference to the documentation of GADF v0.3 follow various sections in, https://gamma-astro-data-formats.readthedocs.io/en/latest/

Functions

add_icrs_position_params(data, source_pos, ...)

Updating data with ICRS position values of reconstructed positions in RA, DEC coordinates and add a column, theta, on separation from the true source position.

create_event_list(data, run_number, ...)

Create the BinTableHDUs for storing DL3 event data.

create_hdu_index_hdu(file_list, hdu_index_file)

Create the hdu index table and write it to the given file.

create_obs_index_hdu(file_list, ...)

Create the obs index table and write it to the given file.

fill_reco_altaz_w_expected_pos(data)

Fill the reconstructed alt, az positions with the expected source positions, for source-dependent analysis.

get_timing_params(data[, epoch])

Retrieve some timing parameters for the DL3 event list as a dict and also have a common time_utc object for use in other functions.

get_pointing_params(data, source_pos, time_utc)

Convert the telescope pointing directions into ICRS frame of reference, and average them for the run (excluding the first events, for which there may be some mispointing due to not-yet-stable tracking).

get_timing_params(data[, epoch])

Retrieve some timing parameters for the DL3 event list as a dict and also have a common time_utc object for use in other functions.

lstchain.high_level.interpolate Module

Functions

check_in_delaunay_triangle(irfs, data_params)

From a given list of IRFs as grid points used for interpolation, retrieve the Delaunay triangulation list of IRFs, where the simplex includes the target points in data_params.

compare_irfs(irfs)

Compare the given list of IRFs with data binning and relevant metadata values.

get_nearest_az_node(irf_params_sel, index, ...)

Check to see if a given IRF node overlaps with another node, in the interpolation parameter space, and to select based on the azimuth angle, the nearest node to the target.

interp_params(params_list, data)

From a given list of angular parameters, to be used for interpolation, take values from a given data dict.

interpolate_cuts(cuts, grid_points, target_point)

Interpolates a grid of GH_CUTS or AL_CUTS tables to a target-point.

interpolate_irf(irfs, data_pars[, interp_method])

Using pyirf functions with a list of IRFs and parameters to compare with data, to interpolate over in a selected interpolation parameter space.

load_irf_grid(irfs, extname, interp_col[, ...])

From a given list of IRFs, load the list of IRF data values that can be interpolated.

lstchain.high_level.significance_calculation Module

Collection of core analysis methods

Functions

analyze_on_off(config)

Extracts the theta2 plot of a dataset taken with ON/OFF observations

analyze_wobble(config)

Extracts the theta2 plot of a dataset taken with wobble observations

setup_logging([verbosity])

Setup logger console and file descriptors