Detailed Description

A metadata class to load tabulated solar axion fluxes. Mass independent.

TRestAxionSolarQCDFlux will use a file in ASCII or binary format to initialize a solar flux table that will describe the solar flux spectrum as a function of the solar radius.

This class may serve to load any generic flux definition that is independent from the axion mass. However, since the design of the class was motivated to reproduce the standard QCD axion flux, therefore the name of the class.

Another scenario arises when the axion or an axion-like particle production mechanism depends on its mass, and we may need to introduce a flux description as the given in the class TRestAxionSolarHiddenPhotonFlux. Both classes are prototyped by a pure base class TRestAxionSolarFlux that defines common methods used to evaluate the flux, and generate Monte-Carlo events inside TRestAxionGeneratorProcess.

Basic use

Once the class has been initialized, the main use of this class will be provided by the method TRestAxionSolarQCDFlux::GetRandomEnergyAndRadius. This method will return a random axion energy and position inside the solar radius following the distributions given by the solar flux tables.

Description of the specific parameters accepted by this metadata class.

fluxDataFile: A table with 100 rows representing the solar ring flux from the center to the corona, and 200 columns representing the flux, measured in cm-2 s-1 keV-1, for the range (0,20)keV in steps of 100eV. The table could be provided in ASCII format, using .dat extension, or it might be a binary table using .N200f extension.
fluxSptFile: A table where each column represents a monochromatic energy. The column contains 101 rows, the first element is the energy of the monochromatic line while the next 100 elements contain the flux, measured in cm-2 s-1, integrated to each solar ring, being the second element the ring in the center of the sun.

Additionally this class will be able to read .flux files that are the original files produced in 3-columns format (inner radius [solar units] / energy [keV] / flux [cm-2 s-1 keV-1]). The .flux files may contain the full information, continuum and spectral components. Those components will be splited into two independent contributions by TRestAxionSolarQCDFlux::ReadFluxFile to be managed internally. Two additional parameters will be required to translate the .flux files into the tables that are understood by this class.

binSize: The energy binning used on the .flux file and inside the histogram used for monochromatic lines identification.
peakSigma: The ratio between the flux provided at the .flux file and the average flux calculated in the peak surroundings. If the flux ratio is higher than this value, the flux at that particular bin will be considered a peak.

Optionally, if we want to consider a different binning on the monochromatic/continuum histogram used internally for the calculation we may specify optionally a new parameter. In that case, fBinSize will be the binning of the internal histogram, while the new parameter will be the binning given inside the .flux file.

fluxBinSize: The bin size used on the .flux table.

Pre-generated solar axion flux tables will be available at the axionlib-data repository. The different RML flux definitions used to load those tables will be found at the fluxes.rml file found at the axionlib-data repository.

Inside a local REST installation, the fluxes.rml file will be found at the REST installation directory, and it will be located automatically by the TRestMetadata::SearchFile method.

A basic RML definition

The following definition integrates an axion-photon component with a continuum spectrum using a Primakoff production model, and a dummy spectrum file that includes two monocrhomatic lines at different solar disk radius positions.

<TRestAxionSolarQCDFlux name="sunPrimakoff" verboseLevel="debug" >
 
    <!-- Common parameters that belong to TRestAxionSolarFlux -->
       <parameter name="couplingType" value="g_ag"/>
       <parameter name="couplingStrength" value="1.e-10"/>
 
    <!-- Specific parameters that belong to TRestAxionSolarQCDFlux -->
       <parameter name="fluxDataFile" value="Primakoff_Gianotti_201904.dat"/>
       <parameter name="fluxSptFile" value="Dummy_Galan_202202.spt"/>
 
</TRestAxionSolarQCDFlux>

Warning: When the flux is loaded manually inside the restRoot interactive shell, or inside a macro or script, after metadata initialization, it is necessary to call the method TRestAxionSolarQCDFlux::LoadTables to trigger the tables initialization.

Performing MonteCarlo tests using pre-loaded tables

In order to test the response of different solar flux definitions we may use the script solarPlot.py found at pipeline/metadata/solarFlux/. This script will generate a number of particles and it will assign to each particle an energy and solar disk location with the help of the method TRestAxionSolarQCDFlux::GetRandomEnergyAndRadius.

python3 solarPlotQCD.py --fluxname LennertHoofABC --N 1000000

By default, it will load the flux definition found at fluxes.rml from the axionlib-data repository, and generate a png image with the resuts from the Monte Carlo execution.

Solar flux distributions MC-generated with TRestAxionSolarQCDFlux.

Reading solar flux tables from `.flux` files

In a similar way we may initialize the class using a .flux file. The .flux described previously will contain a high definition flux table measured in cm-2 s-1 keV-1 as a function of the energy (from 0 to 20keV) and the inner solar radius (from 0 to 1). The binning of this table will be typically between 1eV and 10eV.

The class TRestAxionSolarQCDFlux will be initialized directly using the .flux file provided under the fluxDataFile parameter. During the initialization, the flux will be splitted into two independent flux components. One smooth continuum component integrated in 100 eV steps, and a monochromatic peak components.

