def __init__(self, radec, zabs, vlim, NHI, **kwargs):
"""Standard init
NHI keyword is required
"""
# NHI
if NHI < 20.3:
raise ValueError("This is not a DLA! Try an LLS (or SLLS)")
# vlim
if vlim is None:
vlim = [-500., 500.]*u.km/u.s
# Generate with type
IGMSystem.__init__(self, radec, zabs, vlim, NHI=NHI, abs_type='DLA', **kwargs)
def __init__(self, radec, zabs, vlim, **kwargs):
"""Standard init
NHI keyword is required
Parameters
----------
radec : tuple or coordinate
RA/Dec of the sightline or astropy.coordinate
zabs : float
Absorption redshift
vlim : Quantity array (2)
Velocity limits of the system
Defaulted to +/- 500 km/s if None (see Prochaska et al. 2016 HDLLS)
NHI= : float, required despite being a keyword
log10 of HI column density
**kwargs : keywords
passed to IGMSystem.__init__
"""
# NHI
try:
NHI = kwargs['NHI']
except KeyError:
raise ValueError("NHI must be specified for LLSSystem")
else:
kwargs.pop('NHI')
# vlim
if vlim is None:
vlim = [-500., 500.] * u.km / u.s
# Generate with type
IGMSystem.__init__(self,
radec,
zabs,
vlim,
NHI=NHI,
abs_type='LLS',
**kwargs)
# Set tau_LL
self.tau_LL = (10.**self.NHI) * ltaa.photo_cross(
1, 1, 1 * u.Ry).to('cm**2').value
# Other
self.zpeak = None # Optical depth weighted redshift
self.ZH = 0.
self.metallicity = None # MetallicityPDF class usually
# Subsystems
self.nsub = 0
self.subsys = {}
def __init__(self, radec, zabs, vlim, **kwargs):
"""Standard init
NHI keyword is required
Parameters
----------
radec : tuple or coordinate
RA/Dec of the sightline or astropy.coordinate
zabs : float
Absorption redshift
vlim : Quantity array (2)
Velocity limits of the system
Defaulted to +/- 500 km/s if None (see Prochaska et al. 2016 HDLLS)
NHI= : float, required despite being a keyword
log10 of HI column density
**kwargs : keywords
passed to IGMSystem.__init__
"""
# NHI
try:
NHI = kwargs['NHI']
except KeyError:
raise ValueError("NHI must be specified for LLSSystem")
else:
kwargs.pop('NHI')
# vlim
if vlim is None:
vlim = [-500.,500.]*u.km/u.s
# Generate with type
IGMSystem.__init__(self, 'LLS', radec, zabs, vlim, NHI=NHI, **kwargs)
# Set tau_LL
self.tau_LL = (10.**self.NHI)*ltaa.photo_cross(1, 1, 1*u.Ry).to('cm**2').value
# Other
self.zpeak = None # Optical depth weighted redshift
self.ZH = 0.
self.metallicity = None # MetallicityPDF class usually
# Subsystems
self.nsub = 0
self.subsys = {}
def __init__(self, radec, zabs, vlim, NHI, **kwargs):
"""Standard init
NHI keyword is required
"""
# NHI
if NHI < 20.3:
raise ValueError("This is not a DLA! Try an LLS (or SLLS)")
# vlim
if vlim is None:
vlim = [-500., 500.] * u.km / u.s
# Generate with type
IGMSystem.__init__(self,
radec,
zabs,
vlim,
NHI=NHI,
abs_type='DLA',
**kwargs)
