def __init__(self, N=0.0, b=10.0, D_fraction=1.5e-5): fuf.OneDFit.__init__(self, ["N","b", "Dfrac"], rootName="LyaTransmission") self["b"] = b self["Dfrac"] = D_fraction # (Only) Einstein coefficient relevant for LyA elya = 6.258085e8 # Hydrogen self._absH = VoigtAstroP() self._absH["w0"] = 1215.67 self._absH["b"] = self["b"] self._absH["gamma"] = convertDampingConstant(elya, 1215.67) self._absH["f"] = 0.416 # Deuterium, its width is sqrt(2) times smaller self._absD = VoigtAstroP() self._absD["w0"] = 1215.34 self._absD["b"] = self["b"]/np.sqrt(2.0) self._absD["gamma"] = convertDampingConstant(elya, 1215.67) self._absD["f"] = 0.416
def __init__(self, N=0.0, b=10.0, D_fraction=1.5e-5):
fuf.OneDFit.__init__(self, ["N", "b", "Dfrac"],
rootName="LyaTransmission")
self["b"] = b
self["Dfrac"] = D_fraction
# (Only) Einstein coefficient relevant for LyA
elya = 6.258085e8
# Hydrogen
self._absH = VoigtAstroP()
self._absH["w0"] = 1215.67
self._absH["b"] = self["b"]
self._absH["gamma"] = convertDampingConstant(elya, 1215.67)
self._absH["f"] = 0.416
# Deuterium, its width is sqrt(2) times smaller
self._absD = VoigtAstroP()
self._absD["w0"] = 1215.34
self._absD["b"] = self["b"] / np.sqrt(2.0)
self._absD["gamma"] = convertDampingConstant(elya, 1215.67)
self._absD["f"] = 0.416
def sanity_example(self):
"""
Checking example for converting damping constant into line width.
"""
from PyAstronomy import pyasl
# Einstein coefficient relevant for hydrogen LyA
gLya = 6.258085e8
print "Width of H LyA line at 1215.67 A = %e cm" % \
pyasl.convertDampingConstant(gLya, 1215.67)