In order to help with the identification of peaks we need to define the binSize used in the .flux table and the peakSigma defining the number of sigmas over the average for a bin to be considered a peak.

<TRestAxionSolarQCDFlux name="LennertHoofABC_Flux" verboseLevel="warning" >
    <parameter name="couplingType" value="g_ae"/>
    <parameter name="couplingStrength" value="1.e-13"/>
    <parameter name="fluxDataFile" value="ABC_LennertHoof_202203.flux"/>
 
    <parameter name="binSize" value="10eV" />
    <parameter name="peakSigma" value="10" />
 
    <parameter name="seed" value="137" />
</TRestAxionSolarQCDFlux>

We will be able to load this file as usual, using the following recipe inside restRoot,

TRestAxionSolarQCDFlux *sFlux = new TRestAxionSolarQCDFlux("fluxes.rml", "LennertHoofABC")
sFlux->Initialize()
TCanvas *c = sFlux->DrawSolarFluxes()
c->Print("ABC_FluxTable.png" )

will generate the following figure.

Solar flux distributions generated with TRestAxionSolarQCDFlux::DrawSolarFlux.

Exporting the solar flux tables

On top of that, we will be able to export those tables to the TRestAxionSolarQCDFlux standard format to be used in later occasions.

TRestAxionSolarQCDFlux *sFlux = new TRestAxionSolarQCDFlux("fluxes.rml", "LennertHoofABC")
sFlux->Initialize()
sFlux->ExportTables()

which will produce two files, a binary table .N200f with the continuum flux, and an ASCII table containning the .spt monochromatic lines. The filename root will be extracted from the original .flux file. Optionally we may export the continuum flux to an ASCII file by indicating it at the TRestAxionSolarQCDFlux::ExportTables method call. The files will be placed at the REST user space, at $HOME/.rest/export/ directory.

TODO Implement the method TRestAxionSolarQCDFlux::InitializeSolarTable using a solar model description by TRestAxionSolarModel.

TODO Perhaps it would be interesting to replace fFluxTable for a TH2D

RESTsoft - Software for Rare Event Searches with TPCs

History of developments:

2023-May: Specific methods extracted from TRestAxionSolarFlux Javier Galan

Author: Javier Galan

Definition at line 30 of file TRestAxionSolarQCDFlux.h.

#include <TRestAxionSolarQCDFlux.h>

Inheritance diagram for TRestAxionSolarQCDFlux:

Public Member Functions
	ClassDefOverride (TRestAxionSolarQCDFlux, 1)

virtual TCanvas *	DrawSolarFlux () override
	It draws the contents of a .flux file. This method just receives the. More...

void	ExportTables (Bool_t ascii=false) override
	It will create files with the continuum and spectral flux components to be used in a later ocasion. More...

TH1F *	GetContinuumSpectrum ()
	It builds a histogram with the continuum spectrum component. The flux will be expressed in cm-2 s-1 keV-1. Binned in 100eV steps. More...

TH1F *	GetEnergySpectrum (Double_t m=0) override
	It returns an energy integrated spectrum in cm-2 s-1 keV-1. More...

TH1F *	GetMonochromaticSpectrum ()
	It builds a histogram with the monochromatic spectrum component. The flux will be expressed in cm-2 s-1 eV-1. Binned in 1eV steps. More...

std::pair< Double_t, Double_t >	GetRandomEnergyAndRadius (TVector2 eRange=TVector2(-1, -1), Double_t mass=0) override
	It defines how to generate Monte Carlo energy and radius values to reproduce the flux. More...

Double_t	GetTotalFlux (Double_t mass=0) override
	It returns the total integrated flux at earth in cm-2 s-1. More...

TH1F *	GetTotalSpectrum ()
	It builds a histogram adding the continuum and the monochromatic spectrum component. The flux will be expressed in cm-2 s-1 keV-1. Binned in 1eV steps. More...

void	InitializeSolarTable (TRestAxionSolarModel *model)
	Tables might be loaded using a solar model description by TRestAxionSolarModel. More...

Double_t	IntegrateFluxInRange (TVector2 eRange=TVector2(-1, -1), Double_t mass=0) override
	It returns the integrated flux at earth in cm-2 s-1 for the given energy range. More...

Bool_t	isSolarSpectrumLoaded ()
	It returns true if monochromatic flux spectra was loaded. More...

Bool_t	isSolarTableLoaded ()
	It returns true if continuum flux spectra was loaded. More...

Bool_t	LoadTables () override
	It defines how to read the solar tables at the inhereted class. More...

void	PrintContinuumSolarTable ()
	It prints on screen the table that has been loaded in memory. More...

void	PrintIntegratedRingFlux ()
	It prints on screen the integrated solar flux per solar ring. More...

void	PrintMetadata () override
	Prints on screen the information about the metadata members of TRestAxionSolarQCDFlux. More...

void	PrintMonoChromaticFlux ()
	It prints on screen the spectral lines loaded in memory. More...

	TRestAxionSolarQCDFlux ()
	Default constructor. More...

	TRestAxionSolarQCDFlux (const char *cfgFileName, std::string name="")
	Constructor loading data from a config file. More...

	~TRestAxionSolarQCDFlux ()
	Default destructor. More...

