def add_forest(self, inp, z): '''Add a Lya/Lyb forest line ''' from xastropy.igm.abs_sys.abssys_utils import GenericAbsSystem forest = GenericAbsSystem(zabs=z) # NHI NHI_dict = {'6': 12., '7': 13., '8': 14., '9': 15.} forest.NHI = NHI_dict[inp] # Lines for name in ['HI 1215', 'HI 1025', 'HI 972']: aline = AbsLine(name, linelist=self.llist[self.llist['List']]) # Attributes aline.attrib['N'] = forest.NHI aline.attrib['b'] = 20. * u.km / u.s aline.attrib['z'] = forest.zabs # Append forest.lines.append(aline) # Append to forest lines self.all_forest.append(forest)
def add_forest(self,inp,z): '''Add a Lya/Lyb forest line ''' from xastropy.igm.abs_sys.abssys_utils import GenericAbsSystem forest = GenericAbsSystem((0.*u.deg,0.*u.deg), z, [-300.,300.]*u.km/u.s) # NHI NHI_dict = {'6':12.,'7':13.,'8':14.,'9':15.} forest.NHI=NHI_dict[inp] # Lines for name in ['HI 1215','HI 1025', 'HI 972']: aline = AbsLine(name, linelist=self.llist[self.llist['List']]) # Attributes aline.attrib['N'] = 10**forest.NHI * u.cm**-2 aline.attrib['b'] = 20.*u.km/u.s aline.attrib['z'] = forest.zabs # Append forest.lines.append(aline) # Append to forest lines self.all_forest.append(forest)
def add_forest(self, inp, z): """Add a Lya/Lyb forest line """ from xastropy.igm.abs_sys.abssys_utils import GenericAbsSystem forest = GenericAbsSystem((0.0 * u.deg, 0.0 * u.deg), z, [-300.0, 300.0] * u.km / u.s) # NHI NHI_dict = {"6": 12.0, "7": 13.0, "8": 14.0, "9": 15.0} forest.NHI = NHI_dict[inp] # Lines for name in ["HI 1215", "HI 1025", "HI 972"]: aline = AbsLine(name, linelist=self.llist[self.llist["List"]]) # Attributes aline.attrib["N"] = 10 ** forest.NHI * u.cm ** -2 aline.attrib["b"] = 20.0 * u.km / u.s aline.attrib["z"] = forest.zabs # Append forest.lines.append(aline) # Append to forest lines self.all_forest.append(forest)
def test_generic(): # Start gensurvey = GenericAbsSurvey() # Systems gensys1 = GenericAbsSystem(NHI=16., zabs=1.244) gensys1.coord = SkyCoord(ra=123.1143 * u.deg, dec=-12.4321 * u.deg) gensys1.name = 'Sys1' # gensys2 = GenericAbsSystem(NHI=17., zabs=1.744) gensys2.coord = SkyCoord(ra=223.1143 * u.deg, dec=42.4321 * u.deg) gensys2.name = 'Sys2' # Combine gensurvey.nsys = 2 gensurvey._abs_sys.append(gensys1) gensurvey._abs_sys.append(gensys2) # Attribute aNHI = gensurvey.NHI np.testing.assert_allclose(aNHI, np.array([16., 17.]))
def test_generic(): # Start gensurvey = GenericAbsSurvey() # Systems gensys1 = GenericAbsSystem(NHI=16., zabs=1.244) gensys1.coord = SkyCoord(ra=123.1143*u.deg, dec=-12.4321*u.deg) gensys1.name = 'Sys1' # gensys2 = GenericAbsSystem(NHI=17., zabs=1.744) gensys2.coord = SkyCoord(ra=223.1143*u.deg, dec=42.4321*u.deg) gensys2.name = 'Sys2' # Combine gensurvey.nsys = 2 gensurvey._abs_sys.append(gensys1) gensurvey._abs_sys.append(gensys2) # Attribute aNHI = gensurvey.NHI np.testing.assert_allclose(aNHI, np.array([16.,17.]))