REST-for-Physics
v2.3
Rare Event Searches ToolKit for Physics
|
A basic class to define analytical axion-photon conversion calculations for axion helioscopes.
TRestAxionField is a class used to calculate the axion-photon mixing and determine the probability of the particle being in a photon state after propagating inside a magnetic field.
The class allows to assign a buffer gas using a TRestAxionBufferGas and a magnetic field map loaded through TRestAxionMagneticField. If these objects have been assigned, the methods implemented in this class may consider those - if needed - inside the calculation.
In practice this class provides different versions of the method GammaTransmissionProbability that allows to calculate the axion-photon probability using different strategies.
For calculations inside a constant magnetic field one may simply invoke the following code which will launch the calculation in vacuum, for a field of 2T, a coherence length of 10000mm, and an axion energy and mass of 4.2keV and 0.1eV, respectively.
It is possible to reduce the number of arguments we are giving to this function by assigning few data member values present inside TRestAxionField as follows:
Indeed, the complete version with all arguments will also update the data member values, so that next calls to this method will use the same magnetic field, coherence length and energy. In the following example the second call to GammaTransmissionProbability will use again the same values for magnetic field, coherence length and energy from the preceding call.
The axion-photon probability can also be calculated in the sin of a gaseous medium. For that we need to assign a buffer gas instance containing the relevant gas properties in the form of a TRestAxionBufferGas.
The following piece of code will assign the gas properties to the field calculation as defined inside the bufferGases.rml
official file. The gas provides the absorption and equivalent photon mass as a funtion of the energy.
Once we have assigned a buffer gas to the class we may return back to the vacuum state by just assigning a nullptr to the gas.
There are two main strategies presently implemented inside this class to integrate the axion field along an unhomogeneous magnetic field. The first one uses a simple trapezoidal integration method with the advantage that the computing time is controlled by length of the integration step, but no error calculation is given. The second method advantage is that in principle it allows to fix the accuracy of our calculation and an error of the integration is returned by the GSL integration method used. The main disadvantage is that by fixing the accuracy we lose control on the computation time required to reach such accuracy. Here we describe both methods.
Using a external field profile requires to prepare the field profile and pass it along to TRestAxionField::GammaTransmissionProbability method performing the integration. The profile will be placed in a std::vector<Double_t>
that describes the magnetic field profile with a fixed length step. The method TRestAxionMagneticField::GetTransversalComponentAlongPath can be used to generate such profile and calculate the field using the first method, as in the following example:
In a similar way to the constant field case, it would suffice to assign a buffer gas to the TRestAxionField instance to have the calculation performed in a buffer gas.
A new method has been added in which the integration method is given a pointer to the function that will allow evaluation of the magnetic field along the B-profile parameterized track. Thus, it is the integration method the one requesting the points on the parametric curve that will be necessary to perform the integration. The required field evaluation is encoded inside TRestAxionField::Integrand which accesses the magnetic field instance TRestAxionField::fMagneticField in order to recover a parametric track along the path.
The method name using this technique requires thus direct access to a TRestAxionMagneticField instance that must be assigned beforehand. For computing the integral following this strategy one needs to write the following code where we consider the most generic case including a gas buffer medium.
Where the later call will return a std::pair
with the probability and its corresponding error.
A method may be used to determine the masses and gas densities to achieve a continuous scanning TRestAxionField::GetMassDensityScanning. This method will place a new mass or gas density setting at FWHM/2 till it reaches the maximum axion mass specified in the argument.
The following code recovers the axion masses and density settings required for a continuous scan till 0.2eV.
RESTsoft - Software for Rare Event Searches with TPCs
History of developments:
2019-March: First concept and implementation of TRestAxionField class. Javier Galan 2023-December: Implementing methods to recover the axion mass scanning nodes Fran Candon
Definition at line 30 of file TRestAxionField.h.
#include <TRestAxionField.h>
Public Member Functions | |
void | AssignBufferGas (TRestAxionBufferGas *buffGas) |
It assigns a gas buffer medium to the calculation. More... | |
void | AssignMagneticField (TRestAxionMagneticField *mField) |
It assigns a magnetic field to the calculation. More... | |
Double_t | AxionAbsorptionProbability (Double_t Bmag, Double_t Lcoh, Double_t Ea, Double_t ma, Double_t mg=0, Double_t absLength=0) |
On top of calculating the axion absorption probability it will assign new values for the magnetic field (Bmag/T), coherence length (Lcoh/mm) and axion energy (Ea/keV). More... | |
Double_t | AxionAbsorptionProbability (Double_t ma, Double_t mg=0, Double_t absLength=0) |
Performs the calculation of axion-photon absorption probability using directly equation (18) from van Bibber, Phys Rev D Part Fields. 1989. More... | |
Double_t | BL (Double_t Bmag, Double_t Lcoh) |
Performs the calculation of (BL) factor in natural units. More... | |
Double_t | BLHalfSquared (Double_t Bmag, Double_t Lcoh) |
Performs the calculation of (BL/2)^2 factor in natural units. More... | |
std::pair< Double_t, Double_t > | GammaTransmissionFieldMapProbability (Double_t Ea, Double_t ma, Double_t accuracy=1.e-1, Int_t num_intervals=100, Int_t qawo_levels=20) |
Performs the calculation of axion-photon conversion probability using directly equation (28) from J. Redondo and A. Ringwald, Light shinning through walls. https://arxiv.org/pdf/1011.3741.pdf. More... | |
Double_t | GammaTransmissionFWHM (Double_t step=0.00001) |
Performs the calculation of the FWHM for the axion-photon conversion probability computed in TRestAxionField::GammaTransmissionProbability . More... | |
Double_t | GammaTransmissionProbability (Double_t Bmag, Double_t Lcoh, Double_t Ea, Double_t ma, Double_t mg=0, Double_t absLength=0) |
On top of calculating the gamma transmission probability it will assign new values for the magnetic field (Bmag/T), coherence length (Lcoh/mm) and axion energy (Ea/keV). More... | |
Double_t | GammaTransmissionProbability (Double_t Ea, Double_t ma, Double_t mg=0, Double_t absLength=0) |
Performs the calculation of axion-photon conversion probability using directly equation (11) from van Bibber, Phys Rev D Part Fields. 1989. More... | |
Double_t | GammaTransmissionProbability (std::vector< Double_t > Bmag, Double_t deltaL, Double_t Ea, Double_t ma, Double_t mg=0, Double_t absLength=0) |
Performs the calculation of axion-photon conversion probability using directly equation (28) from J. Redondo and A. Ringwald, Light shinning through walls. https://arxiv.org/pdf/1011.3741.pdf. More... | |
Double_t | GetAxionEnergy () const |
Double_t | GetCoherenceLength () const |
Double_t | GetMagneticField () const |
std::vector< std::pair< Double_t, Double_t > > | GetMassDensityScanning (std::string gasName="He", Double_t maMax=0.15, Double_t rampDown=5.0) |
This method determines the proper densities to be used in an axion helioscope experiment in order to achieve a continuous axion mass scan. More... | |
void | SetAxionEnergy (Double_t e) |
void | SetBufferGas (TRestAxionBufferGas *buffGas) |
It assigns a gas buffer medium to the calculation. More... | |
void | SetCoherenceLength (Double_t l) |
void | SetDebug (Bool_t v) |
It enables/disables debug mode. More... | |
void | SetMagneticField (Double_t b) |
TRestAxionField () | |
Default constructor. More... | |
~TRestAxionField () | |
Default destructor. More... | |
Static Public Member Functions | |
static double | Integrand (double x, void *params) |
Integrand used for axion-photon probability integration. More... | |
Private Member Functions | |
std::pair< Double_t, Double_t > | ComputeOffResonanceIntegral (Double_t q, Double_t Gamma, Double_t accuracy, Int_t num_intervals, Int_t qawo_levels) |
Performs the calculation of axion-photon conversion probability using directly equation (28) from J. Redondo and A. Ringwald, Light shinning through walls. https://arxiv.org/pdf/1011.3741.pdf. More... | |
std::pair< Double_t, Double_t > | ComputeResonanceIntegral (Double_t Gamma, Double_t accuracy, Int_t num_intervals) |
Performs the calculation of axion-photon conversion probability using directly equation (28) from J. Redondo and A. Ringwald, Light shinning through walls. https://arxiv.org/pdf/1011.3741.pdf. More... | |
void | Initialize () |
Initialization of TRestAxionField class. More... | |
Private Attributes | |
Double_t | fBmag = 2.5 |
The magnetic field in Teslas (used for constant field formulas) More... | |
TRestAxionBufferGas * | fBufferGas = nullptr |
A pointer to the buffer gas definition. More... | |
Bool_t | fDebug = false |
Double_t | fEa = 4.2 |
The energy of the axion in keV. More... | |
Double_t | fLcoh = 10000 |
The coherence lenght (in mm) where the magnetic field is defined (for constant field) More... | |
TRestAxionMagneticField * | fMagneticField = nullptr |
A pointer to the magnetic field definition. More... | |
TRestAxionField::TRestAxionField | ( | ) |
Default constructor.
Definition at line 227 of file TRestAxionField.cxx.
TRestAxionField::~TRestAxionField | ( | ) |
Default destructor.
Definition at line 232 of file TRestAxionField.cxx.
|
inline |
It assigns a gas buffer medium to the calculation.
Definition at line 73 of file TRestAxionField.h.
|
inline |
It assigns a magnetic field to the calculation.
Definition at line 76 of file TRestAxionField.h.
Double_t TRestAxionField::AxionAbsorptionProbability | ( | Double_t | Bmag, |
Double_t | Lcoh, | ||
Double_t | Ea, | ||
Double_t | ma, | ||
Double_t | mg = 0 , |
||
Double_t | absLength = 0 |
||
) |
On top of calculating the axion absorption probability it will assign new values for the magnetic field (Bmag/T), coherence length (Lcoh/mm) and axion energy (Ea/keV).
Definition at line 745 of file TRestAxionField.cxx.
Double_t TRestAxionField::AxionAbsorptionProbability | ( | Double_t | ma, |
Double_t | mg = 0 , |
||
Double_t | absLength = 0 |
||
) |
Performs the calculation of axion-photon absorption probability using directly equation (18) from van Bibber, Phys Rev D Part Fields. 1989.
If not provided as argument, m_gamma it will be attempted to obtain the value from the buffer gas definition. If no buffer gas has been assigned the medium will be assumed to be vacuum.
Ea in keV, ma in eV, mgamma in eV, Lcoh in mm, Bmag in T
mg in eV, absLength in cm-1
The returned value is given for g_ag = 10^-10 GeV-1
Definition at line 684 of file TRestAxionField.cxx.
double TRestAxionField::BL | ( | Double_t | Bmag, |
Double_t | Lcoh | ||
) |
Performs the calculation of (BL) factor in natural units.
Lcoh
should be expressed in mm
, and Bmag
in T
. The result will be given for an axion-photon coupling of 10^{-10} GeV^{-1}
Definition at line 251 of file TRestAxionField.cxx.
double TRestAxionField::BLHalfSquared | ( | Double_t | Bmag, |
Double_t | Lcoh | ||
) |
Performs the calculation of (BL/2)^2 factor in natural units.
Lcoh
should be expressed in mm
, and Bmag
in T
. The result will be given for an axion-photon coupling of 10^{-10} GeV^{-1}
Definition at line 268 of file TRestAxionField.cxx.
|
private |
Performs the calculation of axion-photon conversion probability using directly equation (28) from J. Redondo and A. Ringwald, Light shinning through walls. https://arxiv.org/pdf/1011.3741.pdf.
It integrates the Integrand function defined on the header of TRestAxionField.
This method uses the GSL QAWO method for oscillatory functions. That is the case when the q is not zero, in the off-resonance case.
See https://www.gnu.org/software/gsl/doc/html/integration.html for more details on the integration parameters.
The Gamma function should be expressed in mm-1 since the field map accessed in the integrand is evaluated using mm.
Definition at line 611 of file TRestAxionField.cxx.
|
private |
Performs the calculation of axion-photon conversion probability using directly equation (28) from J. Redondo and A. Ringwald, Light shinning through walls. https://arxiv.org/pdf/1011.3741.pdf.
It integrates the Integrand function defined on the header of TRestAxionField.
This method uses the GSL QAG integration method. We use this method when the cosine function vanishes when q = 0.
See https://www.gnu.org/software/gsl/doc/html/integration.html for more details on the integration parameters.
The Gamma function should be expressed in mm-1 since the field map accessed in the integrand is evaluated using mm.
Definition at line 551 of file TRestAxionField.cxx.
std::pair< Double_t, Double_t > TRestAxionField::GammaTransmissionFieldMapProbability | ( | Double_t | Ea, |
Double_t | ma, | ||
Double_t | accuracy = 1.e-1 , |
||
Int_t | num_intervals = 100 , |
||
Int_t | qawo_levels = 20 |
||
) |
Performs the calculation of axion-photon conversion probability using directly equation (28) from J. Redondo and A. Ringwald, Light shinning through walls. https://arxiv.org/pdf/1011.3741.pdf.
m_gamma will be obtainned from the buffer gas definition. If no buffer gas has been assigned then the medium will be assumed to be vacuum.
Ea in keV, ma in eV
The mgamma and absorption lengths are the ones defined by the gas.
The returned value is given for g_ag = 10^-10 GeV-1
Definition at line 484 of file TRestAxionField.cxx.
Double_t TRestAxionField::GammaTransmissionFWHM | ( | Double_t | step = 0.00001 | ) |
Performs the calculation of the FWHM for the axion-photon conversion probability computed in TRestAxionField::GammaTransmissionProbability
.
Ea in keV, ma in eV, mgamma in eV, deltaL in mm, Bmag in T.
Default value for the parameters are: Bmag = 2.5 T, Lcoh = 10000 mm, Ea = 4 keV eV,
IMPORTANT: In the case that the buffer gas is not defined, this method will return the mass at which the probability reaches half of the maximum vacuum probability.
Scanning towards the right (valid also for vacuum)
Definition at line 766 of file TRestAxionField.cxx.
Double_t TRestAxionField::GammaTransmissionProbability | ( | Double_t | Bmag, |
Double_t | Lcoh, | ||
Double_t | Ea, | ||
Double_t | ma, | ||
Double_t | mg = 0 , |
||
Double_t | absLength = 0 |
||
) |
On top of calculating the gamma transmission probability it will assign new values for the magnetic field (Bmag/T), coherence length (Lcoh/mm) and axion energy (Ea/keV).
Definition at line 355 of file TRestAxionField.cxx.
Double_t TRestAxionField::GammaTransmissionProbability | ( | Double_t | Ea, |
Double_t | ma, | ||
Double_t | mg = 0 , |
||
Double_t | absLength = 0 |
||
) |
Performs the calculation of axion-photon conversion probability using directly equation (11) from van Bibber, Phys Rev D Part Fields. 1989.
If m_gamma (mg) is not given as an argument, i.e. it is equal to zero, then m_gamma will be obtainned from the buffer gas definition. If no buffer gas has been assigned then the medium will be assumed to be vacuum.
Ea in keV, ma in eV, mgamma in eV, Lcoh in mm, Bmag in T
mg in eV, absLength in cm-1
The returned value is given for g_ag = 10^-10 GeV-1
Definition at line 294 of file TRestAxionField.cxx.
Double_t TRestAxionField::GammaTransmissionProbability | ( | std::vector< Double_t > | Bmag, |
Double_t | deltaL, | ||
Double_t | Ea, | ||
Double_t | ma, | ||
Double_t | mg = 0 , |
||
Double_t | absLength = 0 |
||
) |
Performs the calculation of axion-photon conversion probability using directly equation (28) from J. Redondo and A. Ringwald, Light shinning through walls. https://arxiv.org/pdf/1011.3741.pdf.
If m_gamma (mg) is not given as an argument, i.e. it is equal to zero, then m_gamma will be obtainned from the buffer gas definition. If no buffer gas has been assigned then the medium will be assumed to be vacuum.
Ea in keV, ma in eV, mgamma in eV, deltaL in mm, Bmag in T
mg in eV, absLength in cm-1
The returned value is given for g_ag = 10^-10 GeV-1
We integrate following the Midpoint rule method. (Other potential options : Trapezoidal, Simpsons)
Definition at line 383 of file TRestAxionField.cxx.
|
inline |
Definition at line 64 of file TRestAxionField.h.
|
inline |
Definition at line 63 of file TRestAxionField.h.
|
inline |
Definition at line 62 of file TRestAxionField.h.
std::vector< std::pair< Double_t, Double_t > > TRestAxionField::GetMassDensityScanning | ( | std::string | gasName = "He" , |
Double_t | maMax = 0.15 , |
||
Double_t | rampDown = 5.0 |
||
) |
This method determines the proper densities to be used in an axion helioscope experiment in order to achieve a continuous axion mass scan.
The first scanning density is placed where the axion-photon vacuum probability reaches half the value, P_ag(max)/2
. Once the first density, or step, has been obtained, the method calculates the FWHM resonance probability for each density/mass and moves the next scanning axion mass by a step of amplitude FWHM/factor
, where the factor is a value that moves from 2 to 1 as the mass increases, and it falls down with a velocity related with the argument rampDown
, where the factor follows this formula:
factor = TMath::Exp( - ma*rampDown ) + 1;
The method stops when the axion mass is bigger than the maximum axion mass given as an argument, maMax
.
Default arguments: gasName="He", ma_max=0.15 eV, rampDown=5
For additional info see PR: https://github.com/rest-for-physics/axionlib/pull/78
Setting mass-density pair for the first step
Definition at line 816 of file TRestAxionField.cxx.
|
private |
Initialization of TRestAxionField class.
MOVED TO TRestAxionFieldPropagationProcess class faxion = SetComplexReal(1, 0); fAem = SetComplexReal(0, 0);
Definition at line 237 of file TRestAxionField.cxx.
|
inlinestatic |
Integrand used for axion-photon probability integration.
Definition at line 100 of file TRestAxionField.h.
|
inline |
Definition at line 60 of file TRestAxionField.h.
|
inline |
It assigns a gas buffer medium to the calculation.
Definition at line 79 of file TRestAxionField.h.
|
inline |
Definition at line 59 of file TRestAxionField.h.
|
inline |
It enables/disables debug mode.
Definition at line 70 of file TRestAxionField.h.
|
inline |
Definition at line 58 of file TRestAxionField.h.
|
private |
The magnetic field in Teslas (used for constant field formulas)
Definition at line 35 of file TRestAxionField.h.
|
private |
A pointer to the buffer gas definition.
Definition at line 46 of file TRestAxionField.h.
|
private |
Definition at line 32 of file TRestAxionField.h.
|
private |
The energy of the axion in keV.
Definition at line 41 of file TRestAxionField.h.
|
private |
The coherence lenght (in mm) where the magnetic field is defined (for constant field)
Definition at line 38 of file TRestAxionField.h.
|
private |
A pointer to the magnetic field definition.
Definition at line 49 of file TRestAxionField.h.