Isotope Explorer user's manual

Isotope Explorer user's manual

The script language

This section describes the script language in build process. The information is only needed if you want to produce your own charts.

General rules

Any text between // and the end of a line is treated as comment
Case is ignored

Commands

BandOut

Output band labels for all selected data sets

ColorAccordingTo(parameter)

Define which parameter should be used for the color coding of the chart

Database(database)

where database is TOI or Moller or databasename

The TOI and Moller databases contain ground state properties of the nuclides. The TOI database contains experimental data, while the Moller database is calculated from theory. The database file name are toi.dat and moller.dat. TOI is the default database
The user can also supply his own database. The format of the database is very simple (see toi.dat for example). The file extension must be .dat and the file name appears in place of databasename

Def(param=...)

See section on predefined functions below

Display(Draw or Table or Plot)

A left click on a nuclide box will show a list of data. If the list is datasets or bands, the user can select a band or a data set by clicking. A display of the given type will be opened

DSOut

Output data set names for selected datasets

EnsdfDB(ENSDF or PENSDF or SENSDF or EHSDF or EDDF)
In order to search for e.g. levels or gammas, the user needs to use an ENSDF database. The default is ENSDF, but there are five to choose from:

ENSDF is the Evaluated Nuclear Structure Data File. It contains complete level information for each nuclide.
PENSDF database is derived from ENSDF by grouping together decay datasets with the same parent
SENSDF database contains data for nuclides with known superdeformed bands
EHSDF database contains only adopted datasets from ENSDF and most of the comments have been removed in order to make the files smaller. There are some high-spin data added to this database
EDDF database contains only decay datasets from ENSDF

FindGamma(...)

Find gammas that fulfil certain criteria
Example:
FindGamma(e=250-260) means find Gamma Energy between 250 and 260 keV
Note: The command GammaOut(quantity) is needed to generate output

FindLevel(...)

Find levels that fulfil certain criteria.
Examples:
FindLevel(e=300-350) means find levels with an excitation energy in the range 300-350 keV
FindLevel(t=1-1.e5) means find levels with half life in the range 1-10⁵ seconds
FindLevel(spin={0+,2+,1-}) means find the lowest 1- level which is higher than a 2+ which is higher than a 0+ *** This works, but it is a bit inconsistent. In FindLevel you require the term spin, but in LevelOut the term JPI. Choose one or the other ***
FindLevel(e=0,t=1-1e7) means find ground states with half lives in the range 1-10⁷ seconds
Note: The command LevelOut(quantity) is needed to generate output

FindOut(quantity)

Extract quantity from a chart database (.DAT file) and store it in the chart. If the command OutToFile(filename) is given, the data is stored in the file filename

GammaOut(E,ELI,ELF,JPII,JPIF,RI)

Output the given quantities to the chart
E is gamma-ray energy
ELI, ELF is initial, final level energy
JPII, JPIF is initial, final level spin/parity.
RI is relative intensity
If uncertainty is required, add D to the key, e.g. ED, RID *** RID picks ut RI again not Delta(RI)***

KeepLines

***What is this????

LevelOut(E,JPI,J(J+1),T)

Output the given quantities to the chart
E is level energy
JPI is spin/parity of the level
J(J+1) is the spin squared
T is the half life
If uncertainty is required, add D to the key, e.g. ED, TD

Nucleus(filter1,filter2...)

where filter can be a=50-, a=50-100, a=-50, z=152, even-a, odd-z, etc. If this command is not given, all nuclides in the Chart section (or if the Chart section is absent, the Data file) will be processed.

OutToFile(FileName)

When a chart is built, it is possible to redirect the output to a file with the name given inside the brackets. If this command is missing, the output will be stored in the chart (memory). To save the output to a disk, select Chart, Save as... to save the chart file.

SelectBand(All or Yrast or GS or SD)

The default database is Ensdf. The default data set is the first data set of the nuclide. ** Now yrast is defined only for even-even and positive parity even spin such as 0+,2+,4+.... ***this is temporary, I assume. Yrast is defined for any spin. There are yrast bands defined for other nuclei, but they are very strange. Try to define the yrast band for all nuclei.***

SelectDS(All or Decay or B-Decay or EC-Decay or NGThermal)

Select data sets of the specified type. This command must be followed by a DSOut command.*** B-Decay does not work. What is it supposed to mean? B--decay or just B-decay?*** ***EC-decay and NGThermal do not work either***

SlopeOut

***What is this????

Examples

