Fourier and machine learning tools for constraining axion-like-particles in astrophysical observations as described in https://arxiv.org/abs/1808.05916 .
Python 2.7, PyXspec 2
for Fourier Analysis: https://github.com/NFFT/nfft
for Machine Learning: tensorflow 1.5
-
Use mathematica notebook PhotonAxionConversionCluster.nb to generate output: survivalProbs_*Bfield.txt
-
Delete all survivalProbs_*Bfield.mod files
-
Use python ALPmod.py to generate survivalProbs_*Bfield.mod table models as input for PyXspec.
-
Run python PyXspec.py to perform PyXspec scan over all *.mod files. Output: g_chisq.txt
-
Run python Analyze.py to analyze g_chsiq.txt. Outputs: histo_*.pdf and g_deltachi2.pdf
-
Residuals can also be fed to the machine learning algorithm described in tf.py.