Flavor Transformations plus Earth matter #258

jpkneller · 2023-06-16T13:43:31Z

A rewrite of the flavor transformation classes.

Includes a new method called get_probabilities which returns a matrix of transition probabilities. This method is meant t replace the 8 prob_xy methods in the previous versions of the classes
To avoid duplicating code, there are now two new intermediate classes called ThreeFlavorTransformation and FourFlavorTransformation between the abstract base class FlavorTransformation and the various flavor transformation classes.
The ThreeFlavorTransformation class has two methods: one called vacuum and the other Earth_matter. These methods return the elements of the D matrix previously calculated in the flavor transformation classes. The Earth_matter method does an earth-matter effect calculation using the Sqa3Earth module. Compiling that module is described elsewhere. The FourFlavorTransformation class only has the vacuum method.
The flavor transformation classes now have an optional argument to the constructor called AltAz. This is an astropy object which gives the altitude and azimuth of the supernova relative to the observer / detector. If not None, the earth matter effect will be done (if the altitude is less than zero)
The previous version of the flavor transformation classes had two optional arguments to the constructor called mixing_angles and mass hierarchy. These have been replaced with the optional argument called mix_params which must be an instance of either the MixingParamters3Flavor or MixingParamters4Flavor classes in neutrino.py. If None, the default is to the use the normal mass hierarchy as before.

A rewrite of the flavor transformation classes. 1) Includes a new method called get_probabilities which returns a matrix of transition probabilities. This method is meant t replace the 8 prob_xy methods in the previous versions of the classes 2) To avoid duplicating code, there are now two new intermediate classes called ThreeFlavorTransformation and FourFlavorTransformation between the abstract base class FlavorTransformation and the various flavor transformation classes. 3) The ThreeFlavorTransformation class has two methods: one called vacuum and the other Earth_matter. These methods return the elements of the D matrix previously calculated in the flavor transformation classes. The Earth_matter method does an earth-matter effect calculation using the Sqa3Earth module. Compiling that module is described elsewhere. The FourFlavorTransformation class only has the vacuum method. 4) The flavor transformation classes now have an optional argument to the constructor called AltAz. This is an astropy object which gives the altitude and azimuth of the supernova relative to the observer / detector. If not None, the earth matter effect will be done (if the altitude is less than zero) 5) The previous version of the flavor transformation classes had two optional arguments to the constructor called mixing_angles and mass hierarchy. These have been replaced with the optional argument called mix_params which must be an instance of either the MixingParamters3Flavor or MixingParamters4Flavor classes in neutrino.py. If None, the default is to the use the normal mass hierarchy as before.

…on class inputs The transformation_type string argument to the generate_* functions has been replaced with a flavor_transformation argument which can be either a string or a flavor transformation class. If it is a string the get_transformation function is called to identify the class.

Added __str__ methods

Now uses the get_probabilities method of the flavor transformation classes rather than the prob_xy methods.

…/adiabatic_basis.cpp

…tl/algorithm2.h

…newpy/SQA/src/mstl/physics/units and constants.cpp

…stl/math2/number.h

…newpy/SQA/src/mstl/math2/spline/discontinuous.cpp

…py/SQA/src/mstl/math2/analysis/algorithm3.cpp

… python/snewpy/SQA/src/mstl/math2/algebra/column and row vectors.cpp

…to python/snewpy/EMEWS/src/mstl/math2/algebra/column and row vectors.cpp

…ewpy/EMEWS/src/mstl/math2/analysis/algorithm3.cpp

…/snewpy/EMEWS/src/mstl/math2/spline/discontinuous.cpp

…/snewpy/EMEWS/src/mstl/physics/units and constants.cpp

Moved EMEWS to separate repository

EMEWS is now optional

jpkneller · 2023-11-17T16:55:13Z

The Earth-matter effect module is now a separate repo and is not required for SNEWPY to work.

Resolving conflicts in flux.py and snowglobes.py

Sheshuk

Hi Jim, thank you for updating this!
Having EMEWS decoupled makes it much easier for users to explore the SNEWPY without having to run additional steps (compilation, installation etc), in case they don't need this functionality.

Also it's great to finally have the transformation matrices in place!

I have two suggestions:

As you can see, the tests are failing, because they try to use mh parameter in the constructors of the FlavorTransformations. So for backward compatibility we should still use mh input, if it is provided (and add a note in the code to clean it up for SNEWPY v2.0)
Move the repeating code, which fills the transformation matrix from p_ee, pbar_ee... to a separate function (I added an example in comments) and use it. Having less repeating code means less space for possible copy-paste errors.

Sheshuk · 2023-11-20T11:41:25Z