Public Member Functions inherited from TRestAxionSolarFlux
Bool_t	AreTablesLoaded ()

TCanvas *	DrawFluxFile (std::string fname, Double_t binSize=0.01)
	It draws the contents of a .flux file. This method just receives the name of the .flux file and it works stand-alone. More...

virtual TCanvas *	DrawSolarFlux ()
	It draws the contents of a .flux file. This method just receives the. More...

virtual void	ExportTables (Bool_t ascii=false)

virtual TH1F *	GetEnergySpectrum (Double_t m=0)=0
	It returns an energy integrated spectrum in cm-2 s-1 keV-1. More...

Double_t	GetFluxAtEnergy (Double_t energy, Double_t m=0)
	It returns a flux in cm-2 s-1 keV-1 at the energy given by argument. More...

TH1F *	GetFluxHistogram (std::string fname, Double_t binSize=0.01)
	It builds a histogram using the contents of the .flux file given in the argument. More...

virtual std::pair< Double_t, Double_t >	GetRandomEnergyAndRadius (TVector2 eRange=TVector2(-1, -1), Double_t mass=0)=0
	It defines how to generate Monte Carlo energy and radius values to reproduce the flux. More...

virtual Double_t	GetTotalFlux (Double_t mass=0)=0
	It returns the total integrated flux at earth in cm-2 s-1. More...

void	Initialize ()
	It is required in order to load solar flux tables into memory. More...

virtual Double_t	IntegrateFluxInRange (TVector2 eRange=TVector2(-1, -1), Double_t mass=0)=0
	It returns the integrated flux at earth in cm-2 s-1 for the given energy range. More...

void	PrintMetadata ()
	Prints on screen the information about the metadata members of TRestAxionSolarFlux. More...

	~TRestAxionSolarFlux ()
	Default destructor. More...

Public Member Functions inherited from TRestMetadata
void	AddLog (std::string log="", bool print=true)
	Add logs to messageBuffer. More...

void	DoNotStore ()
	If this method is called the metadata information will not be stored in disk. More...

TVector2	Get2DVectorParameterWithUnits (std::string parName, TVector2 defaultValue=TVector2(-1, -1))

TVector3	Get3DVectorParameterWithUnits (std::string parName, TVector3 defaultValue=TVector3(-1, -1, -1))

TString	GetCommit ()
	Returns the REST commit value stored in fCommit. More...

std::string	GetConfigBuffer ()
	Returns the config section of this class. More...

std::string	GetDataMemberValue (std::string memberName)
	Get the value of data member as string. More...

std::vector< std::string >	GetDataMemberValues (std::string memberName, Int_t precision=0)
	Get the value of datamember as a vector of strings. More...

TString	GetDataPath ()
	Returns a std::string with the path used for data storage. More...

Double_t	GetDblParameterWithUnits (std::string parName, Double_t defaultValue=PARAMETER_NOT_FOUND_DBL)
	Gets the value of the parameter name parName, after applying unit conversion. More...

Bool_t	GetError () const
	It returns true if an error was identified by a derived metadata class. More...

TString	GetErrorMessage ()
	Returns a std::string containing the error message. More...

TString	GetLibraryVersion ()
	Returns the REST libraty version stored in fLibraryVersion. More...

TString	GetMainDataPath ()
	Gets a std::string with the path used for data storage. More...

Int_t	GetNumberOfErrors () const

Int_t	GetNumberOfWarnings () const

std::string	GetParameter (std::string parName, TString defaultValue=PARAMETER_NOT_FOUND_STR)
	Returns corresponding REST Metadata parameter from multiple sources. More...

std::string	GetSectionName ()
	Returns the section name of this class, defined at the beginning of fSectionName. More...

TRestStringOutput::REST_Verbose_Level	GetVerboseLevel ()
	returns the verboselevel in type of REST_Verbose_Level enumerator More...

TString	GetVerboseLevelString ()
	returns the verbose level in type of TString More...

TString	GetVersion ()
	Returns the REST version stored in fVersion. More...

Int_t	GetVersionCode ()

UInt_t	GetVersionMajor () const

UInt_t	GetVersionMinor () const

UInt_t	GetVersionPatch () const

Bool_t	GetWarning () const
	It returns true if an error was identified by a derived metadata class. More...

TString	GetWarningMessage ()
	Returns a std::string containing the warning message. More...

virtual void	Initialize ()
	Making default settings. More...

TRestMetadata *	InstantiateChildMetadata (int index, std::string pattern="")
	This method will retrieve a new TRestMetadata instance of a child element of the present TRestMetadata instance based on the `index` given by argument, which defines the element order to be retrieved, 0 for first element found, 1 for the second element found, etc. More...

TRestMetadata *	InstantiateChildMetadata (std::string pattern="", std::string name="")
	This method will retrieve a new TRestMetadata instance of a child element of the present TRestMetadata instance based on the `name` given by argument. More...

