pisa

pisa copied to clipboard

Errors running in single precision mode with non-Asimov datasets

14

Running `hypo_testing.py` for NMO studies and getting errors with parameters going to NaN when running in single precision. I modified the error message so it would tell me which parameter...

steven-j-wren

plotter.py: 2D y=log-scaled energy, x=coszen plots end up log(log(E))

3

apparently a bug in linlog logic; removing the imshow section, plots look fine

jllanfranchi

param selections requires all selections in one stage

2

See issue #306 for more discussion

jllanfranchi

This has hooks all the way down into c++ code where the earth model is carried around everywhere, even if vacuum oscillations are being used. Some options would be to...

jllanfranchi

Unify mc/cpu and mc/gpu

1

The two stages mc/cpu.py and mc/gpu.py contain a lot of duplicate code and therefore should be merged into one more general stage.

philippeller

kde_histogramdd, test need to be made more generic and have hacks removed

2

* requires specific binning dimension names * reordering dimensions for PID is performed in ad hoc, arcane fashion * energy units are forced to be GeV or it fails *...

jllanfranchi

`DistributionMaker.get_outputs(return_sum=True)`: first Map name is used (others discarded)

3

First Map name is kept (others discarded) `DistributionMaker.get_outputs(return_sum=True)`. (My first thoughts are that this bug may be in `DistributionMaker` but also might be in `Map` and/or `MapSet`.) Options: 1. rename...

jllanfranchi

Be able to use either the current implementation (Improved Sheather Jones pilot density estimate + Abramov bandwidth selection) or something from Sebastian Schoenen's KDE module.

jllanfranchi

What I have implemented now I always knew might be problematic, so this is a TODO from the code that should be discussed among the group, and particularly after several...

jllanfranchi

Locks for Param, ParamSet, & ParamSelector (used by Stage, Pipeline, DistributionMaker)

1

* registering parent sets one's lock to be parent's lock object * disconnecting parent grabs new lock * `@synchronized` decorator so we don't have to copy-paste and add lots of...

jllanfranchi