Jan Weldert

hypo_testing.py doesn't work with osc.prob3gpu and reco.vbwkde

9

Using the hypo_testing.py script with a pipeline containing osc.prob3gpu and reco.vbwkde delivers the following error: pycuda._driver.LogicError: cuFuncSetBlockShape failed: invalid resource handle This problem is probably related to https://github.com/WIPACrepo/pisa/blob/cake_pre_opensource/pisa/utils/gaussians.py : ```...

Variable binning options

2

When defining a multidimensional binning we combine one dimensional binnings in a homogeneous grid. So if we for example bin our data in energy, coszen, and PID, we have to...

At the moment we only use the correlation between Daemonflux parameters. There is no general method that allows for parameters to be correlated. That would require to pass a covariance...

Supporting Bayesian Analysis

1

At the moment we focus on the frequentist approach. To also support Bayesian analyses we would mainly need some sort of MCMC sampling tool and its integration into the analysis...

Making the bin predictions differentiable was suggested multiple times now. That would for example help minimizers and allow to extract the Fisher-information at the best-fit point. This will be a...

scipy.interpolate.interp1d may be removed in future scipy version

3

We use `from scipy.interpolate import interp1d` in multiple parts of the code, namely - [https://github.com/icecube/pisa/blob/cfeb9b45a38f981059197c293f1a53aa5ab81a39/pisa/core/prior.py#L14](https://github.com/icecube/pisa/blob/cfeb9b45a38f981059197c293f1a53aa5ab81a39/pisa/core/prior.py#L14) - [https://github.com/icecube/pisa/blob/cfeb9b45a38f981059197c293f1a53aa5ab81a39/pisa/stages/background/atm_muons.py#L7](https://github.com/icecube/pisa/blob/cfeb9b45a38f981059197c293f1a53aa5ab81a39/pisa/stages/background/atm_muons.py#L7) - [https://github.com/icecube/pisa/blob/cfeb9b45a38f981059197c293f1a53aa5ab81a39/pisa/utils/spline_smooth.py#L10](https://github.com/icecube/pisa/blob/cfeb9b45a38f981059197c293f1a53aa5ab81a39/pisa/utils/spline_smooth.py#L10) - The new stage introduced in #https://github.com/icecube/pisa/pull/798 At the moment we...

First version of an MCMC sampler

2

First version of an simple MCMC sampler as part of the analysis class. You need to manually install emcee at the moment. This PR is part of #804. I also...

The Prior class uses a lot of lambda functions and not top-module level functions (e.g. llh within the __init_* functions). That makes it impossible to pickle them. That is an...