Bool_t	isCleanState () const

Bool_t	isOfficialRelease () const

Int_t	LoadConfigFromBuffer ()
	Initialize data from a string element buffer. More...

Int_t	LoadConfigFromElement (TiXmlElement eSectional, TiXmlElement eGlobal, std::map< std::string, std::string > envs={})
	Main starter method. More...

Int_t	LoadConfigFromFile (const std::string &configFilename, const std::string &sectionName="")
	Give the file name, find out the corresponding section. Then call the main starter. More...

virtual void	Merge (const TRestMetadata &)

TRestMetadata &	operator= (const TRestMetadata &)

void	Print ()
	Implementing TObject::Print() method. More...

void	PrintConfigBuffer ()
	Print the config xml section stored in the class. More...

void	PrintMessageBuffer ()
	Print the buffered message. More...

virtual void	PrintMetadata ()
	Implemented it in the derived metadata class to print out specific metadata information. More...

void	PrintTimeStamp (Double_t timeStamp)
	Print the current time on local machine. More...

void	SetConfigFile (std::string configFilename)
	set config file path from external More...

void	SetError (std::string message="", bool print=true, int maxPrint=5)
	A metadata class may use this method to signal that something went wrong. More...

void	SetHostmgr (TRestManager *m)
	Set the host manager for this class. More...

void	SetSectionName (std::string sName)
	set the section name, clear the section content More...

void	SetVerboseLevel (TRestStringOutput::REST_Verbose_Level v)
	sets the verbose level More...

void	SetWarning (std::string message="", bool print=true, int maxPrint=5)
	A metadata class may use this method to signal that something went wrong. More...

void	Store ()
	If this method is called the metadata information will be stored in disk. More...

	TRestMetadata (const TRestMetadata &)

virtual void	UpdateMetadataMembers ()
	Method to allow implementation of specific metadata members updates at inherited classes. More...

virtual Int_t	Write (const char *name=nullptr, Int_t option=0, Int_t bufsize=0)
	overwriting the write() method with fStore considered More...

void	WriteConfigBuffer (std::string fName)
	Writes the config buffer to a file in append mode. More...

	~TRestMetadata ()
	TRestMetadata default destructor. More...

Private Member Functions
void	IntegrateSolarFluxes ()
	A helper method to initialize the internal class data members with the integrated flux for each solar ring. It will be called by TRestAxionSolarQCDFlux::Initialize. More...

void	LoadContinuumFluxTable ()
	A helper method to load the data file containning continuum spectra as a function of the solar radius. It will be called by TRestAxionSolarQCDFlux::Initialize. More...

void	LoadMonoChromaticFluxTable ()
	A helper method to load the data file containning monochromatic spectral lines as a function of the solar radius. It will be called by TRestAxionSolarQCDFlux::Initialize. More...

void	ReadFluxFile ()
	It loads a .flux file. It will split continuum and monochromatic peaks, loading both internal flux tables. More...

Private Attributes
Double_t	fBinSize = 0
	It will be used when loading `.flux` files to define the input file energy binsize in eV. More...

TH1F *	fContinuumHist = nullptr
	A pointer to the continuum spectrum histogram. More...

std::string	fFluxDataFile = ""
	The filename containning the solar flux table with continuum spectrum. More...

std::vector< Double_t >	fFluxLineIntegrals
	Accumulative integrated solar flux for each monochromatic energy (renormalized to unity) More...

std::map< Double_t, TH1F * >	fFluxLines
	The tabulated solar flux in cm-2 s-1 for a number of monochromatic energies versus solar radius. More...

Double_t	fFluxRatio = 0
	The ratio between monochromatic and total flux. More...

std::string	fFluxSptFile = ""
	The filename containning the solar flux spectra for monochromatic spectrum. More...

std::vector< TH1F * >	fFluxTable
	The tabulated solar flux continuum spectra TH1F(200,0,20)keV in cm-2 s-1 keV-1 versus solar radius. More...

std::vector< Double_t >	fFluxTableIntegrals
	Accumulative integrated solar flux for each solar ring for continuum spectrum (renormalized to unity) More...

TH1F *	fMonoHist = nullptr
	A pointer to the monochromatic spectrum histogram. More...

Double_t	fPeakSigma = 0
	It will be used when loading `.flux` files to define the threshold for peak identification. More...

Double_t	fTotalContinuumFlux = 0
	Total solar flux for monochromatic contributions. More...

TH1F *	fTotalHist = nullptr
	A pointer to the superposed monochromatic and continuum spectrum histogram. More...

Double_t	fTotalMonochromaticFlux = 0
	Total solar flux for monochromatic contributions. More...

Additional Inherited Members
Protected Member Functions inherited from TRestAxionSolarFlux
virtual Bool_t	LoadTables ()=0
	It defines how to read the solar tables at the inhereted class. More...

	TRestAxionSolarFlux ()
	Default constructor. More...

	TRestAxionSolarFlux (const char *cfgFileName, std::string name="")
	Constructor loading data from a config file. More...

Protected Member Functions inherited from TRestMetadata
std::string	ElementToString (TiXmlElement *ele)
	Convert an TiXmlElement object to string. More...

