pisa
pisa copied to clipboard
icecube
Monte Carlo-based data analysis
See comment from @steven-j-wren here: https://github.com/jllanfranchi/pisa/pull/394#issuecomment-335250536 In summary: "standard `honda.py` systematics should be included in `mceq.py`. I think this could be as simple as importing them from `honda.py` but I'm...
Unit test `test_Data()` fails due to missing data files. Since I'm in the process of removing test, data, and all but example settings files from the PISA "master" branch (i.e....
originally shows up here: https://github.com/jllanfranchi/pisa/blob/37f168b2d98bd51d8ab90a66c06a80e1a179ef40/pisa/core/events_pi.py#L309-L315 where the np.append was being done, but that will operate on a dict but you can't np.concat dicts. I think whoever wrote this got confused...
`simple_data_loader` is becoming increasingly bloated with bespoke behaviour and is no longer really simple. We should make a dedicated data loader for our IceCube PISA case containing all this functionality,...
Just putting this idea here not to forget: For scans (e.g. 1-d nutau norm, 2-d contours, ...) we always start the minimizer from nominal, but it would converge much faster...
Running ```bash export CUDA_VISIBLE_DEVICES=0 export PISA_RESOURCES=/fastio/justin/pisa_resources export PISA_FTYPE=fp32 for t in {1..100} do OUTDIR=/tmp/test${t} $PISA/pisa/scripts/analysis.py discrete_hypo \ --h0-pipeline settings/pipeline/example_mc.cfg \ --h0-param-selections=ih \ --h1-param-selections=nh \ --data-param-selections=nh \ --data-is-mc \ --min-settings settings/minimizer/l-bfgs-b_ftol2e-5_gtol1e-5_eps1e-4_maxiter200.json...
Issues are: * setting of theta23 start values (both mirroring for 2nd fit + shifting in case of 3rd fit) may violate bounds * the outcome of the 2nd fit...
Setting `np.errstate` in https://github.com/jllanfranchi/pisa/blob/dev/pisa/core/map.py#L521 doesn't help when `hist` is an `uncertainties.unumpy.uarray` (I think this has been an open issue for a long time: https://github.com/lebigot/uncertainties/issues/12). This can be fixed by applying...
In `map.py` functions such as `compare`, there are cases where one map is divided by another. If the denominator map has any bins with 0 entries, this fails with a...
It currently forces the user to: * specify some settings (at least a discrete selection, otherwise there's an exception raised since "there's nothing of substance to be done") for the...