Berkeley Nuclear Data Software
Public Member Functions | Public Attributes | List of all members
TALYSUtils Class Reference

#include <TALYSUtils.h>

Collaboration diagram for TALYSUtils:
Collaboration graph
[legend]

Public Member Functions

 TALYSUtils ()
 
 TALYSUtils (int a_ZID, int a_AID, std::string a_talysPath)
 
void setTALYSEnergies (vector< double > a_talysEnergies)
 
void setRunECIS (bool a_runECIS)
 perform ECIS-06 calculation - if yes, will also save ECIS results More...
 
int setLDModel (int a_ldmodel)
 set the Level Density model More...
 
void setTALYSPath (std::string a_talysPath)
 
void setZID (int a_ZID)
 
void setAID (int a_AID)
 
void setBashScriptName (std::string a_fileName)
 
int readTALYSParams (std::string a_paramFileName)
 
int setTALYSParams (std::map< std::string, double > a_params)
 
void readTALYSKeywords (std::string a_fileName)
 reads in TALYS keyword map More...
 
vector< std::string > getTALYSKeywords (std::string a_keyword)
 
vector< std::string > getAllTALYSKeywords ()
 returns string of keys of m_keywords More...
 
void clearAdvancedParamList ()
 clears m_params More...
 
int runDefaultTALYSCalc ()
 
int runTALYSCalc ()
 
int getNumAdvancedParams ()
 returns number of parameters that minimizer will try to fit More...
 
map< std::string, double > getAdvancedParamVals ()
 
vector< double > getAdvancedParamLimits (std::string a_paramName)
 
map< std::string, int > setupTALYSMinimizer (ROOT::Math::Minimizer *a_minimizer, bool a_useScaler=true)
 
int processTALYSOutput (int a_verbose=0, bool a_writeGammaLevelMap=false)
 
Eigen::VectorXf readTALYSCrossSection (int a_initialLevel, int a_finalLevel, double a_thresholdAdder=0)
 
Eigen::VectorXf readTALYSCrossSection (std::string a_initialLevel, std::string a_finalLevel, double a_thresholdAdder=0)
 
Eigen::VectorXf drawGammaProduction (std::string a_initialLevel, std::string a_finalLevel)
 draws a cross section More...
 
Eigen::VectorXf drawGammaProduction (int a_initialLevel, int a_finalLevel)
 
TGraph2D * drawDiscreteInelNeutronAngDist (int a_levelNumber, int a_component=0)
 
Eigen::MatrixXf drawElasticAngDist ()
 2D hist incoming energy vs outgoing angle More...
 
TH2F * drawOutEnergyAngle2D (double a_incidentE, int a_component=0)
 
vector< vector< float > > drawInEnerOutEner (double a_angle, bool a_draw=false, int a_component=0)
 
vector< vector< float > > drawOutEner (double a_angle, double a_inEner, int a_component=0)
 outgoing energy distribution at a given incident energy and outgoing angle More...
 
TH2F * drawInEnerOutAngle (double a_outEnergy, int a_component=0)
 
Eigen::VectorXf drawElasticXSec ()
 
vector< float > getOutEnergies (double a_angle)
 returns the energies More...
 
Eigen::VectorXf drawLevelXSec (int a_level, int a_component=2)
 
void plotChannelIntegratedXSec ()
 
vector< float > getBeamEnergies ()
 
vector< std::string > getInitAndFinalLevel (float a_gammaEnergy)
 
float getGammaProdThreshold (float a_gammaEnergy)
 
Eigen::VectorXf plotExclusiveCrossSection (std::string a_rxnType="(n,2n)")
 

Public Attributes

vector< TALYSoutputInfo > m_outputInfo
 holds information parsed from TALYS "output" file More...
 

Constructor & Destructor Documentation

◆ TALYSUtils() [1/2]

TALYSUtils::TALYSUtils ( )

default constructor Z = 1, A = 1, talysPath = "talysOutput/"

◆ TALYSUtils() [2/2]

TALYSUtils::TALYSUtils ( int  a_ZID,
int  a_AID,
std::string  a_talysPath 
)

constructor that takes: (Z,A) of desired isotope path to directory where TALYS will be run e.g. "talysOutput/"

Member Function Documentation

◆ clearAdvancedParamList()

void TALYSUtils::clearAdvancedParamList ( )

clears m_params

◆ drawDiscreteInelNeutronAngDist()

TGraph2D * TALYSUtils::drawDiscreteInelNeutronAngDist ( int  a_levelNumber,
int  a_component = 0 
)

2D hist incoming energy vs outgoing angle a_component = 0 --> total a_component = 1 --> direct a_component = 2 --> compound