TVector2	Get2DVectorParameterWithUnits (std::string parName, TiXmlElement *e, TVector2 defaultValue=TVector2(-1, -1))

TVector3	Get3DVectorParameterWithUnits (std::string parName, TiXmlElement *e, TVector3 defaultValue=TVector3(-1, -1, -1))

Double_t	GetDblParameterWithUnits (std::string parName, TiXmlElement *e, Double_t defaultVal=PARAMETER_NOT_FOUND_DBL)

TiXmlElement *	GetElement (std::string eleDeclare, TiXmlElement *e=nullptr)
	Get an xml element from a given parent element, according to its declaration. More...

TiXmlElement *	GetElementFromFile (std::string configFilename, std::string NameOrDecalre="")
	Open an xml encoded file and find its element. More...

TiXmlElement *	GetElementWithName (std::string eleDeclare, std::string eleName)
	Get an xml element from the default location, according to its declaration and its field "name". More...

TiXmlElement *	GetElementWithName (std::string eleDeclare, std::string eleName, TiXmlElement *e)
	Get an xml element from a given parent element, according to its declaration and its field "name". More...

std::string	GetFieldValue (std::string fieldName, std::string definition, size_t fromPosition=0)
	Gets field value in an xml element string by parsing it as TiXmlElement. More...

std::string	GetFieldValue (std::string parName, TiXmlElement *e)
	Returns the field value of an xml element which has the specified name. More...

std::string	GetKEYDefinition (std::string keyName)
	Gets the first key definition for keyName found inside buffer starting at fromPosition. More...

std::string	GetKEYDefinition (std::string keyName, size_t &Position)

std::string	GetKEYDefinition (std::string keyName, size_t &Position, std::string buffer)

std::string	GetKEYDefinition (std::string keyName, std::string buffer)

std::string	GetKEYStructure (std::string keyName)
	Gets the first key structure for keyName found inside buffer after fromPosition. More...

std::string	GetKEYStructure (std::string keyName, size_t &Position)

std::string	GetKEYStructure (std::string keyName, size_t &Position, std::string buffer)

std::string	GetKEYStructure (std::string keyName, size_t &Position, TiXmlElement *ele)

std::string	GetKEYStructure (std::string keyName, std::string buffer)

TiXmlElement *	GetNextElement (TiXmlElement *e)
	Get the next sibling xml element of this element, with same eleDeclare. More...

std::string	GetParameter (std::string parName, size_t &pos, std::string inputString)
	Returns the value for the parameter name parName found in inputString. More...

std::string	GetParameter (std::string parName, TiXmlElement *e, TString defaultValue=PARAMETER_NOT_FOUND_STR)
	Returns the value for the parameter named parName in the given section. More...

std::pair< std::string, std::string >	GetParameterAndUnits (std::string parname, TiXmlElement *e=nullptr)
	Returns the unit string of the given parameter of the given xml section. More...

std::map< std::string, std::string >	GetParametersList ()
	It retrieves a map of all parameter:value found in the metadata class. More...

TString	GetSearchPath ()

virtual void	InitFromConfigFile ()
	To make settings from rml file. This method must be implemented in the derived class. More...

virtual void	InitFromRootFile ()
	Method called after the object is retrieved from root file. More...

virtual Int_t	LoadSectionMetadata ()
	This method does some preparation of xml section. More...

void	ReadAllParameters ()
	Reflection methods, Set value of a datamember in class according to TRestMetadata::fElement. More...

void	ReadParametersList (std::map< std::string, std::string > &list)
	It reads a parameter list and associates it to its corresponding metadata member. par0 --> fPar0. More...

std::string	ReplaceConstants (const std::string buffer)
	Identifies "constants" in the input buffer, and replace them with corresponding value. More...

std::string	ReplaceVariables (const std::string buffer)
	Identifies environmental variable replacing marks in the input buffer, and replace them with corresponding value. More...

void	ReSetVersion ()
	Resets the version of TRestRun to REST_RELEASE. Only TRestRun is allowed to update version. More...

std::string	SearchFile (std::string filename)
	Search files in current directory and directories specified in "searchPath" section. More...

void	SetLibraryVersion (TString version)
	Set the library version of this metadata class. More...

TiXmlElement *	StringToElement (std::string definition)
	Parsing a string into TiXmlElement object. More...

	TRestMetadata ()
	TRestMetadata default constructor. More...

	TRestMetadata (const char *configFilename)
	constructor More...

void	UnSetVersion ()
	Resets the version of TRestRun to -1, in case the file is old REST file. Only TRestRun is allowed to update version. More...

Protected Attributes inherited from TRestAxionSolarFlux
TCanvas *	fCanvas = nullptr
	A canvas pointer for drawing. More...

TRandom3 *	fRandom = nullptr
	Random number generator. More...

Bool_t	fTablesLoaded = false
	A metadata member to control if this class has been initialized. More...

Protected Attributes inherited from TRestMetadata
std::string	configBuffer
	The buffer where the corresponding metadata section is stored. Filled only during Write() More...

