The trigger rate application (and concept) has been developed in these PRs: #801, #935, #1000, #1007 according to Konrad's recipe in #642.
Current functionalities
The simtools-calculate-trigger-rate tool calculates the trigger rate for a simulation. It takes as input preferentially the .simtel.zst file. It can also take the histogram file .hdata.zst but in this case, the user has to pass more information about the configuration used for the simulation in the first place;
The energy dependent trigger rate is saved in a ecsv as output and the total trigger rate with uncertainty is printed out in the terminal (if --save_tables is passed);
One might pass more than one input file or a list with input files. Two options are then possible: either the tool calculates the trigger rate individually for each file (default) or adds ups the histograms before estimating an overall trigger rate (by passing --stack_files to the application);
An alternative method for calculating the area thrown, important for the trigger rate estimate, is available by passing --area_from_distribution. In this case, the area thrown is estimated based on the distribution of events on the camera and not the input maximum distance of simulated area.
Missing functionalities
Allow the production of a report (pdf file) with the main results of the trigger rate estimate, i.e. the energy dependent trigger rate, some statistics (percentile) and the total trigger rate;
Calculate trigger rate for other particles (and for the CR all particle spectrum), as explained in #931. To implement this, one would need to first generalize the methods available for the Power Law function to the other functions in ctao-dpps-cosmic-ray-spectra package.
If strictly necessary, implement Konrad's extract_simtel_array.C in simtools and allow the trigger rate to be calculated for subarrays (and ultimately single telescopes), given a production simulation input (see #1001 for details).
It is possible to write a workflow that first simulate protons (or all needed particles) for a single telescope and then calculate directly the trigger rate for this same telescope using both the simtools-simulate-showers-for-trigger-rates and simtools-calculate-trigger-rate tools.
History
The trigger rate application (and concept) has been developed in these PRs: #801, #935, #1000, #1007 according to Konrad's recipe in #642.
Current functionalities
simtools-calculate-trigger-rate
tool calculates the trigger rate for a simulation. It takes as input preferentially the.simtel.zst
file. It can also take the histogram file.hdata.zst
but in this case, the user has to pass more information about the configuration used for the simulation in the first place;--save_tables
is passed);--stack_files
to the application);--area_from_distribution
. In this case, the area thrown is estimated based on the distribution of events on the camera and not the input maximum distance of simulated area.Missing functionalities
extract_simtel_array.C
in simtools and allow the trigger rate to be calculated for subarrays (and ultimately single telescopes), given a production simulation input (see #1001 for details).simtools-simulate-showers-for-trigger-rates
andsimtools-calculate-trigger-rate
tools.