v2.2
What's Changed
GSL Integration implementation, TRestAxionMagneticField ReMap feature and HOT bug fix by @jgalan in #82
-
CRITICAL fix where the axion mass was given in keV inside the TRestAxionFieldPropagationProcess while it should have been expressed in eV. This was leading to wrong results.
-
TRestAxionMagneticField::ReMap
method added. This method allows to remap the field to a new grid that it is necessarily coarser than the original. -
TRestAxionField::GammaTransmissionFieldMapProbability
method added. This method will integrate the axion-field using GSL library integration methods. The advantage is that we define the desired accuracy. This method will return a pair with the probability and the corresponding error. -
TRestAxionField
implements now a pointer tofMagneticField
that needs to be initialized using TRestAxionField::AssignMagneticField so that the new methodGammaTransmissionFieldMapProbability
knows how to evaluate the magnetic field alongx
. -
TRestAxionMagneticField::GetTransversalComponentInParametricTrack(Double_t x)
new method implemented that allows to evaluate the field along a parametric track. The track is calculated using the field boundaries by providing the particle position and direction to the magnetic field, using the new methodSetTrack( const TVector3 &position, const TVector3 &direction);
. -
TRestAxionFieldPropagationProcess
now it integrates the probability using the new GSL integration method. The data members in this class have been updated to be able to set-up the integration method. -
TRestAxionFieldPropagationProcess
implements afReMap
data member that allows the process to re-map the magnetic field to a new grid size.
Full Changelog: v2.1...v2.2