std::string	fConfigFileName
	Full name of the rml file. More...

std::map< std::string, std::string >	fConstants
	Saving a list of rml constants. name-value std::pair. Constants are temporary for this class only. More...

TiXmlElement *	fElement
	Saving the sectional element together with global element. More...

TiXmlElement *	fElementGlobal
	Saving the global element, to be passed to the resident class, if necessary. More...

Bool_t	fError = false
	It can be used as a way to identify that something went wrong using SetError method. More...

TString	fErrorMessage = ""
	A std::string to store an optional error message through method SetError. More...

TRestManager *	fHostmgr
	All metadata classes can be initialized and managed by TRestManager. More...

Int_t	fNErrors = 0
	It counts the number of errors notified. More...

Int_t	fNWarnings = 0
	It counts the number of warnings notified. More...

std::string	fSectionName
	Section name given in the constructor of the derived metadata class. More...

Bool_t	fStore
	This variable is used to determine if the metadata structure should be stored in the ROOT file. More...

std::map< std::string, std::string >	fVariables
	Saving a list of rml variables. name-value std::pair. More...

TRestStringOutput::REST_Verbose_Level	fVerboseLevel
	Verbose level used to print debug info. More...

Bool_t	fWarning = false
	It can be used as a way to identify that something went wrong using SetWarning method. More...

TString	fWarningMessage = ""
	It can be used as a way to identify that something went wrong using SetWarning method. More...

std::string	messageBuffer
	The buffer to store the output message through TRestStringOutput in this class. More...

endl_t	RESTendl
	Termination flag object for TRestStringOutput. More...

Constructor & Destructor Documentation

◆ TRestAxionSolarQCDFlux() [1/2]

TRestAxionSolarQCDFlux::TRestAxionSolarQCDFlux ( )

Default constructor.

Definition at line 220 of file TRestAxionSolarQCDFlux.cxx.

◆ TRestAxionSolarQCDFlux() [2/2]

TRestAxionSolarQCDFlux::TRestAxionSolarQCDFlux	(	const char *	cfgFileName,
		std::string	name = `""`
	)

Constructor loading data from a config file.

If no configuration path is defined using TRestMetadata::SetConfigFilePath the path to the config file must be specified using full path, absolute or relative.

The default behaviour is that the config file must be specified with full path, absolute or relative.

Parameters

cfgFileName	A const char* giving the path to an RML file.
name	The name of the specific metadata. It will be used to find the corresponding TRestAxionMagneticField section inside the RML.

Definition at line 236 of file TRestAxionSolarQCDFlux.cxx.

◆ ~TRestAxionSolarQCDFlux()

TRestAxionSolarQCDFlux::~TRestAxionSolarQCDFlux ( )

Default destructor.

Definition at line 246 of file TRestAxionSolarQCDFlux.cxx.

Member Function Documentation

◆ DrawSolarFlux()

TCanvas * TRestAxionSolarQCDFlux::DrawSolarFlux ( )

overridevirtual

It draws the contents of a .flux file. This method just receives the.

Reimplemented from TRestAxionSolarFlux.

Definition at line 569 of file TRestAxionSolarQCDFlux.cxx.

◆ ExportTables()

void TRestAxionSolarQCDFlux::ExportTables ( Bool_t ascii = false )

overridevirtual

It will create files with the continuum and spectral flux components to be used in a later ocasion.

Reimplemented from TRestAxionSolarFlux.

Definition at line 766 of file TRestAxionSolarQCDFlux.cxx.

◆ GetContinuumSpectrum()

TH1F * TRestAxionSolarQCDFlux::GetContinuumSpectrum ( )

It builds a histogram with the continuum spectrum component. The flux will be expressed in cm-2 s-1 keV-1. Binned in 100eV steps.

Definition at line 482 of file TRestAxionSolarQCDFlux.cxx.

◆ GetEnergySpectrum()

TH1F * TRestAxionSolarQCDFlux::GetEnergySpectrum ( Double_t m = 0 )

inlineoverridevirtual

It returns an energy integrated spectrum in cm-2 s-1 keV-1.

Implements TRestAxionSolarFlux.

Definition at line 102 of file TRestAxionSolarQCDFlux.h.

◆ GetMonochromaticSpectrum()

TH1F * TRestAxionSolarQCDFlux::GetMonochromaticSpectrum ( )

It builds a histogram with the monochromatic spectrum component. The flux will be expressed in cm-2 s-1 eV-1. Binned in 1eV steps.

Definition at line 508 of file TRestAxionSolarQCDFlux.cxx.

◆ GetRandomEnergyAndRadius()

std::pair< Double_t, Double_t > TRestAxionSolarQCDFlux::GetRandomEnergyAndRadius	(	TVector2	eRange = `TVector2(-1, -1)`,
		Double_t	mass = `0`
	)

overridevirtual

It defines how to generate Monte Carlo energy and radius values to reproduce the flux.

It returns a random solar radius position and energy according to the flux distributions defined inside the solar tables loaded in the class.

Implements TRestAxionSolarFlux.

