Skip to content

v2.2

Compare
Choose a tag to compare
@jgalan jgalan released this 12 Feb 13:56
· 160 commits to master since this release
3800d90

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 to fMagneticField that needs to be initialized using TRestAxionField::AssignMagneticField so that the new method GammaTransmissionFieldMapProbability knows how to evaluate the magnetic field along x.

  • 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 method SetTrack( 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 a fReMap data member that allows the process to re-map the magnetic field to a new grid size.

Full Changelog: v2.1...v2.2