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)