The main part of the ray tracer can be found in raytracer2018.nim. To run it, different files like X-Ray reflectivity of different materials are needed but can be requested. The possibility to change from BabyIAXO to CAST is in the very bottom of the file.
The calculation of the gas phase can be found in the folder axionMass and the emission rates from the sun are calculated in readOpacityFile.nim and readSolarModel.nim.
Here is for example the flux of a DFSZ axion with 57Fe and plasmon contributions:
And here the same for a KSVZ axion:
The solar model used for the calculations is the AGSS09 model:
- https://www.imperial.ac.uk/astrophysics/research/the-sun/data/
- https://wwwmpa.mpa-garching.mpg.de/~aldos/solar_main.html (contains multiple & has the data files). The file for the solar model is also found in ReadSolarModel/resources/AGSS09_solar_model.dat)
The calculations for the differential flux further use the Opacity Project (OP) in particular version 3.3:
http://cdsweb.u-strasbg.fr/topbase/TheOP.html
Warning: The server is extremely slow! Use something like aria2c
or similar to download using multiple connections in parallel.
They need to be downloaded manually and placed somewhere accessible
(currently the path is hard coded to be relative to the
readOpacityFile binary, but this will be changed soon).
For version 2.0 a second host is: https://code.osu.edu/wiag/opacities.osc.edu/-/tree/master/OPserver/OPCD_2.0
After downloading and extracting the .tar archive of the OPCD data,
you also need to extract the files in the mono subdirectory:
for f in *.xz; do gzip -d $f; doneFinally, the fm01.mesh file needs to be modified slightly. The first
line:
1.0553E-04 10000
must be modified to:
u
1.0553E-04
After that and fixing the paths in readOpacityFile it should run.