For A=100-200, with database ENSDF, find gammas with energy in the range 1200-1250 keV and output gamma-ray energy, initial level energy, spin/parity of initial and final level to the chart and to the file test.out:


nucleus(a=100-200)

database(ensdf)

outtofile(test.out)

findgamma(e=1200-1250)

gammaout(e,eli,jpii,jpif)



****This is what comes out in test.out, wrong, isn't it?

100Y  (Z=39): 2

     1241.66     1340.74                        

     1240.12     1412.14                    (3+)

100Y 

 2

     1241.66     1340.74                        

     1240.12     1412.14                    (3+)

100Zr (Z=40): 1

     1228.99     1441.44      (1,2+)          2+

100Zr

 1

     1228.99     1441.44      (1,2+)          2+

With database SENSDF, for A=80-200, select superdeformed bands and output the labels to the chart. When a band is selected, a band plot is shown:
```
DataBase(SEnsdf)

Nucleus(A=80-200)

SelectBand(sd)

BandOut

Display(Plot)
```

Predefined functions

Functions

Even(argument), Odd(argument): Return 1 if argument is an even integer, odd integer. Otherwise return 0
Sqrt(argument), Exp(argument), Log(argument), Log10(argument): Return values as the corresponding Fortran functions
Sin(argument), Cos(argument), Tan(argument): Return values as the corresponding Fortran functions
Asin(argument), Acos(argument), Atan(argument), Atan2(argument1,argument2): Return values as the corresponding Fortran functions
Pi(): Returns 3.141592653
NuclA(), NuclZ(), NuclN(): Return A, Z, N for the nuclide
If(condition,YesValue,NoValue): Returns YesValue if condition is not equal to 0. Otherwise returns NoValue
Val(parameter):: Returns value of parameter if it is defined. Otherwise returns 0
IsQM(parameter): Returns 1 if value of parameter contains a '?'. Otherwise returns 0
IsGT(parameter), IsLT(Parameter), IsAP(Parameter), IsEQ(Parameter): Return 1 if value of parameter defined as '>', '<', '~', '='. Otherwise return 0
IsDef(Parameter): Returns 1 if parameter is defined in this nuclide. Otherwise returns 0
EQ(argument1,argument2): If argument1 is equal to argument2 return 1. Otherwise return 0
GE(argument1,argument2): If argument1 is greater than or equal to argument2 return 1. Otherwise return 0
LE(argument1,argument2): If argument1 is smaller than or equal to argument2 return 1. Otherwise return 0
GT(argument1,argument2): If argument1 is greater than argument2 return 1. Otherwise return 0
LT(argument1,argument2): If argument1 is smaller than argument2 return 1. Otherwise return 0

Rules

If a parameter's name contains non-alphanumeric characters it must be embraced by double quotes ("a parameter name","qb-").
Optionally the parameter can be preceded by a N,Z offset or a database name so that the value of the parameter is taken from a neighboring nuclide or another database. For example, if the parameter mass is defined in the database Moller,
"[[moller]][-1,0]mass"
is the mass for the N-1,Z neighbor nuclide from the Moller database
Spaces between operators are ignored
If a line ends with the character '\', the line will be merged with the following line

Examples

If %A is uncertain, set it to 500:


  Def(myvar1=if(isqm("%A"),500,"%a"))

  FindOut(myvar1)

Calculate the nuclide mass from the Weizäcker mass formula (mass0). Get the experimental mass (mass) from the TOI database:


  def(Z=nuclZ())

  def(A=nuclA())

  def(N=nuclN())

  def(av=14.1)

  def(as=13)

  def(ac=0.595)

  def(aa=19)

  def(ap=33.5)

  def(Delta=ap/A^(3/4)*If(Even(A),If(Even(Z),-1,1),0))

  def(mass0=Z*938.767+N*939.549-av*A +as*A^(2/3)+ac*Z^2/A^(1/3)+aa*(A-2*Z)^2/A+Delta)

  FindOut(mass0)

  FindOut(mass)

Chart user-interface

The chart mode can be used as a user interface to select data to be displayed. To activate the chart user-interface, a section like the following should be included in the build section of the chart file (comments after //):


  BUILD:

  DataBase(Ensdf)      //define database

  SelectDS(decay)      //select decay data sets

  dsout                //output list of data sets to chart

  Display(level)       //display level plot of selected data set

  END BUILD

The Data files

Chart data can also be made available to Isotope Explorer as data files with the extension .DAT. See the file toi.dat for the required format.

Comments to: F Chu (program), P Ekström (manual)