SNEWS2 / snewpy

A Python package for working with supernova neutrinos
https://snewpy.readthedocs.io
BSD 3-Clause "New" or "Revised" License
26 stars 19 forks source link

Make `FlavorTransformation` produce a transformation matrix #242

Open Sheshuk opened 1 year ago

Sheshuk commented 1 year ago

Currently FlavorTransformation can only produce prob_ee, prob_ex, prob_xe, prob_xx (and their bar counterparts). I think it would be useful to have a method to output a transformation matrix, for all the four flavors, like this:

M = [
[prob_ee, prob_ex,          0,          0],
[prob_xe, prob_xx,          0,          0],
[      0,       0, prob_eebar, prob_xebar],
[      0,       0, prob_eebar, prob_xebar]]

and its permutation depending on the order of our Flavor enum.