Definition at line 665 of file TRestAxionSolarQCDFlux.cxx.

◆ GetTotalFlux()

Double_t TRestAxionSolarQCDFlux::GetTotalFlux ( Double_t mass = 0 )

inlineoverridevirtual

It returns the total integrated flux at earth in cm-2 s-1.

Implements TRestAxionSolarFlux.

Definition at line 97 of file TRestAxionSolarQCDFlux.h.

◆ GetTotalSpectrum()

TH1F * TRestAxionSolarQCDFlux::GetTotalSpectrum ( )

It builds a histogram adding the continuum and the monochromatic spectrum component. The flux will be expressed in cm-2 s-1 keV-1. Binned in 1eV steps.

Definition at line 535 of file TRestAxionSolarQCDFlux.cxx.

◆ InitializeSolarTable()

void TRestAxionSolarQCDFlux::InitializeSolarTable ( TRestAxionSolarModel * model )

inline

Tables might be loaded using a solar model description by TRestAxionSolarModel.

Definition at line 111 of file TRestAxionSolarQCDFlux.h.

◆ IntegrateFluxInRange()

Double_t TRestAxionSolarQCDFlux::IntegrateFluxInRange	(	TVector2	eRange = `TVector2(-1, -1)`,
		Double_t	mass = `0`
	)

overridevirtual

It returns the integrated flux at earth in cm-2 s-1 for the given energy range.

Implements TRestAxionSolarFlux.

Definition at line 641 of file TRestAxionSolarQCDFlux.cxx.

◆ IntegrateSolarFluxes()

void TRestAxionSolarQCDFlux::IntegrateSolarFluxes ( )

private

A helper method to initialize the internal class data members with the integrated flux for each solar ring. It will be called by TRestAxionSolarQCDFlux::Initialize.

Definition at line 614 of file TRestAxionSolarQCDFlux.cxx.

◆ isSolarSpectrumLoaded()

Bool_t TRestAxionSolarQCDFlux::isSolarSpectrumLoaded ( )

inline

It returns true if monochromatic flux spectra was loaded.

Definition at line 84 of file TRestAxionSolarQCDFlux.h.

◆ isSolarTableLoaded()

Bool_t TRestAxionSolarQCDFlux::isSolarTableLoaded ( )

inline

It returns true if continuum flux spectra was loaded.

Definition at line 81 of file TRestAxionSolarQCDFlux.h.

◆ LoadContinuumFluxTable()

void TRestAxionSolarQCDFlux::LoadContinuumFluxTable ( )

private

A helper method to load the data file containning continuum spectra as a function of the solar radius. It will be called by TRestAxionSolarQCDFlux::Initialize.

Definition at line 275 of file TRestAxionSolarQCDFlux.cxx.

◆ LoadMonoChromaticFluxTable()

void TRestAxionSolarQCDFlux::LoadMonoChromaticFluxTable ( )

private

A helper method to load the data file containning monochromatic spectral lines as a function of the solar radius. It will be called by TRestAxionSolarQCDFlux::Initialize.

Definition at line 326 of file TRestAxionSolarQCDFlux.cxx.

◆ LoadTables()

Bool_t TRestAxionSolarQCDFlux::LoadTables ( )

overridevirtual

It defines how to read the solar tables at the inhereted class.

It will load the tables in memory by using the filename information provided inside the metadata members.

Implements TRestAxionSolarFlux.

Definition at line 252 of file TRestAxionSolarQCDFlux.cxx.

◆ PrintContinuumSolarTable()

void TRestAxionSolarQCDFlux::PrintContinuumSolarTable ( )

It prints on screen the table that has been loaded in memory.

Definition at line 700 of file TRestAxionSolarQCDFlux.cxx.

◆ PrintIntegratedRingFlux()

void TRestAxionSolarQCDFlux::PrintIntegratedRingFlux ( )

It prints on screen the integrated solar flux per solar ring.

Definition at line 715 of file TRestAxionSolarQCDFlux.cxx.

◆ PrintMetadata()

void TRestAxionSolarQCDFlux::PrintMetadata ( )

overridevirtual

Prints on screen the information about the metadata members of TRestAxionSolarQCDFlux.

Reimplemented from TRestAxionSolarFlux.

Definition at line 743 of file TRestAxionSolarQCDFlux.cxx.

◆ PrintMonoChromaticFlux()

void TRestAxionSolarQCDFlux::PrintMonoChromaticFlux ( )

It prints on screen the spectral lines loaded in memory.

Definition at line 728 of file TRestAxionSolarQCDFlux.cxx.

◆ ReadFluxFile()

void TRestAxionSolarQCDFlux::ReadFluxFile ( )

private

It loads a .flux file. It will split continuum and monochromatic peaks, loading both internal flux tables.

Definition at line 373 of file TRestAxionSolarQCDFlux.cxx.

Field Documentation

◆ fBinSize

Double_t TRestAxionSolarQCDFlux::fBinSize = 0

private

It will be used when loading .flux files to define the input file energy binsize in eV.

Definition at line 39 of file TRestAxionSolarQCDFlux.h.

◆ fContinuumHist

