Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add option to calculate errors on the fit params assuming Wilks' theorem. #170

Open
chiarabellenghi opened this issue Aug 30, 2023 · 0 comments
Assignees
Labels
enhancement New feature or request

Comments

@chiarabellenghi
Copy link
Member

skyllh/skyllh/core/analysis.py

Lines 1071 to 1124 in 074c508

def do_trial(
self,
rss,
mean_n_bkg_list=None,
mean_n_sig=0,
bkg_kwargs=None,
sig_kwargs=None,
minimizer_status_dict=None,
tl=None,
**kwargs):
"""This method performs an analysis trial by generating a
pseudo data sample with background events and possible signal events
via the :meth:`generate_pseudo_data` method, and performs the analysis
on that random pseudo data sample by calling the
:meth:`do_trial_with_given_pseudo_data` method.
Parameters
----------
rss : instance of RandomStateService
The instance of RandomStateService instance to use for generating
random numbers.
mean_n_bkg_list : list of float | None
The mean number of background events that should be generated for
each dataset. If set to None (the default), the background
generation method needs to obtain this number itself.
mean_n_sig : float
The mean number of signal events that should be generated for the
trial.
bkg_kwargs : dict | None
Additional keyword arguments for the `generate_events` method of the
background generation method class. An usual keyword argument is
`poisson`.
sig_kwargs : dict | None
Additional keyword arguments for the `generate_signal_events` method
of the `SignalGenerator` class. An usual keyword argument is
`poisson`.
minimizer_status_dict : dict | None
If a dictionary is provided, it will be updated with the minimizer
status dictionary.
tl : instance of TimeLord | None
The optional instance of TimeLord that should be used to time
individual tasks.
**kwargs : dict
Additional keyword arguments are passed to the
:meth:`do_trial_with_given_pseudo_data` method.
Returns
-------
recarray : instance of numpy record ndarray
The numpy record ndarray holding the result of the trial.
See the documentation of the
:py:meth:`~skyllh.core.analysis.Analysis.do_trial_with_given_pseudo_data`
method for further information.
"""

@chiarabellenghi chiarabellenghi self-assigned this Aug 30, 2023
@chiarabellenghi chiarabellenghi added the enhancement New feature or request label Sep 29, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

1 participant