simtools-production-derive-corsika-limits#

Derive CORSIKA configuration limits for energy, core distance, and viewcone radius.

This tool determines configuration limits based on triggered events from broad-range simulations. It supports setting:

  • ERANGE: Derives the lower energy limit; upper limit is user-defined.

  • CSCAT: Derives the upper core distance; lower limit is user-defined.

  • VIEWCONE: Derives the viewcone radius; lower limit is user-defined.

Limits are computed based on a configurable maximum event loss fraction. Results are provided as a table with the following columns:

  • particle_type: Particle type (e.g., gamma, proton, electron).

  • telescope_ids: List of telescope IDs used in the simulation.

  • zenith: Zenith angle.

  • azimuth: Azimuth angle.

  • nsb: Night sky background level

  • layout: Layout of the telescope array used in the simulation.

  • lower_energy_limit: Derived lower energy limit.

  • upper_radius_limit: Derived upper radial distance limit.

  • viewcone_radius: Derived upper viewcone radius limit.

The input event data files are generated using the application simtools-generate-simtel-event-data and are required for each point in the lookup table.

Command line arguments#

event_data_files (str, required)

Path to a file containing event data files derived with ‘simtools-generate-simtel-event-data’.

telescope_ids (str, required)

Path to a file containing telescope configurations.

loss_fraction (float, required)

Maximum event-loss fraction for limit computation.

plot_histograms (bool, optional)

Plot histograms of the event data.

output_file (str, optional)

Path to the output file for the derived limits.

Example#

Derive limits for a given file with a specified loss fraction.

simtools-production-derive-corsika-limits \\
    --event_data_files path/to/event_data_files.yaml \\
    --telescope_ids path/to/telescope_configs.yaml \\
    --loss_fraction 1e-6 \\
    --plot_histograms \\
    --output_file corsika_simulation_limits_lookup.ecsv
production_derive_corsika_limits.main()[source]#

Derive limits for energy, radial distance, and viewcone.