TH1F* TRestAxionSolarQCDFlux::fContinuumHist = nullptr

private

A pointer to the continuum spectrum histogram.

Definition at line 66 of file TRestAxionSolarQCDFlux.h.

◆ fFluxDataFile

std::string TRestAxionSolarQCDFlux::fFluxDataFile = ""

private

The filename containning the solar flux table with continuum spectrum.

Definition at line 33 of file TRestAxionSolarQCDFlux.h.

◆ fFluxLineIntegrals

std::vector<Double_t> TRestAxionSolarQCDFlux::fFluxLineIntegrals

private

Accumulative integrated solar flux for each monochromatic energy (renormalized to unity)

Definition at line 54 of file TRestAxionSolarQCDFlux.h.

◆ fFluxLines

std::map<Double_t, TH1F*> TRestAxionSolarQCDFlux::fFluxLines

private

The tabulated solar flux in cm-2 s-1 for a number of monochromatic energies versus solar radius.

Definition at line 48 of file TRestAxionSolarQCDFlux.h.

◆ fFluxRatio

Double_t TRestAxionSolarQCDFlux::fFluxRatio = 0

private

The ratio between monochromatic and total flux.

Definition at line 63 of file TRestAxionSolarQCDFlux.h.

◆ fFluxSptFile

std::string TRestAxionSolarQCDFlux::fFluxSptFile = ""

private

The filename containning the solar flux spectra for monochromatic spectrum.

Definition at line 36 of file TRestAxionSolarQCDFlux.h.

◆ fFluxTable

std::vector<TH1F*> TRestAxionSolarQCDFlux::fFluxTable

private

The tabulated solar flux continuum spectra TH1F(200,0,20)keV in cm-2 s-1 keV-1 versus solar radius.

Definition at line 45 of file TRestAxionSolarQCDFlux.h.

◆ fFluxTableIntegrals

std::vector<Double_t> TRestAxionSolarQCDFlux::fFluxTableIntegrals

private

Accumulative integrated solar flux for each solar ring for continuum spectrum (renormalized to unity)

Definition at line 51 of file TRestAxionSolarQCDFlux.h.

◆ fMonoHist

TH1F* TRestAxionSolarQCDFlux::fMonoHist = nullptr

private

A pointer to the monochromatic spectrum histogram.

Definition at line 69 of file TRestAxionSolarQCDFlux.h.

◆ fPeakSigma

Double_t TRestAxionSolarQCDFlux::fPeakSigma = 0

private

It will be used when loading .flux files to define the threshold for peak identification.

Definition at line 42 of file TRestAxionSolarQCDFlux.h.

◆ fTotalContinuumFlux

Double_t TRestAxionSolarQCDFlux::fTotalContinuumFlux = 0

private

Total solar flux for monochromatic contributions.

Definition at line 60 of file TRestAxionSolarQCDFlux.h.

◆ fTotalHist

TH1F* TRestAxionSolarQCDFlux::fTotalHist = nullptr

private

A pointer to the superposed monochromatic and continuum spectrum histogram.

Definition at line 72 of file TRestAxionSolarQCDFlux.h.

◆ fTotalMonochromaticFlux

Double_t TRestAxionSolarQCDFlux::fTotalMonochromaticFlux = 0

private

Total solar flux for monochromatic contributions.

Definition at line 57 of file TRestAxionSolarQCDFlux.h.

The documentation for this class was generated from the following files:

Detailed Description

Basic use

A basic RML definition

Performing MonteCarlo tests using pre-loaded tables

Reading solar flux tables from .flux files

Exporting the solar flux tables

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ TRestAxionSolarQCDFlux() [1/2]

◆ TRestAxionSolarQCDFlux() [2/2]

◆ ~TRestAxionSolarQCDFlux()

Member Function Documentation

◆ DrawSolarFlux()

◆ ExportTables()

◆ GetContinuumSpectrum()

◆ GetEnergySpectrum()

◆ GetMonochromaticSpectrum()

◆ GetRandomEnergyAndRadius()

◆ GetTotalFlux()

◆ GetTotalSpectrum()

◆ InitializeSolarTable()

◆ IntegrateFluxInRange()

◆ IntegrateSolarFluxes()

◆ isSolarSpectrumLoaded()

◆ isSolarTableLoaded()

◆ LoadContinuumFluxTable()

◆ LoadMonoChromaticFluxTable()

◆ LoadTables()

◆ PrintContinuumSolarTable()

◆ PrintIntegratedRingFlux()

◆ PrintMetadata()

◆ PrintMonoChromaticFlux()

◆ ReadFluxFile()

Field Documentation

◆ fBinSize

◆ fContinuumHist

◆ fFluxDataFile

◆ fFluxLineIntegrals

◆ fFluxLines

◆ fFluxRatio

◆ fFluxSptFile

◆ fFluxTable

◆ fFluxTableIntegrals

◆ fMonoHist

◆ fPeakSigma

◆ fTotalContinuumFlux

◆ fTotalHist

◆ fTotalMonochromaticFlux

Reading solar flux tables from `.flux` files