◆ drawElasticAngDist()

Eigen::MatrixXf TALYSUtils::drawElasticAngDist ( )

2D hist incoming energy vs outgoing angle

Here is the caller graph for this function:

◆ drawElasticXSec()

Eigen::VectorXf TALYSUtils::drawElasticXSec ( )
Here is the caller graph for this function:

◆ drawGammaProduction() [1/2]

Eigen::VectorXf TALYSUtils::drawGammaProduction ( int  a_initialLevel,
int  a_finalLevel 
)

◆ drawGammaProduction() [2/2]

Eigen::VectorXf TALYSUtils::drawGammaProduction ( std::string  a_initialLevel,
std::string  a_finalLevel 
)

draws a cross section

read in talys production cross section

◆ drawInEnerOutAngle()

TH2F * TALYSUtils::drawInEnerOutAngle ( double  a_outEnergy,
int  a_component = 0 
)

a_component = 0 --> total a_component = 1 --> direct a_component = 4 --> compound

◆ drawInEnerOutEner()

vector< vector< float > > TALYSUtils::drawInEnerOutEner ( double  a_angle,
bool  a_draw = false,
int  a_component = 0 
)

vector {energy in, energy out, cross section} same components as above

Here is the caller graph for this function:

◆ drawLevelXSec()

Eigen::VectorXf TALYSUtils::drawLevelXSec ( int  a_level,
int  a_component = 2 
)

comp 0 = direct comp 1 = compound comp 2 = total

Here is the caller graph for this function:

◆ drawOutEner()

vector< vector< float > > TALYSUtils::drawOutEner ( double  a_angle,
double  a_inEner,
int  a_component = 0 
)

outgoing energy distribution at a given incident energy and outgoing angle

Here is the call graph for this function:

◆ drawOutEnergyAngle2D()

TH2F * TALYSUtils::drawOutEnergyAngle2D ( double  a_incidentE,
int  a_component = 0 
)

a_component = 0 --> total a_component = 1 --> direct a_component = 4 --> compound

◆ getAdvancedParamLimits()

vector< double > TALYSUtils::getAdvancedParamLimits ( std::string  a_paramName)
Here is the caller graph for this function:

◆ getAdvancedParamVals()

map< std::string, double > TALYSUtils::getAdvancedParamVals ( )
Here is the caller graph for this function:

◆ getAllTALYSKeywords()

vector< std::string > TALYSUtils::getAllTALYSKeywords ( )

returns string of keys of m_keywords

◆ getBeamEnergies()

vector< float > TALYSUtils::getBeamEnergies ( )
Here is the caller graph for this function:

◆ getGammaProdThreshold()

float TALYSUtils::getGammaProdThreshold ( float  a_gammaEnergy)

◆ getInitAndFinalLevel()

vector< std::string > TALYSUtils::getInitAndFinalLevel ( float  a_gammaEnergy)

◆ getNumAdvancedParams()

int TALYSUtils::getNumAdvancedParams ( )

returns number of parameters that minimizer will try to fit

Here is the caller graph for this function:

◆ getOutEnergies()

vector< float > TALYSUtils::getOutEnergies ( double  a_angle)

returns the energies

◆ getTALYSKeywords()

vector< std::string > TALYSUtils::getTALYSKeywords ( std::string  a_keyword)

returns string formatted for TALYS input file only handles class 1,2,3,6,7,9 for Z,A keywords, looks at (a_Z,a_A),(a_Z,a_A+1) i.e. a_keyword = T, Z=26, A=56: returns {"T 26 56","T 26 57"} for class 3, only does ground state barrier for parType, only does 'n' for neutron

◆ plotChannelIntegratedXSec()

void TALYSUtils::plotChannelIntegratedXSec ( )

◆ plotExclusiveCrossSection()

Eigen::VectorXf TALYSUtils::plotExclusiveCrossSection ( std::string  a_rxnType = "(n,2n)")

◆ processTALYSOutput()

int TALYSUtils::processTALYSOutput ( int  a_verbose = 0,
bool  a_writeGammaLevelMap = false 
)

parses the "output" file to extract various cross sections W.I.P. -> original focus on secondary neutron energy/angle dists still need to add sections:

  1. : gamma-ray production
  1. : binary non-elastic
  2. : Total particle production
  3. : Residual product production a_verbose = 0 --> no printing as file is parsed a_verbose = 1 --> some printing

total xsecs

end of file

