2024-10-29T12:45:54Z
CHIANTI Lines in LineTAP
CHIANTI consists of a critically evaluated set of up-to-date atomic data
and derived spectral lines. At http://www.chiantidatabase.org/, many
tools for this data set are available. The present resource only
takes the line data and makes it available in standard LineTAP.
This is the state as of CHIANTI 10.1.
atomic-spectroscopy
line-positions
10.1
Dere, K.; Mason, H.; Monsignori-Fossi, B.; Landi, E.;
Young, P.; Landini, M.; Pike, D.; Bromage, G.; Bromage, B.; Del Zanna G.;
Phillips, K.J.H.
2023ApJS..268...52D
Research
Catalog
Infrared
Optical
UV
EUV
X-Ray
Please use the sentence “CHIANTI is a collaborative project
involving George Mason University, the University of Michigan (USA),
University of Cambridge (UK) and NASA Goddard Space Flight Center
(USA).” in your acknowlegement and reference 2019ApJS..241...22D
# the scups files (lines/line strength) we are parsing are
# accompanied by .elvlc files which contain the metadata of
# the levels. We parse these manually into a temporary levels
# tables that the rowmaker can use to pull its metadata from.
levelsDD = data.dd.rd.getById("get_levels")
levelsTable = api.makeData(
levelsDD,
forceSource=sourceToken.replace(".wgfa", ".elvlc")
).getPrimaryTable()
rec= {"levels": levelsTable}
with pyfits.open(sourceToken.replace(".wgfa", ".scups.fits.gz")
) as hdus:
colNames = [c.name for c in hdus[1].data.columns]
global_meta = dict(zip(colNames, hdus[1].data[0]))
rec.update(global_meta)
return rec
LVL1: 1-5
LVL2: 7-10
wavelength: 12-25
gf: 27-40
einstein_a: 42-55
trans_info: 57-
2.1798723611035e-18
@ION_ROMAN.split()[0]
@ION_N-1
int(@LVL1)
int(@LVL2)
abs(float(@wavelength))
@element
@ion_charge
def selectEnergy(row):
if row["observed_energy"] is None:
return "theory", row["theoretical_energy"]
else:
return "observation", row["observed_energy"]
method1, lowerWN = selectEnergy(@levels.getRow(@LVL1))
method2, upperWN = selectEnergy(@levels.getRow(@LVL2))
# if one method is theory, the whole thing is theory
# (we have to pick one for LineTAP)
if "theory" in [method1, method2]:
@method = "theory"
else:
@method = "observation"
if lowerWN==upperWN:
# sometimes two levels have exactly the same energy
# in the observation column. We *could* fall back
# to theory (where they then are different, naturally),
# but all that's so exotic that we simply skip the line
raise IgnoreThisRow(
"Degenerate {element}{ion_charge} {LVL1} {LVL2}".
format(**vars))
if lowerWN>upperWN:
lowerWN, upperWN = upperWN, lowerWN
@upper_energy = PLANCK_H*upperWN*100*LIGHT_C
@lower_energy = PLANCK_H*lowerWN*100*LIGHT_C
Let's make sure that the wavelength is actually
consistent with the energy levels.
BREAKAGE_LOG = open("bad-transitions.txt", "w")
A set of energy levels for internal use
level_index: 1-7
configuration: 8-38
level_label_string: 39-42
two_s: 43-47
l: 48-53
j: 54-59
observed_energy: 60-72
theoretical_energy: 74-87
/tap/sync
row = self.getFirstVOTableRow()
self.assertEqual(row["title"], 'He II 164.053 nm')
self.assertAlmostEqual(row['upper_energy'], 7.749939155470103E-18)
self.assertEqual(row["inchi"], 'InChI=1S/He/q+1')
self.assertEqual(row["ion_charge"], 1)
self.assertEqual(row['method'], 'observation')