excitation functions 5 -> (n,g) then (n,n') 7 -> gamma ray productions

legendre elastic

skip the table header

loop over all the legendre polynomials factors

inelastic legendre

angle elastic

inelastic angle

Binary reactions to discrete levels and continuum

parameters

total outgoing neutron energy spectrum

DDX mode 3

DDX for outgoing neutron energy -> angle

DDX for outoing neutron angle -> energy

DDX for outgoing neutron energy -> angle

Here is the caller graph for this function:

◆ readTALYSCrossSection() [1/2]

Eigen::VectorXf TALYSUtils::readTALYSCrossSection ( int  a_initialLevel,
int  a_finalLevel,
double  a_thresholdAdder = 0 
)

reads gamma-ray production cross section for a given transition from outputs generated with outgamdis y e.g. for transition from 1st excited to g.s.: a_initialLevel = 1, a_finalLevel = 0 a_thresholdAdder increases reaction threshold i.e. if a_thresholdAdder = 1 [MeV] and reaction threshold is 2 MeV first element in output will be cross section at 3 MeV returns vector of cross section

Here is the caller graph for this function:

◆ readTALYSCrossSection() [2/2]

Eigen::VectorXf TALYSUtils::readTALYSCrossSection ( std::string  a_initialLevel,
std::string  a_finalLevel,
double  a_thresholdAdder = 0 
)

same as above e.g. for transition from 1st excited to g.s.: a_initialLevel = "01", a_finalLevel = "00"

read in talys production cross section

◆ readTALYSKeywords()

void TALYSUtils::readTALYSKeywords ( std::string  a_fileName)

reads in TALYS keyword map

◆ readTALYSParams()

int TALYSUtils::readTALYSParams ( std::string  a_paramFileName)

reads in TALYS parameters from json file "basicSetup" includes parameters there are immutable these parameters are included automatically and don't need to be in the json projectile (n) element (m_ZID) mass (A_ZID) ejectiles (g) outgamdis (y) energy (default: energySpecrum.txt) "advancedSetup" includes parameters that you may want to minimize requires an initial value and the range of values the parameter can take which can be found in the TALYS manual "basicSetup" params are stored in m_basicParams "advancedSetup" params are stored in m_params clears previously stored parameters example of the json is misc/Fe56TalysParams.json

Here is the caller graph for this function:

◆ runDefaultTALYSCalc()

int TALYSUtils::runDefaultTALYSCalc ( )

writes new input file only using 'best' and m_ldmodel expecting m_bashScriptName file in m_talysPath directory that at least contains command to run TALYS i.e.

talys <input> output

◆ runTALYSCalc()

int TALYSUtils::runTALYSCalc ( )

writes input file with default + m_params expecting m_bashScriptName file in m_talysPath directory that at least contains command to run TALYS i.e.

talys <input> output

Here is the caller graph for this function:

◆ setAID()

void TALYSUtils::setAID ( int  a_AID)

◆ setBashScriptName()

void TALYSUtils::setBashScriptName ( std::string  a_fileName)

set the name for the bash script that contains the command: save as m_bashScriptName talys <input> output default is runTalys.sh

◆ setLDModel()

int TALYSUtils::setLDModel ( int  a_ldmodel)

set the Level Density model

◆ setRunECIS()

void TALYSUtils::setRunECIS ( bool  a_runECIS)

perform ECIS-06 calculation - if yes, will also save ECIS results

Here is the caller graph for this function:

◆ setTALYSEnergies()

void TALYSUtils::setTALYSEnergies ( vector< double >  a_talysEnergies)

energies to run TALYS calculations for writes file called "energySpectrum.txt" in m_talysPath directory

Here is the caller graph for this function:

◆ setTALYSParams()

int TALYSUtils::setTALYSParams ( std::map< std::string, double >  a_params)

change the value of a parameter in m_params returns 0 if a parameter in this list in not in m_params

Here is the caller graph for this function:

◆ setTALYSPath()

void TALYSUtils::setTALYSPath ( std::string  a_talysPath)
Here is the caller graph for this function:

◆ setupTALYSMinimizer()

map< std::string, int > TALYSUtils::setupTALYSMinimizer ( ROOT::Math::Minimizer *  a_minimizer,
bool  a_useScaler = true 
)

add "advancedSetup" parameters to minimizer as limitied variable sets: parameter name, initial value, and range set a_useScaler to scale entire cross section by a constant returns map of parameter name and ID

Here is the caller graph for this function:

◆ setZID()

void TALYSUtils::setZID ( int  a_ZID)

Member Data Documentation

◆ m_outputInfo

vector<TALYSoutputInfo> TALYSUtils::m_outputInfo

holds information parsed from TALYS "output" file

The documentation for this class was generated from the following files:
Generated by doxygen 1.9.1